MCTruthHelper Namespace Reference

Functions
int	absPdgId (const reco::GenParticle &p)
int	absPdgId (const HepMC::GenParticle &p)
const reco::GenParticle *	daughter (const reco::GenParticle &p, unsigned int idau)
const HepMC::GenParticle *	daughter (const HepMC::GenParticle &p, unsigned int idau)
template<typename P >
void	fillGenStatusFlags (const P &p, reco::GenStatusFlags &statusFlags)
template<typename P >
const P *	findDecayedMother (const P &p)
template<typename P >
const P *	findDecayedMother (const P &p, int abspdgid)
template<typename P >
const P *	firstCopy (const P &p)
template<typename P >
bool	fromHardProcess (const P &p)
template<typename P >
bool	fromHardProcessBeforeFSR (const P &p)
template<typename P >
bool	fromHardProcessDecayed (const P &p)
template<typename P >
bool	fromHardProcessFinalState (const P &p)
template<typename P >
const P *	hardProcessMotherCopy (const P &p)
template<typename P >
bool	isDecayedLeptonHadron (const P &p)
template<typename P >
bool	isDirectHadronDecayProduct (const P &p)
template<typename P >
bool	isDirectHardProcessTauDecayProduct (const P &p)
template<typename P >
bool	isDirectPromptTauDecayProduct (const P &p)
template<typename P >
bool	isDirectTauDecayProduct (const P &p)
template<typename P >
bool	isFirstCopy (const P &p)
template<typename P >
bool	isHadron (const P &p)
template<typename P >
bool	isHardProcess (const P &p)
template<typename P >
bool	isHardProcessTauDecayProduct (const P &p)
template<typename P >
bool	isLastCopy (const P &p)
template<typename P >
bool	isLastCopyBeforeFSR (const P &p)
template<typename P >
bool	isMuonDecayProduct (const P &p)
template<typename P >
bool	isPrompt (const P &p)
template<typename P >
bool	isPromptDecayed (const P &p)
template<typename P >
bool	isPromptFinalState (const P &p)
template<typename P >
bool	isPromptMuonDecayProduct (const P &p)
template<typename P >
bool	isPromptTauDecayProduct (const P &p)
template<typename P >
bool	isTauDecayProduct (const P &p)
template<typename P >
const P *	lastCopy (const P &p)
template<typename P >
const P *	lastCopyBeforeFSR (const P &p)
template<typename P >
const P *	lastDaughterCopyBeforeFSR (const P &p)
const reco::GenParticle *	mother (const reco::GenParticle &p, unsigned int imoth=0)
const HepMC::GenParticle *	mother (const HepMC::GenParticle &p, unsigned int imoth=0)
template<typename P >
const P *	nextCopy (const P &p)
unsigned int	numberOfDaughters (const reco::GenParticle &p)
unsigned int	numberOfDaughters (const HepMC::GenParticle &p)
unsigned int	numberOfMothers (const reco::GenParticle &p)
unsigned int	numberOfMothers (const HepMC::GenParticle &p)
int	pdgId (const reco::GenParticle &p)
int	pdgId (const HepMC::GenParticle &p)
template<typename P >
const P *	previousCopy (const P &p)
template<typename P >
const P *	uniqueMother (const P &p)

Function Documentation

int MCTruthHelper::absPdgId

(

const reco::GenParticle &

p

)

Definition at line 574 of file MCTruthHelper.h.

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

Referenced by MCEmbeddingValidationAnalyzer::analyze(), findDecayedMother(), PseudoTopProducer::isBHadron(), isDecayedLeptonHadron(), PseudoTopProducer::produce(), TauGenJetProducer::produce(), and ParticleReplacerZtautau::testEvent().

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

int MCTruthHelper::absPdgId

(

const HepMC::GenParticle &

p

)

Definition at line 579 of file MCTruthHelper.h.

References funct::abs().

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

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

Definition at line 614 of file MCTruthHelper.h.

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

Referenced by TopDecaySubset::addDaughters(), FBaseSimEvent::addParticles(), TopDecaySubset::addRadiation(), AcceptanceHistoProducer::analyze(), EmbeddingKineReweightNtupleProducer::analyze(), JetTester_HeavyIons::analyze(), JetAnalyzer_HeavyIons::analyze(), MCEmbeddingValidationAnalyzer::analyze(), gen::Pythia6Gun::attachPy6DecaysToGenEvent(), EPOS::IO_EPOS::build_end_vertex(), QGTagger::calcVariables(), DatabasePDG::CheckImpossibleDecays(), compGenParticleP4afterRad(), tauImpactParameter::TrackHelixVertexFitter::computeMotherLorentzVectorPar(), TtDecayChannelSelector::countProngs(), DatabasePDG::DumpData(), reco::modules::JetFlavourIdentifier::fillLeptons(), TopDecaySubset::fillReferences(), JetMCTagUtils::genTauDecayMode(), getSelMuons(), PrimaryVertexAnalyzer::getSimPVs(), PrimaryVertexAnalyzer4PU::getSimPVs(), TtSemiLeptonicEvent::hadronicDecayTop(), reco::tau::xclean::CrossCleanPiZeros< PtrIter >::initialize(), reco::tau::xclean::CrossCleanPtrs< PtrIter >::initialize(), CandMatcherBase< C1, C2 >::initMaps(), lastCopyBeforeFSR(), lastDaughterCopyBeforeFSR(), TtSemiLeptonicEvent::leptonicDecayTop(), nextCopy(), reco::tau::RecoTauBuilderCombinatoricPlugin::operator()(), GenParticlesFromZsSelectorForMCEmbedding::produce(), reco::CentralityProducer::produce(), PFTau3ProngReco::produce(), gen::Cascade2Hadronizer::residualDecay(), gen::Pythia6Hadronizer::residualDecay(), CandCommonVertexFitterBase::set(), PFCandCommonVertexFitterBase::set(), CandKinematicVertexFitter::set(), TtDecayChannelSelector::tauDecay(), ParticleReplacerZtautau::testEvent(), TtFullLeptonicEvent::top(), TtFullHadronicEvent::top(), TtFullLeptonicEvent::topBar(), TtFullHadronicEvent::topBar(), gen::EvtGenInterface::update_particles(), gen::TauolappInterface::update_particles(), and TrajectoryManager::updateWithDaughters().

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

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

Definition at line 619 of file MCTruthHelper.h.

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

template<typename P >

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

Definition at line 625 of file MCTruthHelper.h.

References fromHardProcess(), fromHardProcessBeforeFSR(), isDecayedLeptonHadron(), isDirectHadronDecayProduct(), isDirectHardProcessTauDecayProduct(), isDirectPromptTauDecayProduct(), isDirectTauDecayProduct(), isFirstCopy(), isHardProcess(), isHardProcessTauDecayProduct(), isLastCopy(), isLastCopyBeforeFSR(), isPrompt(), isPromptTauDecayProduct(), isTauDecayProduct(), reco::GenStatusFlags::setFromHardProcess(), reco::GenStatusFlags::setFromHardProcessBeforeFSR(), reco::GenStatusFlags::setIsDecayedLeptonHadron(), reco::GenStatusFlags::setIsDirectHadronDecayProduct(), reco::GenStatusFlags::setIsDirectHardProcessTauDecayProduct(), reco::GenStatusFlags::setIsDirectPromptTauDecayProduct(), reco::GenStatusFlags::setIsDirectTauDecayProduct(), reco::GenStatusFlags::setIsFirstCopy(), reco::GenStatusFlags::setIsHardProcess(), reco::GenStatusFlags::setIsHardProcessTauDecayProduct(), reco::GenStatusFlags::setIsLastCopy(), reco::GenStatusFlags::setIsLastCopyBeforeFSR(), reco::GenStatusFlags::setIsPrompt(), reco::GenStatusFlags::setIsPromptTauDecayProduct(), and reco::GenStatusFlags::setIsTauDecayProduct().

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

)

Definition at line 545 of file MCTruthHelper.h.

References isDecayedLeptonHadron(), mother(), and P.

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

                                         {
    const P *mo = mother(p);
    while (mo && !isDecayedLeptonHadron(*mo)) {
      mo = mother(*mo);
    }
    return mo;
  }  

template<typename P >

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

Definition at line 555 of file MCTruthHelper.h.

References absPdgId(), isDecayedLeptonHadron(), mother(), and P.

                                                       {
    const P *mo = mother(p);
    while (mo && (absPdgId(*mo)!=abspdgid || !isDecayedLeptonHadron(*mo)) ) {
      mo = mother(*mo);
    }
    return mo;
  }

template<typename P >

const P * MCTruthHelper::firstCopy

(

const P &

p

)

Definition at line 420 of file MCTruthHelper.h.

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

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

                                 {
    const P *pcopy = &p;
    while (previousCopy(*pcopy)) {
      pcopy = previousCopy(*pcopy);
    }
    return pcopy;    
  }

template<typename P >

bool MCTruthHelper::fromHardProcess

(

const P &

p

)

Definition at line 332 of file MCTruthHelper.h.

References hardProcessMotherCopy().

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

                                   {
    return hardProcessMotherCopy(p) != 0;
  }

template<typename P >

bool MCTruthHelper::fromHardProcessBeforeFSR

(

const P &

p

)

Definition at line 364 of file MCTruthHelper.h.

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

Referenced by 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

)

Definition at line 344 of file MCTruthHelper.h.

References fromHardProcess(), and isDecayedLeptonHadron().

Referenced by isHardProcessTauDecayProduct().

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

template<typename P >

bool MCTruthHelper::fromHardProcessFinalState

(

const P &

p

)

Definition at line 338 of file MCTruthHelper.h.

References fromHardProcess().

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

template<typename P >

const P * MCTruthHelper::hardProcessMotherCopy

(

const P &

p

)

Definition at line 500 of file MCTruthHelper.h.

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

Referenced by fromHardProcess(), and 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;
    while (previousCopy(*pcopy)) {
      pcopy = previousCopy(*pcopy);
      if (isHardProcess(*pcopy)) return pcopy;
    }
    return 0;
  }

template<typename P >

bool MCTruthHelper::isDecayedLeptonHadron

(

const P &

p

)

Definition at line 235 of file MCTruthHelper.h.

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

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

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

template<typename P >

bool MCTruthHelper::isDirectHadronDecayProduct

(

const P &

p

)

Definition at line 288 of file MCTruthHelper.h.

References isDecayedLeptonHadron(), isHadron(), P, and uniqueMother().

Referenced by fillGenStatusFlags().

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

template<typename P >

bool MCTruthHelper::isDirectHardProcessTauDecayProduct

(

const P &

p

)

Definition at line 357 of file MCTruthHelper.h.

References findDecayedMother(), fromHardProcess(), P, and metsig::tau.

Referenced by 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

)

Definition at line 267 of file MCTruthHelper.h.

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

Referenced by 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

)

Definition at line 259 of file MCTruthHelper.h.

References findDecayedMother(), P, and metsig::tau.

Referenced by 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

)

Definition at line 392 of file MCTruthHelper.h.

References firstCopy().

Referenced by fillGenStatusFlags().

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

template<typename P >

bool MCTruthHelper::isHadron

(

const P &

p

)

Definition at line 295 of file MCTruthHelper.h.

References RecoTau_DiTaus_pt_20-420_cfg::ParticleID, and pdgId().

Referenced by SiStripDigitizerAlgorithm::accumulateSimHits(), isDecayedLeptonHadron(), isDirectHadronDecayProduct(), HiGammaJetSignalDef::IsIsolated(), DigiSimLinkAlgorithm::run(), and InputGenJetsParticleSelector::testPartonChildren().

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

template<typename P >

bool MCTruthHelper::isHardProcess

(

const P &

p

)

Definition at line 302 of file MCTruthHelper.h.

References firstCopy(), mother(), nextCopy(), and P.

Referenced by fillGenStatusFlags(), and 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) {
        bool fromResonance = firstCopy(*um)->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

)

Definition at line 350 of file MCTruthHelper.h.

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

Referenced by fillGenStatusFlags().

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

template<typename P >

bool MCTruthHelper::isLastCopy

(

const P &

p

)

Definition at line 398 of file MCTruthHelper.h.

References lastCopy().

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

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

template<typename P >

bool MCTruthHelper::isLastCopyBeforeFSR

(

const P &

p

)

Definition at line 404 of file MCTruthHelper.h.

References lastCopyBeforeFSR().

Referenced by fillGenStatusFlags().

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

template<typename P >

bool MCTruthHelper::isMuonDecayProduct

(

const P &

p

)

Definition at line 275 of file MCTruthHelper.h.

References findDecayedMother().

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

template<typename P >

bool MCTruthHelper::isPrompt

(

const P &

p

)

Definition at line 221 of file MCTruthHelper.h.

References findDecayedMother().

Referenced by fillGenStatusFlags(), PythiaFilterGammaGamma::filter(), isDirectPromptTauDecayProduct(), isPromptDecayed(), isPromptFinalState(), isPromptMuonDecayProduct(), and 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) == 0;
  }

template<typename P >

bool MCTruthHelper::isPromptDecayed

(

const P &

p

)

Definition at line 240 of file MCTruthHelper.h.

References isDecayedLeptonHadron(), and isPrompt().

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

template<typename P >

bool MCTruthHelper::isPromptFinalState

(

const P &

p

)

Definition at line 230 of file MCTruthHelper.h.

References isPrompt().

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

template<typename P >

bool MCTruthHelper::isPromptMuonDecayProduct

(

const P &

p

)

Definition at line 281 of file MCTruthHelper.h.

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

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

template<typename P >

bool MCTruthHelper::isPromptTauDecayProduct

(

const P &

p

)

Definition at line 252 of file MCTruthHelper.h.

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

Referenced by fillGenStatusFlags().

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

template<typename P >

bool MCTruthHelper::isTauDecayProduct

(

const P &

p

)

Definition at line 246 of file MCTruthHelper.h.

References findDecayedMother().

Referenced by fillGenStatusFlags().

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

template<typename P >

const P * MCTruthHelper::lastCopy

(

const P &

p

)

Definition at line 430 of file MCTruthHelper.h.

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

Referenced by isLastCopy().

                                {
    const P *pcopy = &p;
    while (nextCopy(*pcopy)) {
      pcopy = nextCopy(*pcopy);
    }
    return pcopy;    
  }

template<typename P >

const P * MCTruthHelper::lastCopyBeforeFSR

(

const P &

p

)

Definition at line 440 of file MCTruthHelper.h.

References daughter(), firstCopy(), numberOfDaughters(), P, and pdgId().

Referenced by isLastCopyBeforeFSR().

                                         {
    //start with first copy and then walk down until there is FSR
    const P *pcopy = firstCopy(p);
    bool hasDaughterCopy = true;
    while (hasDaughterCopy) {
      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;
        }
      }
    }
    return pcopy;       
  }

template<typename P >

const P * MCTruthHelper::lastDaughterCopyBeforeFSR

(

const P &

p

)

Definition at line 470 of file MCTruthHelper.h.

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

Referenced by fromHardProcessBeforeFSR().

                                                 {
    //start with this particle and then walk down until there is FSR
    const P *pcopy = &p;
    bool hasDaughterCopy = true;
    while (hasDaughterCopy) {
      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;
        }
      }
    }
    return pcopy;       
  }  

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

Definition at line 594 of file MCTruthHelper.h.

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

