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 (const edm::ParameterSet &parSet)
 
void loadParameters ()
 
 ~GflashKaonMinusShowerProfile ()
 
- Public Member Functions inherited from GflashHadronShowerProfile
std::vector< GflashHit > & getGflashHitList ()
 
GflashShowino * getGflashShowino ()
 
 GflashHadronShowerProfile (const 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, const 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 ( const edm::ParameterSet &  parSet)
inline

Definition at line 12 of file GflashKaonMinusShowerProfile.h.

12  :
13  GflashHadronShowerProfile (parSet) {};
GflashHadronShowerProfile::GflashHadronShowerProfile
GflashHadronShowerProfile(const 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 GflashHadronShowerProfile::depthScale(), GflashHadronShowerProfile::energyScale, create_public_lumi_plots::exp, GflashHadronShowerProfile::fTanh(), GflashShowino::getEnergy(), GflashHadronShowerProfile::getFluctuationVector(), GflashShowino::getPositionAtShower(), GflashShowino::getShowerType(), i, j, Gflash::kENCA, Gflash::kESPM, Gflash::kHB, Gflash::kHE, Gflash::kminus_correl_hadem, Gflash::kminus_emscale, Gflash::kminus_hadscale, Gflash::kminus_par, Gflash::kminus_rho, Gflash::kminus_rpar, GflashHadronShowerProfile::lateralPar, Gflash::LengthCrystalEB, Gflash::LengthCrystalEE, cmsBatch::log, GflashHadronShowerProfile::longEcal, GflashHadronShowerProfile::longHcal, bookConverter::max, Gflash::NPar, Gflash::Nrpar, position, diffTwoXMLs::r1, diffTwoXMLs::r2, Gflash::RFrontCrystalEB, Gflash::Rmin, 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  double r1 = 0.0;
15  double r2 = 0.0;
16 
17  if(showerType == 0 || showerType == 1 || showerType == 4 || showerType == 5) {
18 
19  //@@@ energy dependent energyRho based on tuning with testbeam data
20  double energyRho = fTanh(einc,Gflash::kminus_correl_hadem);
21 
22  if(showerType == 0 || showerType == 1) {
23  do {
24  r1 = CLHEP::RandGaussQ::shoot();
25 
26  energyScale[Gflash::kESPM] = einc*(fTanh(einc,Gflash::kminus_emscale[0]) + fTanh(einc,Gflash::kminus_emscale[1])*r1);
27 
28  //LogNormal mean and sigma of Hcal energy
29  energyMeanHcal = (fTanh(einc,Gflash::kminus_hadscale[0]) +
30  fTanh(einc,Gflash::kminus_hadscale[1])*depthScale(position.getRho(),Gflash::RFrontCrystalEB,Gflash::LengthCrystalEB));
31  energySigmaHcal = (fTanh(einc,Gflash::kminus_hadscale[2]) +
32  fTanh(einc,Gflash::kminus_hadscale[3])*depthScale(position.getRho(),Gflash::RFrontCrystalEB,Gflash::LengthCrystalEB));
33 
34  r2 = CLHEP::RandGaussQ::shoot();
35  energyScale[Gflash::kHB] = exp(energyMeanHcal+energySigmaHcal*(energyRho*r1 + sqrt(1.0- energyRho*energyRho)*r2 ));
36  } while ( energyScale[Gflash::kESPM] < 0 || energyScale[Gflash::kHB] > einc*1.5 );
37  }
38  else {
39  do {
40  r1 = CLHEP::RandGaussQ::shoot();
41  energyScale[Gflash::kENCA] = einc*(fTanh(einc,Gflash::kminus_emscale[0]) + fTanh(einc,Gflash::kminus_emscale[1])*r1);
42 
43  //@@@extend depthScale for HE
44  energyMeanHcal = (fTanh(einc,Gflash::kminus_hadscale[0]) +
45  fTanh(einc,Gflash::kminus_hadscale[1])*depthScale(std::fabs(position.getZ()),Gflash::ZFrontCrystalEE,Gflash::LengthCrystalEE));
46  energySigmaHcal = (fTanh(einc,Gflash::kminus_hadscale[2]) +
47  fTanh(einc,Gflash::kminus_hadscale[3])*depthScale(std::fabs(position.getZ()),Gflash::ZFrontCrystalEE,Gflash::LengthCrystalEE));
48  r2 = CLHEP::RandGaussQ::shoot();
49  energyScale[Gflash::kHE] = exp(energyMeanHcal+energySigmaHcal*(energyRho*r1 + sqrt(1.0- energyRho*energyRho)*r2 ));
50  } while ( energyScale[Gflash::kENCA] < 0 || energyScale[Gflash::kHE] > einc*1.5 );
51  }
52  }
53  else if(showerType == 2 || showerType == 3 || showerType == 6 || showerType == 7) {
54  //Hcal response for mip-like pions (mip)
55 
56  energyMeanHcal = fTanh(einc,Gflash::kminus_hadscale[4]);
57  energySigmaHcal = fTanh(einc,Gflash::kminus_hadscale[5]);
58 
59  double gap_corr = einc*fTanh(einc,Gflash::kminus_hadscale[6]);
60 
61  if(showerType == 2 || showerType == 3) {
62  energyScale[Gflash::kESPM] = 0.0;
63 
64  do {
65  r1 = CLHEP::RandGaussQ::shoot();
66  energyScale[Gflash::kHB] = exp(energyMeanHcal+energySigmaHcal*r1);
67  } while ( energyScale[Gflash::kHB] > einc*1.5 );
68 
69  if(showerType == 2) {
70  energyScale[Gflash::kHE] = std::max(0.0,energyScale[Gflash::kHB]
71  - gap_corr*depthScale(position.getRho(),Gflash::Rmin[Gflash::kHB],28.));
72  }
73  }
74  else if(showerType == 6 || showerType == 7 ) {
75  energyScale[Gflash::kENCA] = 0.0;
76 
77  do {
78  r1 = CLHEP::RandGaussQ::shoot();
79  energyMeanHcal += std::log(1.0-fTanh(einc,Gflash::kminus_hadscale[7]));
80  energyScale[Gflash::kHE] = exp(energyMeanHcal+energySigmaHcal*r1);
81  } while ( energyScale[Gflash::kHE] > einc*1.5 );
82 
83 
84  if(showerType == 6) {
85  energyScale[Gflash::kHE] = std::max(0.0,energyScale[Gflash::kHE]
86  - gap_corr*depthScale(std::fabs(position.getZ()),Gflash::Zmin[Gflash::kHE],66.));
87  }
88  }
89  }
90 
91  // parameters for the longitudinal profiles
92  //@@@check longitudinal profiles of endcaps for possible variations
93 
94  double *rhoHcal = new double [2*Gflash::NPar];
95  double *correlationVectorHcal = new double [Gflash::NPar*(Gflash::NPar+1)/2];
96 
97  //@@@until we have a separate parameterization for Endcap
98 
99  if(showerType>3) {
100  showerType -= 4;
101  }
102  //no separate parameterization before crystal
103  if(showerType==0) showerType = 1;
104 
105  //Hcal parameters are always needed regardless of showerType
106  for(int i = 0 ; i < 2*Gflash::NPar ; i++ ) {
107  rhoHcal[i] = fTanh(einc,Gflash::kminus_rho[i + showerType*2*Gflash::NPar]);
108  }
109 
110  getFluctuationVector(rhoHcal,correlationVectorHcal);
111 
112  double normalZ[Gflash::NPar];
113  for (int i = 0; i < Gflash::NPar ; i++) normalZ[i] = CLHEP::RandGaussQ::shoot();
114 
115  for(int i = 0 ; i < Gflash::NPar ; i++) {
116  double correlationSum = 0.0;
117 
118  for(int j = 0 ; j < i+1 ; j++) {
119  correlationSum += correlationVectorHcal[i*(i+1)/2+j]*normalZ[j];
120  }
121  longHcal[i] = fTanh(einc,Gflash::kminus_par[i+showerType*Gflash::NPar]) +
122  fTanh(einc,Gflash::kminus_par[i+(4+showerType)*Gflash::NPar])*correlationSum;
123  }
124  delete [] rhoHcal;
125  delete [] correlationVectorHcal;
126 
127  // lateral parameters for Hcal
128 
129  for (int i = 0 ; i < Gflash::Nrpar ; i++) {
130  lateralPar[Gflash::kHB][i] = fTanh(einc,Gflash::kminus_rpar[i+showerType*Gflash::Nrpar]);
131  lateralPar[Gflash::kHE][i] = lateralPar[Gflash::kHB][i];
132  }
133 
134  //Ecal parameters are needed if and only if the shower starts inside the crystal
135 
136  if(showerType == 1) {
137 
138  double *rhoEcal = new double [2*Gflash::NPar];
139  double *correlationVectorEcal = new double [Gflash::NPar*(Gflash::NPar+1)/2];
140  for(int i = 0 ; i < 2*Gflash::NPar ; i++ ) rhoEcal[i] = fTanh(einc,Gflash::kminus_rho[i]);
141 
142  getFluctuationVector(rhoEcal,correlationVectorEcal);
143 
144  for(int i = 0 ; i < Gflash::NPar ; i++) normalZ[i] = CLHEP::RandGaussQ::shoot();
145  for(int i = 0 ; i < Gflash::NPar ; i++) {
146  double correlationSum = 0.0;
147 
148  for(int j = 0 ; j < i+1 ; j++) {
149  correlationSum += correlationVectorEcal[i*(i+1)/2+j]*normalZ[j];
150  }
151  longEcal[i] = fTanh(einc,Gflash::kminus_par[i]) +
152  0.5*fTanh(einc,Gflash::kminus_par[i+4*Gflash::NPar])*correlationSum;
153  }
154 
155  delete [] rhoEcal;
156  delete [] correlationVectorEcal;
157 
158  // lateral parameters for Ecal
159 
160  for (int i = 0 ; i < Gflash::Nrpar ; i++) {
161  lateralPar[Gflash::kESPM][i] = fTanh(einc,Gflash::kminus_rpar[i]);
162  lateralPar[Gflash::kENCA][i] = lateralPar[Gflash::kESPM][i];
163  }
164  }
165 
166 }
Gflash::ZFrontCrystalEE
const double ZFrontCrystalEE
Definition: GflashNameSpace.h:36
GflashHadronShowerProfile::fTanh
double fTanh(double einc, const double *par)
Definition: GflashHadronShowerProfile.cc:485
i
int i
Definition: DBlmapReader.cc:9
create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1087
Gflash::kminus_rho
const double kminus_rho[8 *NPar][5]
Definition: GflashNameSpace.h:603
GflashHadronShowerProfile::lateralPar
double lateralPar[Gflash::kNumberCalorimeter][Gflash::Nrpar]
Definition: GflashHadronShowerProfile.h:63
Gflash::kminus_rpar
const double kminus_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:651
GflashShowino::getEnergy
double getEnergy()
Definition: GflashShowino.h:24
GflashShowino::getShowerType
int getShowerType()
Definition: GflashShowino.h:23
Gflash::Nrpar
const int Nrpar
Definition: GflashNameSpace.h:88
Gflash::kminus_hadscale
const double kminus_hadscale[8][5]
Definition: GflashNameSpace.h:534
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
Gflash::kminus_par
const double kminus_par[8 *NPar][5]
Definition: GflashNameSpace.h:553
Gflash::LengthCrystalEE
const double LengthCrystalEE
Definition: GflashNameSpace.h:40
j
int j
Definition: DBlmapReader.cc:9
Gflash::kminus_emscale
const double kminus_emscale[2][5]
Definition: GflashNameSpace.h:529
Gflash::LengthCrystalEB
const double LengthCrystalEB
Definition: GflashNameSpace.h:39
GflashShowino::getPositionAtShower
Gflash3Vector & getPositionAtShower()
Definition: GflashShowino.h:27
Gflash::kminus_correl_hadem
const double kminus_correl_hadem[5]
Definition: GflashNameSpace.h:548
Gflash::NPar
const int NPar
Definition: GflashNameSpace.h:87
Gflash::RFrontCrystalEB
const double RFrontCrystalEB
Definition: GflashNameSpace.h:35
Gflash::Zmin
const double Zmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:27
Gflash3Vector
CLHEP::Hep3Vector Gflash3Vector
Definition: Gflash3Vector.h:6
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
GflashHadronShowerProfile::getFluctuationVector
void getFluctuationVector(double *lowTriangle, double *correlationVector)
Definition: GflashHadronShowerProfile.cc:370
GflashHadronShowerProfile::depthScale
double depthScale(double ssp, double ssp0, double length)
Definition: GflashHadronShowerProfile.cc:497
Gflash::Rmin
const double Rmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:29
cmsBatch.log
tuple log
Definition: cmsBatch.py:341
GflashHadronShowerProfile::energyScale
double energyScale[Gflash::kNumberCalorimeter]
Definition: GflashHadronShowerProfile.h:59
Generated for CMSSW Reference Manual by doxygen 1.8.5