HcalDeterministicFit Class Reference

#include <HcalDeterministicFit.h>

Public Types
enum	FType { shapeLandau, shape205, shape206, shape207 }

Public Member Functions
void	get205Frac (float tStart, float tEnd, float &sum) const
void	get206Frac (float tStart, float tEnd, float &sum) const
void	get207Frac (float tStart, float tEnd, float &sum) const
void	getFrac (float, float, float &, FType) const
void	getLandauFrac (float tStart, float tEnd, float &sum) const
	HcalDeterministicFit ()
void	init (HcalTimeSlew::ParaSource tsParam, HcalTimeSlew::BiasSetting bias, bool iApplyTimeSlew, double respCorr)
void	phase1Apply (const HBHEChannelInfo &channelData, float &reconstructedEnergy, float &reconstructedTime, const HcalTimeSlew *hcalTimeSlew_delay) const
	~HcalDeterministicFit ()

Private Attributes
bool	applyTimeSlew_
PedestalSub	fPedestalSubFxn_
double	frespCorr_
HcalTimeSlew::ParaSource	fTimeSlew_
HcalTimeSlew::BiasSetting	fTimeSlewBias_

Static Private Attributes
static constexpr int	HcalRegion [2] = {16, 17}
static constexpr float	invGpar [3] = {-13.11, 11.29, 5.133}
static constexpr float	landauFrac [125]
static constexpr float	negThresh [2] = {-3., 15.}
static constexpr float	rCorr [2] = {0.95, 0.95}
static constexpr float	rCorrSiPM [2] = {1., 1.}
static constexpr float	siPM205Frac [125]
static constexpr float	siPM206Frac [125]
static constexpr float	siPM207Frac [125]
static constexpr int	tsWidth = 25

Detailed Description

Definition at line 16 of file HcalDeterministicFit.h.

Member Enumeration Documentation

FType

enum HcalDeterministicFit::FType

Enumerator
shapeLandau
shape205
shape206
shape207

Definition at line 21 of file HcalDeterministicFit.h.

{ shapeLandau, shape205, shape206, shape207 };

Constructor & Destructor Documentation

HcalDeterministicFit()

HcalDeterministicFit::HcalDeterministicFit

(

)

Definition at line 14 of file HcalDeterministicFit.cc.

{}

~HcalDeterministicFit()

HcalDeterministicFit::~HcalDeterministicFit

(

)

Definition at line 16 of file HcalDeterministicFit.cc.

{}

Member Function Documentation

get205Frac()

void HcalDeterministicFit::get205Frac	(	float	tStart,
		float	tEnd,
		float &	sum
	)		const

Definition at line 43 of file HcalDeterministicFit.cc.

                                                                                {
  if (std::abs(tStart - tEnd - tsWidth) < 0.1f) {
    sum = 0.f;
    return;
  }
  sum = siPM205Frac[int(ceil(tStart + tsWidth))];
  return;
}

References funct::abs(), reco::ceil(), f, and createfilelist::int.

get206Frac()

void HcalDeterministicFit::get206Frac	(	float	tStart,
		float	tEnd,
		float &	sum
	)		const

Definition at line 53 of file HcalDeterministicFit.cc.

                                                                                {
  if (std::abs(tStart - tEnd - tsWidth) < 0.1f) {
    sum = 0.f;
    return;
  }
  sum = siPM206Frac[int(ceil(tStart + tsWidth))];
  return;
}

References funct::abs(), reco::ceil(), f, and createfilelist::int.

get207Frac()

void HcalDeterministicFit::get207Frac	(	float	tStart,
		float	tEnd,
		float &	sum
	)		const

Definition at line 63 of file HcalDeterministicFit.cc.

                                                                                {
  if (std::abs(tStart - tEnd - tsWidth) < 0.1f) {
    sum = 0.f;
    return;
  }
  sum = siPM207Frac[int(ceil(tStart + tsWidth))];
  return;
}

References funct::abs(), reco::ceil(), f, and createfilelist::int.

getFrac()

void HcalDeterministicFit::getFrac	(	float	tStart,
		float	tEnd,
		float &	sum,
		FType	fType
	)		const

Definition at line 72 of file HcalDeterministicFit.cc.

                                                                                          {
  switch (fType) {
    case shape205:
      get205Frac(tStart, tEnd, sum);
      break;
    case shape206:
      get206Frac(tStart, tEnd, sum);
      break;
    case shape207:
      get207Frac(tStart, tEnd, sum);
      break;
    case shapeLandau:
      getLandauFrac(tStart, tEnd, sum);
      break;
  }
}

getLandauFrac()

void HcalDeterministicFit::getLandauFrac	(	float	tStart,
		float	tEnd,
		float &	sum
	)		const

Definition at line 33 of file HcalDeterministicFit.cc.

                                                                                   {
  if (std::abs(tStart - tEnd - tsWidth) < 0.1f) {
    sum = 0.f;
    return;
  }
  sum = landauFrac[int(ceil(tStart + tsWidth))];
  return;
}

References funct::abs(), reco::ceil(), f, and createfilelist::int.

init()

void HcalDeterministicFit::init	(	HcalTimeSlew::ParaSource	tsParam,
		HcalTimeSlew::BiasSetting	bias,
		bool	iApplyTimeSlew,
		double	respCorr
	)

Definition at line 18 of file HcalDeterministicFit.cc.

                                                 {
  fTimeSlew_ = tsParam;
  fTimeSlewBias_ = bias;
  applyTimeSlew_ = iApplyTimeSlew;
  frespCorr_ = respCorr;
}

phase1Apply()

void HcalDeterministicFit::phase1Apply	(	const HBHEChannelInfo &	channelData,
		float &	reconstructedEnergy,
		float &	reconstructedTime,
		const HcalTimeSlew *	hcalTimeSlew_delay
	)		const

Definition at line 89 of file HcalDeterministicFit.cc.

                                                                                     {
  unsigned int soi = channelData.soi();
 
  std::vector<double> corrCharge;
  corrCharge.reserve(channelData.nSamples());
  std::vector<double> inputCharge;
  inputCharge.reserve(channelData.nSamples());
  std::vector<double> inputPedestal;
  inputPedestal.reserve(channelData.nSamples());
  std::vector<double> inputNoise;
  inputNoise.reserve(channelData.nSamples());
 
  double gainCorr = 0;
  double respCorr = 0;
 
  for (unsigned int ip = 0; ip < channelData.nSamples(); ip++) {
    double charge = channelData.tsRawCharge(ip);
    double ped = channelData.tsPedestal(ip);
    double noise = channelData.tsPedestalWidth(ip);
    double gain = channelData.tsGain(ip);
 
    gainCorr = gain;
    inputCharge.push_back(charge);
    inputPedestal.push_back(ped);
    inputNoise.push_back(noise);
  }
 
  fPedestalSubFxn_.calculate(inputCharge, inputPedestal, inputNoise, corrCharge, soi, channelData.nSamples());
 
  if (fTimeSlew_ == 0)
    respCorr = 1.0;
  else if (fTimeSlew_ == 1)
    channelData.hasTimeInfo() ? respCorr = rCorrSiPM[0] : respCorr = rCorr[0];
  else if (fTimeSlew_ == 2)
    channelData.hasTimeInfo() ? respCorr = rCorrSiPM[1] : respCorr = rCorr[1];
  else if (fTimeSlew_ == 3)
    respCorr = frespCorr_;
 
  float tsShift3, tsShift4, tsShift5;
  tsShift3 = 0.f, tsShift4 = 0.f, tsShift5 = 0.f;
 
  if (applyTimeSlew_) {
    tsShift3 = hcalTimeSlew_delay->delay(inputCharge[soi - 1], fTimeSlew_, fTimeSlewBias_, !channelData.hasTimeInfo());
    tsShift4 = hcalTimeSlew_delay->delay(inputCharge[soi], fTimeSlew_, fTimeSlewBias_, !channelData.hasTimeInfo());
    tsShift5 = hcalTimeSlew_delay->delay(inputCharge[soi + 1], fTimeSlew_, fTimeSlewBias_, !channelData.hasTimeInfo());
  }
 
  float ch3, ch4, ch5, i3, n3, nn3, i4, n4, i5, n5;
  ch4 = 0.f, i3 = 0.f, n3 = 0.f, nn3 = 0.f, i4 = 0.f, n4 = 0.f, i5 = 0.f, n5 = 0.f;
 
  FType fType;
  if (channelData.hasTimeInfo() && channelData.recoShape() == 205)
    fType = shape205;
  else if (channelData.hasTimeInfo() && channelData.recoShape() == 206)
    fType = shape206;
  else if (channelData.hasTimeInfo() && channelData.recoShape() == 207)
    fType = shape207;
  else
    fType = shapeLandau;
 
  getFrac(-tsShift3, -tsShift3 + tsWidth, i3, fType);
  getFrac(-tsShift3 + tsWidth, -tsShift3 + tsWidth * 2, n3, fType);
  getFrac(-tsShift3 + tsWidth * 2, -tsShift3 + tsWidth * 3, nn3, fType);
 
  getFrac(-tsShift4, -tsShift4 + tsWidth, i4, fType);
  getFrac(-tsShift4 + tsWidth, -tsShift4 + tsWidth * 2, n4, fType);
 
  getFrac(-tsShift5, -tsShift5 + tsWidth, i5, fType);
  getFrac(-tsShift5 + tsWidth, -tsShift5 + tsWidth * 2, n5, fType);
 
  if (i3 != 0 && i4 != 0 && i5 != 0) {
    ch3 = corrCharge[soi - 1] / i3;
    ch4 = (i3 * corrCharge[soi] - n3 * corrCharge[soi - 1]) / (i3 * i4);
    ch5 = (n3 * n4 * corrCharge[soi - 1] - i4 * nn3 * corrCharge[soi - 1] - i3 * n4 * corrCharge[soi] +
           i3 * i4 * corrCharge[soi + 1]) /
          (i3 * i4 * i5);
 
    if (ch3 < negThresh[0]) {
      ch3 = negThresh[0];
      ch4 = corrCharge[soi] / i4;
      ch5 = (i4 * corrCharge[soi + 1] - n4 * corrCharge[soi]) / (i4 * i5);
    }
    if (ch5 < negThresh[0] && ch4 > negThresh[1]) {
      double ratio = (corrCharge[soi] - ch3 * i3) / (corrCharge[soi + 1] - negThresh[0] * i5);
      if (ratio < 5 && ratio > 0.5) {
        double invG = invGpar[0] + invGpar[1] * std::sqrt(2 * std::log(invGpar[2] / ratio));
        float iG = 0.f;
        getFrac(-invG, -invG + tsWidth, iG, fType);
        if (iG != 0) {
          ch4 = (corrCharge[soi] - ch3 * n3) / (iG);
          tsShift4 = invG;
        }
      }
    }
  }
 
  if (ch4 < 1) {
    ch4 = 0.f;
  }
 
  reconstructedEnergy = ch4 * gainCorr * respCorr;
  reconstructedTime = tsShift4;
}

References ALCARECOTkAlJpsiMuMu_cff::charge, HcalTimeSlew::delay(), PedestalClient_cfi::gain, HBHEChannelInfo::hasTimeInfo(), testProducerWithPsetDescEmpty_cfi::i3, dqm-mbProfile::log, hgcalDigitizer_cfi::noise, HBHEChannelInfo::nSamples(), particleFlowDisplacedVertex_cfi::ratio, HBHEChannelInfo::recoShape(), HBHEChannelInfo::soi(), mathSSE::sqrt(), HBHEChannelInfo::tsGain(), HBHEChannelInfo::tsPedestal(), HBHEChannelInfo::tsPedestalWidth(), and HBHEChannelInfo::tsRawCharge().

Referenced by SimpleHBHEPhase1Algo::reconstruct().

Member Data Documentation

applyTimeSlew_

bool HcalDeterministicFit::applyTimeSlew_

private

Definition at line 41 of file HcalDeterministicFit.h.

fPedestalSubFxn_

PedestalSub HcalDeterministicFit::fPedestalSubFxn_

private

Definition at line 40 of file HcalDeterministicFit.h.

frespCorr_

double HcalDeterministicFit::frespCorr_

private

Definition at line 42 of file HcalDeterministicFit.h.

fTimeSlew_

HcalTimeSlew::ParaSource HcalDeterministicFit::fTimeSlew_

private

Definition at line 38 of file HcalDeterministicFit.h.

fTimeSlewBias_

HcalTimeSlew::BiasSetting HcalDeterministicFit::fTimeSlewBias_

private

Definition at line 39 of file HcalDeterministicFit.h.

HcalRegion

constexpr int HcalDeterministicFit::HcalRegion = {16, 17}

staticconstexprprivate

Definition at line 44 of file HcalDeterministicFit.h.

invGpar

constexpr float HcalDeterministicFit::invGpar = {-13.11, 11.29, 5.133}

staticconstexprprivate

Definition at line 47 of file HcalDeterministicFit.h.

landauFrac

constexpr float HcalDeterministicFit::landauFrac

staticconstexprprivate

Initial value:

= {
      0,          7.6377e-05, 0.000418655, 0.00153692, 0.00436844, 0.0102076,  0.0204177,   0.0360559,  0.057596,
      0.0848493,  0.117069,   0.153152,    0.191858,   0.23198,    0.272461,   0.312438,    0.351262,   0.388476,
      0.423788,   0.457036,   0.488159,    0.517167,   0.54412,    0.569112,   0.592254,    0.613668,   0.633402,
      0.651391,   0.667242,   0.680131,    0.688868,   0.692188,   0.689122,   0.67928,     0.662924,   0.64087,
      0.614282,   0.584457,   0.552651,    0.51997,    0.487317,   0.455378,   0.424647,    0.395445,   0.367963,
      0.342288,   0.318433,   0.29636,     0.275994,   0.257243,   0.24,       0.224155,    0.2096,     0.196227,
      0.183937,   0.172635,   0.162232,    0.15265,    0.143813,   0.135656,   0.128117,    0.12114,    0.114677,
      0.108681,   0.103113,   0.0979354,   0.0931145,  0.0886206,  0.0844264,  0.0805074,   0.0768411,  0.0734075,
      0.0701881,  0.0671664,  0.0643271,   0.0616564,  0.0591418,  0.0567718,  0.054536,    0.0524247,  0.0504292,
      0.0485414,  0.046754,   0.0450602,   0.0434538,  0.041929,   0.0404806,  0.0391037,   0.0377937,  0.0365465,
      0.0353583,  0.0342255,  0.0331447,   0.032113,   0.0311274,  0.0301854,  0.0292843,   0.0284221,  0.0275964,
      0.0268053,  0.0253052,  0.0238536,   0.0224483,  0.0210872,  0.0197684,  0.0184899,   0.01725,    0.0160471,
      0.0148795,  0.0137457,  0.0126445,   0.0115743,  0.0105341,  0.00952249, 0.00853844,  0.00758086, 0.00664871,
      0.00574103, 0.00485689, 0.00399541,  0.00315576, 0.00233713, 0.00153878, 0.000759962, 0}

Definition at line 51 of file HcalDeterministicFit.h.

negThresh

constexpr float HcalDeterministicFit::negThresh = {-3., 15.}

staticconstexprprivate

Definition at line 46 of file HcalDeterministicFit.h.

rCorr

constexpr float HcalDeterministicFit::rCorr = {0.95, 0.95}

staticconstexprprivate

Definition at line 48 of file HcalDeterministicFit.h.

rCorrSiPM

constexpr float HcalDeterministicFit::rCorrSiPM = {1., 1.}

staticconstexprprivate

Definition at line 49 of file HcalDeterministicFit.h.

siPM205Frac

constexpr float HcalDeterministicFit::siPM205Frac

staticconstexprprivate

Initial value:

= {
      0,         0,         0,         0,         0.00133129, 0.00444633, 0.0115,    0.0243992, 0.0443875, 0.0716386,
      0.105298,  0.143832,  0.185449,  0.228439,  0.271367,   0.31315,    0.353041,  0.390587,  0.425555,  0.45788,
      0.487604,  0.514843,  0.539752,  0.562504,  0.583282,   0.602263,   0.619612,  0.635457,  0.649765,  0.66208,
      0.671249,  0.675509,  0.673048,  0.662709,  0.644394,   0.619024,   0.588194,  0.55375,   0.517448,  0.480768,
      0.444831,  0.410418,  0.378015,  0.347879,  0.320103,   0.294667,   0.271474,  0.250391,  0.231257,  0.213907,
      0.198178,  0.183914,  0.170967,  0.159205,  0.148505,   0.138758,   0.129864,  0.121737,  0.114299,  0.107478,
      0.101214,  0.0954507, 0.0901402, 0.0852385, 0.0807069,  0.0765108,  0.0726194, 0.0690052, 0.0656435, 0.0625123,
      0.0595916, 0.0568637, 0.0543125, 0.0519236, 0.0496838,  0.0475815,  0.0456058, 0.0437472, 0.0419966, 0.0403463,
      0.0387887, 0.0373173, 0.0359259, 0.034609,  0.0333615,  0.0321786,  0.0310561, 0.02999,   0.0289767, 0.0280127,
      0.0270951, 0.0262209, 0.0253875, 0.0245923, 0.0238333,  0.0231082,  0.022415,  0.021752,  0.0211174, 0.0205097,
      0.0199274, 0.0193692, 0.0188336, 0.0183196, 0.017826,   0.0173518,  0.0168959, 0.0164575, 0.0160356, 0.0156296,
      0.0152385, 0.0148617, 0.0144984, 0.0141482, 0.0138103,  0.0134842,  0.0131693, 0.0128652, 0.0125714, 0.0122873,
      0.0120127, 0.011747,  0.01149,   0.0112412, 0.0110002}

Definition at line 67 of file HcalDeterministicFit.h.

siPM206Frac

constexpr float HcalDeterministicFit::siPM206Frac

