CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
HDRShower Class Reference

#include <HDRShower.h>

Public Member Functions

bool computeShower ()
 
float getR ()
 
 HDRShower (const RandomEngineAndDistribution *engine, HDShowerParametrization *myParam, EcalHitMaker *myGrid, HcalHitMaker *myHcalHitMaker, int onECAL, double epart)
 
void setFuncParam ()
 
bool setHit (float espot, float theta)
 
void thetaFunction (int nthetaStep)
 
virtual ~HDRShower ()
 

Private Attributes

float decal
 
float depthECAL
 
float depthGAP
 
double depthStart
 
double e
 
float EcalShift
 
float EgridTable [7]
 
float eHDspot
 
std::vector< float > elastspot
 
float EsCut
 
double lambdaEM
 
double lambdaHD
 
float maxDepth
 
int nthetaStep
 
int onEcal
 
bool qstatus
 
const RandomEngineAndDistribution * random
 
float rpdf [100]
 
EcalHitMaker * theGrid
 
HcalHitMaker * theHcalHitMaker
 
HDShowerParametrization * theParam
 
float Theta1amp [7]
 
float Theta1ampSig [7]
 
float Theta1Lambda [7]
 
float Theta1LambdaSig [7]
 
float ThetaLam21 [7]
 
float ThetaLam21Sig [7]
 
std::vector< int > thetaSpots
 
float thetaStep
 
double x0EM
 
double x0HD
 

Detailed Description

Definition at line 22 of file HDRShower.h.

Constructor & Destructor Documentation

◆ HDRShower()

HDRShower::HDRShower ( const RandomEngineAndDistribution *  engine,
HDShowerParametrization *  myParam,
EcalHitMaker *  myGrid,
HcalHitMaker *  myHcalHitMaker,
int  onECAL,
double  epart 
)

Definition at line 23 of file HDRShower.cc.

29  : theParam(myParam), theGrid(myGrid), theHcalHitMaker(myHcalHitMaker), onEcal(onECAL), e(epart), random(engine) {
30  eHDspot = 0.2;
31  EsCut = 0.050;
32  EcalShift = 0.12;
33  nthetaStep = 10;
34  thetaStep = 0.5 * M_PI / nthetaStep;
35 
36  if (e < 0)
37  e = 0.;
38  setFuncParam();
39 }

References e, EcalShift, eHDspot, EsCut, M_PI, nthetaStep, setFuncParam(), and thetaStep.

◆ ~HDRShower()

virtual HDRShower::~HDRShower ( )
inlinevirtual

Definition at line 31 of file HDRShower.h.

31 { ; }

Member Function Documentation

◆ computeShower()

bool HDRShower::computeShower ( )

Definition at line 41 of file HDRShower.cc.

41  {
42  if (onEcal) {
43  depthECAL = theGrid->ecalTotalL0(); // ECAL depth segment
44  depthGAP = theGrid->ecalHcalGapTotalL0(); // GAP depth segment
45  } else
46  depthECAL = depthGAP = 0;
47 
48  float depthHCAL = theGrid->hcalTotalL0(); // HCAL depth segment
49 
50  // maxDepth = depthECAL + depthGAP + depthHCAL - 1.0;
51  maxDepth = depthECAL + depthHCAL - 0.5;
52  depthStart = log(1. / random->flatShoot()); // starting point lambda unts
53 
54  if (depthStart > maxDepth) {
55  depthStart = maxDepth * random->flatShoot();
56  if (depthStart < 0.)
57  depthStart = 0.;
58  }
59 
60  if (depthStart < EcalShift)
61  depthStart = EcalShift;
62 
63  decal = (depthECAL + depthStart) * 0.5;
64  qstatus = false;
65  if (decal < depthECAL) {
66  qstatus = theGrid->getPads(decal);
67  // if(!qstatus)
68  // cout<<" depth rejected by getQuads(decal="<<decal<<") status="<<qstatus
69  // <<" depthECAL="<<depthECAL<<endl;
70  }
71 
72  thetaFunction(nthetaStep);
73  int maxLoops = 10000;
74  float esum = e;
75  for (int itheta = 0; itheta < nthetaStep; itheta++) {
76  float theta, es;
77  for (int i = 0; i <= thetaSpots[itheta]; i++) {
78  if (i == thetaSpots[itheta])
79  es = elastspot[itheta];
80  else
81  es = eHDspot;
82  float loops = 0;
83  for (int j = 0; j < maxLoops; j++) {
84  theta = (itheta + random->flatShoot()) * thetaStep;
85  if (setHit(es, theta))
86  break;
87  loops++;
88  }
89  esum -= es; // to check only
90  }
91  }
92  return (true);
93 }

