HcalHF_S9S1algorithm.h

Go to the documentation of this file.

#ifndef GUARD_HCALHF_S9S1ALGORITHM_H
#define GUARD_HCALHF_S9S1ALGORITHM_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_S9S1algorithm {
 public:
  HcalHF_S9S1algorithm();
  
  HcalHF_S9S1algorithm(const std::vector<double>& short_optimumSlope, 
               const std::vector<double>& short_Energy, 
               const std::vector<double>& short_ET, 
               const std::vector<double>& long_optimumSlope, 
               const std::vector<double>& long_Energy, 
               const std::vector<double>& long_ET,
               int HcalAcceptSeverityLevel,
               bool isS8S1);

  // Destructor
  ~HcalHF_S9S1algorithm();

  void HFSetFlagFromS9S1(HFRecHit& hf,
            HFRecHitCollection& rec,
            const HcalChannelQuality* myqual,
            const HcalSeverityLevelComputer* mySeverity);
  double CalcSlope(int abs_ieta, const std::vector<double>& params);
  double CalcEnergyThreshold(double abs_energy,const std::vector<double>& params);

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

 private:

  std::vector<double> short_ET_;
  std::vector<double> short_Energy_; 
  std::vector<double> long_ET_;
  std::vector<double> long_Energy_;

  std::vector<double> LongSlopes;
  std::vector<double> ShortSlopes;
  std::vector<double> LongEnergyThreshold;
  std::vector<double> ShortEnergyThreshold;
  std::vector<double> LongETThreshold;
  std::vector<double> ShortETThreshold;
  int HcalAcceptSeverityLevel_;
  bool isS8S1_;
};


#endif