Inheritance diagram for ShowerRecoTools::ShowerNumElectronsEnergy:

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

int	CalculateElement (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerElementHolder) override

Private Member Functions
double	CalculateEnergy (const detinfo::DetectorClocksData &clockData, const detinfo::DetectorPropertiesData &detProp, const std::vector< art::Ptr< recob::Hit >> &hits, const geo::PlaneID::PlaneID_t plane) const

Private Member Functions inherited from ShowerRecoTools::IShowerTool
	IShowerTool (const fhicl::ParameterSet &pset)

virtual	~IShowerTool () noexcept=default

int	RunShowerTool (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder, std::string evd_display_name_append="")

virtual void	InitialiseProducers ()

void	SetPtr (art::ProducesCollector *collector)

void	InitaliseProducerPtr (reco::shower::ShowerProducedPtrsHolder &uniqueproducerPtrs)

virtual int	AddAssociations (const art::Ptr< recob::PFParticle > &pfpPtr, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)

const shower::LArPandoraShowerAlg &	GetLArPandoraShowerAlg () const

template<class T >
art::Ptr< T >	GetProducedElementPtr (std::string Name, reco::shower::ShowerElementHolder &ShowerEleHolder, int iter=-1)

template<class T >
void	InitialiseProduct (std::string Name, std::string InstanceName="")

template<class T , class A , class B >
void	AddSingle (A &a, B &b, std::string Name)

int	GetVectorPtrSize (std::string Name)

void	PrintPtrs ()

void	PrintPtr (std::string Name)

Private Attributes
art::InputTag	fPFParticleLabel

int	fVerbose

std::string	fShowerEnergyOutputLabel

std::string	fShowerBestPlaneOutputLabel

art::ServiceHandle< geo::Geometry >	fGeom

calo::CalorimetryAlg	fCalorimetryAlg

double	fRecombinationFactor

Detailed Description

Definition at line 25 of file ShowerNumElectronsEnergy_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerNumElectronsEnergy::ShowerNumElectronsEnergy ( const fhicl::ParameterSet & pset )

Definition at line 55 of file ShowerNumElectronsEnergy_tool.cc.

     : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
     , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
     , fVerbose(pset.get<int>("Verbose"))
     , fShowerEnergyOutputLabel(pset.get<std::string>("ShowerEnergyOutputLabel"))
     , fShowerBestPlaneOutputLabel(pset.get<std::string>("ShowerBestPlaneOutputLabel"))
     , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
     , fRecombinationFactor(pset.get<double>("RecombinationFactor"))
   {}

Member Function Documentation

int ShowerRecoTools::ShowerNumElectronsEnergy::CalculateElement	(	const art::Ptr< recob::PFParticle > &	pfparticle,
		art::Event &	Event,
		reco::shower::ShowerElementHolder &	ShowerElementHolder
	)

overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 66 of file ShowerNumElectronsEnergy_tool.cc.

   {
 
     if (fVerbose)
       std::cout
         << "~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Shower Reco Energy Tool ~~~~~~~~~~~~~~~~~~~~~~~~~~~~"
         << std::endl;
 
     // Get the assocated pfParicle vertex PFParticles
     auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);
 
     //Get the clusters
     auto const clusHandle = Event.getValidHandle<std::vector<recob::Cluster>>(fPFParticleLabel);
 
     const art::FindManyP<recob::Cluster>& fmc =
       ShowerEleHolder.GetFindManyP<recob::Cluster>(pfpHandle, Event, fPFParticleLabel);
     // art::FindManyP<recob::Cluster> fmc(pfpHandle, Event, fPFParticleLabel);
     std::vector<art::Ptr<recob::Cluster>> clusters = fmc.at(pfparticle.key());
 
     //Get the hit association
     const art::FindManyP<recob::Hit>& fmhc =
       ShowerEleHolder.GetFindManyP<recob::Hit>(clusHandle, Event, fPFParticleLabel);
     // art::FindManyP<recob::Hit> fmhc(clusHandle, Event, fPFParticleLabel);
 
     std::map<geo::PlaneID::PlaneID_t, std::vector<art::Ptr<recob::Hit>>> planeHits;
 
     //Loop over the clusters in the plane and get the hits
     for (auto const& cluster : clusters) {
 
       //Get the hits
       std::vector<art::Ptr<recob::Hit>> hits = fmhc.at(cluster.key());
 
       //Get the plane.
       const geo::PlaneID::PlaneID_t plane(cluster->Plane().Plane);
 
       planeHits[plane].insert(planeHits[plane].end(), hits.begin(), hits.end());
     }
 
     // Calculate the energy for each plane && best plane
     geo::PlaneID::PlaneID_t bestPlane = std::numeric_limits<geo::PlaneID::PlaneID_t>::max();
     unsigned int bestPlaneNumHits = 0;
 
     //Holder for the final product
     std::vector<double> energyVec(fGeom->Nplanes(), -999.);
     std::vector<double> energyError(fGeom->Nplanes(), -999.);
 
     auto const clockData =
       art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(Event);
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(Event, clockData);
 
     for (auto const& [plane, hits] : planeHits) {
 
       unsigned int planeNumHits = hits.size();
 
       //Calculate the Energy for
       double Energy = CalculateEnergy(clockData, detProp, hits, plane);
       // If the energy is negative, leave it at -999
       if (Energy > 0) energyVec.at(plane) = Energy;
 
       if (planeNumHits > bestPlaneNumHits) {
         bestPlane = plane;
         bestPlaneNumHits = planeNumHits;
       }
     }
 
     ShowerEleHolder.SetElement(energyVec, energyError, fShowerEnergyOutputLabel);
     // Only set the best plane if it has some hits in it
     if (bestPlane < fGeom->Nplanes()) {
       // Need to cast as an int for legacy default of -999
       // have to define a new variable as we pass-by-reference when filling
       int bestPlaneVal(bestPlane);
       ShowerEleHolder.SetElement(bestPlaneVal, fShowerBestPlaneOutputLabel);
     }
 
     return 0;
   }

double ShowerRecoTools::ShowerNumElectronsEnergy::CalculateEnergy	(	const detinfo::DetectorClocksData &	clockData,
		const detinfo::DetectorPropertiesData &	detProp,
		const std::vector< art::Ptr< recob::Hit >> &	hits,
		const geo::PlaneID::PlaneID_t	plane
	)		const

private

Definition at line 148 of file ShowerNumElectronsEnergy_tool.cc.

   {
 
     double totalCharge = 0;
     double totalEnergy = 0;
     double correctedtotalCharge = 0;
     double nElectrons = 0;
 
     for (auto const& hit : hits) {
       totalCharge +=
         hit->Integral() *
         fCalorimetryAlg.LifetimeCorrection(
           clockData, detProp, hit->PeakTime()); // obtain charge and correct for lifetime
     }
 
     // correct charge due to recombination
     correctedtotalCharge = totalCharge / fRecombinationFactor;
     // calculate # of electrons and the corresponding energy
     nElectrons = fCalorimetryAlg.ElectronsFromADCArea(correctedtotalCharge, plane);
     totalEnergy = (nElectrons / util::kGeVToElectrons) * 1000; // energy in MeV
     return totalEnergy;
   }