builds a GenJet from the visible daughters of each status 2 tau in the event. More...

#include <TauGenJetProducer.h>

Inheritance diagram for TauGenJetProducer:

Public Member Functions
void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

	TauGenJetProducer (const edm::ParameterSet &)

	~TauGenJetProducer () override

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const bool	includeNeutrinos_
	if yes, neutrinos will be included, for debug purposes More...

const edm::InputTag	inputTagGenParticles_
	Input PFCandidates. More...

const edm::EDGetTokenT< reco::GenParticleCollection >	tokenGenParticles_

const bool	verbose_
	verbose ? More...

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::EDProducerBase
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

builds a GenJet from the visible daughters of each status 2 tau in the event.

Author: Colin Bernet

Date: february 2008

Definition at line 25 of file TauGenJetProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file TauGenJetProducer.cc.

                                                                    :
   inputTagGenParticles_(iConfig.getParameter<InputTag>("GenParticles")),
   tokenGenParticles_(consumes<GenParticleCollection>(inputTagGenParticles_)),
   includeNeutrinos_(iConfig.getParameter<bool>("includeNeutrinos")),
   verbose_(iConfig.getUntrackedParameter<bool>("verbose",false))
 {
   
 
   produces<GenJetCollection>();
 }

TauGenJetProducer::~TauGenJetProducer ( )

override

Definition at line 30 of file TauGenJetProducer.cc.

30 { }

Member Function Documentation

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

override

Definition at line 32 of file TauGenJetProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, gather_cfg::cout, DEFINE_FWK_MODULE, edm::RefVector< C, T, F >::empty(), edm::RefVector< C, T, F >::end(), Exception, GenParticlesHelper::findDescendents(), GenParticlesHelper::findParticles(), runEdmFileComparison::found, GenHFHadronMatcher_cfi::genParticles, edm::Event::getByToken(), includeNeutrinos_, inputTagGenParticles_, metsig::jet, eostools::move(), edm::Event::put(), edm::refToPtr(), reco::LeafCandidate::setCharge(), tokenGenParticles_, and verbose_.

                                                 {
 
   Handle<GenParticleCollection> genParticles;
 
   bool found = iEvent.getByToken( tokenGenParticles_, genParticles);
 
   if ( !found ) {
     std::ostringstream  err;
     err<<" cannot get collection: "
        <<inputTagGenParticles_<<std::endl;
     edm::LogError("TauGenJetProducer")<<err.str();
     throw cms::Exception( "MissingProduct", err.str());
   }
 
   auto pOutVisTaus = std::make_unique<GenJetCollection>();
 
   using namespace GenParticlesHelper;
 
   GenParticleRefVector allStatus2Taus;
   findParticles( *genParticles,
          allStatus2Taus, 15, 2);
 
   for ( IGR iTau=allStatus2Taus.begin(); iTau!=allStatus2Taus.end(); ++iTau ) {
 
     // look for all status 1 (stable) descendents
     GenParticleRefVector descendents;
     findDescendents( *iTau, descendents, 1);
 
     // CV: skip status 2 taus that radiate-off a photon
     //    --> have a status 2 tau lepton in the list of descendents
     GenParticleRefVector status2TauDaughters;
     findDescendents( *iTau, status2TauDaughters, 2, 15 );
     if ( !status2TauDaughters.empty() ) continue;
 
     // loop on descendents, and take all except neutrinos
     math::XYZTLorentzVector sumVisMom;
     Particle::Charge charge = 0;
     Jet::Constituents constituents;
 
     if(verbose_)
       cout<<"tau "<<(*iTau)<<endl;
 
     for(IGR igr = descendents.begin();
     igr!= descendents.end(); ++igr ) {
 
       int absPdgId = abs((*igr)->pdgId());
 
       // neutrinos
       if(!includeNeutrinos_ ) {
     if( absPdgId == 12 ||
         absPdgId == 14 ||
         absPdgId == 16  )
       continue;
       }
 
       if(verbose_)
     cout<<"\t"<<(*igr)<<endl;
 
       charge += (*igr)->charge();
       sumVisMom += (*igr)->p4();
 
       // need to convert the vector of reference to the constituents
       // to a vector of pointers to build the genjet
       constituents.push_back( refToPtr( *igr) );
     }
 
     math::XYZPoint vertex;
     GenJet::Specific specific;
 
     GenJet jet( sumVisMom, vertex, specific, constituents);
 
     if (charge != (*iTau)->charge() )
       std::cout<<" charge of Tau: " << (*iTau) << " not equal to charge of sum of charge of all descendents. " << std::endl;
 
     jet.setCharge(charge);
     pOutVisTaus->push_back( jet );
 
   }
   iEvent.put(std::move(pOutVisTaus) );
 }

Member Data Documentation

const bool TauGenJetProducer::includeNeutrinos_

private

if yes, neutrinos will be included, for debug purposes

Definition at line 41 of file TauGenJetProducer.h.

Referenced by produce().

const edm::InputTag TauGenJetProducer::inputTagGenParticles_

private

Input PFCandidates.

Definition at line 37 of file TauGenJetProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::GenParticleCollection> TauGenJetProducer::tokenGenParticles_

private

Definition at line 38 of file TauGenJetProducer.h.

Referenced by produce().

const bool TauGenJetProducer::verbose_

private

verbose ?

Definition at line 44 of file TauGenJetProducer.h.

Referenced by produce().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation