HcalDeterministicFit.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_HcalRecAlgos_HcalDeterministicFit_h
#define RecoLocalCalo_HcalRecAlgos_HcalDeterministicFit_h

#include <typeinfo>
#include <vector>

#include "CalibCalorimetry/HcalAlgos/interface/HcalTimeSlew.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/PedestalSub.h"

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/HcalRecHit/interface/HBHEChannelInfo.h"
#include "CalibCalorimetry/HcalAlgos/interface/HcalPulseShapes.h"
#include "CalibFormats/HcalObjects/interface/HcalCoder.h"
#include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"

class HcalDeterministicFit {
 public:
  HcalDeterministicFit();
  ~HcalDeterministicFit();

  enum FType {shapeLandau, shape205, shape206, shape207};

  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;

  // This is the CMSSW Implementation of the apply function
  void getLandauFrac(float tStart, float tEnd, float &sum) const;
  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;

 private:
  HcalTimeSlew::ParaSource fTimeSlew_;
  HcalTimeSlew::BiasSetting fTimeSlewBias_;
  PedestalSub fPedestalSubFxn_;
  bool applyTimeSlew_;
  double frespCorr_;
 
  static constexpr int HcalRegion[2] = {16, 17};
  static constexpr int tsWidth = 25;
  static constexpr float negThresh[2] = {-3., 15.};
  static constexpr float invGpar[3] = {-13.11, 11.29, 5.133};
  static constexpr float rCorr[2] = {0.95, 0.95};
  static constexpr float rCorrSiPM[2] = {1., 1.};

  static constexpr float landauFrac[125] = { 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 };

  static constexpr float siPM205Frac[125] = {
    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 };

  static constexpr float siPM206Frac[125] = {
    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};

  static constexpr float siPM207Frac[125] = {
    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};
  
};

#endif // HLTAnalyzer_h