Referenced by FBaseSimEvent::addParticles(), AlignmentParameterStore::alignableFromAlignableDet(), ZLONLOHistogrammer::analyze(), EWKSystUnc::analyze(), zPdfUnc::analyze(), HLTTauRelvalQTester::analyze(), ZMuPtScaleAnalyzer::analyze(), ElectronMcSignalValidator::analyze(), TkAlCaRecoMonitor::analyze(), ElectronSeedAnalyzer::analyze(), SimplePhotonAnalyzer::analyze(), GsfElectronMCAnalyzer::analyze(), TauValidation::analyze(), PhotonValidator::analyze(), GenHFHadronMatcher::analyzeMothers(), HepMC::Pythia8ToHepMCA::append_event(), EPOS::IO_EPOS::build_production_vertex(), PedeSteerer::buildNoHierarchyCollection(), TopGenEvent::candidate(), AlignmentTwoBodyDecayTrackSelector::checkMass(), AlignmentTwoBodyDecayTrackSelector::checkMETMass(), CheckOverlap::checkPV(), KillSecondariesStackingAction::ClassifyNewTrack(), StackingAction::ClassifyNewTrack(), tauImpactParameter::TrackHelixVertexFitter::computeMotherLorentzVectorPar(), DDHCalBarrelAlgo::constructInsideLayers(), DDHCalBarrelAlgo::constructMidLayer(), DDHCalBarrelAlgo::constructSideLayer(), TopGenEvent::daughterQuarkOfWPlus(), gen::TauolappInterface::decay(), DecayChannel::DecayChannel(), PrintGeomInfoAction::dumpPV(), PrintGeomMatInfo::dumpPV(), PrintSensitive::dumpTouch(), PrintGeomInfoAction::dumpTouch(), PrintGeomMatInfo::dumpTouch(), DDLinear::execute(), DDAngular::execute(), DDPixFwdDiskAlgo::execute(), DDTrackerLinearXY::execute(), DDTrackerLinear::execute(), DDTrackerZPosAlgo::execute(), DDTrackerXYZPosAlgo::execute(), DDPixBarLayerUpgradeAlgoCoverage::execute(), DDPixBarStackLinear::execute(), DDPixBarLayerUpgradeAlgo::execute(), DDPixBarStackLayerAlgo::execute(), DDPixBarStackLinearGap::execute(), DDPixBarStackTrigLayerAlgo::execute(), DDPixBarTPGStackLayerAlgo::execute(), DDPixBarLayerAlgo::execute(), DDHCalTBZposAlgo::execute(), DDTECCoolAlgo::execute(), DDTECOptoHybAlgo::execute(), DDTrackerAngularV1::execute(), DDTrackerPhiAlgo::execute(), DDTrackerPhiAltAlgo::execute(), DDHCalLinearXY::execute(), DDHCalFibreBundle::execute(), DDTECAxialCableAlgo::execute(), DDTECPhiAlgo::execute(), DDTECPhiAltAlgo::execute(), DDTIDAxialCableAlgo::execute(), DDTIDRingAlgo::execute(), DDTrackerAngular::execute(), DDPixFwdBlades::execute(), GenParticleProducer::fillDaughters(), MuScleFitPlotter::fillGen(), PythiaFilterGammaGamma::filter(), PythiaFilter::filter(), findDecayedMother(), MuScleFitUtils::findGenMuFromRes(), MuScleFitMuonSelector::findGenMuFromRes(), HiggsValidation::findHiggsDecayChannel(), MuScleFitUtils::findSimMuFromRes(), MuScleFitMuonSelector::findSimMuFromRes(), TauValidation::findTauList(), gen::TauolappInterface::FirstTauInChain(), gen::HijingHadronizer::get_particles(), gen::AMPTHadronizer::get_particles(), gen::HydjetHadronizer::get_particles(), gen::Hydjet2Hadronizer::get_particles(), MEChannel::getAncestor(), MEChannel::getListOfAncestors(), gen::TauolappInterface::GetMother(), TauValidation::GetMother(), HeavyFlavorValidation::getMotherId(), TauValidation::GetMothers(), PrimaryVertexAnalyzer::getSimPVs(), PrimaryVertexAnalyzer4PU::getSimPVs(), PseudoTopProducer::isFromHadron(), isHardProcess(), TtGenEvent::lepton(), TtGenEvent::leptonBar(), FSimEvent::load(), gen::TauolappInterface::MatchedLHESpinUp(), TtGenEvent::neutrino(), TtGenEvent::neutrinoBar(), TopGenEvent::numberOfBQuarks(), TopGenEvent::numberOfLeptons(), ParametersToParametersDerivatives::ParametersToParametersDerivatives(), previousCopy(), FBaseSimEvent::printMCTruth(), ISRGammaWeightProducer::produce(), GenParticles2HepMCConverter::produce(), pat::PATGenCandsFromSimTracksProducer::produce(), MuonTrajectoryUpdator::propagateState(), pf2pat::TopProjectorAlgo< Top, Bottom >::ptrToAncestor(), BtoCharmDecayVertexMerger::resolveBtoDchain(), tnp::TagProbePairMaker::run(), PedeSteerer::selectCoordinateAlis(), StGenEvent::singleLepton(), TtGenEvent::singleLepton(), StGenEvent::singleNeutrino(), TtGenEvent::singleNeutrino(), and uniqueMother().

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

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

Definition at line 599 of file MCTruthHelper.h.

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

template<typename P >

const P * MCTruthHelper::nextCopy

(

const P &

p

)

Definition at line 530 of file MCTruthHelper.h.

References daughter(), numberOfDaughters(), P, and pdgId().

Referenced by isHardProcess(), and 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 0;     
  }  

unsigned int MCTruthHelper::numberOfDaughters

(

const reco::GenParticle &

p

)

Definition at line 604 of file MCTruthHelper.h.

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

Referenced by lastCopyBeforeFSR(), lastDaughterCopyBeforeFSR(), gen.GeneratorAnalyzer.GeneratorAnalyzer::makeMCInfo(), TemplatedSecondaryVertexProducer< IPTI, VTX >::matchSubjets(), JetFlavourClustering::matchSubjets(), nextCopy(), reco::PFClusterJet::print(), FlavorHistoryProducer::produce(), and reco::TrackJet::tracks().

                                                           {
    return p.numberOfDaughters();
  }

unsigned int MCTruthHelper::numberOfDaughters

(

const HepMC::GenParticle &

p

)

Definition at line 609 of file MCTruthHelper.h.

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

unsigned int MCTruthHelper::numberOfMothers

(

const reco::GenParticle &

p

)

Definition at line 584 of file MCTruthHelper.h.

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

Referenced by GenParticleProducer::fillDaughters(), EcalGenEvtSelectorFrag::filter(), and previousCopy().

                                                         {
    return p.numberOfMothers();
  }

unsigned int MCTruthHelper::numberOfMothers

(

const HepMC::GenParticle &

p

)

Definition at line 589 of file MCTruthHelper.h.

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

int MCTruthHelper::pdgId

(

const reco::GenParticle &

p

)

Definition at line 564 of file MCTruthHelper.h.

References reco::LeafCandidate::pdgId().

Referenced by isHadron(), lastCopyBeforeFSR(), lastDaughterCopyBeforeFSR(), nextCopy(), previousCopy(), and uniqueMother().

                                      {
    return p.pdgId();
  }

int MCTruthHelper::pdgId

(

const HepMC::GenParticle &

p

)

Definition at line 569 of file MCTruthHelper.h.

                                       {
    return p.pdg_id();
  }  

template<typename P >

const P * MCTruthHelper::previousCopy

(

const P &

p

)

Definition at line 515 of file MCTruthHelper.h.

References mother(), numberOfMothers(), P, and pdgId().

Referenced by firstCopy(), and 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 0;     
  }   

template<typename P >

const P * MCTruthHelper::uniqueMother

(

const P &

p

)

Definition at line 410 of file MCTruthHelper.h.

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

Referenced by isDirectHadronDecayProduct().

                                    {
    const P *mo = &p;
    while (mo && pdgId(*mo)==pdgId(p)) {
      mo = mother(*mo);
    }
    return mo;
  }