Inheritance diagram for filt::FinalStateParticleFilter:

Public Member Functions
	FinalStateParticleFilter (fhicl::ParameterSet const &)

bool	filter (art::Event &evt)

void	beginJob ()

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 Member Functions
bool	isSubset (std::vector< int > const &a, std::vector< int > const &b) const

Private Attributes
std::string	fGenieModuleLabel

std::vector< int >	fPDG

std::vector< int >	fStatusCode

TH1D *	fSelectedEvents

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 30 of file FinalStateParticleFilter_module.cc.

Constructor & Destructor Documentation

filt::FinalStateParticleFilter::FinalStateParticleFilter ( fhicl::ParameterSet const & pset )

explicit

Definition at line 57 of file FinalStateParticleFilter_module.cc.

     : EDFilter{pset}
   {
     fGenieModuleLabel = pset.get< std::string      >("GenieModuleLabel");
     fPDG              = pset.get< std::vector<int> >("PDG");
   }

Member Function Documentation

void filt::FinalStateParticleFilter::beginJob ( )

virtual

Reimplemented from art::EDFilter.

Definition at line 65 of file FinalStateParticleFilter_module.cc.

   {
     art::ServiceHandle<art::TFileService const> 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
   }

bool filt::FinalStateParticleFilter::filter ( art::Event & evt )

virtual

Implements art::EDFilter.

Definition at line 73 of file FinalStateParticleFilter_module.cc.

   {
 
     //const TDatabasePDG* databasePDG = TDatabasePDG::Instance();
 
     art::Handle< std::vector<simb::MCTruth> > mclist;
     evt.getByLabel(fGenieModuleLabel,mclist);
     art::Ptr<simb::MCTruth> mc(mclist,0);
 
     fTotalEvents->Fill(1);
 
     std::vector<int> finalstateparticles;
 
     //get a vector of final state particles
     for(int i = 0; i < mc->NParticles(); ++i){
       simb::MCParticle part(mc->GetParticle(i));
       if(part.StatusCode()== 1)
         finalstateparticles.push_back(part.PdgCode());
     }
 
     if(isSubset(fPDG, finalstateparticles)){
       fSelectedEvents->Fill(1);
       std::cout << "this is a selected event" << std::endl;
     }
 
     return isSubset(fPDG, finalstateparticles); // returns true if the user-defined fPDG exist(s) in the final state particles
 
   } // bool

bool filt::FinalStateParticleFilter::isSubset	(	std::vector< int > const &	a,
		std::vector< int > const &	b
	)		const

private

Definition at line 105 of file FinalStateParticleFilter_module.cc.

 {
   for (auto const a_int : a) {
     bool found = false;
     for (auto const b_int : b) {
       if (a_int == b_int){
         found = true;
         break;
       }
     }
 
     if (!found){
       return false;
     }
   }
   return true;
 }