HFFilter Class Reference

#include <PhysicsTools/HFFilter/src/HFFilter.cc>

Inheritance diagram for HFFilter:

Public Member Functions
virtual void	endJob ()
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	HFFilter (const edm::ParameterSet &)
	~HFFilter ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Attributes
double	etaMax_
edm::InputTag	genJetsCollName_
double	ptMin_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Description: Filter to see if there are heavy flavor GenJets in this event

Implementation: The implementation is simple, it loops over the GenJets and checks if any constituents have a pdg ID that matches a list. It also has a switch to count objects from a gluon parent, so the user can turn off counting gluon splitting.

Definition at line 44 of file HFFilter.h.

Constructor & Destructor Documentation

HFFilter::HFFilter ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 12 of file HFFilter.cc.

References edm::ParameterSet::getParameter().

{
  genJetsCollName_     = iConfig.getParameter<edm::InputTag>("genJetsCollName");
  ptMin_               = iConfig.getParameter<double>("ptMin");
  etaMax_              = iConfig.getParameter<double>("etaMax");
}

HFFilter::~HFFilter

(

)

Definition at line 20 of file HFFilter.cc.

{
}

Member Function Documentation

void HFFilter::endJob ( void  )

virtual

Reimplemented from edm::EDFilter.

Definition at line 72 of file HFFilter.cc.

                 {
}

bool HFFilter::filter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDFilter.

Definition at line 32 of file HFFilter.cc.

References JetMCTagUtils::decayFromBHadron(), JetMCTagUtils::decayFromCHadron(), end, edm::Event::getByLabel(), and dt_dqm_sourceclient_common_cff::reco.

{

  // Get the GenJetCollection
   using namespace edm;
   using namespace reco;
   Handle<std::vector<GenJet> > hGenJets;
   iEvent.getByLabel(genJetsCollName_,hGenJets);

   // Loop over the GenJetCollection
   vector<GenJet>::const_iterator ijet = hGenJets->begin();
   vector<GenJet>::const_iterator end  = hGenJets->end();
   for ( ; ijet != end; ++ijet ) {

     // Check to make sure the GenJet satisfies kinematic cuts. Ignore those that don't. 
     if ( ijet->pt() < ptMin_ || fabs(ijet->eta()) > etaMax_ ) continue;

     // Get the constituent particles
     vector<const GenParticle*> particles = ijet->getGenConstituents ();
    
     // Loop over the constituent particles
     vector<const GenParticle*>::const_iterator genit = particles.begin();
     vector<const GenParticle*>::const_iterator genend = particles.end();
     for ( ; genit != genend; ++genit ) {

       // See if any of them come from B or C hadrons
       const GenParticle & genitref = **genit;
       if ( JetMCTagUtils::decayFromBHadron( genitref ) ||
        JetMCTagUtils::decayFromCHadron( genitref ) ) {
     return true;
       }
     }// end loop over constituents
   }// end loop over jets


   return false;
}

Member Data Documentation

double HFFilter::etaMax_

private

Definition at line 56 of file HFFilter.h.

edm::InputTag HFFilter::genJetsCollName_

private

Definition at line 54 of file HFFilter.h.

double HFFilter::ptMin_

private

Definition at line 55 of file HFFilter.h.