ShowerQuality Class Reference

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	WorkerType = WorkerT< EDAnalyzer >
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	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...
template<class T >
art::Handle< T >	GetDataOrDie (art::Event const &e, std::string producer)
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	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 38 of file ShowerQuality_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

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

explicit

Definition at line 170 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 )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 332 of file ShowerQuality_module.cc.

References _mc_energy_max, _mc_energy_min, ShowerQuality::TreeParams_t::best_plane_id, btutil::MCMatchAlg::BTAlg(), btutil::MCMatchAlg::BuildMap(), ShowerQuality::TreeParams_t::cluster_eff, ShowerQuality::TreeParams_t::cluster_pur, btutil::MCMatchAlg::ClusterEP(), e, energy, fBTAlg, fMCShowerProducer, fShowerProducer, fSimChannelProducer, fTree, fTreeParams, h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, art::Handle< T >::isValid(), ShowerQuality::TreeParams_t::match_correctness, ShowerQuality::TreeParams_t::mc_containment, ShowerQuality::TreeParams_t::mc_dcosx, ShowerQuality::TreeParams_t::mc_dcosy, ShowerQuality::TreeParams_t::mc_dcosz, ShowerQuality::TreeParams_t::mc_energy, ShowerQuality::TreeParams_t::mc_pdgid, ShowerQuality::TreeParams_t::mc_reco_anglediff, ShowerQuality::TreeParams_t::mc_reco_dist, ShowerQuality::TreeParams_t::mc_x, ShowerQuality::TreeParams_t::mc_y, ShowerQuality::TreeParams_t::mc_z, btutil::MCBTAlg::MCQ(), mDEDX, art::Handle< T >::provenance(), ShowerQuality::TreeParams_t::reco_dcosx, ShowerQuality::TreeParams_t::reco_dcosy, ShowerQuality::TreeParams_t::reco_dcosz, ShowerQuality::TreeParams_t::reco_dedx, ShowerQuality::TreeParams_t::reco_energy, ShowerQuality::TreeParams_t::reco_x, ShowerQuality::TreeParams_t::reco_y, ShowerQuality::TreeParams_t::reco_z, and btutil::MCMatchAlg::ShowerCorrectness().

{
  //auto geo = larutil::Geometry::GetME();
  
  // Retrieve mcshower data product
  auto mcsHandle = GetDataOrDie< std::vector<sim::MCShower  > > (e,fMCShowerProducer);
  auto resHandle = GetDataOrDie< std::vector<recob::Shower  > > (e,fShowerProducer);
  auto schHandle = GetDataOrDie< std::vector<sim::SimChannel> > (e,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.size())) return;

  //auto clsHandle = GetDataOrDie<recob::Cluster>  (e,fClusterProducer);

  //art::FindManyP<recob::Hit>     hit_m     (resHandle, e, fShowerProducer);

  // Get the whole clusters + associated clusters
  art::Handle<std::vector<recob::Cluster> > clsHandle;
  art::FindManyP<recob::Cluster> cluster_m (resHandle, e, fShowerProducer);
  e.get(cluster_m.at(0).front().id(),clsHandle);
  if(!clsHandle.isValid()) throw ::showerreco::ShowerRecoException("Failed to retrieve cluster handle!");
  const std::vector<recob::Cluster>& ev_cluster (*clsHandle);

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

  if(!fBTAlg.BuildMap(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.push_back( ::btutil::WireRange_t( hit_ptr->Channel(),
                                              hit_ptr->StartTick(),
                                              hit_ptr->EndTick() ) 
                       );
      }
    }
    
    auto mcq_v = fBTAlg.BTAlg().MCQ(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 ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 210 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, InitializeAnaTree(), and art::TFileDirectory::make().

{

  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> geo;
  art::ServiceHandle<art::TFileService> 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();
  
}

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

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

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

  {
    return selectors_;
  }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

{
  return rng()->createEngine(placeholder_schedule_id(), seed);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make);
}

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

{
  return rng()->createEngine(
    placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
}

CurrentProcessingContext const * art::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

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

  {
    return current_context_.get();
  }

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

template<class T >

art::Handle<T> ShowerQuality::GetDataOrDie ( art::Event const &  e, std::string  producer  )

inline

Definition at line 72 of file ShowerQuality_module.cc.

References art::DataViewImpl::getByLabel(), and art::Handle< T >::isValid().

  {
    art::Handle<T> h;
    e.getByLabel(producer,h);
    if(!h.isValid()) {
      std::string msg;
      msg += "Could not find a data product by: " + producer;
      std::cout<<msg.c_str()<<std::endl;
      throw showerreco::ShowerRecoException(msg);
    }
    return h;
  }

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

void ShowerQuality::InitializeAnaTree ( )

private

Function to prepare TTree.

Definition at line 592 of file ShowerQuality_module.cc.

References ShowerQuality::TreeParams_t::best_plane_id, ShowerQuality::TreeParams_t::cluster_eff, ShowerQuality::TreeParams_t::cluster_pur, DEFINE_ART_MODULE, fTree, fTreeParams, ShowerQuality::TreeParams_t::match_correctness, ShowerQuality::TreeParams_t::mc_containment, ShowerQuality::TreeParams_t::mc_dcosx, ShowerQuality::TreeParams_t::mc_dcosy, ShowerQuality::TreeParams_t::mc_dcosz, ShowerQuality::TreeParams_t::mc_energy, ShowerQuality::TreeParams_t::mc_pdgid, ShowerQuality::TreeParams_t::mc_reco_anglediff, ShowerQuality::TreeParams_t::mc_reco_dist, ShowerQuality::TreeParams_t::mc_x, ShowerQuality::TreeParams_t::mc_y, ShowerQuality::TreeParams_t::mc_z, ShowerQuality::TreeParams_t::reco_dcosx, ShowerQuality::TreeParams_t::reco_dcosy, ShowerQuality::TreeParams_t::reco_dcosz, ShowerQuality::TreeParams_t::reco_dedx, ShowerQuality::TreeParams_t::reco_dedx_U, ShowerQuality::TreeParams_t::reco_dedx_V, ShowerQuality::TreeParams_t::reco_dedx_Y, ShowerQuality::TreeParams_t::reco_energy, ShowerQuality::TreeParams_t::reco_x, ShowerQuality::TreeParams_t::reco_y, and ShowerQuality::TreeParams_t::reco_z.

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

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.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

{
  static Consumer invalid{InvalidTag{}};
  return &invalid;
}

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

delete

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

delete

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

{
  if (!moduleContext_)
    return;

  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
  for (auto& consumablesPerBranch : consumables_) {
    cet::sort_all(consumablesPerBranch);
  }
}

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

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

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

  {
    return process_name_;
  }

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

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

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

inline

Set maximum energy for MCShowers to be considered.

Definition at line 66 of file ShowerQuality_module.cc.

References _mc_energy_max, and energy.

Referenced by ShowerQuality().

{ _mc_energy_max = energy; }

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

inline

Definition at line 60 of file ShowerQuality_module.cc.

References fMCShowerProducer.

Referenced by ShowerQuality().

  { fMCShowerProducer = name; }

void ShowerQuality::SetMinEnergyCut ( const double  energy )

inline

Set minimum energy for MCShowers to be considered.

Definition at line 69 of file ShowerQuality_module.cc.

References _mc_energy_min, and energy.

Referenced by ShowerQuality().

{ _mc_energy_min = energy; }

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

inline

Setter function for a shower producer name.

Definition at line 58 of file ShowerQuality_module.cc.

References fShowerProducer.

Referenced by ShowerQuality().

  { fShowerProducer = name; }

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

inline

Definition at line 62 of file ShowerQuality_module.cc.

References fSimChannelProducer.

Referenced by ShowerQuality().

  { fSimChannelProducer = name; }

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

double ShowerQuality::_mc_energy_max

private

Maximum MC shower energy cut.

Definition at line 95 of file ShowerQuality_module.cc.

Referenced by analyze(), and SetMaxEnergyCut().

double ShowerQuality::_mc_energy_min

private

Minimum MC shower energy cut.

Definition at line 92 of file ShowerQuality_module.cc.

Referenced by analyze(), and SetMinEnergyCut().

::btutil::MCMatchAlg ShowerQuality::fBTAlg

private

Shower back tracking algorithm.

Definition at line 89 of file ShowerQuality_module.cc.

Referenced by analyze().

std::string ShowerQuality::fMCShowerProducer

private

MCShower Producer's Name.

Definition at line 101 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetMCShowerProducer().

std::string ShowerQuality::fShowerProducer

private

Shower Producer's Name.

Definition at line 98 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetShowerProducer().

std::string ShowerQuality::fSimChannelProducer

private

SimChannel Producer's Name.

Definition at line 104 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetSimChannelProducer().

TTree* ShowerQuality::fTree

private

Analysis TTree.

Definition at line 161 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), InitializeAnaTree(), and ShowerQuality().

struct ShowerQuality::TreeParams_t ShowerQuality::fTreeParams

private

Referenced by analyze(), and InitializeAnaTree().

TH1D* ShowerQuality::h3DAngleDiff

private

Opening angle between reco & MC 3D direction.

Definition at line 117 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hBestPlane

private

Best plane id.

Definition at line 131 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosX

private

Direction unit vector X component difference.

Definition at line 114 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosY

private

Direction unit vector Y component difference.

Definition at line 115 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosZ

private

Direction unit vector Z component difference.

Definition at line 116 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hEnergyAssym

private

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

Definition at line 121 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH2D* ShowerQuality::hEnergyCorr

private

Energy correlation reco (x) vs. MC (y)

Definition at line 119 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hEnergyDiff

private

Energy difference: reco E - MC E.

Definition at line 122 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchCorrectness

private

Matching correctness.

Definition at line 107 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchedClusterEff

private

Matched 3D shower's cluster efficiency (combined across planes)

Definition at line 124 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchedClusterPur

private

Matched 3D shower's cluster purity (combined across planes)

Definition at line 125 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDR

private

3D vtx distance between reco to MC in cm

Definition at line 112 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDX

private

X difference (reco-MC) in cm.

Definition at line 109 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDY

private

Y difference (reco-MC) in cm.

Definition at line 110 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDZ

private

Z difference (reco-MC) in cm.

Definition at line 111 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

std::map<int,TH1D*> ShowerQuality::mDEDX

private

dEdx per particle per PDG code

Definition at line 128 of file ShowerQuality_module.cc.

Referenced by analyze(), and ShowerQuality().

---

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

• larreco/v06_64_02/source/larreco/ShowerFinder/ShowerReco3D/ShowerQuality_module.cc