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(const std::vector<double>& short_R, 
              const std::vector<double>& short_Energy, 
              const std::vector<double>& short_ET, 
              const std::vector<double>& long_R, 
              const std::vector<double>& long_Energy, 
              const std::vector<double>& long_ET,
              int HcalAcceptSeverityLevel,
              // special case for ieta=29
              const std::vector<double>& short_R_29,
              const std::vector<double>& long_R_29);

  // Destructor
  ~HcalHF_PETalgorithm();

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

  void SetShort_R(const std::vector<double>& x){short_R=x;}
  void SetShort_ET_Thresh(const std::vector<double>& x){short_ET_Thresh=x;}
  void SetShort_Energy_Thresh(const std::vector<double>& x){short_Energy_Thresh=x;}
  void SetLong_R(const std::vector<double>& x){long_R=x;}
  void SetLong_ET_Thresh(const std::vector<double>& x){long_ET_Thresh=x;}
  void SetLong_Energy_Thresh(const 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