#include "CalorimetryAlg.h"

Classes
struct	Config

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

	CalorimetryAlg (const Config &config)

	~CalorimetryAlg ()

void	reconfigure (const Config &config)

void	reconfigure (const fhicl::ParameterSet &pset)

double	dEdx_AMP (art::Ptr< recob::Hit > hit, double pitch, double T0=0) const

double	dEdx_AMP (recob::Hit const &hit, double pitch, double T0=0) const

double	dEdx_AMP (double dQ, double time, double pitch, unsigned int plane, double T0=0) const

double	dEdx_AMP (double dQdx, double time, unsigned int plane, double T0=0) const

double	dEdx_AREA (art::Ptr< recob::Hit > hit, double pitch, double T0=0) const

double	dEdx_AREA (recob::Hit const &hit, double pitch, double T0=0) const

double	dEdx_AREA (double dQ, double time, double pitch, unsigned int plane, double T0=0) const

double	dEdx_AREA (double dQdx, double time, unsigned int plane, double T0=0) const

double	ElectronsFromADCPeak (double adc, unsigned short plane) const

double	ElectronsFromADCArea (double area, unsigned short plane) const

double	LifetimeCorrection (double time, double T0=0) const

Private Member Functions
double	dEdx_from_dQdx_e (double dQdx_e, double time, double T0=0) const

Private Attributes
art::ServiceHandle< geo::Geometry >	geom

const detinfo::DetectorProperties *	detprop

std::vector< double >	fCalAmpConstants

std::vector< double >	fCalAreaConstants

bool	fUseModBox

int	fLifeTimeForm

bool	fDoLifeTimeCorrection

Detailed Description

Definition at line 28 of file CalorimetryAlg.h.

Constructor & Destructor Documentation

calo::CalorimetryAlg::CalorimetryAlg ( const fhicl::ParameterSet & pset )

inline

Definition at line 62 of file CalorimetryAlg.h.

                                                       :
                 CalorimetryAlg(fhicl::Table<Config>(pset, {})())
         {}

calo::CalorimetryAlg::CalorimetryAlg ( const Config & config )

Definition at line 26 of file CalorimetryAlg.cxx.

References detprop, and reconfigure().

   {
      this->reconfigure(config);
 
      detprop = art::ServiceHandle<detinfo::DetectorPropertiesService>()->provider();
   }

calo::CalorimetryAlg::~CalorimetryAlg ( )

Definition at line 34 of file CalorimetryAlg.cxx.

   {
     
   }

Member Function Documentation

double calo::CalorimetryAlg::dEdx_AMP	(	art::Ptr< recob::Hit >	hit,
		double	pitch,
		double	T0 = `0`
	)		const

Definition at line 55 of file CalorimetryAlg.cxx.

References recob::Hit::PeakAmplitude(), recob::Hit::PeakTime(), geo::PlaneID::Plane, and recob::Hit::WireID().

