Inheritance diagram for filt::FinalStateParticleFilter:

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

bool	filter (art::Event &evt)

void	beginJob ()

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

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 ( )

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 )

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;
 }