Inheritance diagram for cheat::RecoCheckAna:

Public Types
using	ModuleType = EDAnalyzer

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

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)

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
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

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	analyze (art::Event const &e) override

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

void	CheckReco (detinfo::DetectorClocksData const &clockData, 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 (detinfo::DetectorClocksData const &clockData, 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 const >	fBT
	the back tracker service More...

art::ServiceHandle< cheat::ParticleInventoryService const >	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 50 of file RecoCheckAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 162 of file RecoCheckAna_module.cc.

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

   : EDAnalyzer(p)
   , 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")}
 {}

Member Function Documentation

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

overrideprivatevirtual

Implements art::EDAnalyzer.

Definition at line 178 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::ProductRetriever::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();
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
   this->FillResults(clockData, allhits);
 
   return;
 }

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

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 231 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, and fVertexPurityEfficiency.

 {
   art::ServiceHandle<art::TFileService const> 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;
 }

void cheat::RecoCheckAna::CheckReco	(	detinfo::DetectorClocksData const &	clockData,
		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(clockData, colHits);
 
   geo::View_t view = colHits[0]->View();
 
   std::set<int>::iterator itr = trackIDs.begin();
   while (itr != trackIDs.end()) {
 
     std::set<int> id;
     id.insert(*itr);
 
     // use the cheat::BackTrackerService to find purity and efficiency for these
     // hits
     double purity = fBT->HitCollectionPurity(clockData, id, colHits);
     double efficiency = fBT->HitCollectionEfficiency(clockData, 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 339 of file RecoCheckAna_module.cc.

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

Referenced by analyze().

 {
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
   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(clockData, 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 465 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);
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
   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(clockData, ids, hits);
     double efficiency = fBT->HitCollectionEfficiency(clockData, 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 381 of file RecoCheckAna_module.cc.

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

Referenced by analyze().

 {
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
   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(clockData, 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 360 of file RecoCheckAna_module.cc.

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

Referenced by analyze().

 {
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
   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(clockData, 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 404 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);
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
 
   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(clockData, ids[tv], hits);
       double efficiency = fBT->HitCollectionEfficiency(clockData, 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 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::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::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 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

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

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

void cheat::RecoCheckAna::FillResults	(	detinfo::DetectorClocksData const &	clockData,
		std::vector< art::Ptr< recob::Hit >> const &	allhits
	)

private

Definition at line 566 of file RecoCheckAna_module.cc.

References util::abs(), 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(), if(), simb::MCParticle::P(), part, simb::MCParticle::PdgCode(), 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(clockData, 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 507 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;
 }

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

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

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

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

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

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

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

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

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }

Member Data Documentation

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

private

the back tracker service

Definition at line 98 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 108 of file RecoCheckAna_module.cc.

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

bool cheat::RecoCheckAna::fCheckEvents

private

should we check the reconstruction of events?

Definition at line 112 of file RecoCheckAna_module.cc.

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

bool cheat::RecoCheckAna::fCheckShowers

private

should we check the reconstruction of showers?

Definition at line 109 of file RecoCheckAna_module.cc.

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

bool cheat::RecoCheckAna::fCheckTracks

private

should we check the reconstruction of tracks?

Definition at line 110 of file RecoCheckAna_module.cc.

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

bool cheat::RecoCheckAna::fCheckVertices

private

should we check the reconstruction of vertices?

Definition at line 111 of file RecoCheckAna_module.cc.

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

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

private

cluster efficiencies

Definition at line 148 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

cluster IDs

Definition at line 150 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

cluster purities

Definition at line 149 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fClusterEfficiency

private

histogram of cluster efficiency

Definition at line 115 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 102 of file RecoCheckAna_module.cc.

Referenced by analyze(), and RecoCheckAna().

TH1D* cheat::RecoCheckAna::fClusterPurity

private

histogram of cluster purity

Definition at line 114 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fClusterPurityEfficiency

private

histogram of cluster efficiency times purity

Definition at line 116 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 117 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fevent

private

event number

Definition at line 142 of file RecoCheckAna_module.cc.

Referenced by analyze(), and beginRun().

TH1D* cheat::RecoCheckAna::fEventEfficiency

private

histogram of event efficiency

Definition at line 130 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 106 of file RecoCheckAna_module.cc.

Referenced by analyze(), and RecoCheckAna().

TH1D* cheat::RecoCheckAna::fEventPurity

private

histogram of event purity

Definition at line 129 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoEvents().

TH1D* cheat::RecoCheckAna::fEventPurityEfficiency

private

histogram of event efficiency times purity

Definition at line 131 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 136 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 137 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 138 of file RecoCheckAna_module.cc.

Referenced by CheckRecoTracks(), and FillResults().

double cheat::RecoCheckAna::fhiteff

private

hitfinder efficiency for this particle

Definition at line 146 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 101 of file RecoCheckAna_module.cc.

Referenced by analyze().

int cheat::RecoCheckAna::fnclu

private

number of clusters for this particle

Definition at line 147 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fnshw

private

number of showers for this particle

Definition at line 151 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fntrk

private

number of tracks for this particle

Definition at line 155 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::fpdg

private

particle pdg code

Definition at line 144 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

the back tracker service

Definition at line 99 of file RecoCheckAna_module.cc.

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

double cheat::RecoCheckAna::fpmom

private

particle momentum

Definition at line 145 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::frun

private

run number

Definition at line 141 of file RecoCheckAna_module.cc.

Referenced by analyze(), and beginRun().

TH1D* cheat::RecoCheckAna::fShowerEfficiency

private

histogram of shower efficiency

Definition at line 119 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 103 of file RecoCheckAna_module.cc.

Referenced by analyze(), and RecoCheckAna().

TH1D* cheat::RecoCheckAna::fShowerPurity

private

histogram of shower purity

Definition at line 118 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fShowerPurityEfficiency

private

histogram of shower efficiency times purity

Definition at line 120 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 121 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

shower efficiencies

Definition at line 152 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

shower IDs

Definition at line 154 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

shower purities

Definition at line 153 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fTrackEfficiency

private

histogram of track efficiency

Definition at line 123 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

int cheat::RecoCheckAna::ftrackid

private

geant track ID

Definition at line 143 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 104 of file RecoCheckAna_module.cc.

Referenced by analyze(), and RecoCheckAna().

TH1D* cheat::RecoCheckAna::fTrackPurity

private

histogram of track purity

Definition at line 122 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fTrackPurityEfficiency

private

histogram of track efficiency times purity

Definition at line 124 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 125 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TTree* cheat::RecoCheckAna::fTree

private

TTree to save efficiencies.

Definition at line 140 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

track efficiencies

Definition at line 156 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

track IDs

Definition at line 158 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

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

private

track purities

Definition at line 157 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and FillResults().

TH1D* cheat::RecoCheckAna::fVertexEfficiency

private

histogram of vertex efficiency

Definition at line 127 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 105 of file RecoCheckAna_module.cc.

Referenced by analyze(), and RecoCheckAna().

TH1D* cheat::RecoCheckAna::fVertexPurity

private

histogram of vertex purity

Definition at line 126 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoVertices().

TH1D* cheat::RecoCheckAna::fVertexPurityEfficiency

private

histogram of vertex efficiency times purity

Definition at line 128 of file RecoCheckAna_module.cc.

Referenced by beginRun(), and CheckRecoVertices().

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

larsim/v09_43_00/source/larsim/MCCheater/RecoCheckAna_module.cc

Public Types

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