Inheritance diagram for ShowerQuality:

Classes
struct	TreeParams_t
	For convenience: struct to define a set of parameters per shower to be stored in TTree. More...

Public Types
using	ModuleType = EDAnalyzer

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

	ShowerQuality (ShowerQuality const &)=delete

	ShowerQuality (ShowerQuality &&)=delete

ShowerQuality &	operator= (ShowerQuality const &)=delete

ShowerQuality &	operator= (ShowerQuality &&)=delete

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	beginJob () override

void	SetShowerProducer (const std::string name)

void	SetMCShowerProducer (const std::string name)

void	SetSimChannelProducer (const std::string name)

void	SetMaxEnergyCut (const double energy)
	Set maximum energy for MCShowers to be considered. More...

void	SetMinEnergyCut (const double energy)
	Set minimum energy for MCShowers to be considered. More...

void	InitializeAnaTree ()
	Function to prepare TTree. More...

Private Attributes
::btutil::MCMatchAlg	fBTAlg
	Shower back tracking algorithm. More...

double	_mc_energy_min
	Minimum MC shower energy cut. More...

double	_mc_energy_max
	Maximum MC shower energy cut. More...

std::string	fShowerProducer
	Shower Producer's Name. More...

std::string	fMCShowerProducer
	MCShower Producer's Name. More...

std::string	fSimChannelProducer
	SimChannel Producer's Name. More...

TH1D *	hMatchCorrectness
	Matching correctness. More...

TH1D *	hVtxDX
	X difference (reco-MC) in cm. More...

TH1D *	hVtxDY
	Y difference (reco-MC) in cm. More...

TH1D *	hVtxDZ
	Z difference (reco-MC) in cm. More...

TH1D *	hVtxDR
	3D vtx distance between reco to MC in cm More...

TH1D *	hDCosX
	Direction unit vector X component difference. More...

TH1D *	hDCosY
	Direction unit vector Y component difference. More...

TH1D *	hDCosZ
	Direction unit vector Z component difference. More...

TH1D *	h3DAngleDiff
	Opening angle between reco & MC 3D direction. More...

TH2D *	hEnergyCorr
	Energy correlation reco (x) vs. MC (y) More...

TH1D *	hEnergyAssym
	Energy assym. parameter: (reco E - MC E) / (reco E + MC E) * 2. More...

TH1D *	hEnergyDiff
	Energy difference: reco E - MC E. More...

TH1D *	hMatchedClusterEff
	Matched 3D shower's cluster efficiency (combined across planes) More...

TH1D *	hMatchedClusterPur
	Matched 3D shower's cluster purity (combined across planes) More...

std::map< int, TH1D * >	mDEDX
	dEdx per particle per PDG code More...

TH1D *	hBestPlane
	Best plane id. More...

struct ShowerQuality::TreeParams_t	fTreeParams

TTree *	fTree
	Analysis TTree. More...

Detailed Description

Definition at line 34 of file ShowerQuality_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

ShowerQuality::ShowerQuality ( fhicl::ParameterSet const & p )

explicit

Definition at line 146 of file ShowerQuality_module.cc.

References fTree, fhicl::ParameterSet::get(), h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, mDEDX, SetMaxEnergyCut(), SetMCShowerProducer(), SetMinEnergyCut(), SetShowerProducer(), and SetSimChannelProducer().

                                                        : EDAnalyzer(p) // ,
 // More initializers here.
 {
 
   //fShowerProducer   = "";
   //fMCShowerProducer = "";
   //fSimChannelProducer = "";
   SetShowerProducer(p.get<std::string>("ShowerProducer"));
   SetMCShowerProducer(p.get<std::string>("MCShowerProducer"));
   SetSimChannelProducer(p.get<std::string>("SimChannelProducer"));
   SetMinEnergyCut(p.get<double>("MCShowerEnergyMin"));
   SetMaxEnergyCut(p.get<double>("MCShowerEnergyMax"));
 
   hMatchCorrectness = nullptr;
 
   hVtxDX = nullptr;
   hVtxDY = nullptr;
   hVtxDZ = nullptr;
   hVtxDR = nullptr;
 
   hDCosX = nullptr;
   hDCosY = nullptr;
   hDCosZ = nullptr;
   h3DAngleDiff = nullptr;
 
   hEnergyCorr = nullptr;
   hEnergyAssym = nullptr;
   hEnergyDiff = nullptr;
 
   hMatchedClusterPur = nullptr;
   hMatchedClusterEff = nullptr;
 
   mDEDX.clear();
   hBestPlane = nullptr;
 
   fTree = nullptr;
 }

ShowerQuality::ShowerQuality ( ShowerQuality const & )

delete

ShowerQuality::ShowerQuality ( ShowerQuality && )

delete

Member Function Documentation

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

overrideprivatevirtual

Implements art::EDAnalyzer.

Definition at line 314 of file ShowerQuality_module.cc.

 {
   // Retrieve mcshower data product
   auto mcsHandle = e.getValidHandle<std::vector<sim::MCShower>>(fMCShowerProducer);
   auto resHandle = e.getValidHandle<std::vector<recob::Shower>>(fShowerProducer);
   auto schHandle = e.getValidHandle<std::vector<sim::SimChannel>>(fSimChannelProducer);
   const std::vector<sim::MCShower>& ev_mcs(*mcsHandle);
   const std::vector<recob::Shower>& ev_shower(*resHandle);
   const std::vector<sim::SimChannel>& ev_simch(*schHandle);
 
   if (ev_shower.empty()) return;
 
   // Get the whole clusters + associated clusters
   art::Handle<std::vector<recob::Cluster>> clsHandle;
   art::FindManyP<recob::Cluster> cluster_m(resHandle, e, fShowerProducer);
   const std::vector<recob::Cluster>& ev_cluster = cluster_m.at(0).front().parentAs<std::vector>();
 
   // Make clusters in terms of hit vector to feed into BT algorithm
   art::FindManyP<recob::Hit> hit_m(clsHandle, e, clsHandle.provenance()->moduleLabel());
   std::vector<std::vector<art::Ptr<recob::Hit>>> ev_cluster_hit;
   ev_cluster_hit.reserve(clsHandle->size());
   std::map<art::Ptr<recob::Cluster>, size_t> cluster_ptr_map;
   for (size_t i = 0; i < ev_cluster.size(); ++i) {
     const art::Ptr<recob::Cluster> cluster_ptr(clsHandle, i);
     cluster_ptr_map[cluster_ptr] = ev_cluster_hit.size();
     ev_cluster_hit.push_back(hit_m.at(i));
   }
 
   // Create ass_cluster_v index vector
   std::vector<std::vector<unsigned int>> ass_cluster_v;
   ass_cluster_v.reserve(ev_shower.size());
   for (size_t shower_index = 0; shower_index < ev_shower.size(); ++shower_index) {
     ass_cluster_v.push_back(std::vector<unsigned int>());
     for (auto const& p : cluster_m.at(shower_index))
       ass_cluster_v.back().push_back(cluster_ptr_map[p]);
   }
 
   // Create G4 track ID vector for which we are interested in
   std::vector<std::vector<unsigned int>> g4_trackid_v;
   std::vector<unsigned int> mc_index_v;
   g4_trackid_v.reserve(ev_mcs.size());
   for (size_t mc_index = 0; mc_index < ev_mcs.size(); ++mc_index) {
     auto const& mcs = ev_mcs[mc_index];
     double energy = mcs.DetProfile().E();
     std::vector<unsigned int> id_v;
     id_v.reserve(mcs.DaughterTrackID().size());
     if (_mc_energy_min < energy && energy < _mc_energy_max) {
       for (auto const& id : mcs.DaughterTrackID()) {
         if (id == mcs.TrackID()) continue;
         id_v.push_back(id);
       }
       id_v.push_back(mcs.TrackID());
       g4_trackid_v.push_back(id_v);
       mc_index_v.push_back(mc_index);
     }
   }
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
   if (!fBTAlg.BuildMap(clockData, g4_trackid_v, ev_simch, ev_cluster_hit)) {
     std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> Failed to "
                  "build back-tracking map for MC..."
               << std::endl;
     return;
   }
 
   // Find the best-representative reco-ed Shower given an MCShower
   std::vector<std::vector<double>> shower_mcq_vv(ev_shower.size(),
                                                  std::vector<double>(mc_index_v.size(), 0));
 
   for (size_t shower_index = 0; shower_index < ass_cluster_v.size(); ++shower_index) {
 
     auto const& ass_cluster = ass_cluster_v[shower_index];
 
     std::vector<::btutil::WireRange_t> w_v;
 
     for (auto const& cluster_index : ass_cluster) {
 
       auto const& ass_hit = ev_cluster_hit[cluster_index];
 
       w_v.reserve(ass_hit.size() + w_v.size());
 
       for (auto const& hit_ptr : ass_hit) {
         w_v.emplace_back(hit_ptr->Channel(), hit_ptr->StartTick(), hit_ptr->EndTick());
       }
     }
 
     auto mcq_v = fBTAlg.BTAlg().MCQ(clockData, w_v);
 
     auto& shower_mcq_v = shower_mcq_vv[shower_index];
 
     for (size_t mcs_index = 0; mcs_index < (mcq_v.size() - 1); ++mcs_index) {
 
       shower_mcq_v[mcs_index] = mcq_v[mcs_index];
     }
   }
 
   // Loop over MCShower and inspect corresponding shower quality
   for (size_t mcs_index = 0; mcs_index < mc_index_v.size(); ++mcs_index) {
 
     auto const& mc_shower = ev_mcs[mc_index_v[mcs_index]];
 
     // Search for the best representative shower
     size_t best_shower_index = shower_mcq_vv.size();
     double max_mcq = 0;
     for (size_t shower_index = 0; shower_index < shower_mcq_vv.size(); ++shower_index) {
 
       if (shower_mcq_vv[shower_index][mcs_index] > max_mcq) best_shower_index = shower_index;
     }
 
     if (best_shower_index == shower_mcq_vv.size()) {
       std::string msg;
       std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> "
                 << "Failed to find a corresponding shower for MCShower " << mc_index_v[mcs_index]
                 << std::endl;
       continue;
     }
 
     auto const& reco_shower = ev_shower[best_shower_index];
 
     auto res = fBTAlg.ShowerCorrectness(ass_cluster_v[best_shower_index]);
 
     fTreeParams.match_correctness = res.second;
 
     if (fTreeParams.match_correctness < 0) {
       std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> "
                 << "Failed to find a corresponding MCShower for shower " << best_shower_index
                 << std::endl;
       continue;
     }
 
     // MC Info
     fTreeParams.mc_x = mc_shower.DetProfile().X();
     fTreeParams.mc_y = mc_shower.DetProfile().Y();
     fTreeParams.mc_z = mc_shower.DetProfile().Z();
 
     fTreeParams.mc_energy = mc_shower.DetProfile().E();
     fTreeParams.mc_pdgid = mc_shower.PdgCode();
     fTreeParams.mc_containment = mc_shower.DetProfile().E() / mc_shower.Start().E();
 
     //fTreeParams.mc_dcosx = mc_shower.DetProfile().Px() / fTreeParams.mc_energy;
     //fTreeParams.mc_dcosy = mc_shower.DetProfile().Py() / fTreeParams.mc_energy;
     //fTreeParams.mc_dcosz = mc_shower.DetProfile().Pz() / fTreeParams.mc_energy;
     fTreeParams.mc_dcosx = mc_shower.Start().Px() / mc_shower.Start().E();
     fTreeParams.mc_dcosy = mc_shower.Start().Py() / mc_shower.Start().E();
     fTreeParams.mc_dcosz = mc_shower.Start().Pz() / mc_shower.Start().E();
 
     // Reco vtx
     fTreeParams.reco_x = reco_shower.ShowerStart()[0];
     fTreeParams.reco_y = reco_shower.ShowerStart()[1];
     fTreeParams.reco_z = reco_shower.ShowerStart()[2];
 
     // Reco angle
     fTreeParams.reco_dcosx = reco_shower.Direction()[0];
     fTreeParams.reco_dcosy = reco_shower.Direction()[1];
     fTreeParams.reco_dcosz = reco_shower.Direction()[2];
 
     // Reco - MC angle diff
     fTreeParams.mc_reco_anglediff = acos(fTreeParams.reco_dcosx * fTreeParams.mc_dcosx +
                                          fTreeParams.reco_dcosy * fTreeParams.mc_dcosy +
                                          fTreeParams.reco_dcosz * fTreeParams.mc_dcosz) /
                                     3.14159265359 * 180.;
     // Reco - MC vtx distance
     fTreeParams.mc_reco_dist = sqrt(pow(fTreeParams.reco_x - fTreeParams.mc_x, 2) +
                                     pow(fTreeParams.reco_y - fTreeParams.mc_y, 2) +
                                     pow(fTreeParams.reco_z - fTreeParams.mc_z, 2));
 
     // Reco cluster efficiency & purity
     fTreeParams.cluster_eff = 1.;
     fTreeParams.cluster_pur = 1.;
     for (auto const& cluster_index : ass_cluster_v[best_shower_index]) {
       auto ep = fBTAlg.ClusterEP(cluster_index, mcs_index);
       if (ep.first == 0 && ep.second == 0) continue;
       fTreeParams.cluster_eff *= ep.first;
       fTreeParams.cluster_pur *= ep.second;
     }
 
     // Reco energy & dedx info
     fTreeParams.best_plane_id = reco_shower.best_plane();
 
     /*
       int best_plane_index = -1;
 
       for(size_t i=0; i < ass_cluster_v[best_shower_index].size(); ++i) {
 
         size_t cluster_index = ass_cluster_v[best_shower_index][i];
         //std::cout<<best_plane_index<<" : "<<ev_cluster->at(cluster_index).View()<<std::endl;
         if( ev_cluster->at(cluster_index).View() == reco_shower.best_plane() ) {
           best_plane_index = i;
           break;
         }
       }
 
       if(best_plane_index < 0) {
       throw ::showerreco::ShowerRecoException(Form("Failed to identify the best plane for shower %zu",
       best_shower_index)
       );
       }
     */
 
     fTreeParams.reco_energy = reco_shower.Energy().at(reco_shower.best_plane());
     //fTreeParams.reco_dedx_U   = reco_shower.dEdx().at(0);
     //fTreeParams.reco_dedx_V   = reco_shower.dEdx().at(1);
     //fTreeParams.reco_dedx_Y   = reco_shower.dEdx().at(2);
     fTreeParams.reco_dedx = reco_shower.dEdx().at(reco_shower.best_plane());
 
     //
     // Fill histograms
     //
     hMatchCorrectness->Fill(fTreeParams.match_correctness);
 
     hVtxDX->Fill(fTreeParams.reco_x - fTreeParams.mc_x);
     hVtxDY->Fill(fTreeParams.reco_y - fTreeParams.mc_y);
     hVtxDZ->Fill(fTreeParams.reco_z - fTreeParams.mc_z);
     hVtxDR->Fill(fTreeParams.mc_reco_dist);
 
     // Angular info
     hDCosX->Fill(fTreeParams.reco_dcosx - fTreeParams.mc_dcosx);
     hDCosY->Fill(fTreeParams.reco_dcosy - fTreeParams.mc_dcosy);
     hDCosZ->Fill(fTreeParams.reco_dcosz - fTreeParams.mc_dcosz);
     h3DAngleDiff->Fill(fTreeParams.mc_reco_anglediff);
 
     hEnergyCorr->Fill(fTreeParams.reco_energy, fTreeParams.mc_energy);
 
     hEnergyAssym->Fill((fTreeParams.reco_energy - fTreeParams.mc_energy) /
                        (fTreeParams.reco_energy + fTreeParams.mc_energy) * 2.);
 
     hEnergyDiff->Fill(fTreeParams.mc_energy - fTreeParams.reco_energy);
 
     if (mDEDX.find(fTreeParams.mc_pdgid) == mDEDX.end())
 
       mDEDX.insert(std::make_pair(
         fTreeParams.mc_pdgid,
         new TH1D(Form("hdEdx_PDG_%d", fTreeParams.mc_pdgid),
                  Form("Reco dE/dx for PDG = %d; dE/dx [MeV/cm]; Showers", fTreeParams.mc_pdgid),
                  100,
                  0,
                  50)));
 
     hMatchedClusterEff->Fill(fTreeParams.cluster_eff);
     hMatchedClusterPur->Fill(fTreeParams.cluster_pur);
 
     mDEDX[fTreeParams.mc_pdgid]->Fill(fTreeParams.reco_dedx);
 
     hBestPlane->Fill(fTreeParams.best_plane_id);
 
     // Fill Tree
     fTree->Fill();
   }
 }

