Inheritance diagram for trkf::SeedAna:

Classes
struct	MCHists

struct	RecoHists

Public Types
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

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

virtual	~SeedAna ()

void	analyze (const art::Event &evt)

void	endJob ()

std::string	workerType () const

bool	modifiesEvent () const

void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)

std::string const &	processName () const

bool	wantAllEvents () const

bool	wantEvent (Event const &e)

fhicl::ParameterSetID	selectorConfig () const

art::Handle< art::TriggerResults >	getTriggerResults (Event const &e) 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

detail::CachedProducts &	cachedProducts ()

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

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Attributes
std::string	fSeedModuleLabel

std::string	fMCTrackModuleLabel

int	fDump

double	fMinMCKE

double	fMinMCLen

double	fMatchColinearity

double	fMatchDisp

bool	fIgnoreSign

std::map< int, MCHists >	fMCHistMap

std::map< int, RecoHists >	fRecoHistMap

int	fNumEvent

Detailed Description

Definition at line 233 of file SeedAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

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

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 641 of file SeedAna_module.cc.

References fDump, fMCTrackModuleLabel, fMinMCKE, fMinMCLen, and fSeedModuleLabel.

     : EDAnalyzer(pset)
     , fSeedModuleLabel(pset.get<std::string>("SeedModuleLabel"))
     , fMCTrackModuleLabel(pset.get<std::string>("MCTrackModuleLabel"))
     , fDump(pset.get<int>("Dump"))
     , fMinMCKE(pset.get<double>("MinMCKE"))
     , fMinMCLen(pset.get<double>("MinMCLen"))
     , fMatchColinearity(pset.get<double>("MatchColinearity"))
     , fMatchDisp(pset.get<double>("MatchDisp"))
     , fIgnoreSign(pset.get<bool>("IgnoreSign"))
     , fNumEvent(0)
   {
 
     // Report.
 
     mf::LogInfo("SeedAna") 
       << "SeedAna configured with the following parameters:\n"
       << "  SeedModuleLabel = " << fSeedModuleLabel << "\n"
       << "  MCTrackModuleLabel = " << fMCTrackModuleLabel << "\n"
       << "  Dump = " << fDump << "\n"
       << "  MinMCKE = " << fMinMCKE << "\n"
       << "  MinMCLen = " << fMinMCLen;
   }

trkf::SeedAna::~SeedAna ( )

virtual

Definition at line 670 of file SeedAna_module.cc.

674 {}

Member Function Documentation

void trkf::SeedAna::analyze ( const art::Event & evt )

Definition at line 676 of file SeedAna_module.cc.

   {
     ++fNumEvent;
 
     // Optional dump stream.
 
     std::unique_ptr<mf::LogInfo> pdump;
     if(fDump > 0) {
       --fDump;
       pdump = std::unique_ptr<mf::LogInfo>(new mf::LogInfo("TrackAna"));
     }
 
     // Make sure histograms are booked.
 
     bool mc = !evt.isRealData();
 
     // Get seed handle.
 
     art::Handle< std::vector<recob::Seed> > seedh;
     evt.getByLabel(fSeedModuleLabel, seedh);
 
     const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
     
     // Seed->mc track matching map.
 
     std::map<const recob::Seed*, int> seedmap;
 
     if(mc) {
 
       // Get MCTracks.
 
       art::Handle< std::vector<sim::MCTrack> > mctrackh;
       evt.getByLabel(fMCTrackModuleLabel, mctrackh);
 
       // Dump MCTracks.
 
       if(pdump) {
         *pdump << "MC Tracks\n"
                << "       Id   pdg           x         y         z          dx        dy        dz           p\n"
                << "-------------------------------------------------------------------------------------------\n";
       }
 
       // Loop over mc tracks, and fill histograms that depend only
       // on mc particles.
 
       for(std::vector<sim::MCTrack>::const_iterator imctrk = mctrackh->begin();
           imctrk != mctrackh->end(); ++imctrk) {
         const sim::MCTrack& mctrk = *imctrk;
         int pdg = mctrk.PdgCode();
         if(fIgnoreSign)
           pdg = std::abs(pdg);
 
         // 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
 
           // Apply minimum energy cut.
 
           if(mctrk.Start().E() >= mctrk.Start().Momentum().Mag() + 1000.*fMinMCKE) {
 
             // Calculate the x offset due to nonzero mc particle time.
 
             double mctime = mctrk.Start().T();                                // nsec
             double mcdx = mctime * 1.e-3 * detprop->DriftVelocity();  // cm
 
             // Calculate the length of this mc particle inside the fiducial volume.
 
             TVector3 mcstart;
             TVector3 mcend;
             TVector3 mcstartmom;
             TVector3 mcendmom;
             double plen = length(mctrk, mcdx, mcstart, mcend, mcstartmom, mcendmom);
 
             // Apply minimum fiducial length cut.  Always reject particles that have
             // zero length in the tpc regardless of the configured cut.
 
             if(plen > 0. && plen > fMinMCLen) {
 
               // Dump MC particle information here.
 
               if(pdump) {
                 double pstart = mcstartmom.Mag();
                 double pend = mcendmom.Mag();
                 *pdump << "\nOffset"
                        << std::setw(3) << mctrk.TrackID()
                        << std::setw(6) << mctrk.PdgCode()
                        << "  " 
                        << std::fixed << std::setprecision(2) 
                        << std::setw(10) << mcdx
                        << "\nStart " 
                        << std::setw(3) << mctrk.TrackID()
                        << std::setw(6) << mctrk.PdgCode()
                        << "  " 
                        << std::fixed << std::setprecision(2) 
                        << std::setw(10) << mcstart[0]
                        << std::setw(10) << mcstart[1]
                        << std::setw(10) << mcstart[2];
                 if(pstart > 0.) {
                   *pdump << "  "
                          << std::fixed << std::setprecision(3) 
                          << std::setw(10) << mcstartmom[0]/pstart
                          << std::setw(10) << mcstartmom[1]/pstart
                          << std::setw(10) << mcstartmom[2]/pstart;
                 }
                 else
                   *pdump << std::setw(32) << " ";
                 *pdump << std::setw(12) << std::fixed << std::setprecision(3) << pstart;
                 *pdump << "\nEnd   " 
                        << std::setw(3) << mctrk.TrackID()
                        << std::setw(6) << mctrk.PdgCode()
                        << "  " 
                        << std::fixed << std::setprecision(2)
                        << std::setw(10) << mcend[0]
                        << std::setw(10) << mcend[1]
                        << std::setw(10) << mcend[2];
                 if(pend > 0.01) {
                   *pdump << "  " 
                          << std::fixed << std::setprecision(3) 
                          << std::setw(10) << mcendmom[0]/pend
                          << std::setw(10) << mcendmom[1]/pend
                          << std::setw(10) << mcendmom[2]/pend;
                 }
                 else
                   *pdump << std::setw(32)<< " ";
                 *pdump << std::setw(12) << std::fixed << std::setprecision(3) << pend << "\n";
               }
 
               // Fill histograms.
 
               if(fMCHistMap.count(pdg) == 0) {
                 std::ostringstream ostr;
                 ostr << "MC" << (fIgnoreSign ? "All" : (pdg > 0 ? "Pos" : "Neg")) << std::abs(pdg);
                 fMCHistMap[pdg] = MCHists(ostr.str());
               }
               const MCHists& mchists = fMCHistMap[pdg];
 
               double mctheta_xz = std::atan2(mcstartmom.X(), mcstartmom.Z());
               double mctheta_yz = std::atan2(mcstartmom.Y(), mcstartmom.Z());
 
               mchists.fHmcstartx->Fill(mcstart.X());
               mchists.fHmcstarty->Fill(mcstart.Y());
               mchists.fHmcstartz->Fill(mcstart.Z());
               mchists.fHmcendx->Fill(mcend.X());
               mchists.fHmcendy->Fill(mcend.Y());
               mchists.fHmcendz->Fill(mcend.Z());
               mchists.fHmctheta->Fill(mcstartmom.Theta());
               mchists.fHmcphi->Fill(mcstartmom.Phi());
               mchists.fHmctheta_xz->Fill(mctheta_xz);
               mchists.fHmctheta_yz->Fill(mctheta_yz);
               mchists.fHmcmom->Fill(mcstartmom.Mag());
               mchists.fHmclen->Fill(plen);
 
               // Loop over seeds and do matching.
 
               int nmatch = 0;
               if(seedh.isValid()) {
 
                 // Loop over seeds.
 
                 int nseed = seedh->size();
                 for(int i = 0; i < nseed; ++i) {
                   art::Ptr<recob::Seed> pseed(seedh, i);
                   const recob::Seed& seed = *pseed;
                   if(seedmap.count(&seed) == 0)
                     seedmap[&seed] = -1;
 
                   // Get parameters of this seed.
 
                   TVector3 pos;
                   TVector3 dir;
                   double err[3];
                   seed.GetPoint(&pos[0], err);
                   seed.GetDirection(&dir[0], err);
 
                   // Calculate the global-to-local rotation matrix.
                   // Copied from Track.cxx.
 
                   TMatrixD rot(3,3);
                   double dirmag = dir.Mag();
                   double diryz = std::sqrt(dir.Y()*dir.Y() + dir.Z()*dir.Z());
 
                   double sinth = dir.X() / dirmag;
                   double costh = diryz / dirmag;
                   double sinphi = 0.;
                   double cosphi = 1.;
                   if(diryz != 0) {
                     sinphi = -dir.Y() / diryz;
                     cosphi = dir.Z() / diryz;
                   }
                   rot(0,0) = costh;
                   rot(1,0) = 0.;
                   rot(2,0) = sinth;
                   rot(0,1) = sinth * sinphi;
                   rot(1,1) = cosphi;
                   rot(2,1) = -costh * sinphi;
                   rot(0,2) = -sinth * cosphi;
                   rot(1,2) = sinphi;
                   rot(2,2) = costh * cosphi;
 
                   // Get best matching mc trajectory point.
 
                   int itraj = mcmatch(mctrk, seed);
                   if(itraj >= 0) {
 
                     // Get mc relative position and direction at matching trajectory point.
 
                     TVector3 mcpos = mctrk[itraj].Position().Vect() - pos;
                     TVector3 mcmom = mctrk[itraj].Momentum().Vect();
                     mcpos[0] += mcdx;
 
                     // Rotate the momentum and position to the
                     // seed-local coordinate system.
 
                     TVector3 mcmoml = rot * mcmom;
                     TVector3 mcposl = rot * mcpos;
 
                     if(mcmoml.Z() < 0.)
                       mcmoml = -mcmoml;
                     double costh = mcmoml.Z() / mcmoml.Mag();
 
                     double u = mcposl.X();
                     double v = mcposl.Y();
                     double w = mcposl.Z();
 
                     double pu = mcmoml.X();
                     double pv = mcmoml.Y();
                     double pw = mcmoml.Z();
 
                     double dudw = pu / pw;
                     double dvdw = pv / pw;
 
                     double u0 = u - w * dudw;
                     double v0 = v - w * dvdw;
                     double uv0 = std::sqrt(u0*u0 + v0*v0);
 
                     // Fill matching histograms.
 
                     mchists.fHduvcosth->Fill(costh, uv0);
                     if(std::abs(uv0) < fMatchDisp) {
 
                       // Fill slope matching histograms.
 
                       mchists.fHmcdudw->Fill(dudw);
                       mchists.fHmcdvdw->Fill(dvdw);
                     }
                     mchists.fHcosth->Fill(costh);
                     if(costh > fMatchColinearity) {
 
                       // Fill displacement matching histograms.
 
                       mchists.fHmcu->Fill(u0);
                       mchists.fHmcv->Fill(v0);
                       mchists.fHmcw->Fill(w);
 
                       if(std::abs(uv0) < fMatchDisp) {
 
                         // Now we have passed all matching cuts and we have a matching
                         // mc particle + seed pair.
 
                         ++nmatch;
                         seedmap[&seed] = mctrk.TrackID();
 
                         // Fill matched seed histograms (seed multiplicity).
 
                         mchists.fHmstartx->Fill(mcstart.X());
                         mchists.fHmstarty->Fill(mcstart.Y());
                         mchists.fHmstartz->Fill(mcstart.Z());
                         mchists.fHmendx->Fill(mcend.X());
                         mchists.fHmendy->Fill(mcend.Y());
                         mchists.fHmendz->Fill(mcend.Z());
                         mchists.fHmtheta->Fill(mcstartmom.Theta());
                         mchists.fHmphi->Fill(mcstartmom.Phi());
                         mchists.fHmtheta_xz->Fill(mctheta_xz);
                         mchists.fHmtheta_yz->Fill(mctheta_yz);
                         mchists.fHmmom->Fill(mcstartmom.Mag());
                         mchists.fHmlen->Fill(plen);
                       }
                     }
                   }
                 }
 
                 // If we found at least one matched seed, fill good
                 // particle histograms.
 
                 if(nmatch > 0) {
                   mchists.fHgstartx->Fill(mcstart.X());
                   mchists.fHgstarty->Fill(mcstart.Y());
                   mchists.fHgstartz->Fill(mcstart.Z());
                   mchists.fHgendx->Fill(mcend.X());
                   mchists.fHgendy->Fill(mcend.Y());
                   mchists.fHgendz->Fill(mcend.Z());
                   mchists.fHgtheta->Fill(mcstartmom.Theta());
                   mchists.fHgphi->Fill(mcstartmom.Phi());
                   mchists.fHgtheta_xz->Fill(mctheta_xz);
                   mchists.fHgtheta_yz->Fill(mctheta_yz);
                   mchists.fHgmom->Fill(mcstartmom.Mag());
                   mchists.fHglen->Fill(plen);
                 }
               }
             }
           }
         }
       }
     }
 
     // Loop over seeds and fill reco-only seed histograms.
 
     if(seedh.isValid()) {
 
       // Loop over seeds.
 
       int nseed = seedh->size();
 
       if(nseed > 0 && pdump != 0) {
         *pdump << "\nReconstructed Seeds\n"
                << "           MCid           x         y         z          dx        dy        dz           p\n"
                << "-------------------------------------------------------------------------------------------\n";
       }
 
       for(int i = 0; i < nseed; ++i) {
         art::Ptr<recob::Seed> pseed(seedh, i);
         const recob::Seed& seed = *pseed;
 
         // Fill histograms involving reco seeds only.
 
         TVector3 pos;
         TVector3 dir;
         double err[3];
         seed.GetPoint(&pos[0], err);
         seed.GetDirection(&dir[0], err);
         double mdir = dir.Mag();
         if(mdir != 0.) {
           dir *= (1./mdir);
         }
 
         double dpos = bdist(pos);
         double theta_xz = std::atan2(dir.X(), dir.Z());
         double theta_yz = std::atan2(dir.Y(), dir.Z());
 
         // Dump seed information here.
 
         if(pdump) {
           int mcid = seedmap[&seed];
           *pdump << std::setw(15) << mcid
                  << "  "
                  << std::fixed << std::setprecision(2) 
                  << std::setw(10) << pos[0]
                  << std::setw(10) << pos[1]
                  << std::setw(10) << pos[2]
                  << "  "
                  << std::fixed << std::setprecision(3) 
                  << std::setw(10) << dir[0]
                  << std::setw(10) << dir[1]
                  << std::setw(10) << dir[2]
                  << "\n";
         }
 
         // Fill histograms.
 
         if(fRecoHistMap.count(0) == 0)
           fRecoHistMap[0] = RecoHists("Reco");
         const RecoHists& rhists = fRecoHistMap[0];
 
         rhists.fHx->Fill(pos.X());
         rhists.fHy->Fill(pos.Y());
         rhists.fHz->Fill(pos.Z());
         rhists.fHdist->Fill(dpos);
         rhists.fHtheta->Fill(dir.Theta());
         rhists.fHphi->Fill(dir.Phi());
         rhists.fHtheta_xz->Fill(theta_xz);
         rhists.fHtheta_yz->Fill(theta_yz);
       }
     }
   }

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

Referenced by art::EDAnalyzer::doEvent(), and art::OutputModule::doWriteEvent().

   {
     return selectors_;
   }

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::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

Referenced by art::EDAnalyzer::workerType().

   {
     return current_context_.get();
   }

void trkf::SeedAna::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 1060 of file SeedAna_module.cc.

   {
     // Print summary.
 
     mf::LogInfo("SeedAna") 
       << "SeedAna statistics:\n"
       << "  Number of events = " << fNumEvent;
 
     // Fill multiplicity histograms.
 
     for(std::map<int, MCHists>::const_iterator i = fMCHistMap.begin();
         i != fMCHistMap.end(); ++i) {
       const MCHists& mchists = i->second;
       mulcalc(mchists.fHmstartx, mchists.fHmcstartx, mchists.fHmulstartx);
       mulcalc(mchists.fHmstarty, mchists.fHmcstarty, mchists.fHmulstarty);
       mulcalc(mchists.fHmstartz, mchists.fHmcstartz, mchists.fHmulstartz);
       mulcalc(mchists.fHmendx, mchists.fHmcendx, mchists.fHmulendx);
       mulcalc(mchists.fHmendy, mchists.fHmcendy, mchists.fHmulendy);
       mulcalc(mchists.fHmendz, mchists.fHmcendz, mchists.fHmulendz);
       mulcalc(mchists.fHmtheta, mchists.fHmctheta, mchists.fHmultheta);
       mulcalc(mchists.fHmphi, mchists.fHmcphi, mchists.fHmulphi);
       mulcalc(mchists.fHmtheta_xz, mchists.fHmctheta_xz, mchists.fHmultheta_xz);
       mulcalc(mchists.fHmtheta_yz, mchists.fHmctheta_yz, mchists.fHmultheta_yz);
       mulcalc(mchists.fHmmom, mchists.fHmcmom, mchists.fHmulmom);
       mulcalc(mchists.fHmlen, mchists.fHmclen, mchists.fHmullen);
     }
     // Fill efficiency histograms.
 
     for(std::map<int, MCHists>::const_iterator i = fMCHistMap.begin();
         i != fMCHistMap.end(); ++i) {
       const MCHists& mchists = i->second;
       effcalc(mchists.fHgstartx, mchists.fHmcstartx, mchists.fHestartx);
       effcalc(mchists.fHgstarty, mchists.fHmcstarty, mchists.fHestarty);
       effcalc(mchists.fHgstartz, mchists.fHmcstartz, mchists.fHestartz);
       effcalc(mchists.fHgendx, mchists.fHmcendx, mchists.fHeendx);
       effcalc(mchists.fHgendy, mchists.fHmcendy, mchists.fHeendy);
       effcalc(mchists.fHgendz, mchists.fHmcendz, mchists.fHeendz);
       effcalc(mchists.fHgtheta, mchists.fHmctheta, mchists.fHetheta);
       effcalc(mchists.fHgphi, mchists.fHmcphi, mchists.fHephi);
       effcalc(mchists.fHgtheta_xz, mchists.fHmctheta_xz, mchists.fHetheta_xz);
       effcalc(mchists.fHgtheta_yz, mchists.fHmctheta_yz, mchists.fHetheta_yz);
       effcalc(mchists.fHgmom, mchists.fHmcmom, mchists.fHemom);
       effcalc(mchists.fHglen, mchists.fHmclen, mchists.fHelen);
     }
   }

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();
 }

art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const & e ) const

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

Referenced by art::OutputModule::doWriteEvent().

   {
     return selectors_.getOneTriggerResults(e);
   }

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::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

   {
     return false;
   }

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);
   }
 }

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

   {
     return process_name_;
   }

void art::EventObserverBase::registerProducts	(	MasterProductRegistry &	,
		ProductDescriptions &	,
		ModuleDescription const &
	)

inlineinherited

Definition at line 33 of file EventObserverBase.h.

36 {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

   {
     return selector_config_id_;
   }

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);
 }

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), art::RootOutput::RootOutput(), and art::OutputWorker::wantAllEvents().

   {
     return wantAllEvents_;
   }

bool art::EventObserverBase::wantEvent ( Event const & e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), and art::OutputModule::doWriteEvent().

   {
     return selectors_.wantEvent(e);
   }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

     {
       return "WorkerT<EDAnalyzer>";
     }

Member Data Documentation

int trkf::SeedAna::fDump

private

Definition at line 371 of file SeedAna_module.cc.

Referenced by analyze(), and SeedAna().

bool trkf::SeedAna::fIgnoreSign

private

Definition at line 376 of file SeedAna_module.cc.

Referenced by analyze().

double trkf::SeedAna::fMatchColinearity

private

Definition at line 374 of file SeedAna_module.cc.

Referenced by analyze().

double trkf::SeedAna::fMatchDisp

private

Definition at line 375 of file SeedAna_module.cc.

Referenced by analyze().

std::map<int, MCHists> trkf::SeedAna::fMCHistMap

private

Definition at line 380 of file SeedAna_module.cc.

Referenced by analyze(), and endJob().

std::string trkf::SeedAna::fMCTrackModuleLabel

private

Definition at line 370 of file SeedAna_module.cc.

Referenced by analyze(), and SeedAna().

double trkf::SeedAna::fMinMCKE

private

Definition at line 372 of file SeedAna_module.cc.

Referenced by analyze(), and SeedAna().

double trkf::SeedAna::fMinMCLen

private

Definition at line 373 of file SeedAna_module.cc.

Referenced by analyze(), and SeedAna().

int trkf::SeedAna::fNumEvent

private

Definition at line 385 of file SeedAna_module.cc.

Referenced by analyze(), and endJob().

std::map<int, RecoHists> trkf::SeedAna::fRecoHistMap

private

Definition at line 381 of file SeedAna_module.cc.

Referenced by analyze().

std::string trkf::SeedAna::fSeedModuleLabel

private

Definition at line 369 of file SeedAna_module.cc.

Referenced by analyze(), and SeedAna().

The documentation for this class was generated from the following file:

/cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/TrackFinder/SeedAna_module.cc

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation