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::unique_ptr< reco::GenParticleRefVector > &partons) override
	~Pythia8PartonSelector () override
Public Member Functions inherited from BasePartonSelector
	BasePartonSelector ()
virtual	~BasePartonSelector ()

Detailed Description

Pythia8 parton selector derived from the base parton selector.

Definition at line 10 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 11 of file Pythia8PartonSelector.cc.

{}

Pythia8PartonSelector::~Pythia8PartonSelector ( )

override

Definition at line 13 of file Pythia8PartonSelector.cc.

{}

Member Function Documentation

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

overridevirtual

Reimplemented from BasePartonSelector.

Definition at line 15 of file Pythia8PartonSelector.cc.

References mps_fire::i, CandMCTagUtils::isParton(), and mps_update::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;
}