|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
#ifndef HFShower_H
#define HFShower_H
//FastSimulation Headers
#include "FastSimulation/ShowerDevelopment/interface/HDShowerParametrization.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include <vector>
/**
* \author Salavat Abdullin
* \date: 21-Oct-2004
* \parameterized hadronic shower simulation
* \based on paper G.Gridhammer, M.Rudowicz, S.Peters, SLAC-PUB-5072
*/
class EcalHitMaker;
class HcalHitMaker;
class RandomEngineAndDistribution;
class HFShower {
public:
typedef math::XYZVector XYZPoint;
typedef std::pair<XYZPoint, double> Spot;
typedef std::pair<unsigned int, double> Step;
typedef std::vector<Step> Steps;
typedef Steps::const_iterator step_iterator;
HFShower(const RandomEngineAndDistribution* engine,
HDShowerParametrization* myParam,
EcalHitMaker* myGrid,
HcalHitMaker* myHcalHitMaker,
int onECAL,
double epart);
virtual ~HFShower() { ; }
/// Compute the shower longitudinal and lateral development
bool compute();
private:
// The longitudinal development ersatzt.
double gam(double x, double a) const { return pow(x, a - 1.) * exp(-x); }
// Transverse integrated probability function (for 4R max size)
// integral of the transverse ansatz f(r) = 2rR/(r**2 + R**2)**2 ->
// 1/R - R/(R**2+r**2) | at limits 4R - 0
double transProb(double factor, double R, double r) {
double fsq = factor * factor;
return ((fsq + 1.) / fsq) * r * r / (r * r + R * R);
}
// Compute longE[i] and transR[i] for all nsteps
void makeSteps(int nsteps);
int indexFinder(double x, const std::vector<double>& Fhist);
// The parametrization
HDShowerParametrization* theParam;
// The Calorimeter properties
const ECALProperties* theECALproperties;
const HCALProperties* theHCALproperties;
// Basic parameters of the simulation
double theR1, theR2, theR3;
double alpEM, betEM, alpHD, betHD, part, tgamEM, tgamHD;
// The basic quantities for the shower development.
double lambdaEM, lambdaHD, x0EM, x0HD;
double depthStart;
double aloge;
std::vector<int> detector, nspots;
std::vector<double> eStep, rlamStep;
std::vector<double> x0curr, x0depth;
std::vector<double> lamstep, lamcurr, lamdepth, lamtotal;
int infinity; // big number of cycles if exit is on a special condition
// The crystal grid
EcalHitMaker* theGrid;
// The HCAL
HcalHitMaker* theHcalHitMaker;
// OnECAL flag as an input parameter ...
int onEcal;
// Input energy to distribute
double e;
// HCAL losses option (0-off, 1-on)
int lossesOpt;
// Number of longitudinal steps in HCAL
int nDepthSteps;
// Number of bins in the transverse probability histogram
int nTRsteps;
// Transverse size tuning factor
double transParam;
// Transverse normalization : 1 for HB/HE, 0.5 for HF (narrow showers)
double transFactor;
// HCAL energy spot size
double eSpotSize;
// Longitudinal step size (lambda units)
double depthStep;
// Energy threshold (one depth step -> nDepthSteps);
double criticalEnergy;
// Transverse size cut (in character transverse size units)
double maxTRfactor;
// Balance between ECAL and HCAL "visible" energy (default = 1.)
double balanceEH;
// Regulator of HCAL depth of the shower (to adjust/shrink it to CMS depth)
double hcalDepthFactor;
// Famos Random Engine
const RandomEngineAndDistribution* random;
};
#endif
|