Add in Joe Zennamo's idea of cutting on number of samples a certain distance from pedestal, both in raw digits and recob::wire

Author: tomjunk

sbndcode/Utilities/DigitalNoiseChannelStatus.h

@@ -47,6 +47,11 @@ class sbnd::DigitalNoiseChannelStatus {
   std::string fRawDigitLabel;
   std::string fRecobWireLabel;
 
+  int fNAwayFromPedestalRawDigit;
+  int fDistFromPedestalRawDigit;
+  int fNAwayFromPedestalRecobWire;
+  float fDistFromPedestalRecobWire;
+  
 };
 
 DECLARE_ART_SERVICE(sbnd::DigitalNoiseChannelStatus, LEGACY)

sbndcode/Utilities/DigitalNoiseChannelStatusDefaults.fcl

@@ -7,6 +7,10 @@ DigitalNoiseChannelStatusDefaults :
   RMSCutRawDigit: 100
   RMSCutWire: 100
   NBADCutRawDigit: 5
+  NAwayFromPedestalRawDigit: 100
+  DistFromPedestalRawDigit: 100
+  NAwayFromPedestalRecobWire: 100
+  DistFromPedestalRecobWire: 100
 }
 
 END_PROLOG

sbndcode/Utilities/DigitalNoiseChannelStatus_service.cc

@@ -13,7 +13,11 @@ sbnd::DigitalNoiseChannelStatus::DigitalNoiseChannelStatus(fhicl::ParameterSet c
   fNBADCutRawDigit = p.get<float>("NBADCutRawDigit",5);
   fRawDigitLabel = p.get<std::string>("RawDigitLabel","daq");
   fRecobWireLabel = p.get<std::string>("RecobWireLabel","caldata");
-
+  fNAwayFromPedestalRawDigit = p.get<int>("NAwayFromPedestalRawDigit",100);
+  fDistFromPedestalRawDigit = p.get<int>("DistFromPedestalRawDigit",100);
+  fNAwayFromPedestalRecobWire = p.get<int>("NAwayFromPedestalRecobWire",100);
+  fDistFromPedestalRecobWire = p.get<float>("DistFromPedestalRecobWire",100);
+  
   areg.sPreProcessEvent.watch(this, &sbnd::DigitalNoiseChannelStatus::pub_PrepEvent);
 }
 
@@ -52,15 +56,21 @@ void sbnd::DigitalNoiseChannelStatus::pub_PrepEvent(const art::Event& evt, art::
 	  rawadc.resize(rd.Samples());
 	  raw::Uncompress(rd.ADCs(), rawadc, rd.GetPedestal(), rd.Compression());
 	  bool alleven = true;
+	  int naway = 0;
 	  const short adc0 = rawadc.at(0);
 	  for (size_t i=0; i< rd.Samples(); ++i)
 	    {
 	      const short adc = rawadc.at(i);
 	      if (adc == 0xBAD) ++nhexbad;
-	      alleven &= ( ((adc - adc0) % 2) == 0 ); 
+	      alleven &= ( ((adc - adc0) % 2) == 0 );
+	      if (std::abs(adc - rd.GetPedestal()) > fDistFromPedestalRawDigit)
+		{
+		  ++naway;
+		}
 	    }
 	  chanbad |= (nhexbad > fNBADCutRawDigit);
 	  chanbad |= alleven;
+	  chanbad |= (naway > fNAwayFromPedestalRawDigit);
 	  if (chanbad) fDNChannels.emplace(rd.Channel());
 	}
     }
@@ -73,6 +83,16 @@ void sbnd::DigitalNoiseChannelStatus::pub_PrepEvent(const art::Event& evt, art::
 	  bool chanbad = false;
 	  auto rms = TMath::RMS(rw.NSignal(),rw.Signal().data());
 	  chanbad |= (rms > fRMSCutWire);
+	  auto mns = TMath::Mean(rw.NSignal(),rw.Signal().data());
+	  int naway = 0;
+	  for (size_t i=0; i<rw.NSignal(); ++i)
+	    {
+	      if (std::abs(rw.Signal()[i] - mns) > fDistFromPedestalRecobWire)
+		{
+		  naway ++;
+		}
+	    }
+	  chanbad |= (naway > fNAwayFromPedestalRecobWire);
 	  if (chanbad) fDNChannels.emplace(rw.Channel());
 	}
     }