TestBeamHit.h

Go to the documentation of this file.

 
#ifndef TRIGSCINT_EVENT_TESTBEAMHIT_H
#define TRIGSCINT_EVENT_TESTBEAMHIT_H
 
/*~~~~~~~~~~~*/
/*   Event   */
/*~~~~~~~~~~~*/
#include "TrigScint/Event/TrigScintHit.h"
 
namespace trigscint {
 
class TestBeamHit : public ldmx::TrigScintHit {
 public:
  TestBeamHit() = default;
 
  ~TestBeamHit() = default;
 
  void Clear(Option_t *option = "");
 
  void Print(Option_t *option = "") const;
 
  void setPedestal(const float pedestal) { pedestal_ = pedestal; };
 
 
  float getEarlyPedestal() const { return earlyPedestal_; }
 
  void setEarlyPedestal(const float earlyPed) { earlyPedestal_ = earlyPed; };
 
 
  float getPedestal() const { return pedestal_; }
 
  void setQ(const float q) { pulseQ_ = q; }
 
  float getQ() const { return pulseQ_; }
 
  //  void setPE(const float pe) { PE_ = pe; }
 
  // float getPE() const { return PE_; }
 
  void setStartSample(const int startSample) { startSample_ = startSample; };
 
 
  int getStartSample() const { return startSample_; }
 
  void setSampAbovePed(const int sampAbovePed) {
    sampAbovePed_ = sampAbovePed;
  };
 
  int getSampAbovePed() const { return sampAbovePed_; }
 
  void setSampAboveThr(const int sampAboveThr) {
    sampAboveThr_ = sampAboveThr;
  };
 
  int getSampAboveThr() const { return sampAboveThr_; }
 
  void setHitQuality(const int isClean) { passHitQuality_ = isClean; };
 
  int getHitQuality() const { return passHitQuality_; }
 
  void setQualityFlag(const uint flag) { this->flag_ = flag; };
 
 
  float getQualityFlag() const { return flag_; }
 
  void setPulseWidth(const int pulseWidth) { pulseWidth_ = pulseWidth; };
 
  int getPulseWidth() const { return pulseWidth_; }
 
  bool operator<(const TestBeamHit &rhs) const {
    return this->pulseQ_ < rhs.pulseQ_;
  }
 
 private:
  float pedestal_{
      -999.};  // assumed/average channel pedestal used in subtraction [fC]
  float earlyPedestal_{
      -999.};  // early pedestal (average over first 5 time samples)   [fC]
  float pulseQ_{-999.};  // integrated, ped subtracted charge over pulse [fC]
  // float PE_{-1.};       // number of photoelectrons in the hit (override from
  // trigscinthit which seems unreliable)
  int startSample_{-1};   // start sample
  int pulseWidth_{-1};    // specified pulse width (in number of samples)
  int sampAbovePed_{-1};  // actual number of consecutive samples above pedestal
                          // after start sample
  int sampAboveThr_{-1};  // actual number of consecutive samples above
                          // threshold after start sample
  //  bool passHitQuality_{false};  //track if hit has been checked for and
  //  passed quality criteria
  int passHitQuality_{
      0};  // track if hit has been checked for and passed quality criteria
 
  uint flag_{0};  // more elaborate quality flag (binary sum)
 
  //        - nPulses -- TODO, when running in continuous mode
  /*
  - total Q // this is an event readout property
  - ped subtracted total Q // so is this
  - PE value // already property of trigscinthit
  - max amplitude in pulse // so is this
  - store material assumption? isLYSO? -- no, this is a bar property, not a hit
  property
  */
 
  ClassDef(TestBeamHit, 3);
 
};  // TestBeamHit
 
}  // namespace trigscint
 
#endif  // TRIGSCINT_EVENT_TESTBEAMHIT_H