ADCFilter_module.cc

Go to the documentation of this file.

//
// ADCFilter class:
// Algorithm to ignore events with no ADC values
// above user-defined threshold.
//
// msoderbe@syr.edu
//

#include <algorithm>

//Framework Includes
#include "art/Framework/Core/EDFilter.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/View.h"
#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "fhiclcpp/ParameterSet.h"

//Larsoft Includes
#include "lardataobj/RawData/RawDigit.h"
#include "lardataobj/RawData/raw.h"
#include "larevt/CalibrationDBI/Interface/ChannelStatusProvider.h"
#include "larevt/CalibrationDBI/Interface/ChannelStatusService.h"

namespace filter {

  class ADCFilter : public art::EDFilter {

  public:
    explicit ADCFilter(fhicl::ParameterSet const&);

  private:
    bool filter(art::Event& evt);

    std::string fDigitModuleLabel;
    double fMinADC;
  }; // class ADCFilter

  //-------------------------------------------------
  ADCFilter::ADCFilter(fhicl::ParameterSet const& pset) : EDFilter{pset}
  {
    fDigitModuleLabel = pset.get<std::string>("DigitModuleLabel");
    fMinADC = pset.get<double>("MinADC");
  }

  //-------------------------------------------------
  bool ADCFilter::filter(art::Event& evt)
  {
    //Read in raw data
    art::View<raw::RawDigit> rawdigitView;
    evt.getView(fDigitModuleLabel, rawdigitView);

    if (!rawdigitView.size()) return false;

    lariov::ChannelStatusProvider const& channelFilter =
      art::ServiceHandle<lariov::ChannelStatusService const>()->GetProvider();

    // look through the good channels
    //      for(const raw::RawDigit* digit: filter::SelectGoodChannels(rawdigitView))
    for (const raw::RawDigit* digit : rawdigitView) {
      if (!channelFilter.IsGood(digit->Channel())) continue;
      //get ADC values after decompressing
      std::vector<short> rawadc(digit->Samples());
      raw::Uncompress(digit->ADCs(), rawadc, digit->Compression());
      short max = *std::max_element(rawadc.begin(), rawadc.end()) - digit->GetPedestal();
      if (max >= fMinADC) return true; //found one ADC value above threshold, pass filter
    }

    return false; //didn't find ADC above threshold
  }

  DEFINE_ART_MODULE(ADCFilter)

} //namespace filter