Inheritance diagram for cheat::RecoCheckAna:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Public Member Functions
	RecoCheckAna (fhicl::ParameterSet const &p)

virtual	~RecoCheckAna ()

virtual void	analyze (art::Event const &e)

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

virtual void	beginRun (art::Run const &r)

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 Member Functions
void	CheckReco (int const &colID, std::vector< art::Ptr< recob::Hit > > const &allhits, std::vector< art::Ptr< recob::Hit > > const &colHits, std::map< std::pair< int, int >, std::pair< double, double > > &g4RecoBaseIDToPurityEfficiency)

void	CheckRecoClusters (art::Event const &evt, std::string const &label, art::Handle< std::vector< recob::Cluster > > const &clscol, std::vector< art::Ptr< recob::Hit > > const &allhits)

void	CheckRecoTracks (art::Event const &evt, std::string const &label, art::Handle< std::vector< recob::Track > > const &tcol, std::vector< art::Ptr< recob::Hit > > const &allhits)

void	CheckRecoShowers (art::Event const &evt, std::string const &label, art::Handle< std::vector< recob::Shower > > const &scol, std::vector< art::Ptr< recob::Hit > > const &allhits)

void	CheckRecoVertices (art::Event const &evt, std::string const &label, art::Handle< std::vector< recob::Vertex > > const &vtxcol, std::vector< art::Ptr< recob::Hit > > const &allhits)

void	CheckRecoEvents (art::Event const &evt, std::string const &label, art::Handle< std::vector< recob::Event > > const &evtcol, std::vector< art::Ptr< recob::Hit > > const &allhits)

void	FillResults (std::vector< art::Ptr< recob::Hit > > const &allhits)

void	FlattenMap (std::map< std::pair< int, int >, std::pair< double, double > > const &g4RecoBaseIDToPurityEfficiency, std::map< int, std::vector< std::pair< int, std::pair< double, double > > > > &g4IDToRecoBasePurityEfficiency, TH1D purity, TH1D efficiency, TH1D purityEfficiency, TH2D purityEfficiency2D)

Private Attributes
art::ServiceHandle< cheat::BackTrackerService >	fBT
	the back tracker service More...

art::ServiceHandle< cheat::ParticleInventoryService >	fPI
	the back tracker service More...

std::string	fHitModuleLabel
	label for module making the hits More...

std::string	fClusterModuleLabel
	label for module making the clusters More...

std::string	fShowerModuleLabel
	label for module making the showers More...

std::string	fTrackModuleLabel
	label for module making the tracks More...

std::string	fVertexModuleLabel
	label for module making the vertices More...

std::string	fEventModuleLabel
	label for module making the events More...

bool	fCheckClusters
	should we check the reconstruction of clusters? More...

bool	fCheckShowers
	should we check the reconstruction of showers? More...

bool	fCheckTracks
	should we check the reconstruction of tracks? More...

bool	fCheckVertices
	should we check the reconstruction of vertices? More...

bool	fCheckEvents
	should we check the reconstruction of events? More...

TH1D *	fClusterPurity
	histogram of cluster purity More...

TH1D *	fClusterEfficiency
	histogram of cluster efficiency More...

TH1D *	fClusterPurityEfficiency
	histogram of cluster efficiency times purity More...

TH2D *	fClusterPurityEfficiency2D
	scatter histogram of cluster purity and efficiency More...

TH1D *	fShowerPurity
	histogram of shower purity More...

TH1D *	fShowerEfficiency
	histogram of shower efficiency More...

TH1D *	fShowerPurityEfficiency
	histogram of shower efficiency times purity More...

TH2D *	fShowerPurityEfficiency2D
	scatter histogram of cluster purity and efficiency More...

TH1D *	fTrackPurity
	histogram of track purity More...

TH1D *	fTrackEfficiency
	histogram of track efficiency More...

TH1D *	fTrackPurityEfficiency
	histogram of track efficiency times purity More...

TH2D *	fTrackPurityEfficiency2D
	scatter histogram of cluster purity and efficiency More...

TH1D *	fVertexPurity
	histogram of vertex purity More...

TH1D *	fVertexEfficiency
	histogram of vertex efficiency More...

TH1D *	fVertexPurityEfficiency
	histogram of vertex efficiency times purity More...

TH1D *	fEventPurity
	histogram of event purity More...

TH1D *	fEventEfficiency
	histogram of event efficiency More...

TH1D *	fEventPurityEfficiency
	histogram of event efficiency times purity More...

std::map< std::pair< int, int >, std::pair< double, double > >	fG4ClusterIDToPurityEfficiency

std::map< std::pair< int, int >, std::pair< double, double > >	fG4ShowerIDToPurityEfficiency

std::map< std::pair< int, int >, std::pair< double, double > >	fG4TrackIDToPurityEfficiency

TTree *	fTree
	TTree to save efficiencies. More...

int	frun
	run number More...

int	fevent
	event number More...

int	ftrackid
	geant track ID More...

int	fpdg
	particle pdg code More...

double	fpmom
	particle momentum More...

double	fhiteff
	hitfinder efficiency for this particle More...

int	fnclu
	number of clusters for this particle More...

std::vector< double >	fclueff
	cluster efficiencies More...

std::vector< double >	fclupur
	cluster purities More...

