Use saturation to correctly label strip hits

Author: henrylay97

sbndcode/CRT/CRTReco/CRTStripHitProducer_module.cc

@@ -62,6 +62,7 @@ class sbnd::crt::CRTStripHitProducer : public art::EDProducer {
 
   std::string           fFEBDataModuleLabel;
   uint16_t              fADCThreshold;
+  uint16_t              fADCSaturation;
   std::vector<double>   fErrorCoeff;
   bool                  fAllowFlag1;
   bool                  fApplyTs0Window;
@@ -88,6 +89,7 @@ sbnd::crt::CRTStripHitProducer::CRTStripHitProducer(fhicl::ParameterSet const& p
   : EDProducer{p}
   , fFEBDataModuleLabel(p.get<std::string>("FEBDataModuleLabel"))
   , fADCThreshold(p.get<uint16_t>("ADCThreshold"))
+  , fADCSaturation(p.get<uint16_t>("ADCSaturation"))
   , fErrorCoeff(p.get<std::vector<double>>("ErrorCoeff"))
   , fAllowFlag1(p.get<bool>("AllowFlag1"))
   , fApplyTs0Window(p.get<bool>("ApplyTs0Window"))
@@ -290,7 +292,7 @@ std::vector<sbnd::crt::CRTStripHit> sbnd::crt::CRTStripHitProducer::CreateStripH
           if(pos - err < 0)
             err = pos;
 
-          stripHits.emplace_back(offline_channel_id, t0, t1, ref_time_s, pos, err, adc1, adc2);
+          stripHits.emplace_back(offline_channel_id, t0, t1, ref_time_s, pos, err, adc1, adc2, fADCSaturation);
         }
     }
 

sbndcode/CRT/CRTReco/crtrecoproducers_sbnd.fcl

@@ -6,6 +6,7 @@ crtstriphitproducer_sbnd:
 {
    FEBDataModuleLabel:        "crtsim"
    ADCThreshold:              60
+   ADCSaturation:             @local::sbnd_crtsim.DetSimParams.AdcSaturation
    ErrorCoeff:                [ 0.26, -0.27, 0.025 ]
    AllowFlag1:                false
    ApplyTs0Window:            false