EMECALShowerParametrization Class Reference

#include <EMECALShowerParametrization.h>

Public Member Functions
double	correlationAlphaT (double lny) const
double	correlationAlphaTHom (double lny) const
double	correlationAlphaTSam (double lny) const
const ECALProperties *	ecalProperties () const
	EMECALShowerParametrization (const ECALProperties *ecal, const HCALProperties *hcal, const PreshowerLayer1Properties *layer1, const PreshowerLayer2Properties *layer2, const std::vector< double > &coreIntervals, const std::vector< double > &tailIntervals, double RCFact=1., double RTFact=1.)
const std::vector< double > &	getCoreIntervals () const
const std::vector< double > &	getTailIntervals () const
const HCALProperties *	hcalProperties () const
const PreshowerLayer1Properties *	layer1Properties () const
const PreshowerLayer2Properties *	layer2Properties () const
double	meanAlpha (double lny) const
double	meanAlphaHom (double lny) const
double	meanAlphaSam (double lny) const
double	meanAlphaSpot (double alpha) const
double	meanAlphaSpotHom (double alpha) const
double	meanAlphaSpotSam (double alpha) const
double	meanLnAlpha (double lny) const
double	meanLnAlphaHom (double lny) const
double	meanLnAlphaSam (double lny) const
double	meanLnT (double lny) const
double	meanLnTHom (double lny) const
double	meanLnTSam (double lny) const
double	meanT (double lny) const
double	meanTHom (double lny) const
double	meanTSam (double lny) const
double	meanTSpot (double T) const
double	meanTSpotHom (double T) const
double	meanTSpotSam (double T) const
double	nSpots (double E) const
double	nSpotsHom (double E) const
double	nSpotsSam (double E) const
double	p (double tau, double E) const
double	pHom (double tau, double E) const
double	pSam (double tau, double E) const
double	rC (double tau, double E) const
double	rCHom (double tau, double E) const
double	rCSam (double tau, double E) const
double	rT (double tau, double E) const
double	rTHom (double tau, double E) const
double	rTSam (double tau, double E) const
double	sigmaLnAlpha (double lny) const
double	sigmaLnAlphaHom (double lny) const
double	sigmaLnAlphaSam (double lny) const
double	sigmaLnT (double lny) const
double	sigmaLnTHom (double lny) const
double	sigmaLnTSam (double lny) const
virtual	~EMECALShowerParametrization ()

Private Member Functions
double	k1 () const
double	k2 () const
double	k3 () const
double	k4 (double E) const
double	p1 () const
double	p2 () const
double	p3 (double E) const
double	z1 (double E) const
double	z2 () const

Private Attributes
const std::vector< double > &	theCore
const ECALProperties *	theECAL
const HCALProperties *	theHCAL
const PreshowerLayer1Properties *	theLayer1
const PreshowerLayer2Properties *	theLayer2
double	theRcfactor
double	theRtfactor
const std::vector< double > &	theTail

Detailed Description

Electromagnetic Shower parametrization utilities according to G. Grindhammer and S. Peters, hep-ex/0001020, Appendix A

Author

Patrick Janot

Date

: 25-Jan-2004

Definition at line 19 of file EMECALShowerParametrization.h.

Constructor & Destructor Documentation

EMECALShowerParametrization::EMECALShowerParametrization ( const ECALProperties *  ecal, const HCALProperties *  hcal, const PreshowerLayer1Properties *  layer1, const PreshowerLayer2Properties *  layer2, const std::vector< double > &  coreIntervals, const std::vector< double > &  tailIntervals, double  RCFact = 1., double  RTFact = 1.  )

inline

Definition at line 23 of file EMECALShowerParametrization.h.

                                    : 
    theECAL(ecal),  
    theHCAL(hcal),
    theLayer1(layer1), 
    theLayer2(layer2),
    theCore(coreIntervals),
    theTail(tailIntervals),
    theRcfactor(RCFact),
    theRtfactor(RTFact){}

virtual EMECALShowerParametrization::~EMECALShowerParametrization ( )

inlinevirtual

Definition at line 40 of file EMECALShowerParametrization.h.

{ }

Member Function Documentation

double EMECALShowerParametrization::correlationAlphaT ( double  lny ) const

inline

Definition at line 101 of file EMECALShowerParametrization.h.

References correlationAlphaTHom(), correlationAlphaTSam(), ECALProperties::isHom(), and theECAL.

Referenced by EMShower::EMShower().

                                                    {
    if (theECAL->isHom()) return correlationAlphaTHom(lny); 
    return correlationAlphaTSam(lny); }