std::vector< int >	fcluid
	cluster IDs More...

int	fnshw
	number of showers for this particle More...

std::vector< double >	fshweff
	shower efficiencies More...

std::vector< double >	fshwpur
	shower purities More...

std::vector< int >	fshwid
	shower IDs More...

int	fntrk
	number of tracks for this particle More...

std::vector< double >	ftrkeff
	track efficiencies More...

std::vector< double >	ftrkpur
	track purities More...

std::vector< int >	ftrkid
	track IDs More...

Detailed Description

Definition at line 42 of file RecoCheckAna_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

cheat::RecoCheckAna::RecoCheckAna ( fhicl::ParameterSet const & p )

explicit

Definition at line 156 of file RecoCheckAna_module.cc.

References reconfigure().

   : EDAnalyzer(p)
 {
   this->reconfigure(p);
 }

cheat::RecoCheckAna::~RecoCheckAna ( )

virtual

Definition at line 163 of file RecoCheckAna_module.cc.

 {
   // Clean up dynamic memory and other resources here.
 }

Member Function Documentation

void cheat::RecoCheckAna::analyze ( art::Event const & e )

virtual

Implements art::EDAnalyzer.

Definition at line 169 of file RecoCheckAna_module.cc.

References CheckRecoClusters(), CheckRecoEvents(), CheckRecoShowers(), CheckRecoTracks(), CheckRecoVertices(), art::EventID::event(), art::Handle< T >::failedToGet(), fCheckClusters, fCheckEvents, fCheckShowers, fCheckTracks, fCheckVertices, fClusterModuleLabel, fevent, fEventModuleLabel, fHitModuleLabel, art::fill_ptr_vector(), FillResults(), frun, fShowerModuleLabel, fTrackModuleLabel, fVertexModuleLabel, art::DataViewImpl::getByLabel(), art::Event::id(), art::Event::isRealData(), and art::Event::run().

 {
   // check that this is MC, stop if it isn't
   if(e.isRealData()){
     mf::LogWarning("RecoVetter") << "attempting to run MC truth check on "
                                  << "real data, bail";
     return;
   }
 
   // get all hits in the event to figure out how many there are
   art::Handle< std::vector<recob::Hit> > hithdl;
   e.getByLabel(fHitModuleLabel, hithdl);
   std::vector< art::Ptr<recob::Hit> > allhits;
   art::fill_ptr_vector(allhits, hithdl);
 
   // define variables to hold the reconstructed objects
   art::Handle< std::vector<recob::Cluster> > clscol;
   art::Handle< std::vector<recob::Track> >   trkcol;
   art::Handle< std::vector<recob::Shower> >  shwcol;
   art::Handle< std::vector<recob::Vertex> >  vtxcol;
   art::Handle< std::vector<recob::Event> >   evtcol;
 
   if(fCheckClusters){
     e.getByLabel(fClusterModuleLabel, clscol);
     if( !clscol.failedToGet() ) this->CheckRecoClusters(e, fClusterModuleLabel, clscol, allhits);
   }
   if(fCheckTracks){
     e.getByLabel(fTrackModuleLabel, trkcol);
     if( !trkcol.failedToGet() ) this->CheckRecoTracks(e, fTrackModuleLabel, trkcol, allhits);
   }
   if(fCheckShowers){
     e.getByLabel(fShowerModuleLabel, shwcol);
     if( !shwcol.failedToGet() ) this->CheckRecoShowers(e, fShowerModuleLabel, shwcol, allhits);
   }
   if(fCheckVertices){
     e.getByLabel(fVertexModuleLabel, vtxcol);
     if( !vtxcol.failedToGet() ) this->CheckRecoVertices(e, fVertexModuleLabel, vtxcol, allhits);
   }
   if(fCheckEvents){
     e.getByLabel(fEventModuleLabel, evtcol);
     if( !evtcol.failedToGet() ) this->CheckRecoEvents(e, fEventModuleLabel, evtcol, allhits);
   }
 
   frun = e.run();
   fevent = e.id().event();
 
   this->FillResults(allhits);
 
   return;
  
 }

void cheat::RecoCheckAna::beginRun ( art::Run const & r )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 239 of file RecoCheckAna_module.cc.

