Selects hadrons and partons from a collection of GenParticles. More...

#include <PhysicsTools/JetMCAlgos/plugins/HadronAndPartonSelector.cc>

Inheritance diagram for HadronAndPartonSelector:

Public Member Functions
	HadronAndPartonSelector (const edm::ParameterSet &)

	~HadronAndPartonSelector ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	beginJob ()

virtual void	beginLuminosityBlock (edm::LuminosityBlock &, edm::EventSetup const &)

virtual void	beginRun (edm::Run &, edm::EventSetup const &)

virtual void	endJob ()

virtual void	endLuminosityBlock (edm::LuminosityBlock &, edm::EventSetup const &)

virtual void	endRun (edm::Run &, edm::EventSetup const &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
const edm::EDGetTokenT < reco::GenParticleCollection >	particlesToken_

std::string	partonMode_

PartonSelectorPtr	partonSelector_

const edm::EDGetTokenT < GenEventInfoProduct >	srcToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

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

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Selects hadrons and partons from a collection of GenParticles.

This producer selects hadrons, partons, and leptons from a collection of GenParticles and stores vectors of EDM references to these particles in the event. The following hadrons are selected:

b hadrons that do not have other b hadrons as daughters
c hadrons that do not have other c hadrons as daughters

Older Fortran Monte Carlo generators (Pythia6 and Herwig6) follow the HEPEVT [1] particle status code convention while newer C++ Monte Carlo generators (Pythia8, Herwig++, and Sherpa) follow the HepMC [2] particle status code convention. However, both conventions give considerable freedom in defining the status codes of intermediate particle states. Hence, the parton selection is generator-dependent and is described in each of the parton selectors separately.

Using the provenance information of the GenEventInfoProduct, the producer attempts to automatically determine what generator was used to hadronize events and based on that information decides what parton selection mode to use. It is also possible to enforce any of the supported parton selection modes.

The selected hadrons and partons are finally used by the JetFlavourClustering producer to determine the jet flavour.

The following leptons are selected:

status==1 electrons and muons
status==2 taus

[1] http://cepa.fnal.gov/psm/stdhep/ [2] http://lcgapp.cern.ch/project/simu/HepMC/

Definition at line 79 of file HadronAndPartonSelector.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 111 of file HadronAndPartonSelector.cc.

                                                                                :
 
   srcToken_(mayConsume<GenEventInfoProduct>( iConfig.getParameter<edm::InputTag>("src") )),
   particlesToken_(consumes<reco::GenParticleCollection>( iConfig.getParameter<edm::InputTag>("particles") )),
   partonMode_(iConfig.getParameter<std::string>("partonMode"))
 
 {
    //register your products
    produces<reco::GenParticleRefVector>( "bHadrons" );
    produces<reco::GenParticleRefVector>( "cHadrons" );
    produces<reco::GenParticleRefVector>( "partons" );
    produces<reco::GenParticleRefVector>( "leptons" );
 }

HadronAndPartonSelector::~HadronAndPartonSelector ( )

Definition at line 126 of file HadronAndPartonSelector.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void HadronAndPartonSelector::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 263 of file HadronAndPartonSelector.cc.

264 {

265 }

void HadronAndPartonSelector::beginLuminosityBlock	(	edm::LuminosityBlock &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 286 of file HadronAndPartonSelector.cc.

287 {

288 }

void HadronAndPartonSelector::beginRun	(	edm::Run &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 274 of file HadronAndPartonSelector.cc.

275 {

276 }

void HadronAndPartonSelector::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 269 of file HadronAndPartonSelector.cc.

269 {

270 }

void HadronAndPartonSelector::endLuminosityBlock	(	edm::LuminosityBlock &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 292 of file HadronAndPartonSelector.cc.

293 {

294 }

void HadronAndPartonSelector::endRun	(	edm::Run &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 280 of file HadronAndPartonSelector.cc.

281 {

282 }

void HadronAndPartonSelector::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 298 of file HadronAndPartonSelector.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

void HadronAndPartonSelector::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 141 of file HadronAndPartonSelector.cc.

References edm::hlt::Exception, edm::Event::getByToken(), edm::Event::getProvenance(), CandMCTagUtils::hasBottom(), CandMCTagUtils::hasCharm(), i, edm::HandleBase::id(), reco::isElectron(), reco::isMuon(), reco::isTau(), edm::HandleBase::isValid(), EgammaValidation_Wenu_cff::leptons, edm::moduleName(), particlesToken_, partonMode_, partonSelector_, edm::Event::put(), srcToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
    // determine hadronizer type (done only once per job)
    if( partonMode_=="Auto" )
    {
      edm::Handle<GenEventInfoProduct> genEvtInfoProduct;
      iEvent.getByToken(srcToken_, genEvtInfoProduct);
 
      std::string moduleName = "";
      if( genEvtInfoProduct.isValid() ) {
        const edm::Provenance& prov = iEvent.getProvenance(genEvtInfoProduct.id());
        moduleName = edm::moduleName(prov);
      }
 
      if( moduleName.find("Pythia6")!=std::string::npos )
        partonMode_="Pythia6";
      else if( moduleName.find("Pythia8")!=std::string::npos )
        partonMode_="Pythia8";
      else if( moduleName.find("Herwig6")!=std::string::npos )
        partonMode_="Herwig6";
      else if( moduleName.find("ThePEG")!=std::string::npos )
        partonMode_="Herwig++";
      else if( moduleName.find("Sherpa")!=std::string::npos )
        partonMode_="Sherpa";
      else
        partonMode_="Undefined";
    }
 
    // set the parton selection mode (done only once per job)
    if( !partonSelector_ )
    {
      if ( partonMode_=="Undefined" )
        edm::LogWarning("UndefinedPartonMode") << "Could not automatically determine the hadronizer type and set the correct parton selection mode. Parton-based jet flavour will not be defined.";
      else if ( partonMode_=="Pythia6" )
      {
        partonSelector_ = PartonSelectorPtr( new Pythia6PartonSelector() );
        edm::LogInfo("PartonModeDefined") << "Using Pythia6 parton selection mode.";
      }
      else if ( partonMode_=="Pythia8" )
      {
        partonSelector_ = PartonSelectorPtr( new Pythia8PartonSelector() );
        edm::LogInfo("PartonModeDefined") << "Using Pythia8 parton selection mode.";
      }
      else if ( partonMode_=="Herwig6" )
      {
        partonSelector_ = PartonSelectorPtr( new Herwig6PartonSelector() );
        edm::LogInfo("PartonModeDefined") << "Using Herwig6 parton selection mode.";
      }
      else if ( partonMode_=="Herwig++" )
      {
        partonSelector_ = PartonSelectorPtr( new HerwigppPartonSelector() );
        edm::LogInfo("PartonModeDefined") << "Using Herwig++ parton selection mode.";
      }
      else if ( partonMode_=="Sherpa" )
      {
        partonSelector_ = PartonSelectorPtr( new SherpaPartonSelector() );
        edm::LogInfo("PartonModeDefined") << "Using Sherpa parton selection mode.";
      }
      else
        throw cms::Exception("InvalidPartonMode") <<"Parton selection mode is invalid: " << partonMode_ << ", use Auto | Pythia6 | Pythia8 | Herwig6 | Herwig++ | Sherpa" << std::endl;
    }
 
    edm::Handle<reco::GenParticleCollection> particles;
    iEvent.getByToken(particlesToken_, particles);
 
    std::auto_ptr<reco::GenParticleRefVector> bHadrons ( new reco::GenParticleRefVector );
    std::auto_ptr<reco::GenParticleRefVector> cHadrons ( new reco::GenParticleRefVector );
    std::auto_ptr<reco::GenParticleRefVector> partons  ( new reco::GenParticleRefVector );
    std::auto_ptr<reco::GenParticleRefVector> leptons  ( new reco::GenParticleRefVector );
 
    // loop over particles and select b and c hadrons and leptons
    for(reco::GenParticleCollection::const_iterator it = particles->begin(); it != particles->end(); ++it)
    {
      // if b hadron
      if( CandMCTagUtils::hasBottom( *it ) )
      {
        // check if any of the daughters is also a b hadron
        bool hasbHadronDaughter = false;
        for(size_t i=0; i < it->numberOfDaughters(); ++i)
        {
          if( CandMCTagUtils::hasBottom( *(it->daughter(i)) ) ) { hasbHadronDaughter = true; break; }
        }
        if( hasbHadronDaughter ) continue; // skip excited b hadrons that have other b hadrons as daughters
 
        bHadrons->push_back( reco::GenParticleRef( particles, it - particles->begin() ) );
      }
 
      // if c hadron
      if( CandMCTagUtils::hasCharm( *it ) )
      {
        // check if any of the daughters is also a c hadron
        bool hascHadronDaughter = false;
        for(size_t i=0; i < it->numberOfDaughters(); ++i)
        {
          if( CandMCTagUtils::hasCharm( *(it->daughter(i)) ) ) { hascHadronDaughter = true; break; }
        }
        if( hascHadronDaughter ) continue; // skip excited c hadrons that have other c hadrons as daughters
 
        cHadrons->push_back( reco::GenParticleRef( particles, it - particles->begin() ) );
      }
 
      // status==1 electrons and muons
      if( ( reco::isElectron( *it ) || reco::isMuon( *it ) ) && it->status()==1 )
        leptons->push_back( reco::GenParticleRef( particles, it - particles->begin() ) );
 
      // status==2 taus
      if( reco::isTau( *it ) && it->status()==2 )
        leptons->push_back( reco::GenParticleRef( particles, it - particles->begin() ) );
    }
 
    // select partons
    if ( partonMode_!="Undefined" )
      partonSelector_->run(particles,partons);
 
    iEvent.put( bHadrons, "bHadrons" );
    iEvent.put( cHadrons, "cHadrons" );
    iEvent.put( partons,  "partons" );
    iEvent.put( leptons,  "leptons" );
 }

Member Data Documentation

const edm::EDGetTokenT<reco::GenParticleCollection> HadronAndPartonSelector::particlesToken_

private

Definition at line 98 of file HadronAndPartonSelector.cc.

Referenced by produce().

std::string HadronAndPartonSelector::partonMode_

private

Definition at line 100 of file HadronAndPartonSelector.cc.

Referenced by produce().

PartonSelectorPtr HadronAndPartonSelector::partonSelector_

private

Definition at line 101 of file HadronAndPartonSelector.cc.

Referenced by produce().

const edm::EDGetTokenT<GenEventInfoProduct> HadronAndPartonSelector::srcToken_

private

Definition at line 97 of file HadronAndPartonSelector.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation