/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/RecoLocalCalo/HcalRecAlgos/interface/HcalHFStatusBitFromDigis.h

Go to the documentation of this file.

#ifndef HCALHFSTATUSFROMDIGIS_H
#define HCALHFSTATUSFROMDIGIS_H 1

#include "DataFormats/HcalDigi/interface/HFDataFrame.h"
#include "DataFormats/HcalRecHit/interface/HFRecHit.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalCaloFlagLabels.h"
#include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CalibFormats/HcalObjects/interface/HcalCoderDb.h"

class HcalHFStatusBitFromDigis {
public:
  HcalHFStatusBitFromDigis();
  HcalHFStatusBitFromDigis(int recoFirstSample,
                           int recoSamplesToAdd,
                           const edm::ParameterSet& HFDigiTimeParams,
                           const edm::ParameterSet& HFTimeInWindowParams);
  // Destructor
  ~HcalHFStatusBitFromDigis();

  // The important stuff!  Methods for setting the status flag values
  void hfSetFlagFromDigi(HFRecHit& hf, const HFDataFrame& digi,
                         const HcalCoder& coder,
                         const HcalCalibrations& calib);
  // Hack for 3_11 processing of HF in both old, new runs
  void resetTimeSamples(int firstSample, int samplesToAdd){firstSample_=firstSample; samplesToAdd_=samplesToAdd;}

private:

  // variables for cfg files

  // VARIABLES FOR SETTING HFDigiTime FLAG
  double minthreshold_;
  // Reco Window
  int recoFirstSample_;
  int recoSamplesToAdd_;
  // Special window for Igor's algorithm (not necessarily the same as reco window)
  int firstSample_;
  int samplesToAdd_;
  int expectedPeak_;
  
  // Coefficients used to determine energy ratio threshold:
  // E_peak/(Etotal) > coef0_-exp(coef1_+coef2_*Energy)
  double coef0_;
  double coef1_;
  double coef2_;

  // VARIABLES FOR SETTING HFInTimeWindow FLAG
  double HFlongwindowEthresh_;
  std::vector<double> HFlongwindowMinTime_;
  std::vector<double> HFlongwindowMaxTime_;
  double HFshortwindowEthresh_;
  std::vector<double> HFshortwindowMinTime_;
  std::vector<double> HFshortwindowMaxTime_;
};

#endif