double EMECALShowerParametrization::correlationAlphaTHom ( double  lny ) const

inline

Definition at line 123 of file EMECALShowerParametrization.h.

Referenced by correlationAlphaT().

                                                       {
    return 0.705-0.023*lny; }

double EMECALShowerParametrization::correlationAlphaTSam ( double  lny ) const

inline

Definition at line 141 of file EMECALShowerParametrization.h.

Referenced by correlationAlphaT().

                                                       {
    return 0.784-0.023*lny; }

const ECALProperties* EMECALShowerParametrization::ecalProperties ( ) const

inline

Definition at line 248 of file EMECALShowerParametrization.h.

References theECAL.

Referenced by EMShower::EMShower().

                                                      { 
    return theECAL; 
  }

const std::vector<double>& EMECALShowerParametrization::getCoreIntervals ( ) const

inline

Definition at line 264 of file EMECALShowerParametrization.h.

References theCore.

Referenced by EMShower::setIntervals().

{ return theCore;}

const std::vector<double>& EMECALShowerParametrization::getTailIntervals ( ) const

inline

Definition at line 266 of file EMECALShowerParametrization.h.

References theTail.

Referenced by EMShower::setIntervals().

{ return theTail;}

const HCALProperties* EMECALShowerParametrization::hcalProperties ( ) const

inline

Definition at line 252 of file EMECALShowerParametrization.h.

References theHCAL.

Referenced by EMShower::EMShower().

                                                      { 
    return theHCAL; 
  }

double EMECALShowerParametrization::k1 ( ) const

inlineprivate

Definition at line 288 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by rTHom().

{ return 0.6590-0.00309*theECAL->theZeff(); }

double EMECALShowerParametrization::k2 ( ) const

inlineprivate

Definition at line 289 of file EMECALShowerParametrization.h.

Referenced by rTHom().

{ return 0.6450; }

double EMECALShowerParametrization::k3 ( ) const

inlineprivate

Definition at line 290 of file EMECALShowerParametrization.h.

Referenced by rTHom().

{ return -2.59; }

double EMECALShowerParametrization::k4 ( double  E ) const

inlineprivate

Definition at line 291 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log.

Referenced by rTHom().

{ return 0.3585+0.0421*std::log(E); }

const PreshowerLayer1Properties* EMECALShowerParametrization::layer1Properties ( ) const

inline

Definition at line 256 of file EMECALShowerParametrization.h.

References theLayer1.

Referenced by EMShower::EMShower().

                                                                   { 
    return theLayer1; 
  }

const PreshowerLayer2Properties* EMECALShowerParametrization::layer2Properties ( ) const

inline

Definition at line 260 of file EMECALShowerParametrization.h.

References theLayer2.

Referenced by EMShower::EMShower().

                                                                   { 
    return theLayer2; 
  }

double EMECALShowerParametrization::meanAlpha ( double  lny ) const

inline

Definition at line 51 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanAlphaHom(), meanAlphaSam(), and theECAL.

Referenced by EMShower::EMShower().

                                            { 
    if (theECAL->isHom()) return meanAlphaHom(lny);
    return meanAlphaSam(lny); }

double EMECALShowerParametrization::meanAlphaHom ( double  lny ) const

inline

Definition at line 61 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by meanAlpha(), and meanAlphaSam().

                                               { 
    return 0.21+(0.492+2.38/theECAL->theZeff())*lny; }

double EMECALShowerParametrization::meanAlphaSam ( double  lny ) const

inline

Definition at line 71 of file EMECALShowerParametrization.h.

References meanAlphaHom(), theECAL, and ECALProperties::theFs().

Referenced by meanAlpha().

                                               { 
    return meanAlphaHom(lny) - 0.444/theECAL->theFs(); }

double EMECALShowerParametrization::meanAlphaSpot ( double  alpha ) const

inline

Definition at line 214 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanAlphaSpotHom(), meanAlphaSpotSam(), and theECAL.

Referenced by EMShower::EMShower().

                                                  {
    if (theECAL->isHom()) return meanAlphaSpotHom(alpha); 
    return meanAlphaSpotSam(alpha); }

double EMECALShowerParametrization::meanAlphaSpotHom ( double  alpha ) const

inline

Definition at line 229 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by meanAlphaSpot(), and meanAlphaSpotSam().

                                                     {
    return alpha*(0.639+0.00334*theECAL->theZeff()); }

double EMECALShowerParametrization::meanAlphaSpotSam ( double  alpha ) const

inline

Definition at line 241 of file EMECALShowerParametrization.h.

References meanAlphaSpotHom(), theECAL, and ECALProperties::theZeff().

Referenced by meanAlphaSpot().

                                                     {
    return meanAlphaSpotHom(alpha)*(0.844+0.0026*theECAL->theZeff()); }

double EMECALShowerParametrization::meanLnAlpha ( double  lny ) const

inline

Definition at line 91 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanLnAlphaHom(), meanLnAlphaSam(), and theECAL.

Referenced by EMShower::EMShower().

                                              {
    if (theECAL->isHom()) return meanLnAlphaHom(lny); 
    return meanLnAlphaSam(lny); }

double EMECALShowerParametrization::meanLnAlphaHom ( double  lny ) const

inline

Definition at line 117 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log, theECAL, and ECALProperties::theZeff().

Referenced by meanLnAlpha(), and meanLnAlphaSam().

                                                 {
    return std::log(0.81+(0.458+2.26/theECAL->theZeff())*lny); }

double EMECALShowerParametrization::meanLnAlphaSam ( double  lny ) const

inline

Definition at line 135 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::exp, create_public_lumi_plots::log, meanLnAlphaHom(), theECAL, and ECALProperties::theZeff().

Referenced by meanLnAlpha().

                                                 {
    return log( std::exp(meanLnAlphaHom(lny)) - 0.476/theECAL->theZeff() ); }

double EMECALShowerParametrization::meanLnT ( double  lny ) const

inline

Definition at line 82 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanLnTHom(), meanLnTSam(), and theECAL.

Referenced by EMShower::EMShower().

                                          {
    if (theECAL->isHom()) return meanLnTHom(lny); 
    return meanLnTSam(lny); }

double EMECALShowerParametrization::meanLnTHom ( double  lny ) const

inline

Definition at line 110 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log.

Referenced by meanLnT(), and meanLnTSam().

                                             {
    return std::log(lny-0.812); }

double EMECALShowerParametrization::meanLnTSam ( double  lny ) const

inline

Definition at line 129 of file EMECALShowerParametrization.h.

References ECALProperties::ehat(), create_public_lumi_plots::exp, create_public_lumi_plots::log, meanLnTHom(), theECAL, and ECALProperties::theFs().

Referenced by meanLnT().

                                             {
    return log( std::exp(meanLnTHom(lny)) - 0.55/theECAL->theFs() - 0.69*(1-theECAL->ehat()) ); }

double EMECALShowerParametrization::meanT ( double  lny ) const

inline

Definition at line 47 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanTHom(), meanTSam(), and theECAL.

Referenced by EMShower::EMShower().

                                        { 
    if (theECAL->isHom()) return meanTHom(lny);
    return meanTSam(lny); }

double EMECALShowerParametrization::meanTHom ( double  lny ) const

inline

Definition at line 58 of file EMECALShowerParametrization.h.

Referenced by meanT(), and meanTSam().

                                           { 
    return lny-0.858; }

double EMECALShowerParametrization::meanTSam ( double  lny ) const

inline

Definition at line 68 of file EMECALShowerParametrization.h.

References ECALProperties::ehat(), meanTHom(), theECAL, and ECALProperties::theFs().

Referenced by meanT().

                                           { 
    return meanTHom(lny) - 0.59/theECAL->theFs() - 0.53*(1.-theECAL->ehat()); }

double EMECALShowerParametrization::meanTSpot ( double  T ) const

inline

Definition at line 210 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), meanTSpotHom(), meanTSpotSam(), and theECAL.

Referenced by EMShower::EMShower().

                                          {
    if (theECAL->isHom()) return meanTSpotHom(T); 
    return meanTSpotSam(T); }

double EMECALShowerParametrization::meanTSpotHom ( double  T ) const

inline

Definition at line 226 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by meanTSpot(), and meanTSpotSam().

                                             {
    return T*(0.698+0.00212*theECAL->theZeff()); }

double EMECALShowerParametrization::meanTSpotSam ( double  T ) const

inline

Definition at line 238 of file EMECALShowerParametrization.h.

References meanTSpotHom(), theECAL, and ECALProperties::theZeff().

Referenced by meanTSpot().

                                             {
    return meanTSpotHom(T)*(0.813+0.0019*theECAL->theZeff()); }

double EMECALShowerParametrization::nSpots ( double  E ) const

inline

Definition at line 206 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), nSpotsHom(), nSpotsSam(), and theECAL.

Referenced by EMShower::EMShower().

                                       {
    if (theECAL->isHom()) return nSpotsHom(E); 
    return nSpotsSam(E); }

double EMECALShowerParametrization::nSpotsHom ( double  E ) const

inline

Definition at line 223 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log, funct::pow(), theECAL, and ECALProperties::theZeff().

Referenced by nSpots().

                                          {
    return 93.*std::log(theECAL->theZeff()) * std::pow(E,0.876); }

double EMECALShowerParametrization::nSpotsSam ( double  E ) const

inline

Definition at line 235 of file EMECALShowerParametrization.h.

References funct::pow(), ECALProperties::resE(), and theECAL.

Referenced by nSpots().

                                          {
    return 10.3/theECAL->resE()*std::pow(E, 0.959); }

double EMECALShowerParametrization::p ( double  tau, double  E  ) const

inline

Definition at line 162 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), pHom(), pSam(), and theECAL.

Referenced by EMShower::compute().

                                              {
    if (theECAL->isHom()) return pHom(tau, E); 
    return pSam(tau, E); }

double EMECALShowerParametrization::p1 ( ) const

inlineprivate

Definition at line 281 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by pHom().

{ return 2.632-0.00094*theECAL->theZeff(); }

double EMECALShowerParametrization::p2 ( ) const

inlineprivate

Definition at line 282 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by pHom().

{ return 0.401+0.00187*theECAL->theZeff(); }

double EMECALShowerParametrization::p3 ( double  E ) const

inlineprivate

Definition at line 283 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log.

Referenced by pHom().

{ return 1.313-0.0686*std::log(E); }

double EMECALShowerParametrization::pHom ( double  tau, double  E  ) const

inline

Definition at line 179 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::exp, p1(), p2(), p3(), and metsig::tau.

Referenced by p(), and pSam().

                                                 {
    double arg = (p2()-tau)/p3(E);
    return p1()* std::exp(arg-std::exp(arg));
  }

double EMECALShowerParametrization::pSam ( double  tau, double  E  ) const

inline

Definition at line 195 of file EMECALShowerParametrization.h.

References ECALProperties::ehat(), create_public_lumi_plots::exp, pHom(), funct::pow(), theECAL, and ECALProperties::theFs().

Referenced by p().

                                                 {
    return pHom(tau, E) + (1-theECAL->ehat())*(0.348-0.642/theECAL->theFs()*std::exp(-1.*std::pow((tau-1),2) ) );
  }

double EMECALShowerParametrization::rC ( double  tau, double  E  ) const

inline

Definition at line 154 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), rCHom(), rCSam(), and theECAL.

Referenced by EMShower::compute().

                                               {
    if (theECAL->isHom()) return rCHom(tau, E); 
    return rCSam(tau, E); }

double EMECALShowerParametrization::rCHom ( double  tau, double  E  ) const

inline

Definition at line 170 of file EMECALShowerParametrization.h.

References metsig::tau, theRcfactor, z1(), and z2().

Referenced by rC(), and rCSam().

                                                  {
    return theRcfactor*(z1(E) + z2()*tau);
  }  

double EMECALShowerParametrization::rCSam ( double  tau, double  E  ) const

inline

Definition at line 187 of file EMECALShowerParametrization.h.

References ECALProperties::ehat(), create_public_lumi_plots::exp, rCHom(), theECAL, and ECALProperties::theFs().

Referenced by rC().

                                                  {
    return rCHom(tau, E) - 0.0203*(1-theECAL->ehat()) + 0.0397/theECAL->theFs()*std::exp(-1.*tau);
  }

double EMECALShowerParametrization::rT ( double  tau, double  E  ) const

inline

Definition at line 158 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), rTHom(), rTSam(), and theECAL.

Referenced by EMShower::compute().

                                               {
    if (theECAL->isHom()) return rTHom(tau, E); 
    return rTSam(tau, E); }

double EMECALShowerParametrization::rTHom ( double  tau, double  E  ) const

inline

Definition at line 174 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::exp, k1(), k2(), k3(), k4(), and theRtfactor.

Referenced by rT(), and rTSam().

                                                 {
    return theRtfactor*k1() * ( std::exp(k3()*(tau-k2()))+
                std::exp(k4(E)*(tau-k2())) );
  }

double EMECALShowerParametrization::rTSam ( double  tau, double  E  ) const

inline

Definition at line 191 of file EMECALShowerParametrization.h.

References ECALProperties::ehat(), create_public_lumi_plots::exp, rTHom(), theECAL, and ECALProperties::theFs().

