#include <MCTruthHelper.h>

Public Member Functions
int	absPdgId (const reco::GenParticle &p) const

int	absPdgId (const HepMC::GenParticle &p) const

const reco::GenParticle *	daughter (const reco::GenParticle &p, unsigned int idau) const

const HepMC::GenParticle *	daughter (const HepMC::GenParticle &p, unsigned int idau) const

void	fillGenStatusFlags (const P &p, reco::GenStatusFlags &statusFlags) const

const P *	findDecayedMother (const P &p) const

const P *	findDecayedMother (const P &p, int abspdgid) const

const P *	firstCopy (const P &p) const

bool	fromHardProcess (const P &p) const

bool	fromHardProcessBeforeFSR (const P &p) const

bool	fromHardProcessDecayed (const P &p) const

bool	fromHardProcessFinalState (const P &p) const

const P *	hardProcessMotherCopy (const P &p) const

bool	isDecayedLeptonHadron (const P &p) const

bool	isDirectHadronDecayProduct (const P &p) const

bool	isDirectHardProcessTauDecayProduct (const P &p) const

bool	isDirectPromptTauDecayProduct (const P &p) const

bool	isDirectTauDecayProduct (const P &p) const

bool	isFirstCopy (const P &p) const

bool	isHadron (const P &p) const

bool	isHardProcess (const P &p) const

bool	isHardProcessTauDecayProduct (const P &p) const

bool	isLastCopy (const P &p) const

bool	isLastCopyBeforeFSR (const P &p) const

bool	isMuonDecayProduct (const P &p) const

bool	isPrompt (const P &p) const

bool	isPromptDecayed (const P &p) const

bool	isPromptFinalState (const P &p) const

bool	isPromptMuonDecayProduct (const P &p) const

bool	isPromptTauDecayProduct (const P &p) const

bool	isTauDecayProduct (const P &p) const

const P *	lastCopy (const P &p) const

const P *	lastCopyBeforeFSR (const P &p) const

const P *	lastDaughterCopyBeforeFSR (const P &p) const

const reco::GenParticle *	mother (const reco::GenParticle &p, unsigned int imoth=0) const

const HepMC::GenParticle *	mother (const HepMC::GenParticle &p, unsigned int imoth=0) const

const P *	nextCopy (const P &p) const

unsigned int	numberOfDaughters (const reco::GenParticle &p) const

unsigned int	numberOfDaughters (const HepMC::GenParticle &p) const

unsigned int	numberOfMothers (const reco::GenParticle &p) const

unsigned int	numberOfMothers (const HepMC::GenParticle &p) const

int	pdgId (const reco::GenParticle &p) const

int	pdgId (const HepMC::GenParticle &p) const

const P *	previousCopy (const P &p) const

const P *	uniqueMother (const P &p) const

Detailed Description

template<typename P>
class MCTruthHelper

Definition at line 14 of file MCTruthHelper.h.

Member Function Documentation

template<typename P >

int MCTruthHelper::absPdgId ( const reco::GenParticle & p ) const

Definition at line 573 of file MCTruthHelper.h.

References funct::abs(), and reco::LeafCandidate::pdgId().

Referenced by MCTruthHelper::findDecayedMother(), and MCTruthHelper::isDecayedLeptonHadron().

                                                              {
   return std::abs(p.pdgId());
 }

template<typename P >

int MCTruthHelper::absPdgId ( const HepMC::GenParticle & p ) const

Definition at line 579 of file MCTruthHelper.h.

References funct::abs().

                                                               {
   return std::abs(p.pdg_id());
 }

template<typename P >

const reco::GenParticle * MCTruthHelper< P >::daughter	(	const reco::GenParticle &	p,
		unsigned int	idau
	)		const

Definition at line 621 of file MCTruthHelper.h.

References reco::CompositeRefCandidateT< D >::daughter().

Referenced by MCTruthHelper::lastCopyBeforeFSR(), MCTruthHelper::lastDaughterCopyBeforeFSR(), and MCTruthHelper::nextCopy().

                                                                                                    {
   return static_cast<const reco::GenParticle*>(p.daughter(idau));
 }

