#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

◆ ShowerDepth()

hgcal::ShowerDepth::ShowerDepth ( )

inline

Definition at line 29 of file ShowerDepth.h.

29 {}

◆ ~ShowerDepth()

hgcal::ShowerDepth::~ShowerDepth ( )

inline

Definition at line 30 of file ShowerDepth.h.

30 {}

Member Function Documentation

◆ getClusterDepthCompatibility()

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

Definition at line 6 of file ShowerDepth.cc.

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