Herwig++ parton selector derived from the base parton selector. More...

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

Inheritance diagram for HerwigppPartonSelector:

Public Member Functions
	HerwigppPartonSelector ()

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

	~HerwigppPartonSelector () override

Public Member Functions inherited from BasePartonSelector
	BasePartonSelector ()

virtual	~BasePartonSelector ()

Detailed Description

Herwig++ parton selector derived from the base parton selector.

Definition at line 10 of file HerwigppPartonSelector.h.

Constructor & Destructor Documentation

◆ HerwigppPartonSelector()

HerwigppPartonSelector::HerwigppPartonSelector ( )

This is a Herwig++-specific parton selector that selects all status==2 partons. This is likely a temporary choice since Herwig++ status codes in CMSSW currently break the HepMC convention. Once the status codes are fixed, the selector will be updated.

Definition at line 11 of file HerwigppPartonSelector.cc.

11 {}

◆ ~HerwigppPartonSelector()

HerwigppPartonSelector::~HerwigppPartonSelector ( )

override

Definition at line 13 of file HerwigppPartonSelector.cc.

13 {}

Member Function Documentation

◆ run()

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

overridevirtual

Reimplemented from BasePartonSelector.

Definition at line 15 of file HerwigppPartonSelector.cc.

References CandMCTagUtils::isParton(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, primaryVertexAssociation_cfi::particles, and dqmAnalyzer_cff::partons.

                                                                                    {
   // loop over particles and select partons
   for (reco::GenParticleCollection::const_iterator it = particles->begin(); it != particles->end(); ++it) {
     if (it->status() != 2)
       continue;  // only accept status==2 particles
     if (!CandMCTagUtils::isParton(*it))
       continue;  // skip particle if not a parton
 
     partons->push_back(reco::GenParticleRef(particles, it - particles->begin()));
   }
 
   return;
 }