Inheritance diagram for trkf::TrackKalmanCheater:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	TrackKalmanCheater (fhicl::ParameterSet const &pset)

virtual	~TrackKalmanCheater ()
	Destructor. More...

virtual void	reconfigure (fhicl::ParameterSet const &pset)

virtual void	produce (art::Event &e)

virtual void	beginJob ()
	Begin job method. More...

virtual void	endJob ()
	End job method. More...

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () const

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Attributes
bool	fHist
	Make histograms. More...

KalmanFilterAlg	fKFAlg
	Kalman filter algorithm. More...

SpacePointAlg	fSpacePointAlg
	Space point algorithm. More...

bool	fUseClusterHits
	Use cluster hits or all hits? More...

std::string	fHitModuleLabel
	Unclustered Hits. More...

std::string	fClusterModuleLabel
	Clustered Hits. More...

double	fMaxTcut
	Maximum delta ray energy in MeV for restricted dE/dx. More...

const Propagator *	fProp
	Propagator. More...

TH1F *	fHIncChisq
	Incremental chisquare. More...

TH1F *	fHPull
	Hit pull. More...

int	fNumEvent
	Number of events seen. More...

int	fNumTrack
	Number of tracks produced. More...

Detailed Description

Definition at line 54 of file TrackKalmanCheater_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

trkf::TrackKalmanCheater::TrackKalmanCheater ( fhicl::ParameterSet const & pset )

explicit

Constructor.

Arguments:

p - Fcl parameters.

Definition at line 107 of file TrackKalmanCheater_module.cc.

References fClusterModuleLabel, fHitModuleLabel, fUseClusterHits, and reconfigure().

   : fHist(false)
   , fKFAlg(pset.get<fhicl::ParameterSet>("KalmanFilterAlg"))
   , fSpacePointAlg(pset.get<fhicl::ParameterSet>("SpacePointAlg"))
   , fUseClusterHits(false)
   , fMaxTcut(0.)
   , fProp(0)
   , fHIncChisq(0)
   , fHPull(0)
   , fNumEvent(0)
   , fNumTrack(0)
 {
   reconfigure(pset);
   produces<std::vector<recob::Track>                   >();
   produces<std::vector<recob::SpacePoint>              >();
   produces<art::Assns<recob::Track, recob::Hit>        >();
   produces<art::Assns<recob::Track, recob::SpacePoint> >();
   produces<art::Assns<recob::SpacePoint, recob::Hit>   >();
 
   // Report.
 
   mf::LogInfo("TrackKalmanCheater") 
     << "TrackKalmanCheater configured with the following parameters:\n"
     << "  UseClusterHits = " << fUseClusterHits << "\n"
     << "  HitModuleLabel = " << fHitModuleLabel << "\n"
     << "  ClusterModuleLabel = " << fClusterModuleLabel;
 }

trkf::TrackKalmanCheater::~TrackKalmanCheater ( )

virtual

Destructor.

Definition at line 137 of file TrackKalmanCheater_module.cc.

138 {}

Member Function Documentation

void trkf::TrackKalmanCheater::beginJob ( )

virtual

Begin job method.

Reimplemented from art::EDProducer.

Definition at line 163 of file TrackKalmanCheater_module.cc.

References dir, fHIncChisq, fHist, fHPull, art::TFileDirectory::make(), and art::TFileDirectory::mkdir().

 {
   if(fHist) {
 
     // Book histograms.
 
     art::ServiceHandle<art::TFileService> tfs;
     art::TFileDirectory dir = tfs->mkdir("hitkalman", "TrackKalmanCheater histograms");
 
     fHIncChisq = dir.make<TH1F>("IncChisq", "Incremental Chisquare", 100, 0., 20.);
     fHPull = dir.make<TH1F>("Pull", "Hit Pull", 100, -10., 10.);
   }
 }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

void trkf::TrackKalmanCheater::endJob ( )

virtual

End job method.

Reimplemented from art::EDProducer.

Definition at line 498 of file TrackKalmanCheater_module.cc.

References fNumEvent, and fNumTrack.

 {
   mf::LogInfo("TrackKalmanCheater") 
     << "TrackKalmanCheater statistics:\n"
     << "  Number of events = " << fNumEvent << "\n"
     << "  Number of tracks created = " << fNumTrack;
 }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void trkf::TrackKalmanCheater::produce ( art::Event & evt )

virtual

Produce method.

Arguments:

e - Art event.

This method extracts Hit from the event and produces and adds Track objects.

Implements art::EDProducer.

Definition at line 187 of file TrackKalmanCheater_module.cc.

References trkf::Surface::BACKWARD, art::PtrVector< T >::begin(), sim::ParticleList::begin(), trkf::KalmanFilterAlg::buildTrack(), trkf::SpacePointAlg::clearHitMap(), util::CreateAssn(), dir, art::PtrVector< T >::empty(), art::PtrVector< T >::end(), sim::ParticleList::end(), fClusterModuleLabel, fHIncChisq, fHist, fHitModuleLabel, fHPull, trkf::KHitContainerWireX::fill(), trkf::KGTrack::fillHits(), trkf::SpacePointAlg::fillSpacePoints(), trkf::KGTrack::fillTrack(), trkf::KalmanFilterAlg::fitMomentum(), fKFAlg, fNumEvent, fNumTrack, trkf::Surface::FORWARD, trkf::Propagator::FORWARD, fProp, fSpacePointAlg, fUseClusterHits, trkf::SpacePointAlg::getAssociatedHits(), art::DataViewImpl::getByLabel(), trkf::KHitTrack::getHit(), trkf::KHit< N >::getResError(), trkf::KHit< N >::getResVector(), trkf::KGTrack::getTrackMap(), hits(), cheat::BackTrackerService::HitToTrackIDEs(), ipart, mf::isDebugEnabled(), art::Handle< T >::isValid(), simb::MCParticle::Momentum(), part, cheat::ParticleInventoryService::ParticleList(), simb::MCParticle::PdgCode(), geo::PlaneID::Plane, simb::MCParticle::Position(), art::Event::productGetter(), art::PtrVector< T >::push_back(), art::Event::put(), art::PtrVector< T >::reserve(), trkf::KalmanFilterAlg::setPlane(), art::PtrVector< T >::size(), trkf::KalmanFilterAlg::smoothTrackIter(), simb::MCParticle::StatusCode(), simb::MCParticle::TrackId(), trkf::KGTrack::TrackMap(), trkf::KalmanFilterAlg::updateMomentum(), lar::dump::vector(), recob::Hit::WireID(), x, y, and z.

 {
   ++fNumEvent;
 
   // Make a collection of tracks, plus associations, that will
   // eventually be inserted into the event.
 
   std::unique_ptr<std::vector<recob::Track> > tracks(new std::vector<recob::Track>);
   std::unique_ptr< art::Assns<recob::Track, recob::Hit> > th_assn(new art::Assns<recob::Track, recob::Hit>);
   std::unique_ptr< art::Assns<recob::Track, recob::SpacePoint> > tsp_assn(new art::Assns<recob::Track, recob::SpacePoint>);
 
   // Make a collection of space points, plus associations, that will
   // be inserted into the event.
 
   std::unique_ptr<std::vector<recob::SpacePoint> > spts(new std::vector<recob::SpacePoint>);
   std::unique_ptr< art::Assns<recob::SpacePoint, recob::Hit> > sph_assn(new art::Assns<recob::SpacePoint, recob::Hit>);
 
   // Make a collection of KGTracks where we will save our results.
 
   std::deque<KGTrack> kalman_tracks;
 
   // Get Services.
 
   art::ServiceHandle<geo::Geometry> geom;
   art::ServiceHandle<cheat::BackTrackerService> bt_serv;
   art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
 
   // Reset space point algorithm.
 
   fSpacePointAlg.clearHitMap();
 
 
   // Get Hits.
 
   art::PtrVector<recob::Hit> hits;
 
   if(fUseClusterHits) {
 
     // Get clusters.
 
     art::Handle< std::vector<recob::Cluster> > clusterh;
     evt.getByLabel(fClusterModuleLabel, clusterh);
 
     // Get hits from all clusters.
     art::FindManyP<recob::Hit> fm(clusterh, evt, fClusterModuleLabel);
 
     if(clusterh.isValid()) {
       int nclus = clusterh->size();
 
       for(int i = 0; i < nclus; ++i) {
         art::Ptr<recob::Cluster> pclus(clusterh, i);
         std::vector< art::Ptr<recob::Hit> > clushits = fm.at(i);
         int nhits = clushits.size();
         hits.reserve(hits.size() + nhits);
 
         for(std::vector< art::Ptr<recob::Hit> >::const_iterator ihit = clushits.begin();
             ihit != clushits.end(); ++ihit) {
           hits.push_back(*ihit);
         }
       }
     }
   }
   else {
 
     // Get unclustered hits.
 
     art::Handle< std::vector<recob::Hit> > hith;
     evt.getByLabel(fHitModuleLabel, hith);
     if(hith.isValid()) {
       int nhits = hith->size();
       hits.reserve(nhits);
 
       for(int i = 0; i < nhits; ++i)
         hits.push_back(art::Ptr<recob::Hit>(hith, i));
     }
   }
 
   // Sort hits into separate PtrVectors based on track id.
 
   std::map<int, art::PtrVector<recob::Hit> > hitmap;
 
   // Loop over hits.
 
   for(art::PtrVector<recob::Hit>::const_iterator ihit = hits.begin();
         ihit != hits.end(); ++ihit) {
     //const recob::Hit& hit = **ihit;
 
     // Get track ids for this hit.
 
     std::vector<sim::TrackIDE> tids = bt_serv->HitToTrackIDEs(*ihit);
 
     // Loop over track ids.
 
     for(std::vector<sim::TrackIDE>::const_iterator itid = tids.begin();
         itid != tids.end(); ++itid) {
       int trackID = itid->trackID;
 
       // Add hit to PtrVector corresponding to this track id.
 
       hitmap[trackID].push_back(*ihit);
     }
   }
 
   // Extract geant mc particles.
 
   sim::ParticleList const& plist = pi_serv->ParticleList();
 
   // Loop over geant particles.
 
   {
     // use mf::LogDebug instead of LOG_DEBUG because we reuse it in many lines;
     // insertions are protected by mf::isDebugEnabled()
     mf::LogDebug log("TrackKalmanCheater");
     for(sim::ParticleList::const_iterator ipart = plist.begin();
         ipart != plist.end(); ++ipart) {
       const simb::MCParticle* part = (*ipart).second;
       if (!part) throw cet::exception("TrackKalmanCheater") << "no particle!\n";
       int pdg = part->PdgCode();
 
       // Ignore everything except stable charged nonshowering particles.
 
       int apdg = std::abs(pdg);
       if(apdg == 13 ||     // Muon
          apdg == 211 ||    // Charged pion
          apdg == 321 ||    // Charged kaon
          apdg == 2212) {   // (Anti)proton
 
         int trackid = part->TrackId();
         int stat = part->StatusCode();
         int nhit = 0;
         if(hitmap.count(trackid) != 0)
           nhit = hitmap[trackid].size();
         const TLorentzVector& pos = part->Position();
         const TLorentzVector& mom = part->Momentum();
         
         if (mf::isDebugEnabled()) {
           log << "Trackid=" << trackid 
               << ", pdgid=" << pdg 
               << ", status=" << stat
               << ", NHit=" << nhit << "\n"
               << "  x = " << pos.X()
               << ", y = " << pos.Y()
               << ", z = " << pos.Z() << "\n"
               << "  px= " << mom.Px()
               << ", py= " << mom.Py()
               << ", pz= " << mom.Pz() << "\n";
         } // if debugging
         
         double x = pos.X();
         double y = pos.Y();
         double z = pos.Z();
         double px = mom.Px();
         double py = mom.Py();
         double pz = mom.Pz();
         double p = std::sqrt(px*px + py*py + pz*pz);
 
         if(nhit > 0 && pz != 0.) {
           const art::PtrVector<recob::Hit>& trackhits = hitmap[trackid];
           if (trackhits.empty())
             throw cet::exception("TrackKalmanCheater") << "No hits in track\n";
 
           // Make a seed track (KTrack).
 
           std::shared_ptr<const Surface> psurf(new SurfYZPlane(0., y, z, 0.));
           TrackVector vec(5);
           vec(0) = x;
           vec(1) = 0.;
           vec(2) = px / pz;
           vec(3) = py / pz;
           vec(4) = 1. / p;
           Surface::TrackDirection dir = (pz > 0. ? Surface::FORWARD : Surface::BACKWARD);
           KTrack trk(psurf, vec, dir, pdg);
 
           // Fill KHitContainer with hits.
 
           KHitContainerWireX cont;
           cont.fill(trackhits, -1);
 
           // Count hits in each plane.  Set the preferred plane to be
           // the one with the most hits.
 
           std::vector<unsigned int> planehits(3, 0);
           for(art::PtrVector<recob::Hit>::const_iterator ih = trackhits.begin();
               ih != trackhits.end(); ++ih) {
             const recob::Hit& hit = **ih;
             unsigned int plane = hit.WireID().Plane;
 
             if (plane >= planehits.size()) {
               throw cet::exception("TrackKalmanCheater") << "plane " << plane << "...\n"; }
             ++planehits[plane];
           }
           unsigned int prefplane = 0;
           for(unsigned int i=0; i<planehits.size(); ++i) {
             if(planehits[i] > planehits[prefplane])
               prefplane = i;
           }
           if (mf::isDebugEnabled())
             log << "Preferred plane = " << prefplane << "\n";
 
           // Build and smooth track.
 
           KGTrack trg(prefplane);
           fKFAlg.setPlane(prefplane);
           bool ok = fKFAlg.buildTrack(trk, trg, fProp, Propagator::FORWARD, cont, false);
           if(ok) {
             ok = fKFAlg.smoothTrackIter(5, trg, fProp);
             if(ok) {
               KETrack tremom;
               bool pok = fKFAlg.fitMomentum(trg, fProp, tremom);
               if(pok)
                 fKFAlg.updateMomentum(tremom, fProp, trg);
               ++fNumTrack;
               kalman_tracks.push_back(trg);
             }
           }
           if (mf::isDebugEnabled())
             log << (ok? "Build track succeeded.": "Build track failed.") << "\n";
           
         }
       }
     }
   }
 
   // Fill histograms.
 
   if(fHist) {
 
     // First loop over tracks.
 
     for(std::deque<KGTrack>::const_iterator k = kalman_tracks.begin();
         k != kalman_tracks.end(); ++k) {
       const KGTrack& trg = *k;
 
       // Loop over measurements in this track.
 
       const std::multimap<double, KHitTrack>& trackmap = trg.getTrackMap();
       for(std::multimap<double, KHitTrack>::const_iterator ih = trackmap.begin();
           ih != trackmap.end(); ++ih) {
         const KHitTrack& trh = (*ih).second;
         const std::shared_ptr<const KHitBase>& hit = trh.getHit();
         double chisq = hit->getChisq();
         fHIncChisq->Fill(chisq);
         const KHit<1>* ph1 = dynamic_cast<const KHit<1>*>(&*hit);
         if(ph1 != 0) {
           double pull = ph1->getResVector()(0) / std::sqrt(ph1->getResError()(0, 0));
           fHPull->Fill(pull);
         }
       }
     }
   }
 
   // Process Kalman filter tracks into persistent objects.
 
   tracks->reserve(kalman_tracks.size());
   for(std::deque<KGTrack>::const_iterator k = kalman_tracks.begin();
       k != kalman_tracks.end(); ++k) {
     const KGTrack& kalman_track = *k;
 
     // Add Track object to collection.
 
     tracks->push_back(recob::Track());
     kalman_track.fillTrack(tracks->back(), tracks->size() - 1);
 
     // Make Track to Hit associations.  
 
     art::PtrVector<recob::Hit> trhits;
     std::vector<unsigned int> hittpindex;
     kalman_track.fillHits(hits, hittpindex);
     util::CreateAssn(*this, evt, *tracks, trhits, *th_assn, tracks->size()-1);
 
     // Make space points from this track.
 
     int nspt = spts->size();
     fSpacePointAlg.fillSpacePoints(*spts, kalman_track.TrackMap());
 
     // Associate newly created space points with hits.
     // Also associate track with newly created space points.
 
     art::PtrVector<recob::SpacePoint> sptvec;
 
     // Loop over newly created space points.
       
     for(unsigned int ispt = nspt; ispt < spts->size(); ++ispt) {
       const recob::SpacePoint& spt = (*spts)[ispt];
       art::ProductID sptid = getProductID<std::vector<recob::SpacePoint> >();
       art::Ptr<recob::SpacePoint> sptptr(sptid, ispt, evt.productGetter(sptid));
       sptvec.push_back(sptptr);
 
       // Make space point to hit associations.
 
       const art::PtrVector<recob::Hit>& sphits = 
         fSpacePointAlg.getAssociatedHits(spt);
       util::CreateAssn(*this, evt, *spts, sphits, *sph_assn, ispt);
     }
 
     // Make track to space point associations.
 
     util::CreateAssn(*this, evt, *tracks, sptvec, *tsp_assn, tracks->size()-1);
   }
 
   // Add tracks and associations to event.
 
   evt.put(std::move(tracks));
   evt.put(std::move(spts));
   evt.put(std::move(th_assn));
   evt.put(std::move(tsp_assn));
   evt.put(std::move(sph_assn));
 }

void trkf::TrackKalmanCheater::reconfigure ( fhicl::ParameterSet const & pset )

virtual

Reconfigure method.

Arguments:

p - Fcl parameter set.

Definition at line 147 of file TrackKalmanCheater_module.cc.

References fClusterModuleLabel, fHist, fHitModuleLabel, fKFAlg, fMaxTcut, fProp, fSpacePointAlg, fUseClusterHits, fhicl::ParameterSet::get(), trkf::KalmanFilterAlg::reconfigure(), and trkf::SpacePointAlg::reconfigure().

Referenced by TrackKalmanCheater().

 {
   fHist = pset.get<bool>("Hist");
   fKFAlg.reconfigure(pset.get<fhicl::ParameterSet>("KalmanFilterAlg"));
   fSpacePointAlg.reconfigure(pset.get<fhicl::ParameterSet>("SpacePointAlg"));
   fUseClusterHits = pset.get<bool>("UseClusterHits");
   fHitModuleLabel = pset.get<std::string>("HitModuleLabel");
   fClusterModuleLabel = pset.get<std::string>("ClusterModuleLabel");
   fMaxTcut = pset.get<double>("MaxTcut");
   if(fProp != 0)
     delete fProp;
   fProp = new PropYZPlane(fMaxTcut, true);
 }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }