GenParticlesFromZsSelectorForMCEmbedding Class Reference

#include <GenParticlesFromZsSelectorForMCEmbedding.h>

Inheritance diagram for GenParticlesFromZsSelectorForMCEmbedding:

Public Member Functions
	GenParticlesFromZsSelectorForMCEmbedding (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &)
	~GenParticlesFromZsSelectorForMCEmbedding ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	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
std::vector< ConsumesInfo >	consumesInfo () const
	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
enum	{ kBeforeFSR, kAfterFSR }
typedef std::vector< int >	vint

Private Attributes
int	before_or_afterFSR_
int	maxDaughters_
int	minDaughters_
vint	pdgIdsDaughters_
vint	pdgIdsMothers_
edm::InputTag	src_
int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Select tau leptons, muons and electrons produced in Z/gamma* –> l+ l- processes

NOTE: Can handle case that virtual Z/gamma is missing in GenEVT record

CV: This class is a copy of the copy of the code in TauAnalysis/GenSimTools/plugins/GenParticlesFromZsSelector.h TauAnalysis/GenSimTools/plugins/GenParticlesFromZsSelector.cc The copy was created to avoid package dependencies (of either TauAnalysis/MCEmbeddingTools on TauAnalysis/GenSimTools or of TauAnalysis/GenSimTools on TauAnalysis/MCEmbeddingTools... I know this is not nice)

Author

Christian Veelken, LLR

Version

Revision:

1.2

Id:

GenParticlesFromZsSelectorForMCEmbedding.h,v 1.2 2013/01/31 09:07:18 veelken Exp

Definition at line 33 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Member Typedef Documentation

typedef std::vector<int> GenParticlesFromZsSelectorForMCEmbedding::vint

private

Definition at line 47 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
kBeforeFSR
kAfterFSR

Definition at line 54 of file GenParticlesFromZsSelectorForMCEmbedding.h.

{ kBeforeFSR, kAfterFSR };

Constructor & Destructor Documentation

GenParticlesFromZsSelectorForMCEmbedding::GenParticlesFromZsSelectorForMCEmbedding ( const edm::ParameterSet &  cfg )

explicit

Definition at line 13 of file GenParticlesFromZsSelectorForMCEmbedding.cc.

References before_or_afterFSR_, Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), kAfterFSR, kBeforeFSR, maxDaughters_, minDaughters_, pdgIdsDaughters_, pdgIdsMothers_, src_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

{
  src_ = cfg.getParameter<edm::InputTag>("src");

  pdgIdsMothers_ = cfg.getParameter<vint>("pdgIdsMothers");
  pdgIdsDaughters_ = cfg.getParameter<vint>("pdgIdsDaughters");

  maxDaughters_ = cfg.getParameter<int>("maxDaughters");
  minDaughters_ = cfg.getParameter<int>("minDaughters");

  std::string before_or_afterFSR_string = cfg.getParameter<std::string>("before_or_afterFSR");
  if      ( before_or_afterFSR_string == "beforeFSR" ) before_or_afterFSR_ = kBeforeFSR;
  else if ( before_or_afterFSR_string == "afterFSR"  ) before_or_afterFSR_ = kAfterFSR;
  else throw cms::Exception("Configuration")
    << " Invalid Configuration Parameter 'before_or_afterFSR' = " << before_or_afterFSR_string << " !!\n";

  verbosity_ = ( cfg.exists("verbosity") ) ? 
    cfg.getParameter<int>("verbosity") : 0;

  produces<reco::GenParticleCollection>("");
}

GenParticlesFromZsSelectorForMCEmbedding::~GenParticlesFromZsSelectorForMCEmbedding

(

)

Definition at line 35 of file GenParticlesFromZsSelectorForMCEmbedding.cc.

{ 
// nothing to be done yet...
}

Member Function Documentation

void GenParticlesFromZsSelectorForMCEmbedding::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDProducer.

Definition at line 98 of file GenParticlesFromZsSelectorForMCEmbedding.cc.

References assert(), before_or_afterFSR_, reco::Candidate::begin(), gather_cfg::cout, reco::Candidate::end(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByLabel(), customizeTrackingMonitorSeedNumber::idx, kAfterFSR, kBeforeFSR, maxDaughters_, minDaughters_, pdgIdsDaughters_, pdgIdsMothers_, reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), edm::Event::put(), src_, and verbosity_.

{
  if ( verbosity_ ) {
    std::cout << "<GenParticlesFromZsSelectorForMCEmbedding::produce>:" << std::endl;
    std::cout << " src = " << src_ << std::endl;
  }

  edm::Handle<reco::GenParticleCollection> genParticles;
  evt.getByLabel(src_, genParticles);
  
  std::vector<const reco::GenParticle*> genParticlesFromZs_tmp;

//--- check if HepEVT record contains any Z/gamma* --> l+ l- entries
//
//    NOTE: iteration over mothers in the outer loop
//          gives Z --> l+ l- decays priority over gamma --> e+ e- 
//
  for ( vint::const_iterator pdgIdMother = pdgIdsMothers_.begin();
    pdgIdMother != pdgIdsMothers_.end(); ++pdgIdMother ) {
    for ( vint::const_iterator pdgIdDaughter = pdgIdsDaughters_.begin();
      pdgIdDaughter != pdgIdsDaughters_.end(); ++pdgIdDaughter ) {
      if ( (int)genParticlesFromZs_tmp.size() < maxDaughters_ || maxDaughters_ == -1 )
    findGenParticles(*genParticles, *pdgIdMother, *pdgIdDaughter, genParticlesFromZs_tmp);
    }
  }
  
//--- check if HepEVT record contains l+ l- entries without Z/gamma* parent
//    
//    NOTE: in order to avoid ambiguities, give preference to l+ l- pairs
//          which are subsequent in HepEVT record
//
  if ( !((int)genParticlesFromZs_tmp.size() >= minDaughters_) ) {
    for ( vint::const_iterator pdgIdDaughter = pdgIdsDaughters_.begin();
      pdgIdDaughter != pdgIdsDaughters_.end(); ++pdgIdDaughter ) {
      findGenParticles(*genParticles, *pdgIdDaughter, minDaughters_, true, genParticlesFromZs_tmp);
    }
  }
  if ( !((int)genParticlesFromZs_tmp.size() >= minDaughters_) ) {
    for ( vint::const_iterator pdgIdDaughter = pdgIdsDaughters_.begin();
      pdgIdDaughter != pdgIdsDaughters_.end(); ++pdgIdDaughter ) {
      findGenParticles(*genParticles, *pdgIdDaughter, minDaughters_, false, genParticlesFromZs_tmp);
    }
  }
  
  std::auto_ptr<reco::GenParticleCollection> genParticlesFromZs(new reco::GenParticleCollection());
  
  int idx = 0;
  for ( std::vector<const reco::GenParticle*>::const_iterator genParticleFromZ_beforeFSR = genParticlesFromZs_tmp.begin();
    genParticleFromZ_beforeFSR != genParticlesFromZs_tmp.end(); ++genParticleFromZ_beforeFSR ) {
    if ( before_or_afterFSR_ == kBeforeFSR ) {
      genParticlesFromZs->push_back(**genParticleFromZ_beforeFSR);
    } else if ( before_or_afterFSR_ == kAfterFSR ) {
      std::vector<const reco::GenParticle*> daughters;
      findDaughters(*genParticleFromZ_beforeFSR, daughters, -1);
      const reco::GenParticle* genParticleFromZ_afterFSR = (*genParticleFromZ_beforeFSR);
      for ( std::vector<const reco::GenParticle*>::const_iterator daughter = daughters.begin();
        daughter != daughters.end(); ++daughter ) {
    if ( (*daughter)->pdgId() == (*genParticleFromZ_beforeFSR)->pdgId() && 
         (*daughter)->energy() < genParticleFromZ_afterFSR->energy() ) genParticleFromZ_afterFSR = (*daughter);
      }
      if ( verbosity_ ) {
    std::cout << "genParticleFromZ #" << idx << " (beforeFSR): Pt = " << (*genParticleFromZ_beforeFSR)->pt() << "," 
          << " eta = " << (*genParticleFromZ_beforeFSR)->eta() << ", phi = " << (*genParticleFromZ_beforeFSR)->phi() << std::endl;
    std::cout << "genParticleFromZ #" << idx << " (afterFSR): Pt = " << genParticleFromZ_afterFSR->pt() << "," 
          << " eta = " << genParticleFromZ_afterFSR->eta() << ", phi = " << genParticleFromZ_afterFSR->phi() << std::endl;
      }
      genParticlesFromZs->push_back(*genParticleFromZ_afterFSR);
    } else assert(0);
    ++idx;
  }
  
  evt.put(genParticlesFromZs);
}

Member Data Documentation

int GenParticlesFromZsSelectorForMCEmbedding::before_or_afterFSR_

private

Definition at line 55 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

int GenParticlesFromZsSelectorForMCEmbedding::maxDaughters_

private

Definition at line 51 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

int GenParticlesFromZsSelectorForMCEmbedding::minDaughters_

private

Definition at line 52 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

vint GenParticlesFromZsSelectorForMCEmbedding::pdgIdsDaughters_

private

Definition at line 49 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

vint GenParticlesFromZsSelectorForMCEmbedding::pdgIdsMothers_

private

Definition at line 48 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

edm::InputTag GenParticlesFromZsSelectorForMCEmbedding::src_

private

Definition at line 45 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().

int GenParticlesFromZsSelectorForMCEmbedding::verbosity_

private

Definition at line 57 of file GenParticlesFromZsSelectorForMCEmbedding.h.

Referenced by GenParticlesFromZsSelectorForMCEmbedding(), and produce().