Referenced by dEdx_AMP(), trkf::BezierTrack::GetCalorimetryObject(), shwf::ShowerReco::LongTransEnergy(), calo::GeneralCalorimetry::produce(), calo::Calorimetry::produce(), and showerreco::ShowerRecoAlg::RecoOneShower().

   {
     return dEdx_AMP(hit->PeakAmplitude()/pitch, hit->PeakTime(), hit->WireID().Plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AMP	(	recob::Hit const &	hit,
		double	pitch,
		double	T0 = `0`
	)		const

Definition at line 61 of file CalorimetryAlg.cxx.

References dEdx_AMP(), recob::Hit::PeakAmplitude(), recob::Hit::PeakTime(), geo::PlaneID::Plane, and recob::Hit::WireID().

   {
     return dEdx_AMP(hit.PeakAmplitude()/pitch, hit.PeakTime(), hit.WireID().Plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AMP	(	double	dQ,
		double	time,
		double	pitch,
		unsigned int	plane,
		double	T0 = `0`
	)		const

Todo:: The plane argument should really be for a view instead

Definition at line 68 of file CalorimetryAlg.cxx.

References dEdx_AMP().

   {
     double dQdx   = dQ/pitch;           // in ADC/cm
     return dEdx_AMP(dQdx, time, plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AMP	(	double	dQdx,
		double	time,
		unsigned int	plane,
		double	T0 = `0`
	)		const

Definition at line 75 of file CalorimetryAlg.cxx.

References dEdx_from_dQdx_e(), and fCalAmpConstants.

   {
     double fADCtoEl=1.;
     
     fADCtoEl = fCalAmpConstants[plane];
     
     double dQdx_e = dQdx/fADCtoEl;  // Conversion from ADC/cm to e/cm
     return dEdx_from_dQdx_e(dQdx_e,time, T0);
   }

double calo::CalorimetryAlg::dEdx_AREA	(	art::Ptr< recob::Hit >	hit,
		double	pitch,
		double	T0 = `0`
	)		const

Definition at line 88 of file CalorimetryAlg.cxx.

References recob::Hit::Integral(), recob::Hit::PeakTime(), geo::PlaneID::Plane, and recob::Hit::WireID().

Referenced by calo::TrackCalorimetryAlg::AnalyzeHit(), mvapid::MVAAlg::CalcSegmentdEdxDist(), ems::EMShower3D::ConvertFrom(), ems::EMShower3D::ConvertFrom2(), dEdx_AREA(), tca::FilldEdx(), shower::EMShowerAlg::FinddEdx(), shwf::ShowerReco::LongTransEnergy(), shower::EMShowerAlg::MakeShower(), calo::Calorimetry::produce(), and showerreco::ShowerRecoAlg::RecoOneShower().

   {
     return dEdx_AREA(hit->Integral()/pitch, hit->PeakTime(), hit->WireID().Plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AREA	(	recob::Hit const &	hit,
		double	pitch,
		double	T0 = `0`
	)		const

Definition at line 94 of file CalorimetryAlg.cxx.

References dEdx_AREA(), recob::Hit::Integral(), recob::Hit::PeakTime(), geo::PlaneID::Plane, and recob::Hit::WireID().

   {
     return dEdx_AREA(hit.Integral()/pitch, hit.PeakTime(), hit.WireID().Plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AREA	(	double	dQ,
		double	time,
		double	pitch,
		unsigned int	plane,
		double	T0 = `0`
	)		const

Definition at line 100 of file CalorimetryAlg.cxx.

References dEdx_AREA().

   {
     double dQdx   = dQ/pitch;           // in ADC/cm
     return dEdx_AREA(dQdx, time, plane, T0);
   }

double calo::CalorimetryAlg::dEdx_AREA	(	double	dQdx,
		double	time,
		unsigned int	plane,
		double	T0 = `0`
	)		const

Definition at line 107 of file CalorimetryAlg.cxx.

References dEdx_from_dQdx_e(), and fCalAreaConstants.

   {
     double fADCtoEl=1.;
     
     fADCtoEl = fCalAreaConstants[plane];
     
     double dQdx_e = dQdx/fADCtoEl;  // Conversion from ADC/cm to e/cm
     return dEdx_from_dQdx_e(dQdx_e, time, T0);
   }

double calo::CalorimetryAlg::dEdx_from_dQdx_e	(	double	dQdx_e,
		double	time,
		double	T0 = `0`
	)		const

private

Definition at line 118 of file CalorimetryAlg.cxx.

References detinfo::DetectorProperties::BirksCorrection(), detprop, fDoLifeTimeCorrection, fUseModBox, LifetimeCorrection(), and detinfo::DetectorProperties::ModBoxCorrection().

Referenced by dEdx_AMP(), and dEdx_AREA().

   {
     if (fDoLifeTimeCorrection)
       dQdx_e *= LifetimeCorrection(time, T0);   // Lifetime Correction (dQdx_e in e/cm)
     if(fUseModBox) {
       return detprop->ModBoxCorrection(dQdx_e);
     } else {
       return detprop->BirksCorrection(dQdx_e);
     }
   }

double calo::CalorimetryAlg::ElectronsFromADCArea	(	double	area,
		unsigned short	plane
	)		const

inline

Definition at line 87 of file CalorimetryAlg.h.

Referenced by showerreco::ShowerRecoAlg::RecoOneShower().

88 { return area / fCalAreaConstants[plane]; }

calo::CalorimetryAlg::fCalAreaConstants

std::vector< double > fCalAreaConstants

Definition: CalorimetryAlg.h:100

double calo::CalorimetryAlg::ElectronsFromADCPeak	(	double	adc,
		unsigned short	plane
	)		const

inline

Definition at line 84 of file CalorimetryAlg.h.

Referenced by img::DataProviderAlg::DataProviderAlg(), and showerreco::ShowerRecoAlg::RecoOneShower().

85 { return adc / fCalAmpConstants[plane]; }

calo::CalorimetryAlg::fCalAmpConstants

std::vector< double > fCalAmpConstants

Definition: CalorimetryAlg.h:99

double calo::CalorimetryAlg::LifetimeCorrection	(	double	time,
		double	T0 = `0`
	)		const

Definition at line 133 of file CalorimetryAlg.cxx.

References detprop, e, detinfo::DetectorProperties::ElectronLifetime(), fLifeTimeForm, detinfo::DetectorProperties::SamplingRate(), and detinfo::DetectorProperties::TriggerOffset().

Referenced by dEdx_from_dQdx_e(), img::DataProviderAlg::LifetimeCorrection(), showerreco::ShowerRecoAlg::RecoOneShower(), and img::DataProviderAlg::resizeView().

   {  
     float t = time;
 
     double timetick = detprop->SamplingRate()*1.e-3;    //time sample in microsec
     double presamplings = detprop->TriggerOffset();
     
     t -= presamplings;
     time = t * timetick - T0*1e-3;  //  (in microsec)
 
     if (fLifeTimeForm==0){
       //Exponential form
       double tau = detprop->ElectronLifetime();
     
       double correction = exp(time/tau);
       return correction;
     }
     else if (fLifeTimeForm==1){
       //Exponential+constant form
       const lariov::ElectronLifetimeProvider& elifetime_provider = art::ServiceHandle<lariov::ElectronLifetimeService>()->GetProvider();
       double correction = elifetime_provider.Lifetime(time);
       //std::cout<<correction<<std::endl;
       return correction;
     }
     else{
       throw cet::exception("CalorimetryAlg") << "Unknow CaloLifeTimeForm "<<fLifeTimeForm<<std::endl;
     }
   }

void calo::CalorimetryAlg::reconfigure ( const Config & config )

Definition at line 40 of file CalorimetryAlg.cxx.

References calo::CalorimetryAlg::Config::CalAmpConstants, calo::CalorimetryAlg::Config::CalAreaConstants, calo::CalorimetryAlg::Config::CaloDoLifeTimeCorrection, calo::CalorimetryAlg::Config::CaloLifeTimeForm, calo::CalorimetryAlg::Config::CaloUseModBox, fCalAmpConstants, fCalAreaConstants, fDoLifeTimeCorrection, fLifeTimeForm, and fUseModBox.

Referenced by CalorimetryAlg(), img::DataProviderAlg::DataProviderAlg(), and calo::TrackCalorimetryAlg::reconfigure().

   {
     
     fCalAmpConstants    = config.CalAmpConstants();
     fCalAreaConstants   = config.CalAreaConstants();
     fUseModBox          = config.CaloUseModBox();
     fLifeTimeForm       = config.CaloLifeTimeForm();
     fDoLifeTimeCorrection = config.CaloDoLifeTimeCorrection();
 
     return;
   }