References fCheckClusters, fCheckEvents, fCheckShowers, fCheckTracks, fCheckVertices, fclueff, fcluid, fclupur, fClusterEfficiency, fClusterPurity, fClusterPurityEfficiency, fClusterPurityEfficiency2D, fevent, fEventEfficiency, fEventPurity, fEventPurityEfficiency, fhiteff, fnclu, fnshw, fntrk, fpdg, fpmom, frun, fShowerEfficiency, fShowerPurity, fShowerPurityEfficiency, fShowerPurityEfficiency2D, fshweff, fshwid, fshwpur, fTrackEfficiency, ftrackid, fTrackPurity, fTrackPurityEfficiency, fTrackPurityEfficiency2D, fTree, ftrkeff, ftrkid, ftrkpur, fVertexEfficiency, fVertexPurity, fVertexPurityEfficiency, and art::TFileDirectory::make().

 {
   art::ServiceHandle<art::TFileService> tfs;
 
   if(fCheckEvents){
     fEventPurity             = tfs->make<TH1D>("eventPurity",       ";Purity;Events",       100, 0., 1.1);
     fEventEfficiency         = tfs->make<TH1D>("eventEfficiency",   ";Efficiency;Events",   100, 0., 1.1);
     fEventPurityEfficiency   = tfs->make<TH1D>("eventPurityEfficiency", ";purityEfficiency;Events", 110, 0., 1.1);
   }
   if(fCheckVertices){
     fVertexPurity            = tfs->make<TH1D>("vertexPurity",      ";Purity;Vertices",     100, 0., 1.1);
     fVertexEfficiency        = tfs->make<TH1D>("vertexEfficiency",  ";Efficiency;Vertices", 100, 0., 1.1);
     fVertexPurityEfficiency  = tfs->make<TH1D>("vertexPurityEfficiency", ";purityEfficiency;Vertex", 110, 0., 1.1);
   }
   if(fCheckTracks){
     fTrackPurity             = tfs->make<TH1D>("trackPurity",       ";Purity;Tracks",       100, 0., 1.1);
     fTrackEfficiency         = tfs->make<TH1D>("trackEfficiency",   ";Efficiency;Tracks",   100, 0., 1.1);
     fTrackPurityEfficiency   = tfs->make<TH1D>("trackPurityEfficiency", ";purityEfficiency;Tracks", 110, 0., 1.1);
     fTrackPurityEfficiency2D = tfs->make<TH2D>("trackPurityEfficiency2D", ";purity;efficiency", 110, 0., 1.1, 110, 0., 1.1);
   }
   if(fCheckShowers){
     fShowerPurity            = tfs->make<TH1D>("showerPurity",      ";Purity;Showers",      100, 0., 1.1);
     fShowerEfficiency        = tfs->make<TH1D>("showerEfficiency",  ";Efficiency;Showers",  100, 0., 1.1);
     fShowerPurityEfficiency  = tfs->make<TH1D>("showerPurityEfficiency", ";purityEfficiency;Showers", 110, 0., 1.1);
     fShowerPurityEfficiency2D= tfs->make<TH2D>("showerPurityEfficiency2D", ";purity;efficiency", 110, 0., 1.1, 110, 0., 1.1);
   }
   if(fCheckClusters){
     fClusterPurity              = tfs->make<TH1D>("clusterPurity",     ";Purity;Clusters",     110, 0., 1.1);
     fClusterEfficiency          = tfs->make<TH1D>("clusterEfficiency", ";Efficiency;Clusters", 110, 0., 1.1);
     fClusterPurityEfficiency    = tfs->make<TH1D>("clusterPurityEfficiency", ";purityEfficiency;Clusters", 110, 0., 1.1);
     fClusterPurityEfficiency2D  = tfs->make<TH2D>("clusterPurityEfficiency2D", ";purity;efficiency", 110, 0., 1.1, 110, 0., 1.1);
   }
 
   fTree = tfs->make<TTree>("cheatertree","cheater tree");
   fTree->Branch("run",     &frun,     "run/I"    );
   fTree->Branch("event",   &fevent,   "event/I"  );
   fTree->Branch("trackid", &ftrackid, "trackid/I");
   fTree->Branch("pdg",     &fpdg,     "pdg/I"    );
   fTree->Branch("pmom",    &fpmom,    "pmom/D"   );
   fTree->Branch("hiteff",  &fhiteff,  "hiteff/D" );
   fTree->Branch("nclu",    &fnclu,    "nclu/I"   );
   fTree->Branch("clueff",  &fclueff              );
   fTree->Branch("clupur",  &fclupur              );
   fTree->Branch("cluid",   &fcluid               );
   fTree->Branch("nshw",    &fnshw,    "nshw/I"   );
   fTree->Branch("shweff",  &fshweff              );
   fTree->Branch("shwpur",  &fshwpur              );
   fTree->Branch("shwid",   &fshwid               );
   fTree->Branch("ntrk",    &fntrk,    "ntrk/I"   );
   fTree->Branch("trkeff",  &ftrkeff              );
   fTree->Branch("trkpur",  &ftrkpur              );
   fTree->Branch("trkid",   &ftrkid               );
 
   return;
 }

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

void cheat::RecoCheckAna::CheckReco	(	int const &	colID,
		std::vector< art::Ptr< recob::Hit > > const &	allhits,
		std::vector< art::Ptr< recob::Hit > > const &	colHits,
		std::map< std::pair< int, int >, std::pair< double, double > > &	g4RecoBaseIDToPurityEfficiency
	)

private

Definition at line 298 of file RecoCheckAna_module.cc.

References fBT, cheat::BackTrackerService::GetSetOfTrackIds(), cheat::BackTrackerService::HitCollectionEfficiency(), and cheat::BackTrackerService::HitCollectionPurity().

Referenced by CheckRecoClusters(), CheckRecoShowers(), and CheckRecoTracks().

 {
 
   // grab the set of track IDs for these hits
   std::set<int> trackIDs = fBT->GetSetOfTrackIds(colHits);
 
   geo::View_t view = colHits[0]->View();
     
   std::set<int>::iterator itr = trackIDs.begin();
   while( itr != trackIDs.end() ){
 
     //std::cout << "*itr: " << *itr << std::endl;
 
     std::set<int> id;
     id.insert(*itr);
 
     // use the cheat::BackTrackerService to find purity and efficiency for these hits
     double purity     = fBT->HitCollectionPurity(id, colHits);
     double efficiency = fBT->HitCollectionEfficiency(id, colHits, allhits, view);
     
     // make the purity and efficiency pair
     std::pair<double, double> pe(purity, efficiency);
     
     // make the pair of the RecoBase object id to the pair of purity/efficiency
     std::pair<int, int> g4reco(*itr, colID);
     
     // insert idpe into the map
     g4RecoBaseIDToPurityEfficiency[g4reco] = pe;
 
     itr++;
 
   } // end loop over eveIDs
 
   return;
 }

void cheat::RecoCheckAna::CheckRecoClusters	(	art::Event const &	evt,
		std::string const &	label,
		art::Handle< std::vector< recob::Cluster > > const &	clscol,
		std::vector< art::Ptr< recob::Hit > > const &	allhits
	)

private

Definition at line 338 of file RecoCheckAna_module.cc.

References CheckReco(), fG4ClusterIDToPurityEfficiency, and hits().

Referenced by analyze().

 {
 
   art::FindManyP<recob::Hit> fmh(clscol, evt, label);
 
   for(size_t c = 0; c < clscol->size(); ++c){
 
     // get the hits associated with this event
     std::vector< art::Ptr< recob::Hit > > hits = fmh.at(c);
 
     this->CheckReco(clscol->at(c).ID(), allhits, hits, fG4ClusterIDToPurityEfficiency);
       
   }// end loop over clusters
 
   return;
 }

void cheat::RecoCheckAna::CheckRecoEvents	(	art::Event const &	evt,
		std::string const &	label,
		art::Handle< std::vector< recob::Event > > const &	evtcol,
		std::vector< art::Ptr< recob::Hit > > const &	allhits
	)

private

Todo:: need to divy it up in the case where there is more than 1 true interaction in a spill

Definition at line 461 of file RecoCheckAna_module.cc.

References d, simb::MCParticle::Daughter(), fBT, fEventEfficiency, fEventPurity, fEventPurityEfficiency, fPI, cheat::BackTrackerService::HitCollectionEfficiency(), cheat::BackTrackerService::HitCollectionPurity(), hits(), geo::k3D, simb::MCParticle::NumberDaughters(), part, and cheat::ParticleInventoryService::ParticleList().

Referenced by analyze().

 {
   const sim::ParticleList& plist = fPI->ParticleList();
 
   // loop over all primaries in the plist and grab them and their daughters to put into 
   // the set of track ids to pass on to the back tracker
   std::set<int> ids;
   for (const auto& PartPair: plist) {
     auto trackID = PartPair.first;
     if (!plist.IsPrimary(trackID)) continue;
     const simb::MCParticle& part = *(PartPair.second);
     ids.insert(trackID);
     for(int d = 0; d < part.NumberDaughters(); ++d)
       ids.insert(part.Daughter(d));
   }// end loop over primaries
 
 
 
   art::FindManyP<recob::Hit> fmh(evtcol, evt, label);
 
   for(size_t ev = 0; ev < evtcol->size(); ++ev){
 
     // get the hits associated with this event
     std::vector< art::Ptr< recob::Hit > > hits = fmh.at(ev);
 
     // use the cheat::BackTrackerService to find purity and efficiency for these hits
     double purity     = fBT->HitCollectionPurity(ids, hits);
     double efficiency = fBT->HitCollectionEfficiency(ids, hits, allhits, geo::k3D);
 
     fEventPurity    ->Fill(purity);
     fEventEfficiency->Fill(efficiency);
     fEventPurityEfficiency->Fill(purity*efficiency);
 
   }// end loop over events
 
   return;
 }

void cheat::RecoCheckAna::CheckRecoShowers	(	art::Event const &	evt,
		std::string const &	label,
		art::Handle< std::vector< recob::Shower > > const &	scol,
		std::vector< art::Ptr< recob::Hit > > const &	allhits
	)

private

Definition at line 380 of file RecoCheckAna_module.cc.

References CheckReco(), fG4ShowerIDToPurityEfficiency, and hits().

Referenced by analyze().

 {
 
   art::FindManyP<recob::Hit> fmh(scol, evt, label);
 
   for(size_t p = 0; p < scol->size(); ++p){
 
     // get the hits associated with this event
     std::vector< art::Ptr< recob::Hit > > hits = fmh.at(p);
 
     this->CheckReco(scol->at(p).ID(), allhits, hits, fG4ShowerIDToPurityEfficiency);
 
   }// end loop over events
 
   return;
 }

void cheat::RecoCheckAna::CheckRecoTracks	(	art::Event const &	evt,
		std::string const &	label,
		art::Handle< std::vector< recob::Track > > const &	tcol,
		std::vector< art::Ptr< recob::Hit > > const &	allhits
	)

private

Definition at line 359 of file RecoCheckAna_module.cc.

References CheckReco(), fG4TrackIDToPurityEfficiency, and hits().

