#include <ShowerDepth.h>

Public Member Functions
float	getClusterDepthCompatibility (float measuredDepth, float emEnergy, float &expectedDepth, float &expectedSigma) const

	ShowerDepth ()

	~ShowerDepth ()

Static Private Attributes
static constexpr float	corrlnalphalnt0_ { 0.7858 }

static constexpr float	corrlnalphalnt1_ { -0.0232 }

static constexpr float	criticalEnergy_ = 0.00536

static constexpr float	meanalpha0_ { -0.0433 }

static constexpr float	meanalpha1_ { 0.540 }

static constexpr float	meant0_ { -1.396 }

static constexpr float	meant1_ { 1.007 }

static constexpr float	radiationLength_ = 0.968

static constexpr float	sigmalnalpha0_ { -0.08442 }

static constexpr float	sigmalnalpha1_ { 0.7904 }

static constexpr float	sigmalnt0_ { -2.506 }

static constexpr float	sigmalnt1_ { 1.245 }

Detailed Description

Definition at line 27 of file ShowerDepth.h.

Constructor & Destructor Documentation

hgcal::ShowerDepth::ShowerDepth ( )

inline

Definition at line 29 of file ShowerDepth.h.

29 {}

hgcal::ShowerDepth::~ShowerDepth ( )

inline

Definition at line 30 of file ShowerDepth.h.

References getClusterDepthCompatibility().

30 {}

Member Function Documentation

float ShowerDepth::getClusterDepthCompatibility	(	float	measuredDepth,
		float	emEnergy,
		float &	expectedDepth,
		float &	expectedSigma
	)		const

Definition at line 6 of file ShowerDepth.cc.

References corrlnalphalnt0_, corrlnalphalnt1_, criticalEnergy_, cmsBatch::log, meanalpha0_, meanalpha1_, meant0_, meant1_, radiationLength_, sigmalnalpha0_, sigmalnalpha1_, sigmalnt0_, sigmalnt1_, and mathSSE::sqrt().

Referenced by hgcal::EGammaPCAHelper::clusterDepthCompatibility(), and ~ShowerDepth().

                                                                                                                                      {
     float lny = (emEnergy > criticalEnergy_) ? std::log(emEnergy/criticalEnergy_) : 0.;
 
     // inject here parametrization results
     float meantmax = meant0_ + meant1_*lny;
     float meanalpha = meanalpha0_ + meanalpha1_*lny;
     if ( meanalpha < 1. ) meanalpha = 1.1; // no poles
     float sigmalntmax = 1. / (sigmalnt0_+sigmalnt1_*lny);
     float sigmalnalpha = 1. / (sigmalnalpha0_+sigmalnalpha1_*lny);
     float corrlnalphalntmax = corrlnalphalnt0_+corrlnalphalnt1_*lny;
 
     float invbeta = meantmax/(meanalpha-1.);
     float predictedDepth = meanalpha*invbeta;
     predictedDepth *= radiationLength_;
 
     float predictedSigma = sigmalnalpha*sigmalnalpha/((meanalpha-1.)*(meanalpha-1.));
     predictedSigma += sigmalntmax*sigmalntmax;
     predictedSigma -= 2*sigmalnalpha*sigmalntmax*corrlnalphalntmax/(meanalpha-1.);
     if ( predictedSigma < 0. ) predictedSigma = 1.e10; // say we can't predict anything
     predictedSigma = predictedDepth*std::sqrt(predictedSigma);
 
     expectedDepth = predictedDepth;
     expectedSigma = predictedSigma;
     return (measuredDepth-predictedDepth)/predictedSigma;
 }