EventFilter_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       EventFilter
// Module Type: filter
// File:        EventFilter_module.cc
//
// Generated at Wed Jan 28 16:48:29 2015 by Matthew Tamsett using artmod
// from cetpkgsupport v1_07_00.
////////////////////////////////////////////////////////////////////////

#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"
//// this was working at least up to S17-02-21     
//#include "art/Utilities/InputTag.h"
#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

#include <memory>
#include <iostream>
#include <fstream>
#include <string>


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

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

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

private:

  // Declare member data here.
  std::ifstream eventList;
  std::string _InputFile;

  std::vector<art::SubRunNumber_t> run;
  std::vector<art::RunNumber_t> subrun;
  std::vector<art::EventNumber_t> event;
};
//----------------------------------------------------------------------
EventFilter::EventFilter(fhicl::ParameterSet const & p)
  : _InputFile  (p.get<std::string >("InputFile"  ))
// :
// Initialize member data here.
{
  // Call appropriate produces<>() functions here.
}
//----------------------------------------------------------------------
void EventFilter::beginJob()
{
  run.clear();
  subrun.clear();
  event.clear();

  eventList.open(_InputFile.c_str());
  int fRun, fSubrun, fEvent, fSlice;
  double slicetime;

  while(!eventList.eof())
    {
      // Change this once the release is updated
      eventList >> fRun >> fSubrun >> fEvent >> fSlice >> slicetime;
      //eventList >> fRun >> fSubrun >> fEvent >> fSlice;
      run.push_back(fRun);
      subrun.push_back(fSubrun);
      event.push_back(fEvent);
    }
  eventList.close();
}

bool EventFilter::filter(art::Event& evt)
{

  if(std::find(run.begin(),    run.end(),    evt.id().run())    != run.end() &&
     std::find(subrun.begin(), subrun.end(), evt.id().subRun()) != subrun.end() &&
     std::find(event.begin(),  event.end(),  evt.id().event())  != event.end())
    return true;
  
  else return false;
}
//----------------------------------------------------------------------
DEFINE_ART_MODULE(EventFilter)