Referenced by analyze().

 {
 
   art::FindManyP<recob::Hit> fmh(tcol, evt, label);
 
   for(size_t p = 0; p < tcol->size(); ++p){
 
     // get the hits associated with this event
     std::vector< art::Ptr< recob::Hit > > hits = fmh.at(p);
 
     this->CheckReco(tcol->at(p).ID(), allhits, hits, fG4TrackIDToPurityEfficiency);
 
   }// end loop over tracks
 
   return;
 }

void cheat::RecoCheckAna::CheckRecoVertices	(	art::Event const &	evt,
		std::string const &	label,
		art::Handle< std::vector< recob::Vertex > > const &	vtxcol,
		std::vector< art::Ptr< recob::Hit > > const &	allhits
	)

private

Definition at line 403 of file RecoCheckAna_module.cc.

References d, simb::MCParticle::Daughter(), fBT, fPI, fVertexEfficiency, fVertexPurity, fVertexPurityEfficiency, cheat::BackTrackerService::HitCollectionEfficiency(), cheat::BackTrackerService::HitCollectionPurity(), hits(), geo::k3D, simb::MCParticle::NumberDaughters(), part, and cheat::ParticleInventoryService::ParticleList().

Referenced by analyze().

 {
   const sim::ParticleList& plist = fPI->ParticleList();
 
   std::vector< std::set<int> > ids(1);
   // loop over all primary particles and put their ids into the first set of the 
   // vector.  add another set for each primary particle that also has daughters
   // and put those daughters into the new set
   // PartPair is a (track ID, particle pointer) pair
   for (const auto& PartPair: plist) {
     auto trackID = PartPair.first;
     if (!plist.IsPrimary(trackID)) continue;
     const simb::MCParticle& part = *(PartPair.second);
     ids[0].insert(trackID);
     if(part.NumberDaughters() > 0){
       std::set<int> dv;
       for(int d = 0; d < part.NumberDaughters(); ++d)
         dv.insert(part.Daughter(d));
       ids.push_back(std::move(dv));
     }//end if this primary particle has daughters
   }// end loop over primaries
 
   art::FindManyP<recob::Hit> fmh(vtxcol, evt, label);
 
   for(size_t v = 0; v < vtxcol->size(); ++v){
 
     // get the hits associated with this event
     std::vector< art::Ptr< recob::Hit > > hits = fmh.at(v);
 
     double maxPurity     = -1.;
     double maxEfficiency = -1.;
 
     for(size_t tv = 0; tv < ids.size(); ++tv){
 
       // use the cheat::BackTrackerService to find purity and efficiency for these hits
       double purity     = fBT->HitCollectionPurity(ids[tv], hits);
       double efficiency = fBT->HitCollectionEfficiency(ids[tv], hits, allhits, geo::k3D);
 
       if(purity     > maxPurity    ) maxPurity     = purity;
       if(efficiency > maxEfficiency) maxEfficiency = efficiency;
     }
 
     fVertexPurity    ->Fill(maxPurity);
     fVertexEfficiency->Fill(maxEfficiency);
     fVertexPurityEfficiency->Fill(maxPurity*maxEfficiency);
 
   }// end loop over vertices
 
   return;
 }

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 cheat::RecoCheckAna::FillResults ( std::vector< art::Ptr< recob::Hit > > const & allhits )

private

Definition at line 560 of file RecoCheckAna_module.cc.

References DEFINE_ART_MODULE, e, energy, fBT, fCheckClusters, fCheckShowers, fCheckTracks, fclueff, fcluid, fclupur, fClusterEfficiency, fClusterPurity, fClusterPurityEfficiency, fClusterPurityEfficiency2D, fG4ClusterIDToPurityEfficiency, fG4ShowerIDToPurityEfficiency, fG4TrackIDToPurityEfficiency, fhiteff, FlattenMap(), fnclu, fnshw, fntrk, fpdg, fPI, fpmom, fShowerEfficiency, fShowerPurity, fShowerPurityEfficiency, fShowerPurityEfficiency2D, fshweff, fshwid, fshwpur, fTrackEfficiency, ftrackid, fTrackPurity, fTrackPurityEfficiency, fTrackPurityEfficiency2D, fTree, ftrkeff, ftrkid, ftrkpur, cheat::BackTrackerService::GetSetOfTrackIds(), cheat::BackTrackerService::HitToTrackIDEs(), simb::MCParticle::P(), part, simb::MCParticle::PdgCode(), s, cheat::ParticleInventoryService::TrackIdToParticle_P(), and cheat::BackTrackerService::TrackIdToSimIDEs_Ps().

Referenced by analyze().

 {
   // map the g4 track id to energy deposited in a hit
   std::map<int, double> g4IDToHitEnergy;
   for(size_t h = 0; h < allhits.size(); ++h){
     const std::vector<sim::TrackIDE> hitTrackIDs = fBT->HitToTrackIDEs(allhits[h]);
     for(size_t e = 0; e < hitTrackIDs.size(); ++e){
       g4IDToHitEnergy[hitTrackIDs[e].trackID] += hitTrackIDs[e].energy;
     }
   } // end loop over hits to fill map
 
   // flatten the G4RecoBaseIDToPurityEfficiency maps to have just the g4ID as the key and the 
   // rest of the information in vector form
   std::map<int, std::vector< std::pair<int, std::pair<double, double> > > > g4IDToClusterPurityEfficiency;
   std::map<int, std::vector< std::pair<int, std::pair<double, double> > > > g4IDToShowerPurityEfficiency;
   std::map<int, std::vector< std::pair<int, std::pair<double, double> > > > g4IDToTrackPurityEfficiency;
   std::map<int, std::vector< std::pair<int, std::pair<double, double> > > >::iterator g4peItr;
 
   if(fCheckClusters) this->FlattenMap(fG4ClusterIDToPurityEfficiency, g4IDToClusterPurityEfficiency, fClusterPurity, fClusterEfficiency, fClusterPurityEfficiency, fClusterPurityEfficiency2D);
   if(fCheckShowers)  this->FlattenMap(fG4ShowerIDToPurityEfficiency,  g4IDToShowerPurityEfficiency,  fShowerPurity,  fShowerEfficiency, fShowerPurityEfficiency, fShowerPurityEfficiency2D);
   if(fCheckTracks)   this->FlattenMap(fG4TrackIDToPurityEfficiency,   g4IDToTrackPurityEfficiency,   fTrackPurity,   fTrackEfficiency, fTrackPurityEfficiency, fTrackPurityEfficiency2D);
 
   // fill the tree vectors
   // get all the eveIDs from this event
   std::set<int> trackIDs = fBT->GetSetOfTrackIds();
   std::set<int>::const_iterator trackItr = trackIDs.begin();
 
   // loop over them
   while( trackItr != trackIDs.end() ){
 
     const simb::MCParticle* part = fPI->TrackIdToParticle_P(*trackItr);
 
     ftrackid = std::abs(*trackItr);
     fpdg     = part->PdgCode();
     fpmom    = part->P();
 
     // figure out how much of the energy deposited from this particle is stored in hits
     std::vector<const sim::IDE*> ides = fBT->TrackIdToSimIDEs_Ps(*trackItr);
     double totalDep = 0.;
     for(size_t i = 0; i < ides.size(); ++i) totalDep += ides[i]->energy;
 
     if(totalDep > 0.)
       fhiteff  = g4IDToHitEnergy[*trackItr]/totalDep;
 
     std::vector< std::pair<int, std::pair<double, double> > > clVec;
     std::vector< std::pair<int, std::pair<double, double> > > shVec;
     std::vector< std::pair<int, std::pair<double, double> > > trVec;
    
 
     if( g4IDToClusterPurityEfficiency.find(*trackItr) != g4IDToClusterPurityEfficiency.end() )
       clVec = g4IDToClusterPurityEfficiency.find(*trackItr)->second;
 
     if( g4IDToShowerPurityEfficiency.find(*trackItr) != g4IDToShowerPurityEfficiency.end() )
       shVec = g4IDToShowerPurityEfficiency.find(*trackItr)->second;
 
     if( g4IDToTrackPurityEfficiency.find(*trackItr) != g4IDToTrackPurityEfficiency.end() )
       trVec = g4IDToTrackPurityEfficiency.find(*trackItr)->second;
 
     fnclu = clVec.size();
     fnshw = shVec.size();
     fntrk = trVec.size();
 
     for(size_t c = 0; c < clVec.size(); ++c){
       fcluid .push_back(clVec[c].first);
       fclupur.push_back(clVec[c].second.first);
       fclueff.push_back(clVec[c].second.second);
     }      
 
     for(size_t s = 0; s < shVec.size(); ++s){
       fshwid .push_back(shVec[s].first);
       fshwpur.push_back(shVec[s].second.first);
       fshweff.push_back(shVec[s].second.second);
     }      
 
     for(size_t t = 0; t < trVec.size(); ++t){
       ftrkid .push_back(trVec[t].first);
       ftrkpur.push_back(trVec[t].second.first);
       ftrkeff.push_back(trVec[t].second.second);
     }      
 
     fTree->Fill();
 
     trackItr++;
   }
 
   // clean up for the next event
 
   // clear the maps of G4 track id to efficiency and purity for 
   // various RecoBase objects
   fG4ClusterIDToPurityEfficiency.clear();
   fG4ShowerIDToPurityEfficiency .clear();
   fG4TrackIDToPurityEfficiency  .clear(); 
 
   // clear the vectors hooked up to the tree
   fclueff.clear();
   fclupur.clear();
   fcluid .clear();
   ftrkeff.clear();
   ftrkpur.clear();
   ftrkid .clear();
   fshweff.clear();
   fshwpur.clear();
   fshwid .clear();
 
   return;
 }

void cheat::RecoCheckAna::FlattenMap	(	std::map< std::pair< int, int >, std::pair< double, double > > const &	g4RecoBaseIDToPurityEfficiency,
		std::map< int, std::vector< std::pair< int, std::pair< double, double > > > > &	g4IDToRecoBasePurityEfficiency,
		TH1D *	purity,
		TH1D *	efficiency,
		TH1D *	purityEfficiency,
		TH2D *	purityEfficiency2D
	)

private

Definition at line 503 of file RecoCheckAna_module.cc.

