Go to the documentation of this file.00001 #include <cmath>
00002
00003 #include "SimCalorimetry/EcalSimAlgos/interface/EBShape.h"
00004
00005 EBShape::~EBShape()
00006 {
00007 }
00008
00009 EBShape::EBShape() :
00010 EcalShapeBase( true )
00011 {
00012 buildMe() ;
00013 }
00014
00015 double
00016 EBShape::threshold() const
00017 {
00018 return 0.00013 ;
00019 }
00020
00021 void
00022 EBShape::fillShape( EcalShapeBase::DVec& aVec ) const
00023 {
00024
00025 for( unsigned int i ( k1NSecBins ) ; i != k1NSecBinsTotal ; ++i )
00026 {
00027 aVec[i] = exp(2.39735 - 0.0151053* ((double)i+1.0));
00028
00029 }
00030
00031
00032 aVec[0] = 6.94068e-05 ;
00033 aVec[1] = -5.03304e-05 ;
00034 aVec[2] = -2.13404e-05 ;
00035 aVec[3] = 6.017e-05 ;
00036 aVec[4] = 2.01697e-05 ;
00037 aVec[5] = 0.000114845 ;
00038 aVec[6] = 2.13998e-05 ;
00039 aVec[7] = 2.74476e-05 ;
00040 aVec[8] = 5.2824e-05 ;
00041 aVec[9] = 8.754e-05 ;
00042 aVec[10] = 2.95346e-06 ;
00043 aVec[11] = -7.58699e-05 ;
00044 aVec[12] = -2.72224e-05 ;
00045 aVec[13] = 3.10997e-06 ;
00046 aVec[14] = -3.97771e-05 ;
00047 aVec[15] = -1.06916e-05 ;
00048 aVec[16] = -0.000113865 ;
00049 aVec[17] = 6.05044e-05 ;
00050 aVec[18] = -5.81202e-05 ;
00051 aVec[19] = -6.58974e-06 ;
00052 aVec[20] = 5.37494e-05 ;
00053 aVec[21] = -0.000123729 ;
00054 aVec[22] = 7.50938e-06 ;
00055 aVec[23] = -1.35628e-05 ;
00056 aVec[24] = 8.33725e-05 ;
00057 aVec[25] = 3.19299e-05 ;
00058 aVec[26] = -3.09232e-05 ;
00059 aVec[27] = -7.0086e-05 ;
00060 aVec[28] = 1.78937e-06 ;
00061 aVec[29] = -2.20365e-05 ;
00062 aVec[30] = 7.68054e-05 ;
00063 aVec[31] = -2.5368e-05 ;
00064 aVec[32] = 5.67291e-06 ;
00065 aVec[33] = 5.87096e-05 ;
00066 aVec[34] = -2.62771e-06 ;
00067 aVec[35] = 4.31832e-05 ;
00068 aVec[36] = 8.33616e-06 ;
00069 aVec[37] = 7.27813e-05 ;
00070 aVec[38] = 7.6159e-05 ;
00071 aVec[39] = -1.60446e-05 ;
00072 aVec[40] = -4.12127e-06 ;
00073 aVec[41] = -5.93381e-05 ;
00074 aVec[42] = 1.61444e-05 ;
00075 aVec[43] = -5.49559e-05 ;
00076 aVec[44] = 5.55254e-05 ;
00077 aVec[45] = 3.32251e-05 ;
00078 aVec[46] = -3.15897e-05 ;
00079 aVec[47] = 7.86588e-05 ;
00080 aVec[48] = -2.9704e-05 ;
00081 aVec[49] = 5.66838e-05 ;
00082 aVec[50] = 2.85281e-05 ;
00083 aVec[51] = -3.02436e-05 ;
00084 aVec[52] = -4.16265e-05 ;
00085 aVec[53] = -1.63191e-05 ;
00086 aVec[54] = 6.61193e-05 ;
00087 aVec[55] = 9.23766e-05 ;
00088 aVec[56] = 6.68903e-05 ;
00089 aVec[57] = -3.20994e-05 ;
00090 aVec[58] = 0.00011082 ;
00091 aVec[59] = -4.07997e-05 ;
00092 aVec[60] = -8.29046e-06 ;
00093 aVec[61] = -7.42197e-05 ;
00094 aVec[62] = -1.64386e-05 ;
00095 aVec[63] = 1.02508e-05 ;
00096 aVec[64] = 7.10995e-06 ;
00097 aVec[65] = -5.87486e-05 ;
00098 aVec[66] = -0.000101201 ;
00099 aVec[67] = 1.62003e-05 ;
00100 aVec[68] = -2.53093e-05 ;
00101 aVec[69] = 2.65239e-05 ;
00102 aVec[70] = -2.68722e-05 ;
00103 aVec[71] = -4.02001e-05 ;
00104 aVec[72] = 5.0674e-05 ;
00105 aVec[73] = -1.75884e-05 ;
00106 aVec[74] = 4.7902e-05 ;
00107 aVec[75] = -1.01079e-05 ;
00108 aVec[76] = 1.08427e-05 ;
00109 aVec[77] = -0.000112906 ;
00110 aVec[78] = 3.33076e-05 ;
00111 aVec[79] = 0.000181201 ;
00112 aVec[80] = 0.000426875 ;
00113 aVec[81] = 0.00114222 ;
00114 aVec[82] = 0.00237804 ;
00115 aVec[83] = 0.00541858 ;
00116 aVec[84] = 0.0089021 ;
00117 aVec[85] = 0.0149157 ;
00118 aVec[86] = 0.0231397 ;
00119 aVec[87] = 0.0344671 ;
00120 aVec[88] = 0.0471013 ;
00121 aVec[89] = 0.0625517 ;
00122 aVec[90] = 0.0857351 ;
00123 aVec[91] = 0.108561 ;
00124 aVec[92] = 0.133481 ;
00125 aVec[93] = 0.163557 ;
00126 aVec[94] = 0.200243 ;
00127 aVec[95] = 0.225919 ;
00128 aVec[96] = 0.269213 ;
00129 aVec[97] = 0.302929 ;
00130 aVec[98] = 0.342722 ;
00131 aVec[99] = 0.378522 ;
00132 aVec[100] = 0.436563 ;
00133 aVec[101] = 0.467581 ;
00134 aVec[102] = 0.510133 ;
00135 aVec[103] = 0.550063 ;
00136 aVec[104] = 0.583509 ;
00137 aVec[105] = 0.619187 ;
00138 aVec[106] = 0.653245 ;
00139 aVec[107] = 0.686101 ;
00140 aVec[108] = 0.721178 ;
00141 aVec[109] = 0.745129 ;
00142 aVec[110] = 0.774163 ;
00143 aVec[111] = 0.799011 ;
00144 aVec[112] = 0.822177 ;
00145 aVec[113] = 0.838315 ;
00146 aVec[114] = 0.858847 ;
00147 aVec[115] = 0.875559 ;
00148 aVec[116] = 0.891294 ;
00149 aVec[117] = 0.90537 ;
00150 aVec[118] = 0.919617 ;
00151 aVec[119] = 0.930632 ;
00152 aVec[120] = 0.936216 ;
00153 aVec[121] = 0.947739 ;
00154 aVec[122] = 0.955306 ;
00155 aVec[123] = 0.961876 ;
00156 aVec[124] = 0.968124 ;
00157 aVec[125] = 0.97327 ;
00158 aVec[126] = 0.977513 ;
00159 aVec[127] = 0.984885 ;
00160 aVec[128] = 0.986497 ;
00161 aVec[129] = 0.990039 ;
00162 aVec[130] = 0.994798 ;
00163 aVec[131] = 0.994884 ;
00164 aVec[132] = 0.99795 ;
00165 aVec[133] = 0.99834 ;
00166 aVec[134] = 0.999607 ;
00167 aVec[135] = 1 ;
00168 aVec[136] = 0.999047 ;
00169 aVec[137] = 0.998745 ;
00170 aVec[138] = 0.999219 ;
00171 aVec[139] = 0.99814 ;
00172 aVec[140] = 0.995082 ;
00173 aVec[141] = 0.992449 ;
00174 aVec[142] = 0.990418 ;
00175 aVec[143] = 0.985032 ;
00176 aVec[144] = 0.982308 ;
00177 aVec[145] = 0.978696 ;
00178 aVec[146] = 0.975656 ;
00179 aVec[147] = 0.971027 ;
00180 aVec[148] = 0.964811 ;
00181 aVec[149] = 0.959428 ;
00182 aVec[150] = 0.95096 ;
00183 aVec[151] = 0.947428 ;
00184 aVec[152] = 0.9419 ;
00185 aVec[153] = 0.933223 ;
00186 aVec[154] = 0.926482 ;
00187 aVec[155] = 0.922172 ;
00188 aVec[156] = 0.912777 ;
00189 aVec[157] = 0.907388 ;
00190 aVec[158] = 0.897289 ;
00191 aVec[159] = 0.891889 ;
00192 aVec[160] = 0.882056 ;
00193 aVec[161] = 0.873382 ;
00194 aVec[162] = 0.865442 ;
00195 aVec[163] = 0.860032 ;
00196 aVec[164] = 0.85202 ;
00197 aVec[165] = 0.841013 ;
00198 aVec[166] = 0.833802 ;
00199 aVec[167] = 0.825259 ;
00200 aVec[168] = 0.815013 ;
00201 aVec[169] = 0.807465 ;
00202 aVec[170] = 0.799428 ;
00203 aVec[171] = 0.792165 ;
00204 aVec[172] = 0.783088 ;
00205 aVec[173] = 0.773392 ;
00206 aVec[174] = 0.764982 ;
00207 aVec[175] = 0.752174 ;
00208 aVec[176] = 0.746487 ;
00209 aVec[177] = 0.737678 ;
00210 aVec[178] = 0.727396 ;
00211 aVec[179] = 0.718692 ;
00212 aVec[180] = 0.712737 ;
00213 aVec[181] = 0.702738 ;
00214 aVec[182] = 0.69559 ;
00215 aVec[183] = 0.684389 ;
00216 aVec[184] = 0.677989 ;
00217 aVec[185] = 0.667643 ;
00218 aVec[186] = 0.659009 ;
00219 aVec[187] = 0.650217 ;
00220 aVec[188] = 0.644479 ;
00221 aVec[189] = 0.636017 ;
00222 aVec[190] = 0.625257 ;
00223 aVec[191] = 0.618507 ;
00224 aVec[192] = 0.609798 ;
00225 aVec[193] = 0.600097 ;
00226 aVec[194] = 0.592788 ;
00227 aVec[195] = 0.584895 ;
00228 aVec[196] = 0.578228 ;
00229 aVec[197] = 0.569299 ;
00230 aVec[198] = 0.560576 ;
00231 aVec[199] = 0.552404 ;
00232 aVec[200] = 0.541405 ;
00233 aVec[201] = 0.536271 ;
00234 aVec[202] = 0.528734 ;
00235 aVec[203] = 0.519813 ;
00236 aVec[204] = 0.512264 ;
00237 aVec[205] = 0.507001 ;
00238 aVec[206] = 0.49828 ;
00239 aVec[207] = 0.492416 ;
00240 aVec[208] = 0.483181 ;
00241 aVec[209] = 0.477907 ;
00242 aVec[210] = 0.469623 ;
00243 aVec[211] = 0.462528 ;
00244 aVec[212] = 0.455099 ;
00245 aVec[213] = 0.45055 ;
00246 aVec[214] = 0.443576 ;
00247 aVec[215] = 0.435364 ;
00248 aVec[216] = 0.429789 ;
00249 aVec[217] = 0.422724 ;
00250 aVec[218] = 0.415621 ;
00251 aVec[219] = 0.409469 ;
00252 aVec[220] = 0.40401 ;
00253 aVec[221] = 0.398121 ;
00254 aVec[222] = 0.391079 ;
00255 aVec[223] = 0.384414 ;
00256 aVec[224] = 0.378214 ;
00257 aVec[225] = 0.369851 ;
00258 aVec[226] = 0.365966 ;
00259 aVec[227] = 0.359865 ;
00260 aVec[228] = 0.353505 ;
00261 aVec[229] = 0.347899 ;
00262 aVec[230] = 0.343829 ;
00263 aVec[231] = 0.337585 ;
00264 aVec[232] = 0.333089 ;
00265 aVec[233] = 0.326289 ;
00266 aVec[234] = 0.322249 ;
00267 aVec[235] = 0.316079 ;
00268 aVec[236] = 0.31061 ;
00269 aVec[237] = 0.305426 ;
00270 aVec[238] = 0.301885 ;
00271 aVec[239] = 0.296753 ;
00272 aVec[240] = 0.290931 ;
00273 aVec[241] = 0.286877 ;
00274 aVec[242] = 0.281831 ;
00275 aVec[243] = 0.276633 ;
00276 aVec[244] = 0.272283 ;
00277 aVec[245] = 0.268069 ;
00278 aVec[246] = 0.26399 ;
00279 aVec[247] = 0.258457 ;
00280 aVec[248] = 0.253549 ;
00281 aVec[249] = 0.249493 ;
00282 }
1.7.3