/data/refman/pasoursint/CMSSW_5_3_10_patch2/src/SimCalorimetry/HcalSimAlgos/src/HcalSiPM.cc

Go to the documentation of this file.

#include "SimCalorimetry/HcalSimAlgos/interface/HcalSiPM.h"
#include "FWCore/Utilities/interface/RandomNumberGenerator.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include <cmath>

using std::vector;

HcalSiPM::HcalSiPM(int nCells) :
  theCellCount(nCells), theRndGauss(0) {

  assert(theCellCount>0);

  defaultErrInit();
}

HcalSiPM::~HcalSiPM() {
  delete theRndGauss;
}

int HcalSiPM::hitCells(int photons, int integral) const {
  //don't need to do zero or negative photons.
  if (photons < 1) return 0;
  if (integral < 0) integral = 0;
  if (integral >= theCellCount) return 0;

  if (theRndGauss == 0) {
    //random number generator setup
    edm::Service<edm::RandomNumberGenerator> rng;
    if ( ! rng.isAvailable()) {
      throw cms::Exception("Configuration")
        << "HcalSiPM requires the RandomNumberGeneratorService\n"
        "which is not present in the configuration file.  "
        "You must add the service\n"
        "in the configuration file or remove the modules that require it.";
    }

    CLHEP::HepRandomEngine& engine = rng->getEngine();
    theRndGauss = new CLHEP::RandGaussQ(engine);
  }

  //normalize by theCellCount to remove dependency on SiPM size and pixel density.
  double x = double(photons)/double(theCellCount);
  double prehit = double(integral)/double(theCellCount);

  //calculate the width and mean of the distribution for a given x
  double width = errOnX(x, prehit);
  double mean = 1. - std::exp(-x);

  //you can't hit more than everything.
  if (mean > 1.) mean = 1.;

  //convert back to absolute pixels
  mean *= (1-prehit)*theCellCount;
  width *= (1-prehit)*theCellCount;

  double npe;
  while (true) {
    npe = theRndGauss->fire(mean, width);
    if ((npe > -0.5) && (npe <= theCellCount-integral) && (npe <= photons))
      return int(npe + 0.5);
  }
}

void HcalSiPM::setNCells(int nCells) {
  assert(nCells>0);
  theCellCount = nCells;
}

void HcalSiPM::initRandomEngine(CLHEP::HepRandomEngine& engine) {
  if(theRndGauss) delete theRndGauss;
  theRndGauss = new CLHEP::RandGaussQ(engine);
}

double HcalSiPM::errOnX(double x, double prehit) const {
  //interpolate between the points in the x,prehit space to determine the
  //width of the distribution of x
  int before, after;
  getBeforeAndAfter(x, before, after, theXSamples);
  int i, j;
  getBeforeAndAfter(prehit, i, j, thePrehitSamples);
  double m1 = (theErrSamples[i][after]-theErrSamples[i][before])/
    (theXSamples[after]-theXSamples[before]);
  double err1 = theErrSamples[i][before]+m1*(x-theXSamples[before]);
  double m2 = (theErrSamples[j][after]-theErrSamples[j][before])/
    (theXSamples[after]-theXSamples[before]);
  double err2 = theErrSamples[j][before]+m2*(x-theXSamples[before]);

  return err1 + (err2-err1)/(thePrehitSamples[j]-thePrehitSamples[i])*
    (prehit - thePrehitSamples[i]);
}

void HcalSiPM::getBeforeAndAfter(double val, int& before, int& after, 
                                 const vector<double>& vec) const {
  //binary like search to indeces points around an arbitrary value.
  before = 0;
  after = vec.size()-1;
  int mid = vec.size()/2;
  bool go = true;
  do {
    if (vec[mid]>val) after = mid;
    else before = mid;
    mid = (after+before)/2;
    go = (after-before > 1);
  } while (go);
}