Referenced by FillResults().

 {
 
   std::map<std::pair<int, int>, std::pair<double, double> >::const_iterator rbItr = g4RecoBaseIDToPurityEfficiency.begin();
 
   // map of key cluster ID to pair of purity, efficiency
   std::map<int, std::pair<double, double> > recoBIDToPurityEfficiency;
   std::map<int, std::pair<double, double> >::iterator rbpeItr;
 
   while( rbItr != g4RecoBaseIDToPurityEfficiency.end() ){
     
 
     // trackID, cluster ID
     std::pair<int, int>       g4cl = rbItr->first;
     // purity, efficiency
     std::pair<double, double> pe   = rbItr->second;
 
     // add the efficiency and purity values for clusters corresponding
     // to the current g4 id to the map
     // pair of cluster id, pair of purity, efficiency
     std::pair<int, std::pair<double, double> > clpe(g4cl.second, pe);
     // g4IDToRecoBasePurityEfficiency is a map with key of trackID of a vector of clusterIDs of pairs of purity and efficiency
     g4IDToRecoBasePurityEfficiency[g4cl.first].push_back(clpe);
 
     // now find the maximum purity to determine the purity and efficiency
     // for this RecoBase object
     rbpeItr = recoBIDToPurityEfficiency.find(g4cl.second);
     if( rbpeItr != recoBIDToPurityEfficiency.end() ){
       std::pair<double, double> curpe = rbpeItr->second;
       if(pe.first > curpe.first) recoBIDToPurityEfficiency[g4cl.second] = pe;
     }
     else
       recoBIDToPurityEfficiency[g4cl.second] = pe;
 
     rbItr++;
   }
 
   rbpeItr = recoBIDToPurityEfficiency.begin();
 
   // now fill the histograms, 
   while(rbpeItr != recoBIDToPurityEfficiency.end() ){
     purity    ->Fill(rbpeItr->second.first);
     efficiency->Fill(rbpeItr->second.second);
     purityEfficiency->Fill(rbpeItr->second.first*rbpeItr->second.second);
     purityEfficiency2D->Fill(rbpeItr->second.first,rbpeItr->second.second);
     rbpeItr++;
   }
 
   return;
 }

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 cheat::RecoCheckAna::reconfigure ( fhicl::ParameterSet const & p )

virtual

Definition at line 222 of file RecoCheckAna_module.cc.

References fCheckClusters, fCheckEvents, fCheckShowers, fCheckTracks, fCheckVertices, fClusterModuleLabel, fEventModuleLabel, fHitModuleLabel, fShowerModuleLabel, fTrackModuleLabel, fVertexModuleLabel, and fhicl::ParameterSet::get().

Referenced by RecoCheckAna().

 {
   fHitModuleLabel     = p.get< std::string >("HitModuleLabel");
   fClusterModuleLabel = p.get< std::string >("ClusterModuleLabel");
   fShowerModuleLabel  = p.get< std::string >("ShowerModuleLabel"  );
   fTrackModuleLabel   = p.get< std::string >("TrackModuleLabel" );
   fVertexModuleLabel  = p.get< std::string >("VertexModuleLabel" );
   fEventModuleLabel   = p.get< std::string >("EventModuleLabel"  );
 
   fCheckClusters      = p.get< bool        >("CheckClusters");
   fCheckShowers       = p.get< bool        >("CheckShowers" );
   fCheckTracks        = p.get< bool        >("CheckTracks"  );
   fCheckVertices      = p.get< bool        >("CheckVertices");
   fCheckEvents        = p.get< bool        >("CheckEvents"  );
 }

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

art::ServiceHandle<cheat::BackTrackerService> cheat::RecoCheckAna::fBT

private

the back tracker service

Definition at line 91 of file RecoCheckAna_module.cc.

Referenced by CheckReco(), CheckRecoEvents(), CheckRecoVertices(), and FillResults().

bool cheat::RecoCheckAna::fCheckClusters

private

should we check the reconstruction of clusters?

Definition at line 101 of file RecoCheckAna_module.cc.

Referenced by analyze(), beginRun(), FillResults(), and reconfigure().

bool cheat::RecoCheckAna::fCheckEvents

private

should we check the reconstruction of events?

Definition at line 105 of file RecoCheckAna_module.cc.

Referenced by analyze(), beginRun(), and reconfigure().

bool cheat::RecoCheckAna::fCheckShowers

private

should we check the reconstruction of showers?

Definition at line 102 of file RecoCheckAna_module.cc.

Referenced by analyze(), beginRun(), FillResults(), and reconfigure().

bool cheat::RecoCheckAna::fCheckTracks

private

should we check the reconstruction of tracks?

Definition at line 103 of file RecoCheckAna_module.cc.

Referenced by analyze(), beginRun(), FillResults(), and reconfigure().

bool cheat::RecoCheckAna::fCheckVertices

private

should we check the reconstruction of vertices?

Definition at line 104 of file RecoCheckAna_module.cc.

Referenced by analyze(), beginRun(), and reconfigure().

std::vector<double> cheat::RecoCheckAna::fclueff

private

cluster efficiencies

Definition at line 141 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<int> cheat::RecoCheckAna::fcluid

private

cluster IDs

Definition at line 143 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<double> cheat::RecoCheckAna::fclupur

private

cluster purities

Definition at line 142 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fClusterEfficiency

private

histogram of cluster efficiency

Definition at line 108 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::string cheat::RecoCheckAna::fClusterModuleLabel

private

label for module making the clusters

Definition at line 95 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* cheat::RecoCheckAna::fClusterPurity

private