References decal, depthECAL, depthGAP, depthStart, e, EcalHitMaker::ecalHcalGapTotalL0(), EcalShift, EcalHitMaker::ecalTotalL0(), eHDspot, elastspot, RandomEngineAndDistribution::flatShoot(), EcalHitMaker::getPads(), EcalHitMaker::hcalTotalL0(), mps_fire::i, dqmiolumiharvest::j, dqm-mbProfile::log, maxDepth, AlignmentProducer_cff::maxLoops, nthetaStep, onEcal, qstatus, random, setHit(), theGrid, theta(), thetaFunction(), thetaSpots, and thetaStep.

Referenced by CalorimetryManager::HDShowerSimulation().

◆ getR()

float HDRShower::getR ( )

Definition at line 123 of file HDRShower.cc.

123  {
124  float p = random->flatShoot();
125  unsigned int i = 1;
126  while (rpdf[i] < p && i < R_range - 1) {
127  i++;
128  }
129  float r;
130  float dr = rpdf[i] - rpdf[i - 1];
131  if (dr != 0.0)
132  r = (float(i) + (p - rpdf[i - 1]) / dr) / lambdaHD;
133  else
134  r = float(i) / lambdaHD;
135  return (r);
136 }

References flavorHistoryFilter_cfi::dr, RandomEngineAndDistribution::flatShoot(), dqmMemoryStats::float, mps_fire::i, lambdaHD, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, R_range, random, and rpdf.

Referenced by setHit().

◆ setFuncParam()

void HDRShower::setFuncParam ( )

Definition at line 208 of file HDRShower.cc.

208  {
209  lambdaHD = theParam->hcalProperties()->interactionLength();
210  x0HD = theParam->hcalProperties()->radLenIncm();
211  if (onEcal)
212  lambdaEM = theParam->ecalProperties()->interactionLength();
213  else
214  lambdaEM = lambdaHD;
215 
216  if (debug)
217  LogDebug("FastCalorimetry") << "setFuncParam-> lambdaEM=" << lambdaEM << " lambdaHD=" << lambdaHD << endl;
218 
219  float _EgridTable[NEnergyScan] = {10, 20, 30, 50, 100, 300, 500};
220  float _Theta1amp[NEnergyScan] = {1.57, 2.05, 2.27, 2.52, 2.66, 2.76, 2.76};
221  float _Theta1ampSig[NEnergyScan] = {2.40, 1.50, 1.25, 1.0, 0.8, 0.52, 0.52};
222 
223  float _Theta1Lambda[NEnergyScan] = {0.086, 0.092, 0.88, 0.80, 0.0713, 0.0536, 0.0536};
224  float _Theta1LambdaSig[NEnergyScan] = {0.038, 0.037, 0.027, 0.03, 0.023, 0.018, 0.018};
225 
226  float _ThetaLam21[NEnergyScan] = {2.8, 2.44, 2.6, 2.77, 3.16, 3.56, 3.56};
227  float _ThetaLam21Sig[NEnergyScan] = {1.8, 0.97, 0.87, 0.77, 0.7, 0.49, 0.49};
228 
229  for (int i = 0; i < NEnergyScan; i++) {
230  EgridTable[i] = _EgridTable[i];
231  Theta1amp[i] = _Theta1amp[i];
232  Theta1ampSig[i] = _Theta1ampSig[i];
233  Theta1Lambda[i] = _Theta1Lambda[i];
234  Theta1LambdaSig[i] = _Theta1LambdaSig[i];
235  ThetaLam21[i] = _ThetaLam21[i];
236  ThetaLam21Sig[i] = _ThetaLam21Sig[i];
237  }
238 
239 #define lambdafit 15.05
240  float R_alfa = -0.0993 + 0.1114 * log(e);
241  float R_p = 0.589191 + 0.0463392 * log(e);
242  float R_beta_lam = (0.54134 - 0.00011148 * e) / 4.0 * lambdafit; //was fitted in 4cmbin
243  float LamOverX0 = lambdaHD / x0HD; // 10.52
244  // int R_range = 100; // 7 lambda
245  // rpdf.erase(rpdf.begin(),rpdf.end());
246 
247  rpdf[0] = 0.;
248  for (int i = 1; i < R_range; i++) {
249  float x = (float(i)) / lambdaHD;
250  float r = pow(x, R_alfa) * (R_p * exp(-R_beta_lam * x) + (1 - R_p) * exp(-LamOverX0 * R_beta_lam * x));
251  rpdf[i] = r;
252  // rpdf.push_back(r);
253  }
254 
255  for (int i = 1; i < R_range; i++)
256  rpdf[i] += rpdf[i - 1];
257  for (int i = 0; i < R_range; i++)
258  rpdf[i] /= rpdf[R_range - 1];
259 }

