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::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.

11 {}

◆ ~Pythia8PartonSelector()

Pythia8PartonSelector::~Pythia8PartonSelector ( )

override

Definition at line 13 of file Pythia8PartonSelector.cc.

13 {}

Member Function Documentation

◆ run()

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(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, ecalTrigSettings_cff::particles, dqmAnalyzer_cff::partons, 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;
 }