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

#include <HDShower.h>

Public Types

typedef std::pair< XYZPoint, double > Spot
 
typedef std::pair< unsigned int, double > Step
 
typedef Steps::const_iterator step_iterator
 
typedef std::vector< Step > Steps
 
typedef math::XYZVector XYZPoint
 

Public Member Functions

bool compute ()
 Compute the shower longitudinal and lateral development. More...
 
int getmip ()
 
 HDShower (const RandomEngineAndDistribution *engine, HDShowerParametrization *myParam, EcalHitMaker *myGrid, HcalHitMaker *myHcalHitMaker, int onECAL, double epart, double pmip)
 
virtual ~HDShower ()
 

Private Member Functions

double gam (double x, double a) const
 
int indexFinder (double x, const std::vector< double > &Fhist)
 
void makeSteps (int nsteps)
 
double transProb (double factor, double R, double r)
 

Private Attributes

double aloge
 
double alpEM
 
double alpHD
 
double balanceEH
 
double betEM
 
double betHD
 
double criticalEnergy
 
double depthECAL
 
double depthGAP
 
double depthGAPx0
 
double depthHCAL
 
double depthStart
 
double depthStep
 
double depthToHCAL
 
std::vector< int > detector
 
double e
 
double eSpotSize
 
std::vector< double > eStep
 
double hcalDepthFactor
 
int infinity
 
double lambdaEM
 
double lambdaHD
 
std::vector< double > lamcurr
 
std::vector< double > lamdepth
 
std::vector< double > lamstep
 
std::vector< double > lamtotal
 
int lossesOpt
 
double maxTRfactor
 
int mip
 
int nDepthSteps
 
std::vector< int > nspots
 
int nTRsteps
 
int onEcal
 
double part
 
const RandomEngineAndDistribution * random
 
std::vector< double > rlamStep
 
double tgamEM
 
double tgamHD
 
const ECALProperties * theECALproperties
 
EcalHitMaker * theGrid
 
HcalHitMaker * theHcalHitMaker
 
const HCALProperties * theHCALproperties
 
HDShowerParametrization * theParam
 
double theR1
 
double theR2
 
double theR3
 
double transFactor
 
double transParam
 
std::vector< double > x0curr
 
std::vector< double > x0depth
 
double x0EM
 
double x0HD
 

Detailed Description

Definition at line 22 of file HDShower.h.

Member Typedef Documentation

◆ Spot

typedef std::pair<XYZPoint, double> HDShower::Spot

Definition at line 26 of file HDShower.h.

◆ Step

typedef std::pair<unsigned int, double> HDShower::Step

Definition at line 27 of file HDShower.h.

◆ step_iterator

typedef Steps::const_iterator HDShower::step_iterator

Definition at line 29 of file HDShower.h.

◆ Steps

typedef std::vector<Step> HDShower::Steps

Definition at line 28 of file HDShower.h.

◆ XYZPoint

typedef math::XYZVector HDShower::XYZPoint

Definition at line 24 of file HDShower.h.

Constructor & Destructor Documentation

◆ HDShower()

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

Definition at line 21 of file HDShower.cc.

28  : theParam(myParam),
29  theGrid(myGrid),
30  theHcalHitMaker(myHcalHitMaker),
31  onEcal(onECAL),
32  e(epart),
33  // pmip(pmip),
34  random(engine) {
35  // To get an access to constants read in FASTCalorimeter
36  // FASTCalorimeter * myCalorimeter= FASTCalorimeter::instance();
37 
38  // Values taken from FamosGeneric/FamosCalorimeter/src/FASTCalorimeter.cc
39  lossesOpt = myParam->hsParameters()->getHDlossesOpt();
40  nDepthSteps = myParam->hsParameters()->getHDnDepthSteps();
41  nTRsteps = myParam->hsParameters()->getHDnTRsteps();
42  transParam = myParam->hsParameters()->getHDtransParam();
43  eSpotSize = myParam->hsParameters()->getHDeSpotSize();
44  depthStep = myParam->hsParameters()->getHDdepthStep();
45  criticalEnergy = myParam->hsParameters()->getHDcriticalEnergy();
46  maxTRfactor = myParam->hsParameters()->getHDmaxTRfactor();
47  balanceEH = myParam->hsParameters()->getHDbalanceEH();
48  hcalDepthFactor = myParam->hsParameters()->getHDhcalDepthFactor();
49 
50  // Special tr.size fluctuations
51  transParam *= (1. + random->flatShoot());
52 
53  // Special long. fluctuations
54  hcalDepthFactor += 0.05 * (2. * random->flatShoot() - 1.);
55 
56  transFactor = 1.; // normally 1, in HF - might be smaller
57  // to take into account
58  // a narrowness of the HF shower (Cherenkov light)
59 
60  // simple protection ...
61  if (e < 0)
62  e = 0.;
63 
64  // Get the Famos Histos pointer
65  // myHistos = FamosHistos::instance();
66  // std::cout << " Hello FamosShower " << std::endl;
67 
68  theECALproperties = theParam->ecalProperties();
69  theHCALproperties = theParam->hcalProperties();
70 
71  double emax = theParam->emax();
72  double emid = theParam->emid();
73  double emin = theParam->emin();
74  double effective = e;
75 
76  if (e < emid) {
77  theParam->setCase(1);
78  // avoid "underflow" below Emin (for parameters calculation only)
79  if (e < emin)
80  effective = emin;
81  } else
82  theParam->setCase(2);
83 
84  // Avoid "overflow" beyond Emax (for parameters)
85  if (effective > 0.5 * emax) {
86  eSpotSize *= 2.5;
87  if (effective > emax) {
88  effective = emax;
89  eSpotSize *= 4.;
90  depthStep *= 2.;
91  if (effective > 1000.)
92  eSpotSize *= 2.;
93  }
94  }
95 
96  if (debug == 2)
97  LogInfo("FastCalorimetry") << " HDShower : " << std::endl
98  << " Energy " << e << std::endl
99  << " lossesOpt " << lossesOpt << std::endl
100  << " nDepthSteps " << nDepthSteps << std::endl
101  << " nTRsteps " << nTRsteps << std::endl
102  << " transParam " << transParam << std::endl
103  << " eSpotSize " << eSpotSize << std::endl
104  << " criticalEnergy " << criticalEnergy << std::endl
105  << " maxTRfactor " << maxTRfactor << std::endl
106  << " balanceEH " << balanceEH << std::endl
107  << "hcalDepthFactor " << hcalDepthFactor << std::endl;
108 
109  double alpEM1 = theParam->alpe1();
110  double alpEM2 = theParam->alpe2();
111 
112  double betEM1 = theParam->bete1();
113  double betEM2 = theParam->bete2();
114 
115  double alpHD1 = theParam->alph1();
116  double alpHD2 = theParam->alph2();
117 
118  double betHD1 = theParam->beth1();
119  double betHD2 = theParam->beth2();
120 
121  double part1 = theParam->part1();
122  double part2 = theParam->part2();
123 
124  aloge = std::log(effective);
125 
126  double edpar = (theParam->e1() + aloge * theParam->e2()) * effective;
127  double aedep = std::log(edpar);
128 
129  if (debug == 2)
130  LogInfo("FastCalorimetry") << " HDShower : " << std::endl
131  << " edpar " << edpar << " aedep " << aedep << std::endl
132  << " alpEM1 " << alpEM1 << std::endl
133  << " alpEM2 " << alpEM2 << std::endl
134  << " betEM1 " << betEM1 << std::endl
135  << " betEM2 " << betEM2 << std::endl
136  << " alpHD1 " << alpHD1 << std::endl
137  << " alpHD2 " << alpHD2 << std::endl
138  << " betHD1 " << betHD1 << std::endl
139  << " betHD2 " << betHD2 << std::endl
140  << " part1 " << part1 << std::endl
141  << " part2 " << part2 << std::endl;
142 
143  // private members to set
144  theR1 = theParam->r1();
145  theR2 = theParam->r2();
146  theR3 = theParam->r3();
147 
148  alpEM = alpEM1 + alpEM2 * aedep;
149  tgamEM = tgamma(alpEM);
150  betEM = betEM1 - betEM2 * aedep;
151  alpHD = alpHD1 + alpHD2 * aedep;
152  tgamHD = tgamma(alpHD);
153  betHD = betHD1 - betHD2 * aedep;
154  part = part1 - part2 * aedep;
155  if (part > 1.)
156  part = 1.; // protection - just in case of
157 
158  if (debug == 2)
159  LogInfo("FastCalorimetry") << " HDShower : " << std::endl
160  << " alpEM " << alpEM << std::endl
161  << " tgamEM " << tgamEM << std::endl
162  << " betEM " << betEM << std::endl
163  << " alpHD " << alpHD << std::endl
164  << " tgamHD " << tgamHD << std::endl
165  << " betHD " << betHD << std::endl
166  << " part " << part << std::endl;
167 
168  if (onECAL) {
169  lambdaEM = theParam->ecalProperties()->interactionLength();
170  x0EM = theParam->ecalProperties()->radLenIncm();
171  } else {
172  lambdaEM = 0.;
173  x0EM = 0.;
174  }
175  lambdaHD = theParam->hcalProperties()->interactionLength();
176  x0HD = theParam->hcalProperties()->radLenIncm();
177 
178  if (debug == 2)
179  LogInfo("FastCalorimetry") << " HDShower e " << e << std::endl
180  << " x0EM = " << x0EM << std::endl
181  << " x0HD = " << x0HD << std::endl
182  << " lamEM = " << lambdaEM << std::endl
183  << " lamHD = " << lambdaHD << std::endl;
184 
185  // Starting point of the shower
186  // try first with ECAL lambda
187 
188  double sum1 = 0.; // lambda path from the ECAL/HF entrance;
189  double sum2 = 0.; // lambda path from the interaction point;
190  double sum3 = 0.; // x0 path from the interaction point;
191  int nsteps = 0; // full number of longitudinal steps (counter);
192 
193  int nmoresteps; // how many longitudinal steps in addition to
194  // one (if interaction happens there) in ECAL
195 
196  mip = 1; // just to initiate particle as MIP in ECAL
197 
198  if (e < criticalEnergy)
199  nmoresteps = 1;
200  else
201  nmoresteps = nDepthSteps;
202 
203  depthECAL = 0.;
204  depthGAP = 0.;
205  depthGAPx0 = 0.;
206  if (onECAL) {
207  depthECAL = theGrid->ecalTotalL0(); // ECAL depth segment
208  //depthECAL = theGrid->ecalTotalL0() + theGrid->ps1TotalL0() + theGrid->ps2TotalL0() + theGrid->ps2eeTotalL0(); //TEST: include preshower
209  depthGAP = theGrid->ecalHcalGapTotalL0(); // GAP depth segment
210  depthGAPx0 = theGrid->ecalHcalGapTotalX0(); // GAP depth x0
211  }
212 
213  depthHCAL = theGrid->hcalTotalL0(); // HCAL depth segment
214  depthToHCAL = depthECAL + depthGAP;
215 
216  //---------------------------------------------------------------------------
217  // Depth simulation & various protections, among them
218  // if too deep - get flat random in the allowed region
219  // if no HCAL material behind - force to deposit in ECAL
220  double maxDepth = depthToHCAL + depthHCAL - 1.1 * depthStep;
221  double depthStart = std::log(1. / random->flatShoot()); // starting point lambda unts
222  //std::cout<<"generated depth "<<depthStart<<std::endl;
223  if (pmip == 1) {
224  depthStart = depthToHCAL;
225  } else {
226  depthStart = depthStart * 0.9 * depthECAL / std::log(1. / pmip);
227  }
228  // std::cout<<"modified depth "<<depthStart<<std::endl;
229 
230  if (e < emin) {
231  if (debug)
232  LogInfo("FastCalorimetry") << " FamosHDShower : e <emin -> depthStart = 0" << std::endl;
233  depthStart = 0.;
234  }
235 
236  if (depthStart > maxDepth) {
237  if (debug)
238  LogInfo("FastCalorimetry") << " FamosHDShower : depthStart too big ... = " << depthStart << std::endl;
239  depthStart = maxDepth * random->flatShoot();
240  if (depthStart < 0.)
241  depthStart = 0.;
242  if (debug)
243  LogInfo("FastCalorimetry") << " FamosHDShower : depthStart re-calculated = " << depthStart << std::endl;
244  }
245 
246  if (onECAL && e < emid) {
247  if (depthECAL > depthStep && (depthECAL - depthStart) / depthECAL > 0.2) {
248  depthStart = 0.5 * depthECAL * random->flatShoot();
249  if (debug)
250  LogInfo("FastCalorimetry") << " FamosHDShower : small energy, "
251  << " depthStart reduced to = " << depthStart << std::endl;
252  }
253  }
254 
255  if (depthHCAL < depthStep) {
256  if (debug)
257  LogInfo("FastCalorimetry") << " FamosHDShower : depthHCAL too small ... = " << depthHCAL
258  << " depthStart -> forced to 0 !!!" << std::endl;
259  depthStart = 0.;
260  nmoresteps = 0;
261 
262  if (depthECAL < depthStep) {
263  nsteps = -1;
264  LogInfo("FastCalorimetry") << " FamosHDShower : too small ECAL and HCAL depths - "
265  << " particle is lost !!! " << std::endl;
266  }
267  }
268 
269  if (debug)
270  LogInfo("FastCalorimetry") << " FamosHDShower depths(lam) - " << std::endl
271  << " ECAL = " << depthECAL << std::endl
272  << " GAP = " << depthGAP << std::endl
273  << " HCAL = " << depthHCAL << std::endl
274  << " starting point = " << depthStart << std::endl;
275 
276  if (onEcal) {
277  if (debug)
278  LogInfo("FastCalorimetry") << " FamosHDShower : onECAL" << std::endl;
279  if (depthStart < depthECAL) {
280  if (debug)
281  LogInfo("FastCalorimetry") << " FamosHDShower : depthStart < depthECAL" << std::endl;
282  if (depthECAL > depthStep && (depthECAL - depthStart) / depthECAL > 0.1) {
283  if (debug)
284  LogInfo("FastCalorimetry") << " FamosHDShower : enough space to make ECAL step" << std::endl;
285 
286  // ECAL - one step
287  nsteps++;
288  sum1 += depthECAL; // at the end of step
289  sum2 += depthECAL - depthStart;
290  sum3 += sum2 * lambdaEM / x0EM;
291  lamtotal.push_back(sum1);
292  lamdepth.push_back(sum2);
293  lamcurr.push_back(lambdaEM);
294  lamstep.push_back(depthECAL - depthStart);
295  x0depth.push_back(sum3);
296  x0curr.push_back(x0EM);
297  detector.push_back(1);
298  mip = 0;
299 
300  if (debug)
301  LogInfo("FastCalorimetry") << " FamosHDShower : "
302  << " in ECAL sum1, sum2 " << sum1 << " " << sum2 << std::endl;
303 
304  // Gap - no additional step after ECAL
305  // just move further to HCAL over the gap
306  sum1 += depthGAP;
307  sum2 += depthGAP;
308  sum3 += depthGAPx0;
309  } else { // Just shift starting point to HCAL
310  // cout << " FamosHDShower : not enough space to make ECAL step" << std::endl;
311  if (debug)
312  LogInfo("FastCalorimetry") << " FamosHDShower : goto HCAL" << std::endl;
313 
314  depthStart = depthToHCAL;
315  sum1 += depthStart;
316  }
317  } else { // GAP or HCAL
318  if (depthStart >= depthECAL && depthStart < depthToHCAL) {
319  depthStart = depthToHCAL; // just a shift to HCAL for simplicity
320  }
321  sum1 += depthStart;
322  if (debug)
323  LogInfo("FastCalorimetry") << " FamosHDShower : goto HCAL" << std::endl;
324  }
325  } else { // Forward
326  if (debug)
327  LogInfo("FastCalorimetry") << " FamosHDShower : forward" << std::endl;
328  sum1 += depthStart;
329  transFactor = 0.5; // makes narower tresverse size of shower
330  }
331 
332  for (int i = 0; i < nmoresteps; i++) {
333  sum1 += depthStep;
334  if (sum1 > (depthECAL + depthGAP + depthHCAL))
335  break;
336  sum2 += depthStep;
337  sum3 += sum2 * lambdaHD / x0HD;
338  lamtotal.push_back(sum1);
339  lamdepth.push_back(sum2);
340  lamcurr.push_back(lambdaHD);
341  lamstep.push_back(depthStep);
342  x0depth.push_back(sum3);
343  x0curr.push_back(x0HD);
344  detector.push_back(3);
345  nsteps++;
346  }
347 
348  // Make fractions of energy and transverse radii at each step
349 
350  if (nsteps > 0) {
351  makeSteps(nsteps);
352  }
353 }

References aloge, HDShowerParametrization::alpe1(), HDShowerParametrization::alpe2(), alpEM, HDShowerParametrization::alph1(), HDShowerParametrization::alph2(), alpHD, balanceEH, HDShowerParametrization::bete1(), HDShowerParametrization::bete2(), betEM, HDShowerParametrization::beth1(), HDShowerParametrization::beth2(), betHD, criticalEnergy, debug, depthECAL, depthGAP, depthGAPx0, depthHCAL, depthStart, depthStep, depthToHCAL, detector, e, HDShowerParametrization::e1(), HDShowerParametrization::e2(), EcalHitMaker::ecalHcalGapTotalL0(), EcalHitMaker::ecalHcalGapTotalX0(), HDShowerParametrization::ecalProperties(), EcalHitMaker::ecalTotalL0(), HDShowerParametrization::emax(), HDShowerParametrization::emid(), HDShowerParametrization::emin(), eSpotSize, RandomEngineAndDistribution::flatShoot(), HSParameters::getHDbalanceEH(), HSParameters::getHDcriticalEnergy(), HSParameters::getHDdepthStep(), HSParameters::getHDeSpotSize(), HSParameters::getHDhcalDepthFactor(), HSParameters::getHDlossesOpt(), HSParameters::getHDmaxTRfactor(), HSParameters::getHDnDepthSteps(), HSParameters::getHDnTRsteps(), HSParameters::getHDtransParam(), hcalDepthFactor, HDShowerParametrization::hcalProperties(), EcalHitMaker::hcalTotalL0(), HDShowerParametrization::hsParameters(), mps_fire::i, ECALProperties::interactionLength(), HCALProperties::interactionLength(), lambdaEM, lambdaHD, lamcurr, lamdepth, lamstep, lamtotal, dqm-mbProfile::log, lossesOpt, makeSteps(), HLT_FULL_cff::maxDepth, maxTRfactor, mip, nDepthSteps, nTRsteps, onEcal, HDShowerParametrization::part1(), HDShowerParametrization::part2(), HDShowerParametrization::r1(), HDShowerParametrization::r2(), HDShowerParametrization::r3(), ECALProperties::radLenIncm(), HCALProperties::radLenIncm(), random, HDShowerParametrization::setCase(), combinedConstraintHelpers::sum2(), tgamEM, tgamHD, theECALproperties, theGrid, theHCALproperties, theParam, theR1, theR2, theR3, transFactor, transParam, x0curr, x0depth, x0EM, and x0HD.

◆ ~HDShower()

virtual HDShower::~HDShower ( )
inlinevirtual

Definition at line 41 of file HDShower.h.

41 { ; }

Member Function Documentation

◆ compute()

bool HDShower::compute ( )

Compute the shower longitudinal and lateral development.

Definition at line 453 of file HDShower.cc.

453  {
454  // TimeMe theT("FamosHDShower::compute");
455 
456  bool status = false;
457  int numLongit = eStep.size();
458  if (debug)
459  LogInfo("FastCalorimetry") << " FamosHDShower::compute - "
460  << " N_long.steps required : " << numLongit << std::endl;
461 
462  if (numLongit > 0) {
463  status = true;
464  // Prepare the trsanverse probability function
465  std::vector<double> Fhist;
466  std::vector<double> rhist;
467  for (int j = 0; j < nTRsteps + 1; j++) {
468  rhist.push_back(maxTRfactor * j / nTRsteps);
469  Fhist.push_back(transProb(maxTRfactor, 1., rhist[j]));
470  if (debug == 3)
471  LogInfo("FastCalorimetry") << "indexFinder - i, Fhist[i] = " << j << " " << Fhist[j] << std::endl;
472  }
473 
474  // Longitudinal steps
475  for (int i = 0; i < numLongit; i++) {
476  double currentDepthL0 = lamtotal[i] - 0.5 * lamstep[i];
477  // vary the longitudinal profile if needed
478  if (detector[i] != 1)
479  currentDepthL0 *= hcalDepthFactor;
480  if (debug)
481  LogInfo("FastCalorimetry") << " FamosHDShower::compute - detector = " << detector[i]
482  << " currentDepthL0 = " << currentDepthL0 << std::endl;
483 
484  double maxTRsize = maxTRfactor * rlamStep[i]; // in lambda units
485  double rbinsize = maxTRsize / nTRsteps;
486  double espot = eStep[i] / (double)nspots[i]; // re-adjust espot
487 
488  if (espot > 2. || espot < 0.)
489  LogInfo("FastCalorimetry") << " FamosHDShower::compute - unphysical espot = " << espot << std::endl;
490 
491  int ecal = 0;
492  if (detector[i] != 1) {
493  bool setHDdepth = theHcalHitMaker->setDepth(currentDepthL0);
494 
495  if (debug)
496  LogInfo("FastCalorimetry") << " FamosHDShower::compute - status of "
497  << " theHcalHitMaker->setDepth(currentDepthL0) is " << setHDdepth << std::endl;
498 
499  if (!setHDdepth) {
500  currentDepthL0 -= lamstep[i];
501  setHDdepth = theHcalHitMaker->setDepth(currentDepthL0);
502  }
503 
504  if (!setHDdepth)
505  continue;
506 
507  theHcalHitMaker->setSpotEnergy(espot);
508 
509  //fill hcal longitudinal distribution histogram
510  } else {
511  ecal = 1;
512  bool status = theGrid->getPads(currentDepthL0);
513 
514  if (debug)
515  LogInfo("FastCalorimetry") << " FamosHDShower::compute - status of Grid = " << status << std::endl;
516 
517  if (!status)
518  continue;
519 
520  int ntry = nspots[i] * 10;
521  if (ntry >= infinity) { // use max allowed in case of too many spots
522  nspots[i] = 0.5 * infinity;
523  espot *= 0.1 * (double)ntry / double(nspots[i]);
524  } else {
525  espot *= 0.1; // fine-grain energy spots in ECAL
526  // to avoid false ECAL clustering
527  nspots[i] = ntry;
528  }
529 
530  theGrid->setSpotEnergy(espot);
531 
532  //fill ecal longitudinal distribution histogram
533  }
534 
535  // Transverse distribition
536  int nok = 0; // counter of OK
537  int count = 0;
538  int inf = infinity;
539  if (lossesOpt)
540  inf = nspots[i]; // if losses are enabled, otherwise
541  // only OK points are counted ...
542  if (nspots[i] > inf)
543  std::cout << " FamosHDShower::compute - at long.step " << i << " too many spots required : " << nspots[i]
544  << " !!! " << std::endl;
545 
546  for (int j = 0; j < inf; j++) {
547  if (nok == nspots[i])
548  break;
549  count++;
550 
551  double prob = random->flatShoot();
552  int index = indexFinder(prob, Fhist);
553  double radius = rlamStep[i] * rhist[index] + random->flatShoot() * rbinsize; // in-bin
554  double phi = 2. * M_PI * random->flatShoot();
555 
556  if (debug == 2)
557  LogInfo("FastCalorimetry") << std::endl
558  << " FamosHDShower::compute "
559  << " r = " << radius << " phi = " << phi << std::endl;
560 
561  bool result;
562  if (ecal) {
563  result = theGrid->addHit(radius, phi, 0);
564 
565  if (debug == 2)
566  LogInfo("FastCalorimetry") << " FamosHDShower::compute - "
567  << " theGrid->addHit result = " << result << std::endl;
568 
569  //fill ecal transverse distribution histogram
570  } else {
571  result = theHcalHitMaker->addHit(radius, phi, 0);
572 
573  if (debug == 2)
574  LogInfo("FastCalorimetry") << " FamosHDShower::compute - "
575  << " theHcalHitMaker->addHit result = " << result << std::endl;
576 
577  //fill hcal transverse distribution histogram
578  }
579 
580  if (result)
581  nok++;
582 
583  } // end of tranverse simulation
584 
585  if (count == infinity) {
586  if (debug)
587  LogInfo("FastCalorimetry") << " FamosHDShower::compute "
588  << " maximum number of"
589  << " transverse points " << count << " is used !!!" << std::endl;
590  }
591 
592  if (debug)
593  LogInfo("FastCalorimetry") << " FamosHDShower::compute "
594  << " long.step No." << i << " Ntry, Nok = " << count << " " << nok << std::endl;
595  } // end of longitudinal steps
596  } // end of no steps
597 
598  return status;
599 }

References HcalHitMaker::addHit(), EcalHitMaker::addHit(), submitPVResolutionJobs::count, gather_cfg::cout, debug, detector, eStep, RandomEngineAndDistribution::flatShoot(), EcalHitMaker::getPads(), hcalDepthFactor, mps_fire::i, indexFinder(), dqmiodatasetharvest::inf, infinity, dqmiolumiharvest::j, lamstep, lamtotal, lossesOpt, M_PI, maxTRfactor, nspots, nTRsteps, phi, TtFullHadEvtBuilder_cfi::prob, CosmicsPD_Skims::radius, random, mps_fire::result, rlamStep, HcalHitMaker::setDepth(), HcalHitMaker::setSpotEnergy(), EcalHitMaker::setSpotEnergy(), mps_update::status, theGrid, theHcalHitMaker, and transProb().

Referenced by CalorimetryManager::HDShowerSimulation().

◆ gam()

double HDShower::gam ( double  x,
double  a 
) const
inlineprivate

Definition at line 48 of file HDShower.h.

48 { return pow(x, a - 1.) * exp(-x); }

References a, JetChargeProducer_cfi::exp, funct::pow(), and x.

Referenced by makeSteps().

◆ getmip()

int HDShower::getmip ( )
inline

Definition at line 39 of file HDShower.h.

39 { return mip; }

References mip.

Referenced by CalorimetryManager::HDShowerSimulation().

◆ indexFinder()

int HDShower::indexFinder ( double  x,
const std::vector< double > &  Fhist 
)
private

Definition at line 601 of file HDShower.cc.

601  {
602  // binary search in the vector of doubles
603  int size = Fhist.size();
604 
605  int curr = size / 2;
606  int step = size / 4;
607  int iter;
608  int prevdir = 0;
609  int actudir = 0;
610 
611  for (iter = 0; iter < size; iter++) {
612  if (curr >= size || curr < 1)
613  LogWarning("FastCalorimetry") << " FamosHDShower::indexFinder - wrong current index = " << curr << " !!!"
614  << std::endl;
615 
616  if ((x <= Fhist[curr]) && (x > Fhist[curr - 1]))
617  break;
618  prevdir = actudir;
619  if (x > Fhist[curr]) {
620  actudir = 1;
621  } else {
622  actudir = -1;
623  }
624  if (prevdir * actudir < 0) {
625  if (step > 1)
626  step /= 2;
627  }
628  curr += actudir * step;
629  if (curr > size)
630  curr = size;
631  else {
632  if (curr < 1) {
633  curr = 1;
634  }
635  }
636 
637  if (debug == 3)
638  LogInfo("FastCalorimetry") << " indexFinder - end of iter." << iter << " curr, F[curr-1], F[curr] = " << curr
639  << " " << Fhist[curr - 1] << " " << Fhist[curr] << std::endl;
640  }
641 
642  if (debug == 3)
643  LogInfo("FastCalorimetry") << " indexFinder x = " << x << " found index = " << curr - 1 << std::endl;
644 
645  return curr - 1;
646 }

References debug, findQualityFiles::size, and x.

Referenced by compute().

◆ makeSteps()

void HDShower::makeSteps ( int  nsteps)
private

Definition at line 355 of file HDShower.cc.

355  {
356  double sumes = 0.;
357  double sum = 0.;
358  std::vector<double> temp;
359 
360  if (debug)
361  LogInfo("FastCalorimetry") << " FamosHDShower::makeSteps - "
362  << " nsteps required : " << nsteps << std::endl;
363 
364  int count = 0;
365  for (int i = 0; i < nsteps; i++) {
366  double deplam = lamdepth[i] - 0.5 * lamstep[i];
367  double depx0 = x0depth[i] - 0.5 * lamstep[i] / x0curr[i];
368  double x = betEM * depx0;
369  double y = betHD * deplam;
370 
371  if (debug == 2)
372  LogInfo("FastCalorimetry") << " FamosHDShower::makeSteps "
373  << " - step " << i << " depx0, x = " << depx0 << ", " << x
374  << " deplam, y = " << deplam << ", " << y << std::endl;
375 
376  double est = (part * betEM * gam(x, alpEM) * lamcurr[i] / (x0curr[i] * tgamEM) +
377  (1. - part) * betHD * gam(y, alpHD) / tgamHD) *
378  lamstep[i];
379 
380  // protection ...
381  if (est < 0.) {
382  LogInfo("FastCalorimetry") << "*** FamosHDShower::makeSteps "
383  << " - negative step energy !!!" << std::endl;
384  est = 0.;
385  break;
386  }
387 
388  // for estimates only
389  sum += est;
390  int nPest = (int)(est * e / sum / eSpotSize);
391 
392  if (debug == 2)
393  LogInfo("FastCalorimetry") << " FamosHDShower::makeSteps - nPoints estimate = " << nPest << std::endl;
394 
395  if (nPest <= 1 && count != 0)
396  break;
397 
398  // good step - to proceed
399 
400  temp.push_back(est);
401  sumes += est;
402 
403  rlamStep.push_back(transParam * (theR1 + (theR2 - theR3 * aloge)) * deplam * transFactor);
404  count++;
405  }
406 
407  // fluctuations in ECAL and re-distribution of remaining energy in HCAL
408  if (detector[0] == 1 && count > 1) {
409  double oldECALenergy = temp[0];
410  double oldHCALenergy = sumes - oldECALenergy;
411  double newECALenergy = 2. * sumes;
412  for (int i = 0; newECALenergy > sumes && i < infinity; i++)
413  newECALenergy = 2. * balanceEH * random->flatShoot() * oldECALenergy;
414 
415  if (debug == 2)
416  LogInfo("FastCalorimetry") << "*** FamosHDShower::makeSteps "
417  << " ECAL fraction : old/new - " << oldECALenergy / sumes << "/"
418  << newECALenergy / sumes << std::endl;
419 
420  temp[0] = newECALenergy;
421  double newHCALenergy = sumes - newECALenergy;
422  double newHCALreweight = newHCALenergy / oldHCALenergy;
423 
424  for (int i = 1; i < count; i++) {
425  temp[i] *= newHCALreweight;
426  }
427  }
428 
429  // final re-normalization of the energy fractions
430  for (int i = 0; i < count; i++) {
431  eStep.push_back(temp[i] * e / sumes);
432  nspots.push_back((int)(eStep[i] / eSpotSize) + 1);
433 
434  if (debug)
435  LogInfo("FastCalorimetry") << " step " << i << " det: " << detector[i]
436  << " xO and lamdepth at the end of step = " << x0depth[i] << " " << lamdepth[i]
437  << " Estep func = " << eStep[i] << " Rstep = " << rlamStep[i]
438  << " Nspots = " << nspots[i] << " espot = " << eStep[i] / (double)nspots[i]
439  << std::endl;
440  }
441 
442  // The only step is in ECAL - let's make the size bigger ...
443  if (count == 1 and detector[0] == 1)
444  rlamStep[0] *= 2.;
445 
446  if (debug) {
447  if (eStep[0] > 0.95 * e && detector[0] == 1)
448  LogInfo("FastCalorimetry") << " FamosHDShower::makeSteps - "
449  << "ECAL energy = " << eStep[0] << " out of total = " << e << std::endl;
450  }
451 }

References aloge, alpEM, alpHD, balanceEH, betEM, betHD, submitPVResolutionJobs::count, debug, detector, e, eSpotSize, eStep, RandomEngineAndDistribution::flatShoot(), gam(), mps_fire::i, infinity, createfilelist::int, lamcurr, lamdepth, lamstep, nspots, random, rlamStep, groupFilesInBlocks::temp, tgamEM, tgamHD, theR1, theR2, theR3, transFactor, transParam, x, x0curr, x0depth, and y.

Referenced by HDShower().

◆ transProb()

double HDShower::transProb ( double  factor,
double  R,
double  r 
)
inlineprivate

Definition at line 53 of file HDShower.h.

53  {
54  double fsq = factor * factor;
55  return ((fsq + 1.) / fsq) * r * r / (r * r + R * R);
56  }

References DQMScaleToClient_cfi::factor, dttmaxenums::R, and alignCSCRings::r.

Referenced by compute().

Member Data Documentation

◆ aloge

double HDShower::aloge
private

Definition at line 76 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ alpEM

double HDShower::alpEM
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ alpHD

double HDShower::alpHD
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ balanceEH

double HDShower::balanceEH
private

Definition at line 119 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ betEM

double HDShower::betEM
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ betHD

double HDShower::betHD
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ criticalEnergy

double HDShower::criticalEnergy
private

Definition at line 115 of file HDShower.h.

Referenced by HDShower().

◆ depthECAL

double HDShower::depthECAL
private

Definition at line 127 of file HDShower.h.

Referenced by HDShower().

◆ depthGAP

double HDShower::depthGAP
private

Definition at line 127 of file HDShower.h.

Referenced by HDShower().

◆ depthGAPx0

double HDShower::depthGAPx0
private

Definition at line 127 of file HDShower.h.

Referenced by HDShower().

◆ depthHCAL

double HDShower::depthHCAL
private

Definition at line 127 of file HDShower.h.

Referenced by HDShower().

◆ depthStart

double HDShower::depthStart
private

Definition at line 75 of file HDShower.h.

Referenced by HDShower().

◆ depthStep

double HDShower::depthStep
private

Definition at line 113 of file HDShower.h.

Referenced by HDShower().

◆ depthToHCAL

double HDShower::depthToHCAL
private

Definition at line 127 of file HDShower.h.

Referenced by HDShower().

◆ detector

std::vector<int> HDShower::detector
private

Definition at line 78 of file HDShower.h.

Referenced by compute(), HDShower(), and makeSteps().

◆ e

double HDShower::e
private

Definition at line 98 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ eSpotSize

double HDShower::eSpotSize
private

Definition at line 111 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ eStep

std::vector<double> HDShower::eStep
private

Definition at line 79 of file HDShower.h.

Referenced by compute(), and makeSteps().

◆ hcalDepthFactor

double HDShower::hcalDepthFactor
private

Definition at line 121 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ infinity

int HDShower::infinity
private

Definition at line 83 of file HDShower.h.

Referenced by compute(), and makeSteps().

◆ lambdaEM

double HDShower::lambdaEM
private

Definition at line 74 of file HDShower.h.

Referenced by HDShower().

◆ lambdaHD

double HDShower::lambdaHD
private

Definition at line 74 of file HDShower.h.

Referenced by HDShower().

◆ lamcurr

std::vector<double> HDShower::lamcurr
private

Definition at line 81 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ lamdepth

std::vector<double> HDShower::lamdepth
private

Definition at line 81 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ lamstep

std::vector<double> HDShower::lamstep
private

Definition at line 81 of file HDShower.h.

Referenced by compute(), HDShower(), and makeSteps().

◆ lamtotal

std::vector<double> HDShower::lamtotal
private

Definition at line 81 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ lossesOpt

int HDShower::lossesOpt
private

Definition at line 101 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ maxTRfactor

double HDShower::maxTRfactor
private

Definition at line 117 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ mip

int HDShower::mip
private

Definition at line 95 of file HDShower.h.

Referenced by getmip(), and HDShower().

◆ nDepthSteps

int HDShower::nDepthSteps
private

Definition at line 103 of file HDShower.h.

Referenced by HDShower().

◆ nspots

std::vector<int> HDShower::nspots
private

Definition at line 78 of file HDShower.h.

Referenced by compute(), and makeSteps().

◆ nTRsteps

int HDShower::nTRsteps
private

Definition at line 105 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ onEcal

int HDShower::onEcal
private

Definition at line 92 of file HDShower.h.

Referenced by HDShower().

◆ part

double HDShower::part
private

Definition at line 71 of file HDShower.h.

◆ random

const RandomEngineAndDistribution* HDShower::random
private

Definition at line 124 of file HDShower.h.

Referenced by compute(), HDShower(), and makeSteps().