void ShowerQuality::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 184 of file ShowerQuality_module.cc.

References fMCShowerProducer, fShowerProducer, fSimChannelProducer, fTree, h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, and InitializeAnaTree().

 {
 
   if (fShowerProducer.empty() || fMCShowerProducer.empty() || fSimChannelProducer.empty()) {
     std::string msg;
     msg += "\033[93m[ERROR]\033[00m <<";
     msg += __FUNCTION__;
     msg += ">> Producer's name not set!";
     std::cout << msg.c_str() << std::endl;
     throw ::showerreco::ShowerRecoException(msg.c_str());
   }
 
   art::ServiceHandle<geo::Geometry const> geo;
   art::ServiceHandle<art::TFileService const> tfs;
 
   if (fTree) delete fTree;
   fTree = tfs->make<TTree>("fShowerQualityTree", "");
 
   //
   // Matching correctness histogram initialization
   //
   if (hMatchCorrectness) delete hMatchCorrectness;
   hMatchCorrectness =
     tfs->make<TH1D>("hMatchCorrectness",
                     "Shower 2D Cluster Matching Correctness; Correctness; Showers",
                     101,
                     -0.005,
                     1.005);
 
   //
   // 3D Vtx (start point) MC/Reco comparison histogram initialization
   //
   if (hVtxDX) delete hVtxDX;
   if (hVtxDY) delete hVtxDY;
   if (hVtxDZ) delete hVtxDZ;
   if (hVtxDR) delete hVtxDR;
 
   hVtxDX = tfs->make<TH1D>(
     "hVtxDX", "Reco - MC Start X [cm] Displacement; #DeltaX [cm]; Showers", 200, -100, 100);
 
   hVtxDY = tfs->make<TH1D>(
     "hVtxDY", "Reco - MC Start Y [cm] Displacement; #DeltaY [cm]; Showers", 200, -100, 100);
 
   hVtxDZ = tfs->make<TH1D>(
     "hVtxDZ", "Reco - MC Start Z [cm] Displacement; #DeltaZ [cm]; Showers", 200, -100, 100);
 
   hVtxDR = tfs->make<TH1D>(
     "hVtxDR", "Reco - MC Start 3D Vtx Displacement; #DeltaR [cm]; Showers", 200, -100, 100);
 
   //
   // 3D Angular MC/Reco comparison histogram initialization
   //
   if (hDCosX) delete hDCosX;
   if (hDCosY) delete hDCosY;
   if (hDCosZ) delete hDCosZ;
   if (h3DAngleDiff) delete h3DAngleDiff;
 
   hDCosX = tfs->make<TH1D>(
     "hDCosX", "Direction Unit Vector Reco - MC #DeltaX; #DeltaCosX; Showers", 100, -2, 2);
 
   hDCosY = tfs->make<TH1D>(
     "hDCosY", "Direction Unit Vector Reco - MC #DeltaY; #DeltaCosY; Showers", 100, -2, 2);
 
   hDCosZ = tfs->make<TH1D>(
     "hDCosZ", "Direction Unit Vector Reco - MC #DeltaZ; #DeltaCosZ; Showers", 100, -2, 2);
 
   h3DAngleDiff =
     tfs->make<TH1D>("h3DAngleDiff",
                     "3D Opening Angle Between Reco & MC; Opening Angle [degrees]; Showers",
                     181,
                     -0.5,
                     180.5);
 
   //
   // Energy MC/Reco comparison histogram initialization
   //
   if (hEnergyCorr) delete hEnergyCorr;
   if (hEnergyAssym) delete hEnergyAssym;
   if (hEnergyDiff) delete hEnergyDiff;
 
   hEnergyCorr =
     tfs->make<TH2D>("hEnergyCorr",
                     "Reco (x) vs. MC (y) Energy Comparison; Reco Energy [MeV]; MC Energy [MeV]",
                     200,
                     0,
                     1000,
                     200,
                     0,
                     1000);
 
   hEnergyAssym = tfs->make<TH1D>("hEnergyAssym",
                                  "MC - Reco Energy Fractional Difference; Assymetry; Showers",
                                  201,
                                  -1.005,
                                  1.005);
 
   hEnergyDiff = tfs->make<TH1D>(
     "hEnergyDiff", "MC - Reco Energy Difference; Energy Difference [MeV]; Showers", 200, 0, 1000);
 
   //
   // Shower cluster purity & efficiency histograms initialization
   //
   if (hMatchedClusterEff) delete hMatchedClusterEff;
   if (hMatchedClusterPur) delete hMatchedClusterPur;
 
   hMatchedClusterEff =
     tfs->make<TH1D>("hMatchedClusterEff_PlaneCombo",
                     "Matched Shower Cluster's Charge Efficiency; Efficiency; Clusters",
                     101,
                     -0.005,
                     1.005);
 
   hMatchedClusterPur = tfs->make<TH1D>("hMatchedClusterPur_PlaneCombo",
                                        "Matched Shower Cluster's Charge Purity; Purity; Clusters",
                                        101,
                                        -0.005,
                                        1.005);
 
   //
   // Best plane ID histogram initialization
   //
   hBestPlane = tfs->make<TH1D>("hBestPlane",
                                "Best Plane (for energy & dE/dx estimate); Plane ID; Showers",
                                geo->Nplanes(),
                                -0.5,
                                geo->Nplanes() - 0.5);
 
   InitializeAnaTree();
 }

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

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

void ShowerQuality::InitializeAnaTree ( )

private

Function to prepare TTree.

Definition at line 564 of file ShowerQuality_module.cc.

Referenced by beginJob().

 {
 
   fTree->Branch("reco_x", &fTreeParams.reco_x, "reco_x/D");
   fTree->Branch("reco_y", &fTreeParams.reco_y, "reco_y/D");
   fTree->Branch("reco_z", &fTreeParams.reco_z, "reco_z/D");
   fTree->Branch("reco_dcosx", &fTreeParams.reco_dcosx, "reco_dcosx/D");
   fTree->Branch("reco_dcosy", &fTreeParams.reco_dcosy, "reco_dcosy/D");
   fTree->Branch("reco_dcosz", &fTreeParams.reco_dcosz, "reco_dcosz/D");
   fTree->Branch("reco_energy", &fTreeParams.reco_energy, "reco_energy/D");
 
   fTree->Branch("best_plane_id", &fTreeParams.best_plane_id, "best_plane_id/i");
 
   fTree->Branch("mc_x", &fTreeParams.mc_x, "mc_x/D");
   fTree->Branch("mc_y", &fTreeParams.mc_y, "mc_y/D");
   fTree->Branch("mc_z", &fTreeParams.mc_z, "mc_z/D");
   fTree->Branch("mc_dcosx", &fTreeParams.mc_dcosx, "mc_dcosx/D");
   fTree->Branch("mc_dcosy", &fTreeParams.mc_dcosy, "mc_dcosy/D");
   fTree->Branch("mc_dcosz", &fTreeParams.mc_dcosz, "mc_dcosz/D");
   fTree->Branch("mc_energy", &fTreeParams.mc_energy, "mc_energy/D");
 
   fTree->Branch("reco_dedx", &fTreeParams.reco_dedx, "reco_dedx_/D");
   fTree->Branch("reco_dedx_U", &fTreeParams.reco_dedx_U, "reco_dedx_U/D");
   fTree->Branch("reco_dedx_V", &fTreeParams.reco_dedx_V, "reco_dedx_V/D");
   fTree->Branch("reco_dedx_Y", &fTreeParams.reco_dedx_Y, "reco_dedx_Y/D");
   fTree->Branch("mc_pdgid", &fTreeParams.mc_pdgid, "mc_pdgid/i");
 
   fTree->Branch("mc_reco_anglediff", &fTreeParams.mc_reco_anglediff, "mc_reco_anglediff/D");
   fTree->Branch("mc_reco_dist", &fTreeParams.mc_reco_dist, "mc_reco_dist/D");
 
   fTree->Branch("mc_containment", &fTreeParams.mc_containment, "mc_containment/D");
 
   fTree->Branch("match_correctness", &fTreeParams.match_correctness, "match_correctness/D");
   fTree->Branch("cluster_eff", &fTreeParams.cluster_eff, "cluster_eff/D");
   fTree->Branch("cluster_pur", &fTreeParams.cluster_pur, "cluster_pur/D");
 }

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

