HcalReconConditions.h

 
#ifndef HCAL_HCALRECONCONDTIONS_H_
#define HCAL_HCALRECONCONDTIONS_H_
 
#include "Conditions/SimpleTableCondition.h"
#include "DetDescr/HcalDigiID.h"
 
namespace hcal {
 
class HcalReconConditions : public framework::ConditionsObject {
 public:
  static const unsigned int i_m_adc_i = 0;
  static const unsigned int i_cut_point_tot = 1;
  static const unsigned int i_high_slope = 2;
  static const unsigned int i_high_offset = 3;
  static const unsigned int i_low_slope = 4;
  static const unsigned int i_low_power = 5;
  static const unsigned int i_lower_offset = 6;
  static const unsigned int i_tot_not = 7;
  static const unsigned int i_channel = 8;
  static const unsigned int i_flagged = 9;
  static const std::string CONDITIONS_NAME;
 
  HcalReconConditions(const conditions::DoubleTableCondition& adc_ped,
                      const conditions::DoubleTableCondition& adc_gain,
                      const conditions::DoubleTableCondition& tot_calib,
                      const conditions::DoubleTableCondition& toa_calib);
 
  double adcPedestal(const ldmx::HcalDigiID& id, int idx = 0) const {
    return adc_pedestals_.get(id.raw(), idx);
  }
 
  double adcGain(const ldmx::HcalDigiID& id, int idx = 0) const {
    return adc_gains_.get(id.raw(), idx);
  }
 
  bool is_adc(const ldmx::HcalDigiID& id, double sum_tot) const;
 
  double linearize(const ldmx::HcalDigiID& id, double sum_tot) const;
 
  double totCalib(const ldmx::HcalDigiID& id, int idx = 0) const {
    return tot_calibs_.get(id.raw(), idx);
  }
 
  double toaCalib(const ldmx::HcalDigiID& id, int idx = 0) const {
    return toa_calibs_.get(id.raw(), idx);
  }
 
 private:
  const conditions::DoubleTableCondition& adc_pedestals_;
  const conditions::DoubleTableCondition& adc_gains_;
  const conditions::DoubleTableCondition& tot_calibs_;
  const conditions::DoubleTableCondition& toa_calibs_;
};  // HcalReconConditions
 
}  // namespace hcal
 
#endif  // HCAL_HCALRECONCONDITIONS_H_