CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEShape.cc
Go to the documentation of this file.
1 #include <cmath>
2 
3 #include "SimCalorimetry/EcalSimAlgos/interface/EEShape.h"
4 
5 EEShape::~EEShape()
6 {
7 }
8 
9 EEShape::EEShape() :
10  EcalShapeBase( true )
11 {
12  buildMe() ;
13 }
14 
15 double
16 EEShape::threshold() const
17 {
18  return 0.00025 ; // endcap threhsold (largest pedestal val is 0.00019)
19 }
20 
21 void
22 EEShape::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  aVec[i] = exp(2.63163e-01 - 1.37795e-02*(i - 151.0));
29 
30  }
31 
32  aVec[0] = 9.09091e-05;
33  aVec[1] = 9.96016e-05;
34  aVec[2] = 8.08219e-05;
35  aVec[3] = 0.000119685;
36  aVec[4] = 9.56522e-05;
37  aVec[5] = 0.000143969;
38  aVec[6] = 0.000101639;
39  aVec[7] = 0.000115789;
40  aVec[8] = 8.31461e-05;
41  aVec[9] = 0.000117091;
42  aVec[10] = 0.000141224;
43  aVec[11] = 5.38745e-05;
44  aVec[12] = 0.000158079;
45  aVec[13] = 0.000104906;
46  aVec[14] = 7.91946e-05;
47  aVec[15] = 8e-05;
48  aVec[16] = 0.000119703;
49  aVec[17] = 9.96516e-05;
50  aVec[18] = 8.82591e-05;
51  aVec[19] = 0.000146988;
52  aVec[20] = 0.000120588;
53  aVec[21] = 9.89474e-05;
54  aVec[22] = 0.000116788;
55  aVec[23] = 0.000151292;
56  aVec[24] = 0.000108434;
57  aVec[25] = 0.000108651;
58  aVec[26] = 8.5259e-05;
59  aVec[27] = 0.000146575;
60  aVec[28] = 8.66142e-05;
61  aVec[29] = 0.000103679;
62  aVec[30] = 0.000153307;
63  aVec[31] = 7.80328e-05;
64  aVec[32] = 8.42105e-05;
65  aVec[33] = 0.000114607;
66  aVec[34] = 6.90909e-05;
67  aVec[35] = 9.22449e-05;
68  aVec[36] = 8.92989e-05;
69  aVec[37] = 0.000105677;
70  aVec[38] = 9.58491e-05;
71  aVec[39] = 0.000107383;
72  aVec[40] = 8.85714e-05;
73  aVec[41] = 0.000119703;
74  aVec[42] = 0.000141463;
75  aVec[43] = 0.000122267;
76  aVec[44] = 6.0241e-05;
77  aVec[45] = 0.000120588;
78  aVec[46] = 0.000115789;
79  aVec[47] = 0.000121168;
80  aVec[48] = 8.04428e-05;
81  aVec[49] = 9.39759e-05;
82  aVec[50] = 3.391e-05;
83  aVec[51] = -7.72908e-05;
84  aVec[52] = 5.61644e-05;
85  aVec[53] = 0.000181102;
86  aVec[54] = -3.67893e-05;
87  aVec[55] = 6.45914e-05;
88  aVec[56] = 2.68852e-05;
89  aVec[57] = -1.57895e-05;
90  aVec[58] = -2.47191e-05;
91  aVec[59] = -5.30909e-05;
92  aVec[60] = -1.06122e-05;
93  aVec[61] = -3.69004e-06;
94  aVec[62] = 5.85153e-05;
95  aVec[63] = -3.54717e-05;
96  aVec[64] = -9.39597e-06;
97  aVec[65] = 5e-05;
98  aVec[66] = -3.19703e-05;
99  aVec[67] = 3.69338e-05;
100  aVec[68] = -0.000169231;
101  aVec[69] = -4.09639e-05;
102  aVec[70] = -1.61765e-05;
103  aVec[71] = -8.21053e-05;
104  aVec[72] = 6.86131e-05;
105  aVec[73] = -6.12546e-05;
106  aVec[74] = -6.50602e-05;
107  aVec[75] = 2.56055e-05;
108  aVec[76] = -3.90438e-05;
109  aVec[77] = -5.89041e-05;
110  aVec[78] = 0.000195276;
111  aVec[79] = -7.29097e-05;
112  aVec[80] = 1.32296e-05;
113  aVec[81] = 3.27869e-06;
114  aVec[82] = -6.31579e-05;
115  aVec[83] = -3.37079e-05;
116  aVec[84] = 5.16364e-05;
117  aVec[85] = -4.97959e-05;
118  aVec[86] = 8.92989e-05;
119  aVec[87] = 0.000561572;
120  aVec[88] = 0.00171698;
121  aVec[89] = 0.0035906;
122  aVec[90] = 0.00742143;
123  aVec[91] = 0.01249;
124  aVec[92] = 0.020922;
125  aVec[93] = 0.030336;
126  aVec[94] = 0.043588;
127  aVec[95] = 0.0593529;
128  aVec[96] = 0.0764947;
129  aVec[97] = 0.0958978;
130  aVec[98] = 0.11567;
131  aVec[99] = 0.137087;
132  aVec[100] = 0.160643;
133  aVec[101] = 0.188973;
134  aVec[102] = 0.218062;
135  aVec[103] = 0.240682;
136  aVec[104] = 0.265773;
137  aVec[105] = 0.293327;
138  aVec[106] = 0.320832;
139  aVec[107] = 0.348689;
140  aVec[108] = 0.374411;
141  aVec[109] = 0.403575;
142  aVec[110] = 0.431147;
143  aVec[111] = 0.458644;
144  aVec[112] = 0.488217;
145  aVec[113] = 0.515254;
146  aVec[114] = 0.53902;
147  aVec[115] = 0.566819;
148  aVec[116] = 0.589752;
149  aVec[117] = 0.615587;
150  aVec[118] = 0.63808;
151  aVec[119] = 0.664137;
152  aVec[120] = 0.687959;
153  aVec[121] = 0.707379;
154  aVec[122] = 0.726826;
155  aVec[123] = 0.747043;
156  aVec[124] = 0.762952;
157  aVec[125] = 0.781342;
158  aVec[126] = 0.799515;
159  aVec[127] = 0.818559;
160  aVec[128] = 0.832301;
161  aVec[129] = 0.846658;
162  aVec[130] = 0.860791;
163  aVec[131] = 0.874009;
164  aVec[132] = 0.886242;
165  aVec[133] = 0.896982;
166  aVec[134] = 0.908713;
167  aVec[135] = 0.9189;
168  aVec[136] = 0.9283;
169  aVec[137] = 0.938453;
170  aVec[138] = 0.945995;
171  aVec[139] = 0.95186;
172  aVec[140] = 0.958876;
173  aVec[141] = 0.963283;
174  aVec[142] = 0.968989;
175  aVec[143] = 0.972846;
176  aVec[144] = 0.977824;
177  aVec[145] = 0.981341;
178  aVec[146] = 0.982743;
179  aVec[147] = 0.984082;
180  aVec[148] = 0.986609;
181  aVec[149] = 0.985925;
182  aVec[150] = 0.984707;
183  aVec[151] = 0.981384;
184  aVec[152] = 0.980941;
185  aVec[153] = 0.980132;
186  aVec[154] = 0.978176;
187  aVec[155] = 0.975987;
188  aVec[156] = 0.973396;
189  aVec[157] = 0.970395;
190  aVec[158] = 0.967422;
191  aVec[159] = 0.96346;
192  aVec[160] = 0.959622;
193  aVec[161] = 0.955396;
194  aVec[162] = 0.950846;
195  aVec[163] = 0.945995;
196  aVec[164] = 0.940835;
197  aVec[165] = 0.935373;
198  aVec[166] = 0.929544;
199  aVec[167] = 0.923836;
200  aVec[168] = 0.91787;
201  aVec[169] = 0.910812;
202  aVec[170] = 0.903924;
203  aVec[171] = 0.897669;
204  aVec[172] = 0.89127;
205  aVec[173] = 0.884613;
206  aVec[174] = 0.878051;
207  aVec[175] = 0.867734;
208  aVec[176] = 0.856756;
209  aVec[177] = 0.848477;
210  aVec[178] = 0.842189;
211  aVec[179] = 0.833787;
212  aVec[180] = 0.825903;
213  aVec[181] = 0.817931;
214  aVec[182] = 0.8102;
215  aVec[183] = 0.802771;
216  aVec[184] = 0.79388;
217  aVec[185] = 0.785348;
218  aVec[186] = 0.777295;
219  aVec[187] = 0.767854;
220  aVec[188] = 0.759619;
221  aVec[189] = 0.752358;
222  aVec[190] = 0.743394;
223  aVec[191] = 0.736674;
224  aVec[192] = 0.726982;
225  aVec[193] = 0.719321;
226  aVec[194] = 0.708937;
227  aVec[195] = 0.700003;
228  aVec[196] = 0.693202;
229  aVec[197] = 0.685961;
230  aVec[198] = 0.676889;
231  aVec[199] = 0.670513;
232  aVec[200] = 0.658788;
233  aVec[201] = 0.648343;
234  aVec[202] = 0.639381;
235  aVec[203] = 0.632761;
236  aVec[204] = 0.623614;
237  aVec[205] = 0.616164;
238  aVec[206] = 0.607884;
239  aVec[207] = 0.601437;
240  aVec[208] = 0.593944;
241  aVec[209] = 0.585817;
242  aVec[210] = 0.577248;
243  aVec[211] = 0.570328;
244  aVec[212] = 0.56147;
245  aVec[213] = 0.553509;
246  aVec[214] = 0.548765;
247  aVec[215] = 0.539964;
248  aVec[216] = 0.53343;
249  aVec[217] = 0.525787;
250  aVec[218] = 0.519868;
251  aVec[219] = 0.510817;
252  aVec[220] = 0.503274;
253  aVec[221] = 0.498478;
254  aVec[222] = 0.493058;
255  aVec[223] = 0.484176;
256  aVec[224] = 0.479896;
257  aVec[225] = 0.469998;
258  aVec[226] = 0.462893;
259  aVec[227] = 0.45576;
260  aVec[228] = 0.450375;
261  aVec[229] = 0.442683;
262  aVec[230] = 0.437887;
263  aVec[231] = 0.430736;
264  aVec[232] = 0.427216;
265  aVec[233] = 0.420775;
266  aVec[234] = 0.41485;
267  aVec[235] = 0.407769;
268  aVec[236] = 0.403546;
269  aVec[237] = 0.397028;
270  aVec[238] = 0.389534;
271  aVec[239] = 0.388569;
272  aVec[240] = 0.380296;
273  aVec[241] = 0.375312;
274  aVec[242] = 0.369457;
275  aVec[243] = 0.365884;
276  aVec[244] = 0.358889;
277  aVec[245] = 0.353344;
278  aVec[246] = 0.350364;
279  aVec[247] = 0.347136;
280  aVec[248] = 0.33872;
281  aVec[249] = 0.336706;
282 
283 }
i
int i
Definition: DBlmapReader.cc:9
EEShape::threshold
virtual double threshold() const
Definition: EEShape.cc:16
EEShape::EEShape
EEShape()
Definition: EEShape.cc:9
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
EEShape::~EEShape
virtual ~EEShape()
Definition: EEShape.cc:5
EcalShapeBase::DVec
std::vector< double > DVec
Definition: EcalShapeBase.h:16
EEShape::fillShape
virtual void fillShape(EcalShapeBase::DVec &aVec) const
Definition: EEShape.cc:22
Generated for CMSSW Reference Manual by doxygen 1.8.5