histogram of cluster purity

Definition at line 107 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fClusterPurityEfficiency

private

histogram of cluster efficiency times purity

Definition at line 109 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH2D* cheat::RecoCheckAna::fClusterPurityEfficiency2D

private

scatter histogram of cluster purity and efficiency

Definition at line 110 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fevent

private

event number

Definition at line 135 of file RecoCheckAna_module.cc.

Referenced by analyze(), and beginRun().

TH1D* cheat::RecoCheckAna::fEventEfficiency

private

histogram of event efficiency

Definition at line 123 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoEvents().

std::string cheat::RecoCheckAna::fEventModuleLabel

private

label for module making the events

Definition at line 99 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* cheat::RecoCheckAna::fEventPurity

private

histogram of event purity

Definition at line 122 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoEvents().

TH1D* cheat::RecoCheckAna::fEventPurityEfficiency

private

histogram of event efficiency times purity

Definition at line 124 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoEvents().

std::map< std::pair<int, int>, std::pair<double, double> > cheat::RecoCheckAna::fG4ClusterIDToPurityEfficiency

private

Definition at line 129 of file RecoCheckAna_module.cc.

Referenced by CheckRecoClusters(), and FillResults().

std::map< std::pair<int, int>, std::pair<double, double> > cheat::RecoCheckAna::fG4ShowerIDToPurityEfficiency

private

Definition at line 130 of file RecoCheckAna_module.cc.

Referenced by CheckRecoShowers(), and FillResults().

std::map< std::pair<int, int>, std::pair<double, double> > cheat::RecoCheckAna::fG4TrackIDToPurityEfficiency

private

Definition at line 131 of file RecoCheckAna_module.cc.

Referenced by CheckRecoTracks(), and FillResults().

double cheat::RecoCheckAna::fhiteff

private

hitfinder efficiency for this particle

Definition at line 139 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::string cheat::RecoCheckAna::fHitModuleLabel

private

label for module making the hits

Definition at line 94 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

int cheat::RecoCheckAna::fnclu

private

number of clusters for this particle

Definition at line 140 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fnshw

private

number of showers for this particle

Definition at line 144 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fntrk

private

number of tracks for this particle

Definition at line 148 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fpdg

private

particle pdg code

Definition at line 137 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

art::ServiceHandle<cheat::ParticleInventoryService> cheat::RecoCheckAna::fPI

private

the back tracker service

Definition at line 92 of file RecoCheckAna_module.cc.

Referenced by CheckRecoEvents(), CheckRecoVertices(), and FillResults().

double cheat::RecoCheckAna::fpmom

private

particle momentum

Definition at line 138 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::frun

private

run number

Definition at line 134 of file RecoCheckAna_module.cc.

Referenced by analyze(), and beginRun().

TH1D* cheat::RecoCheckAna::fShowerEfficiency

private

histogram of shower efficiency

Definition at line 112 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::string cheat::RecoCheckAna::fShowerModuleLabel

private

label for module making the showers

Definition at line 96 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* cheat::RecoCheckAna::fShowerPurity

private

histogram of shower purity

Definition at line 111 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fShowerPurityEfficiency

private

histogram of shower efficiency times purity

Definition at line 113 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH2D* cheat::RecoCheckAna::fShowerPurityEfficiency2D

private

scatter histogram of cluster purity and efficiency

Definition at line 114 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<double> cheat::RecoCheckAna::fshweff

private

shower efficiencies

Definition at line 145 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<int> cheat::RecoCheckAna::fshwid

private

shower IDs

Definition at line 147 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<double> cheat::RecoCheckAna::fshwpur

private

shower purities

Definition at line 146 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fTrackEfficiency

private

histogram of track efficiency

Definition at line 116 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::ftrackid

private

geant track ID

Definition at line 136 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::string cheat::RecoCheckAna::fTrackModuleLabel

private

label for module making the tracks

Definition at line 97 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* cheat::RecoCheckAna::fTrackPurity

private

histogram of track purity

Definition at line 115 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fTrackPurityEfficiency

private

histogram of track efficiency times purity

Definition at line 117 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH2D* cheat::RecoCheckAna::fTrackPurityEfficiency2D

private

scatter histogram of cluster purity and efficiency

Definition at line 118 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TTree* cheat::RecoCheckAna::fTree

private

TTree to save efficiencies.

Definition at line 133 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<double> cheat::RecoCheckAna::ftrkeff

private

track efficiencies

Definition at line 149 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<int> cheat::RecoCheckAna::ftrkid

private

track IDs

Definition at line 151 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

std::vector<double> cheat::RecoCheckAna::ftrkpur

private

track purities

Definition at line 150 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fVertexEfficiency

private

histogram of vertex efficiency

Definition at line 120 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoVertices().

std::string cheat::RecoCheckAna::fVertexModuleLabel

private

label for module making the vertices

Definition at line 98 of file RecoCheckAna_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* cheat::RecoCheckAna::fVertexPurity

private

histogram of vertex purity

Definition at line 119 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoVertices().

TH1D* cheat::RecoCheckAna::fVertexPurityEfficiency

private

histogram of vertex efficiency times purity

Definition at line 121 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoVertices().

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

larsim/v06_53_00/source/larsim/MCCheater/RecoCheckAna_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation