Inheritance diagram for trkf::TrackKalmanCheater:

Public Types
using	ModuleType = EDProducer

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

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

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
ConsumesCollector &	consumesCollector ()

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

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

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

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

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

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

Private Member Functions
void	produce (art::Event &e) override

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

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

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...

TH1F *	fHIncChisq {nullptr}
	Incremental chisquare. More...

TH1F *	fHPull {nullptr}
	Hit pull. More...

int	fNumEvent {0}
	Number of events seen. More...

int	fNumTrack {0}
	Number of tracks produced. More...

Detailed Description

Definition at line 65 of file TrackKalmanCheater_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

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

explicit

Constructor.

Arguments:

p - Fcl parameters.

Definition at line 106 of file TrackKalmanCheater_module.cc.

References fClusterModuleLabel, fHist, fHitModuleLabel, fKFAlg, fMaxTcut, fSpacePointAlg, and fUseClusterHits.

   : EDProducer{pset}
   , fHist(pset.get<bool>("Hist"))
   , fKFAlg(pset.get<fhicl::ParameterSet>("KalmanFilterAlg"))
   , fSpacePointAlg(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"))
 {
   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;
 }

Member Function Documentation

void trkf::TrackKalmanCheater::beginJob ( )

overrideprivatevirtual

Begin job method.

Reimplemented from art::EDProducer.

Definition at line 133 of file TrackKalmanCheater_module.cc.

References dir, fHIncChisq, fHist, and fHPull.

 {
   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 BT>

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

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

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

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::doEvent	(	EventPrincipal &	ep,
		ModuleContext const &	mc,
		std::atomic< std::size_t > &	counts_run,
		std::atomic< std::size_t > &	counts_passed,
		std::atomic< std::size_t > &	counts_failed
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void trkf::TrackKalmanCheater::endJob ( )

overrideprivatevirtual

End job method.

Reimplemented from art::EDProducer.

Definition at line 449 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;
 }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

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

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

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

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

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

overrideprivatevirtual

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 157 of file TrackKalmanCheater_module.cc.

References util::abs(), 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::SpacePointAlg::fillSpacePoints(), trkf::KalmanFilterAlg::fitMomentum(), fKFAlg, fMaxTcut, fNumEvent, fNumTrack, trkf::Surface::FORWARD, trkf::Propagator::FORWARD, fSpacePointAlg, fUseClusterHits, trkf::SpacePointAlg::getAssociatedHits(), art::ProductRetriever::getByLabel(), trkf::KHitTrack::getHit(), art::ProductRetriever::getProductID(), trkf::KHit< N >::getResError(), trkf::KHit< N >::getResVector(), 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::ProductRetriever::productGetter(), art::PtrVector< T >::push_back(), art::Event::put(), art::PtrVector< T >::reserve(), trkf::KalmanFilterAlg::setPlane(), art::PtrVector< T >::size(), util::size(), trkf::KalmanFilterAlg::smoothTrackIter(), simb::MCParticle::StatusCode(), simb::MCParticle::TrackId(), 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 const> geom;
   art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
   art::ServiceHandle<cheat::ParticleInventoryService const> pi_serv;
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
 
   // 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(clockData, *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();
 
   auto const det_prop =
     art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt);
 
   PropYZPlane const prop{det_prop, fMaxTcut, true};
 
   // Loop over geant particles.
 
   {
     // use mf::LogDebug instead of MF_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 (not accepted_particle(apdg)) { continue; }
 
       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::hypot(px, py, 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(det_prop, 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);
         if (fKFAlg.buildTrack(trk, trg, prop, Propagator::FORWARD, cont, false)) {
           if (fKFAlg.smoothTrackIter(5, trg, prop)) {
             KETrack tremom;
             if (fKFAlg.fitMomentum(trg, prop, tremom)) { fKFAlg.updateMomentum(tremom, prop, trg); }
             ++fNumTrack;
             kalman_tracks.push_back(trg);
           }
         }
       }
     }
   }
 
   // Fill histograms.
 
   if (fHist) {
 
     // First loop over tracks.
 
     for (auto const& trg : kalman_tracks) {
 
       // 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 (auto const& kalman_track : kalman_tracks) {
     tracks->push_back(recob::Track());
     kalman_track.fillTrack(det_prop, 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(evt, *tracks, trhits, *th_assn, tracks->size() - 1);
 
     // Make space points from this track.
 
     int nspt = spts->size();
     fSpacePointAlg.fillSpacePoints(det_prop, *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 = evt.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(evt, *spts, sphits, *sph_assn, ispt);
     }
 
     // Make track to space point associations.
 
     util::CreateAssn(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 art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }