HcalHF_PETalgorithm.h

Go to the documentation of this file.

#ifndef GUARD_HCALHF_PETALGORITHM_H
#define GUARD_HCALHF_PETALGORITHM_H 1

#include "DataFormats/HcalRecHit/interface/HFRecHit.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalCaloFlagLabels.h"

// Forward declarations
class HcalChannelQuality;
class HcalSeverityLevelComputer;

class HcalHF_PETalgorithm {
 public:
  HcalHF_PETalgorithm();
  
  HcalHF_PETalgorithm(std::vector<double> short_R, 
              std::vector<double> short_Energy, 
              std::vector<double> short_ET, 
              std::vector<double> long_R, 
              std::vector<double> long_Energy, 
              std::vector<double> long_ET,
              int HcalAcceptSeverityLevel,
              // special case for ieta=29
              std::vector<double> short_R_29,
              std::vector<double> long_R_29);

  // Destructor
  ~HcalHF_PETalgorithm();

  void HFSetFlagFromPET(HFRecHit& hf,
            HFRecHitCollection& rec,
            HcalChannelQuality* myqual,
            const HcalSeverityLevelComputer* mySeverity);
  double CalcThreshold(double abs_energy,std::vector<double> params);

  void SetShort_R(std::vector<double> x){short_R=x;}
  void SetShort_ET_Thresh(std::vector<double> x){short_ET_Thresh=x;}
  void SetShort_Energy_Thresh(std::vector<double> x){short_Energy_Thresh=x;}
  void SetLong_R(std::vector<double> x){long_R=x;}
  void SetLong_ET_Thresh(std::vector<double> x){long_ET_Thresh=x;}
  void SetLong_Energy_Thresh(std::vector<double> x){long_Energy_Thresh=x;}

  std::vector<double> GetShort_R(){return short_R;}
  std::vector<double> GetShort_ET_Thresh(){return short_ET_Thresh;}
  std::vector<double> GetShort_Energy_Thresh(){return short_Energy_Thresh;}
  std::vector<double> GetLong_R(){return long_R;}
  std::vector<double> GetLong_ET_Thresh(){return long_ET_Thresh;}
  std::vector<double> GetLong_Energy_Thresh(){return long_Energy_Thresh;}

  double bit(){return HcalCaloFlagLabels::HFLongShort;}

 private:
  std::vector<double> short_R;
  std::vector<double> short_ET_Thresh;
  std::vector<double> short_Energy_Thresh; 

  std::vector<double> long_R;
  std::vector<double> long_ET_Thresh;
  std::vector<double> long_Energy_Thresh;
  int HcalAcceptSeverityLevel_;
  std::vector<double> short_R_29;
  std::vector<double> long_R_29;
};


#endif