Referenced by rT().

                                                 {
    return rTHom(tau, E) -0.14*(1-theECAL->ehat()) - 0.495/theECAL->theFs()*std::exp(-1.*tau);
  }

double EMECALShowerParametrization::sigmaLnAlpha ( double  lny ) const

inline

Definition at line 96 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), sigmaLnAlphaHom(), sigmaLnAlphaSam(), and theECAL.

Referenced by EMShower::EMShower().

                                               {
    if (theECAL->isHom()) return sigmaLnAlphaHom(lny); 
    return sigmaLnAlphaSam(lny); }

double EMECALShowerParametrization::sigmaLnAlphaHom ( double  lny ) const

inline

Definition at line 120 of file EMECALShowerParametrization.h.

Referenced by sigmaLnAlpha().

                                                  {
    return 1./(-0.58+0.86*lny); }

double EMECALShowerParametrization::sigmaLnAlphaSam ( double  lny ) const

inline

Definition at line 138 of file EMECALShowerParametrization.h.

Referenced by sigmaLnAlpha().

                                                  {
    return 1./(-0.82+0.79*lny); }

double EMECALShowerParametrization::sigmaLnT ( double  lny ) const

inline

Definition at line 86 of file EMECALShowerParametrization.h.

References ECALProperties::isHom(), sigmaLnTHom(), sigmaLnTSam(), and theECAL.

Referenced by EMShower::EMShower().

                                           {
    if (theECAL->isHom()) return sigmaLnTHom(lny); 
    return sigmaLnTSam(lny); }

double EMECALShowerParametrization::sigmaLnTHom ( double  lny ) const

inline

Definition at line 114 of file EMECALShowerParametrization.h.

Referenced by sigmaLnT().

                                              {
    return 1./(-1.4+1.26*lny); }

double EMECALShowerParametrization::sigmaLnTSam ( double  lny ) const

inline

Definition at line 132 of file EMECALShowerParametrization.h.

Referenced by sigmaLnT().

                                              {
    return 1/(-2.5 + 1.25*lny); }

double EMECALShowerParametrization::z1 ( double  E ) const

inlineprivate

Definition at line 285 of file EMECALShowerParametrization.h.

References create_public_lumi_plots::log.

Referenced by rCHom().

{ return 0.0251+0.00319*std::log(E); }

double EMECALShowerParametrization::z2 ( ) const

inlineprivate

Definition at line 286 of file EMECALShowerParametrization.h.

References theECAL, and ECALProperties::theZeff().

Referenced by rCHom().

{ return 0.1162-0.000381*theECAL->theZeff(); }

Member Data Documentation

const std::vector<double>& EMECALShowerParametrization::theCore

private

Definition at line 275 of file EMECALShowerParametrization.h.

Referenced by getCoreIntervals().

const ECALProperties* EMECALShowerParametrization::theECAL

private

Definition at line 270 of file EMECALShowerParametrization.h.

Referenced by correlationAlphaT(), ecalProperties(), k1(), meanAlpha(), meanAlphaHom(), meanAlphaSam(), meanAlphaSpot(), meanAlphaSpotHom(), meanAlphaSpotSam(), meanLnAlpha(), meanLnAlphaHom(), meanLnAlphaSam(), meanLnT(), meanLnTSam(), meanT(), meanTSam(), meanTSpot(), meanTSpotHom(), meanTSpotSam(), nSpots(), nSpotsHom(), nSpotsSam(), p(), p1(), p2(), pSam(), rC(), rCSam(), rT(), rTSam(), sigmaLnAlpha(), sigmaLnT(), and z2().

const HCALProperties* EMECALShowerParametrization::theHCAL

private

Definition at line 271 of file EMECALShowerParametrization.h.

Referenced by hcalProperties().

const PreshowerLayer1Properties* EMECALShowerParametrization::theLayer1

private

Definition at line 272 of file EMECALShowerParametrization.h.

Referenced by layer1Properties().

const PreshowerLayer2Properties* EMECALShowerParametrization::theLayer2

private

Definition at line 273 of file EMECALShowerParametrization.h.

Referenced by layer2Properties().

double EMECALShowerParametrization::theRcfactor

private

Definition at line 278 of file EMECALShowerParametrization.h.

Referenced by rCHom().

double EMECALShowerParametrization::theRtfactor

private

Definition at line 279 of file EMECALShowerParametrization.h.

Referenced by rTHom().

const std::vector<double>& EMECALShowerParametrization::theTail

private

Definition at line 276 of file EMECALShowerParametrization.h.

Referenced by getTailIntervals().