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
virtual void	produce (edm::Event &, const edm::EventSetup &) override

	TauGenJetProducer (const edm::ParameterSet &)

	~TauGenJetProducer ()

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 (const std::string &iProcessName, std::vector< const char * > &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 Attributes
bool	includeNeutrinos_
	if yes, neutrinos will be included, for debug purposes More...

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

edm::EDGetTokenT < reco::GenParticleCollection >	tokenGenParticles_

bool	verbose_
	verbose ? More...

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

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.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

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

TauGenJetProducer::~TauGenJetProducer ( )

Definition at line 35 of file TauGenJetProducer.cc.

35 { }

Member Function Documentation

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

overridevirtual

Implements edm::EDProducer.

Definition at line 37 of file TauGenJetProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), RecoTauCleanerPlugins::charge, gather_cfg::cout, edm::RefVector< C, T, F >::end(), edm::hlt::Exception, GenParticlesHelper::findDescendents(), GenParticlesHelper::findParticles(), newFWLiteAna::found, genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByToken(), metsig::jet, edm::Event::put(), edm::refToPtr(), reco::LeafCandidate::setCharge(), edm::RefVector< C, T, F >::size(), and timingPdfMaker::specific.

                                           {
 
   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());
   }
 
   std::auto_ptr<GenJetCollection>
     pOutVisTaus(new 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.size() > 0 ) 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::cout;
 
     jet.setCharge(charge);
     pOutVisTaus->push_back( jet );
 
   }
   iEvent.put( pOutVisTaus );
 }

Member Data Documentation

bool TauGenJetProducer::includeNeutrinos_

private

if yes, neutrinos will be included, for debug purposes

Definition at line 41 of file TauGenJetProducer.h.

edm::InputTag TauGenJetProducer::inputTagGenParticles_

private

Input PFCandidates.

Definition at line 37 of file TauGenJetProducer.h.

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

private

Definition at line 38 of file TauGenJetProducer.h.

bool TauGenJetProducer::verbose_

private

verbose ?

Definition at line 44 of file TauGenJetProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation