---

TauGenJetProducer Class Reference

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

#include <PhysicsTools/JetMCAlgos/plugins/TauGenJetProducer.h>

Inheritance diagram for TauGenJetProducer:

List of all members.

Public Member Functions
virtual void	beginJob (const edm::EventSetup &c)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	TauGenJetProducer (const edm::ParameterSet &)
	~TauGenJetProducer ()
Private Attributes
bool	includeNeutrinos_
	if yes, neutrinos will be included, for debug purposes
edm::InputTag	inputTagGenParticles_
	Input PFCandidates.
bool	verbose_
	verbose ?

---

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 23 of file TauGenJetProducer.h.

---

Constructor & Destructor Documentation

TauGenJetProducer::TauGenJetProducer

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 22 of file TauGenJetProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), includeNeutrinos_, inputTagGenParticles_, and verbose_.

                                                                   {
  

  inputTagGenParticles_ 
    = iConfig.getParameter<InputTag>("GenParticles");

  includeNeutrinos_ 
    = iConfig.getParameter<bool>("includeNeutrinos");

  verbose_ = 
    iConfig.getUntrackedParameter<bool>("verbose",false);



  produces<GenJetCollection>();
}

TauGenJetProducer::~TauGenJetProducer

(

)

Definition at line 41 of file TauGenJetProducer.cc.

{ }

---

Member Function Documentation

void TauGenJetProducer::beginJob

(

const edm::EventSetup &

c

)

[virtual]

Reimplemented from edm::EDProducer.

Definition at line 45 of file TauGenJetProducer.cc.

{ }

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

Implements edm::EDProducer.

Definition at line 48 of file TauGenJetProducer.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), GenMuonPlsPt100GeV_cfg::cout, edm::RefVector< C, T, F >::end(), lat::endl(), err, Exception, GenParticlesHelper::findDescendents(), GenParticlesHelper::findParticles(), genParticles_cfi::genParticles, edm::Event::getByLabel(), includeNeutrinos_, inputTagGenParticles_, metsig::jet, edm::Event::put(), edm::refToPtr(), and verbose_.

                                                          {
  
  Handle<GenParticleCollection> genParticles;

  bool found = iEvent.getByLabel( inputTagGenParticles_, 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);
    
    // 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); 
    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 40 of file TauGenJetProducer.h.

Referenced by produce(), and TauGenJetProducer().

edm::InputTag TauGenJetProducer::inputTagGenParticles_ [private]

Input PFCandidates.

Definition at line 37 of file TauGenJetProducer.h.

Referenced by produce(), and TauGenJetProducer().

bool TauGenJetProducer::verbose_ [private]

verbose ?

Definition at line 43 of file TauGenJetProducer.h.

Referenced by produce(), and TauGenJetProducer().

---

The documentation for this class was generated from the following files:

• PhysicsTools/JetMCAlgos/plugins/TauGenJetProducer.h
• PhysicsTools/JetMCAlgos/plugins/TauGenJetProducer.cc

---