Reference Manual

Go to the documentation of this file.
00001 #include "SimCalorimetry/HcalSimAlgos/interface/HcalSiPM.h"
00002 #include "FWCore/Utilities/interface/RandomNumberGenerator.h"
00003 #include "FWCore/ServiceRegistry/interface/Service.h"
00004 
00005 #include <cmath>
00006 
00007 using std::vector;
00008 
00009 HcalSiPM::HcalSiPM(int nCells) :
00010   theCellCount(nCells), theRndGauss(0) {
00011 
00012   assert(theCellCount>0);
00013 
00014   defaultErrInit();
00015 }
00016 
00017 HcalSiPM::~HcalSiPM() {
00018   delete theRndGauss;
00019 }
00020 
00021 int HcalSiPM::hitCells(int photons, int integral) const {
00022   //don't need to do zero or negative photons.
00023   if (photons < 1) return 0;
00024   if (integral < 0) integral = 0;
00025   if (integral >= theCellCount) return 0;
00026 
00027   if (theRndGauss == 0) {
00028     //random number generator setup
00029     edm::Service<edm::RandomNumberGenerator> rng;
00030     if ( ! rng.isAvailable()) {
00031       throw cms::Exception("Configuration")
00032         << "HcalSiPM requires the RandomNumberGeneratorService\n"
00033         "which is not present in the configuration file.  "
00034         "You must add the service\n"
00035         "in the configuration file or remove the modules that require it.";
00036     }
00037 
00038     CLHEP::HepRandomEngine& engine = rng->getEngine();
00039     theRndGauss = new CLHEP::RandGaussQ(engine);
00040   }
00041 
00042   //normalize by theCellCount to remove dependency on SiPM size and pixel density.
00043   double x = double(photons)/double(theCellCount);
00044   double prehit = double(integral)/double(theCellCount);
00045 
00046   //calculate the width and mean of the distribution for a given x
00047   double width = errOnX(x, prehit);
00048   double mean = 1. - std::exp(-x);
00049 
00050   //you can't hit more than everything.
00051   if (mean > 1.) mean = 1.;
00052 
00053   //convert back to absolute pixels
00054   mean *= (1-prehit)*theCellCount;
00055   width *= (1-prehit)*theCellCount;
00056 
00057   double npe;
00058   while (true) {
00059     npe = theRndGauss->fire(mean, width);
00060     if ((npe > -0.5) && (npe <= theCellCount-integral) && (npe <= photons))
00061       return int(npe + 0.5);
00062   }
00063 }
00064 
00065 void HcalSiPM::setNCells(int nCells) {
00066   assert(nCells>0);
00067   theCellCount = nCells;
00068 }
00069 
00070 void HcalSiPM::initRandomEngine(CLHEP::HepRandomEngine& engine) {
00071   if(theRndGauss) delete theRndGauss;
00072   theRndGauss = new CLHEP::RandGaussQ(engine);
00073 }
00074 
00075 double HcalSiPM::errOnX(double x, double prehit) const {
00076   //interpolate between the points in the x,prehit space to determine the
00077   //width of the distribution of x
00078   int before, after;
00079   getBeforeAndAfter(x, before, after, theXSamples);
00080   int i, j;
00081   getBeforeAndAfter(prehit, i, j, thePrehitSamples);
00082   double m1 = (theErrSamples[i][after]-theErrSamples[i][before])/
00083     (theXSamples[after]-theXSamples[before]);
00084   double err1 = theErrSamples[i][before]+m1*(x-theXSamples[before]);
00085   double m2 = (theErrSamples[j][after]-theErrSamples[j][before])/
00086     (theXSamples[after]-theXSamples[before]);
00087   double err2 = theErrSamples[j][before]+m2*(x-theXSamples[before]);
00088 
00089   return err1 + (err2-err1)/(thePrehitSamples[j]-thePrehitSamples[i])*
00090     (prehit - thePrehitSamples[i]);
00091 }
00092 
00093 void HcalSiPM::getBeforeAndAfter(double val, int& before, int& after, 
00094                                  const vector<double>& vec) const {
00095   //binary like search to indeces points around an arbitrary value.
00096   before = 0;
00097   after = vec.size()-1;
00098   int mid = vec.size()/2;
00099   bool go = true;
00100   do {
00101     if (vec[mid]>val) after = mid;
00102     else before = mid;
00103     mid = (after+before)/2;
00104     go = (after-before > 1);
00105   } while (go);
00106 }
00107 
00108 void HcalSiPM::defaultErrInit() {
00109   //load up values for the interpolation for the errs on x used to generate
00110   //the number of photons.  These were produced from a toy SiPM model.
00111   theXSamples.clear();
00112   thePrehitSamples.clear();
00113   theErrSamples.clear();
00114 
00115   theXSamples.push_back(7.14286e-05);
00116   theXSamples.push_back(0.000142857);
00117   theXSamples.push_back(0.000428571);
00118   theXSamples.push_back(0.00107143);
00119   theXSamples.push_back(0.00221429);
00120   theXSamples.push_back(0.004);
00121   theXSamples.push_back(0.00657143);
00122   theXSamples.push_back(0.0100714);
00123   theXSamples.push_back(0.0146429);
00124   theXSamples.push_back(0.0204286);
00125   theXSamples.push_back(0.0275714);
00126   theXSamples.push_back(0.0362143);
00127   theXSamples.push_back(0.0465);
00128   theXSamples.push_back(0.0585714);
00129   theXSamples.push_back(0.0725714);
00130   theXSamples.push_back(0.0886429);
00131   theXSamples.push_back(0.106929);
00132   theXSamples.push_back(0.127571);
00133   theXSamples.push_back(0.150714);
00134   theXSamples.push_back(0.1765);
00135   theXSamples.push_back(0.205071);
00136   theXSamples.push_back(0.236571);
00137   theXSamples.push_back(0.271143);
00138   theXSamples.push_back(0.308929);
00139   theXSamples.push_back(0.350071);
00140   theXSamples.push_back(0.394714);
00141   theXSamples.push_back(0.443);
00142   theXSamples.push_back(0.495071);
00143   theXSamples.push_back(0.551071);
00144   theXSamples.push_back(0.611143);
00145   theXSamples.push_back(0.675429);
00146   theXSamples.push_back(0.744071);
00147   theXSamples.push_back(0.817214);
00148   theXSamples.push_back(0.895);
00149   theXSamples.push_back(0.977571);
00150   theXSamples.push_back(1.06507);
00151   theXSamples.push_back(1.15764);
00152   theXSamples.push_back(1.25543);
00153   theXSamples.push_back(1.35857);
00154   theXSamples.push_back(1.46721);
00155   theXSamples.push_back(1.5815);
00156   theXSamples.push_back(1.70157);
00157   theXSamples.push_back(1.82757);
00158   theXSamples.push_back(1.95964);
00159   theXSamples.push_back(2.09793);
00160   theXSamples.push_back(2.24257);
00161   theXSamples.push_back(2.39371);
00162   theXSamples.push_back(2.5515);
00163   theXSamples.push_back(2.71607);
00164   theXSamples.push_back(2.88757);
00165   theXSamples.push_back(3.06614);
00166   theXSamples.push_back(3.25193);
00167   theXSamples.push_back(3.44507);
00168   theXSamples.push_back(3.64571);
00169   theXSamples.push_back(3.854);
00170   theXSamples.push_back(4.07007);
00171   theXSamples.push_back(4.29407);
00172   theXSamples.push_back(4.52614);
00173   theXSamples.push_back(4.76643);
00174   theXSamples.push_back(5.01507);
00175   theXSamples.push_back(5.27221);
00176   theXSamples.push_back(5.538);
00177   theXSamples.push_back(5.81257);
00178   theXSamples.push_back(6.09607);
00179   theXSamples.push_back(6.38864);
00180   theXSamples.push_back(6.69043);
00181   theXSamples.push_back(7.00157);
00182   theXSamples.push_back(7.32221);
00183   theXSamples.push_back(7.6525);
00184   theXSamples.push_back(7.99257);
00185   theXSamples.push_back(8.34257);
00186   theXSamples.push_back(8.70264);
00187   theXSamples.push_back(9.07293);
00188   theXSamples.push_back(9.45357);
00189   theXSamples.push_back(9.84471);
00190   theXSamples.push_back(10.2465);
00191 
00192   for(int i=0; i<16; i++) 
00193     thePrehitSamples.push_back(i*0.02);
00194 
00195   vector<double> temp;
00196 
00197   temp.push_back(0);
00198   temp.push_back(0);
00199   temp.push_back(0);
00200   temp.push_back(8.40058e-06);
00201   temp.push_back(1.17956e-05);
00202   temp.push_back(2.16323e-05);
00203   temp.push_back(3.82308e-05);
00204   temp.push_back(6.09053e-05);
00205   temp.push_back(8.38503e-05);
00206   temp.push_back(0.000119927);
00207   temp.push_back(0.000162094);
00208   temp.push_back(0.000214873);
00209   temp.push_back(0.000263888);
00210   temp.push_back(0.000358384);
00211   temp.push_back(0.000405752);
00212   temp.push_back(0.0004813);
00213   temp.push_back(0.000578313);
00214   temp.push_back(0.000708949);
00215   temp.push_back(0.000806313);
00216   temp.push_back(0.000947174);
00217   temp.push_back(0.00103474);
00218   temp.push_back(0.00117471);
00219   temp.push_back(0.00128605);
00220   temp.push_back(0.00137934);
00221   temp.push_back(0.00154853);
00222   temp.push_back(0.00177649);
00223   temp.push_back(0.00183272);
00224   temp.push_back(0.00186753);
00225   temp.push_back(0.0021082);
00226   temp.push_back(0.00224273);
00227   temp.push_back(0.00230361);
00228   temp.push_back(0.00237872);
00229   temp.push_back(0.00245845);
00230   temp.push_back(0.00265224);
00231   temp.push_back(0.00271459);
00232   temp.push_back(0.00256564);
00233   temp.push_back(0.00266094);
00234   temp.push_back(0.00280213);
00235   temp.push_back(0.00267646);
00236   temp.push_back(0.00263196);
00237   temp.push_back(0.00256547);
00238   temp.push_back(0.00248559);
00239   temp.push_back(0.00254747);
00240   temp.push_back(0.00246776);
00241   temp.push_back(0.0024417);
00242   temp.push_back(0.00233673);
00243   temp.push_back(0.00221526);
00244   temp.push_back(0.00212567);
00245   temp.push_back(0.00179446);
00246   temp.push_back(0.00174539);
00247   temp.push_back(0.00163708);
00248   temp.push_back(0.00152897);
00249   temp.push_back(0.00137436);
00250   temp.push_back(0.0012264);
00251   temp.push_back(0.00115666);
00252   temp.push_back(0.00104865);
00253   temp.push_back(0.000931932);
00254   temp.push_back(0.000867489);
00255   temp.push_back(0.000773467);
00256   temp.push_back(0.000659164);
00257   temp.push_back(0.00060111);
00258   temp.push_back(0.00054868);
00259   temp.push_back(0.000450312);
00260   temp.push_back(0.000393679);
00261   temp.push_back(0.000341696);
00262   temp.push_back(0.000305829);
00263   temp.push_back(0.000238913);
00264   temp.push_back(0.000216915);
00265   temp.push_back(0.000181357);
00266   temp.push_back(0.000164413);
00267   temp.push_back(0.000133342);
00268   temp.push_back(0.000113031);
00269   temp.push_back(9.45459e-05);
00270   temp.push_back(7.48834e-05);
00271   temp.push_back(5.68479e-05);
00272   temp.push_back(4.7456e-05);
00273   theErrSamples.push_back(temp);
00274   temp.clear();
00275 
00276   temp.push_back(1.07019e-05);
00277   temp.push_back(1.4297e-05);
00278   temp.push_back(2.47542e-05);
00279   temp.push_back(3.98772e-05);
00280   temp.push_back(6.1143e-05);
00281   temp.push_back(7.55069e-05);
00282   temp.push_back(0.000101314);
00283   temp.push_back(0.000126813);
00284   temp.push_back(0.000169132);
00285   temp.push_back(0.00022141);
00286   temp.push_back(0.000261971);
00287   temp.push_back(0.000337033);
00288   temp.push_back(0.000353671);
00289   temp.push_back(0.000455462);
00290   temp.push_back(0.000529259);
00291   temp.push_back(0.00065141);
00292   temp.push_back(0.0007213);
00293   temp.push_back(0.000755735);
00294   temp.push_back(0.000905117);
00295   temp.push_back(0.000983389);
00296   temp.push_back(0.00117113);
00297   temp.push_back(0.00123377);
00298   temp.push_back(0.00140575);
00299   temp.push_back(0.0014624);
00300   temp.push_back(0.00158725);
00301   temp.push_back(0.0017951);
00302   temp.push_back(0.0019545);
00303   temp.push_back(0.00203318);
00304   temp.push_back(0.00228151);
00305   temp.push_back(0.00233437);
00306   temp.push_back(0.00236975);
00307   temp.push_back(0.00248279);
00308   temp.push_back(0.00257292);
00309   temp.push_back(0.00265839);
00310   temp.push_back(0.00291408);
00311   temp.push_back(0.00287581);
00312   temp.push_back(0.0027742);
00313   temp.push_back(0.0028192);
00314   temp.push_back(0.00290722);
00315   temp.push_back(0.00290954);
00316   temp.push_back(0.00289248);
00317   temp.push_back(0.0028919);
00318   temp.push_back(0.00274249);
00319   temp.push_back(0.00275067);
00320   temp.push_back(0.00254015);
00321   temp.push_back(0.00263388);
00322   temp.push_back(0.00246535);
00323   temp.push_back(0.00230695);
00324   temp.push_back(0.00219084);
00325   temp.push_back(0.00201822);
00326   temp.push_back(0.00195889);
00327   temp.push_back(0.00191163);
00328   temp.push_back(0.00195144);
00329   temp.push_back(0.0017406);
00330   temp.push_back(0.00171877);
00331   temp.push_back(0.00161755);
00332   temp.push_back(0.00161152);
00333   temp.push_back(0.00142289);
00334   temp.push_back(0.00142294);
00335   temp.push_back(0.00131059);
00336   temp.push_back(0.00131986);
00337   temp.push_back(0.00130397);
00338   temp.push_back(0.00131082);
00339   temp.push_back(0.00127312);
00340   temp.push_back(0.00127991);
00341   temp.push_back(0.00125962);
00342   temp.push_back(0.00123696);
00343   temp.push_back(0.00127488);
00344   temp.push_back(0.0012194);
00345   temp.push_back(0.00124351);
00346   temp.push_back(0.00130209);
00347   temp.push_back(0.00123718);
00348   temp.push_back(0.00126401);
00349   temp.push_back(0.0012762);
00350   temp.push_back(0.00118605);
00351   temp.push_back(0.00121322);
00352   theErrSamples.push_back(temp);
00353   temp.clear();
00354 
00355   temp.push_back(1.52753e-05);
00356   temp.push_back(1.9252e-05);
00357   temp.push_back(3.51239e-05);
00358   temp.push_back(5.42772e-05);
00359   temp.push_back(8.42695e-05);
00360   temp.push_back(0.000112201);
00361   temp.push_back(0.000147315);
00362   temp.push_back(0.000181674);
00363   temp.push_back(0.000224435);
00364   temp.push_back(0.000284745);
00365   temp.push_back(0.000327223);
00366   temp.push_back(0.000395938);
00367   temp.push_back(0.000472805);
00368   temp.push_back(0.000519986);
00369   temp.push_back(0.000635296);
00370   temp.push_back(0.000669197);
00371   temp.push_back(0.000775362);
00372   temp.push_back(0.000874534);
00373   temp.push_back(0.00102977);
00374   temp.push_back(0.00113617);
00375   temp.push_back(0.00121311);
00376   temp.push_back(0.00141158);
00377   temp.push_back(0.00154837);
00378   temp.push_back(0.00169857);
00379   temp.push_back(0.00188519);
00380   temp.push_back(0.00191032);
00381   temp.push_back(0.00206575);
00382   temp.push_back(0.00222722);
00383   temp.push_back(0.00227325);
00384   temp.push_back(0.00245166);
00385   temp.push_back(0.00259628);
00386   temp.push_back(0.00249795);
00387   temp.push_back(0.00272387);
00388   temp.push_back(0.00297299);
00389   temp.push_back(0.0030973);
00390   temp.push_back(0.00299177);
00391   temp.push_back(0.00288901);
00392   temp.push_back(0.00304009);
00393   temp.push_back(0.00300757);
00394   temp.push_back(0.003216);
00395   temp.push_back(0.00301877);
00396   temp.push_back(0.00305196);
00397   temp.push_back(0.00305521);
00398   temp.push_back(0.00264599);
00399   temp.push_back(0.00289825);
00400   temp.push_back(0.00273878);
00401   temp.push_back(0.0028623);
00402   temp.push_back(0.00262132);
00403   temp.push_back(0.00244983);
00404   temp.push_back(0.00238771);
00405   temp.push_back(0.00241464);
00406   temp.push_back(0.0022464);
00407   temp.push_back(0.00216141);
00408   temp.push_back(0.00216927);
00409   temp.push_back(0.00202911);
00410   temp.push_back(0.00201579);
00411   temp.push_back(0.00206435);
00412   temp.push_back(0.00189058);
00413   temp.push_back(0.0019032);
00414   temp.push_back(0.00180546);
00415   temp.push_back(0.00181397);
00416   temp.push_back(0.00177982);
00417   temp.push_back(0.00169705);
00418   temp.push_back(0.00175611);
00419   temp.push_back(0.00183977);
00420   temp.push_back(0.00181481);
00421   temp.push_back(0.0017491);
00422   temp.push_back(0.00167498);
00423   temp.push_back(0.00176056);
00424   temp.push_back(0.00168206);
00425   temp.push_back(0.00178997);
00426   temp.push_back(0.00160131);
00427   temp.push_back(0.00175127);
00428   temp.push_back(0.00175144);
00429   temp.push_back(0.00172983);
00430   temp.push_back(0.00176049);
00431   theErrSamples.push_back(temp);
00432   temp.clear();
00433 
00434   temp.push_back(1.56003e-05);
00435   temp.push_back(2.56961e-05);
00436   temp.push_back(3.98157e-05);
00437   temp.push_back(7.51062e-05);
00438   temp.push_back(0.000104808);
00439   temp.push_back(0.000137126);
00440   temp.push_back(0.000181945);
00441   temp.push_back(0.00023428);
00442   temp.push_back(0.000276685);
00443   temp.push_back(0.000340639);
00444   temp.push_back(0.000365219);
00445   temp.push_back(0.000455067);
00446   temp.push_back(0.000517277);
00447   temp.push_back(0.000653087);
00448   temp.push_back(0.000704082);
00449   temp.push_back(0.00081798);
00450   temp.push_back(0.000919737);
00451   temp.push_back(0.000997528);
00452   temp.push_back(0.00110852);
00453   temp.push_back(0.00127015);
00454   temp.push_back(0.00135952);
00455   temp.push_back(0.00157116);
00456   temp.push_back(0.00169714);
00457   temp.push_back(0.00181507);
00458   temp.push_back(0.00188281);
00459   temp.push_back(0.00187945);
00460   temp.push_back(0.00224784);
00461   temp.push_back(0.00234832);
00462   temp.push_back(0.00246759);
00463   temp.push_back(0.00258435);
00464   temp.push_back(0.00287211);
00465   temp.push_back(0.00282387);
00466   temp.push_back(0.00301207);
00467   temp.push_back(0.00301706);
00468   temp.push_back(0.00314772);
00469   temp.push_back(0.00313356);
00470   temp.push_back(0.00330207);
00471   temp.push_back(0.00321868);
00472   temp.push_back(0.00340633);
00473   temp.push_back(0.00316687);
00474   temp.push_back(0.0030367);
00475   temp.push_back(0.00346355);
00476   temp.push_back(0.00308358);
00477   temp.push_back(0.00313356);
00478   temp.push_back(0.00322158);
00479   temp.push_back(0.00286642);
00480   temp.push_back(0.0030063);
00481   temp.push_back(0.0028074);
00482   temp.push_back(0.00276659);
00483   temp.push_back(0.00284753);
00484   temp.push_back(0.00260215);
00485   temp.push_back(0.00256294);
00486   temp.push_back(0.00264208);
00487   temp.push_back(0.00247397);
00488   temp.push_back(0.00249761);
00489   temp.push_back(0.00232896);
00490   temp.push_back(0.00237722);
00491   temp.push_back(0.00223927);
00492   temp.push_back(0.00218478);
00493   temp.push_back(0.00219469);
00494   temp.push_back(0.00223608);
00495   temp.push_back(0.00225412);
00496   temp.push_back(0.00218023);
00497   temp.push_back(0.00223298);
00498   temp.push_back(0.0021489);
00499   temp.push_back(0.00216245);
00500   temp.push_back(0.00222762);
00501   temp.push_back(0.0021608);
00502   temp.push_back(0.00226409);
00503   temp.push_back(0.00215114);
00504   temp.push_back(0.00209354);
00505   temp.push_back(0.00220478);
00506   temp.push_back(0.00212817);
00507   temp.push_back(0.00216157);
00508   temp.push_back(0.00210899);
00509   temp.push_back(0.00221498);
00510   theErrSamples.push_back(temp);
00511   temp.clear();
00512 
00513   temp.push_back(1.98294e-05);
00514   temp.push_back(2.74883e-05);
00515   temp.push_back(4.97341e-05);
00516   temp.push_back(7.8167e-05);
00517   temp.push_back(0.000115158);
00518   temp.push_back(0.000159637);
00519   temp.push_back(0.000198413);
00520   temp.push_back(0.000249918);
00521   temp.push_back(0.000307694);
00522   temp.push_back(0.00037631);
00523   temp.push_back(0.000420037);
00524   temp.push_back(0.000530557);
00525   temp.push_back(0.000574406);
00526   temp.push_back(0.000678883);
00527   temp.push_back(0.000783538);
00528   temp.push_back(0.000828001);
00529   temp.push_back(0.000945383);
00530   temp.push_back(0.0010672);
00531   temp.push_back(0.00119183);
00532   temp.push_back(0.00134696);
00533   temp.push_back(0.00155348);
00534   temp.push_back(0.00167101);
00535   temp.push_back(0.00181066);
00536   temp.push_back(0.00187759);
00537   temp.push_back(0.00199899);
00538   temp.push_back(0.00225269);
00539   temp.push_back(0.0023458);
00540   temp.push_back(0.00253072);
00541   temp.push_back(0.00270998);
00542   temp.push_back(0.00263773);
00543   temp.push_back(0.00294804);
00544   temp.push_back(0.00294624);
00545   temp.push_back(0.00321668);
00546   temp.push_back(0.00332147);
00547   temp.push_back(0.00334766);
00548   temp.push_back(0.00330149);
00549   temp.push_back(0.00340882);
00550   temp.push_back(0.00347742);
00551   temp.push_back(0.00352589);
00552   temp.push_back(0.0035523);
00553   temp.push_back(0.0035411);
00554   temp.push_back(0.00340469);
00555   temp.push_back(0.00335693);
00556   temp.push_back(0.00345872);
00557   temp.push_back(0.00316091);
00558   temp.push_back(0.0031948);
00559   temp.push_back(0.00307664);
00560   temp.push_back(0.00320621);
00561   temp.push_back(0.00304138);
00562   temp.push_back(0.00299107);
00563   temp.push_back(0.00281006);
00564   temp.push_back(0.00278658);
00565   temp.push_back(0.0028095);
00566   temp.push_back(0.00268299);
00567   temp.push_back(0.0027215);
00568   temp.push_back(0.00264144);
00569   temp.push_back(0.00260714);
00570   temp.push_back(0.00255874);
00571   temp.push_back(0.00266137);
00572   temp.push_back(0.00248823);
00573   temp.push_back(0.0026194);
00574   temp.push_back(0.0026251);
00575   temp.push_back(0.00245108);
00576   temp.push_back(0.00254617);
00577   temp.push_back(0.00264651);
00578   temp.push_back(0.00247705);
00579   temp.push_back(0.00252758);
00580   temp.push_back(0.00262686);
00581   temp.push_back(0.00248693);
00582   temp.push_back(0.00254178);
00583   temp.push_back(0.00258846);
00584   temp.push_back(0.0025353);
00585   temp.push_back(0.00240568);
00586   temp.push_back(0.00250232);
00587   temp.push_back(0.00251452);
00588   temp.push_back(0.00258121);
00589   theErrSamples.push_back(temp);
00590   temp.clear();
00591 
00592   temp.push_back(2.362e-05);
00593   temp.push_back(3.24457e-05);
00594   temp.push_back(5.93625e-05);
00595   temp.push_back(9.71924e-05);
00596   temp.push_back(0.000126279);
00597   temp.push_back(0.000176152);
00598   temp.push_back(0.000224722);
00599   temp.push_back(0.000290262);
00600   temp.push_back(0.000351171);
00601   temp.push_back(0.000439513);
00602   temp.push_back(0.000497545);
00603   temp.push_back(0.000596431);
00604   temp.push_back(0.000686665);
00605   temp.push_back(0.000716486);
00606   temp.push_back(0.000835977);
00607   temp.push_back(0.000943229);
00608   temp.push_back(0.00104521);
00609   temp.push_back(0.00120528);
00610   temp.push_back(0.00140133);
00611   temp.push_back(0.00144132);
00612   temp.push_back(0.00164603);
00613   temp.push_back(0.00179182);
00614   temp.push_back(0.00188172);
00615   temp.push_back(0.00203418);
00616   temp.push_back(0.00224695);
00617   temp.push_back(0.00224972);
00618   temp.push_back(0.00259959);
00619   temp.push_back(0.00249426);
00620   temp.push_back(0.00280589);
00621   temp.push_back(0.00286528);
00622   temp.push_back(0.00309813);
00623   temp.push_back(0.00321952);
00624   temp.push_back(0.00324915);
00625   temp.push_back(0.00312759);
00626   temp.push_back(0.00338569);
00627   temp.push_back(0.00351915);
00628   temp.push_back(0.00336895);
00629   temp.push_back(0.00358675);
00630   temp.push_back(0.00361005);
00631   temp.push_back(0.00356783);
00632   temp.push_back(0.00381571);
00633   temp.push_back(0.00358167);
00634   temp.push_back(0.00355341);
00635   temp.push_back(0.00369961);
00636   temp.push_back(0.00345098);
00637   temp.push_back(0.00345749);
00638   temp.push_back(0.00366643);
00639   temp.push_back(0.00328805);
00640   temp.push_back(0.00317392);
00641   temp.push_back(0.00325794);
00642   temp.push_back(0.00325624);
00643   temp.push_back(0.00298589);
00644   temp.push_back(0.00313121);
00645   temp.push_back(0.00307129);
00646   temp.push_back(0.00303754);
00647   temp.push_back(0.00299004);
00648   temp.push_back(0.00283794);
00649   temp.push_back(0.0029663);
00650   temp.push_back(0.00278505);
00651   temp.push_back(0.00294565);
00652   temp.push_back(0.00295431);
00653   temp.push_back(0.0029553);
00654   temp.push_back(0.0027788);
00655   temp.push_back(0.00272655);
00656   temp.push_back(0.00267875);
00657   temp.push_back(0.00288197);
00658   temp.push_back(0.00284249);
00659   temp.push_back(0.00281758);
00660   temp.push_back(0.00285672);
00661   temp.push_back(0.00260368);
00662   temp.push_back(0.0027489);
00663   temp.push_back(0.00270647);
00664   temp.push_back(0.00300406);
00665   temp.push_back(0.00298123);
00666   temp.push_back(0.00284949);
00667   temp.push_back(0.00275038);
00668   theErrSamples.push_back(temp);
00669   temp.clear();
00670 
00671   temp.push_back(2.62119e-05);
00672   temp.push_back(4.06637e-05);
00673   temp.push_back(6.45721e-05);
00674   temp.push_back(9.8698e-05);
00675   temp.push_back(0.000145844);
00676   temp.push_back(0.000209191);
00677   temp.push_back(0.000267177);
00678   temp.push_back(0.000312714);
00679   temp.push_back(0.000397544);
00680   temp.push_back(0.000436342);
00681   temp.push_back(0.000537476);
00682   temp.push_back(0.000588913);
00683   temp.push_back(0.000733794);
00684   temp.push_back(0.000834061);
00685   temp.push_back(0.000884576);
00686   temp.push_back(0.001049);
00687   temp.push_back(0.00109329);
00688   temp.push_back(0.00128614);
00689   temp.push_back(0.00139434);
00690   temp.push_back(0.00155489);
00691   temp.push_back(0.00169101);
00692   temp.push_back(0.00192918);
00693   temp.push_back(0.00192531);
00694   temp.push_back(0.00211887);
00695   temp.push_back(0.00227517);
00696   temp.push_back(0.00250976);
00697   temp.push_back(0.00268224);
00698   temp.push_back(0.00276841);
00699   temp.push_back(0.00287698);
00700   temp.push_back(0.00299388);
00701   temp.push_back(0.00306973);
00702   temp.push_back(0.00337693);
00703   temp.push_back(0.00332126);
00704   temp.push_back(0.00344551);
00705   temp.push_back(0.00364214);
00706   temp.push_back(0.00354195);
00707   temp.push_back(0.00371903);
00708   temp.push_back(0.0037893);
00709   temp.push_back(0.00381634);
00710   temp.push_back(0.00377593);
00711   temp.push_back(0.00382319);
00712   temp.push_back(0.00378659);
00713   temp.push_back(0.00379056);
00714   temp.push_back(0.00392088);
00715   temp.push_back(0.00396608);
00716   temp.push_back(0.00361685);
00717   temp.push_back(0.00356004);
00718   temp.push_back(0.00377212);
00719   temp.push_back(0.00359799);
00720   temp.push_back(0.00349761);
00721   temp.push_back(0.00345083);
00722   temp.push_back(0.00341753);
00723   temp.push_back(0.0031536);
00724   temp.push_back(0.00329225);
00725   temp.push_back(0.00318691);
00726   temp.push_back(0.00342211);
00727   temp.push_back(0.00324309);
00728   temp.push_back(0.00343059);
00729   temp.push_back(0.00313786);
00730   temp.push_back(0.00309648);
00731   temp.push_back(0.00308665);
00732   temp.push_back(0.00315973);
00733   temp.push_back(0.00320034);
00734   temp.push_back(0.00327231);
00735   temp.push_back(0.00315402);
00736   temp.push_back(0.00306078);
00737   temp.push_back(0.00307093);
00738   temp.push_back(0.00304372);
00739   temp.push_back(0.0032217);
00740   temp.push_back(0.00315268);
00741   temp.push_back(0.0031173);
00742   temp.push_back(0.00308137);
00743   temp.push_back(0.0031778);
00744   temp.push_back(0.0030366);
00745   temp.push_back(0.00304068);
00746   temp.push_back(0.0033214);
00747   theErrSamples.push_back(temp);
00748   temp.clear();
00749 
00750   temp.push_back(2.74013e-05);
00751   temp.push_back(4.22899e-05);
00752   temp.push_back(6.78705e-05);
00753   temp.push_back(0.000108305);
00754   temp.push_back(0.000155136);
00755   temp.push_back(0.000215087);
00756   temp.push_back(0.000267491);
00757   temp.push_back(0.000324307);
00758   temp.push_back(0.000450428);
00759   temp.push_back(0.000525255);
00760   temp.push_back(0.000593983);
00761   temp.push_back(0.000670897);
00762   temp.push_back(0.000773512);
00763   temp.push_back(0.000886391);
00764   temp.push_back(0.000944693);
00765   temp.push_back(0.0012177);
00766   temp.push_back(0.00122543);
00767   temp.push_back(0.00137642);
00768   temp.push_back(0.00154086);
00769   temp.push_back(0.00165634);
00770   temp.push_back(0.00191029);
00771   temp.push_back(0.00190192);
00772   temp.push_back(0.00211654);
00773   temp.push_back(0.00212919);
00774   temp.push_back(0.00228748);
00775   temp.push_back(0.00254164);
00776   temp.push_back(0.00280712);
00777   temp.push_back(0.00302764);
00778   temp.push_back(0.00306686);
00779   temp.push_back(0.00327206);
00780   temp.push_back(0.00331136);
00781   temp.push_back(0.0036489);
00782   temp.push_back(0.00342041);
00783   temp.push_back(0.003781);
00784   temp.push_back(0.00378694);
00785   temp.push_back(0.00386973);
00786   temp.push_back(0.00369683);
00787   temp.push_back(0.0038999);
00788   temp.push_back(0.00395901);
00789   temp.push_back(0.00396228);
00790   temp.push_back(0.00405743);
00791   temp.push_back(0.00410244);
00792   temp.push_back(0.00413032);
00793   temp.push_back(0.00373123);
00794   temp.push_back(0.00395574);
00795   temp.push_back(0.00378686);
00796   temp.push_back(0.00406811);
00797   temp.push_back(0.00395806);
00798   temp.push_back(0.00379701);
00799   temp.push_back(0.00386635);
00800   temp.push_back(0.00361041);
00801   temp.push_back(0.00378335);
00802   temp.push_back(0.00353608);
00803   temp.push_back(0.00346058);
00804   temp.push_back(0.00361939);
00805   temp.push_back(0.00350638);
00806   temp.push_back(0.00351027);
00807   temp.push_back(0.00356889);
00808   temp.push_back(0.00353562);
00809   temp.push_back(0.00338339);
00810   temp.push_back(0.00358931);
00811   temp.push_back(0.0033404);
00812   temp.push_back(0.0034386);
00813   temp.push_back(0.00356428);
00814   temp.push_back(0.00342609);
00815   temp.push_back(0.00351769);
00816   temp.push_back(0.0034948);
00817   temp.push_back(0.0032928);
00818   temp.push_back(0.00354759);
00819   temp.push_back(0.0035006);
00820   temp.push_back(0.00344599);
00821   temp.push_back(0.00360147);
00822   temp.push_back(0.00350432);
00823   temp.push_back(0.00328072);
00824   temp.push_back(0.00352725);
00825   temp.push_back(0.00344304);
00826   theErrSamples.push_back(temp);
00827   temp.clear();
00828 
00829   temp.push_back(3.19736e-05);
00830   temp.push_back(4.2274e-05);
00831   temp.push_back(7.42979e-05);
00832   temp.push_back(0.000115308);
00833   temp.push_back(0.000181278);
00834   temp.push_back(0.000238536);
00835   temp.push_back(0.000295333);
00836   temp.push_back(0.000372787);
00837   temp.push_back(0.000450952);
00838   temp.push_back(0.000549762);
00839   temp.push_back(0.000587952);
00840   temp.push_back(0.000717376);
00841   temp.push_back(0.000823059);
00842   temp.push_back(0.000929103);
00843   temp.push_back(0.00101405);
00844   temp.push_back(0.00120059);
00845   temp.push_back(0.00129663);
00846   temp.push_back(0.00138426);
00847   temp.push_back(0.0016377);
00848   temp.push_back(0.00177567);
00849   temp.push_back(0.00193567);
00850   temp.push_back(0.00193153);
00851   temp.push_back(0.00217634);
00852   temp.push_back(0.00226176);
00853   temp.push_back(0.00269487);
00854   temp.push_back(0.00260071);
00855   temp.push_back(0.0029404);
00856   temp.push_back(0.00314849);
00857   temp.push_back(0.00329682);
00858   temp.push_back(0.00359669);
00859   temp.push_back(0.0033118);
00860   temp.push_back(0.00346258);
00861   temp.push_back(0.00376314);
00862   temp.push_back(0.00361881);
00863   temp.push_back(0.00386849);
00864   temp.push_back(0.00410369);
00865   temp.push_back(0.0039706);
00866   temp.push_back(0.0041219);
00867   temp.push_back(0.00407418);
00868   temp.push_back(0.00410784);
00869   temp.push_back(0.00417039);
00870   temp.push_back(0.00396376);
00871   temp.push_back(0.00402818);
00872   temp.push_back(0.00408969);
00873   temp.push_back(0.00398906);
00874   temp.push_back(0.00407158);
00875   temp.push_back(0.00407663);
00876   temp.push_back(0.00429465);
00877   temp.push_back(0.00399424);
00878   temp.push_back(0.00398511);
00879   temp.push_back(0.00413987);
00880   temp.push_back(0.00402974);
00881   temp.push_back(0.00396823);
00882   temp.push_back(0.00406268);
00883   temp.push_back(0.00397054);
00884   temp.push_back(0.00406304);
00885   temp.push_back(0.00381928);
00886   temp.push_back(0.00373153);
00887   temp.push_back(0.0038722);
00888   temp.push_back(0.00372542);
00889   temp.push_back(0.00371033);
00890   temp.push_back(0.00370851);
00891   temp.push_back(0.00370273);
00892   temp.push_back(0.00375404);
00893   temp.push_back(0.00368417);
00894   temp.push_back(0.00370411);
00895   temp.push_back(0.00360586);
00896   temp.push_back(0.00374387);
00897   temp.push_back(0.00370439);
00898   temp.push_back(0.00362516);
00899   temp.push_back(0.00377766);
00900   temp.push_back(0.00360445);
00901   temp.push_back(0.00360142);
00902   temp.push_back(0.00368594);
00903   temp.push_back(0.00358465);
00904   temp.push_back(0.00361813);
00905   theErrSamples.push_back(temp);
00906   temp.clear();
00907 
00908   temp.push_back(3.33533e-05);
00909   temp.push_back(4.71255e-05);
00910   temp.push_back(7.95104e-05);
00911   temp.push_back(0.000125354);
00912   temp.push_back(0.000184391);
00913   temp.push_back(0.000241585);
00914   temp.push_back(0.000312129);
00915   temp.push_back(0.000416606);
00916   temp.push_back(0.00048679);
00917   temp.push_back(0.000590734);
00918   temp.push_back(0.000680798);
00919   temp.push_back(0.00074592);
00920   temp.push_back(0.000925237);
00921   temp.push_back(0.000945055);
00922   temp.push_back(0.00111491);
00923   temp.push_back(0.00124566);
00924   temp.push_back(0.00130333);
00925   temp.push_back(0.00159285);
00926   temp.push_back(0.00160697);
00927   temp.push_back(0.00173291);
00928   temp.push_back(0.00205625);
00929   temp.push_back(0.0020842);
00930   temp.push_back(0.00234011);
00931   temp.push_back(0.00257964);
00932   temp.push_back(0.00271621);
00933   temp.push_back(0.00274754);
00934   temp.push_back(0.00300009);
00935   temp.push_back(0.00313262);
00936   temp.push_back(0.0033966);
00937   temp.push_back(0.00370825);
00938   temp.push_back(0.00365397);
00939   temp.push_back(0.00380956);
00940   temp.push_back(0.00406663);
00941   temp.push_back(0.00390965);
00942   temp.push_back(0.00400766);
00943   temp.push_back(0.0040838);
00944   temp.push_back(0.00430457);
00945   temp.push_back(0.00430605);
00946   temp.push_back(0.00452344);
00947   temp.push_back(0.00444883);
00948   temp.push_back(0.00462339);
00949   temp.push_back(0.00441296);
00950   temp.push_back(0.0044215);
00951   temp.push_back(0.00437965);
00952   temp.push_back(0.00450718);
00953   temp.push_back(0.00411557);
00954   temp.push_back(0.0044587);
00955   temp.push_back(0.00428602);
00956   temp.push_back(0.00415902);
00957   temp.push_back(0.00433486);
00958   temp.push_back(0.0043004);
00959   temp.push_back(0.00407849);
00960   temp.push_back(0.00423246);
00961   temp.push_back(0.00394595);
00962   temp.push_back(0.00402799);
00963   temp.push_back(0.00390578);
00964   temp.push_back(0.00403325);
00965   temp.push_back(0.00401427);
00966   temp.push_back(0.00413036);
00967   temp.push_back(0.00396507);
00968   temp.push_back(0.00397724);
00969   temp.push_back(0.00405608);
00970   temp.push_back(0.00410989);
00971   temp.push_back(0.00390445);
00972   temp.push_back(0.00408041);
00973   temp.push_back(0.00374423);
00974   temp.push_back(0.00380687);
00975   temp.push_back(0.00396012);
00976   temp.push_back(0.00403351);
00977   temp.push_back(0.00397022);
00978   temp.push_back(0.00420906);
00979   temp.push_back(0.00395898);
00980   temp.push_back(0.00414565);
00981   temp.push_back(0.00419401);
00982   temp.push_back(0.00398256);
00983   temp.push_back(0.00405866);
00984   theErrSamples.push_back(temp);
00985   temp.clear();
00986 
00987   temp.push_back(3.60154e-05);
00988   temp.push_back(4.78073e-05);
00989   temp.push_back(8.81868e-05);
00990   temp.push_back(0.00014119);
00991   temp.push_back(0.000199197);
00992   temp.push_back(0.000271447);
00993   temp.push_back(0.000348848);
00994   temp.push_back(0.000414794);
00995   temp.push_back(0.000517428);
00996   temp.push_back(0.000613861);
00997   temp.push_back(0.000706405);
00998   temp.push_back(0.000830507);
00999   temp.push_back(0.000954556);
01000   temp.push_back(0.00102745);
01001   temp.push_back(0.00115016);
01002   temp.push_back(0.00130855);
01003   temp.push_back(0.00141543);
01004   temp.push_back(0.00154704);
01005   temp.push_back(0.00183892);
01006   temp.push_back(0.00182328);
01007   temp.push_back(0.00211743);
01008   temp.push_back(0.00238445);
01009   temp.push_back(0.00251855);
01010   temp.push_back(0.00259965);
01011   temp.push_back(0.00276547);
01012   temp.push_back(0.00279785);
01013   temp.push_back(0.00310652);
01014   temp.push_back(0.00332247);
01015   temp.push_back(0.00348554);
01016   temp.push_back(0.00340047);
01017   temp.push_back(0.00379133);
01018   temp.push_back(0.00394976);
01019   temp.push_back(0.00411556);
01020   temp.push_back(0.00395808);
01021   temp.push_back(0.00413605);
01022   temp.push_back(0.00434804);
01023   temp.push_back(0.00454974);
01024   temp.push_back(0.00448955);
01025   temp.push_back(0.00475876);
01026   temp.push_back(0.00457681);
01027   temp.push_back(0.00444923);
01028   temp.push_back(0.00441159);
01029   temp.push_back(0.00433255);
01030   temp.push_back(0.0044053);
01031   temp.push_back(0.00467223);
01032   temp.push_back(0.00444452);
01033   temp.push_back(0.00429894);
01034   temp.push_back(0.00466592);
01035   temp.push_back(0.004516);
01036   temp.push_back(0.00446295);
01037   temp.push_back(0.00435999);
01038   temp.push_back(0.00467332);
01039   temp.push_back(0.00423821);
01040   temp.push_back(0.00443203);
01041   temp.push_back(0.00434006);
01042   temp.push_back(0.00429902);
01043   temp.push_back(0.0042231);
01044   temp.push_back(0.00423514);
01045   temp.push_back(0.00451164);
01046   temp.push_back(0.00436348);
01047   temp.push_back(0.00436551);
01048   temp.push_back(0.00425205);
01049   temp.push_back(0.00429035);
01050   temp.push_back(0.00411415);
01051   temp.push_back(0.00427361);
01052   temp.push_back(0.0039277);
01053   temp.push_back(0.00405848);
01054   temp.push_back(0.00429804);
01055   temp.push_back(0.00408766);
01056   temp.push_back(0.00384011);
01057   temp.push_back(0.00435199);
01058   temp.push_back(0.004093);
01059   temp.push_back(0.00406005);
01060   temp.push_back(0.00431062);
01061   temp.push_back(0.00411773);
01062   temp.push_back(0.00435653);
01063   theErrSamples.push_back(temp);
01064   temp.clear();
01065 
01066   temp.push_back(3.94789e-05);
01067   temp.push_back(5.3769e-05);
01068   temp.push_back(9.3841e-05);
01069   temp.push_back(0.000149774);
01070   temp.push_back(0.000215262);
01071   temp.push_back(0.000267695);
01072   temp.push_back(0.000370277);
01073   temp.push_back(0.000473963);
01074   temp.push_back(0.000548657);
01075   temp.push_back(0.000625464);
01076   temp.push_back(0.000737716);
01077   temp.push_back(0.000895534);
01078   temp.push_back(0.00100207);
01079   temp.push_back(0.00111602);
01080   temp.push_back(0.00123458);
01081   temp.push_back(0.00138359);
01082   temp.push_back(0.00158882);
01083   temp.push_back(0.00176246);
01084   temp.push_back(0.00177764);
01085   temp.push_back(0.00199677);
01086   temp.push_back(0.00211274);
01087   temp.push_back(0.00240314);
01088   temp.push_back(0.00238892);
01089   temp.push_back(0.00266706);
01090   temp.push_back(0.00292668);
01091   temp.push_back(0.00314009);
01092   temp.push_back(0.00330102);
01093   temp.push_back(0.00354671);
01094   temp.push_back(0.00371555);
01095   temp.push_back(0.00383521);
01096   temp.push_back(0.00392384);
01097   temp.push_back(0.00415042);
01098   temp.push_back(0.00417837);
01099   temp.push_back(0.00410063);
01100   temp.push_back(0.00436123);
01101   temp.push_back(0.00453869);
01102   temp.push_back(0.00456717);
01103   temp.push_back(0.00432658);
01104   temp.push_back(0.00474763);
01105   temp.push_back(0.00479594);
01106   temp.push_back(0.00495067);
01107   temp.push_back(0.00486827);
01108   temp.push_back(0.00501981);
01109   temp.push_back(0.00467952);
01110   temp.push_back(0.00496726);
01111   temp.push_back(0.00471759);
01112   temp.push_back(0.00487664);
01113   temp.push_back(0.00477905);
01114   temp.push_back(0.00457876);
01115   temp.push_back(0.00478556);
01116   temp.push_back(0.00470292);
01117   temp.push_back(0.00474277);
01118   temp.push_back(0.00474074);
01119   temp.push_back(0.00476266);
01120   temp.push_back(0.00459758);
01121   temp.push_back(0.00465773);
01122   temp.push_back(0.00469381);
01123   temp.push_back(0.004336);
01124   temp.push_back(0.00439833);
01125   temp.push_back(0.00434989);
01126   temp.push_back(0.00452046);
01127   temp.push_back(0.00444761);
01128   temp.push_back(0.00452204);
01129   temp.push_back(0.00433173);
01130   temp.push_back(0.00428761);
01131   temp.push_back(0.00454501);
01132   temp.push_back(0.00466898);
01133   temp.push_back(0.00449207);
01134   temp.push_back(0.00478147);
01135   temp.push_back(0.00448187);
01136   temp.push_back(0.00427879);
01137   temp.push_back(0.0045072);
01138   temp.push_back(0.00446825);
01139   temp.push_back(0.00442717);
01140   temp.push_back(0.00450486);
01141   temp.push_back(0.004514);
01142   theErrSamples.push_back(temp);
01143   temp.clear();
01144 
01145   temp.push_back(4.07374e-05);
01146   temp.push_back(5.84427e-05);
01147   temp.push_back(0.000102379);
01148   temp.push_back(0.000154485);
01149   temp.push_back(0.00022167);
01150   temp.push_back(0.000291138);
01151   temp.push_back(0.000366195);
01152   temp.push_back(0.00049284);
01153   temp.push_back(0.000555038);
01154   temp.push_back(0.000718216);
01155   temp.push_back(0.000815782);
01156   temp.push_back(0.000929635);
01157   temp.push_back(0.00105284);
01158   temp.push_back(0.00117775);
01159   temp.push_back(0.00128659);
01160   temp.push_back(0.0014848);
01161   temp.push_back(0.00161139);
01162   temp.push_back(0.00186044);
01163   temp.push_back(0.00193622);
01164   temp.push_back(0.0021136);
01165   temp.push_back(0.00218413);
01166   temp.push_back(0.00254649);
01167   temp.push_back(0.00257878);
01168   temp.push_back(0.00286384);
01169   temp.push_back(0.00277422);
01170   temp.push_back(0.00337728);
01171   temp.push_back(0.0034849);
01172   temp.push_back(0.00340672);
01173   temp.push_back(0.00376872);
01174   temp.push_back(0.00397498);
01175   temp.push_back(0.00407663);
01176   temp.push_back(0.00416444);
01177   temp.push_back(0.00447438);
01178   temp.push_back(0.00421777);
01179   temp.push_back(0.00434779);
01180   temp.push_back(0.00457815);
01181   temp.push_back(0.00510487);
01182   temp.push_back(0.004957);
01183   temp.push_back(0.0050608);
01184   temp.push_back(0.00482047);
01185   temp.push_back(0.00479584);
01186   temp.push_back(0.00506816);
01187   temp.push_back(0.00485049);
01188   temp.push_back(0.00502603);
01189   temp.push_back(0.00496904);
01190   temp.push_back(0.0048041);
01191   temp.push_back(0.00498567);
01192   temp.push_back(0.00505111);
01193   temp.push_back(0.00495093);
01194   temp.push_back(0.00493909);
01195   temp.push_back(0.00518096);
01196   temp.push_back(0.00475333);
01197   temp.push_back(0.00509383);
01198   temp.push_back(0.00455324);
01199   temp.push_back(0.00501535);
01200   temp.push_back(0.00477645);
01201   temp.push_back(0.00467969);
01202   temp.push_back(0.00482058);
01203   temp.push_back(0.0049352);
01204   temp.push_back(0.00479813);
01205   temp.push_back(0.00456841);
01206   temp.push_back(0.00495553);
01207   temp.push_back(0.00475263);
01208   temp.push_back(0.00464536);
01209   temp.push_back(0.0045788);
01210   temp.push_back(0.00482641);
01211   temp.push_back(0.00468534);
01212   temp.push_back(0.00458609);
01213   temp.push_back(0.00469523);
01214   temp.push_back(0.00471367);
01215   temp.push_back(0.00480113);
01216   temp.push_back(0.00477874);
01217   temp.push_back(0.00487724);
01218   temp.push_back(0.00448248);
01219   temp.push_back(0.00457135);
01220   temp.push_back(0.00463623);
01221   theErrSamples.push_back(temp);
01222   temp.clear();
01223 
01224   temp.push_back(4.19494e-05);
01225   temp.push_back(5.81133e-05);
01226   temp.push_back(0.000106598);
01227   temp.push_back(0.000167045);
01228   temp.push_back(0.000245601);
01229   temp.push_back(0.000309417);
01230   temp.push_back(0.000402556);
01231   temp.push_back(0.000531311);
01232   temp.push_back(0.000599414);
01233   temp.push_back(0.000724794);
01234   temp.push_back(0.000843293);
01235   temp.push_back(0.00102089);
01236   temp.push_back(0.00110244);
01237   temp.push_back(0.00119442);
01238   temp.push_back(0.0014845);
01239   temp.push_back(0.00144809);
01240   temp.push_back(0.00163904);
01241   temp.push_back(0.00179288);
01242   temp.push_back(0.00208349);
01243   temp.push_back(0.00218156);
01244   temp.push_back(0.00241731);
01245   temp.push_back(0.00260733);
01246   temp.push_back(0.0028223);
01247   temp.push_back(0.00303541);
01248   temp.push_back(0.00317796);
01249   temp.push_back(0.00318073);
01250   temp.push_back(0.00354012);
01251   temp.push_back(0.0038304);
01252   temp.push_back(0.00394061);
01253   temp.push_back(0.00388992);
01254   temp.push_back(0.00415685);
01255   temp.push_back(0.004271);
01256   temp.push_back(0.00454353);
01257   temp.push_back(0.00436867);
01258   temp.push_back(0.00471855);
01259   temp.push_back(0.00486456);
01260   temp.push_back(0.00458562);
01261   temp.push_back(0.00467012);
01262   temp.push_back(0.00488229);
01263   temp.push_back(0.0051219);
01264   temp.push_back(0.00525306);
01265   temp.push_back(0.00520996);
01266   temp.push_back(0.00512736);
01267   temp.push_back(0.00506308);
01268   temp.push_back(0.00509758);
01269   temp.push_back(0.00513977);
01270   temp.push_back(0.00506956);
01271   temp.push_back(0.0050327);
01272   temp.push_back(0.00528484);
01273   temp.push_back(0.00541921);
01274   temp.push_back(0.00486394);
01275   temp.push_back(0.00532431);
01276   temp.push_back(0.00491166);
01277   temp.push_back(0.0054032);
01278   temp.push_back(0.00515276);
01279   temp.push_back(0.00499206);
01280   temp.push_back(0.00507303);
01281   temp.push_back(0.00528015);
01282   temp.push_back(0.00502935);
01283   temp.push_back(0.00498204);
01284   temp.push_back(0.00482074);
01285   temp.push_back(0.0050209);
01286   temp.push_back(0.00506188);
01287   temp.push_back(0.00501215);
01288   temp.push_back(0.00523337);
01289   temp.push_back(0.00501484);
01290   temp.push_back(0.00482134);
01291   temp.push_back(0.00510347);
01292   temp.push_back(0.00500609);
01293   temp.push_back(0.00508656);
01294   temp.push_back(0.00505679);
01295   temp.push_back(0.00526603);
01296   temp.push_back(0.00491982);
01297   temp.push_back(0.00505811);
01298   temp.push_back(0.00494347);
01299   temp.push_back(0.00506382);
01300   theErrSamples.push_back(temp);
01301   temp.clear();
01302 
01303   temp.push_back(4.58461e-05);
01304   temp.push_back(6.56593e-05);
01305   temp.push_back(0.000106249);
01306   temp.push_back(0.000177259);
01307   temp.push_back(0.000242624);
01308   temp.push_back(0.000318801);
01309   temp.push_back(0.000429442);
01310   temp.push_back(0.000524763);
01311   temp.push_back(0.000627853);
01312   temp.push_back(0.000767861);
01313   temp.push_back(0.000914838);
01314   temp.push_back(0.00101068);
01315   temp.push_back(0.00119111);
01316   temp.push_back(0.00124833);
01317   temp.push_back(0.00147471);
01318   temp.push_back(0.00151766);
01319   temp.push_back(0.00175504);
01320   temp.push_back(0.0019513);
01321   temp.push_back(0.00225088);
01322   temp.push_back(0.00232528);
01323   temp.push_back(0.002477);
01324   temp.push_back(0.0026629);
01325   temp.push_back(0.00282126);
01326   temp.push_back(0.00305865);
01327   temp.push_back(0.00328165);
01328   temp.push_back(0.00342969);
01329   temp.push_back(0.00368517);
01330   temp.push_back(0.00412023);
01331   temp.push_back(0.00408107);
01332   temp.push_back(0.00442514);
01333   temp.push_back(0.00419995);
01334   temp.push_back(0.00445839);
01335   temp.push_back(0.00480241);
01336   temp.push_back(0.00469063);
01337   temp.push_back(0.00488558);
01338   temp.push_back(0.00505738);
01339   temp.push_back(0.00534401);
01340   temp.push_back(0.00529189);
01341   temp.push_back(0.00503914);
01342   temp.push_back(0.00547271);
01343   temp.push_back(0.00532523);
01344   temp.push_back(0.0056474);
01345   temp.push_back(0.00528356);
01346   temp.push_back(0.00537977);
01347   temp.push_back(0.00528494);
01348   temp.push_back(0.00555228);
01349   temp.push_back(0.00522019);
01350   temp.push_back(0.00551423);
01351   temp.push_back(0.00532466);
01352   temp.push_back(0.00562846);
01353   temp.push_back(0.00531636);
01354   temp.push_back(0.00534452);
01355   temp.push_back(0.00546899);
01356   temp.push_back(0.00540307);
01357   temp.push_back(0.00522696);
01358   temp.push_back(0.00543343);
01359   temp.push_back(0.00496315);
01360   temp.push_back(0.0054688);
01361   temp.push_back(0.0054534);
01362   temp.push_back(0.00529457);
01363   temp.push_back(0.00527906);
01364   temp.push_back(0.00505608);
01365   temp.push_back(0.00532934);
01366   temp.push_back(0.00540523);
01367   temp.push_back(0.00546589);
01368   temp.push_back(0.00541721);
01369   temp.push_back(0.00518411);
01370   temp.push_back(0.0050658);
01371   temp.push_back(0.00538019);
01372   temp.push_back(0.00523502);
01373   temp.push_back(0.00518119);
01374   temp.push_back(0.00541954);
01375   temp.push_back(0.00556926);
01376   temp.push_back(0.00516485);
01377   temp.push_back(0.00543988);
01378   temp.push_back(0.00534757);
01379   theErrSamples.push_back(temp);
01380   temp.clear();
01381 
01382   temp.push_back(4.87191e-05);
01383   temp.push_back(6.79557e-05);
01384   temp.push_back(0.00011585);
01385   temp.push_back(0.000170012);
01386   temp.push_back(0.000249995);
01387   temp.push_back(0.000358082);
01388   temp.push_back(0.000441529);
01389   temp.push_back(0.000573455);
01390   temp.push_back(0.000686492);
01391   temp.push_back(0.000768395);
01392   temp.push_back(0.000922129);
01393   temp.push_back(0.0010261);
01394   temp.push_back(0.00121379);
01395   temp.push_back(0.00132394);
01396   temp.push_back(0.00147006);
01397   temp.push_back(0.0016026);
01398   temp.push_back(0.00185945);
01399   temp.push_back(0.00205459);
01400   temp.push_back(0.00224833);
01401   temp.push_back(0.00239438);
01402   temp.push_back(0.00258948);
01403   temp.push_back(0.00285399);
01404   temp.push_back(0.00293392);
01405   temp.push_back(0.00310816);
01406   temp.push_back(0.0035526);
01407   temp.push_back(0.00353642);
01408   temp.push_back(0.003711);
01409   temp.push_back(0.00392946);
01410   temp.push_back(0.00417796);
01411   temp.push_back(0.00440017);
01412   temp.push_back(0.00474434);
01413   temp.push_back(0.00469919);
01414   temp.push_back(0.00476225);
01415   temp.push_back(0.00474708);
01416   temp.push_back(0.00515252);
01417   temp.push_back(0.0053985);
01418   temp.push_back(0.00515269);
01419   temp.push_back(0.00535194);
01420   temp.push_back(0.00519113);
01421   temp.push_back(0.00531732);
01422   temp.push_back(0.00545661);
01423   temp.push_back(0.00562353);
01424   temp.push_back(0.00547828);
01425   temp.push_back(0.00590485);
01426   temp.push_back(0.0054847);
01427   temp.push_back(0.00575722);
01428   temp.push_back(0.00578941);
01429   temp.push_back(0.0055163);
01430   temp.push_back(0.00565193);
01431   temp.push_back(0.00559088);
01432   temp.push_back(0.00578905);
01433   temp.push_back(0.00568463);
01434   temp.push_back(0.00538114);
01435   temp.push_back(0.00569508);
01436   temp.push_back(0.00534805);
01437   temp.push_back(0.00521921);
01438   temp.push_back(0.00568289);
01439   temp.push_back(0.0055866);
01440   temp.push_back(0.00561974);
01441   temp.push_back(0.00548487);
01442   temp.push_back(0.00558498);
01443   temp.push_back(0.00536631);
01444   temp.push_back(0.00541084);
01445   temp.push_back(0.00577058);
01446   temp.push_back(0.00540762);
01447   temp.push_back(0.00545929);
01448   temp.push_back(0.00571948);
01449   temp.push_back(0.00546965);
01450   temp.push_back(0.00557277);
01451   temp.push_back(0.00554671);
01452   temp.push_back(0.00521175);
01453   temp.push_back(0.00546106);
01454   temp.push_back(0.00541317);
01455   temp.push_back(0.00569926);
01456   temp.push_back(0.0055058);
01457   temp.push_back(0.00537131);
01458   theErrSamples.push_back(temp);
01459   temp.clear();
01460 
01461 }
/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/SimCalorimetry/HcalSimAlgos/src/HcalSiPM.cc