ShowerQuality& ShowerQuality::operator= ( ShowerQuality const & )

delete

ShowerQuality& ShowerQuality::operator= ( ShowerQuality && )

delete

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 ShowerQuality::SetMaxEnergyCut ( const double energy )

inlineprivate

Set maximum energy for MCShowers to be considered.

Definition at line 60 of file ShowerQuality_module.cc.

References _mc_energy_max, and energy.

Referenced by ShowerQuality().

60 { _mc_energy_max = energy; }

energy

double energy

Definition: plottest35.C:25

ShowerQuality::_mc_energy_max

double _mc_energy_max

Maximum MC shower energy cut.

Definition: ShowerQuality_module.cc:72

void ShowerQuality::SetMCShowerProducer ( const std::string name )

inlineprivate

Definition at line 56 of file ShowerQuality_module.cc.

References fMCShowerProducer.

Referenced by ShowerQuality().

56 { fMCShowerProducer = name; }

ShowerQuality::fMCShowerProducer

std::string fMCShowerProducer

MCShower Producer's Name.

Definition: ShowerQuality_module.cc:78

void ShowerQuality::SetMinEnergyCut ( const double energy )

inlineprivate

Set minimum energy for MCShowers to be considered.

Definition at line 63 of file ShowerQuality_module.cc.

References _mc_energy_min, and energy.

Referenced by ShowerQuality().

63 { _mc_energy_min = energy; }

energy

double energy

Definition: plottest35.C:25

ShowerQuality::_mc_energy_min

double _mc_energy_min

Minimum MC shower energy cut.

Definition: ShowerQuality_module.cc:69

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

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void ShowerQuality::SetShowerProducer ( const std::string name )

inlineprivate

Setter function for a shower producer name.

Definition at line 55 of file ShowerQuality_module.cc.

References fShowerProducer.

Referenced by ShowerQuality().

55 { fShowerProducer = name; }

ShowerQuality::fShowerProducer

std::string fShowerProducer

Shower Producer's Name.

Definition: ShowerQuality_module.cc:75

void ShowerQuality::SetSimChannelProducer ( const std::string name )

inlineprivate

Definition at line 57 of file ShowerQuality_module.cc.

References fSimChannelProducer.

Referenced by ShowerQuality().

57 { fSimChannelProducer = name; }

ShowerQuality::fSimChannelProducer

std::string fSimChannelProducer

SimChannel Producer's Name.

Definition: ShowerQuality_module.cc:81

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