Go to the documentation of this file.00001 #include <cmath>
00002
00003 #include "SimCalorimetry/EcalSimAlgos/interface/EEShape.h"
00004
00005 EEShape::~EEShape()
00006 {
00007 }
00008
00009 EEShape::EEShape() :
00010 EcalShapeBase( true )
00011 {
00012 buildMe() ;
00013 }
00014
00015 double
00016 EEShape::threshold() const
00017 {
00018 return 0.00025 ;
00019 }
00020
00021 void
00022 EEShape::fillShape( EcalShapeBase::DVec& aVec ) const
00023 {
00024
00025 for( unsigned int i ( k1NSecBins ) ; i != k1NSecBinsTotal ; ++i )
00026 {
00027
00028 aVec[i] = exp(2.63163e-01 - 1.37795e-02*(i - 151.0));
00029
00030 }
00031
00032 aVec[0] = 9.09091e-05;
00033 aVec[1] = 9.96016e-05;
00034 aVec[2] = 8.08219e-05;
00035 aVec[3] = 0.000119685;
00036 aVec[4] = 9.56522e-05;
00037 aVec[5] = 0.000143969;
00038 aVec[6] = 0.000101639;
00039 aVec[7] = 0.000115789;
00040 aVec[8] = 8.31461e-05;
00041 aVec[9] = 0.000117091;
00042 aVec[10] = 0.000141224;
00043 aVec[11] = 5.38745e-05;
00044 aVec[12] = 0.000158079;
00045 aVec[13] = 0.000104906;
00046 aVec[14] = 7.91946e-05;
00047 aVec[15] = 8e-05;
00048 aVec[16] = 0.000119703;
00049 aVec[17] = 9.96516e-05;
00050 aVec[18] = 8.82591e-05;
00051 aVec[19] = 0.000146988;
00052 aVec[20] = 0.000120588;
00053 aVec[21] = 9.89474e-05;
00054 aVec[22] = 0.000116788;
00055 aVec[23] = 0.000151292;
00056 aVec[24] = 0.000108434;
00057 aVec[25] = 0.000108651;
00058 aVec[26] = 8.5259e-05;
00059 aVec[27] = 0.000146575;
00060 aVec[28] = 8.66142e-05;
00061 aVec[29] = 0.000103679;
00062 aVec[30] = 0.000153307;
00063 aVec[31] = 7.80328e-05;
00064 aVec[32] = 8.42105e-05;
00065 aVec[33] = 0.000114607;
00066 aVec[34] = 6.90909e-05;
00067 aVec[35] = 9.22449e-05;
00068 aVec[36] = 8.92989e-05;
00069 aVec[37] = 0.000105677;
00070 aVec[38] = 9.58491e-05;
00071 aVec[39] = 0.000107383;
00072 aVec[40] = 8.85714e-05;
00073 aVec[41] = 0.000119703;
00074 aVec[42] = 0.000141463;
00075 aVec[43] = 0.000122267;
00076 aVec[44] = 6.0241e-05;
00077 aVec[45] = 0.000120588;
00078 aVec[46] = 0.000115789;
00079 aVec[47] = 0.000121168;
00080 aVec[48] = 8.04428e-05;
00081 aVec[49] = 9.39759e-05;
00082 aVec[50] = 3.391e-05;
00083 aVec[51] = -7.72908e-05;
00084 aVec[52] = 5.61644e-05;
00085 aVec[53] = 0.000181102;
00086 aVec[54] = -3.67893e-05;
00087 aVec[55] = 6.45914e-05;
00088 aVec[56] = 2.68852e-05;
00089 aVec[57] = -1.57895e-05;
00090 aVec[58] = -2.47191e-05;
00091 aVec[59] = -5.30909e-05;
00092 aVec[60] = -1.06122e-05;
00093 aVec[61] = -3.69004e-06;
00094 aVec[62] = 5.85153e-05;
00095 aVec[63] = -3.54717e-05;
00096 aVec[64] = -9.39597e-06;
00097 aVec[65] = 5e-05;
00098 aVec[66] = -3.19703e-05;
00099 aVec[67] = 3.69338e-05;
00100 aVec[68] = -0.000169231;
00101 aVec[69] = -4.09639e-05;
00102 aVec[70] = -1.61765e-05;
00103 aVec[71] = -8.21053e-05;
00104 aVec[72] = 6.86131e-05;
00105 aVec[73] = -6.12546e-05;
00106 aVec[74] = -6.50602e-05;
00107 aVec[75] = 2.56055e-05;
00108 aVec[76] = -3.90438e-05;
00109 aVec[77] = -5.89041e-05;
00110 aVec[78] = 0.000195276;
00111 aVec[79] = -7.29097e-05;
00112 aVec[80] = 1.32296e-05;
00113 aVec[81] = 3.27869e-06;
00114 aVec[82] = -6.31579e-05;
00115 aVec[83] = -3.37079e-05;
00116 aVec[84] = 5.16364e-05;
00117 aVec[85] = -4.97959e-05;
00118 aVec[86] = 8.92989e-05;
00119 aVec[87] = 0.000561572;
00120 aVec[88] = 0.00171698;
00121 aVec[89] = 0.0035906;
00122 aVec[90] = 0.00742143;
00123 aVec[91] = 0.01249;
00124 aVec[92] = 0.020922;
00125 aVec[93] = 0.030336;
00126 aVec[94] = 0.043588;
00127 aVec[95] = 0.0593529;
00128 aVec[96] = 0.0764947;
00129 aVec[97] = 0.0958978;
00130 aVec[98] = 0.11567;
00131 aVec[99] = 0.137087;
00132 aVec[100] = 0.160643;
00133 aVec[101] = 0.188973;
00134 aVec[102] = 0.218062;
00135 aVec[103] = 0.240682;
00136 aVec[104] = 0.265773;
00137 aVec[105] = 0.293327;
00138 aVec[106] = 0.320832;
00139 aVec[107] = 0.348689;
00140 aVec[108] = 0.374411;
00141 aVec[109] = 0.403575;
00142 aVec[110] = 0.431147;
00143 aVec[111] = 0.458644;
00144 aVec[112] = 0.488217;
00145 aVec[113] = 0.515254;
00146 aVec[114] = 0.53902;
00147 aVec[115] = 0.566819;
00148 aVec[116] = 0.589752;
00149 aVec[117] = 0.615587;
00150 aVec[118] = 0.63808;
00151 aVec[119] = 0.664137;
00152 aVec[120] = 0.687959;
00153 aVec[121] = 0.707379;
00154 aVec[122] = 0.726826;
00155 aVec[123] = 0.747043;
00156 aVec[124] = 0.762952;
00157 aVec[125] = 0.781342;
00158 aVec[126] = 0.799515;
00159 aVec[127] = 0.818559;
00160 aVec[128] = 0.832301;
00161 aVec[129] = 0.846658;
00162 aVec[130] = 0.860791;
00163 aVec[131] = 0.874009;
00164 aVec[132] = 0.886242;
00165 aVec[133] = 0.896982;
00166 aVec[134] = 0.908713;
00167 aVec[135] = 0.9189;
00168 aVec[136] = 0.9283;
00169 aVec[137] = 0.938453;
00170 aVec[138] = 0.945995;
00171 aVec[139] = 0.95186;
00172 aVec[140] = 0.958876;
00173 aVec[141] = 0.963283;
00174 aVec[142] = 0.968989;
00175 aVec[143] = 0.972846;
00176 aVec[144] = 0.977824;
00177 aVec[145] = 0.981341;
00178 aVec[146] = 0.982743;
00179 aVec[147] = 0.984082;
00180 aVec[148] = 0.986609;
00181 aVec[149] = 0.985925;
00182 aVec[150] = 0.984707;
00183 aVec[151] = 0.981384;
00184 aVec[152] = 0.980941;
00185 aVec[153] = 0.980132;
00186 aVec[154] = 0.978176;
00187 aVec[155] = 0.975987;
00188 aVec[156] = 0.973396;
00189 aVec[157] = 0.970395;
00190 aVec[158] = 0.967422;
00191 aVec[159] = 0.96346;
00192 aVec[160] = 0.959622;
00193 aVec[161] = 0.955396;
00194 aVec[162] = 0.950846;
00195 aVec[163] = 0.945995;
00196 aVec[164] = 0.940835;
00197 aVec[165] = 0.935373;
00198 aVec[166] = 0.929544;
00199 aVec[167] = 0.923836;
00200 aVec[168] = 0.91787;
00201 aVec[169] = 0.910812;
00202 aVec[170] = 0.903924;
00203 aVec[171] = 0.897669;
00204 aVec[172] = 0.89127;
00205 aVec[173] = 0.884613;
00206 aVec[174] = 0.878051;
00207 aVec[175] = 0.867734;
00208 aVec[176] = 0.856756;
00209 aVec[177] = 0.848477;
00210 aVec[178] = 0.842189;
00211 aVec[179] = 0.833787;
00212 aVec[180] = 0.825903;
00213 aVec[181] = 0.817931;
00214 aVec[182] = 0.8102;
00215 aVec[183] = 0.802771;
00216 aVec[184] = 0.79388;
00217 aVec[185] = 0.785348;
00218 aVec[186] = 0.777295;
00219 aVec[187] = 0.767854;
00220 aVec[188] = 0.759619;
00221 aVec[189] = 0.752358;
00222 aVec[190] = 0.743394;
00223 aVec[191] = 0.736674;
00224 aVec[192] = 0.726982;
00225 aVec[193] = 0.719321;
00226 aVec[194] = 0.708937;
00227 aVec[195] = 0.700003;
00228 aVec[196] = 0.693202;
00229 aVec[197] = 0.685961;
00230 aVec[198] = 0.676889;
00231 aVec[199] = 0.670513;
00232 aVec[200] = 0.658788;
00233 aVec[201] = 0.648343;
00234 aVec[202] = 0.639381;
00235 aVec[203] = 0.632761;
00236 aVec[204] = 0.623614;
00237 aVec[205] = 0.616164;
00238 aVec[206] = 0.607884;
00239 aVec[207] = 0.601437;
00240 aVec[208] = 0.593944;
00241 aVec[209] = 0.585817;
00242 aVec[210] = 0.577248;
00243 aVec[211] = 0.570328;
00244 aVec[212] = 0.56147;
00245 aVec[213] = 0.553509;
00246 aVec[214] = 0.548765;
00247 aVec[215] = 0.539964;
00248 aVec[216] = 0.53343;
00249 aVec[217] = 0.525787;
00250 aVec[218] = 0.519868;
00251 aVec[219] = 0.510817;
00252 aVec[220] = 0.503274;
00253 aVec[221] = 0.498478;
00254 aVec[222] = 0.493058;
00255 aVec[223] = 0.484176;
00256 aVec[224] = 0.479896;
00257 aVec[225] = 0.469998;
00258 aVec[226] = 0.462893;
00259 aVec[227] = 0.45576;
00260 aVec[228] = 0.450375;
00261 aVec[229] = 0.442683;
00262 aVec[230] = 0.437887;
00263 aVec[231] = 0.430736;
00264 aVec[232] = 0.427216;
00265 aVec[233] = 0.420775;
00266 aVec[234] = 0.41485;
00267 aVec[235] = 0.407769;
00268 aVec[236] = 0.403546;
00269 aVec[237] = 0.397028;
00270 aVec[238] = 0.389534;
00271 aVec[239] = 0.388569;
00272 aVec[240] = 0.380296;
00273 aVec[241] = 0.375312;
00274 aVec[242] = 0.369457;
00275 aVec[243] = 0.365884;
00276 aVec[244] = 0.358889;
00277 aVec[245] = 0.353344;
00278 aVec[246] = 0.350364;
00279 aVec[247] = 0.347136;
00280 aVec[248] = 0.33872;
00281 aVec[249] = 0.336706;
00282
00283 }
1.7.3