Pythia8PartonSelector Class Reference

Pythia8 parton selector derived from the base parton selector. More...

#include <PhysicsTools/JetMCAlgos/interface/Pythia8PartonSelector.h>

Inheritance diagram for Pythia8PartonSelector:

Public Member Functions
	Pythia8PartonSelector ()
void	run (const edm::Handle< reco::GenParticleCollection > &particles, std::auto_ptr< reco::GenParticleRefVector > &partons)
virtual	~Pythia8PartonSelector ()
Public Member Functions inherited from BasePartonSelector
	BasePartonSelector ()
	~BasePartonSelector ()

Detailed Description

Pythia8 parton selector derived from the base parton selector.

Definition at line 11 of file Pythia8PartonSelector.h.

Constructor & Destructor Documentation

Pythia8PartonSelector::Pythia8PartonSelector

(

)

This is a Pythia8-specific parton selector that selects all partons that don't have other partons as daughters, i.e., partons from the end of the parton showering sequence. An explanation of the particle status codes returned by Pythia8 can be found in Pythia8 online manual (http://home.thep.lu.se/~torbjorn/pythia81html/ParticleProperties.html).

Definition at line 12 of file Pythia8PartonSelector.cc.

{
}

Pythia8PartonSelector::~Pythia8PartonSelector ( )

virtual

Definition at line 16 of file Pythia8PartonSelector.cc.

{
}

Member Function Documentation

void Pythia8PartonSelector::run ( const edm::Handle< reco::GenParticleCollection > &  particles, std::auto_ptr< reco::GenParticleRefVector > &  partons  )

virtual

Reimplemented from BasePartonSelector.

Definition at line 21 of file Pythia8PartonSelector.cc.

References i, CandMCTagUtils::isParton(), and ntuplemaker::status.

{
   // loop over particles and select partons
   for(reco::GenParticleCollection::const_iterator it = particles->begin(); it != particles->end(); ++it)
   {
     int status = it->status();
     if( status==1 ) continue;                         // skip stable particles
     if( status==2 ) continue;                         // skip decayed Standard Model hadrons and leptons
     if( !CandMCTagUtils::isParton( *it ) ) continue;  // skip particle if not a parton

     // check if the parton has other partons as daughters
     int nparton_daughters = 0;
     for(unsigned i=0; i<it->numberOfDaughters(); ++i)
     {
       if( CandMCTagUtils::isParton( *(it->daughter(i)) ) )
         ++nparton_daughters;
     }

     if( nparton_daughters==0 )
       partons->push_back( reco::GenParticleRef( particles, it - particles->begin() ) );
   }

   return;
}