staticconstexprprivate

Initial value:

= {
      0,         0,         0,         4.55043e-06, 0.00292008, 0.0147851, 0.0374319, 0.0688652, 0.105913,  0.145714,
      0.186153,  0.225892,  0.264379,  0.30145,     0.337074,   0.371247,  0.403973,  0.43526,   0.465115,  0.493554,
      0.520596,  0.546269,  0.570605,  0.59364,     0.615418,   0.635984,  0.655384,  0.673669,  0.690889,  0.707091,
      0.719418,  0.721882,  0.7127,    0.693915,    0.668746,   0.640097,  0.610129,  0.580218,  0.550958,  0.522548,
      0.495058,  0.468521,  0.442967,  0.418419,    0.394896,   0.37241,   0.350965,  0.330559,  0.31118,   0.292812,
      0.275432,  0.259013,  0.243523,  0.228928,    0.215193,   0.202281,  0.190154,  0.178773,  0.1681,    0.158098,
      0.148729,  0.139959,  0.131751,  0.124074,    0.116894,   0.110182,  0.103907,  0.0980423, 0.0925613, 0.0874393,
      0.0826525, 0.078179,  0.0739978, 0.0700894,   0.0664353,  0.0630185, 0.0598226, 0.0568328, 0.0540348, 0.0514156,
      0.0489628, 0.046665,  0.0445115, 0.0424924,   0.0405984,  0.038821,  0.037152,  0.0355841, 0.0341104, 0.0327243,
      0.0314199, 0.0301916, 0.0290343, 0.0279431,   0.0269136,  0.0259417, 0.0250235, 0.0241554, 0.0233341, 0.0225566,
      0.0218199, 0.0211214, 0.0204587, 0.0198294,   0.0192313,  0.0186626, 0.0181214, 0.0176059, 0.0171146, 0.016646,
      0.0161986, 0.0157713, 0.0153627, 0.0149719,   0.0145977,  0.0142393, 0.0138956, 0.0135659, 0.0132493, 0.0129451,
      0.0126528, 0.0123715, 0.0121007, 0.0118399,   0.0115885}

Definition at line 82 of file HcalDeterministicFit.h.

siPM207Frac

constexpr float HcalDeterministicFit::siPM207Frac

staticconstexprprivate

Initial value:

= {
      8.79768e-07, 1.0741e-05, 7.99305e-05, 0.000401863, 0.00148019, 0.00425688, 0.0100505, 0.0202671, 0.0360159,
      0.0578184,   0.0855186,  0.118382,    0.155299,    0.195005,   0.236264,   0.277979,  0.319254,  0.359407,
      0.397953,    0.434579,   0.469106,    0.50146,     0.531646,   0.559717,   0.585765,  0.609901,  0.632248,
      0.652923,    0.672002,   0.689414,    0.70477,     0.717237,   0.725596,   0.728535,  0.725028,  0.714632,
      0.697574,    0.674653,   0.647036,    0.616039,    0.58295,    0.548907,   0.514845,  0.481483,  0.449338,
      0.418754,    0.389935,   0.36298,     0.337908,    0.314684,   0.293236,   0.27347,   0.255279,  0.23855,
      0.22317,     0.209031,   0.196027,    0.184061,    0.173041,   0.162884,   0.153513,  0.144857,  0.136854,
      0.129445,    0.122577,   0.116205,    0.110283,    0.104775,   0.0996449,  0.0948609, 0.0903944, 0.0862196,
      0.0823128,   0.0786528,  0.0752202,   0.0719973,   0.0689682,  0.0661182,  0.0634341, 0.0609038, 0.058516,
      0.0562608,   0.0541287,  0.0521113,   0.0502007,   0.0483898,  0.046672,   0.0450411, 0.0434916, 0.0420183,
      0.0406164,   0.0392815,  0.0380094,   0.0367964,   0.035639,   0.0345339,  0.0334781, 0.0324687, 0.0315032,
      0.030579,    0.0296939,  0.0288458,   0.0280326,   0.0272526,  0.0265039,  0.025785,  0.0250943, 0.0244304,
      0.0237919,   0.0231777,  0.0225864,   0.022017,    0.0214684,  0.0209397,  0.0204298, 0.0199381, 0.0194635,
      0.0190053,   0.0185629,  0.0181354,   0.0177223,   0.0173229,  0.0169366,  0.0165628, 0.0162011}

Definition at line 97 of file HcalDeterministicFit.h.

tsWidth

constexpr int HcalDeterministicFit::tsWidth = 25

staticconstexprprivate

Definition at line 45 of file HcalDeterministicFit.h.