CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions
GflashKaonMinusShowerProfile Class Reference

#include <GflashKaonMinusShowerProfile.h>

Inheritance diagram for GflashKaonMinusShowerProfile:
GflashHadronShowerProfile

Public Member Functions

 GflashKaonMinusShowerProfile (edm::ParameterSet parSet)
 
void loadParameters ()
 
 ~GflashKaonMinusShowerProfile ()
 
- Public Member Functions inherited from GflashHadronShowerProfile
std::vector< GflashHit > & getGflashHitList ()
 
GflashShowino * getGflashShowino ()
 
 GflashHadronShowerProfile (edm::ParameterSet parSet)
 
void hadronicParameterization ()
 
void initialize (int showerType, double energy, double globalTime, double charge, Gflash3Vector &position, Gflash3Vector &momentum)
 
virtual ~GflashHadronShowerProfile ()
 

Additional Inherited Members

- Protected Member Functions inherited from GflashHadronShowerProfile
double depthScale (double ssp, double ssp0, double length)
 
void doCholeskyReduction (double **cc, double **vv, const int ndim)
 
double fLnE1 (double einc, const double *par)
 
double fTanh (double einc, const double *par)
 
double gammaProfile (double alpha, double beta, double depth, double lengthUnit)
 
void getFluctuationVector (double *lowTriangle, double *correlationVector)
 
int getNumberOfSpots (Gflash::CalorimeterNumber kCalor)
 
double hoProfile (double pathLength, double refDepth)
 
Gflash3Vector locateHitPosition (GflashTrajectoryPoint &point, double lateralArm)
 
double longitudinalProfile ()
 
double medianLateralArm (double depth, Gflash::CalorimeterNumber kCalor)
 
void setEnergyScale (double einc, Gflash3Vector ssp)
 
double twoGammaProfile (double *par, double depth, Gflash::CalorimeterNumber kIndex)
 
- Protected Attributes inherited from GflashHadronShowerProfile
double averageSpotEnergy [Gflash::kNumberCalorimeter]
 
double energyScale [Gflash::kNumberCalorimeter]
 
double lateralPar [Gflash::kNumberCalorimeter][Gflash::Nrpar]
 
double longEcal [Gflash::NPar]
 
double longHcal [Gflash::NPar]
 
double theBField
 
bool theGflashHcalOuter
 
std::vector< GflashHit > theGflashHitList
 
GflashHistogram * theHisto
 
edm::ParameterSet theParSet
 
GflashShowino * theShowino
 

Detailed Description

Definition at line 6 of file GflashKaonMinusShowerProfile.h.

Constructor & Destructor Documentation

GflashKaonMinusShowerProfile::GflashKaonMinusShowerProfile ( edm::ParameterSet  parSet)
inline

Definition at line 12 of file GflashKaonMinusShowerProfile.h.

12  :
13  GflashHadronShowerProfile (parSet) {};
GflashHadronShowerProfile::GflashHadronShowerProfile
GflashHadronShowerProfile(edm::ParameterSet parSet)
Definition: GflashHadronShowerProfile.cc:17
GflashKaonMinusShowerProfile::~GflashKaonMinusShowerProfile ( )
inline

Definition at line 14 of file GflashKaonMinusShowerProfile.h.

14 {};

Member Function Documentation

void GflashKaonMinusShowerProfile::loadParameters ( )
virtual

Reimplemented from GflashHadronShowerProfile.

Definition at line 4 of file GflashKaonMinusShowerProfile.cc.

References abs, Gflash::correl_hadem, GflashHadronShowerProfile::depthScale(), Gflash::emcorr_pid, Gflash::emscale, GflashHadronShowerProfile::energyScale, funct::exp(), GflashHadronShowerProfile::fLnE1(), GflashHadronShowerProfile::fTanh(), GflashShowino::getEnergy(), GflashHadronShowerProfile::getFluctuationVector(), GflashShowino::getPositionAtShower(), GflashShowino::getShowerType(), Gflash::hadcorr_pid, Gflash::hadscale, i, j, Gflash::kENCA, Gflash::kESPM, Gflash::kHB, Gflash::kHE, GflashHadronShowerProfile::lateralPar, Gflash::LengthCrystalEB, Gflash::LengthCrystalEE, GflashHadronShowerProfile::longEcal, GflashHadronShowerProfile::longHcal, max(), Gflash::mipcorr_pid, Gflash::NPar, Gflash::Nrpar, Gflash::par, position, Gflash::RFrontCrystalEB, Gflash::rho, Gflash::Rmin, Gflash::rpar, mathSSE::sqrt(), GflashHadronShowerProfile::theShowino, Gflash::ZFrontCrystalEE, and Gflash::Zmin.

5 {
6  double einc = theShowino->getEnergy();
7  Gflash3Vector position = theShowino->getPositionAtShower();
8  int showerType = theShowino->getShowerType();
9 
10  // energy scale
11  double energyMeanHcal = 0.0;
12  double energySigmaHcal = 0.0;
13 
14  if(showerType == 0 || showerType == 1 || showerType == 4 || showerType == 5) {
15 
16  double r1 = 0.0;
17  double r2 = 0.0;
18 
19  //@@@ energy dependent energyRho based on tuning with testbeam data
20  double energyRho = fTanh(einc,Gflash::correl_hadem);
21 
22  r1 = CLHEP::RandGaussQ::shoot();
23 
24  if(showerType == 0 || showerType == 1) {
25  energyScale[Gflash::kESPM] = einc*(fTanh(einc,Gflash::emscale[0]) + fTanh(einc,Gflash::emcorr_pid[2])
26  +(fTanh(einc,Gflash::emscale[2]) + fTanh(einc,Gflash::emscale[3])*depthScale(position.getRho(),151.,22.) )*r1);
27 
28  energyMeanHcal = (fTanh(einc,Gflash::hadscale[0]) +
29  fTanh(einc,Gflash::hadscale[1])*depthScale(position.getRho(),Gflash::RFrontCrystalEB,Gflash::LengthCrystalEB));
30  energySigmaHcal = (fTanh(einc,Gflash::hadscale[2]) +
31  fTanh(einc,Gflash::hadscale[3])*depthScale(position.getRho(),Gflash::RFrontCrystalEB,Gflash::LengthCrystalEB));
32 
33  r2 = CLHEP::RandGaussQ::shoot();
34  energyScale[Gflash::kHB] = std::max(0.0,
35  einc*fTanh(einc,Gflash::hadcorr_pid[2]) +
36  exp(energyMeanHcal+energySigmaHcal*(energyRho*r1 + sqrt(1.0- energyRho*energyRho)*r2 ))-0.05*einc);
37  }
38  else {
39  energyScale[Gflash::kENCA] = einc*(fTanh(einc,Gflash::emscale[0]) + fTanh(einc,Gflash::emcorr_pid[3])
40  +(fTanh(einc,Gflash::emscale[2]) + fTanh(einc,Gflash::emscale[3])*depthScale(std::fabs(position.getZ()),338.0,21.) )*r1);
41 
42  //@@@extend depthScale for HE
43  energyMeanHcal = (fTanh(einc,Gflash::hadscale[0]) +
44  fTanh(einc,Gflash::hadscale[1])*depthScale(std::fabs(position.getZ()),Gflash::ZFrontCrystalEE,Gflash::LengthCrystalEE));
45  energySigmaHcal = (fTanh(einc,Gflash::hadscale[2]) +
46  fTanh(einc,Gflash::hadscale[3])*depthScale(std::fabs(position.getZ()),Gflash::ZFrontCrystalEE,Gflash::LengthCrystalEE));
47  r2 = CLHEP::RandGaussQ::shoot();
48  energyScale[Gflash::kHE] = std::max(0.0,
49  einc*fTanh(einc,Gflash::hadcorr_pid[3]) +
50  exp(energyMeanHcal+energySigmaHcal*(energyRho*r1 + sqrt(1.0- energyRho*energyRho)*r2 ))-0.05*einc);
51  }
52  }
53  else if(showerType == 2 || showerType == 6 || showerType == 3 || showerType == 7) {
54  //Hcal response for mip-like pions (mip)
55 
56  energyMeanHcal = fTanh(einc,Gflash::hadscale[4]);
57  energySigmaHcal = fTanh(einc,Gflash::hadscale[5]);
58  double gap_corr = einc * fTanh(einc,Gflash::hadscale[6]);
59 
60  if(showerType == 2 || showerType == 3) {
61  energyScale[Gflash::kHB] = std::max(0.0, einc*fTanh(einc,Gflash::mipcorr_pid[2]) +
62  exp(energyMeanHcal+energySigmaHcal*CLHEP::RandGaussQ::shoot())-2.0);
63  if(showerType == 2) {
64  energyScale[Gflash::kHB] = std::max(0.0,energyScale[Gflash::kHB]
65  - gap_corr*depthScale(position.getRho(),Gflash::Rmin[Gflash::kHB],28.));
66  }
67  }
68  else if(showerType == 6 || showerType == 7 ) {
69  energyScale[Gflash::kHE] = std::max(0.0, einc*fTanh(einc,Gflash::mipcorr_pid[3]) +
70  exp(energyMeanHcal+energySigmaHcal*CLHEP::RandGaussQ::shoot())-2.0);
71  if(showerType == 6) {
72  energyScale[Gflash::kHE] = std::max(0.0,energyScale[Gflash::kHE]
73  - gap_corr*depthScale(std::fabs(position.getZ()),Gflash::Zmin[Gflash::kHE],66.));
74  }
75  }
76  }
77 
78  // parameters for the longitudinal profiles
79  //@@@check longitudinal profiles of endcaps for possible variations
80  double *rhoHcal = new double [2*Gflash::NPar];
81  double *correlationVectorHcal = new double [Gflash::NPar*(Gflash::NPar+1)/2];
82 
83  //@@@until we have a separate parameterization for Endcap
84  bool isEndcap = false;
85  if(showerType>3) {
86  showerType -= 4;
87  isEndcap = true;
88  }
89  //no separate parameterization before crystal
90  if(showerType==0) showerType = 1;
91 
92  //Hcal parameters are always needed regardless of showerType
93  for(int i = 0 ; i < 2*Gflash::NPar ; i++ ) {
94  rhoHcal[i] = fTanh(einc,Gflash::rho[i + showerType*2*Gflash::NPar]);
95  }
96 
97  getFluctuationVector(rhoHcal,correlationVectorHcal);
98 
99  double normalZ[Gflash::NPar];
100  for (int i = 0; i < Gflash::NPar ; i++) normalZ[i] = CLHEP::RandGaussQ::shoot();
101 
102  for(int i = 0 ; i < Gflash::NPar ; i++) {
103  double correlationSum = 0.0;
104  for(int j = 0 ; j < i+1 ; j++) {
105  correlationSum += correlationVectorHcal[i*(i+1)/2+j]*normalZ[j];
106  }
107  longHcal[i] = fTanh(einc,Gflash::par[i+showerType*Gflash::NPar]) +
108  fTanh(einc,Gflash::par[i+(4+showerType)*Gflash::NPar])*correlationSum;
109  }
110 
111  delete [] rhoHcal;
112  delete [] correlationVectorHcal;
113 
114  // lateral parameters for Hcal
115 
116  for (int i = 0 ; i < Gflash::Nrpar ; i++) {
117  lateralPar[Gflash::kHB][i] = fLnE1(einc,Gflash::rpar[i+showerType*Gflash::Nrpar]);
118 
119  //tuning for pure hadronic response: +10%
120  if(showerType==3 && i == 0) lateralPar[Gflash::kHB][i] *= 1.1;
121  lateralPar[Gflash::kHE][i] = lateralPar[Gflash::kHB][i];
122  }
123 
124  //Ecal parameters are needed if and only if the shower starts inside the crystal
125 
126  if(showerType == 1) {
127  //A depth dependent correction for the core term of R in Hcal is the linear in
128  //the shower start point while for the spread term is nearly constant
129 
130  if(!isEndcap) lateralPar[Gflash::kHB][0] -= 2.3562e-01*(position.getRho()-131.0);
131  else lateralPar[Gflash::kHE][0] -= 2.3562e-01*(std::abs(position.getZ())-332.0);
132 
133  double *rhoEcal = new double [2*Gflash::NPar];
134  double *correlationVectorEcal = new double [Gflash::NPar*(Gflash::NPar+1)/2];
135  for(int i = 0 ; i < 2*Gflash::NPar ; i++ ) rhoEcal[i] = fTanh(einc,Gflash::rho[i]);
136 
137  getFluctuationVector(rhoEcal,correlationVectorEcal);
138 
139  for (int i = 0; i < Gflash::NPar ; i++) normalZ[i] = CLHEP::RandGaussQ::shoot();
140  for(int i = 0 ; i < Gflash::NPar ; i++) {
141  double correlationSum = 0.0;
142  for(int j = 0 ; j < i+1 ; j++) {
143  correlationSum += correlationVectorEcal[i*(i+1)/2+j]*normalZ[j];
144  }
145  longEcal[i] = fTanh(einc,Gflash::par[i]) +
146  fTanh(einc,Gflash::par[i+4*Gflash::NPar])*correlationSum;
147  }
148 
149  delete [] rhoEcal;
150  delete [] correlationVectorEcal;
151 
152  // lateral parameters for Ecal
153 
154  for (int i = 0 ; i < Gflash::Nrpar ; i++) {
155  lateralPar[Gflash::kESPM][i] = fLnE1(einc,Gflash::rpar[i]);
156  lateralPar[Gflash::kENCA][i] = lateralPar[Gflash::kESPM][i];
157  }
158  }
159 
160 }
Gflash::ZFrontCrystalEE
const double ZFrontCrystalEE
Definition: GflashNameSpace.h:36
GflashHadronShowerProfile::fTanh
double fTanh(double einc, const double *par)
Definition: GflashHadronShowerProfile.cc:472
i
int i
Definition: DBlmapReader.cc:9
GflashHadronShowerProfile::lateralPar
double lateralPar[Gflash::kNumberCalorimeter][Gflash::Nrpar]
Definition: GflashHadronShowerProfile.h:63
Gflash::correl_hadem
const double correl_hadem[4]
Definition: GflashNameSpace.h:113
GflashShowino::getEnergy
double getEnergy()
Definition: GflashShowino.h:24
GflashShowino::getShowerType
int getShowerType()
Definition: GflashShowino.h:23
abs
#define abs(x)
Definition: mlp_lapack.h:159
Gflash::emcorr_pid
const double emcorr_pid[10][4]
Definition: GflashNameSpace.h:136
Gflash::Nrpar
const int Nrpar
Definition: GflashNameSpace.h:274
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
Gflash::emscale
const double emscale[4][4]
Definition: GflashNameSpace.h:93
position
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
Gflash::hadcorr_pid
const double hadcorr_pid[10][4]
Definition: GflashNameSpace.h:149
Gflash::rpar
const double rpar[4 *Nrpar][2]
Definition: GflashNameSpace.h:275
Gflash::mipcorr_pid
const double mipcorr_pid[10][4]
Definition: GflashNameSpace.h:162
max
const T & max(const T &a, const T &b)
Definition: MaterialBudgetTrackerHistos.cc:4
GflashHadronShowerProfile::fLnE1
double fLnE1(double einc, const double *par)
Definition: GflashHadronShowerProfile.cc:478
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
Gflash::hadscale
const double hadscale[7][4]
Definition: GflashNameSpace.h:100
Gflash::LengthCrystalEE
const double LengthCrystalEE
Definition: GflashNameSpace.h:40
j
int j
Definition: DBlmapReader.cc:9
Gflash::LengthCrystalEB
const double LengthCrystalEB
Definition: GflashNameSpace.h:39
GflashShowino::getPositionAtShower
Gflash3Vector & getPositionAtShower()
Definition: GflashShowino.h:27
Gflash::NPar
const int NPar
Definition: GflashNameSpace.h:86
Gflash::RFrontCrystalEB
const double RFrontCrystalEB
Definition: GflashNameSpace.h:35
Gflash::Zmin
const double Zmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:27
Gflash::rho
const double rho[8 *NPar][4]
Definition: GflashNameSpace.h:227
Gflash3Vector
CLHEP::Hep3Vector Gflash3Vector
Definition: Gflash3Vector.h:6
GflashHadronShowerProfile::getFluctuationVector
void getFluctuationVector(double *lowTriangle, double *correlationVector)
Definition: GflashHadronShowerProfile.cc:357
GflashHadronShowerProfile::depthScale
double depthScale(double ssp, double ssp0, double length)
Definition: GflashHadronShowerProfile.cc:484
Gflash::Rmin
const double Rmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:29
Gflash::par
const double par[8 *NPar][4]
Definition: GflashNameSpace.h:177
GflashHadronShowerProfile::energyScale
double energyScale[Gflash::kNumberCalorimeter]
Definition: GflashHadronShowerProfile.h:59
Generated for CMSSW Reference Manual by doxygen 1.8.5