CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GflashNameSpace.h
Go to the documentation of this file.
1 #ifndef GflashNameSpace_H
2 #define GflashNameSpace_H
3 
4 #include "SimGeneral/GFlash/interface/Gflash3Vector.h"
5 
6 namespace Gflash {
7 
8  //unit convention: energy in GeV, length in cm, time in nanosecond
9 
10  enum CalorimeterNumber {
11  kNULL = -1, // ECAL Barrel - ESPM
12  kESPM, // ECAL Barrel - ESPM
13  kHB, // HCAL Barrel - HB
14  kENCA, // ECAL Endcap - ENCA
15  kHE, // HCAL Endcap - HE
16  kHO, // HCAL Outer - HO
17  kHF, // HCAL Forward - HF
18  kNumberCalorimeter
19  };
20 
21  CalorimeterNumber getCalorimeterNumber(const Gflash3Vector& position);
22  double rhoBackEB(const Gflash3Vector& position);
23  double zBackEE(const Gflash3Vector& position);
24  int findShowerType(const Gflash3Vector& position);
25 
26  // EB HB EE HE HO, HF
27  const double Zmin[kNumberCalorimeter] = {0.000, 0.000, 317.2, 391.95, 0.000, 1110.0}; // in cm
28  const double Zmax[kNumberCalorimeter] = {317.0, 433.2, 399.1, 554.10, 661.0, 1275.0};
29  const double Rmin[kNumberCalorimeter] = {123.8, 177.5, 31.6, 31.6, 382.0, 12.5};
30  const double Rmax[kNumberCalorimeter] = {177.5, 287.7, 171.1, 263.9, 407.0, 140.0};
31 
32  const double EtaMin[kNumberCalorimeter] = {0.000, 0.000, 1.479, 1.479, 0.000, 3.000};
33  const double EtaMax[kNumberCalorimeter] = {1.479, 1.479, 3.000, 3.000, 1.262, 5.000};
34 
35  const double RFrontCrystalEB = 129.0; // rho in cm
36  const double ZFrontCrystalEE = 320.9; // z in cm
37  const double ROffCrystalEB = 5.2; // rho in cm
38  const double ZOffCrystalEE = 3.7; // z in cm
39  const double LengthCrystalEB = 23.0;
40  const double LengthCrystalEE = 22.0;
41 
42  //constants needed for GflashHadronShowerModel
43  const double energyCutOff = 1.0;
44 
45  //cut value for quasi-elastic like interactions on the secondary energy,
46  //ratio = (Leading Energy of Secondaries)/(Total Energy)
47  const double QuasiElasticLike = 0.95;
48 
49  //minimum distance to the back of parameterized envelopes in [cm]
50  const double MinDistanceToOut = 10.0;
51 
52  //constants needed for GflashHadronShowerProfile
53 
54  const double criticalEnergy = 0.0086155;
55  const double rMoliere[kNumberCalorimeter] = { 2.19, 2.19, 2.19, 2.19, 2.19, 1.72}; // in cm
56  const double radLength[kNumberCalorimeter] = { 0.89, 1.49, 0.89, 1.49, 1.49, 1.76}; // in cm
57  const double intLength[kNumberCalorimeter] = { 22.4,16.42, 22.4, 16.42, 16.42, 16.77}; // in cm
58  const double Z[kNumberCalorimeter] = { 68.36, 68.36, 68.36, 68.36, 68.36, 55.845}; // mass of molicule
59 
60  // The step size of showino along the helix trajectory in cm unit
61  const double divisionStep = 1.0;
62 
63  //maximum shower depth for the hadronic lateral in interaction length unit
64  const double maxShowerDepthforR50 = 6.0;
65 
66  //maximum size of the lateral arm for hadronic showers = 12 * interaction length
67  const double maxLateralArmforR50 = 197.0;
68 
69  //minimum particle energy to trigger HO parameterization in [GeV]
70  const double MinEnergyCutOffForHO = 2.5;
71 
72  //additional energy scale for the Hcal sensitive detector
73  const double scaleSensitiveHB = 0.180;
74  // const double scaleSensitiveHE = 0.220;
75  const double scaleSensitiveHE = 0.180;
76 
77  // properties for different sub-detectors (calorimeters)
78 
79  const double SAMHAD[3][kNumberCalorimeter] = {{0.0,0.89,0.0,0.89,0.89,0.0},
80  {0.0,0.00,0.0,0.00,0.00,0.0},
81  {0.0,0.00,0.0,0.00,0.00,0.0}};
82  const double RLTHAD[kNumberCalorimeter] = {32.7,23.7,32.7,23.7,23.7,23.7};
83 
84 
85  //parameters for logitudinal/lateral parameterizations and energy responses
86 
87  const int NPar = 5;
88  const int Nrpar = 4;
89 
90  //Pion(-/+) parameters
91 
92  const double pion_emscale[2][5] = {
93  { 8.4338e-01, 1.8663e-01, -1.1906e+00, 6.8340e-01, -3.1237e-02 },
94  { 1.8918e-01, 2.2620e-02, 4.2555e+00, 1.8317e+00, -6.1358e-03 }
95  };
96 
97  const double pion_hadscale[8][5] = {
98  //based on LogNormal of Had Energy - RegionI
99  { -4.4752e+00, 3.8163e+00, 2.7962e-01, -7.6887e-01, 9.6721e-01 },
100  { 1.1228e+00, 6.5039e-01, 4.5585e+00, -8.0282e+00, -1.7630e-01 },
101  { 2.8223e+00, 5.9029e+00, -1.5873e-01, -8.7078e-01, 3.8598e-01 },
102  { -8.2726e-01, 1.0359e-01, 2.2052e+01, 2.9719e+00, 6.4103e-02 },
103  //LogNormal of Had Energy for mip-like response
104  { -8.1787e-02, 8.3243e+00, 3.7488e-03, 3.2653e+00, 9.5662e-01 },
105  { 1.0567e-01, -2.0485e-01, 2.3345e-01, 4.5563e-01, 2.3694e-02 },
106  //correction of response if ssp is between Ecal and Hcal
107  { 1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00 },
108  { 5.3155e-02, 9.3601e-02, 8.8143e-01, 2.3251e+00, -2.1999e-02 }
109  };
110 
111  const double pion_correl_hadem[5] =
112  { -7.3137e-01, 1.2367e-01, -8.5607e-01, 1.0989e+00, -1.2858e-02 };
113 
114  //pion longitudinal parameters
115 
116  const double pion_par[8*NPar][5] = {
117  { 4.4600e-01, 9.8040e-02, 1.4319e+00, 2.3920e+00, 5.8971e-02 }, //Mean Ecal ShowerType1
118  { 3.2392e+00, 5.4986e+00, 1.5979e-01, 3.2165e+00, -6.7426e-01 },
119  { -2.7502e-01, -2.2450e+00, -2.0636e-01, 2.7225e+00, -3.2656e-01 },
120  { 5.3763e-01, 3.1128e-01, 1.5000e+00, 2.7038e+00, 9.7550e-02 },
121  { 5.8269e+00, 5.9630e+00, 2.3568e-01, 3.2485e+00, -1.2201e+00 },
122 
123  { 1.6324e+00, 1.1370e+00, 3.0111e-01, 4.4501e+00, -2.5984e-01 }, //Mean Hcal ShowerType1
124  { 1.2732e+00, -3.8571e-01, 1.2127e+00, 9.7707e-01, -1.4888e-02 },
125  { -4.0078e-01, -1.0428e+00, 6.6820e-01, 2.3706e-01, -1.3898e-01 },
126  { 1.9612e+00, -1.4451e+00, 5.2179e-01, 7.7221e-01, 4.1455e-02 },
127  { 2.8888e+00, -3.2074e+00, 3.4082e-01, 3.3507e-01, 2.7346e-02 },
128 
129  { 3.5029e+00, 5.5470e+00, 1.5174e-01, 3.9215e+00, -7.5815e-01 }, //Mean Hcal ShowerType2
130  { 1.2035e+00, -4.3918e-01, 9.8431e-01, 5.3485e-01, -3.4481e-02 },
131  { 2.4186e+01, 2.5653e+01, -6.8909e-01, -2.7461e+00, -1.5702e-01 },
132  { 1.0827e+01, 1.0659e+01, -3.3366e-01, -3.1979e+00, 1.1158e-01 },
133  { 1.5010e+01, -1.5988e+01, 2.2580e-01, -4.3233e+00, 1.1448e-01 },
134 
135  { 4.0761e-01, 5.0999e-02, 1.7222e+00, 2.1085e+00, 5.4258e-02 }, //Mean Hcal ShowerType3
136  { 1.0394e+01, -9.5700e+00, 1.1167e+00, -1.6884e+00, -3.3264e-02 },
137  { -1.1882e+00, 1.6460e-01, -2.5897e+00, 5.4484e-01, -1.5845e-01 },
138  { 2.7227e+01, 2.6921e+01, -3.5664e-01, -4.8584e+00, 1.1698e-01 },
139  { 2.1544e+01, -2.2048e+01, 2.5227e-01, -5.2186e+00, 6.7578e-02 },
140 
141  { 3.6152e-02, 1.5236e-01, -5.4226e-01, 3.4917e+00, 3.1429e-02 }, //Sigma Ecal ShowerType1
142  { 4.9213e-01, 1.0093e-01, -7.8460e-01, 2.2417e+00, -1.9762e-02 },
143  { 3.6555e-01, 2.2955e-02, -3.3972e+01, 2.7112e+00, -2.0412e-02 },
144  { 7.1526e-01, 3.1537e-01, -4.8117e-01, 1.3809e+00, -1.1832e-02 },
145  { 5.3792e-01, 1.7556e-01, -5.9724e-01, 3.9307e-01, -2.1343e-02 },
146 
147  { 2.7662e-01, 5.6892e-02, -1.0211e+00, 1.5592e+00, -2.0545e-02 }, //Sigma Hcal ShowerType1
148  { 2.1221e+00, 1.7139e+00, 2.9685e-01, 3.0293e+00, -4.9449e-01 },
149  { 9.8715e-01, 1.7377e-01, -1.5715e+00, 1.1948e+00, -7.4157e-02 },
150  { 1.2424e+00, 5.8446e-01, -7.3994e-01, 1.6481e+00, -6.4211e-02 },
151  { 1.2797e+00, 9.5722e-01, -4.1961e-01, 1.5458e+00, -1.9276e-02 },
152 
153  { 5.6503e-01, 3.2484e-01, -6.8870e-01, -9.6941e-01, -2.4937e-02 }, //Sigma Hcal ShowerType2
154  { 7.2561e-01, 5.3369e-02, 9.1553e+00, 1.9652e+00, -5.4515e-02 },
155  { 3.3058e+01, 8.7333e+01, -2.4388e-02, -1.5815e+01, 1.6538e+00 },
156  { 1.1171e+00, 6.8039e-01, -6.0860e-01, 1.2764e+00, -2.1545e-02 },
157  { 2.2385e+00, 2.4309e+00, -2.4251e-01, -8.3314e-01, 4.4941e-02 },
158 
159  { 2.5595e-01, 2.0105e-02, -1.7804e+00, 2.6746e+00, -2.6966e-02 }, //Sigma Hcal ShowerType3
160  { 3.0882e-01, 3.1336e-01, -5.1932e-01, 3.3669e+00, 6.9642e-02 },
161  { 6.1781e-01, 1.1997e-01, 9.2596e-01, -1.3085e+01, -5.1192e-02 },
162  { 1.4680e+00, -1.1644e+00, 4.1119e-01, -2.3604e-01, 3.1277e-03 },
163  { 3.9081e+00, 4.2913e+00, -1.7321e-01, -3.9169e+00, 6.0504e-02 }
164  };
165 
166  const double pion_rho[8*NPar][5] = {
167  { 1.2318e-01, 3.0000e-01, 8.8725e-01, 1.6928e+00, 0.0000e+00 }, //Ecal ShowerType1
168  { 4.1546e-01, 1.6375e-01, 2.0719e+00, 1.2931e+00, -2.3695e-03 },
169  { -9.4588e-01, 1.7680e+00, 1.1180e+00, -1.7623e+00, 6.6812e-03 },
170  { -1.8668e-01, 1.5000e-01, -5.0000e+00, 2.4000e+00, 6.0000e-02 },
171  { 1.1945e-01, 1.5393e-01, -1.6983e+00, 2.2457e+00, 5.1813e-02 },
172  { 3.4322e-02, 2.0564e-01, -1.4297e+00, 2.1601e+00, 4.0797e-02 },
173  { -1.0894e-01, 1.4305e-01, -1.5543e+00, 2.5491e+00, 5.6509e-02 },
174  { 8.1763e-02, 8.0722e-02, -5.8023e+00, 2.4338e+00, 5.4228e-02 },
175  { 9.4796e-02, 1.5509e-01, -1.8610e+00, 2.3388e+00, 3.9771e-02 },
176  { 7.8380e-01, 1.4165e-01, 7.7385e-01, 1.9248e+00, 4.2539e-03 },
177 
178  { 4.2393e-01, 2.7909e-01, 9.2577e-01, 2.7916e+00, -6.8012e-02 }, //Hcal ShowerType1
179  { 5.6826e-01, 2.9822e-01, 8.1516e-01, 2.3029e+00, -9.2674e-02 },
180  { 8.0676e-01, 1.0045e+00, -4.2949e-03, -1.1090e+01, -9.6997e-04 },
181  { -1.8671e-01, 7.4935e-01, 3.6716e-04, 4.3540e+00, 8.4585e-02 },
182  { 3.7378e-01, 1.8905e+00, 1.3206e-01, 8.3958e-01, -1.4278e-01 },
183  { 2.4541e-01, 1.4579e-01, 1.3386e+00, 1.1209e+00, 1.5902e-03 },
184  { 4.4394e-02, 1.1718e-01, 2.0436e+00, 1.0125e+00, -2.0101e-02 },
185  { 1.5076e-01, 8.8556e-02, -1.8958e+00, 2.2495e+00, 4.6895e-02 },
186  { 5.0817e-01, 1.2550e-01, 2.3684e+00, 1.0669e+00, -1.3888e-02 },
187  { 8.3803e-01, 5.6501e-02, -1.2531e+00, 2.7326e+00, -4.1174e-02 },
188 
189  { 3.0277e-01, 2.1400e-01, 1.2867e+00, 2.3086e+00, -8.7281e-02 }, //Hcal ShowerType2
190  { 5.2508e-01, 2.6391e-01, 1.3357e+00, 1.9172e+00, -9.2390e-02 },
191  { 7.5387e-01, -3.2142e-02, 4.7437e+00, 8.5346e-01, -1.1566e-02 },
192  { -7.5633e-02, 8.8550e-02, 2.2340e+00, 9.7441e-01, 3.2514e-02 },
193  { 2.0314e-01, 1.2932e-01, 9.9352e-01, 6.6276e-01, 3.2726e-02 },
194  { -2.6996e+00, 3.9944e+00, -2.3836e-01, 3.4896e+00, 8.5793e-01 },
195  { 1.2888e-01, 2.4009e-01, 9.8654e-01, 9.8778e-01, -4.3265e-02 },
196  { 2.2281e-01, 8.4319e-02, 4.8402e-04, -2.0111e-01, 1.9263e-02 },
197  { -9.3521e+00, 1.0283e+01, 3.5236e-01, -4.8949e+00, -7.7675e-02 },
198  { 7.6844e-01, -9.0198e-02, 1.0364e+00, 2.3443e+00, -1.4107e-02 },
199 
200  { 3.0154e-01, 3.6790e-01, 7.5730e-01, 1.7676e+00, -9.7870e-02 }, //Hcal ShowerType3
201  { 4.2110e-01, 3.6328e-01, 9.7741e-01, 1.1211e+00, -6.4038e-02 },
202  { 7.6092e-01, 3.5780e-02, 9.7086e-01, 1.9511e+00, -1.2363e-02 },
203  { -1.4874e+00, 2.4273e+00, -1.0552e-01, 7.3591e+00, 1.7302e-01 },
204  { -1.6966e+00, 2.9215e+00, -2.2195e-01, 3.3990e+00, 6.0529e-01 },
205  { 1.6144e-01, 2.1615e-01, -7.4824e-01, 3.1003e+00, 0.0000e+00 },
206  { -2.2588e+00, 3.9894e+01, 5.5166e-02, -1.1180e+00, -2.1166e+00 },
207  { 6.8954e-01, 1.1032e+00, 2.4371e-01, 1.8565e+00, -1.9127e-01 },
208  { 4.7802e-02, 3.0961e+00, 1.4566e-01, -1.0462e+00, -3.6505e-01 },
209  { 6.3815e-01, -5.1494e-01, 3.3634e-01, 1.4517e+00, 8.8357e-02 }
210  };
211 
212  //pion lateral parameters
213 
214  const double pion_rpar[4*Nrpar][5] = {
215  { 1.7239e+00, 9.3646e-01, -7.6784e-01, 1.6319e+00, 0.0000e+00 }, //Ecal for showerType = 1
216  { 8.0152e+00, 1.3573e+00, -2.2490e+00, 2.0799e+00, -9.1649e-01 },
217  { 7.3494e-01, 8.7912e-02, -1.9433e+00, 1.2044e+00, -3.4421e-02 },
218  { -2.2025e-01, 7.3525e-02, 3.1836e+00, 1.4077e+00, 2.3818e-02 },
219 
220  { 2.0617e+01, 2.1354e+00, -1.3979e+00, 2.4424e+00, -2.0928e+00 }, //Hcal for showerType = 1
221  { 1.4732e+02, 1.4804e+02, -1.0534e+00, -2.1180e+00, 1.9513e-02 },
222  { 5.2939e-01, 2.3488e-01, 4.2874e-01, 1.4435e+00, -6.7974e-02 },
223  { -2.2405e-02, 1.7781e-01, 5.8629e-03, -2.6258e+01, 2.0595e-03 },
224 
225  { 1.3184e+01, 5.3308e+00, -6.5926e-01, 1.4287e+00, -9.1948e-01 }, //Hcal for showerType = 2
226  { 6.0658e+00, 7.6593e+00, -6.9664e-01, -1.7824e+00, 1.5619e-01 },
227  { 5.5211e-01, 0.0000e+00, 0.0000e+00, 0.0000e+00, -2.9011e-02 },
228  { 2.1551e-02, -5.4231e-02, 1.6615e+00, 2.7025e-01, 1.1807e-02 },
229 
230  { 9.0230e-01, 8.7071e-01, -5.9984e-01, 2.7973e+00, 2.4553e-01 }, //Hcal for showerType = 3
231  { 7.5508e+00, 3.6556e+00, -8.3028e-01, 1.1326e+00, -6.4124e-01 },
232  { 2.5737e+01, 2.5165e+01, -5.6146e-01, -4.0345e+00, -2.9592e-02 },
233  { -3.1200e-01, 2.0000e-01, 2.8915e+00, 7.5519e-01, 1.5230e-02 }
234  };
235 
236  //Proton paramters
237 
238  const double pro_emscale[2][5] = {
239  { 6.2614e-01, 1.7633e-01, -1.2677e+00, 1.3193e+00, 2.4973e-03 },
240  { 1.6413e-01, 5.2570e-03, 1.8876e+02, 1.2887e+00, -4.0266e-03 }
241  };
242 
243  const double pro_hadscale[8][5] = {
244  //based on LogNormal of Had Energy
245  { -1.2767e+02, 1.2691e+02, 2.6892e-01, -8.4621e+00, 9.7956e-01 },
246  { 1.2283e+00, 1.6307e-01, -3.6863e+00, 1.8677e+00, -4.5990e-02 },
247  { 1.8917e+00, 2.0947e+00, -3.6371e-01, 3.7861e-01, 8.8760e-02 },
248  { -7.6830e-01, 3.1319e-01, 8.2844e-01, 1.5865e+00, 2.4735e-02 },
249  //LogNormal of Had Energy for mip-like response
250  { -1.0927e+01, 1.0727e+01, 2.9922e+00, -8.2073e-01, 9.8722e-01 },
251  { 7.2797e-02, 2.3601e-02, -1.4900e+00, 1.8989e+00, 3.4294e-04 },
252  //correction of response if ssp is between Ecal and Hcal
253  { 1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00 },
254  { 4.7858e-01, 5.8552e-01, 3.9112e-01, 2.6945e+00, -1.7791e-01 }
255  };
256 
257  const double pro_correl_hadem[5] =
258  { -6.7711e-01, 2.5760e-01, -4.9826e-01, 3.4791e-01, 0.0000e+00 };
259 
260  //proton longitudinal parameters
261 
262  const double pro_par[8*NPar][5] = {
263  { 3.5805e-01, 8.2073e-02, 2.4040e+00, 2.8935e+00, 7.1200e-02 }, //Mean Ecal ShowerType1
264  { 8.1351e-01, 9.3288e-02, 4.5959e+00, 4.1930e+00, 8.2825e-02 },
265  { -1.2342e+00, 1.3826e-01, 8.1966e-01, 3.7141e+00, 2.7493e-02 },
266  { 5.6140e-01, 2.1720e-01, 2.0061e+00, 2.8141e+00, 1.1165e-01 },
267  { 3.9409e+00, 2.6083e+00, 3.0382e-01, 3.2679e+00, -6.3072e-01 },
268 
269  { 1.5450e+00, 1.1047e+00, 2.9732e-01, 4.5325e+00, -2.4096e-01 }, //Mean Hcal ShowerType1
270  { 1.5098e+00, -6.8331e-01, 8.7004e-01, 5.0915e-01, -9.3536e-03 },
271  { 6.9904e+00, -8.9057e+00, 3.8745e-01, -2.9089e+00, -6.0353e-02 },
272  { 1.8465e+00, -1.1548e+00, 7.1662e-01, 9.2467e-01, 7.0536e-03 },
273  { 2.1807e+00, -1.4464e+00, 7.2033e-01, 1.0696e+00, -1.3923e-01 },
274 
275  { 3.9185e+00, 6.1651e+00, 1.4738e-01, 4.2079e+00, -8.1423e-01 }, //Mean Hcal ShowerType2
276  { 8.3119e+00, 7.5628e+00, -4.9561e-01, -2.9352e+00, -3.0061e-02 },
277  { 2.7150e+00, 4.5234e+00, -4.6494e-01, -1.9846e+00, -9.9347e-02 },
278  { 2.2762e+00, 1.7668e+00, -5.3421e-01, -6.4613e-03, 5.3533e-02 },
279  { 2.6187e+00, -2.3479e+00, 4.9602e-01, 1.8344e-01, -5.3695e-02 },
280 
281  { 1.9666e+00, 3.0810e+00, 1.7874e-01, 3.2441e+00, -4.4566e-01 }, //Mean Hcal ShowerType3
282  { 1.6223e+01, 1.7175e+01, -1.3214e-01, -1.0103e+01, 1.9030e-01 },
283  { 8.1636e+01, 8.3368e+01, -5.0505e-01, -5.2201e+00, -1.0340e-01 },
284  { 5.5674e+01, 5.5241e+01, -3.8856e-01, -5.5563e+00, 9.2886e-02 },
285  { 6.1932e+01, 6.2078e+01, -3.2123e-01, -5.9991e+00, 2.2640e-02 },
286 
287  { 7.9984e-02, 1.0177e-01, -4.9195e-01, 3.6549e+00, 1.3986e-02 }, //Sigma Ecal ShowerType1
288  { 4.6155e-01, 1.1266e-01, -1.1120e+00, 2.8973e+00, -7.7394e-03 },
289  { 3.6587e-01, 3.0463e-02, -1.8552e+00, 2.9503e+00, -1.9858e-02 },
290  { 7.0777e-01, 1.5344e-01, -9.7573e-01, 2.2658e+00, -4.3013e-02 },
291  { 3.7488e-01, 3.3388e-01, -3.2414e-01, 1.7099e+00, 2.7456e-02 },
292 
293  { 2.3330e-01, 8.3848e-02, -1.4841e+00, 1.6539e+00, -8.9738e-03 }, //Sigma Hcal ShowerType1
294  { 9.0862e-01, 1.9284e-01, -2.9385e+00, 5.5159e-01, -4.4707e-02 },
295  { 1.0615e+00, 4.0817e-01, -9.2562e-01, 7.7349e-01, -5.2196e-02 },
296  { 1.3044e+00, 6.9523e-01, -8.0370e-01, 1.3130e+00, -6.0052e-02 },
297  { 1.3648e+00, 6.7204e-01, -7.7137e-01, 1.3261e+00, -8.1459e-02 },
298 
299  { 2.3378e-01, 8.4699e-02, -1.2727e+00, 1.3170e+00, -9.0939e-03 }, //Sigma Hcal ShowerType2
300  { 6.7311e-01, 2.6483e-02, -2.3889e+01, 3.3071e+00, -3.8701e-02 },
301  { 1.9157e+01, 1.8483e+01, -6.2208e-01, -3.5409e+00, -6.3584e-02 },
302  { 1.2720e+00, 7.7241e-01, -6.6997e-01, 8.5122e-01, -3.7883e-02 },
303  { 3.2209e+00, 2.8976e+00, -3.4428e-01, -1.5930e+00, -2.2328e-02 },
304 
305  { 2.3064e+00, 2.1994e+00, -2.4570e-01, -6.2275e+00, -6.5378e-03 }, //Sigma Hcal ShowerType3
306  { 6.2645e-01, 1.0763e-01, -2.4655e+00, 3.0694e+00, -2.2801e-02 },
307  { 3.9143e+00, -3.1821e+00, 2.2346e+00, -8.5515e-01, -7.0756e-02 },
308  { 8.8440e+00, 8.6952e+00, -2.8769e-01, -4.4932e+00, 2.0374e-02 },
309  { 7.1745e+00, 7.0713e+00, -2.5061e-01, -4.6686e+00, 1.0573e-02 }
310  };
311 
312  const double pro_rho[8*NPar][5] = {
313  { 1.9717e+00, 6.6780e+00, 1.4425e-01, 2.2917e+00, -8.1369e-01 }, //Ecal ShowerType1
314  { -6.3565e+00, 7.0824e+00, 1.9663e-01, -7.8033e+00, -2.3854e-02 },
315  { 8.0629e-01, 6.5672e-02, 1.4501e+00, 3.3136e-01, -5.6348e-03 },
316  { -1.6000e-01, 2.3904e-01, 1.1818e+00, 6.1474e-01, -1.0088e-02 },
317  { 1.0966e-01, 1.8834e-01, -1.4410e+00, 2.6445e+00, 9.0217e-02 },
318  { 1.0565e-02, 2.4379e-01, -1.3677e+00, 2.7338e+00, 7.6137e-02 },
319  { -3.1898e-02, 1.1508e-01, 1.5916e+00, 6.2431e-01, -9.8921e-03 },
320  { 1.5641e-01, 1.8429e-01, 1.8207e+00, 4.9841e-01, 1.3907e-03 },
321  { -1.5532e-01, 3.4349e-01, -8.5195e-01, 2.8581e+00, 1.4698e-01 },
322  { 6.9661e-01, 2.0247e-01, 7.0948e-01, 6.1266e-01, 1.0035e-02 },
323 
324  { 2.8900e-01, 1.9972e-01, 1.3997e+00, 2.8253e+00, -3.0914e-02 }, //Hcal ShowerType1
325  { 3.2869e-01, 8.6843e-02, 6.0682e+01, 2.1983e+00, -1.0146e-02 },
326  { 7.8200e-01, 6.8963e-02, -1.9113e+00, 7.5866e-01, 5.9595e-03 },
327  { 9.5356e-01, 1.1469e+00, 3.0959e-01, 1.0238e+01, 7.6437e-02 },
328  { 3.4378e-01, 1.6417e-01, 1.2061e+00, 1.6405e+00, 2.3441e-02 },
329  { 2.8005e-01, 2.0804e-01, 1.3118e+00, 1.1119e+00, -1.9884e-02 },
330  { -2.7624e-02, 1.0661e-01, 2.7942e+00, 5.8531e-01, -8.9599e-03 },
331  { 2.3749e-01, 2.6046e-02, 1.8467e+01, 6.2294e-01, 9.8564e-03 },
332  { 4.6006e-01, 2.0859e-01, 1.2562e+00, 7.6745e-01, -1.9148e-02 },
333  { 8.1573e-01, 6.8433e-02, -2.1358e+00, 2.2568e+00, -3.8336e-02 },
334 
335  { 1.1195e-01, 1.1879e-01, 2.9551e+00, 1.9034e+00, -2.0542e-02 }, //Hcal ShowerType2
336  { 4.1285e-01, 4.6288e-01, 5.0506e-01, 1.4100e+00, -1.2918e-01 },
337  { 7.7938e-01, -4.7534e-02, 1.8256e+01, 6.8910e-02, -8.8481e-03 },
338  { 2.1819e-01, 2.1512e-01, 7.0238e-01, 3.6517e+00, -2.4946e-02 },
339  { 3.5413e-01, 1.8842e-01, 1.0927e+00, 1.2173e+00, 3.4229e-04 },
340  { -1.9019e+00, 2.8987e+00, -2.5328e-01, 3.4091e+00, 6.5499e-01 },
341  { -7.1433e-02, 7.6591e-02, -6.8698e+00, 2.1165e+00, 5.6828e-02 },
342  { 2.2651e-01, 7.3121e-01, 6.9897e-04, 3.5656e+00, 2.7391e-02 },
343  { -9.1867e+00, 1.0087e+01, 3.4287e-01, -5.2458e+00, -7.1044e-02 },
344  { 7.7621e-01, -9.5476e-02, 1.0892e+00, 1.8212e+00, -1.6587e-02 },
345 
346  { 1.4812e-01, 3.7207e-01, 6.7394e-01, 1.6932e+00, -8.9577e-02 }, //Hcal ShowerType3
347  { -1.0980e+01, 1.1814e+01, 2.9231e-01, -5.5518e+00, -9.7033e-02 },
348  { 7.8865e-01, -4.6499e-02, 1.0071e+01, -1.7465e-01, -7.8028e-03 },
349  { -4.8628e-03, 6.8062e-02, 4.0379e-04, -2.9541e+00, 1.4925e-02 },
350  { 3.2583e-01, 2.2048e-01, 1.3796e+00, 7.7019e-01, -1.9870e-02 },
351  { -4.1067e+00, 8.5741e+00, -1.8161e-01, 3.0661e+00, 1.4121e+00 },
352  { 1.3314e-01, 1.3727e-01, -2.5181e-04, -4.8352e-01, -7.1482e-03 },
353  { 1.8448e-01, -1.0892e-01, 2.0336e-04, -3.5803e+01, 6.2453e-02 },
354  { -4.3394e+00, 6.7933e+00, 1.2905e-01, -6.7132e+00, -2.7930e-01 },
355  { 7.6150e-01, -1.3613e-01, 9.9347e-01, 1.1661e+00, 5.4822e-03 }
356  };
357 
358  //proton lateral parameters
359 
360  const double pro_rpar[4*Nrpar][5] = {
361  { 2.0473e+01, 1.8293e+01, -2.7258e-03, -1.4304e+03, -2.7740e-01 }, //Ecal for showerType = 1
362  { 4.3030e+00, 3.1796e+01, -1.4979e-01, 1.2390e+00, 2.8661e+00 },
363  { -7.1506e-01, 3.4815e+00, -1.8282e-01, 2.4101e+00, 5.4402e-01 },
364  { 2.5011e+00, 5.3549e+00, 1.7772e-01, 2.9511e+00, -8.8951e-01 },
365 
366  { 2.0252e+01, 1.6754e+00, -3.3596e+00, 2.1412e+00, -2.1051e+00 }, //Hcal for showerType = 1
367  { -1.6270e+00, 5.0493e+00, -1.6916e-01, 1.5961e+00, 7.2367e-01 },
368  { 5.0396e-01, 7.1958e-02, 1.4915e+00, 1.4498e+00, -4.0554e-02 },
369  { -1.3402e-01, 1.3657e-01, 2.3093e+00, -3.1556e+00, 5.2742e-03 },
370 
371  { 1.2771e+01, 2.8591e+00, -7.5931e-01, 1.5374e+00, -1.2623e+00 }, //Hcal for showerType = 2
372  { -4.1051e+00, 3.1735e+00, -3.5820e-01, 3.0133e+00, 1.0538e+00 },
373  { 4.9758e-01, 5.9756e-02, -1.9593e+00, 2.9959e+00, -1.3451e-02 },
374  { 2.0201e-01, 1.0967e-01, 1.2310e+00, 3.3239e+00, -6.2154e-02 },
375 
376  { -4.8143e+00, 1.0479e+01, -2.1974e-01, 3.0933e+00, 2.0855e+00 }, //Hcal for showerType = 3
377  { 2.5276e+01, 2.1056e+01, -4.4710e-01, -1.7116e+00, -6.7854e-01 },
378  { -5.8886e-01, 2.6766e+00, -1.9974e-01, 2.6381e+00, 4.3146e-01 },
379  { -1.7710e-02, 9.4732e-02, 4.2625e+00, 3.3545e+00, -2.9552e-02 }
380  };
381 
382  //Anti-Proton paramters
383 
384  const double pbar_emscale[2][5] = {
385  { 1.9246e+01, 1.8870e+01, -3.3127e-01, -4.8509e+00, -2.3785e-02 },
386  { 3.5744e+00, 3.3661e+00, -4.9667e-01, -4.1388e+00, -8.8824e-03 }
387  };
388 
389  const double pbar_hadscale[8][5] = {
390  //based on LogNormal of Had Energy
391  { -1.8857e+00, 9.4549e-01, 1.0173e+00, 2.3654e+00, 1.0381e+00 },
392  { 1.4276e+00, 8.3835e-01, -9.1466e-01, 2.3945e+00, 0.0000e+00 },
393  { 1.0005e+00, 5.0454e-01, -9.6777e-01, 2.4216e+00, -3.0214e-02 },
394  { -7.1494e-01, 3.7825e-01, 1.1627e+00, 2.7411e+00, 5.5459e-03 },
395  //LogNormal of Had Energy for mip-like response
396  { 1.5770e+00, 1.1215e+00, 4.2471e-01, 3.8321e+00, 5.2852e-01 },
397  { -6.5624e-02, 2.7895e-01, -3.2639e-01, 2.3605e+00, 6.3871e-02 },
398  //correction of response if ssp is between Ecal and Hcal
399  { 1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00 },
400  { 9.1334e-02, 8.4321e-02, 9.1042e-01, 2.7686e+00, -3.0344e-02 }
401  };
402 
403  const double pbar_correl_hadem[5] =
404  { -5.8253e-01, 3.6344e-01, -4.6371e-01, 2.8284e-01, 0.0000e+00 };
405 
406  //anti-proton longitudinal parameters
407 
408  const double pbar_par[8*NPar][5] = {
409  { 2.2962e+00, 1.8237e+00, 1.5856e-01, 1.1724e+01, -2.5137e-02 }, //Mean Ecal ShowerType1
410  { 5.1688e-01, 2.9636e-02, 3.9787e+01, 3.1512e+00, 1.7428e-01 },
411  { -1.3058e+00, 4.9530e-02, 1.5425e+00, 4.1539e+00, 5.7473e-02 },
412  { 2.9921e-01, 2.0357e-01, 2.9853e+00, 1.8306e+00, 1.5198e-01 },
413  { 3.8358e+00, 5.0343e+00, 1.3900e-01, 3.2366e+00, -6.1138e-01 },
414 
415  { 9.2918e-01, 3.4905e-01, 4.4705e-01, 4.3199e+00, -1.0055e-01 }, //Mean Hcal ShowerType1
416  { 1.1081e+00, -3.4156e-01, 1.2140e+00, 1.0438e+00, 2.1996e-02 },
417  { 2.5596e+01, -4.7613e+01, 3.8236e-02, -1.6258e+01, 8.3889e-01 },
418  { 1.4731e+00, -8.7282e-01, 7.6302e-01, 1.4135e+00, 4.0020e-02 },
419  { 1.5854e+00, -1.5158e+00, 5.3554e-01, 1.7953e+00, -2.1115e-02 },
420 
421  { 1.6397e+00, 1.4348e+00, 1.8597e-01, 5.0391e+00, -2.1112e-01 }, //Mean Hcal ShowerType2
422  { 1.9561e+00, 9.1382e-01, 5.2437e-01, 3.2069e+00, -3.6565e-01 },
423  { 8.4435e-01, 2.0674e+00, -4.2496e-01, -3.2228e+00, -2.0347e-01 },
424  { 1.3139e+00, 6.8898e-01, -9.0400e-01, 1.1612e+00, 5.2261e-02 },
425  { 1.3894e+00, -2.3487e+00, 3.5770e-01, 1.2953e+00, 1.4271e-01 },
426 
427  { 6.4737e-01, 7.3992e-02, 9.2454e-01, 4.8070e+00, -3.3204e-03 }, //Mean Hcal ShowerType3
428  { 9.6156e-01, 1.5021e-01, 2.1087e+00, 2.7135e+00, -5.6608e-02 },
429  { -2.9797e+00, 3.7631e+00, -1.8261e-01, 2.8053e+00, 4.3358e-01 },
430  { 1.0332e+00, 4.4820e-01, -1.2826e+00, 9.7747e-01, 9.4973e-02 },
431  { 1.6588e+00, 3.6401e+00, -2.4077e-01, 3.4642e-01, 2.7807e-01 },
432 
433  { 7.9348e-01, 6.0791e-01, -5.2480e-01, -4.1979e+00, -9.7594e-03 }, //Sigma Ecal ShowerType1
434  { 6.7623e-01, 1.1835e-01, 2.5074e+01, 2.9602e+00, -7.6399e-02 },
435  { 4.3953e-01, 5.0370e-02, 4.3043e+01, 2.9560e+00, -3.5965e-02 },
436  { 6.7091e-01, 1.1626e-01, -5.7986e-01, 2.9029e-01, -2.8586e-02 },
437  { 4.2821e-01, 8.0538e-02, -9.3285e-01, 1.2197e+00, -1.3245e-02 },
438 
439  { 2.3715e-01, 6.1849e-02, -7.9327e-01, 2.4567e+00, -1.5569e-02 }, //Sigma Hcal ShowerType1
440  { 7.7529e-01, 8.4253e-02, -2.8453e+00, 6.6960e-01, -2.3133e-02 },
441  { 9.7261e-01, 1.1777e-01, -1.2671e+00, 1.3352e+00, -9.0527e-02 },
442  { 1.1170e+00, 5.7366e-01, -7.9215e-01, 1.8696e+00, -3.2848e-02 },
443  { 1.0912e+00, 5.3449e-01, -8.3799e-01, 2.2648e+00, -5.3589e-02 },
444 
445  { 2.0364e-01, 1.1109e-01, -5.4911e-01, 2.2371e+00, -2.7723e-03 }, //Sigma Hcal ShowerType2
446  { 5.0058e-01, 8.3324e-02, -1.5426e+00, 4.3092e+00, 1.7710e-02 },
447  { 4.8258e+00, -5.3363e+00, -6.9399e-02, 1.3639e+01, -3.1232e-01 },
448  { 1.0654e+00, 4.3313e-01, -8.1864e-01, 1.5794e+00, -5.0301e-02 },
449  { 9.7937e-01, 4.9011e-01, -7.5635e-01, 2.0913e+00, -4.3251e-02 },
450 
451  { 2.4410e-01, 6.9464e-02, -5.7121e-01, 1.5400e+00, -1.8889e-02 }, //Sigma Hcal ShowerType3
452  { 6.8356e-01, 4.4101e-02, 2.0167e+01, 2.9570e+00, -4.3937e-02 },
453  { 2.9627e-01, -4.3021e-01, 1.3628e-01, 1.3540e+01, -5.0582e-02 },
454  { 9.2673e-01, 3.5636e-01, -6.8759e-01, 1.1528e+00, -3.6487e-02 },
455  { 8.1179e-01, 5.0238e-01, -5.7945e-01, 1.8664e+00, -9.9021e-03 }
456  };
457 
458  const double pbar_rho[8*NPar][5] = {
459  { 1.6115e-01, 7.7749e-02, 1.3576e+01, 6.4878e-01, -3.2481e-03 }, //Ecal ShowerType1
460  { 4.1746e-01, 6.4980e-02, -1.8986e+01, 2.9812e+00, 8.2744e-03 },
461  { 8.3730e-01, 5.2519e-02, 1.9080e+01, 2.9734e+00, -7.2551e-04 },
462  { 4.7186e-01, 2.3447e+00, 1.5313e-01, 2.1195e+00, -2.7598e-01 },
463  { 2.6387e-01, 1.2494e-01, 2.0479e+01, 2.9541e+00, 3.9418e-03 },
464  { 1.4487e-01, 1.7347e-01, 2.0829e+01, 2.9708e+00, -1.7575e-02 },
465  { -3.2116e-02, 1.0593e-01, -3.6256e+00, 2.0405e+00, 3.9701e-02 },
466  { 3.2368e-01, 3.4897e-01, 4.0152e-01, 1.9683e+00, -4.3362e-02 },
467  { 3.3095e-01, 8.8705e-02, 3.0619e+01, 3.2302e+00, -2.6705e-02 },
468  { 6.8646e-01, 1.3220e-01, 1.8172e+00, 1.8252e+00, 2.5609e-02 },
469 
470  { 2.5579e-01, 1.3158e-01, 2.0117e+00, 2.6834e+00, -1.8803e-02 }, //Hcal ShowerType1
471  { 5.4665e-01, 3.5282e-01, 7.0197e-01, 1.8240e+00, -1.2513e-01 },
472  { 7.8407e-01, 4.6267e-02, -1.9313e+00, 1.2055e+00, 3.1770e-03 },
473  { 1.9774e-01, 1.4336e+00, 1.0833e-01, 2.5249e+00, -6.1467e-02 },
474  { 4.8669e-01, 2.5554e-01, 7.8888e-01, 2.7417e+00, 0.0000e+00 },
475  { 2.5270e-01, 1.3412e-01, 7.7256e-01, 1.3281e+00, 1.2302e-02 },
476  { 4.3116e-02, 1.3977e-01, 1.2806e+00, 8.8424e-01, -2.2153e-02 },
477  { -1.1513e-01, 1.1712e+00, -1.9271e-01, 1.8406e+00, 2.3827e-01 },
478  { 4.8509e-01, 1.9358e-01, 6.8277e-01, 7.8174e-01, -1.2847e-02 },
479  { 6.8735e-01, 1.7148e-01, -8.8434e-01, 2.6390e+00, 1.8336e-02 },
480 
481  { 2.6738e-01, 1.6698e-01, 1.6242e+00, 2.3194e+00, -5.4439e-02 }, //Hcal ShowerType2
482  { 5.1565e-01, 3.1755e-01, 7.6866e-01, 1.4314e+00, -1.2355e-01 },
483  { 1.4488e+00, -1.0738e+00, 1.5373e-01, -4.4945e+00, 5.0825e-02 },
484  { -5.0531e-02, 5.4899e-01, 6.4605e-04, 6.2060e+00, 5.0480e-02 },
485  { 5.0045e-01, 2.8609e-01, 7.7954e-01, 2.4851e+00, -1.7178e-02 },
486  { 2.4183e-01, 1.1418e-01, 1.1548e+00, 6.8739e-01, 1.4071e-03 },
487  { -2.7267e-01, 2.9060e-01, -9.2639e-01, 3.3275e+00, 1.3254e-01 },
488  { 3.6510e-01, 6.8755e-02, 1.7529e+00, 3.3932e+00, 3.5968e-03 },
489  { -2.2256e+00, 4.0953e+00, 1.2237e-01, -6.3486e+00, -1.6487e-01 },
490  { 6.7014e-01, -1.8220e-01, 8.3413e-01, 2.3705e+00, 3.1919e-02 },
491 
492  { -9.2169e-02, 2.7435e-01, -1.4954e+00, 3.4936e+00, 5.4085e-02 }, //Hcal ShowerType3
493  { -5.7426e-02, 4.2696e-01, -9.1241e-01, 3.5629e+00, 1.0904e-01 },
494  { 7.4201e-01, 4.7628e-02, 1.5747e+00, -4.6694e+00, -1.4296e-02 },
495  { 1.8949e-02, 1.0836e-01, -1.2242e+00, 2.7149e+00, 1.3031e-02 },
496  { -1.5217e+00, 2.4872e+00, 8.1875e-02, -1.0666e+01, -1.4661e-03 },
497  { 2.2412e-01, 9.4590e-02, -2.4942e+00, 1.8584e+00, 1.6461e-02 },
498  { -5.4793e-01, 8.2034e-01, -4.2542e-01, 3.6711e+00, 1.9032e-01 },
499  { -2.5079e-01, -8.8412e-01, 1.5082e-01, 5.3762e+00, 1.3426e-01 },
500  { -7.3299e+00, 8.1099e+00, 4.0558e-01, -5.0694e+00, -5.9030e-02 },
501  { 5.6633e-01, -4.7278e-01, 4.1026e-01, 1.7565e+00, 1.1369e-01 }
502  };
503 
504  //anti-proton lateral parameters
505  const double pbar_rpar[4*Nrpar][5] = {
506  { 1.8776e+00, 2.8426e-01, 2.3530e+01, 3.0047e+00, -1.7332e-01 }, //Ecal for showerType = 1
507  { 9.1705e+00, 4.0028e+00, -4.4536e-01, 6.4847e-01, -5.9678e-01 },
508  { 5.9148e-01, 7.5552e-02, -3.0172e+01, 1.3417e+00, 2.0513e-02 },
509  { -1.3812e-01, 8.0542e-02, 5.5994e+01, 1.2251e+00, -1.6587e-02 },
510 
511  { 1.7618e+01, 4.3991e+00, -9.7543e-01, 2.7789e+00, -1.2115e+00 }, //Hcal for showerType = 1
512  { -7.6558e-02, 1.3021e+00, -1.4519e+00, 1.3966e+00, 1.6778e-01 },
513  { 5.9899e-01, 1.6407e-01, 8.4000e-01, 2.1795e+00, -7.0504e-02 },
514  { -1.4042e-01, 1.3416e-01, 2.3339e+00, -2.7260e+00, 5.8295e-03 },
515 
516  { 1.2340e+01, 5.3893e+00, -4.9219e-01, 1.5722e+00, -6.8433e-01 }, //Hcal for showerType = 2
517  { 1.0814e-01, 1.6610e+00, -1.1174e+00, -5.0663e+00, 1.9605e-01 },
518  { 4.9994e-01, 0.0000e+00, 0.0000e+00, 0.0000e+00, -2.0081e-02 },
519  { 4.6168e-01, 3.5308e-01, 5.4512e-01, 3.2036e+00, -1.4086e-01 },
520 
521  { -1.6306e+00, 3.5609e+00, -1.8903e-01, 6.6580e+00, 4.6936e-01 }, //Hcal for showerType = 3
522  { 6.5560e+00, 2.3432e+00, -9.8700e-01, 1.3115e+00, -6.5919e-01 },
523  { 7.4461e-01, 6.6885e-02, -8.9140e-01, -1.2050e+01, -4.5248e-02 },
524  { -2.4275e-01, 1.6289e-01, -2.9514e-04, -1.1208e+00, 4.3240e-02 }
525  };
526 
527  //Kaon-Minus paramters
528 
529  const double kminus_emscale[2][5] = {
530  { 8.3935e-01, 2.4180e-01, -1.4484e+00, 1.0397e+00, -2.2546e-02 },
531  { 2.3334e-01, 3.6810e-02, 1.0247e+01, 1.8556e+00, -1.9288e-02 }
532  };
533 
534  const double kminus_hadscale[8][5] = {
535  //based on LogNormal of Had Energy
536  { -5.9393e+01, 6.4013e+01, 1.0843e-01, -1.2441e+01, 4.2589e-01 },
537  { 1.3649e+01, 1.2266e+01, -4.9381e-01, -3.6921e+00, -1.1441e-01 },
538  { 1.2102e+00, 4.7458e-01, -9.4128e-01, 2.1482e+00, -5.1031e-02 },
539  { -6.2018e-01, 2.0288e-01, 1.1271e+00, 2.8535e+00, 0.0000e+00 },
540  //LogNormal of Had Energy for mip-like response
541  { 1.7106e-01, 1.3766e-01, 7.3972e-01, 4.1351e+00, 9.0284e-01 },
542  { 5.7520e-01, -1.2499e+00, 8.0351e-02, -4.7481e+00, 6.0937e-02 },
543  //correction of response if ssp is between Ecal and Hcal
544  { 1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00 },
545  { 4.2896e-02, 6.7448e-02, 7.4893e-01, 2.8837e+00, -1.1509e-02 }
546  };
547 
548  const double kminus_correl_hadem[5] =
549  { -7.3480e-01, 1.8541e-01, -7.6211e-01, 9.4878e-01, 2.5519e-03 };
550 
551  //kaon-minus longitudinal parameters
552 
553  const double kminus_par[8*NPar][5] = {
554  { 4.5986e-01, 8.0360e-02, 1.9103e+00, 2.7576e+00, 5.8024e-02 }, //Mean Ecal ShowerType1
555  { 5.6306e-01, 8.6175e-02, 7.5699e+00, 1.9142e+00, 1.4928e-01 },
556  { -1.4029e+00, 8.6442e-02, 4.7447e+00, 2.0026e+00, 7.8132e-02 },
557  { 5.1022e-01, 3.1763e-01, 1.4969e+00, 2.4374e+00, 1.0470e-01 },
558  { 4.4840e+00, 4.0858e+00, 2.4837e-01, 3.0941e+00, -8.4496e-01 },
559 
560  { 2.0880e+00, 1.9024e+00, 2.2650e-01, 4.5397e+00, -3.5337e-01 }, //Mean Hcal ShowerType1
561  { 1.2654e+00, -3.2827e-01, 1.4255e+00, 9.2278e-01, -2.8213e-02 },
562  { -8.5073e-01, 4.8949e-01, -1.1802e+00, 8.3284e-01, -1.6258e-01 },
563  { 2.4486e+00, -2.0627e+00, 4.1771e-01, 2.8534e-02, 5.9279e-02 },
564  { 1.0307e+01, -1.2125e+01, 1.8949e-01, -3.5864e+00, 1.9567e-01 },
565 
566  { 3.0833e+00, 4.5830e+00, 1.6380e-01, 3.8496e+00, -6.6617e-01 }, //Mean Hcal ShowerType2
567  { 1.5334e+01, -1.4598e+01, 4.5394e-01, -4.3694e+00, -3.0904e-02 },
568  { 1.1765e+01, 1.3252e+01, -6.1629e-01, -2.8795e+00, -1.5726e-01 },
569  { 1.1343e+01, 1.1615e+01, -2.4147e-01, -4.5126e+00, 1.6850e-01 },
570  { 1.5856e+01, -1.6981e+01, 2.0512e-01, -5.0463e+00, 1.2586e-01 },
571 
572  { 4.3696e-01, 5.9373e-02, 1.5520e+00, 2.8488e+00, 4.4038e-02 }, //Mean Hcal ShowerType3
573  { 5.4176e-01, 4.6653e+00, -1.1337e-01, 6.3126e-01, 4.5058e-01 },
574  { -1.3576e+00, 1.3274e-01, -2.0240e+00, 7.4848e-01, -1.3586e-01 },
575  { 8.4917e+00, 1.0264e+01, -1.5633e-01, -5.1903e+00, 3.7948e-01 },
576  { 2.4254e+01, -2.5119e+01, 2.1971e-01, -6.2339e+00, 1.1753e-01 },
577 
578  { 6.6849e-02, 1.1461e-01, -6.0693e-01, 3.5849e+00, 2.1510e-02 }, //Sigma Ecal ShowerType1
579  { 4.4751e-01, 6.7414e-02, -1.7351e+00, 3.1800e+00, -1.4204e-02 },
580  { 2.3277e-01, 9.8312e-02, -8.1359e-01, 3.3423e+00, 2.0069e-02 },
581  { 6.0942e-01, 2.5435e-01, -5.7408e-01, 2.2588e+00, -2.7305e-03 },
582  { 4.7227e-01, 6.9822e-01, -1.8614e-01, 5.0224e-01, 5.3256e-02 },
583 
584  { 2.7710e-01, 5.4557e-02, -9.1426e-01, 1.3652e+00, -2.1545e-02 }, //Sigma Hcal ShowerType1
585  { 8.5475e-01, 9.0583e-02, 3.6811e+00, 2.1694e+00, -8.1316e-02 },
586  { 9.9153e-01, 1.7054e-01, -2.3509e+00, 1.0596e+00, -7.7110e-02 },
587  { 1.2402e+00, 6.9760e-01, -5.8527e-01, 1.3975e+00, -4.4764e-02 },
588  { 1.1041e+00, 9.4012e-01, -4.0793e-01, 1.7098e+00, 6.9772e-03 },
589 
590  { 4.5757e-01, 2.2138e-01, -8.3016e-01, -6.5474e-01, -2.5193e-02 }, //Sigma Hcal ShowerType2
591  { 4.6612e-01, 1.1899e-01, -8.7496e-01, 4.0150e+00, 2.3778e-02 },
592  { 3.2115e+01, 3.1746e+01, -2.1227e-01, -1.1088e+01, 0.0000e+00 },
593  { 1.4220e+00, 1.3834e+00, -3.3056e-01, 2.2772e-01, 3.5460e-02 },
594  { 1.0770e+00, 2.3134e+00, -2.0993e-01, 7.6896e-01, 1.7177e-01 },
595 
596  { 1.0333e+00, 8.4267e-01, -2.2255e-01, -6.2499e+00, -2.0389e-02 }, //Sigma Hcal ShowerType3
597  { 4.5615e-01, 1.4719e-01, -1.0516e+00, 3.2503e+00, 2.9697e-02 },
598  { 7.1504e-01, 4.0454e-01, 6.8555e-04, -1.7412e+01, -5.0839e-02 },
599  { 3.4149e+00, 3.8890e+00, -1.8545e-01, -3.4381e+00, 9.7396e-02 },
600  { 3.7946e+00, 4.0784e+00, -1.7896e-01, -3.9979e+00, 5.1072e-02 }
601  };
602 
603  const double kminus_rho[8*NPar][5] = {
604  { -4.5951e-02, 3.6373e-01, 4.9181e-01, 6.9207e-01, 1.4429e-02 }, //Ecal ShowerType1
605  { 1.6996e-01, 4.2946e-01, 6.2046e-01, -3.7396e-02, -1.4572e-02 },
606  { 8.3137e-01, 2.9433e-02, 3.2221e+00, 1.6200e+00, -8.8843e-04 },
607  { -1.3483e-01, 9.3124e-02, -4.7555e+01, 1.9746e+00, 4.2961e-02 },
608  { 1.4051e-01, 9.3445e-02, -3.3320e+00, 2.3752e+00, 3.7645e-02 },
609  { 7.8256e-02, 1.2162e-01, -2.9812e+00, 2.1452e+00, 1.7049e-02 },
610  { -8.5479e-03, 9.6577e-02, -2.5794e+00, 2.1738e+00, 2.0528e-02 },
611  { 1.9550e-01, 9.8667e-02, 2.3510e+00, 9.3278e-01, -1.0451e-02 },
612  { 1.3706e-01, -1.1500e-01, 2.4339e+00, 2.3990e+00, 2.4792e-02 },
613  { 7.8427e-01, 1.6946e-01, 6.7162e-01, 1.6939e+00, -6.8641e-04 },
614 
615  { 3.0252e-01, 2.0969e-01, 1.2233e+00, 2.5744e+00, -3.5561e-02 }, //Hcal ShowerType1
616  { 5.0366e-01, 2.6734e-01, 9.1620e-01, 2.1240e+00, -7.5807e-02 },
617  { 7.6809e-01, 5.2611e-02, -3.3231e+00, 1.0126e+00, 6.1410e-03 },
618  { -1.7444e-01, 4.5874e-01, 1.0344e-03, 2.5811e+00, 7.8505e-02 },
619  { 2.5192e-01, 6.2636e-02, 1.3799e+00, 1.4204e+00, 6.1673e-02 },
620  { 2.9799e-01, 1.5329e-01, 1.2920e+00, 1.3530e+00, -1.0006e-02 },
621  { -1.8200e-01, 4.7254e-01, 5.0000e-01, -3.6152e-01, -4.6662e-02 },
622  { 1.6888e-01, -6.7814e-02, 6.5507e+00, 2.1306e+00, 4.1608e-02 },
623  { 4.8725e-01, 1.3573e-01, 1.7021e+00, 8.4936e-01, -1.0376e-02 },
624  { 8.3064e-01, 4.7621e-02, -1.4663e+00, 2.5843e+00, -4.1124e-02 },
625 
626  { 2.5820e-01, 2.1436e-01, 1.3064e+00, 2.1257e+00, -7.9308e-02 }, //Hcal ShowerType2
627  { 4.9118e-01, 2.3914e-01, 2.0019e+00, 1.8543e+00, -8.3488e-02 },
628  { 9.7119e-01, -2.4913e-01, 2.1162e+00, -2.2685e-01, -1.1841e-02 },
629  { -3.1042e-02, 8.5169e-02, 1.6265e+00, 1.3944e+00, 1.6725e-02 },
630  { 2.8430e-01, 1.0624e-01, 1.2654e+00, 8.8490e-01, 2.1335e-02 },
631  { -1.7883e+00, 2.5865e+00, -2.6200e-01, 3.5317e+00, 5.8945e-01 },
632  { -3.4648e-01, 3.0452e-01, -8.5428e-01, 3.7070e+00, 1.4478e-01 },
633  { 9.6655e-01, 1.1316e+00, -1.9048e-01, -3.5009e+00, 8.1655e-02 },
634  { 5.3247e-01, 1.1047e-01, 2.1538e+00, 7.0882e-01, -2.0301e-02 },
635  { 7.7207e-01, -8.5138e-02, 8.8082e-01, 2.0751e+00, -1.7783e-02 },
636 
637  { 4.4074e-02, 1.1797e-01, 9.8745e-01, 1.2754e+00, 0.0000e+00 }, //Hcal ShowerType3
638  { 3.0829e-01, 3.7222e-01, 6.4042e-01, 7.1567e-01, -4.8572e-02 },
639  { 7.2438e-01, 1.3993e-02, 2.8178e+01, 1.1987e+00, -1.8113e-03 },
640  { 3.3772e-02, 5.2144e-02, -2.1020e+01, 2.2369e+00, 0.0000e+00 },
641  { -8.7740e+01, 8.8231e+01, 5.0000e-01, -6.4849e+00, -2.5865e-02 },
642  { 1.6919e-01, 1.5093e-01, -1.6970e+00, 3.1185e+00, 0.0000e+00 },
643  { -2.7414e-01, 4.2637e-01, -4.5968e-01, 4.0680e+00, 1.1120e-01 },
644  { 3.5779e-01, 1.3748e+00, 1.5201e-01, 4.5045e-01, -1.4935e-01 },
645  { -4.4323e+01, 4.5858e+01, 1.4803e-01, -1.5183e+01, -1.7599e-01 },
646  { 4.4308e-01, -1.9361e+00, 1.7101e-01, 1.2267e+00, 2.6132e-01 }
647  };
648 
649  //kaon-minus lateral parameters
650 
651  const double kminus_rpar[4*Nrpar][5] = {
652  { 1.5255e+00, 7.1460e-01, -9.6109e-01, 2.1485e+00, 0.0000e+00 }, //Ecal for showerType = 1
653  { 8.1702e+00, 1.1618e+00, -3.2615e+00, 1.9965e+00, -9.8151e-01 },
654  { 8.0394e-01, 1.1379e-01, -5.7818e+01, 1.8315e+00, -3.1074e-02 },
655  { -3.4073e-01, 1.5998e-01, 5.6939e+01, 1.8305e+00, 1.5185e-02 },
656 
657  { 1.8546e+01, 3.9927e+00, -6.4386e-01, 2.5076e+00, -1.3842e+00 }, //Hcal for showerType = 1
658  { -2.9410e-01, 6.9922e-01, -7.1734e+01, 1.3073e+00, 8.2441e-02 },
659  { 4.8144e-01, 6.6668e-02, 1.4955e+00, 1.5032e+00, -2.8181e-02 },
660  { -5.4174e-03, 1.7362e-02, 1.2282e+01, 2.4035e-01, 6.1898e-04 },
661 
662  { 2.0741e+01, 3.9654e+01, -1.8029e-01, -4.2558e-01, 2.5375e+00 }, //Hcal for showerType = 2
663  { -1.1271e+00, 5.6222e-01, -5.1185e+01, 9.9351e-01, 1.9139e-01 },
664  { 4.8514e-01, 2.6816e-02, -3.7361e+01, 2.7567e+00, -8.2626e-03 },
665  { -9.0682e-02, 7.2902e-02, 2.3265e+00, -6.5459e+00, 8.8080e-03 },
666 
667  { -2.9248e+00, 5.6213e+00, -3.2437e-01, 3.0222e+00, 1.5066e+00 }, //Hcal for showerType = 3
668  { 7.4285e+00, 3.1786e+00, -9.1437e-01, 1.1013e+00, -7.0831e-01 },
669  { 9.6057e-01, 3.6162e-01, -3.6673e-01, -8.9691e-01, -3.2462e-02 },
670  { -2.9505e+01, 2.9362e+01, 1.0665e+00, -2.3946e+00, 1.6595e-02 }
671  };
672 
673  //Kaon-Plus paramters
674 
675  const double kplus_emscale[2][5] = {
676  { 7.8775e-01, 2.7022e-01, -1.4291e+00, 1.3980e+00, -1.2795e-02 },
677  { 2.2224e-01, 1.7480e-02, 1.1984e+01, 1.8265e+00, -1.6718e-02 }
678  };
679 
680  const double kplus_hadscale[8][5] = {
681  //based on LogNormal of Had Energy
682  { -1.1629e+02, 1.1506e+02, 4.8753e-01, -4.4057e+00, 1.0940e+00 },
683  { 2.2018e+00, 1.1327e+00, -6.4208e-01, 1.2240e-01, -7.2109e-02 },
684  { 1.5217e+00, 1.9949e+00, -4.0047e-01, 1.0897e+00, 1.4482e-01 },
685  { -1.2902e+01, 1.5504e+01, 1.1257e-01, -8.4732e+00, -3.0451e-01 },
686  //LogNormal of Had Energy for mip-like response
687  { 3.4887e-01, 1.9201e-01, 8.5205e-01, 3.8457e+00, 8.5097e-01 },
688  { 5.1252e-01, 7.7619e-01, -1.3001e-01, -4.1296e+00, 3.8339e-02 },
689  //correction of response if ssp is between Ecal and Hcal
690  { 1.8736e-01, 1.7391e-01, 1.0991e+00, 3.3557e+00, 0.0000e+00 },
691  { 8.2313e-01, 1.3027e+00, 2.5263e-01, 2.8203e+00, -2.8548e-01 }
692  };
693 
694  const double kplus_correl_hadem[5] =
695  { -6.9227e-01, 1.8762e-01, -7.1395e-01, 7.1711e-01, -7.3016e-03 };
696 
697  //kaon-plus longitudinal parameters
698 
699  const double kplus_par[8*NPar][5] = {
700  { 3.8772e-01, 7.5923e-02, 2.3857e+00, 2.8135e+00, 7.0675e-02 }, //Mean Ecal ShowerType1
701  { 1.4815e+00, 9.8518e-01, 3.4748e-01, 3.3609e+00, -1.2142e-01 },
702  { -7.4688e-01, 7.4632e-01, 3.7387e-01, 3.1767e+00, -1.2898e-01 },
703  { 6.4780e-01, 3.3263e-01, 1.6380e+00, 2.8163e+00, 7.5372e-02 },
704  { 4.8648e+00, 4.4851e+00, 2.5096e-01, 3.1485e+00, -9.4875e-01 },
705 
706  { 1.6983e+00, 1.3064e+00, 2.8039e-01, 4.3610e+00, -2.8089e-01 }, //Mean Hcal ShowerType1
707  { 1.4443e+00, -5.0550e-01, 1.0680e+00, 5.6710e-01, -1.5739e-02 },
708  { -5.3822e-01, 1.0584e+00, -6.4612e-01, 3.5104e-01, -1.1463e-01 },
709  { 2.3248e+00, -2.2851e+00, 3.6709e-01, 1.3197e-01, 1.2563e-01 },
710  { 2.1837e+00, -2.2043e+00, 4.6015e-01, 9.2583e-01, -7.4904e-03 },
711 
712  { 3.3203e+00, 5.0310e+00, 1.6001e-01, 3.9441e+00, -7.1274e-01 }, //Mean Hcal ShowerType2
713  { 1.1513e+00, -1.8085e-01, 3.3413e+00, 7.0993e-01, -6.3645e-02 },
714  { 1.3898e+01, 1.5638e+01, -4.1890e-01, -3.9674e+00, -1.1401e-01 },
715  { 2.0012e+01, 1.9866e+01, -2.9212e-01, -5.0006e+00, 1.1974e-01 },
716  { 2.4295e+00, -2.2834e+00, 4.8288e-01, 2.9977e-01, -3.0573e-02 },
717 
718  { 3.8408e-01, 6.1344e-02, 1.0667e+00, 3.0242e+00, 5.0278e-02 }, //Mean Hcal ShowerType3
719  { 1.8964e+00, 5.1610e+00, -8.8169e-02, -1.9272e+00, 3.1513e-01 },
720  { 1.4200e+01, 1.6127e+01, -2.8959e-01, -6.1679e+00, -7.4159e-02 },
721  { 2.7497e+01, 2.7062e+01, -3.3540e-01, -5.5799e+00, 1.0200e-01 },
722  { 2.3118e+01, -2.3752e+01, 2.6134e-01, -5.2177e+00, 9.7844e-02 },
723 
724  { -3.5273e-02, 2.6492e-01, -2.7154e-01, 4.1599e+00, 3.7807e-02 }, //Sigma Ecal ShowerType1
725  { 5.4462e-01, 4.9507e-02, -1.8513e+01, 2.9250e+00, -3.6791e-02 },
726  { 3.7860e-01, 1.8549e-02, -2.9438e+00, 2.4898e+00, -2.2790e-02 },
727  { 6.6597e-01, 6.5499e-02, -2.1924e+01, 2.9525e+00, -4.7980e-02 },
728  { 1.0753e+00, 9.2583e-01, -1.4431e-01, -5.0882e+00, 2.9234e-03 },
729 
730  { 2.6008e-01, 6.1267e-02, -1.7484e+00, 1.6088e+00, -1.7668e-02 }, //Sigma Hcal ShowerType1
731  { 8.6750e-01, 8.7283e-02, 2.5881e+00, 2.5618e+00, -8.5318e-02 },
732  { 1.0368e+00, 2.5737e-01, -1.3413e+00, 8.7206e-01, -7.0973e-02 },
733  { 1.2517e+00, 6.4973e-01, -7.3033e-01, 1.4095e+00, -5.3997e-02 },
734  { 1.9193e+00, 1.9482e+00, -3.0639e-01, 6.6016e-02, 2.9074e-02 },
735 
736  { 2.7361e-01, 5.8267e-02, -1.5294e+00, 1.0185e+00, -2.2741e-02 }, //Sigma Hcal ShowerType2
737  { 5.3653e-01, 1.0764e-01, -6.6807e-01, 4.7074e+00, -3.4054e-03 },
738  { 5.8520e+01, 5.8162e+01, -2.3280e-01, -1.1136e+01, 0.0000e+00 },
739  { 1.2192e+00, 6.9207e-01, -6.8090e-01, 8.9059e-01, -3.8146e-02 },
740  { 3.2972e+00, 2.9315e+00, -3.4433e-01, -1.6434e+00, -2.8242e-02 },
741 
742  { 1.4349e+00, 1.4443e+00, -1.2105e-01, -9.0252e+00, 2.3150e-03 }, //Sigma Hcal ShowerType3
743  { 4.9465e-01, 1.8568e-01, -9.3747e-01, 3.2665e+00, 2.3940e-02 },
744  { 2.3407e+00, -3.8601e+00, 2.9374e-02, -1.4549e+01, 1.4344e-02 },
745  { 4.8451e+00, 4.6687e+00, -2.7170e-01, -3.7684e+00, 1.8616e-02 },
746  { 5.0127e+00, 5.0025e+00, -2.3093e-01, -4.0513e+00, 2.3821e-02 }
747  };
748 
749  const double kplus_rho[8*NPar][5] = {
750  { 3.8772e-01, 7.5923e-02, 2.3857e+00, 2.8135e+00, 7.0675e-02 }, //Ecal ShowerType1
751  { 1.4815e+00, 9.8518e-01, 3.4748e-01, 3.3609e+00, -1.2142e-01 },
752  { -7.4688e-01, 7.4632e-01, 3.7387e-01, 3.1767e+00, -1.2898e-01 },
753  { 6.4780e-01, 3.3263e-01, 1.6380e+00, 2.8163e+00, 7.5372e-02 },
754  { 4.8648e+00, 4.4851e+00, 2.5096e-01, 3.1485e+00, -9.4875e-01 },
755  { -3.5273e-02, 2.6492e-01, -2.7154e-01, 4.1599e+00, 3.7807e-02 },
756  { 5.4462e-01, 4.9507e-02, -1.8513e+01, 2.9250e+00, -3.6791e-02 },
757  { 3.7860e-01, 1.8549e-02, -2.9438e+00, 2.4898e+00, -2.2790e-02 },
758  { 6.6597e-01, 6.5499e-02, -2.1924e+01, 2.9525e+00, -4.7980e-02 },
759  { 1.0753e+00, 9.2583e-01, -1.4431e-01, -5.0882e+00, 2.9234e-03 },
760 
761  { 1.6983e+00, 1.3064e+00, 2.8039e-01, 4.3610e+00, -2.8089e-01 }, //Hcal ShowerType1
762  { 1.4443e+00, -5.0550e-01, 1.0680e+00, 5.6710e-01, -1.5739e-02 },
763  { -5.3822e-01, 1.0584e+00, -6.4612e-01, 3.5104e-01, -1.1463e-01 },
764  { 2.3248e+00, -2.2851e+00, 3.6709e-01, 1.3197e-01, 1.2563e-01 },
765  { 2.1837e+00, -2.2043e+00, 4.6015e-01, 9.2583e-01, -7.4904e-03 },
766  { 2.6008e-01, 6.1267e-02, -1.7484e+00, 1.6088e+00, -1.7668e-02 },
767  { 8.6750e-01, 8.7283e-02, 2.5881e+00, 2.5618e+00, -8.5318e-02 },
768  { 1.0368e+00, 2.5737e-01, -1.3413e+00, 8.7206e-01, -7.0973e-02 },
769  { 1.2517e+00, 6.4973e-01, -7.3033e-01, 1.4095e+00, -5.3997e-02 },
770  { 1.9193e+00, 1.9482e+00, -3.0639e-01, 6.6016e-02, 2.9074e-02 },
771 
772  { 2.4423e-01, 2.2650e-01, 1.3877e+00, 1.8939e+00, -7.5389e-02 }, //Hcal ShowerType2
773  { 4.6760e-01, 2.6361e-01, 1.4374e+00, 1.6434e+00, -9.9844e-02 },
774  { 7.4796e-01, -4.4661e-02, 1.7292e+01, 7.0879e-01, 2.4763e-03 },
775  { -6.4806e+00, 6.4231e+00, 1.0400e+00, -2.2650e+00, 4.7778e-02 },
776  { 3.4440e-01, 1.1933e-01, 2.1685e+00, 1.3074e+00, 1.5225e-02 },
777  { -1.1990e+00, 1.8273e+00, -2.2560e-01, 4.3776e+00, 3.3393e-01 },
778  { -1.3851e+00, 1.7655e+00, -2.6978e-01, 3.6071e+00, 4.2187e-01 },
779  { 2.8248e-01, 4.8579e-02, 1.0095e+00, 3.5771e+00, 1.1249e-02 },
780  { -4.8075e+00, 5.5604e+00, 5.3060e-01, -3.1000e+00, -4.5298e-02 },
781  { 7.2139e-01, -1.2988e-01, 1.3460e+00, 1.7526e+00, 8.9762e-03 },
782 
783  { 4.7434e-02, 3.1872e+00, 1.6477e-01, 3.0236e-01, -4.2486e-01 }, //Hcal ShowerType3
784  { 9.4033e-02, 5.7792e-01, 6.8058e-01, 2.1194e-01, -8.2046e-02 },
785  { 7.4700e-01, 2.0561e-02, 4.2268e+01, 1.6399e+00, -1.1949e-02 },
786  { 5.2255e-02, 4.6034e-02, -1.7110e+00, 1.8247e+00, 0.0000e+00 },
787  { -1.0684e+02, 1.0739e+02, 5.0000e-01, -6.3340e+00, -2.6291e-02 },
788  { 2.0912e-01, 1.3043e-01, -4.9417e+00, 2.9813e+00, 0.0000e+00 },
789  { -1.8831e+00, 3.4253e+00, -1.8757e-01, 3.5963e+00, 5.6274e-01 },
790  { 1.4591e-03, 1.6858e+00, 1.2556e-01, -7.6405e-01, -1.1943e-01 },
791  { -6.2423e+00, 7.1383e+00, 3.8814e-01, -4.3872e+00, -9.4735e-02 },
792  { 7.2200e-01, -1.1568e-01, 1.3604e+00, 1.3338e+00, 1.6450e-02 }
793  };
794 
795  //kaon-plus lateral parameters
796 
797  const double kplus_rpar[4*Nrpar][5] = {
798  { 2.0917e+00, 1.2528e+00, -6.9494e-01, 1.7600e+00, 0.0000e+00 }, //Ecal for showerType = 1
799  { 7.8866e+00, 6.3777e-01, -4.3657e+00, 2.0492e+00, -1.0208e+00 },
800  { 7.3150e-01, 1.6104e-01, -8.4279e+00, 2.0670e+00, -3.3696e-03 },
801  { -2.4312e-01, 1.8371e-01, 1.5265e+01, 2.0458e+00, -1.6358e-02 },
802 
803  { 2.0008e+01, 2.1622e+00, -3.4949e+00, 1.9954e+00, -2.0569e+00 }, //Hcal for showerType = 1
804  { -4.1668e-02, 5.4366e-01, -1.3780e+01, 7.0957e-01, -1.1259e-02 },
805  { 5.0855e-01, 1.0979e-01, 1.2233e+00, 1.4132e+00, -4.5390e-02 },
806  { -5.9436e-02, -1.0396e-01, 3.7419e-01, 1.9661e+00, 3.6075e-02 },
807 
808  { 1.2970e+01, 2.5925e+00, -1.5404e+00, 1.6391e+00, -1.3978e+00 }, //Hcal for showerType = 2
809  { 1.9183e-01, 1.3373e+00, -1.1090e+00, -7.6909e+00, 6.3587e-02 },
810  { 4.3892e-01, 7.0247e-02, -1.8419e+01, 2.9160e+00, 1.2592e-02 },
811  { 6.2448e-02, 6.3318e-02, 6.1143e-01, 2.4479e+00, -2.4459e-02 },
812 
813  { 2.0954e+00, 4.0916e-01, -6.9750e+01, 1.9870e+00, -5.3104e-02 }, //Hcal for showerType = 3
814  { 6.4162e+00, 1.8264e+00, -9.6305e-01, 1.5012e+00, -7.5696e-01 },
815  { 4.3430e+00, 1.3126e+01, -2.5342e-02, -1.0939e+01, 2.1736e-01 },
816  { -5.0784e+01, 5.0870e+01, 2.0000e-01, -1.3669e+01, -1.2081e-02 }
817  };
818 
819  //HcalOuter parameters
820  const double ho_nonzero[5] =
821  {4.79943e-01,4.61158e-01,7.09011e-01,4.86440e+00,0.0000e+00};
822 }
823 
824 #endif
Gflash::pion_rpar
const double pion_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:214
Gflash::pro_hadscale
const double pro_hadscale[8][5]
Definition: GflashNameSpace.h:243
Gflash::Z
const double Z[kNumberCalorimeter]
Definition: GflashNameSpace.h:58
Gflash::ZFrontCrystalEE
const double ZFrontCrystalEE
Definition: GflashNameSpace.h:36
Gflash::MinDistanceToOut
const double MinDistanceToOut
Definition: GflashNameSpace.h:50
Gflash::RLTHAD
const double RLTHAD[kNumberCalorimeter]
Definition: GflashNameSpace.h:82
Gflash::kplus_hadscale
const double kplus_hadscale[8][5]
Definition: GflashNameSpace.h:680
Gflash::kminus_rho
const double kminus_rho[8 *NPar][5]
Definition: GflashNameSpace.h:603
Gflash::pbar_par
const double pbar_par[8 *NPar][5]
Definition: GflashNameSpace.h:408
Gflash::divisionStep
const double divisionStep
Definition: GflashNameSpace.h:61
Gflash::kminus_rpar
const double kminus_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:651
Gflash::Zmax
const double Zmax[kNumberCalorimeter]
Definition: GflashNameSpace.h:28
Gflash::pro_par
const double pro_par[8 *NPar][5]
Definition: GflashNameSpace.h:262
Gflash::zBackEE
double zBackEE(const Gflash3Vector &position)
Definition: GflashNameSpace.cc:89
Gflash::pion_hadscale
const double pion_hadscale[8][5]
Definition: GflashNameSpace.h:97
Gflash::kplus_correl_hadem
const double kplus_correl_hadem[5]
Definition: GflashNameSpace.h:694
Gflash::pbar_hadscale
const double pbar_hadscale[8][5]
Definition: GflashNameSpace.h:389
Gflash::pion_par
const double pion_par[8 *NPar][5]
Definition: GflashNameSpace.h:116
Gflash::Nrpar
const int Nrpar
Definition: GflashNameSpace.h:88
Gflash::EtaMax
const double EtaMax[kNumberCalorimeter]
Definition: GflashNameSpace.h:33
Gflash::ZOffCrystalEE
const double ZOffCrystalEE
Definition: GflashNameSpace.h:38
Gflash::kplus_rho
const double kplus_rho[8 *NPar][5]
Definition: GflashNameSpace.h:749
Gflash::rMoliere
const double rMoliere[kNumberCalorimeter]
Definition: GflashNameSpace.h:55
Gflash::ho_nonzero
const double ho_nonzero[5]
Definition: GflashNameSpace.h:820
Gflash::pro_rho
const double pro_rho[8 *NPar][5]
Definition: GflashNameSpace.h:312
Gflash::getCalorimeterNumber
CalorimeterNumber getCalorimeterNumber(const Gflash3Vector &position)
Definition: GflashNameSpace.cc:7
Gflash::kminus_hadscale
const double kminus_hadscale[8][5]
Definition: GflashNameSpace.h:534
Gflash::kminus_par
const double kminus_par[8 *NPar][5]
Definition: GflashNameSpace.h:553
Gflash::LengthCrystalEE
const double LengthCrystalEE
Definition: GflashNameSpace.h:40
Gflash::MinEnergyCutOffForHO
const double MinEnergyCutOffForHO
Definition: GflashNameSpace.h:70
Gflash::pro_correl_hadem
const double pro_correl_hadem[5]
Definition: GflashNameSpace.h:257
Gflash::energyCutOff
const double energyCutOff
Definition: GflashNameSpace.h:43
Gflash::kminus_emscale
const double kminus_emscale[2][5]
Definition: GflashNameSpace.h:529
Gflash::pion_correl_hadem
const double pion_correl_hadem[5]
Definition: GflashNameSpace.h:111
Gflash::EtaMin
const double EtaMin[kNumberCalorimeter]
Definition: GflashNameSpace.h:32
Gflash::LengthCrystalEB
const double LengthCrystalEB
Definition: GflashNameSpace.h:39
Gflash::criticalEnergy
const double criticalEnergy
Definition: GflashNameSpace.h:54
Gflash::kminus_correl_hadem
const double kminus_correl_hadem[5]
Definition: GflashNameSpace.h:548
Gflash::pro_emscale
const double pro_emscale[2][5]
Definition: GflashNameSpace.h:238
Gflash::ROffCrystalEB
const double ROffCrystalEB
Definition: GflashNameSpace.h:37
Gflash::pbar_correl_hadem
const double pbar_correl_hadem[5]
Definition: GflashNameSpace.h:403
Gflash::NPar
const int NPar
Definition: GflashNameSpace.h:87
Gflash::kplus_rpar
const double kplus_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:797
Gflash::intLength
const double intLength[kNumberCalorimeter]
Definition: GflashNameSpace.h:57
Gflash::rhoBackEB
double rhoBackEB(const Gflash3Vector &position)
Definition: GflashNameSpace.cc:82
Gflash::findShowerType
int findShowerType(const Gflash3Vector &position)
Definition: GflashNameSpace.cc:44
Gflash::RFrontCrystalEB
const double RFrontCrystalEB
Definition: GflashNameSpace.h:35
Gflash::maxShowerDepthforR50
const double maxShowerDepthforR50
Definition: GflashNameSpace.h:64
Gflash::Zmin
const double Zmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:27
Gflash::radLength
const double radLength[kNumberCalorimeter]
Definition: GflashNameSpace.h:56
Gflash::QuasiElasticLike
const double QuasiElasticLike
Definition: GflashNameSpace.h:47
Gflash::scaleSensitiveHE
const double scaleSensitiveHE
Definition: GflashNameSpace.h:75
Gflash3Vector
CLHEP::Hep3Vector Gflash3Vector
Definition: Gflash3Vector.h:6
Gflash::kplus_emscale
const double kplus_emscale[2][5]
Definition: GflashNameSpace.h:675
Gflash::pro_rpar
const double pro_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:360
Gflash::kplus_par
const double kplus_par[8 *NPar][5]
Definition: GflashNameSpace.h:699
Gflash::pbar_rho
const double pbar_rho[8 *NPar][5]
Definition: GflashNameSpace.h:458
Gflash::pion_rho
const double pion_rho[8 *NPar][5]
Definition: GflashNameSpace.h:166
Gflash::pbar_emscale
const double pbar_emscale[2][5]
Definition: GflashNameSpace.h:384
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
Gflash::maxLateralArmforR50
const double maxLateralArmforR50
Definition: GflashNameSpace.h:67
Gflash::Rmax
const double Rmax[kNumberCalorimeter]
Definition: GflashNameSpace.h:30
Gflash::pbar_rpar
const double pbar_rpar[4 *Nrpar][5]
Definition: GflashNameSpace.h:505
Gflash::scaleSensitiveHB
const double scaleSensitiveHB
Definition: GflashNameSpace.h:73
Gflash::Rmin
const double Rmin[kNumberCalorimeter]
Definition: GflashNameSpace.h:29
Gflash::pion_emscale
const double pion_emscale[2][5]
Definition: GflashNameSpace.h:92
Gflash::SAMHAD
const double SAMHAD[3][kNumberCalorimeter]
Definition: GflashNameSpace.h:79
Generated for CMSSW Reference Manual by doxygen 1.8.5