◆ rlamStep

std::vector<double> HDShower::rlamStep
private

Definition at line 79 of file HDShower.h.

Referenced by compute(), and makeSteps().

◆ tgamEM

double HDShower::tgamEM
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ tgamHD

double HDShower::tgamHD
private

Definition at line 71 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ theECALproperties

const ECALProperties* HDShower::theECALproperties
private

Definition at line 66 of file HDShower.h.

Referenced by HDShower().

◆ theGrid

EcalHitMaker* HDShower::theGrid
private

Definition at line 86 of file HDShower.h.

Referenced by compute(), and HDShower().

◆ theHcalHitMaker

HcalHitMaker* HDShower::theHcalHitMaker
private

Definition at line 89 of file HDShower.h.

Referenced by compute().

◆ theHCALproperties

const HCALProperties* HDShower::theHCALproperties
private

Definition at line 67 of file HDShower.h.

Referenced by HDShower().

◆ theParam

HDShowerParametrization* HDShower::theParam
private

Definition at line 63 of file HDShower.h.

Referenced by HDShower().

◆ theR1

double HDShower::theR1
private

Definition at line 70 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ theR2

double HDShower::theR2
private

Definition at line 70 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ theR3

double HDShower::theR3
private

Definition at line 70 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ transFactor

double HDShower::transFactor
private

Definition at line 109 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ transParam

double HDShower::transParam
private

Definition at line 107 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ x0curr

std::vector<double> HDShower::x0curr
private

Definition at line 80 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ x0depth

std::vector<double> HDShower::x0depth
private

Definition at line 80 of file HDShower.h.

Referenced by HDShower(), and makeSteps().

◆ x0EM

double HDShower::x0EM
private

Definition at line 74 of file HDShower.h.

Referenced by HDShower().

◆ x0HD

double HDShower::x0HD
private

Definition at line 74 of file HDShower.h.

Referenced by HDShower().

debug
#define debug
Definition: HDShower.cc:17
HDShower::theParam
HDShowerParametrization * theParam
Definition: HDShower.h:63
HSParameters::getHDnDepthSteps
int getHDnDepthSteps() const
Definition: HSParameters.h:20
DDAxes::y
HDShower::nspots
std::vector< int > nspots
Definition: HDShower.h:78
mps_fire.i
i
Definition: mps_fire.py:428
HDShower::lamtotal
std::vector< double > lamtotal
Definition: HDShower.h:81
HDShower::betHD
double betHD
Definition: HDShower.h:71
HDShower::theHcalHitMaker
HcalHitMaker * theHcalHitMaker
Definition: HDShower.h:89
HDShower::transFactor
double transFactor
Definition: HDShower.h:109
step
step
Definition: StallMonitor.cc:94
HDShower::gam
double gam(double x, double a) const
Definition: HDShower.h:48
HDShower::detector
std::vector< int > detector
Definition: HDShower.h:78
HDShower::aloge
double aloge
Definition: HDShower.h:76
HSParameters::getHDmaxTRfactor
double getHDmaxTRfactor() const
Definition: HSParameters.h:26
mps_update.status
status
Definition: mps_update.py:68
HDShower::lamstep
std::vector< double > lamstep
Definition: HDShower.h:81
gather_cfg.cout
cout
Definition: gather_cfg.py:144
HCALProperties::radLenIncm
double radLenIncm() const override
Radiation length in cm.
Definition: HCALProperties.h:37
HDShower::balanceEH
double balanceEH
Definition: HDShower.h:119
HDShowerParametrization::beth1
double beth1() const
Definition: HDShowerParametrization.h:81
HDShower::depthHCAL
double depthHCAL
Definition: HDShower.h:127
HDShowerParametrization::bete2
double bete2() const
Definition: HDShowerParametrization.h:63
DDAxes::x
HSParameters::getHDdepthStep
double getHDdepthStep() const
Definition: HSParameters.h:24
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
HDShower::lamdepth
std::vector< double > lamdepth
Definition: HDShower.h:81
HDShowerParametrization::part1
double part1() const
Definition: HDShowerParametrization.h:93
ecal
Definition: ElectronicsMappingGPU.h:13
HDShowerParametrization::r2
double r2() const
Definition: HDShowerParametrization.h:106
groupFilesInBlocks.temp
list temp
Definition: groupFilesInBlocks.py:142
EcalHitMaker::getPads
bool getPads(double depth, bool inCm=false)
Definition: EcalHitMaker.cc:820
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HDShowerParametrization::alph1
double alph1() const
Definition: HDShowerParametrization.h:69
HLT_FULL_cff.maxDepth
maxDepth
Definition: HLT_FULL_cff.py:8665
EcalHitMaker::addHit
bool addHit(double r, double phi, unsigned layer=0) override
Definition: EcalHitMaker.cc:184
HDShower::depthECAL
double depthECAL
Definition: HDShower.h:127
HSParameters::getHDtransParam
double getHDtransParam() const
Definition: HSParameters.h:22
HDShower::criticalEnergy
double criticalEnergy
Definition: HDShower.h:115
HDShower::theR1
double theR1
Definition: HDShower.h:70
HDShower::theR3
double theR3
Definition: HDShower.h:70
HDShowerParametrization::emax
double emax() const
Definition: HDShowerParametrization.h:41
part
part
Definition: HCALResponse.h:20
HDShowerParametrization::r1
double r1() const
Definition: HDShowerParametrization.h:105
HDShower::x0HD
double x0HD
Definition: HDShower.h:74
HDShower::eStep
std::vector< double > eStep
Definition: HDShower.h:79
HDShowerParametrization::part2
double part2() const
Definition: HDShowerParametrization.h:99
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
submitPVResolutionJobs.count
count
Definition: submitPVResolutionJobs.py:352
DQMScaleToClient_cfi.factor
factor
Definition: DQMScaleToClient_cfi.py:8
HDShower::nTRsteps
int nTRsteps
Definition: HDShower.h:105
HDShowerParametrization::setCase
void setCase(int choice)
Definition: HDShowerParametrization.h:29
HDShower::mip
int mip
Definition: HDShower.h:95
HDShower::transParam
double transParam
Definition: HDShower.h:107
HDShower::alpHD
double alpHD
Definition: HDShower.h:71
HSParameters::getHDbalanceEH
double getHDbalanceEH() const
Definition: HSParameters.h:27
HDShower::depthStep
double depthStep
Definition: HDShower.h:113
HSParameters::getHDlossesOpt
int getHDlossesOpt() const
Definition: HSParameters.h:19
HDShower::depthGAP
double depthGAP
Definition: HDShower.h:127
HDShower::tgamHD
double tgamHD
Definition: HDShower.h:71
a
double a
Definition: hdecay.h:119
HSParameters::getHDhcalDepthFactor
double getHDhcalDepthFactor() const
Definition: HSParameters.h:28
HDShower::theHCALproperties
const HCALProperties * theHCALproperties
Definition: HDShower.h:67
HDShower::tgamEM
double tgamEM
Definition: HDShower.h:71
EcalHitMaker::setSpotEnergy
void setSpotEnergy(double e) override
Definition: EcalHitMaker.h:112
HDShower::x0curr
std::vector< double > x0curr
Definition: HDShower.h:80
HDShower::lambdaEM
double lambdaEM
Definition: HDShower.h:74
HDShower::rlamStep
std::vector< double > rlamStep
Definition: HDShower.h:79
HDShower::lamcurr
std::vector< double > lamcurr
Definition: HDShower.h:81
HDShowerParametrization::hsParameters
const HSParameters * hsParameters() const
Definition: HDShowerParametrization.h:26
createfilelist.int
int
Definition: createfilelist.py:10
HDShowerParametrization::e1
double e1() const
Definition: HDShowerParametrization.h:43
HDShower::infinity
int infinity
Definition: HDShower.h:83
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
HDShower::theECALproperties
const ECALProperties * theECALproperties
Definition: HDShower.h:66
HDShower::lambdaHD
double lambdaHD
Definition: HDShower.h:74
HDShower::e
double e
Definition: HDShower.h:98
HDShower::x0EM
double x0EM
Definition: HDShower.h:74
combinedConstraintHelpers::sum2
void sum2(T &x, T y)
Definition: CombinedKinematicConstraintT.h:74
HDShowerParametrization::hcalProperties
const HCALProperties * hcalProperties() const
Definition: HDShowerParametrization.h:24
HDShower::maxTRfactor
double maxTRfactor
Definition: HDShower.h:117
HDShower::makeSteps
void makeSteps(int nsteps)
Definition: HDShower.cc:355
HDShower::eSpotSize
double eSpotSize
Definition: HDShower.h:111
RandomEngineAndDistribution::flatShoot
double flatShoot(double xmin=0.0, double xmax=1.0) const
Definition: RandomEngineAndDistribution.h:27
HDShower::x0depth
std::vector< double > x0depth
Definition: HDShower.h:80
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
HDShowerParametrization::alph2
double alph2() const
Definition: HDShowerParametrization.h:75
HDShowerParametrization::alpe1
double alpe1() const
Definition: HDShowerParametrization.h:45
DDAxes::phi
HDShower::betEM
double betEM
Definition: HDShower.h:71
HDShowerParametrization::alpe2
double alpe2() const
Definition: HDShowerParametrization.h:51
HSParameters::getHDcriticalEnergy
double getHDcriticalEnergy() const
Definition: HSParameters.h:25
ECALProperties::radLenIncm
double radLenIncm() const override
Radiation length in cm.
Definition: ECALProperties.h:32
HDShowerParametrization::ecalProperties
const ECALProperties * ecalProperties() const
Definition: HDShowerParametrization.h:22
HDShower::depthStart
double depthStart
Definition: HDShower.h:75
dqmiodatasetharvest.inf
inf
Definition: dqmiodatasetharvest.py:38
HDShower::transProb
double transProb(double factor, double R, double r)
Definition: HDShower.h:53
HDShower::hcalDepthFactor
double hcalDepthFactor
Definition: HDShower.h:121
HDShower::depthGAPx0
double depthGAPx0
Definition: HDShower.h:127
HSParameters::getHDnTRsteps
int getHDnTRsteps() const
Definition: HSParameters.h:21
EcalHitMaker::ecalHcalGapTotalX0
double ecalHcalGapTotalX0() const
ECAL-HCAL transition.
Definition: EcalHitMaker.h:70
CosmicsPD_Skims.radius
radius
Definition: CosmicsPD_Skims.py:135
HDShower::indexFinder
int indexFinder(double x, const std::vector< double > &Fhist)
Definition: HDShower.cc:601
HDShowerParametrization::bete1
double bete1() const
Definition: HDShowerParametrization.h:57
HDShower::depthToHCAL
double depthToHCAL
Definition: HDShower.h:127
HDShower::random
const RandomEngineAndDistribution * random
Definition: HDShower.h:124
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
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
HDShowerParametrization::e2
double e2() const
Definition: HDShowerParametrization.h:44
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
HDShower::alpEM
double alpEM
Definition: HDShower.h:71
HDShowerParametrization::emin
double emin() const
Definition: HDShowerParametrization.h:37
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
HSParameters::getHDeSpotSize
double getHDeSpotSize() const
Definition: HSParameters.h:23
HDShowerParametrization::emid
double emid() const
Definition: HDShowerParametrization.h:39
EcalHitMaker::hcalTotalL0
double hcalTotalL0() const
in the HCAL
Definition: EcalHitMaker.h:88
HDShower::lossesOpt
int lossesOpt
Definition: HDShower.h:101
HcalHitMaker::setSpotEnergy
void setSpotEnergy(double e) override
Set the spot energy.
Definition: HcalHitMaker.h:26
HDShower::theGrid
EcalHitMaker * theGrid
Definition: HDShower.h:86
dttmaxenums::R
Definition: DTTMax.h:29
HDShowerParametrization::r3
double r3() const
Definition: HDShowerParametrization.h:107
TtFullHadEvtBuilder_cfi.prob
prob
Definition: TtFullHadEvtBuilder_cfi.py:33
HDShower::theR2
double theR2
Definition: HDShower.h:70
HDShower::onEcal
int onEcal
Definition: HDShower.h:92
HDShower::nDepthSteps
int nDepthSteps
Definition: HDShower.h:103
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
HDShowerParametrization::beth2
double beth2() const
Definition: HDShowerParametrization.h:87
Generated for CMSSW Reference Manual by doxygen 1.8.16