#include "ShowerEnergyAlg.h"

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

double	ShowerEnergy (std::vector< art::Ptr< recob::Hit > > const &hits, int plane)
	Finds the total energy deposited by the shower in this view. More...

Private Attributes
double	fUGradient

double	fUIntercept

double	fVGradient

double	fVIntercept

double	fZGradient

double	fZIntercept

detinfo::DetectorProperties const *	detprop = nullptr

Detailed Description

Definition at line 30 of file ShowerEnergyAlg.h.

Constructor & Destructor Documentation

shower::ShowerEnergyAlg::ShowerEnergyAlg ( fhicl::ParameterSet const & pset )

Definition at line 11 of file ShowerEnergyAlg.cxx.

References fUGradient, fUIntercept, fVGradient, fVIntercept, fZGradient, fZIntercept, and fhicl::ParameterSet::get().

   : detprop(lar::providerFrom<detinfo::DetectorPropertiesService>())
 {
   fUGradient  = pset.get<double>("UGradient");
   fUIntercept = pset.get<double>("UIntercept");
   fVGradient  = pset.get<double>("VGradient");
   fVIntercept = pset.get<double>("VIntercept");
   fZGradient  = pset.get<double>("ZGradient");
   fZIntercept = pset.get<double>("ZIntercept");
 }

Member Function Documentation

double shower::ShowerEnergyAlg::ShowerEnergy	(	std::vector< art::Ptr< recob::Hit > > const &	hits,
		int	plane
	)

Finds the total energy deposited by the shower in this view.

Definition at line 22 of file ShowerEnergyAlg.cxx.

References detprop, detinfo::DetectorProperties::ElectronLifetime(), fUGradient, fUIntercept, fVGradient, fVIntercept, fZGradient, fZIntercept, hits(), and detinfo::DetectorProperties::SamplingRate().

Referenced by shower::EMShowerAlg::MakeShower(), and ems::MergeEMShower3D::produce().

                                                                                               {
 
   double totalCharge = 0, totalEnergy = 0;
 
   for (art::PtrVector<recob::Hit>::const_iterator hit = hits.begin(); hit != hits.end(); ++hit){
     if (int((*hit)->WireID().Plane)!=plane) continue;
     totalCharge += ( (*hit)->Integral() * TMath::Exp( (detprop->SamplingRate() * (*hit)->PeakTime()) / (detprop->ElectronLifetime()*1e3) ) );
   }
 
   switch (plane) {
   case 0:
     totalEnergy = (totalCharge * fUGradient) + fUIntercept;
     break;
   case 1:
     totalEnergy = (totalCharge * fVGradient) + fVIntercept;
     break;
   case 2:
     totalEnergy = (totalCharge * fZGradient) + fZIntercept;
     break;
   }
 
   return totalEnergy;
 
 }