HcalFlagHFDigiTimeParam.h

Go to the documentation of this file.

#ifndef GUARD_HcalFlagHFDigiTimeParam
#define GUARD_HcalFlagHFDigiTimeParam

#include <boost/cstdint.hpp>
#include <vector>

class HcalFlagHFDigiTimeParam{

  // Default constructor sets parameters according to 6-TS digis
  // (Is there a way to pass in default vector arguments directly, rather than by using push_back?) 
 public:HcalFlagHFDigiTimeParam():mId(0),mHFdigiflagFirstSample(1),mHFdigiflagSamplesToAdd(3),mHFdigiflagExpectedPeak(2),mHFdigiflagMinEthreshold(40) {
    mHFdigiflagCoefficients.clear();
    mHFdigiflagCoefficients.push_back(0.93);
    mHFdigiflagCoefficients.push_back(-0.38275);
    mHFdigiflagCoefficients.push_back(-0.012667);
  }
  
  // Set parameters based on input
  HcalFlagHFDigiTimeParam(unsigned long fId, 
              unsigned int fFirstSample, 
              unsigned int fSamplesToAdd, 
              unsigned int fExpectedPeak, 
              double fminEThreshold, 
              const std::vector<double>& fcoef): 
    mId(fId), 
    mHFdigiflagFirstSample(fFirstSample), 
    mHFdigiflagSamplesToAdd(fSamplesToAdd), 
    mHFdigiflagExpectedPeak(fExpectedPeak), 
    mHFdigiflagMinEthreshold(fminEThreshold),
    mHFdigiflagCoefficients(fcoef)
      {      }
    
    uint32_t rawId () const {return mId;}
  

    // Return flag parameter values
    uint32_t HFdigiflagFirstSample()   const  {return mHFdigiflagFirstSample;}
    uint32_t HFdigiflagSamplesToAdd()  const  {return mHFdigiflagSamplesToAdd;}
    uint32_t HFdigiflagExpectedPeak()  const  {return mHFdigiflagExpectedPeak;}
    double    HFdigiflagMinEThreshold() const  {return mHFdigiflagMinEthreshold;}
    std::vector<double> HFdigiflagCoefficients()  const {return mHFdigiflagCoefficients;}

 private:
    uint32_t mId; // detector ID
    uint32_t mHFdigiflagFirstSample;         // first sample used in NTS calculation
    uint32_t mHFdigiflagSamplesToAdd;        // # of sampels to use in NTS calculation
    uint32_t mHFdigiflagExpectedPeak;        // expected peak position; used for calculating TS(peak)
    double    mHFdigiflagMinEthreshold;       // minimum energy for flagged rechit
    std::vector<double> mHFdigiflagCoefficients; // coefficients used to parameterize TS(peak)/NTS threshold:  [0]-exp([1]+[2]*E+....)
};

#endif