void HcalSiPM::defaultErrInit() {
  //load up values for the interpolation for the errs on x used to generate
  //the number of photons.  These were produced from a toy SiPM model.
  theXSamples.clear();
  thePrehitSamples.clear();
  theErrSamples.clear();

  theXSamples.push_back(7.14286e-05);
  theXSamples.push_back(0.000142857);
  theXSamples.push_back(0.000428571);
  theXSamples.push_back(0.00107143);
  theXSamples.push_back(0.00221429);
  theXSamples.push_back(0.004);
  theXSamples.push_back(0.00657143);
  theXSamples.push_back(0.0100714);
  theXSamples.push_back(0.0146429);
  theXSamples.push_back(0.0204286);
  theXSamples.push_back(0.0275714);
  theXSamples.push_back(0.0362143);
  theXSamples.push_back(0.0465);
  theXSamples.push_back(0.0585714);
  theXSamples.push_back(0.0725714);
  theXSamples.push_back(0.0886429);
  theXSamples.push_back(0.106929);
  theXSamples.push_back(0.127571);
  theXSamples.push_back(0.150714);
  theXSamples.push_back(0.1765);
  theXSamples.push_back(0.205071);
  theXSamples.push_back(0.236571);
  theXSamples.push_back(0.271143);
  theXSamples.push_back(0.308929);
  theXSamples.push_back(0.350071);
  theXSamples.push_back(0.394714);
  theXSamples.push_back(0.443);
  theXSamples.push_back(0.495071);
  theXSamples.push_back(0.551071);
  theXSamples.push_back(0.611143);
  theXSamples.push_back(0.675429);
  theXSamples.push_back(0.744071);
  theXSamples.push_back(0.817214);
  theXSamples.push_back(0.895);
  theXSamples.push_back(0.977571);
  theXSamples.push_back(1.06507);
  theXSamples.push_back(1.15764);
  theXSamples.push_back(1.25543);
  theXSamples.push_back(1.35857);
  theXSamples.push_back(1.46721);
  theXSamples.push_back(1.5815);
  theXSamples.push_back(1.70157);
  theXSamples.push_back(1.82757);
  theXSamples.push_back(1.95964);
  theXSamples.push_back(2.09793);
  theXSamples.push_back(2.24257);
  theXSamples.push_back(2.39371);
  theXSamples.push_back(2.5515);
  theXSamples.push_back(2.71607);
  theXSamples.push_back(2.88757);
  theXSamples.push_back(3.06614);
  theXSamples.push_back(3.25193);
  theXSamples.push_back(3.44507);
  theXSamples.push_back(3.64571);
  theXSamples.push_back(3.854);
  theXSamples.push_back(4.07007);
  theXSamples.push_back(4.29407);
  theXSamples.push_back(4.52614);
  theXSamples.push_back(4.76643);
  theXSamples.push_back(5.01507);
  theXSamples.push_back(5.27221);
  theXSamples.push_back(5.538);
  theXSamples.push_back(5.81257);
  theXSamples.push_back(6.09607);
  theXSamples.push_back(6.38864);
  theXSamples.push_back(6.69043);
  theXSamples.push_back(7.00157);
  theXSamples.push_back(7.32221);
  theXSamples.push_back(7.6525);
  theXSamples.push_back(7.99257);
  theXSamples.push_back(8.34257);
  theXSamples.push_back(8.70264);
  theXSamples.push_back(9.07293);
  theXSamples.push_back(9.45357);
  theXSamples.push_back(9.84471);
  theXSamples.push_back(10.2465);

  for(int i=0; i<16; i++) 
    thePrehitSamples.push_back(i*0.02);

  vector<double> temp;

  temp.push_back(0);
  temp.push_back(0);
  temp.push_back(0);
  temp.push_back(8.40058e-06);
  temp.push_back(1.17956e-05);
  temp.push_back(2.16323e-05);
  temp.push_back(3.82308e-05);
  temp.push_back(6.09053e-05);
  temp.push_back(8.38503e-05);
  temp.push_back(0.000119927);
  temp.push_back(0.000162094);
  temp.push_back(0.000214873);
  temp.push_back(0.000263888);
  temp.push_back(0.000358384);
  temp.push_back(0.000405752);
  temp.push_back(0.0004813);
  temp.push_back(0.000578313);
  temp.push_back(0.000708949);
  temp.push_back(0.000806313);
  temp.push_back(0.000947174);
  temp.push_back(0.00103474);
  temp.push_back(0.00117471);
  temp.push_back(0.00128605);
  temp.push_back(0.00137934);
  temp.push_back(0.00154853);
  temp.push_back(0.00177649);
  temp.push_back(0.00183272);
  temp.push_back(0.00186753);
  temp.push_back(0.0021082);
  temp.push_back(0.00224273);
  temp.push_back(0.00230361);
  temp.push_back(0.00237872);
  temp.push_back(0.00245845);
  temp.push_back(0.00265224);
  temp.push_back(0.00271459);
  temp.push_back(0.00256564);
  temp.push_back(0.00266094);
  temp.push_back(0.00280213);
  temp.push_back(0.00267646);
  temp.push_back(0.00263196);
  temp.push_back(0.00256547);
  temp.push_back(0.00248559);
  temp.push_back(0.00254747);
  temp.push_back(0.00246776);
  temp.push_back(0.0024417);
  temp.push_back(0.00233673);
  temp.push_back(0.00221526);
  temp.push_back(0.00212567);
  temp.push_back(0.00179446);
  temp.push_back(0.00174539);
  temp.push_back(0.00163708);
  temp.push_back(0.00152897);
  temp.push_back(0.00137436);
  temp.push_back(0.0012264);
  temp.push_back(0.00115666);
  temp.push_back(0.00104865);
  temp.push_back(0.000931932);
  temp.push_back(0.000867489);
  temp.push_back(0.000773467);
  temp.push_back(0.000659164);
  temp.push_back(0.00060111);
  temp.push_back(0.00054868);
  temp.push_back(0.000450312);
  temp.push_back(0.000393679);
  temp.push_back(0.000341696);
  temp.push_back(0.000305829);
  temp.push_back(0.000238913);
  temp.push_back(0.000216915);
  temp.push_back(0.000181357);
  temp.push_back(0.000164413);
  temp.push_back(0.000133342);
  temp.push_back(0.000113031);
  temp.push_back(9.45459e-05);
  temp.push_back(7.48834e-05);
  temp.push_back(5.68479e-05);
  temp.push_back(4.7456e-05);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(1.07019e-05);
  temp.push_back(1.4297e-05);
  temp.push_back(2.47542e-05);
  temp.push_back(3.98772e-05);
  temp.push_back(6.1143e-05);
  temp.push_back(7.55069e-05);
  temp.push_back(0.000101314);
  temp.push_back(0.000126813);
  temp.push_back(0.000169132);
  temp.push_back(0.00022141);
  temp.push_back(0.000261971);
  temp.push_back(0.000337033);
  temp.push_back(0.000353671);
  temp.push_back(0.000455462);
  temp.push_back(0.000529259);
  temp.push_back(0.00065141);
  temp.push_back(0.0007213);
  temp.push_back(0.000755735);
  temp.push_back(0.000905117);
  temp.push_back(0.000983389);
  temp.push_back(0.00117113);
  temp.push_back(0.00123377);
  temp.push_back(0.00140575);
  temp.push_back(0.0014624);
  temp.push_back(0.00158725);
  temp.push_back(0.0017951);
  temp.push_back(0.0019545);
  temp.push_back(0.00203318);
  temp.push_back(0.00228151);
  temp.push_back(0.00233437);
  temp.push_back(0.00236975);
  temp.push_back(0.00248279);
  temp.push_back(0.00257292);
  temp.push_back(0.00265839);
  temp.push_back(0.00291408);
  temp.push_back(0.00287581);
  temp.push_back(0.0027742);
  temp.push_back(0.0028192);
  temp.push_back(0.00290722);
  temp.push_back(0.00290954);
  temp.push_back(0.00289248);
  temp.push_back(0.0028919);
  temp.push_back(0.00274249);
  temp.push_back(0.00275067);
  temp.push_back(0.00254015);
  temp.push_back(0.00263388);
  temp.push_back(0.00246535);
  temp.push_back(0.00230695);
  temp.push_back(0.00219084);
  temp.push_back(0.00201822);
  temp.push_back(0.00195889);
  temp.push_back(0.00191163);
  temp.push_back(0.00195144);
  temp.push_back(0.0017406);
  temp.push_back(0.00171877);
  temp.push_back(0.00161755);
  temp.push_back(0.00161152);
  temp.push_back(0.00142289);
  temp.push_back(0.00142294);
  temp.push_back(0.00131059);
  temp.push_back(0.00131986);
  temp.push_back(0.00130397);
  temp.push_back(0.00131082);
  temp.push_back(0.00127312);
  temp.push_back(0.00127991);
  temp.push_back(0.00125962);
  temp.push_back(0.00123696);
  temp.push_back(0.00127488);
  temp.push_back(0.0012194);
  temp.push_back(0.00124351);
  temp.push_back(0.00130209);
  temp.push_back(0.00123718);
  temp.push_back(0.00126401);
  temp.push_back(0.0012762);
  temp.push_back(0.00118605);
  temp.push_back(0.00121322);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(1.52753e-05);
  temp.push_back(1.9252e-05);
  temp.push_back(3.51239e-05);
  temp.push_back(5.42772e-05);
  temp.push_back(8.42695e-05);
  temp.push_back(0.000112201);
  temp.push_back(0.000147315);
  temp.push_back(0.000181674);
  temp.push_back(0.000224435);
  temp.push_back(0.000284745);
  temp.push_back(0.000327223);
  temp.push_back(0.000395938);
  temp.push_back(0.000472805);
  temp.push_back(0.000519986);
  temp.push_back(0.000635296);
  temp.push_back(0.000669197);
  temp.push_back(0.000775362);
  temp.push_back(0.000874534);
  temp.push_back(0.00102977);
  temp.push_back(0.00113617);
  temp.push_back(0.00121311);
  temp.push_back(0.00141158);
  temp.push_back(0.00154837);
  temp.push_back(0.00169857);
  temp.push_back(0.00188519);
  temp.push_back(0.00191032);
  temp.push_back(0.00206575);
  temp.push_back(0.00222722);
  temp.push_back(0.00227325);
  temp.push_back(0.00245166);
  temp.push_back(0.00259628);
  temp.push_back(0.00249795);
  temp.push_back(0.00272387);
  temp.push_back(0.00297299);
  temp.push_back(0.0030973);
  temp.push_back(0.00299177);
  temp.push_back(0.00288901);
  temp.push_back(0.00304009);
  temp.push_back(0.00300757);
  temp.push_back(0.003216);
  temp.push_back(0.00301877);
  temp.push_back(0.00305196);
  temp.push_back(0.00305521);
  temp.push_back(0.00264599);
  temp.push_back(0.00289825);
  temp.push_back(0.00273878);
  temp.push_back(0.0028623);
  temp.push_back(0.00262132);
  temp.push_back(0.00244983);
  temp.push_back(0.00238771);
  temp.push_back(0.00241464);
  temp.push_back(0.0022464);
  temp.push_back(0.00216141);
  temp.push_back(0.00216927);
  temp.push_back(0.00202911);
  temp.push_back(0.00201579);
  temp.push_back(0.00206435);
  temp.push_back(0.00189058);
  temp.push_back(0.0019032);
  temp.push_back(0.00180546);
  temp.push_back(0.00181397);
  temp.push_back(0.00177982);
  temp.push_back(0.00169705);
  temp.push_back(0.00175611);
  temp.push_back(0.00183977);
  temp.push_back(0.00181481);
  temp.push_back(0.0017491);
  temp.push_back(0.00167498);
  temp.push_back(0.00176056);
  temp.push_back(0.00168206);
  temp.push_back(0.00178997);
  temp.push_back(0.00160131);
  temp.push_back(0.00175127);
  temp.push_back(0.00175144);
  temp.push_back(0.00172983);
  temp.push_back(0.00176049);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(1.56003e-05);
  temp.push_back(2.56961e-05);
  temp.push_back(3.98157e-05);
  temp.push_back(7.51062e-05);
  temp.push_back(0.000104808);
  temp.push_back(0.000137126);
  temp.push_back(0.000181945);
  temp.push_back(0.00023428);
  temp.push_back(0.000276685);
  temp.push_back(0.000340639);
  temp.push_back(0.000365219);
  temp.push_back(0.000455067);
  temp.push_back(0.000517277);
  temp.push_back(0.000653087);
  temp.push_back(0.000704082);
  temp.push_back(0.00081798);
  temp.push_back(0.000919737);
  temp.push_back(0.000997528);
  temp.push_back(0.00110852);
  temp.push_back(0.00127015);
  temp.push_back(0.00135952);
  temp.push_back(0.00157116);
  temp.push_back(0.00169714);
  temp.push_back(0.00181507);
  temp.push_back(0.00188281);
  temp.push_back(0.00187945);
  temp.push_back(0.00224784);
  temp.push_back(0.00234832);
  temp.push_back(0.00246759);
  temp.push_back(0.00258435);
  temp.push_back(0.00287211);
  temp.push_back(0.00282387);
  temp.push_back(0.00301207);
  temp.push_back(0.00301706);
  temp.push_back(0.00314772);
  temp.push_back(0.00313356);
  temp.push_back(0.00330207);
  temp.push_back(0.00321868);
  temp.push_back(0.00340633);
  temp.push_back(0.00316687);
  temp.push_back(0.0030367);
  temp.push_back(0.00346355);
  temp.push_back(0.00308358);
  temp.push_back(0.00313356);
  temp.push_back(0.00322158);
  temp.push_back(0.00286642);
  temp.push_back(0.0030063);
  temp.push_back(0.0028074);
  temp.push_back(0.00276659);
  temp.push_back(0.00284753);
  temp.push_back(0.00260215);
  temp.push_back(0.00256294);
  temp.push_back(0.00264208);
  temp.push_back(0.00247397);
  temp.push_back(0.00249761);
  temp.push_back(0.00232896);
  temp.push_back(0.00237722);
  temp.push_back(0.00223927);
  temp.push_back(0.00218478);
  temp.push_back(0.00219469);
  temp.push_back(0.00223608);
  temp.push_back(0.00225412);
  temp.push_back(0.00218023);
  temp.push_back(0.00223298);
  temp.push_back(0.0021489);
  temp.push_back(0.00216245);
  temp.push_back(0.00222762);
  temp.push_back(0.0021608);
  temp.push_back(0.00226409);
  temp.push_back(0.00215114);
  temp.push_back(0.00209354);
  temp.push_back(0.00220478);
  temp.push_back(0.00212817);
  temp.push_back(0.00216157);
  temp.push_back(0.00210899);
  temp.push_back(0.00221498);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(1.98294e-05);
  temp.push_back(2.74883e-05);
  temp.push_back(4.97341e-05);
  temp.push_back(7.8167e-05);
  temp.push_back(0.000115158);
  temp.push_back(0.000159637);
  temp.push_back(0.000198413);
  temp.push_back(0.000249918);
  temp.push_back(0.000307694);
  temp.push_back(0.00037631);
  temp.push_back(0.000420037);
  temp.push_back(0.000530557);
  temp.push_back(0.000574406);
  temp.push_back(0.000678883);
  temp.push_back(0.000783538);
  temp.push_back(0.000828001);
  temp.push_back(0.000945383);
  temp.push_back(0.0010672);
  temp.push_back(0.00119183);
  temp.push_back(0.00134696);
  temp.push_back(0.00155348);
  temp.push_back(0.00167101);
  temp.push_back(0.00181066);
  temp.push_back(0.00187759);
  temp.push_back(0.00199899);
  temp.push_back(0.00225269);
  temp.push_back(0.0023458);
  temp.push_back(0.00253072);
  temp.push_back(0.00270998);
  temp.push_back(0.00263773);
  temp.push_back(0.00294804);
  temp.push_back(0.00294624);
  temp.push_back(0.00321668);
  temp.push_back(0.00332147);
  temp.push_back(0.00334766);
  temp.push_back(0.00330149);
  temp.push_back(0.00340882);
  temp.push_back(0.00347742);
  temp.push_back(0.00352589);
  temp.push_back(0.0035523);
  temp.push_back(0.0035411);
  temp.push_back(0.00340469);
  temp.push_back(0.00335693);
  temp.push_back(0.00345872);
  temp.push_back(0.00316091);
  temp.push_back(0.0031948);
  temp.push_back(0.00307664);
  temp.push_back(0.00320621);
  temp.push_back(0.00304138);
  temp.push_back(0.00299107);
  temp.push_back(0.00281006);
  temp.push_back(0.00278658);
  temp.push_back(0.0028095);
  temp.push_back(0.00268299);
  temp.push_back(0.0027215);
  temp.push_back(0.00264144);
  temp.push_back(0.00260714);
  temp.push_back(0.00255874);
  temp.push_back(0.00266137);
  temp.push_back(0.00248823);
  temp.push_back(0.0026194);
  temp.push_back(0.0026251);
  temp.push_back(0.00245108);
  temp.push_back(0.00254617);
  temp.push_back(0.00264651);
  temp.push_back(0.00247705);
  temp.push_back(0.00252758);
  temp.push_back(0.00262686);
  temp.push_back(0.00248693);
  temp.push_back(0.00254178);
  temp.push_back(0.00258846);
  temp.push_back(0.0025353);
  temp.push_back(0.00240568);
  temp.push_back(0.00250232);
  temp.push_back(0.00251452);
  temp.push_back(0.00258121);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(2.362e-05);
  temp.push_back(3.24457e-05);
  temp.push_back(5.93625e-05);
  temp.push_back(9.71924e-05);
  temp.push_back(0.000126279);
  temp.push_back(0.000176152);
  temp.push_back(0.000224722);
  temp.push_back(0.000290262);
  temp.push_back(0.000351171);
  temp.push_back(0.000439513);
  temp.push_back(0.000497545);
  temp.push_back(0.000596431);
  temp.push_back(0.000686665);
  temp.push_back(0.000716486);
  temp.push_back(0.000835977);
  temp.push_back(0.000943229);
  temp.push_back(0.00104521);
  temp.push_back(0.00120528);
  temp.push_back(0.00140133);
  temp.push_back(0.00144132);
  temp.push_back(0.00164603);
  temp.push_back(0.00179182);
  temp.push_back(0.00188172);
  temp.push_back(0.00203418);
  temp.push_back(0.00224695);
  temp.push_back(0.00224972);
  temp.push_back(0.00259959);
  temp.push_back(0.00249426);
  temp.push_back(0.00280589);
  temp.push_back(0.00286528);
  temp.push_back(0.00309813);
  temp.push_back(0.00321952);
  temp.push_back(0.00324915);
  temp.push_back(0.00312759);
  temp.push_back(0.00338569);
  temp.push_back(0.00351915);
  temp.push_back(0.00336895);
  temp.push_back(0.00358675);
  temp.push_back(0.00361005);
  temp.push_back(0.00356783);
  temp.push_back(0.00381571);
  temp.push_back(0.00358167);
  temp.push_back(0.00355341);
  temp.push_back(0.00369961);
  temp.push_back(0.00345098);
  temp.push_back(0.00345749);
  temp.push_back(0.00366643);
  temp.push_back(0.00328805);
  temp.push_back(0.00317392);
  temp.push_back(0.00325794);
  temp.push_back(0.00325624);
  temp.push_back(0.00298589);
  temp.push_back(0.00313121);
  temp.push_back(0.00307129);
  temp.push_back(0.00303754);
  temp.push_back(0.00299004);
  temp.push_back(0.00283794);
  temp.push_back(0.0029663);
  temp.push_back(0.00278505);
  temp.push_back(0.00294565);
  temp.push_back(0.00295431);
  temp.push_back(0.0029553);
  temp.push_back(0.0027788);
  temp.push_back(0.00272655);
  temp.push_back(0.00267875);
  temp.push_back(0.00288197);
  temp.push_back(0.00284249);
  temp.push_back(0.00281758);
  temp.push_back(0.00285672);
  temp.push_back(0.00260368);
  temp.push_back(0.0027489);
  temp.push_back(0.00270647);
  temp.push_back(0.00300406);
  temp.push_back(0.00298123);
  temp.push_back(0.00284949);
  temp.push_back(0.00275038);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(2.62119e-05);
  temp.push_back(4.06637e-05);
  temp.push_back(6.45721e-05);
  temp.push_back(9.8698e-05);
  temp.push_back(0.000145844);
  temp.push_back(0.000209191);
  temp.push_back(0.000267177);
  temp.push_back(0.000312714);
  temp.push_back(0.000397544);
  temp.push_back(0.000436342);
  temp.push_back(0.000537476);
  temp.push_back(0.000588913);
  temp.push_back(0.000733794);
  temp.push_back(0.000834061);
  temp.push_back(0.000884576);
  temp.push_back(0.001049);
  temp.push_back(0.00109329);
  temp.push_back(0.00128614);
  temp.push_back(0.00139434);
  temp.push_back(0.00155489);
  temp.push_back(0.00169101);
  temp.push_back(0.00192918);
  temp.push_back(0.00192531);
  temp.push_back(0.00211887);
  temp.push_back(0.00227517);
  temp.push_back(0.00250976);
  temp.push_back(0.00268224);
  temp.push_back(0.00276841);
  temp.push_back(0.00287698);
  temp.push_back(0.00299388);
  temp.push_back(0.00306973);
  temp.push_back(0.00337693);
  temp.push_back(0.00332126);
  temp.push_back(0.00344551);
  temp.push_back(0.00364214);
  temp.push_back(0.00354195);
  temp.push_back(0.00371903);
  temp.push_back(0.0037893);
  temp.push_back(0.00381634);
  temp.push_back(0.00377593);
  temp.push_back(0.00382319);
  temp.push_back(0.00378659);
  temp.push_back(0.00379056);
  temp.push_back(0.00392088);
  temp.push_back(0.00396608);
  temp.push_back(0.00361685);
  temp.push_back(0.00356004);
  temp.push_back(0.00377212);
  temp.push_back(0.00359799);
  temp.push_back(0.00349761);
  temp.push_back(0.00345083);
  temp.push_back(0.00341753);
  temp.push_back(0.0031536);
  temp.push_back(0.00329225);
  temp.push_back(0.00318691);
  temp.push_back(0.00342211);
  temp.push_back(0.00324309);
  temp.push_back(0.00343059);
  temp.push_back(0.00313786);
  temp.push_back(0.00309648);
  temp.push_back(0.00308665);
  temp.push_back(0.00315973);
  temp.push_back(0.00320034);
  temp.push_back(0.00327231);
  temp.push_back(0.00315402);
  temp.push_back(0.00306078);
  temp.push_back(0.00307093);
  temp.push_back(0.00304372);
  temp.push_back(0.0032217);
  temp.push_back(0.00315268);
  temp.push_back(0.0031173);
  temp.push_back(0.00308137);
  temp.push_back(0.0031778);
  temp.push_back(0.0030366);
  temp.push_back(0.00304068);
  temp.push_back(0.0033214);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(2.74013e-05);
  temp.push_back(4.22899e-05);
  temp.push_back(6.78705e-05);
  temp.push_back(0.000108305);
  temp.push_back(0.000155136);
  temp.push_back(0.000215087);
  temp.push_back(0.000267491);
  temp.push_back(0.000324307);
  temp.push_back(0.000450428);
  temp.push_back(0.000525255);
  temp.push_back(0.000593983);
  temp.push_back(0.000670897);
  temp.push_back(0.000773512);
  temp.push_back(0.000886391);
  temp.push_back(0.000944693);
  temp.push_back(0.0012177);
  temp.push_back(0.00122543);
  temp.push_back(0.00137642);
  temp.push_back(0.00154086);
  temp.push_back(0.00165634);
  temp.push_back(0.00191029);
  temp.push_back(0.00190192);
  temp.push_back(0.00211654);
  temp.push_back(0.00212919);
  temp.push_back(0.00228748);
  temp.push_back(0.00254164);
  temp.push_back(0.00280712);
  temp.push_back(0.00302764);
  temp.push_back(0.00306686);
  temp.push_back(0.00327206);
  temp.push_back(0.00331136);
  temp.push_back(0.0036489);
  temp.push_back(0.00342041);
  temp.push_back(0.003781);
  temp.push_back(0.00378694);
  temp.push_back(0.00386973);
  temp.push_back(0.00369683);
  temp.push_back(0.0038999);
  temp.push_back(0.00395901);
  temp.push_back(0.00396228);
  temp.push_back(0.00405743);
  temp.push_back(0.00410244);
  temp.push_back(0.00413032);
  temp.push_back(0.00373123);
  temp.push_back(0.00395574);
  temp.push_back(0.00378686);
  temp.push_back(0.00406811);
  temp.push_back(0.00395806);
  temp.push_back(0.00379701);
  temp.push_back(0.00386635);
  temp.push_back(0.00361041);
  temp.push_back(0.00378335);
  temp.push_back(0.00353608);
  temp.push_back(0.00346058);
  temp.push_back(0.00361939);
  temp.push_back(0.00350638);
  temp.push_back(0.00351027);
  temp.push_back(0.00356889);
  temp.push_back(0.00353562);
  temp.push_back(0.00338339);
  temp.push_back(0.00358931);
  temp.push_back(0.0033404);
  temp.push_back(0.0034386);
  temp.push_back(0.00356428);
  temp.push_back(0.00342609);
  temp.push_back(0.00351769);
  temp.push_back(0.0034948);
  temp.push_back(0.0032928);
  temp.push_back(0.00354759);
  temp.push_back(0.0035006);
  temp.push_back(0.00344599);
  temp.push_back(0.00360147);
  temp.push_back(0.00350432);
  temp.push_back(0.00328072);
  temp.push_back(0.00352725);
  temp.push_back(0.00344304);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(3.19736e-05);
  temp.push_back(4.2274e-05);
  temp.push_back(7.42979e-05);
  temp.push_back(0.000115308);
  temp.push_back(0.000181278);
  temp.push_back(0.000238536);
  temp.push_back(0.000295333);
  temp.push_back(0.000372787);
  temp.push_back(0.000450952);
  temp.push_back(0.000549762);
  temp.push_back(0.000587952);
  temp.push_back(0.000717376);
  temp.push_back(0.000823059);
  temp.push_back(0.000929103);
  temp.push_back(0.00101405);
  temp.push_back(0.00120059);
  temp.push_back(0.00129663);
  temp.push_back(0.00138426);
  temp.push_back(0.0016377);
  temp.push_back(0.00177567);
  temp.push_back(0.00193567);
  temp.push_back(0.00193153);
  temp.push_back(0.00217634);
  temp.push_back(0.00226176);
  temp.push_back(0.00269487);
  temp.push_back(0.00260071);
  temp.push_back(0.0029404);
  temp.push_back(0.00314849);
  temp.push_back(0.00329682);
  temp.push_back(0.00359669);
  temp.push_back(0.0033118);
  temp.push_back(0.00346258);
  temp.push_back(0.00376314);
  temp.push_back(0.00361881);
  temp.push_back(0.00386849);
  temp.push_back(0.00410369);
  temp.push_back(0.0039706);
  temp.push_back(0.0041219);
  temp.push_back(0.00407418);
  temp.push_back(0.00410784);
  temp.push_back(0.00417039);
  temp.push_back(0.00396376);
  temp.push_back(0.00402818);
  temp.push_back(0.00408969);
  temp.push_back(0.00398906);
  temp.push_back(0.00407158);
  temp.push_back(0.00407663);
  temp.push_back(0.00429465);
  temp.push_back(0.00399424);
  temp.push_back(0.00398511);
  temp.push_back(0.00413987);
  temp.push_back(0.00402974);
  temp.push_back(0.00396823);
  temp.push_back(0.00406268);
  temp.push_back(0.00397054);
  temp.push_back(0.00406304);
  temp.push_back(0.00381928);
  temp.push_back(0.00373153);
  temp.push_back(0.0038722);
  temp.push_back(0.00372542);
  temp.push_back(0.00371033);
  temp.push_back(0.00370851);
  temp.push_back(0.00370273);
  temp.push_back(0.00375404);
  temp.push_back(0.00368417);
  temp.push_back(0.00370411);
  temp.push_back(0.00360586);
  temp.push_back(0.00374387);
  temp.push_back(0.00370439);
  temp.push_back(0.00362516);
  temp.push_back(0.00377766);
  temp.push_back(0.00360445);
  temp.push_back(0.00360142);
  temp.push_back(0.00368594);
  temp.push_back(0.00358465);
  temp.push_back(0.00361813);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(3.33533e-05);
  temp.push_back(4.71255e-05);
  temp.push_back(7.95104e-05);
  temp.push_back(0.000125354);
  temp.push_back(0.000184391);
  temp.push_back(0.000241585);
  temp.push_back(0.000312129);
  temp.push_back(0.000416606);
  temp.push_back(0.00048679);
  temp.push_back(0.000590734);
  temp.push_back(0.000680798);
  temp.push_back(0.00074592);
  temp.push_back(0.000925237);
  temp.push_back(0.000945055);
  temp.push_back(0.00111491);
  temp.push_back(0.00124566);
  temp.push_back(0.00130333);
  temp.push_back(0.00159285);
  temp.push_back(0.00160697);
  temp.push_back(0.00173291);
  temp.push_back(0.00205625);
  temp.push_back(0.0020842);
  temp.push_back(0.00234011);
  temp.push_back(0.00257964);
  temp.push_back(0.00271621);
  temp.push_back(0.00274754);
  temp.push_back(0.00300009);
  temp.push_back(0.00313262);
  temp.push_back(0.0033966);
  temp.push_back(0.00370825);
  temp.push_back(0.00365397);
  temp.push_back(0.00380956);
  temp.push_back(0.00406663);
  temp.push_back(0.00390965);
  temp.push_back(0.00400766);
  temp.push_back(0.0040838);
  temp.push_back(0.00430457);
  temp.push_back(0.00430605);
  temp.push_back(0.00452344);
  temp.push_back(0.00444883);
  temp.push_back(0.00462339);
  temp.push_back(0.00441296);
  temp.push_back(0.0044215);
  temp.push_back(0.00437965);
  temp.push_back(0.00450718);
  temp.push_back(0.00411557);
  temp.push_back(0.0044587);
  temp.push_back(0.00428602);
  temp.push_back(0.00415902);
  temp.push_back(0.00433486);
  temp.push_back(0.0043004);
  temp.push_back(0.00407849);
  temp.push_back(0.00423246);
  temp.push_back(0.00394595);
  temp.push_back(0.00402799);
  temp.push_back(0.00390578);
  temp.push_back(0.00403325);
  temp.push_back(0.00401427);
  temp.push_back(0.00413036);
  temp.push_back(0.00396507);
  temp.push_back(0.00397724);
  temp.push_back(0.00405608);
  temp.push_back(0.00410989);
  temp.push_back(0.00390445);
  temp.push_back(0.00408041);
  temp.push_back(0.00374423);
  temp.push_back(0.00380687);
  temp.push_back(0.00396012);
  temp.push_back(0.00403351);
  temp.push_back(0.00397022);
  temp.push_back(0.00420906);
  temp.push_back(0.00395898);
  temp.push_back(0.00414565);
  temp.push_back(0.00419401);
  temp.push_back(0.00398256);
  temp.push_back(0.00405866);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(3.60154e-05);
  temp.push_back(4.78073e-05);
  temp.push_back(8.81868e-05);
  temp.push_back(0.00014119);
  temp.push_back(0.000199197);
  temp.push_back(0.000271447);
  temp.push_back(0.000348848);
  temp.push_back(0.000414794);
  temp.push_back(0.000517428);
  temp.push_back(0.000613861);
  temp.push_back(0.000706405);
  temp.push_back(0.000830507);
  temp.push_back(0.000954556);
  temp.push_back(0.00102745);
  temp.push_back(0.00115016);
  temp.push_back(0.00130855);
  temp.push_back(0.00141543);
  temp.push_back(0.00154704);
  temp.push_back(0.00183892);
  temp.push_back(0.00182328);
  temp.push_back(0.00211743);
  temp.push_back(0.00238445);
  temp.push_back(0.00251855);
  temp.push_back(0.00259965);
  temp.push_back(0.00276547);
  temp.push_back(0.00279785);
  temp.push_back(0.00310652);
  temp.push_back(0.00332247);
  temp.push_back(0.00348554);
  temp.push_back(0.00340047);
  temp.push_back(0.00379133);
  temp.push_back(0.00394976);
  temp.push_back(0.00411556);
  temp.push_back(0.00395808);
  temp.push_back(0.00413605);
  temp.push_back(0.00434804);
  temp.push_back(0.00454974);
  temp.push_back(0.00448955);
  temp.push_back(0.00475876);
  temp.push_back(0.00457681);
  temp.push_back(0.00444923);
  temp.push_back(0.00441159);
  temp.push_back(0.00433255);
  temp.push_back(0.0044053);
  temp.push_back(0.00467223);
  temp.push_back(0.00444452);
  temp.push_back(0.00429894);
  temp.push_back(0.00466592);
  temp.push_back(0.004516);
  temp.push_back(0.00446295);
  temp.push_back(0.00435999);
  temp.push_back(0.00467332);
  temp.push_back(0.00423821);
  temp.push_back(0.00443203);
  temp.push_back(0.00434006);
  temp.push_back(0.00429902);
  temp.push_back(0.0042231);
  temp.push_back(0.00423514);
  temp.push_back(0.00451164);
  temp.push_back(0.00436348);
  temp.push_back(0.00436551);
  temp.push_back(0.00425205);
  temp.push_back(0.00429035);
  temp.push_back(0.00411415);
  temp.push_back(0.00427361);
  temp.push_back(0.0039277);
  temp.push_back(0.00405848);
  temp.push_back(0.00429804);
  temp.push_back(0.00408766);
  temp.push_back(0.00384011);
  temp.push_back(0.00435199);
  temp.push_back(0.004093);
  temp.push_back(0.00406005);
  temp.push_back(0.00431062);
  temp.push_back(0.00411773);
  temp.push_back(0.00435653);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(3.94789e-05);
  temp.push_back(5.3769e-05);
  temp.push_back(9.3841e-05);
  temp.push_back(0.000149774);
  temp.push_back(0.000215262);
  temp.push_back(0.000267695);
  temp.push_back(0.000370277);
  temp.push_back(0.000473963);
  temp.push_back(0.000548657);
  temp.push_back(0.000625464);
  temp.push_back(0.000737716);
  temp.push_back(0.000895534);
  temp.push_back(0.00100207);
  temp.push_back(0.00111602);
  temp.push_back(0.00123458);
  temp.push_back(0.00138359);
  temp.push_back(0.00158882);
  temp.push_back(0.00176246);
  temp.push_back(0.00177764);
  temp.push_back(0.00199677);
  temp.push_back(0.00211274);
  temp.push_back(0.00240314);
  temp.push_back(0.00238892);
  temp.push_back(0.00266706);
  temp.push_back(0.00292668);
  temp.push_back(0.00314009);
  temp.push_back(0.00330102);
  temp.push_back(0.00354671);
  temp.push_back(0.00371555);
  temp.push_back(0.00383521);
  temp.push_back(0.00392384);
  temp.push_back(0.00415042);
  temp.push_back(0.00417837);
  temp.push_back(0.00410063);
  temp.push_back(0.00436123);
  temp.push_back(0.00453869);
  temp.push_back(0.00456717);
  temp.push_back(0.00432658);
  temp.push_back(0.00474763);
  temp.push_back(0.00479594);
  temp.push_back(0.00495067);
  temp.push_back(0.00486827);
  temp.push_back(0.00501981);
  temp.push_back(0.00467952);
  temp.push_back(0.00496726);
  temp.push_back(0.00471759);
  temp.push_back(0.00487664);
  temp.push_back(0.00477905);
  temp.push_back(0.00457876);
  temp.push_back(0.00478556);
  temp.push_back(0.00470292);
  temp.push_back(0.00474277);
  temp.push_back(0.00474074);
  temp.push_back(0.00476266);
  temp.push_back(0.00459758);
  temp.push_back(0.00465773);
  temp.push_back(0.00469381);
  temp.push_back(0.004336);
  temp.push_back(0.00439833);
  temp.push_back(0.00434989);
  temp.push_back(0.00452046);
  temp.push_back(0.00444761);
  temp.push_back(0.00452204);
  temp.push_back(0.00433173);
  temp.push_back(0.00428761);
  temp.push_back(0.00454501);
  temp.push_back(0.00466898);
  temp.push_back(0.00449207);
  temp.push_back(0.00478147);
  temp.push_back(0.00448187);
  temp.push_back(0.00427879);
  temp.push_back(0.0045072);
  temp.push_back(0.00446825);
  temp.push_back(0.00442717);
  temp.push_back(0.00450486);
  temp.push_back(0.004514);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(4.07374e-05);
  temp.push_back(5.84427e-05);
  temp.push_back(0.000102379);
  temp.push_back(0.000154485);
  temp.push_back(0.00022167);
  temp.push_back(0.000291138);
  temp.push_back(0.000366195);
  temp.push_back(0.00049284);
  temp.push_back(0.000555038);
  temp.push_back(0.000718216);
  temp.push_back(0.000815782);
  temp.push_back(0.000929635);
  temp.push_back(0.00105284);
  temp.push_back(0.00117775);
  temp.push_back(0.00128659);
  temp.push_back(0.0014848);
  temp.push_back(0.00161139);
  temp.push_back(0.00186044);
  temp.push_back(0.00193622);
  temp.push_back(0.0021136);
  temp.push_back(0.00218413);
  temp.push_back(0.00254649);
  temp.push_back(0.00257878);
  temp.push_back(0.00286384);
  temp.push_back(0.00277422);
  temp.push_back(0.00337728);
  temp.push_back(0.0034849);
  temp.push_back(0.00340672);
  temp.push_back(0.00376872);
  temp.push_back(0.00397498);
  temp.push_back(0.00407663);
  temp.push_back(0.00416444);
  temp.push_back(0.00447438);
  temp.push_back(0.00421777);
  temp.push_back(0.00434779);
  temp.push_back(0.00457815);
  temp.push_back(0.00510487);
  temp.push_back(0.004957);
  temp.push_back(0.0050608);
  temp.push_back(0.00482047);
  temp.push_back(0.00479584);
  temp.push_back(0.00506816);
  temp.push_back(0.00485049);
  temp.push_back(0.00502603);
  temp.push_back(0.00496904);
  temp.push_back(0.0048041);
  temp.push_back(0.00498567);
  temp.push_back(0.00505111);
  temp.push_back(0.00495093);
  temp.push_back(0.00493909);
  temp.push_back(0.00518096);
  temp.push_back(0.00475333);
  temp.push_back(0.00509383);
  temp.push_back(0.00455324);
  temp.push_back(0.00501535);
  temp.push_back(0.00477645);
  temp.push_back(0.00467969);
  temp.push_back(0.00482058);
  temp.push_back(0.0049352);
  temp.push_back(0.00479813);
  temp.push_back(0.00456841);
  temp.push_back(0.00495553);
  temp.push_back(0.00475263);
  temp.push_back(0.00464536);
  temp.push_back(0.0045788);
  temp.push_back(0.00482641);
  temp.push_back(0.00468534);
  temp.push_back(0.00458609);
  temp.push_back(0.00469523);
  temp.push_back(0.00471367);
  temp.push_back(0.00480113);
  temp.push_back(0.00477874);
  temp.push_back(0.00487724);
  temp.push_back(0.00448248);
  temp.push_back(0.00457135);
  temp.push_back(0.00463623);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(4.19494e-05);
  temp.push_back(5.81133e-05);
  temp.push_back(0.000106598);
  temp.push_back(0.000167045);
  temp.push_back(0.000245601);
  temp.push_back(0.000309417);
  temp.push_back(0.000402556);
  temp.push_back(0.000531311);
  temp.push_back(0.000599414);
  temp.push_back(0.000724794);
  temp.push_back(0.000843293);
  temp.push_back(0.00102089);
  temp.push_back(0.00110244);
  temp.push_back(0.00119442);
  temp.push_back(0.0014845);
  temp.push_back(0.00144809);
  temp.push_back(0.00163904);
  temp.push_back(0.00179288);
  temp.push_back(0.00208349);
  temp.push_back(0.00218156);
  temp.push_back(0.00241731);
  temp.push_back(0.00260733);
  temp.push_back(0.0028223);
  temp.push_back(0.00303541);
  temp.push_back(0.00317796);
  temp.push_back(0.00318073);
  temp.push_back(0.00354012);
  temp.push_back(0.0038304);
  temp.push_back(0.00394061);
  temp.push_back(0.00388992);
  temp.push_back(0.00415685);
  temp.push_back(0.004271);
  temp.push_back(0.00454353);
  temp.push_back(0.00436867);
  temp.push_back(0.00471855);
  temp.push_back(0.00486456);
  temp.push_back(0.00458562);
  temp.push_back(0.00467012);
  temp.push_back(0.00488229);
  temp.push_back(0.0051219);
  temp.push_back(0.00525306);
  temp.push_back(0.00520996);
  temp.push_back(0.00512736);
  temp.push_back(0.00506308);
  temp.push_back(0.00509758);
  temp.push_back(0.00513977);
  temp.push_back(0.00506956);
  temp.push_back(0.0050327);
  temp.push_back(0.00528484);
  temp.push_back(0.00541921);
  temp.push_back(0.00486394);
  temp.push_back(0.00532431);
  temp.push_back(0.00491166);
  temp.push_back(0.0054032);
  temp.push_back(0.00515276);
  temp.push_back(0.00499206);
  temp.push_back(0.00507303);
  temp.push_back(0.00528015);
  temp.push_back(0.00502935);
  temp.push_back(0.00498204);
  temp.push_back(0.00482074);
  temp.push_back(0.0050209);
  temp.push_back(0.00506188);
  temp.push_back(0.00501215);
  temp.push_back(0.00523337);
  temp.push_back(0.00501484);
  temp.push_back(0.00482134);
  temp.push_back(0.00510347);
  temp.push_back(0.00500609);
  temp.push_back(0.00508656);
  temp.push_back(0.00505679);
  temp.push_back(0.00526603);
  temp.push_back(0.00491982);
  temp.push_back(0.00505811);
  temp.push_back(0.00494347);
  temp.push_back(0.00506382);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(4.58461e-05);
  temp.push_back(6.56593e-05);
  temp.push_back(0.000106249);
  temp.push_back(0.000177259);
  temp.push_back(0.000242624);
  temp.push_back(0.000318801);
  temp.push_back(0.000429442);
  temp.push_back(0.000524763);
  temp.push_back(0.000627853);
  temp.push_back(0.000767861);
  temp.push_back(0.000914838);
  temp.push_back(0.00101068);
  temp.push_back(0.00119111);
  temp.push_back(0.00124833);
  temp.push_back(0.00147471);
  temp.push_back(0.00151766);
  temp.push_back(0.00175504);
  temp.push_back(0.0019513);
  temp.push_back(0.00225088);
  temp.push_back(0.00232528);
  temp.push_back(0.002477);
  temp.push_back(0.0026629);
  temp.push_back(0.00282126);
  temp.push_back(0.00305865);
  temp.push_back(0.00328165);
  temp.push_back(0.00342969);
  temp.push_back(0.00368517);
  temp.push_back(0.00412023);
  temp.push_back(0.00408107);
  temp.push_back(0.00442514);
  temp.push_back(0.00419995);
  temp.push_back(0.00445839);
  temp.push_back(0.00480241);
  temp.push_back(0.00469063);
  temp.push_back(0.00488558);
  temp.push_back(0.00505738);
  temp.push_back(0.00534401);
  temp.push_back(0.00529189);
  temp.push_back(0.00503914);
  temp.push_back(0.00547271);
  temp.push_back(0.00532523);
  temp.push_back(0.0056474);
  temp.push_back(0.00528356);
  temp.push_back(0.00537977);
  temp.push_back(0.00528494);
  temp.push_back(0.00555228);
  temp.push_back(0.00522019);
  temp.push_back(0.00551423);
  temp.push_back(0.00532466);
  temp.push_back(0.00562846);
  temp.push_back(0.00531636);
  temp.push_back(0.00534452);
  temp.push_back(0.00546899);
  temp.push_back(0.00540307);
  temp.push_back(0.00522696);
  temp.push_back(0.00543343);
  temp.push_back(0.00496315);
  temp.push_back(0.0054688);
  temp.push_back(0.0054534);
  temp.push_back(0.00529457);
  temp.push_back(0.00527906);
  temp.push_back(0.00505608);
  temp.push_back(0.00532934);
  temp.push_back(0.00540523);
  temp.push_back(0.00546589);
  temp.push_back(0.00541721);
  temp.push_back(0.00518411);
  temp.push_back(0.0050658);
  temp.push_back(0.00538019);
  temp.push_back(0.00523502);
  temp.push_back(0.00518119);
  temp.push_back(0.00541954);
  temp.push_back(0.00556926);
  temp.push_back(0.00516485);
  temp.push_back(0.00543988);
  temp.push_back(0.00534757);
  theErrSamples.push_back(temp);
  temp.clear();

  temp.push_back(4.87191e-05);
  temp.push_back(6.79557e-05);
  temp.push_back(0.00011585);
  temp.push_back(0.000170012);
  temp.push_back(0.000249995);
  temp.push_back(0.000358082);
  temp.push_back(0.000441529);
  temp.push_back(0.000573455);
  temp.push_back(0.000686492);
  temp.push_back(0.000768395);
  temp.push_back(0.000922129);
  temp.push_back(0.0010261);
  temp.push_back(0.00121379);
  temp.push_back(0.00132394);
  temp.push_back(0.00147006);
  temp.push_back(0.0016026);
  temp.push_back(0.00185945);
  temp.push_back(0.00205459);
  temp.push_back(0.00224833);
  temp.push_back(0.00239438);
  temp.push_back(0.00258948);
  temp.push_back(0.00285399);
  temp.push_back(0.00293392);
  temp.push_back(0.00310816);
  temp.push_back(0.0035526);
  temp.push_back(0.00353642);
  temp.push_back(0.003711);
  temp.push_back(0.00392946);
  temp.push_back(0.00417796);
  temp.push_back(0.00440017);
  temp.push_back(0.00474434);
  temp.push_back(0.00469919);
  temp.push_back(0.00476225);
  temp.push_back(0.00474708);
  temp.push_back(0.00515252);
  temp.push_back(0.0053985);
  temp.push_back(0.00515269);
  temp.push_back(0.00535194);
  temp.push_back(0.00519113);
  temp.push_back(0.00531732);
  temp.push_back(0.00545661);
  temp.push_back(0.00562353);
  temp.push_back(0.00547828);
  temp.push_back(0.00590485);
  temp.push_back(0.0054847);
  temp.push_back(0.00575722);
  temp.push_back(0.00578941);
  temp.push_back(0.0055163);
  temp.push_back(0.00565193);
  temp.push_back(0.00559088);
  temp.push_back(0.00578905);
  temp.push_back(0.00568463);
  temp.push_back(0.00538114);
  temp.push_back(0.00569508);
  temp.push_back(0.00534805);
  temp.push_back(0.00521921);
  temp.push_back(0.00568289);
  temp.push_back(0.0055866);
  temp.push_back(0.00561974);
  temp.push_back(0.00548487);
  temp.push_back(0.00558498);
  temp.push_back(0.00536631);
  temp.push_back(0.00541084);
  temp.push_back(0.00577058);
  temp.push_back(0.00540762);
  temp.push_back(0.00545929);
  temp.push_back(0.00571948);
  temp.push_back(0.00546965);
  temp.push_back(0.00557277);
  temp.push_back(0.00554671);
  temp.push_back(0.00521175);
  temp.push_back(0.00546106);
  temp.push_back(0.00541317);
  temp.push_back(0.00569926);
  temp.push_back(0.0055058);
  temp.push_back(0.00537131);
  theErrSamples.push_back(temp);
  temp.clear();

}