CMS 3D CMS Logo
EBShape.cc
Go to the documentation of this file.
1 #include <cmath>
2 
3 #include "SimCalorimetry/EcalSimAlgos/interface/EBShape.h"
4 
5 EBShape::~EBShape()
6 {
7 }
8 
9 EBShape::EBShape() :
10  EcalShapeBase( true )
11 {
12  buildMe() ;
13 }
14 
15 double
16 EBShape::threshold() const
17 {
18  return 0.00013 ; // magic number for barrel
19 }
20 
21 void
22 EBShape::fillShape( EcalShapeBase::DVec& aVec ) const
23 {
24 
25  for( unsigned int i ( k1NSecBins ) ; i != k1NSecBinsTotal ; ++i )
26  {
27  aVec[i] = exp(2.39735 - 0.0151053* ((double)i+1.0));
28 // shapeArray[i] = 10.82918*exp(-i/66.202 );
29  }
30 
31 
32  aVec[0] = 6.94068e-05 ;
33  aVec[1] = -5.03304e-05 ;
34  aVec[2] = -2.13404e-05 ;
35  aVec[3] = 6.017e-05 ;
36  aVec[4] = 2.01697e-05 ;
37  aVec[5] = 0.000114845 ;
38  aVec[6] = 2.13998e-05 ;
39  aVec[7] = 2.74476e-05 ;
40  aVec[8] = 5.2824e-05 ;
41  aVec[9] = 8.754e-05 ;
42  aVec[10] = 2.95346e-06 ;
43  aVec[11] = -7.58699e-05 ;
44  aVec[12] = -2.72224e-05 ;
45  aVec[13] = 3.10997e-06 ;
46  aVec[14] = -3.97771e-05 ;
47  aVec[15] = -1.06916e-05 ;
48  aVec[16] = -0.000113865 ;
49  aVec[17] = 6.05044e-05 ;
50  aVec[18] = -5.81202e-05 ;
51  aVec[19] = -6.58974e-06 ;
52  aVec[20] = 5.37494e-05 ;
53  aVec[21] = -0.000123729 ;
54  aVec[22] = 7.50938e-06 ;
55  aVec[23] = -1.35628e-05 ;
56  aVec[24] = 8.33725e-05 ;
57  aVec[25] = 3.19299e-05 ;
58  aVec[26] = -3.09232e-05 ;
59  aVec[27] = -7.0086e-05 ;
60  aVec[28] = 1.78937e-06 ;
61  aVec[29] = -2.20365e-05 ;
62  aVec[30] = 7.68054e-05 ;
63  aVec[31] = -2.5368e-05 ;
64  aVec[32] = 5.67291e-06 ;
65  aVec[33] = 5.87096e-05 ;
66  aVec[34] = -2.62771e-06 ;
67  aVec[35] = 4.31832e-05 ;
68  aVec[36] = 8.33616e-06 ;
69  aVec[37] = 7.27813e-05 ;
70  aVec[38] = 7.6159e-05 ;
71  aVec[39] = -1.60446e-05 ;
72  aVec[40] = -4.12127e-06 ;
73  aVec[41] = -5.93381e-05 ;
74  aVec[42] = 1.61444e-05 ;
75  aVec[43] = -5.49559e-05 ;
76  aVec[44] = 5.55254e-05 ;
77  aVec[45] = 3.32251e-05 ;
78  aVec[46] = -3.15897e-05 ;
79  aVec[47] = 7.86588e-05 ;
80  aVec[48] = -2.9704e-05 ;
81  aVec[49] = 5.66838e-05 ;
82  aVec[50] = 2.85281e-05 ;
83  aVec[51] = -3.02436e-05 ;
84  aVec[52] = -4.16265e-05 ;
85  aVec[53] = -1.63191e-05 ;
86  aVec[54] = 6.61193e-05 ;
87  aVec[55] = 9.23766e-05 ;
88  aVec[56] = 6.68903e-05 ;
89  aVec[57] = -3.20994e-05 ;
90  aVec[58] = 0.00011082 ;
91  aVec[59] = -4.07997e-05 ;
92  aVec[60] = -8.29046e-06 ;
93  aVec[61] = -7.42197e-05 ;
94  aVec[62] = -1.64386e-05 ;
95  aVec[63] = 1.02508e-05 ;
96  aVec[64] = 7.10995e-06 ;
97  aVec[65] = -5.87486e-05 ;
98  aVec[66] = -0.000101201 ;
99  aVec[67] = 1.62003e-05 ;
100  aVec[68] = -2.53093e-05 ;
101  aVec[69] = 2.65239e-05 ;
102  aVec[70] = -2.68722e-05 ;
103  aVec[71] = -4.02001e-05 ;
104  aVec[72] = 5.0674e-05 ;
105  aVec[73] = -1.75884e-05 ;
106  aVec[74] = 4.7902e-05 ;
107  aVec[75] = -1.01079e-05 ;
108  aVec[76] = 1.08427e-05 ;
109  aVec[77] = -0.000112906 ;
110  aVec[78] = 3.33076e-05 ;
111  aVec[79] = 0.000181201 ;
112  aVec[80] = 0.000426875 ;
113  aVec[81] = 0.00114222 ;
114  aVec[82] = 0.00237804 ;
115  aVec[83] = 0.00541858 ;
116  aVec[84] = 0.0089021 ;
117  aVec[85] = 0.0149157 ;
118  aVec[86] = 0.0231397 ;
119  aVec[87] = 0.0344671 ;
120  aVec[88] = 0.0471013 ;
121  aVec[89] = 0.0625517 ;
122  aVec[90] = 0.0857351 ;
123  aVec[91] = 0.108561 ;
124  aVec[92] = 0.133481 ;
125  aVec[93] = 0.163557 ;
126  aVec[94] = 0.200243 ;
127  aVec[95] = 0.225919 ;
128  aVec[96] = 0.269213 ;
129  aVec[97] = 0.302929 ;
130  aVec[98] = 0.342722 ;
131  aVec[99] = 0.378522 ;
132  aVec[100] = 0.436563 ;
133  aVec[101] = 0.467581 ;
134  aVec[102] = 0.510133 ;
135  aVec[103] = 0.550063 ;
136  aVec[104] = 0.583509 ;
137  aVec[105] = 0.619187 ;
138  aVec[106] = 0.653245 ;
139  aVec[107] = 0.686101 ;
140  aVec[108] = 0.721178 ;
141  aVec[109] = 0.745129 ;
142  aVec[110] = 0.774163 ;
143  aVec[111] = 0.799011 ;
144  aVec[112] = 0.822177 ;
145  aVec[113] = 0.838315 ;
146  aVec[114] = 0.858847 ;
147  aVec[115] = 0.875559 ;
148  aVec[116] = 0.891294 ;
149  aVec[117] = 0.90537 ;
150  aVec[118] = 0.919617 ;
151  aVec[119] = 0.930632 ;
152  aVec[120] = 0.936216 ;
153  aVec[121] = 0.947739 ;
154  aVec[122] = 0.955306 ;
155  aVec[123] = 0.961876 ;
156  aVec[124] = 0.968124 ;
157  aVec[125] = 0.97327 ;
158  aVec[126] = 0.977513 ;
159  aVec[127] = 0.984885 ;
160  aVec[128] = 0.986497 ;
161  aVec[129] = 0.990039 ;
162  aVec[130] = 0.994798 ;
163  aVec[131] = 0.994884 ;
164  aVec[132] = 0.99795 ;
165  aVec[133] = 0.99834 ;
166  aVec[134] = 0.999607 ;
167  aVec[135] = 1 ;
168  aVec[136] = 0.999047 ;
169  aVec[137] = 0.998745 ;
170  aVec[138] = 0.999219 ;
171  aVec[139] = 0.99814 ;
172  aVec[140] = 0.995082 ;
173  aVec[141] = 0.992449 ;
174  aVec[142] = 0.990418 ;
175  aVec[143] = 0.985032 ;
176  aVec[144] = 0.982308 ;
177  aVec[145] = 0.978696 ;
178  aVec[146] = 0.975656 ;
179  aVec[147] = 0.971027 ;
180  aVec[148] = 0.964811 ;
181  aVec[149] = 0.959428 ;
182  aVec[150] = 0.95096 ;
183  aVec[151] = 0.947428 ;
184  aVec[152] = 0.9419 ;
185  aVec[153] = 0.933223 ;
186  aVec[154] = 0.926482 ;
187  aVec[155] = 0.922172 ;
188  aVec[156] = 0.912777 ;
189  aVec[157] = 0.907388 ;
190  aVec[158] = 0.897289 ;
191  aVec[159] = 0.891889 ;
192  aVec[160] = 0.882056 ;
193  aVec[161] = 0.873382 ;
194  aVec[162] = 0.865442 ;
195  aVec[163] = 0.860032 ;
196  aVec[164] = 0.85202 ;
197  aVec[165] = 0.841013 ;
198  aVec[166] = 0.833802 ;
199  aVec[167] = 0.825259 ;
200  aVec[168] = 0.815013 ;
201  aVec[169] = 0.807465 ;
202  aVec[170] = 0.799428 ;
203  aVec[171] = 0.792165 ;
204  aVec[172] = 0.783088 ;
205  aVec[173] = 0.773392 ;
206  aVec[174] = 0.764982 ;
207  aVec[175] = 0.752174 ;
208  aVec[176] = 0.746487 ;
209  aVec[177] = 0.737678 ;
210  aVec[178] = 0.727396 ;
211  aVec[179] = 0.718692 ;
212  aVec[180] = 0.712737 ;
213  aVec[181] = 0.702738 ;
214  aVec[182] = 0.69559 ;
215  aVec[183] = 0.684389 ;
216  aVec[184] = 0.677989 ;
217  aVec[185] = 0.667643 ;
218  aVec[186] = 0.659009 ;
219  aVec[187] = 0.650217 ;
220  aVec[188] = 0.644479 ;
221  aVec[189] = 0.636017 ;
222  aVec[190] = 0.625257 ;
223  aVec[191] = 0.618507 ;
224  aVec[192] = 0.609798 ;
225  aVec[193] = 0.600097 ;
226  aVec[194] = 0.592788 ;
227  aVec[195] = 0.584895 ;
228  aVec[196] = 0.578228 ;
229  aVec[197] = 0.569299 ;
230  aVec[198] = 0.560576 ;
231  aVec[199] = 0.552404 ;
232  aVec[200] = 0.541405 ;
233  aVec[201] = 0.536271 ;
234  aVec[202] = 0.528734 ;
235  aVec[203] = 0.519813 ;
236  aVec[204] = 0.512264 ;
237  aVec[205] = 0.507001 ;
238  aVec[206] = 0.49828 ;
239  aVec[207] = 0.492416 ;
240  aVec[208] = 0.483181 ;
241  aVec[209] = 0.477907 ;
242  aVec[210] = 0.469623 ;
243  aVec[211] = 0.462528 ;
244  aVec[212] = 0.455099 ;
245  aVec[213] = 0.45055 ;
246  aVec[214] = 0.443576 ;
247  aVec[215] = 0.435364 ;
248  aVec[216] = 0.429789 ;
249  aVec[217] = 0.422724 ;
250  aVec[218] = 0.415621 ;
251  aVec[219] = 0.409469 ;
252  aVec[220] = 0.40401 ;
253  aVec[221] = 0.398121 ;
254  aVec[222] = 0.391079 ;
255  aVec[223] = 0.384414 ;
256  aVec[224] = 0.378214 ;
257  aVec[225] = 0.369851 ;
258  aVec[226] = 0.365966 ;
259  aVec[227] = 0.359865 ;
260  aVec[228] = 0.353505 ;
261  aVec[229] = 0.347899 ;
262  aVec[230] = 0.343829 ;
263  aVec[231] = 0.337585 ;
264  aVec[232] = 0.333089 ;
265  aVec[233] = 0.326289 ;
266  aVec[234] = 0.322249 ;
267  aVec[235] = 0.316079 ;
268  aVec[236] = 0.31061 ;
269  aVec[237] = 0.305426 ;
270  aVec[238] = 0.301885 ;
271  aVec[239] = 0.296753 ;
272  aVec[240] = 0.290931 ;
273  aVec[241] = 0.286877 ;
274  aVec[242] = 0.281831 ;
275  aVec[243] = 0.276633 ;
276  aVec[244] = 0.272283 ;
277  aVec[245] = 0.268069 ;
278  aVec[246] = 0.26399 ;
279  aVec[247] = 0.258457 ;
280  aVec[248] = 0.253549 ;
281  aVec[249] = 0.249493 ;
282 }
EBShape::fillShape
void fillShape(EcalShapeBase::DVec &aVec) const override
Definition: EBShape.cc:22
mps_fire.i
i
Definition: mps_fire.py:269
EBShape::~EBShape
~EBShape() override
Definition: EBShape.cc:5
EcalShapeBase::DVec
std::vector< double > DVec
Definition: EcalShapeBase.h:16
EBShape::threshold
double threshold() const override
Definition: EBShape.cc:16
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
EBShape::EBShape
EBShape()
Definition: EBShape.cc:9
Generated for CMSSW Reference Manual by doxygen 1.8.11