ActiveDCMFilter_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       ActiveDCMFilter
// Plugin Type: filter (art v1_17_06)
// File:        ActiveDCMFilter_module.cc
//
// Generated at Tue Mar 27 03:39:21 2018 by Andrew S. Hirsch using cetskelgen
// from cetlib version v1_15_03.
////////////////////////////////////////////////////////////////////////

#include "art/Framework/Core/EDFilter.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"
#include "art/Framework/Principal/Run.h"
#include "art/Framework/Principal/SubRun.h"
#include "canvas/Utilities/InputTag.h"
#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

#include "DDTBaseDataProducts/DAQHit.h"
#include "DDTBaseDataProducts/HitList.h"

#include <memory>

namespace novaddt {
  class ActiveDCMFilter;
}


class novaddt::ActiveDCMFilter : public art::EDFilter {
public:
  explicit ActiveDCMFilter(fhicl::ParameterSet const & p);
  // The compiler-generated destructor is fine for non-base
  // classes without bare pointers or other resource use.

  // Plugins should not be copied or assigned.
  ActiveDCMFilter(ActiveDCMFilter const &) = delete;
  ActiveDCMFilter(ActiveDCMFilter &&) = delete;
  ActiveDCMFilter & operator = (ActiveDCMFilter const &) = delete;
  ActiveDCMFilter & operator = (ActiveDCMFilter &&) = delete;

  // Required functions.
  bool filter(art::Event & e) override;

private:

  // Declare member data here.
  art::InputTag fHitsTag; //where to read the hits for noise activity
  std::map<novaddt::DCM, uint64_t> fDCMOccupancy;
  size_t fNumDCMMinimum;

};


novaddt::ActiveDCMFilter::ActiveDCMFilter(fhicl::ParameterSet const & p):
fHitsTag(p.get<std::string>("hits_tag")),
fNumDCMMinimum(p.get<int>("num_DCM_minimum"))
{
    mf::LogInfo("ActiveDCMFilter")<<"Configured:\n"
        <<"Reading hits from: "<<fHitsTag<<"\n "
        <<"Discard events with: N active DCMs < "<<fNumDCMMinimum;
}

bool novaddt::ActiveDCMFilter::filter(art::Event & e)
{
    //clear map of DCM activities
    fDCMOccupancy.clear();
    //fill the map with hits
    auto hits=e.getValidHandle<novaddt::HitList>(fHitsTag);
    for(const auto &hit: *hits) fDCMOccupancy[hit]+=1;
    
    LOG_DEBUG("Active DCMs")<<fDCMOccupancy.size();
    //filter out if number of active DCMs is less than a given Minimum
    return (fDCMOccupancy.size() >= fNumDCMMinimum);
}

DEFINE_ART_MODULE(novaddt::ActiveDCMFilter)