TauGenJetProducer Class Reference

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

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.

{ }

Member Function Documentation

void TauGenJetProducer::beginJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 37 of file TauGenJetProducer.cc.

{ }

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

virtual

Implements edm::EDProducer.

Definition at line 39 of file TauGenJetProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), DeDxDiscriminatorTools::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 43 of file TauGenJetProducer.h.

edm::InputTag TauGenJetProducer::inputTagGenParticles_

private

Input PFCandidates.

Definition at line 39 of file TauGenJetProducer.h.

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

private

Definition at line 40 of file TauGenJetProducer.h.

bool TauGenJetProducer::verbose_

private

verbose ?

Definition at line 46 of file TauGenJetProducer.h.