Inheritance diagram for protoana::ProtoDUNETruthBeamFilter:

Public Member Functions
	ProtoDUNETruthBeamFilter (fhicl::ParameterSet const &pset)

virtual	~ProtoDUNETruthBeamFilter ()

void	beginJob ()

bool	filter (art::Event &evt)

void	reconfigure (fhicl::ParameterSet const &pset)

void	endJob ()

Public Member Functions inherited from art::EDFilter
	EDFilter (fhicl::ParameterSet const &pset)

template<typename Config >
	EDFilter (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Filter
virtual	~Filter () noexcept

	Filter (fhicl::ParameterSet const &)

	Filter (Filter const &)=delete

	Filter (Filter &&)=delete

Filter &	operator= (Filter const &)=delete

Filter &	operator= (Filter &&)=delete

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
bool	fDebug

std::vector< int >	fPDG
	Controls the verbosity. More...

TH1D *	fSelectedEvents
	List of particle PDGs we want to keep. More...

TH1D *	fTotalEvents

Additional Inherited Members
Public Types inherited from art::EDFilter
using	ModuleType = EDFilter

using	WorkerType = WorkerT< EDFilter >

Public Types inherited from art::detail::Filter
template<typename UserConfig >
using	Table = Modifier::Table< UserConfig >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Attributes inherited from art::detail::Filter
static constexpr bool	Pass {true}

static constexpr bool	Fail {false}

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 45 of file ProtoDUNETruthBeamFilter_module.cc.

Constructor & Destructor Documentation

protoana::ProtoDUNETruthBeamFilter::ProtoDUNETruthBeamFilter ( fhicl::ParameterSet const & pset )

explicit

Definition at line 72 of file ProtoDUNETruthBeamFilter_module.cc.

 : EDFilter(pset)
 {
   this->reconfigure(pset);
 
 }

protoana::ProtoDUNETruthBeamFilter::~ProtoDUNETruthBeamFilter ( )

virtual

Definition at line 80 of file ProtoDUNETruthBeamFilter_module.cc.

80 {}

Member Function Documentation

void protoana::ProtoDUNETruthBeamFilter::beginJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 83 of file ProtoDUNETruthBeamFilter_module.cc.

                                                 {
 
         art::ServiceHandle<art::TFileService> tfs;
     fSelectedEvents = tfs->make<TH1D>("fSelectedEvents", "Number of Selected Events", 3, 0, 3); //counts the number of selected events 
     fTotalEvents = tfs->make<TH1D>("fTotalEvents", "Total Events", 3, 0, 3); //counts the initial number of events in the unfiltered root input file
   
 }

void protoana::ProtoDUNETruthBeamFilter::endJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 147 of file ProtoDUNETruthBeamFilter_module.cc.

147 {}

bool protoana::ProtoDUNETruthBeamFilter::filter ( art::Event & evt )

virtual

Implements art::EDFilter.

Definition at line 101 of file ProtoDUNETruthBeamFilter_module.cc.

                                                           {
 
 
 
   if(fDebug) std::cout << "Reading Event" << std::endl;
     
 
   bool found = false;  
   // Access truth info
   std::vector< art::Ptr<simb::MCTruth> > mclist;
   auto mctruthhandle = evt.getHandle< std::vector<simb::MCTruth> >("generator");
   if (mctruthhandle){
     art::fill_ptr_vector(mclist,mctruthhandle);
   }
   else{
     mf::LogError("protoana::ProtoDUNETruthBeamParticle::analyze") << "Requested protoDUNE beam generator information does not exist!";
     return found;
   }
   fTotalEvents->Fill(1); //count total events
 
   art::Ptr<simb::MCTruth> mctruth;
   mctruth = mclist[0];
   
 
   for(int i = 0; i < mctruth->NParticles(); ++i){
     const simb::MCParticle& part(mctruth->GetParticle(i));
   
           bool isgoodparticle   = (part.Process() == "primary" && mctruth->Origin()==4); //primary particle and beam origin
           if(!isgoodparticle) continue;
     // Select partciles with the chosen pdg
     for(unsigned int j = 0; j < fPDG.size(); j++){
                   if(part.PdgCode() == fPDG[j]){
                   found = true;
                   if(fDebug) std::cout << "this is a selected event with PDG "<< part.PdgCode() << std::endl;
                   fSelectedEvents->Fill(1); //count selected events
                     break;
             }
     }
     
   }   
   
   return found;
   
 }

void protoana::ProtoDUNETruthBeamFilter::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 92 of file ProtoDUNETruthBeamFilter_module.cc.

                                                                                {
 
   fPDG             = pset.get< std::vector<int> >("PDG");
   fDebug             = pset.get<bool>("IsVerbose"); 
 
 }