TrigEnergySum.h

#ifndef TRIGGER_EVENT_TRIGENERGYSUM_H_
#define TRIGGER_EVENT_TRIGENERGYSUM_H_
 
// ldmx-sw
#include <stdint.h>  //uint32_t
 
// ROOT
#include "TObject.h"  //For ClassDef
 
namespace trigger {
 
// Forward declaration needed by typedef
class TrigEnergySum;
 
typedef std::vector<TrigEnergySum> TrigEnergySumCollection;
 
class TrigEnergySum {
 public:
  TrigEnergySum() = default;
 
  TrigEnergySum(int layer, int hwEnergy = 0);
  TrigEnergySum(int layer, int module, float energy);
 
  virtual ~TrigEnergySum() = default;
 
  bool operator<(const TrigEnergySum &sum) { return hwEnergy_ < sum.hwEnergy_; }
 
  void Clear() {
    energy_ = 0;
    hwEnergy_ = 0;
    layer_ = 0;
    module_ = 0;
  }
 
  void setLayer(int layer) { layer_ = layer; }
  int layer() const { return layer_; }
 
  void setModule(int module) { module_ = module; }
  int module() const { return module_; }
 
  void setEnergy(float energy) { energy_ = energy; }
  float energy() const { return energy_; }
 
  void setHwEnergy(int hwEnergy) { hwEnergy_ = hwEnergy; }
  int hwEnergy() const { return hwEnergy_; }
  // inline float pe() const { return hwEnergy_*pe_per_adc_; }
  // inline float hadEnergy() const { return
  // hwEnergy_*pe_per_adc_*mev_per_pe_*had_sample_frac_; }
 
  void Print() const;
 
  friend std::ostream &operator<<(std::ostream &o, const TrigEnergySum &d);
 
  friend std::ostream &operator<<(std::ostream &o,
                                  const TrigEnergySumCollection &c);
 
 private:
  int layer_{0};
  int module_{0};
  int hwEnergy_{0};
  float energy_{0};
  /* const float pe_per_adc_{1.2/5}; // gain * pe/mV */
  /* const float mev_per_pe_{4.66/68}; // MeV/PE (MIP) */
  /* const float had_sample_frac_{0.109813*0.818731}; */
  // const float had_samp_frac = (20/77.07)/(20/77.07 + 25/16.77); // 0.148266
  // const float em_samp_frac = (20/41.31)/(20/41.31 + 25/1.757); // 0.032906
  // const float samp_frac = (em_samp_frac + 2*had_samp_frac)/3; // 0.109813
  // const float attn = exp(-1/5.); // 0.818731 5m attenuation length, 1m
  // half-bar
 
  ClassDef(TrigEnergySum, 1);
};
}  // namespace trigger
 
#endif  // TRIGGER_EVENT_CALOTRIGPRIM_H_