template<typename P >

const HepMC::GenParticle * MCTruthHelper< P >::daughter	(	const HepMC::GenParticle &	p,
		unsigned int	idau
	)		const

Definition at line 627 of file MCTruthHelper.h.

                                                                                                      {
   return *(p.end_vertex()->particles_out_const_begin() + idau);
 }

template<typename P>

void MCTruthHelper< P >::fillGenStatusFlags	(	const P &	p,
		reco::GenStatusFlags &	statusFlags
	)		const

Definition at line 633 of file MCTruthHelper.h.

Referenced by GenParticleProducer::convertParticle().

                                                                                            {
   statusFlags.setIsPrompt(isPrompt(p));
   statusFlags.setIsDecayedLeptonHadron(isDecayedLeptonHadron(p));
   statusFlags.setIsTauDecayProduct(isTauDecayProduct(p));
   statusFlags.setIsPromptTauDecayProduct(isPromptTauDecayProduct(p));
   statusFlags.setIsDirectTauDecayProduct(isDirectTauDecayProduct(p));
   statusFlags.setIsDirectPromptTauDecayProduct(isDirectPromptTauDecayProduct(p));
   statusFlags.setIsDirectHadronDecayProduct(isDirectHadronDecayProduct(p));
   statusFlags.setIsHardProcess(isHardProcess(p));
   statusFlags.setFromHardProcess(fromHardProcess(p));
   statusFlags.setIsHardProcessTauDecayProduct(isHardProcessTauDecayProduct(p));
   statusFlags.setIsDirectHardProcessTauDecayProduct(isDirectHardProcessTauDecayProduct(p));
   statusFlags.setFromHardProcessBeforeFSR(fromHardProcessBeforeFSR(p));
   statusFlags.setIsFirstCopy(isFirstCopy(p));
   statusFlags.setIsLastCopy(isLastCopy(p));
   statusFlags.setIsLastCopyBeforeFSR(isLastCopyBeforeFSR(p));
 }  

template<typename P>

const P * MCTruthHelper::findDecayedMother ( const P & p ) const

Definition at line 535 of file MCTruthHelper.h.

References MCTruthHelper::isDecayedLeptonHadron(), and MCTruthHelper::mother().

Referenced by MCTruthHelper::isDirectHardProcessTauDecayProduct(), MCTruthHelper::isDirectPromptTauDecayProduct(), MCTruthHelper::isDirectTauDecayProduct(), MCTruthHelper::isHardProcessTauDecayProduct(), MCTruthHelper::isMuonDecayProduct(), MCTruthHelper::isPrompt(), MCTruthHelper::isPromptMuonDecayProduct(), MCTruthHelper::isPromptTauDecayProduct(), and MCTruthHelper::isTauDecayProduct().

                                                              {
   const P *mo = mother(p);
   std::unordered_set<const P*> dupCheck;
   while (mo && !isDecayedLeptonHadron(*mo)) {
     dupCheck.insert(mo);
     mo = mother(*mo);
     if (dupCheck.count(mo)) return nullptr;
   }
   return mo;
 }  

template<typename P>

const P * MCTruthHelper< P >::findDecayedMother	(	const P &	p,
		int	abspdgid
	)		const

Definition at line 548 of file MCTruthHelper.h.

References MCTruthHelper::absPdgId(), MCTruthHelper::isDecayedLeptonHadron(), and MCTruthHelper::mother().

                                                                            {
   const P *mo = mother(p);
   std::unordered_set<const P*> dupCheck;
   while (mo && (absPdgId(*mo)!=abspdgid || !isDecayedLeptonHadron(*mo)) ) {
     dupCheck.insert(mo);
     mo = mother(*mo);
     if (dupCheck.count(mo)) return nullptr;
   }
   return mo;
 }

template<typename P>

const P * MCTruthHelper::firstCopy ( const P & p ) const

Definition at line 394 of file MCTruthHelper.h.

References AlCaHLTBitMon_ParallelJobs::p, and MCTruthHelper::previousCopy().

Referenced by MCTruthHelper::isFirstCopy(), MCTruthHelper::isHardProcess(), and MCTruthHelper::lastCopyBeforeFSR().

                                                      {
   const P *pcopy = &p;
   std::unordered_set<const P*> dupCheck;
   while (previousCopy(*pcopy)) {
     dupCheck.insert(pcopy);
     pcopy = previousCopy(*pcopy);
     if (dupCheck.count(pcopy)) return nullptr;
   }
   return pcopy;    
 }

template<typename P>

bool MCTruthHelper::fromHardProcess ( const P & p ) const

Definition at line 303 of file MCTruthHelper.h.

References MCTruthHelper::hardProcessMotherCopy().

Referenced by MCTruthHelper::fillGenStatusFlags(), MCTruthHelper::fromHardProcessDecayed(), MCTruthHelper::fromHardProcessFinalState(), and MCTruthHelper::isDirectHardProcessTauDecayProduct().

                                                        {
   return hardProcessMotherCopy(p) != nullptr;
 }

template<typename P>

bool MCTruthHelper::fromHardProcessBeforeFSR ( const P & p ) const

Definition at line 335 of file MCTruthHelper.h.

References MCTruthHelper::hardProcessMotherCopy(), MCTruthHelper::isLastCopy(), MCTruthHelper::lastDaughterCopyBeforeFSR(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                 {
   //pythia 6 documentation line roughly corresponds to this condition
   if (p.status()==3) return true;
   
   //check hard process mother properties
   const P *hpc = hardProcessMotherCopy(p);
   if (!hpc) return false;
       
   //for incoming partons in pythia8, more useful information is not
   //easily available, so take only the incoming parton itself
   if (hpc->status()==21 && (&p)==hpc) return true;
   
   //for intermediate particles in pythia 8, just take the last copy
   if (hpc->status()==22 && isLastCopy(p)) return true;
   
   //for outgoing particles in pythia 8, explicitly find the last copy
   //before FSR starting from the hardProcess particle, and take only
   //this one
   if ( (hpc->status()==23 || hpc->status()==1) && (&p)==lastDaughterCopyBeforeFSR(*hpc) ) return true;
   
   
   //didn't satisfy any of the conditions
   return false;
   
 }

template<typename P>

bool MCTruthHelper::fromHardProcessDecayed ( const P & p ) const

Definition at line 315 of file MCTruthHelper.h.

References MCTruthHelper::fromHardProcess(), and MCTruthHelper::isDecayedLeptonHadron().

Referenced by MCTruthHelper::isHardProcessTauDecayProduct().

                                                               {
   return isDecayedLeptonHadron(p) && fromHardProcess(p);
 }

template<typename P>

bool MCTruthHelper::fromHardProcessFinalState ( const P & p ) const

Definition at line 309 of file MCTruthHelper.h.

References MCTruthHelper::fromHardProcess().

                                                                  {
   return p.status()==1 && fromHardProcess(p);
 }

template<typename P>

const P * MCTruthHelper::hardProcessMotherCopy ( const P & p ) const

Definition at line 487 of file MCTruthHelper.h.

References MCTruthHelper::isHardProcess(), AlCaHLTBitMon_ParallelJobs::p, and MCTruthHelper::previousCopy().

Referenced by MCTruthHelper::fromHardProcess(), and MCTruthHelper::fromHardProcessBeforeFSR().

                                                                  {
   //is particle itself is hard process particle
   if (isHardProcess(p)) return &p;
   
   //check if any other copies are hard process particles
   const P *pcopy = &p;
   std::unordered_set<const P*> dupCheck;
   while (previousCopy(*pcopy)) {
     dupCheck.insert(pcopy);
     pcopy = previousCopy(*pcopy);
     if (isHardProcess(*pcopy)) return pcopy;
     if (dupCheck.count(pcopy)) break;
   }
   return nullptr;
 }

template<typename P>

bool MCTruthHelper::isDecayedLeptonHadron ( const P & p ) const

Definition at line 205 of file MCTruthHelper.h.

References MCTruthHelper::absPdgId(), MCTruthHelper::isHadron(), and MCTruthHelper::isLastCopy().

Referenced by MCTruthHelper::fillGenStatusFlags(), MCTruthHelper::findDecayedMother(), MCTruthHelper::fromHardProcessDecayed(), MCTruthHelper::isDirectHadronDecayProduct(), and MCTruthHelper::isPromptDecayed().

                                                              {
   return p.status()==2  && (isHadron(p) || absPdgId(p)==13 || absPdgId(p)==15) && isLastCopy(p);
 }

template<typename P>

bool MCTruthHelper::isDirectHadronDecayProduct ( const P & p ) const

Definition at line 258 of file MCTruthHelper.h.

References MCTruthHelper::isDecayedLeptonHadron(), MCTruthHelper::isHadron(), and MCTruthHelper::uniqueMother().

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                   {
   const P *um = uniqueMother(p);
   return um && isHadron(*um) && isDecayedLeptonHadron(*um);
 }  

template<typename P>

bool MCTruthHelper::isDirectHardProcessTauDecayProduct ( const P & p ) const

Definition at line 328 of file MCTruthHelper.h.

References symbols::dm, MCTruthHelper::findDecayedMother(), MCTruthHelper::fromHardProcess(), and metsig::tau.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                           {
   const P *tau = findDecayedMother(p,15);
   const P *dm = findDecayedMother(p);
   return tau && tau==dm && fromHardProcess(*tau);    
 }  

template<typename P>

bool MCTruthHelper::isDirectPromptTauDecayProduct ( const P & p ) const

Definition at line 237 of file MCTruthHelper.h.

References symbols::dm, MCTruthHelper::findDecayedMother(), MCTruthHelper::isPrompt(), and metsig::tau.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                      {
   const P *tau = findDecayedMother(p,15);
   const P *dm = findDecayedMother(p);
   return tau && tau==dm && isPrompt(*tau);
 }

template<typename P>

bool MCTruthHelper::isDirectTauDecayProduct ( const P & p ) const

Definition at line 229 of file MCTruthHelper.h.

References symbols::dm, MCTruthHelper::findDecayedMother(), and metsig::tau.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                {
   const P *tau = findDecayedMother(p,15);
   const P *dm = findDecayedMother(p);
   return tau && tau==dm;
 }

template<typename P>

bool MCTruthHelper::isFirstCopy ( const P & p ) const

Definition at line 363 of file MCTruthHelper.h.

References MCTruthHelper::firstCopy().

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                    {
   return &p == firstCopy(p);
 }  

template<typename P>

bool MCTruthHelper::isHadron ( const P & p ) const

Definition at line 265 of file MCTruthHelper.h.

References source_particleGun_cfi::ParticleID, and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::isDecayedLeptonHadron(), and MCTruthHelper::isDirectHadronDecayProduct().

                                                 {
   HepPDT::ParticleID heppdtid(pdgId(p));
   return heppdtid.isHadron();
 }

template<typename P>

bool MCTruthHelper::isHardProcess ( const P & p ) const

Definition at line 272 of file MCTruthHelper.h.

References MCTruthHelper::firstCopy(), MCTruthHelper::mother(), and MCTruthHelper::nextCopy().

Referenced by MCTruthHelper::fillGenStatusFlags(), and MCTruthHelper::hardProcessMotherCopy().

                                                      {
   
   //status 3 in pythia6 means hard process;
   if (p.status()==3) return true;
   
   //hard process codes for pythia8 are 21-29 inclusive (currently 21,22,23,24 are used)
   if (p.status()>20 && p.status()<30) return true;
   
   //if this is a final state or decayed particle,
   //check if direct mother is a resonance decay in pythia8 but exclude FSR branchings
   //(In pythia8 if a resonance decay product did not undergo any further branchings
   //it will be directly stored as status 1 or 2 without any status 23 copy)
   if (p.status()==1 || p.status()==2) {
     const P *um = mother(p);
     if (um) {
       const P *firstcopy = firstCopy(*um);
       bool fromResonance = firstcopy && firstcopy->status()==22;
       
       const P *umNext = nextCopy(*um);
       bool fsrBranching = umNext && umNext->status()>50 && umNext->status()<60;
       
       if (fromResonance && !fsrBranching) return true;
     }
   }
   
   return false;
   
 }

template<typename P>

bool MCTruthHelper::isHardProcessTauDecayProduct ( const P & p ) const

Definition at line 321 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother(), MCTruthHelper::fromHardProcessDecayed(), and metsig::tau.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                     {
   const P *tau = findDecayedMother(p,15);
   return tau && fromHardProcessDecayed(*tau);    
 }

template<typename P>

bool MCTruthHelper::isLastCopy ( const P & p ) const

Definition at line 369 of file MCTruthHelper.h.

References MCTruthHelper::lastCopy().

Referenced by MCTruthHelper::fillGenStatusFlags(), MCTruthHelper::fromHardProcessBeforeFSR(), and MCTruthHelper::isDecayedLeptonHadron().

                                                   {
   return &p == lastCopy(p);
 }

template<typename P>

bool MCTruthHelper::isLastCopyBeforeFSR ( const P & p ) const

Definition at line 375 of file MCTruthHelper.h.

References MCTruthHelper::lastCopyBeforeFSR().

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                            {
   return &p == lastCopyBeforeFSR(p);
 }  

template<typename P>

bool MCTruthHelper::isMuonDecayProduct ( const P & p ) const

Definition at line 245 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother().

                                                           {
   return findDecayedMother(p,13) != 0;
 }

template<typename P>

bool MCTruthHelper::isPrompt ( const P & p ) const

Definition at line 191 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother().

Referenced by MCTruthHelper::fillGenStatusFlags(), MCTruthHelper::isDirectPromptTauDecayProduct(), MCTruthHelper::isPromptDecayed(), MCTruthHelper::isPromptFinalState(), MCTruthHelper::isPromptMuonDecayProduct(), and MCTruthHelper::isPromptTauDecayProduct().

                                                 {
   //particle from hadron/muon/tau decay -> not prompt
   //checking all the way up the chain treats properly the radiated photon
   //case as well
   return findDecayedMother(p) == nullptr;
 }

template<typename P>

bool MCTruthHelper::isPromptDecayed ( const P & p ) const

Definition at line 210 of file MCTruthHelper.h.

References MCTruthHelper::isDecayedLeptonHadron(), and MCTruthHelper::isPrompt().

                                                        {
   return isDecayedLeptonHadron(p) && isPrompt(p);
 }

template<typename P>

bool MCTruthHelper::isPromptFinalState ( const P & p ) const

Definition at line 200 of file MCTruthHelper.h.

References MCTruthHelper::isPrompt().

                                                           {
   return p.status()==1 && isPrompt(p);
 }

template<typename P>

bool MCTruthHelper::isPromptMuonDecayProduct ( const P & p ) const

Definition at line 251 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother(), MCTruthHelper::isPrompt(), and RPCpg::mu.

                                                                 {
   const P *mu = findDecayedMother(p,13);
   return mu && isPrompt(*mu);
 }  

template<typename P>

bool MCTruthHelper::isPromptTauDecayProduct ( const P & p ) const

Definition at line 222 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother(), MCTruthHelper::isPrompt(), and metsig::tau.

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                                {
   const P *tau = findDecayedMother(p,15);
   return tau && isPrompt(*tau);
 }

template<typename P>

bool MCTruthHelper::isTauDecayProduct ( const P & p ) const

Definition at line 216 of file MCTruthHelper.h.

References MCTruthHelper::findDecayedMother().

Referenced by MCTruthHelper::fillGenStatusFlags().

                                                          {
   return findDecayedMother(p,15) != nullptr;
 }

template<typename P>

const P * MCTruthHelper::lastCopy ( const P & p ) const

Definition at line 407 of file MCTruthHelper.h.

References MCTruthHelper::nextCopy(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by MCTruthHelper::isLastCopy().

                                                     {
   const P *pcopy = &p;
   std::unordered_set<const P*> dupCheck;
   while (nextCopy(*pcopy)) {
     dupCheck.insert(pcopy);
     pcopy = nextCopy(*pcopy);
     if (dupCheck.count(pcopy)) return nullptr;
   }
   return pcopy;    
 }

template<typename P>

const P * MCTruthHelper::lastCopyBeforeFSR ( const P & p ) const

Definition at line 420 of file MCTruthHelper.h.

References MCTruthHelper::daughter(), MCTruthHelper::firstCopy(), MCTruthHelper::numberOfDaughters(), and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::isLastCopyBeforeFSR().

                                                              {
   //start with first copy and then walk down until there is FSR
   const P *pcopy = firstCopy(p);
   if (!pcopy) return nullptr;
   bool hasDaughterCopy = true;
   std::unordered_set<const P*> dupCheck;
   while (hasDaughterCopy) {
     dupCheck.insert(pcopy);
     hasDaughterCopy = false;
     const unsigned int ndau = numberOfDaughters(*pcopy);
     //look for FSR
     for (unsigned int idau = 0; idau<ndau; ++idau) {
       const P *dau = daughter(*pcopy,idau);
       if (pdgId(*dau)==21 || pdgId(*dau)==22) {
         //has fsr (or else decayed and is the last copy by construction)
         return pcopy;
       }        
     }
     //look for daughter copy
     for (unsigned int idau = 0; idau<ndau; ++idau) {
       const P *dau = daughter(*pcopy,idau);
       if (pdgId(*dau)==pdgId(p)) {
         pcopy = dau;
         hasDaughterCopy = true;
         break;
       }
     }
     if (dupCheck.count(pcopy)) return nullptr;
   }
   return pcopy;       
 }

template<typename P>

const P * MCTruthHelper::lastDaughterCopyBeforeFSR ( const P & p ) const

Definition at line 454 of file MCTruthHelper.h.

References MCTruthHelper::daughter(), MCTruthHelper::numberOfDaughters(), AlCaHLTBitMon_ParallelJobs::p, and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::fromHardProcessBeforeFSR().

                                                                      {
   //start with this particle and then walk down until there is FSR
   const P *pcopy = &p;
   bool hasDaughterCopy = true;
   std::unordered_set<const P*> dupCheck;
   while (hasDaughterCopy) {
     dupCheck.insert(pcopy);
     hasDaughterCopy = false;
     const unsigned int ndau = numberOfDaughters(*pcopy);
     //look for FSR
     for (unsigned int idau = 0; idau<ndau; ++idau) {
       const P *dau = daughter(*pcopy,idau);
       if (pdgId(*dau)==21 || pdgId(*dau)==22) {
         //has fsr (or else decayed and is the last copy by construction)
         return pcopy;
       }        
     }
     //look for daughter copy
     for (unsigned int idau = 0; idau<ndau; ++idau) {
       const P *dau = daughter(*pcopy,idau);
       if (pdgId(*dau)==pdgId(p)) {
         pcopy = dau;
         hasDaughterCopy = true;
         break;
       }
     }
     if (dupCheck.count(pcopy)) return nullptr;
   }
   return pcopy;       
 }  

template<typename P >

const reco::GenParticle * MCTruthHelper< P >::mother	(	const reco::GenParticle &	p,
		unsigned int	imoth = `0`
	)		const

Definition at line 597 of file MCTruthHelper.h.

References reco::CompositeRefCandidateT< D >::mother().

Referenced by MCTruthHelper::findDecayedMother(), MCTruthHelper::isHardProcess(), MCTruthHelper::previousCopy(), and MCTruthHelper::uniqueMother().

                                                                                                   {
   return static_cast<const reco::GenParticle*>(p.mother(imoth));
 }

template<typename P >

const HepMC::GenParticle * MCTruthHelper< P >::mother	(	const HepMC::GenParticle &	p,
		unsigned int	imoth = `0`
	)		const

Definition at line 603 of file MCTruthHelper.h.

                                                                                                     {
   return p.production_vertex() && p.production_vertex()->particles_in_size() ? *(p.production_vertex()->particles_in_const_begin() + imoth) : nullptr;
 }

template<typename P>

const P * MCTruthHelper::nextCopy ( const P & p ) const

Definition at line 520 of file MCTruthHelper.h.

References MCTruthHelper::daughter(), MCTruthHelper::numberOfDaughters(), and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::isHardProcess(), and MCTruthHelper::lastCopy().

                                                     {
   
   const unsigned int ndau = numberOfDaughters(p);
   for (unsigned int idau = 0; idau<ndau; ++idau) {
     const P *dau = daughter(p,idau);
     if (pdgId(*dau)==pdgId(p)) {
       return dau;
     }
   }
   
   return nullptr;     
 }  

template<typename P >

unsigned int MCTruthHelper::numberOfDaughters ( const reco::GenParticle & p ) const

Definition at line 609 of file MCTruthHelper.h.

References reco::CompositeRefCandidateT< D >::numberOfDaughters().

Referenced by MCTruthHelper::lastCopyBeforeFSR(), MCTruthHelper::lastDaughterCopyBeforeFSR(), and MCTruthHelper::nextCopy().

                                                                                {
   return p.numberOfDaughters();
 }

template<typename P >

unsigned int MCTruthHelper::numberOfDaughters ( const HepMC::GenParticle & p ) const

Definition at line 615 of file MCTruthHelper.h.

                                                                                 {
   return p.end_vertex() ? p.end_vertex()->particles_out_size() : 0;
 }

template<typename P >

unsigned int MCTruthHelper::numberOfMothers ( const reco::GenParticle & p ) const

Definition at line 585 of file MCTruthHelper.h.

References reco::CompositeRefCandidateT< D >::numberOfMothers().

Referenced by MCTruthHelper::previousCopy().

                                                                              {
   return p.numberOfMothers();
 }

template<typename P >

unsigned int MCTruthHelper::numberOfMothers ( const HepMC::GenParticle & p ) const

Definition at line 591 of file MCTruthHelper.h.

                                                                               {
   return p.production_vertex() ? p.production_vertex()->particles_in_size() : 0;
 }

template<typename P >

int MCTruthHelper::pdgId ( const reco::GenParticle & p ) const

Definition at line 561 of file MCTruthHelper.h.

References reco::LeafCandidate::pdgId().

Referenced by Particle.Particle::__str__(), MCTruthHelper::isHadron(), MCTruthHelper::lastCopyBeforeFSR(), MCTruthHelper::lastDaughterCopyBeforeFSR(), MCTruthHelper::nextCopy(), MCTruthHelper::previousCopy(), and MCTruthHelper::uniqueMother().

                                                           {
   return p.pdgId();
 }

template<typename P >

int MCTruthHelper::pdgId ( const HepMC::GenParticle & p ) const

Definition at line 567 of file MCTruthHelper.h.

Referenced by Particle.Particle::__str__().

                                                            {
   return p.pdg_id();
 }  

template<typename P>

const P * MCTruthHelper::previousCopy ( const P & p ) const

Definition at line 505 of file MCTruthHelper.h.

References MCTruthHelper::mother(), MCTruthHelper::numberOfMothers(), and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::firstCopy(), and MCTruthHelper::hardProcessMotherCopy().

                                                         {
   
   const unsigned int nmoth = numberOfMothers(p);
   for (unsigned int imoth = 0; imoth<nmoth; ++imoth) {
     const P *moth = mother(p,imoth);
     if (pdgId(*moth)==pdgId(p)) {
       return moth;
     }
   }
   
   return nullptr;     
 }   

template<typename P>

const P * MCTruthHelper::uniqueMother ( const P & p ) const

Definition at line 381 of file MCTruthHelper.h.

References MCTruthHelper::mother(), AlCaHLTBitMon_ParallelJobs::p, and MCTruthHelper::pdgId().

Referenced by MCTruthHelper::isDirectHadronDecayProduct().

                                                         {
   const P *mo = &p;
   std::unordered_set<const P*> dupCheck;
   while (mo && pdgId(*mo)==pdgId(p)) {
     dupCheck.insert(mo);
     mo = mother(*mo);
     if (dupCheck.count(mo)) return nullptr;
   }
   return mo;
 }

Public Member Functions

Detailed Description

template<typename P> class MCTruthHelper< P >

Member Function Documentation

template<typename P>
class MCTruthHelper< P >