References debug, e, HDShowerParametrization::ecalProperties(), EgridTable, JetChargeProducer_cfi::exp, dqmMemoryStats::float, HDShowerParametrization::hcalProperties(), mps_fire::i, ECALProperties::interactionLength(), HCALProperties::interactionLength(), lambdaEM, lambdafit, lambdaHD, dqm-mbProfile::log, LogDebug, NEnergyScan, onEcal, funct::pow(), alignCSCRings::r, R_range, HCALProperties::radLenIncm(), rpdf, theParam, Theta1amp, Theta1ampSig, Theta1Lambda, Theta1LambdaSig, ThetaLam21, ThetaLam21Sig, x, and x0HD.

Referenced by HDRShower().

◆ setHit()

bool HDRShower::setHit ( float  espot,
float  theta 
)

Definition at line 95 of file HDRShower.cc.

95  {
96  float phi = 2. * M_PI * random->flatShoot(); // temporary: 1st approximation
97  float rshower = getR(); // temporary: 1st approximation
98 
99  float d = depthStart + rshower * cos(theta);
100  if (d + depthGAP > maxDepth)
101  return (false);
102 
103  // Commented (F.B) to remove a warning. Not used anywhere ?
104  // bool inHcal = !onEcal || d>depthECAL || !qstatus;
105  bool result = false;
106  if (!onEcal || d > depthECAL || !qstatus) { // in HCAL (HF or HB, HE)
107  d += depthGAP;
108  bool setHDdepth = theHcalHitMaker->setDepth(d);
109  if (setHDdepth) {
110  theHcalHitMaker->setSpotEnergy(espot);
111  result = theHcalHitMaker->addHit(rshower * sin(theta), phi, 0);
112  } else
113  LogWarning("FastCalorimetry") << " setHit in HCAL failed d=" << d << " maxDepth=" << maxDepth << " onEcal'"
114  << onEcal << endl;
115  } else {
116  // bool status = theGrid->getQuads(d);
117  theGrid->setSpotEnergy(espot);
118  result = theGrid->addHit(rshower * sin(theta), phi, 0);
119  }
120  return (result);
121 }

References HcalHitMaker::addHit(), EcalHitMaker::addHit(), funct::cos(), ztail::d, depthECAL, depthGAP, depthStart, RandomEngineAndDistribution::flatShoot(), getR(), M_PI, maxDepth, onEcal, phi, qstatus, random, mps_fire::result, HcalHitMaker::setDepth(), HcalHitMaker::setSpotEnergy(), EcalHitMaker::setSpotEnergy(), funct::sin(), theGrid, theHcalHitMaker, and theta().

Referenced by computeShower().

◆ thetaFunction()

void HDRShower::thetaFunction ( int  nthetaStep)

Definition at line 138 of file HDRShower.cc.

138  {
139  unsigned int i = 0;
140  while (EgridTable[i] < e && i < NEnergyScan - 1) {
141  i++;
142  }
143 
144  float amean, asig, lambda1, lambda1sig, lam21, lam21sig;
145  amean = Theta1amp[i];
146  asig = Theta1ampSig[i];
147  lambda1 = Theta1Lambda[i];
148  lambda1sig = Theta1LambdaSig[i];
149  lam21 = ThetaLam21[i];
150  lam21sig = ThetaLam21Sig[i];
151  if (i == 0)
152  i = 1; //extrapolation to the left
153  float c = (e - EgridTable[i - 1]) / (EgridTable[i] - EgridTable[i - 1]);
154 
155  amean += (Theta1amp[i] - Theta1amp[i - 1]) * c;
156  asig += (Theta1ampSig[i] - Theta1ampSig[i - 1]) * c;
157  lambda1 += (Theta1Lambda[i] - Theta1Lambda[i - 1]) * c;
158  lambda1sig += (Theta1LambdaSig[i] - Theta1LambdaSig[i - 1]) * c;
159  lam21 += (ThetaLam21[i] - ThetaLam21[i - 1]) * c;
160  lam21sig += (ThetaLam21Sig[i] - ThetaLam21Sig[i - 1]) * c;
161 
162  float a = exp(amean + asig * random->gaussShoot());
163  float L1 = lambda1 + lambda1sig * random->gaussShoot();
164  if (L1 < 0.02)
165  L1 = 0.02;
166  float L2 = L1 * (lam21 + lam21sig * random->gaussShoot());
167 
168  vector<double> pdf;
169  pdf.erase(pdf.begin(), pdf.end());
170  thetaSpots.erase(thetaSpots.begin(), thetaSpots.end());
171  elastspot.erase(elastspot.begin(), elastspot.end());
172  double sum = 0;
173  for (int it = 0; it < nthetaStep; it++) {
174  float theta = it * thetaStep;
175  float p = a * exp(L1 * theta) + exp(L2 * theta);
176  sum += p;
177  pdf.push_back(p);
178  }
179  float ntot = e / eHDspot;
180  float esum = 0;
181  for (int it = 0; it < nthetaStep; it++) {
182  float fn = ntot * pdf[it] / sum;
183  thetaSpots.push_back(int(fn));
184  elastspot.push_back((fn - int(fn)) * eHDspot);
185  }
186 
187  for (int it = 0; it < nthetaStep; it++)
188  if (elastspot[it] < EsCut) {
189  esum += elastspot[it];
190  elastspot[it] = 0;
191  }
192 
193  float en = esum / EsCut;
194  int n = int(en);
195  en = esum - n * EsCut;
196 
197  for (int ie = 0; ie <= n; ie++) {
198  int k = int(nthetaStep * random->flatShoot());
199  if (k < 0 || k > nthetaStep - 1)
200  k = k % nthetaStep;
201  if (ie == n)
202  elastspot[k] += en;
203  else
204  elastspot[k] += EsCut;
205  }
206 }

References a, c, e, EgridTable, eHDspot, elastspot, EsCut, JetChargeProducer_cfi::exp, RandomEngineAndDistribution::flatShoot(), personalPlayback::fn, RandomEngineAndDistribution::gaussShoot(), mps_fire::i, createfilelist::int, dqmdumpme::k, PixelRegions::L1, PixelRegions::L2, dqmiodumpmetadata::n, NEnergyScan, nthetaStep, gpuPixelDoublets::ntot, AlCaHLTBitMon_ParallelJobs::p, random, theta(), Theta1amp, Theta1ampSig, Theta1Lambda, Theta1LambdaSig, ThetaLam21, ThetaLam21Sig, thetaSpots, and thetaStep.

Referenced by computeShower().

Member Data Documentation

◆ decal

float HDRShower::decal
private

Definition at line 64 of file HDRShower.h.

Referenced by computeShower().

◆ depthECAL

float HDRShower::depthECAL
private

Definition at line 59 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ depthGAP

float HDRShower::depthGAP
private

Definition at line 59 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ depthStart

double HDRShower::depthStart
private

Definition at line 52 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ e

double HDRShower::e
private

Definition at line 45 of file HDRShower.h.

Referenced by computeShower(), HDRShower(), setFuncParam(), and thetaFunction().

◆ EcalShift

float HDRShower::EcalShift
private

Definition at line 55 of file HDRShower.h.

Referenced by computeShower(), and HDRShower().

◆ EgridTable

float HDRShower::EgridTable[7]
private

Definition at line 66 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ eHDspot

float HDRShower::eHDspot
private

Definition at line 53 of file HDRShower.h.

Referenced by computeShower(), HDRShower(), and thetaFunction().

◆ elastspot

std::vector<float> HDRShower::elastspot
private

Definition at line 61 of file HDRShower.h.

Referenced by computeShower(), and thetaFunction().

◆ EsCut

float HDRShower::EsCut
private

Definition at line 54 of file HDRShower.h.

Referenced by HDRShower(), and thetaFunction().

◆ lambdaEM

double HDRShower::lambdaEM
private

Definition at line 51 of file HDRShower.h.

Referenced by setFuncParam().

◆ lambdaHD

double HDRShower::lambdaHD
private

Definition at line 51 of file HDRShower.h.

Referenced by getR(), and setFuncParam().

◆ maxDepth

float HDRShower::maxDepth
private

Definition at line 59 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ nthetaStep

int HDRShower::nthetaStep
private

Definition at line 56 of file HDRShower.h.

Referenced by computeShower(), HDRShower(), and thetaFunction().

◆ onEcal

int HDRShower::onEcal
private

Definition at line 44 of file HDRShower.h.

Referenced by computeShower(), setFuncParam(), and setHit().

◆ qstatus

bool HDRShower::qstatus
private

Definition at line 63 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ random

const RandomEngineAndDistribution* HDRShower::random
private

Definition at line 75 of file HDRShower.h.

Referenced by computeShower(), getR(), setHit(), and thetaFunction().

◆ rpdf

float HDRShower::rpdf[100]
private

Definition at line 62 of file HDRShower.h.

Referenced by getR(), and setFuncParam().

◆ theGrid

EcalHitMaker* HDRShower::theGrid
private

Definition at line 42 of file HDRShower.h.

Referenced by computeShower(), and setHit().

◆ theHcalHitMaker

HcalHitMaker* HDRShower::theHcalHitMaker
private

Definition at line 43 of file HDRShower.h.

Referenced by setHit().

◆ theParam

HDShowerParametrization* HDRShower::theParam
private

Definition at line 41 of file HDRShower.h.

Referenced by setFuncParam().

◆ Theta1amp

float HDRShower::Theta1amp[7]
private

Definition at line 67 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ Theta1ampSig

float HDRShower::Theta1ampSig[7]
private

Definition at line 68 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ Theta1Lambda

float HDRShower::Theta1Lambda[7]
private

Definition at line 69 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ Theta1LambdaSig

float HDRShower::Theta1LambdaSig[7]
private

Definition at line 70 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ ThetaLam21

float HDRShower::ThetaLam21[7]
private

Definition at line 71 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ ThetaLam21Sig

float HDRShower::ThetaLam21Sig[7]
private

Definition at line 72 of file HDRShower.h.

Referenced by setFuncParam(), and thetaFunction().

◆ thetaSpots

std::vector<int> HDRShower::thetaSpots
private

Definition at line 60 of file HDRShower.h.

Referenced by computeShower(), and thetaFunction().

◆ thetaStep

float HDRShower::thetaStep
private

Definition at line 58 of file HDRShower.h.

Referenced by computeShower(), HDRShower(), and thetaFunction().

◆ x0EM

double HDRShower::x0EM
private

Definition at line 51 of file HDRShower.h.

◆ x0HD

double HDRShower::x0HD
private

Definition at line 51 of file HDRShower.h.

Referenced by setFuncParam().

PixelRegions::L1
Definition: PixelRegionContainers.h:32
HDRShower::maxDepth
float maxDepth
Definition: HDRShower.h:59
HDRShower::qstatus
bool qstatus
Definition: HDRShower.h:63
mps_fire.i
i
Definition: mps_fire.py:428
PixelRegions::L2
Definition: PixelRegionContainers.h:32
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
HDRShower::Theta1Lambda
float Theta1Lambda[7]
Definition: HDRShower.h:69
HDRShower::EsCut
float EsCut
Definition: HDRShower.h:54
HDRShower::setHit
bool setHit(float espot, float theta)
Definition: HDRShower.cc:95
HDRShower::theParam
HDShowerParametrization * theParam
Definition: HDRShower.h:41
HDRShower::getR
float getR()
Definition: HDRShower.cc:123
HDRShower::depthGAP
float depthGAP
Definition: HDRShower.h:59
AlignmentProducer_cff.maxLoops
maxLoops
Definition: AlignmentProducer_cff.py:15
HDRShower::thetaSpots
std::vector< int > thetaSpots
Definition: HDRShower.h:60
HCALProperties::radLenIncm
double radLenIncm() const override
Radiation length in cm.
Definition: HCALProperties.h:37
HDRShower::Theta1amp
float Theta1amp[7]
Definition: HDRShower.h:67
HDRShower::setFuncParam
void setFuncParam()
Definition: HDRShower.cc:208
gpuPixelDoublets::ntot
__shared__ uint32_t ntot
Definition: gpuPixelDoubletsAlgos.h:67
HDRShower::ThetaLam21Sig
float ThetaLam21Sig[7]
Definition: HDRShower.h:72
HDRShower::EgridTable
float EgridTable[7]
Definition: HDRShower.h:66
HDRShower::rpdf
float rpdf[100]
Definition: HDRShower.h:62
DDAxes::x
HDRShower::theGrid
EcalHitMaker * theGrid
Definition: HDRShower.h:42
EcalHitMaker::getPads
bool getPads(double depth, bool inCm=false)
Definition: EcalHitMaker.cc:820
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
EcalHitMaker::addHit
bool addHit(double r, double phi, unsigned layer=0) override
Definition: EcalHitMaker.cc:184
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
HDRShower::x0HD
double x0HD
Definition: HDRShower.h:51
debug
#define debug
Definition: HDRShower.cc:19
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
HDRShower::thetaFunction
void thetaFunction(int nthetaStep)
Definition: HDRShower.cc:138
HDRShower::depthECAL
float depthECAL
Definition: HDRShower.h:59
HCALProperties::interactionLength
double interactionLength() const override
Interaction length in cm.
Definition: HCALProperties.h:56
EcalHitMaker::ecalHcalGapTotalL0
double ecalHcalGapTotalL0() const
ECAL-HCAL transition.
Definition: EcalHitMaker.h:91
HDRShower::random
const RandomEngineAndDistribution * random
Definition: HDRShower.h:75
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:150
HDRShower::eHDspot
float eHDspot
Definition: HDRShower.h:53
HDRShower::decal
float decal
Definition: HDRShower.h:64
dqmdumpme.k
k
Definition: dqmdumpme.py:60
R_range
#define R_range
Definition: HDRShower.h:16
HDRShower::elastspot
std::vector< float > elastspot
Definition: HDRShower.h:61
HDRShower::theHcalHitMaker
HcalHitMaker * theHcalHitMaker
Definition: HDRShower.h:43
HDRShower::nthetaStep
int nthetaStep
Definition: HDRShower.h:56
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
a
double a
Definition: hdecay.h:119
HDRShower::e
double e
Definition: HDRShower.h:45
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
EcalHitMaker::setSpotEnergy
void setSpotEnergy(double e) override
Definition: EcalHitMaker.h:112
createfilelist.int
int
Definition: createfilelist.py:10
lambdafit
#define lambdafit
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
HDRShower::lambdaEM
double lambdaEM
Definition: HDRShower.h:51
HDRShower::Theta1ampSig
float Theta1ampSig[7]
Definition: HDRShower.h:68
HDShowerParametrization::hcalProperties
const HCALProperties * hcalProperties() const
Definition: HDShowerParametrization.h:24
HDRShower::onEcal
int onEcal
Definition: HDRShower.h:44
RandomEngineAndDistribution::flatShoot
double flatShoot(double xmin=0.0, double xmax=1.0) const
Definition: RandomEngineAndDistribution.h:27
HcalHitMaker::setDepth
bool setDepth(double, bool inCm=false)
set the depth in X0 or Lambda0 units depending on showerType
Definition: HcalHitMaker.cc:114
alignCSCRings.r
r
Definition: alignCSCRings.py:93
NEnergyScan
#define NEnergyScan
Definition: HDRShower.h:14
DDAxes::phi
HDRShower::thetaStep
float thetaStep
Definition: HDRShower.h:58
HDShowerParametrization::ecalProperties
const ECALProperties * ecalProperties() const
Definition: HDShowerParametrization.h:22
HDRShower::EcalShift
float EcalShift
Definition: HDRShower.h:55
HDRShower::Theta1LambdaSig
float Theta1LambdaSig[7]
Definition: HDRShower.h:70
HDRShower::lambdaHD
double lambdaHD
Definition: HDRShower.h:51
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
personalPlayback.fn
fn
Definition: personalPlayback.py:515
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
ECALProperties::interactionLength
double interactionLength() const override
Interaction length in cm: 18.5 for Standard ECAL.
Definition: ECALProperties.h:49
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
ztail.d
d
Definition: ztail.py:151
mps_fire.result
result
Definition: mps_fire.py:311
HcalHitMaker::addHit
bool addHit(double r, double phi, unsigned layer=0) override
add the hit in the HCAL in local coordinates
Definition: HcalHitMaker.cc:30
HDRShower::ThetaLam21
float ThetaLam21[7]
Definition: HDRShower.h:71
HDRShower::depthStart
double depthStart
Definition: HDRShower.h:52
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:56
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
EcalHitMaker::ecalTotalL0
double ecalTotalL0() const
in the ECAL
Definition: EcalHitMaker.h:85
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
EcalHitMaker::hcalTotalL0
double hcalTotalL0() const
in the HCAL
Definition: EcalHitMaker.h:88
HcalHitMaker::setSpotEnergy
void setSpotEnergy(double e) override
Set the spot energy.
Definition: HcalHitMaker.h:26
RandomEngineAndDistribution::gaussShoot
double gaussShoot(double mean=0.0, double sigma=1.0) const
Definition: RandomEngineAndDistribution.h:29
Generated for CMSSW Reference Manual by doxygen 1.8.16