#include <BCToEFilterAlgo.h>

Public Member Functions
	BCToEFilterAlgo (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

bool	filter (const edm::Event &iEvent)

bool	hasBCAncestors (const reco::GenParticle &gp)

	~BCToEFilterAlgo ()

Private Member Functions
bool	isBCBaryon (const reco::GenParticle &gp)

bool	isBCHadron (const reco::GenParticle &gp)

bool	isBCMeson (const reco::GenParticle &gp)

Private Attributes
float	eTThreshold_

float	FILTER_ETA_MAX_

edm::EDGetTokenT< reco::GenParticleCollection >	genParSource_

Detailed Description

BCToEFilterAlgo returns true for events that have an electron, above configurable eT threshold and within |eta|<2.5, that has an ancestor of a b or c quark

Author: J Lamb, UCSB

Definition at line 26 of file BCToEFilterAlgo.h.

Constructor & Destructor Documentation

BCToEFilterAlgo::BCToEFilterAlgo	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 7 of file BCToEFilterAlgo.cc.

References eTThreshold_, FILTER_ETA_MAX_, objects.autophobj::float, genParSource_, and edm::ParameterSet::getParameter().

                                                                                            {
 
   //set constants
   FILTER_ETA_MAX_=2.5;
   eTThreshold_=(float)iConfig.getParameter<double>("eTThreshold");
   genParSource_=iC.consumes<reco::GenParticleCollection>(iConfig.getParameter<edm::InputTag>("genParSource"));
 
 }

BCToEFilterAlgo::~BCToEFilterAlgo ( )

Definition at line 16 of file BCToEFilterAlgo.cc.

16 {

17 }

Member Function Documentation

bool BCToEFilterAlgo::filter ( const edm::Event & iEvent )

Definition at line 22 of file BCToEFilterAlgo.cc.

References funct::abs(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), eTThreshold_, FILTER_ETA_MAX_, genParSource_, edm::Event::getByToken(), runTauDisplay::gp, hasBCAncestors(), reco::LeafCandidate::pdgId(), mps_fire::result, and reco::LeafCandidate::status().

                                                     {
 
   bool result=false;
 
   
   
   edm::Handle<reco::GenParticleCollection> genParsHandle;
   iEvent.getByToken(genParSource_,genParsHandle);
   reco::GenParticleCollection genPars=*genParsHandle;
 
   for (uint32_t ig=0;ig<genPars.size();ig++) {
     reco::GenParticle gp=genPars.at(ig);
     if (gp.status()==1 && abs(gp.pdgId())==11 && gp.et()>eTThreshold_ && fabs(gp.eta())<FILTER_ETA_MAX_) {
       if (hasBCAncestors(gp)) {
     result=true;
       }
     }
   }
   return result;
 
 }

bool BCToEFilterAlgo::hasBCAncestors ( const reco::GenParticle & gp )

Definition at line 49 of file BCToEFilterAlgo.cc.

References isBCHadron(), reco::CompositeRefCandidateT< D >::motherRef(), and reco::CompositeRefCandidateT< D >::numberOfMothers().

Referenced by filter().

                                                               {
 
   //stopping condition: this particle is a b or c hadron
   if (isBCHadron(gp)) return true;
   //stopping condition: this particle has no mothers
   if (gp.numberOfMothers()==0) return false;
   //otherwise continue
   bool retval=false;
   for (uint32_t im=0;im<gp.numberOfMothers();im++) {
     retval=retval || hasBCAncestors(*gp.motherRef(im));
   }
   return retval;
   
 }

bool BCToEFilterAlgo::isBCBaryon ( const reco::GenParticle & gp )

private

Definition at line 80 of file BCToEFilterAlgo.cc.

References funct::abs(), BPhysicsValidation_cfi::pdgid, and reco::LeafCandidate::pdgId().

Referenced by isBCHadron().

                                                           {
   
   uint32_t pdgid=abs(gp.pdgId());
   uint32_t thousands=pdgid%10000;
   if (thousands>=4000 && thousands <6000) {
     return true;
   } else {
     return false;
   }
 
 }

bool BCToEFilterAlgo::isBCHadron ( const reco::GenParticle & gp )

private

Definition at line 64 of file BCToEFilterAlgo.cc.

References isBCBaryon(), and isBCMeson().

Referenced by hasBCAncestors().

                                                           {
   return isBCMeson(gp) || isBCBaryon(gp);
 }

bool BCToEFilterAlgo::isBCMeson ( const reco::GenParticle & gp )

private

Definition at line 68 of file BCToEFilterAlgo.cc.

References funct::abs(), BPhysicsValidation_cfi::pdgid, and reco::LeafCandidate::pdgId().

Referenced by isBCHadron().

                                                          {
   
   uint32_t pdgid=abs(gp.pdgId());
   uint32_t hundreds=pdgid%1000;
   if (hundreds>=400 && hundreds<600) {
     return true;
   } else {
     return false;
   }
 
 }