DerivationFramework::MenuTruthThinning Class Reference

#include <MenuTruthThinning.h>

Inheritance diagram for DerivationFramework::MenuTruthThinning:

Collaboration diagram for DerivationFramework::MenuTruthThinning:

Public Member Functions
	MenuTruthThinning (const std::string &t, const std::string &n, const IInterface *p)
virtual	~MenuTruthThinning ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual StatusCode	doThinning () const override
bool	isAccepted (const xAOD::TruthParticle *) const
bool	matchHadronIncTau (const xAOD::TruthParticle *part) const
bool	matchQuarkIncTau (const xAOD::TruthParticle *part) const
bool	isOrphanIncTau (const xAOD::TruthParticle *part) const
bool	matchGenParticle (const xAOD::TruthParticle *part, std::vector< int > &targetIDs, std::vector< int > &intermediateIDs, bool targetsAreRange) const
bool	isLeptonFromTau (const xAOD::TruthParticle *) const
bool	isFromTau (const xAOD::TruthParticle *, std::unordered_set< int > &barcode_trace) const
bool	isFsrFromLepton (const xAOD::TruthParticle *) const
bool	parentIsLongLived (const xAOD::TruthParticle *) const

Static Public Member Functions
static bool	isttHFHadron (const xAOD::TruthParticle *)

Private Attributes
StringProperty	m_streamName { this, "StreamName", "", "Name of the stream being thinned" }
SG::ThinningHandleKey< xAOD::TruthParticleContainer >	m_particlesKey { this, "ParticlesKey", "TruthParticles", "TruthParticle container name" }
SG::ThinningHandleKey< xAOD::TruthVertexContainer >	m_verticesKey { this, "VerticesKey", "TruthVertices", "TruthVertex container name" }
bool	m_writePartons
	Parameter: Keep partons? More...
bool	m_writeHadrons
	Parameter: Keep hadrons? More...
bool	m_writeBHadrons
	Parameter: Keep b-hadrons? More...
bool	m_writeCHadrons
	Parameter: Keep c-hadrons? More...
bool	m_writeGeant
	Parameter: Keep geant particles? More...
float	m_geantPhotonPtThresh
	Parameter: Write Geant photons with Pt above this threshold. More...
bool	m_writeTauHad
	Parameter: Keep hadronic tau decays? More...
bool	m_writeBSM
	Parameter: Keep BSM particles? More...
bool	m_writeBosons
	Parameter: Keep bosons? More...
float	m_photonPtCut
float	m_partonPtThresh
	Parameter: Write partons with Pt above this threshold. More...
bool	m_writeEverything
	Parameter: Write absolutely everything. More...
bool	m_writeBosonProducts
	Parameter: Write boson decay products. More...
bool	m_writeBSMProducts
	Parameter: Write BSM decay products. More...
bool	m_writeTopAndDecays
	Parameter: Write top and decay products. More...
bool	m_writeAllLeptons
	Parameter: Write all leptons. More...
bool	m_writeLeptonsNotFromHadrons
	Parameter: Write all leptons. More...
bool	m_writeAllStable
bool	m_writeNotPhysical
	Parameter: Write particles with status code 3. More...
bool	m_writettHFHadrons
	Parameter: Write particles for tt+HF classification. More...
int	m_writeFirstN
	Parameter: First N particles to write. More...
bool	m_preserveDescendants
	Parameter: preserve descendant/ancestor graph completeness. More...
bool	m_preserveGeneratorDescendants
bool	m_preserveAncestors
bool	m_preserveImmediate
bool	m_preserveHadVtx
std::vector< int >	m_pdgIdsToKeep
	Parameter: List of PDG IDs to always keep. More...
std::vector< int >	m_longLivedPdgIds
	Parameter: List of PDG IDs of long lived particles so that one can keep their children. More...
std::atomic< unsigned int >	m_totpart
std::atomic< unsigned int >	m_removedpart
std::atomic< int >	m_eventCount {}

Detailed Description

Definition at line 27 of file MenuTruthThinning.h.

Constructor & Destructor Documentation

MenuTruthThinning()

DerivationFramework::MenuTruthThinning::MenuTruthThinning	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 30 of file MenuTruthThinning.cxx.

                                                                                :
base_class(t,n,p),
m_writeFirstN(-1),
m_totpart(0),
m_removedpart(0)
{
    declareProperty ("WritePartons",
                     m_writePartons = true,
                     "Keep partons?");
    
    declareProperty ("WriteHadrons",
                     m_writeHadrons = false,
                     "Keep hadrons?");
    
    declareProperty ("WriteBHadrons",
                     m_writeBHadrons = true,
                     "Keep b-hadrons?");
 
    declareProperty ("WriteCHadrons",
                     m_writeCHadrons = true,
                     "Keep c-hadrons?");
 
    declareProperty ("WriteGeant",
                     m_writeGeant = false,
                     "Keep Geant particles?");
    
    declareProperty ("GeantPhotonPtThresh",
                     m_geantPhotonPtThresh = 0.5*GeV,
                     "Write Geant photons with Pt above this threshold.  "
                     "Set to < 0 to not write any.");
    
    declareProperty ("WriteTauHad",
                     m_writeTauHad = false,
                     "Keep hadronic taus?");
    
    declareProperty ("PartonPtThresh",
                     m_partonPtThresh = -1,
                     "Write partons with Pt above this threshold.");
    
    declareProperty ("WriteBSM",
                     m_writeBSM = false,
                     "Keep BSM particles?");
    
    declareProperty ("WriteBosons",
                     m_writeBosons = false,
                     "Keep bosons?");
    
    declareProperty ("PhotonPtThresh",
                     m_photonPtCut = 3000.,
                     "Photon pt cut with WriteBosons");
    
    declareProperty ("WriteBSMProducts",
                     m_writeBSMProducts = false,
                     "Keep BSM particle decay products?");
    
    declareProperty ("WriteBosonProducts",
                     m_writeBosonProducts = false,
                     "Keep boson decay products?");
    
    declareProperty ("WriteTopAndDecays",
                     m_writeTopAndDecays = false,
                     "Keep top partons and immediate decay products?");
    
    declareProperty("WriteEverything",
                    m_writeEverything = false,
                    "Keep absolutely everything (overrides all other flags)");
    
    declareProperty("WriteAllLeptons",
                    m_writeAllLeptons = false,
                    "Keep absolutely all leptons");
    
    declareProperty("WriteLeptonsNotFromHadrons",
                    m_writeLeptonsNotFromHadrons = false,
                    "Keep leptons not from hadron decays");
    
    declareProperty("WriteAllStable",
                    m_writeAllStable = false,
                    "Keep all stable particles");
    
    declareProperty("WriteNotPhysical",
                    m_writeNotPhysical = false,
                    "Save also non-physical particles");
    
    declareProperty("WriteFirstN",
                    m_writeFirstN = -1,
                    "Keep first N particles in record");
    
    declareProperty("PreserveDescendants",
                    m_preserveDescendants = false,
                    "Preserve descendants of retained particles");
    
    declareProperty("PreserveGeneratorDescendants",
                    m_preserveGeneratorDescendants = false,
                    "Preserve descendants of retained particles excluding Geant particles");
    
    declareProperty("PreserveAncestors",
                    m_preserveAncestors = false,
                    "Preserve ancestors of retained particles");
    
    declareProperty ("PreserveParentsSiblingsChildren",
                     m_preserveImmediate = false,
                     "Preserve the parents, siblings and children of retained particles");
    
    declareProperty ("PreserveHadronizationVertices",
                     m_preserveHadVtx = false,
                     "Preserve hadronization vertices with parents/children.");
    
    declareProperty ("WritettHFHadrons",
                     m_writettHFHadrons = false,
                     "Keep tt+HF hadrons?");
 
    declareProperty ("PDGIDsToKeep",
                     m_pdgIdsToKeep={},
                     "List of PDG IDs to always keep");
 
    declareProperty ("LongLivedPDGIDs",
                     m_longLivedPdgIds={},
                     "List of PDG IDs of long lived particles so that one can store their children");
}

~MenuTruthThinning()

DerivationFramework::MenuTruthThinning::~MenuTruthThinning ( )

virtual

Definition at line 153 of file MenuTruthThinning.cxx.

                                                         {
}

Member Function Documentation

doThinning()

StatusCode DerivationFramework::MenuTruthThinning::doThinning ( ) const

overridevirtual

Definition at line 185 of file MenuTruthThinning.cxx.

{
    const EventContext& ctx = Gaudi::Hive::currentContext();
 
    // Retrieve the truth collections
    SG::ThinningHandle<xAOD::TruthParticleContainer> importedTruthParticles
      (m_particlesKey, ctx);
    m_totpart += importedTruthParticles->size();
 
    SG::ThinningHandle<xAOD::TruthVertexContainer> importedTruthVertices
      (m_verticesKey, ctx);
 
    // Print events
    //if( m_eventCount < 20 ){
    //    printxAODTruth(m_eventCount, importedTruthParticles);
    //}
    ++m_eventCount;
    
    // Set up a mask with the same number of entries as the full TruthParticle collection
    std::vector<bool> particleMask, vertexMask;
    int nTruthParticles = importedTruthParticles->size();
    int nTruthVertices = importedTruthVertices->size();
    particleMask.assign(nTruthParticles,true); // default: keep all particles
    vertexMask.assign(nTruthVertices,false); // throw all vertices: to be discussed
    
    // Standard particle loop
    for (int particleCounter = 0; particleCounter < nTruthParticles; ++particleCounter) {
        const xAOD::TruthParticle* particle = (*importedTruthParticles)[particleCounter];
        // Keep first N particles
        if ( (particleCounter > m_writeFirstN) && (!isAccepted(particle)) ){
            // Particle removal - note mask is true by default
            particleMask[particleCounter] = false;
            ++m_removedpart;
        }
    }
    
    // If user requested preservation of descendants/ancestors:
    // - loop over the masks and work out which particles need to be descended/ascended from
    // - do the recursive loop
    // - update the masks including the descendants/ancestors
    // To ensure graph completeness, this  over-rides anything set by the special treatment
    // of taus in the section above
    DerivationFramework::DecayGraphHelper decayHelper;
    std::unordered_set<int> encounteredUniqueIDs; // For loop handling
    if (m_preserveDescendants || m_preserveGeneratorDescendants || m_preserveAncestors) {
        for (int i=0; i<nTruthParticles; ++i) {
            bool toKeep = particleMask[i];
            if (!toKeep) continue;
            const xAOD::TruthParticle* particle = (*importedTruthParticles)[i];
            encounteredUniqueIDs.clear();
            if (m_preserveDescendants) decayHelper.descendants(particle,particleMask,vertexMask,encounteredUniqueIDs,true);
            encounteredUniqueIDs.clear();
            if (m_preserveGeneratorDescendants) decayHelper.descendants(particle,particleMask,vertexMask,encounteredUniqueIDs,false);
            encounteredUniqueIDs.clear();
            if (m_preserveAncestors) decayHelper.ancestors(particle,particleMask,vertexMask,encounteredUniqueIDs);
            encounteredUniqueIDs.clear();
        }
    }
    // User only wants to keep parents, siblings and children of retained states
    // Much simpler case - no recursion so no need for uniqueIDs etc
    if (m_preserveImmediate) {
        // Make a copy of the particle mask to avoid changes being propagated
        // down the chain
        std::vector<bool> particleMaskCopy = particleMask;
        for (int i=0; i<nTruthParticles; ++i) {
            bool toKeep = particleMask[i]; // still loop over the orginal list, not the copy
            if (!toKeep) continue;
            const xAOD::TruthParticle* particle = (*importedTruthParticles)[i];
            decayHelper.immediateRelatives(particle,particleMaskCopy,vertexMask,
                                           m_preserveHadVtx); // but only update the copy
        }
        particleMask = particleMaskCopy; // Now update the original list in one go
    }
    
    // Execute the thinning service based on the mask. Finish.
    importedTruthParticles.keep (particleMask);
    importedTruthVertices.keep  (vertexMask);
    
    return StatusCode::SUCCESS;
}

finalize()

StatusCode DerivationFramework::MenuTruthThinning::finalize ( )

overridevirtual

Definition at line 176 of file MenuTruthThinning.cxx.

{
    ATH_MSG_VERBOSE("finalize() ...");
    ATH_MSG_INFO("Total of " << m_totpart << " truth particles");
    ATH_MSG_INFO("Removed "<< m_removedpart << " truth particles");
    return StatusCode::SUCCESS;
}

initialize()

StatusCode DerivationFramework::MenuTruthThinning::initialize ( )

overridevirtual

Definition at line 157 of file MenuTruthThinning.cxx.

{
    if (m_preserveDescendants && m_preserveGeneratorDescendants) {
        ATH_MSG_FATAL("You are asking to keep both all descendants, and only those from the event generator. Please check your job options.");
        return StatusCode::FAILURE;
    }
    if (m_preserveImmediate && (m_preserveDescendants || m_preserveGeneratorDescendants || m_preserveAncestors)) {
        ATH_MSG_FATAL("You are asking to preserve only parents/children/siblings of retained states, and also more distant relatives. Please check your job options.");
        return StatusCode::FAILURE;
    }
    m_totpart = 0;
    m_removedpart = 0;
    m_eventCount = 0;
 
    ATH_CHECK( m_particlesKey.initialize (m_streamName) );
    ATH_CHECK( m_verticesKey.initialize (m_streamName) );
    return StatusCode::SUCCESS;
}

isAccepted()

bool DerivationFramework::MenuTruthThinning::isAccepted

(

const xAOD::TruthParticle *

p

)

const

Definition at line 267 of file MenuTruthThinning.cxx.

{
    bool ok = false;
    
    int pdg_id = p->absPdgId();
 
    // All explicitly requested PDG IDs of long lived particles, this is needed
    // because their childrens uniqueIDs can be above the cut off m_geantOffset
    if(!m_longLivedPdgIds.empty() && parentIsLongLived(p)) ok = true; 
 
 
    if (HepMC::is_simulation_particle(p) && !m_writeGeant && !m_writeEverything && !ok) {
        if (!(MC::isPhoton(pdg_id) && m_geantPhotonPtThresh >= 0 && p->pt() > m_geantPhotonPtThresh) )
            return false;
    }
    
    // Do we want to save everything?
    if (m_writeEverything) ok = true;
    
    // Save NotPhysical particles
    if (m_writeNotPhysical && !MC::isPhysical(p)) ok = true;
    
    // OK if we select partons and are at beginning of event record
    if( m_writePartons && !MC::isHadron(pdg_id) &&
       (m_partonPtThresh<0 || p->pt()>m_partonPtThresh) )
        ok = true;
    
    //  OK if we should select hadrons and are in hadron range
    // JRC: cut changed from PHOTOSMIN to m_geantOffset
    if( m_writeHadrons && MC::isHadron (pdg_id) && !HepMC::is_simulation_particle(p) ) ok = true;
    
    // OK if we should select b hadrons and are in hadron range
    // JRC: cut changed from PHOTOSMIN to m_geantOffset
    if( m_writeBHadrons &&  !HepMC::is_simulation_particle(p) && MC::isBottomHadron (pdg_id)) ok = true;
 
    // OK if we should select c hadrons and are in hadron range
    // JRC: cut changed from PHOTOSMIN to m_geantOffset
    if( m_writeCHadrons &&  !HepMC::is_simulation_particle(p) && MC::isCharmHadron (pdg_id)) ok = true;
 
    // PHOTOS range: check whether photons come from parton range or
    // hadron range
    int motherPDGID = 999999999;
    // AV: dropped the HepMC::PHOTOSMIN condition
    if ( !HepMC::is_simulation_particle(p) && p->hasProdVtx() )
    {
        const xAOD::TruthVertex* vprod = p->prodVtx();
        if (vprod->nIncomingParticles() > 0) {
            const xAOD::TruthParticle* mother = vprod->incomingParticle(0);
            if (mother) motherPDGID = mother->pdgId();
        }
        if( m_writePartons && !MC::isHadron( motherPDGID ) ) ok = true;
        if( m_writeHadrons && MC::isHadron( motherPDGID ) ) ok = true;
    }
    
    // OK if we should select G4 particles and are in G4 range
    if( m_writeGeant && HepMC::is_simulation_particle(p) ) ok = true;
    
    if(isLeptonFromTau(p)) ok = true;
    
    if(isFsrFromLepton(p)) ok = true;
    
    // Hadronic tau decays
    if(m_writeTauHad){
        std::unordered_set<int> uniqueID_trace;
        if (isFromTau(p, uniqueID_trace))
            ok = true;
    }
    
    // Bosons
    if(m_writeBosons && (MC::isZ(p)||MC::isW(p)||MC::isHiggs(p) || (MC::isPhoton(p) && p->pt()>m_photonPtCut) )) ok = true;
    
    // BSM particles
    if(m_writeBSM && MC::isBSM(p)) ok = true;
    
    // tt+HF hadrons
    if (m_writettHFHadrons && isttHFHadron(p)) ok = true;
    
    // Top quarks
    if(m_writeTopAndDecays && MC::isTop(p)) ok = true;
    
    // All leptons
    if(m_writeAllLeptons && MC::isLepton(p))ok = true; // Include 4th generation...
 
    // All stable
    if (m_writeAllStable && MC::isStable(p) && !HepMC::is_simulation_particle(p)) ok = true;
    
    // All leptons not from hadron decays
    if(!ok && m_writeLeptonsNotFromHadrons && MC::isLepton(p) && MC::isStable(p)) {// Include 4th generation...
        ok = !(matchHadronIncTau(p) || matchQuarkIncTau(p) || isOrphanIncTau(p));
    }
    
    if ((m_writeBSMProducts || m_writeBosonProducts || m_writeTopAndDecays) && p->hasProdVtx()){ // Either way we have to go through parents
        const xAOD::TruthVertex* prodVtx = p->prodVtx();
        unsigned int nIncoming = prodVtx->nIncomingParticles();
        for(unsigned int itr=0; itr<nIncoming; ++itr) {
            const xAOD::TruthParticle* incomingParticle = prodVtx->incomingParticle(itr);
            if (!incomingParticle) continue;
            if ((m_writeBSMProducts && MC::isBSM(incomingParticle)) ||
                (m_writeBosonProducts && (MC::isZ(incomingParticle)||MC::isW(incomingParticle)||MC::isHiggs(incomingParticle) || (MC::isPhoton(incomingParticle) && incomingParticle->pt()>m_photonPtCut) )) ||
                (m_writeTopAndDecays && MC::isTop(incomingParticle)) ){
                ok = true;
                break;
            }
        }
        
    }
 
    // All explicitly requested PDG IDs
    for (const auto id : m_pdgIdsToKeep){
      if (pdg_id==id) ok=true;
    } // Loop over PDG IDs
 
    return ok;
}

isFromTau()

bool DerivationFramework::MenuTruthThinning::isFromTau	(	const xAOD::TruthParticle *	part,
		std::unordered_set< int > &	barcode_trace
	)		const

Definition at line 540 of file MenuTruthThinning.cxx.

                                                                                                    {
    
    int pdg = part->pdgId();
    
    const xAOD::TruthVertex* prod = part->prodVtx();
    if(!prod) return false; // no parent.
    
    // Simple loop catch
    if (prod==part->decayVtx()) return false;
    
    // More complex loop catch
    std::unordered_set<int>::const_iterator foundVtx = uniqueID_trace.find( HepMC::uniqueID(prod) );
    if( foundVtx != uniqueID_trace.end() ) {
      ATH_MSG_DEBUG( "Found a loop (a la Sherpa sample).  Backing out." );
      return false;
    }
    uniqueID_trace.insert(HepMC::uniqueID(prod));
    
    // Loop over the parents of this particle.
    unsigned int nIncoming = prod->nIncomingParticles();
    for (unsigned int pitr = 0; pitr<nIncoming; ++pitr){
        const xAOD::TruthParticle* itrParent = prod->incomingParticle(pitr);
        if (!itrParent) continue;
        if (pitr>2) break; // No point in trying - this vertex does not have a quantum meaning...
 
        if( MC::isTau(itrParent) ) {
            // Check if one of the children of this parent was a tau - if it is, then it is
            //   photon radiation, and we already cover that under FSR
            bool has_fsr = false;
            if ( itrParent->hasDecayVtx() ){
                unsigned int nChildren = itrParent->decayVtx()->nOutgoingParticles();
                for (unsigned int citr = 0; citr<nChildren; ++citr) {
                    const xAOD::TruthParticle* itrChild = itrParent->decayVtx()->outgoingParticle(citr);
                    if (!itrChild) continue;
                    if (MC::isTau(itrChild)){
                        has_fsr = true;
                        break;
                    } // Caught FSR check
                } // loop over immediate children
            } // Parent has an end vertex
            if (has_fsr) return false;
            ATH_MSG_DEBUG("Particle with pdgId = " << pdg << ", matched to tau");
            return true; // Has tau parent
        }
        
        // Go up the generator record until a tau is found or not.
        if(isFromTau(itrParent, uniqueID_trace)) {
            return true;
        }
    }
    
    return false;
}

isFsrFromLepton()

bool DerivationFramework::MenuTruthThinning::isFsrFromLepton

(

const xAOD::TruthParticle *

part

)

const

Definition at line 616 of file MenuTruthThinning.cxx.

                                                                                              {
    int pdg = part->pdgId();
    if(!MC::isPhoton(part)) return false;
    if(HepMC::is_simulation_particle(part) ) return false; // Geant photon
    const xAOD::TruthVertex* prod = part->prodVtx();
    if(!prod) return false; // no parent.
    // Simple loop check
    if (prod==part->decayVtx()) return false;
    // Loop over the parents of this particle.
    unsigned int nIncoming = prod->nIncomingParticles();
    for(unsigned int pitr=0; pitr<nIncoming; ++pitr){
        const xAOD::TruthParticle* itrParent = prod->incomingParticle(pitr);
        int parentId = itrParent->pdgId();
        if(MC::isElectron(parentId) ||
           MC::isMuon(parentId) ||
           MC::isTau(parentId)) {
          ATH_MSG_DEBUG("Photon with uniqueID " << HepMC::uniqueID(part) << " matched to particle with pdgId = " << parentId );
            return true; // Has lepton parent
        }
        if(parentId == pdg) { // Same particle just a different MC status
            // Go up the generator record until a tau is found or not.
            if(isFsrFromLepton(itrParent)) return true;
        }
    }
    return false;
}

isLeptonFromTau()

bool DerivationFramework::MenuTruthThinning::isLeptonFromTau

(

const xAOD::TruthParticle *

part

)

const

Definition at line 506 of file MenuTruthThinning.cxx.

                                                                                             {
    
    int pdg = part->pdgId();
    
    if(!MC::isSMLepton(part)) return false; // all leptons including tau.
    
    const xAOD::TruthVertex* prod = part->prodVtx();
    if(!prod) return false; // no parent.
    
    // Simple loop catch
    if (prod==part->decayVtx()) return false;
    
    // Loop over the parents of this particle.
    unsigned int nIncoming = prod->nIncomingParticles();
    for(unsigned int itr = 0; itr<nIncoming; ++itr){
        int parentId = prod->incomingParticle(itr)->pdgId();
        if( MC::isTau(parentId) ) {
            ATH_MSG_DEBUG("Particle with pdgId = " << pdg << ", matched to tau");
            return true; // Has tau parent
        }
        
        if(parentId == pdg) { // Same particle just a different MC status
            // Go up the generator record until a tau is found or not.
            // Note that this requires a connected *lepton* chain, while calling
            //  isFromTau would allow leptons from hadrons from taus
            if(isLeptonFromTau(prod->incomingParticle(itr))) {
                return true;
            }
        }
    }
    
    return false;
}

isOrphanIncTau()

bool DerivationFramework::MenuTruthThinning::isOrphanIncTau

(

const xAOD::TruthParticle *

part

)

const

Definition at line 398 of file MenuTruthThinning.cxx.

                                                                                             {
    //check if part descends from nothing, possibly passing through a tau decay
    //(indicating a broken truth record)
    //ported from TopFiducial
    
    int pdgId = part->pdgId();
    
    if (!part->hasProdVtx()) return true;
    
    const xAOD::TruthVertex* prodVtx = part->prodVtx();
    unsigned int nIncoming = prodVtx->nIncomingParticles();
    unsigned int itrParent = 0;
    
    // Simple loop catch
    if (prodVtx==part->decayVtx()) return false;
    
    //0 incoming-> orphan
    if(nIncoming == 0) return true;
    
    unsigned int nParents = 0;
    while(nParents==0 && itrParent<nIncoming) {
        const xAOD::TruthParticle* incomingParticle = prodVtx->incomingParticle(itrParent);
        
        // If the parent is itself migrate up the generator record
        if(incomingParticle->pdgId() == pdgId) {
            if(!isOrphanIncTau(incomingParticle)) nParents++;
        }
        else if(MC::isTau(incomingParticle)) {
            if(!isOrphanIncTau(incomingParticle)) nParents++;
        }
        else {
            nParents++;
        }
        ++itrParent;
    }
    return (nParents==0);
}

isttHFHadron()

bool DerivationFramework::MenuTruthThinning::isttHFHadron ( const xAOD::TruthParticle *  part )

static

Definition at line 597 of file MenuTruthThinning.cxx.

                                                                                     {
    
    int ttHFClassification=6;
 
    static const SG::ConstAccessor<int> TopHadronOriginFlagAcc("TopHadronOriginFlag");
    if (TopHadronOriginFlagAcc.isAvailable(*part)){
        ttHFClassification = TopHadronOriginFlagAcc(*part);
    }
    else{
        return false;
    }
    
    if (ttHFClassification < 6 ) // useful Hadrons
        return true;
    
    return false;
}

matchGenParticle()

bool DerivationFramework::MenuTruthThinning::matchGenParticle	(	const xAOD::TruthParticle *	part,
		std::vector< int > &	targetIDs,
		std::vector< int > &	intermediateIDs,
		bool	targetsAreRange
	)		const

Definition at line 436 of file MenuTruthThinning.cxx.

                                                                                                                              {
    //check if part descends from a particle, possibly passing through an intermediate decay
    //ported from TopFiducial
    
    int pdgId = part->pdgId();
    
    bool found = false;
    
    // Iterators for the target pdg or the intermediate pdg ids
    std::vector<int>::const_iterator itrPdgId, itrPdgIdEnd;
    
    if (!part->hasProdVtx()) return false;
    
    // Loop over the parents
    const xAOD::TruthVertex* prodVtx = part->prodVtx();
    unsigned int nIncoming = prodVtx->nIncomingParticles();
    
    // Simple loop catch
    if (prodVtx==part->decayVtx()) return false;
    
    unsigned int itrParent = 0;
    while(!found && itrParent<nIncoming) {
        const xAOD::TruthParticle* incomingParticle = prodVtx->incomingParticle(itrParent);
        
        if(!targetsAreRange) {
            
            // If the parent is one of the particles in the target pdg list
            itrPdgId = targetPdgIds.begin();
            itrPdgIdEnd = targetPdgIds.end();
            for(;itrPdgId != itrPdgIdEnd; ++itrPdgId) {
                if(incomingParticle->pdgId() == (*itrPdgId)) return true;
            }
        }
        else {
            int absPdgId = incomingParticle->absPdgId();
            
            // If the parent is within the range given in the target pdg list
            if(targetPdgIds.size() == 1) {
                if(absPdgId >= targetPdgIds.at(0)) return true;
            }
            else if(targetPdgIds.size() >= 2) {
                if(absPdgId >= targetPdgIds.at(0) &&
                   absPdgId <= targetPdgIds.at(1)) return true;
            }
        }
        
        // If the parent is itself migrate up the generator record
        if(incomingParticle->pdgId() == pdgId) {
            found = matchGenParticle(incomingParticle, targetPdgIds, intermediatePdgIds, targetsAreRange);
        }
        else {
            // If the parent is in the intermediatePdgIds list migrate up the generator record
            itrPdgId = intermediatePdgIds.begin();
            itrPdgIdEnd = intermediatePdgIds.end();
            bool foundIntermediate = false;
            while(!foundIntermediate && itrPdgId != itrPdgIdEnd) {
                if(incomingParticle->pdgId() == (*itrPdgId)) foundIntermediate = true;
                else ++itrPdgId;
            }
            if(foundIntermediate) {
                found = matchGenParticle(incomingParticle, targetPdgIds, intermediatePdgIds, targetsAreRange);
            }
        }
        
        ++itrParent;
    }
    return found;
}

matchHadronIncTau()

bool DerivationFramework::MenuTruthThinning::matchHadronIncTau

(

const xAOD::TruthParticle *

part

)

const

Definition at line 382 of file MenuTruthThinning.cxx.

                                                                                                {
    //check if part descends from a hadron, possibly passing through a tau decay
    //ported from TopFiducial
    std::vector<int> hadrons = std::vector<int>{111};
    std::vector<int> taus = std::vector<int>{15,-15};
    return matchGenParticle(part, hadrons, taus, true);
}

matchQuarkIncTau()

bool DerivationFramework::MenuTruthThinning::matchQuarkIncTau

(

const xAOD::TruthParticle *

part

)

const

Definition at line 390 of file MenuTruthThinning.cxx.

                                                                                               {
    //check if part descends from a hadron, possibly passing through a tau decay
    //ported from TopFiducial
    std::vector<int> quarks = std::vector<int>{1,6};
    std::vector<int> taus = std::vector<int>{15,-15};
    return matchGenParticle(part, quarks, taus, true);
}

parentIsLongLived()

bool DerivationFramework::MenuTruthThinning::parentIsLongLived

(

const xAOD::TruthParticle *

part

)

const

Definition at line 643 of file MenuTruthThinning.cxx.

                                                                                                {
  //loop over the parents of the truth particle, if the parent is in the list of long lived particles
  //store this truth particle.
  for(size_t parent_itr = 0; parent_itr < part->nParents(); parent_itr++){
    if(!part->parent(parent_itr)) continue;
    const xAOD::TruthParticle* parent = part->parent(parent_itr);
    const int parent_abs_pdgid = abs(parent->pdgId());
    if(std::find(m_longLivedPdgIds.begin(), m_longLivedPdgIds.end(), parent_abs_pdgid) != m_longLivedPdgIds.end() ){
      return true;
    }
  }
  return false;
}

Member Data Documentation

m_eventCount

std::atomic<int> DerivationFramework::MenuTruthThinning::m_eventCount {}

mutableprivate

Definition at line 138 of file MenuTruthThinning.h.

m_geantPhotonPtThresh

float DerivationFramework::MenuTruthThinning::m_geantPhotonPtThresh

private

Parameter: Write Geant photons with Pt above this threshold.

Set to < 0 to not write any.

Definition at line 77 of file MenuTruthThinning.h.

m_longLivedPdgIds

std::vector<int> DerivationFramework::MenuTruthThinning::m_longLivedPdgIds

private

Parameter: List of PDG IDs of long lived particles so that one can keep their children.

Definition at line 133 of file MenuTruthThinning.h.

m_particlesKey

SG::ThinningHandleKey<xAOD::TruthParticleContainer> DerivationFramework::MenuTruthThinning::m_particlesKey { this, "ParticlesKey", "TruthParticles", "TruthParticle container name" }

private

Definition at line 55 of file MenuTruthThinning.h.

m_partonPtThresh

float DerivationFramework::MenuTruthThinning::m_partonPtThresh

private

Parameter: Write partons with Pt above this threshold.

Definition at line 90 of file MenuTruthThinning.h.

m_pdgIdsToKeep

std::vector<int> DerivationFramework::MenuTruthThinning::m_pdgIdsToKeep

private

Parameter: List of PDG IDs to always keep.

Definition at line 130 of file MenuTruthThinning.h.

m_photonPtCut

float DerivationFramework::MenuTruthThinning::m_photonPtCut

private

Definition at line 87 of file MenuTruthThinning.h.

m_preserveAncestors

bool DerivationFramework::MenuTruthThinning::m_preserveAncestors

private

Definition at line 125 of file MenuTruthThinning.h.

m_preserveDescendants

bool DerivationFramework::MenuTruthThinning::m_preserveDescendants

private

Parameter: preserve descendant/ancestor graph completeness.

Definition at line 123 of file MenuTruthThinning.h.

m_preserveGeneratorDescendants

bool DerivationFramework::MenuTruthThinning::m_preserveGeneratorDescendants

private

Definition at line 124 of file MenuTruthThinning.h.

m_preserveHadVtx

bool DerivationFramework::MenuTruthThinning::m_preserveHadVtx

private

Definition at line 127 of file MenuTruthThinning.h.

m_preserveImmediate

bool DerivationFramework::MenuTruthThinning::m_preserveImmediate

private

Definition at line 126 of file MenuTruthThinning.h.

m_removedpart

std::atomic<unsigned int> DerivationFramework::MenuTruthThinning::m_removedpart

mutableprivate

Definition at line 137 of file MenuTruthThinning.h.

m_streamName

StringProperty DerivationFramework::MenuTruthThinning::m_streamName { this, "StreamName", "", "Name of the stream being thinned" }

private

Definition at line 53 of file MenuTruthThinning.h.

m_totpart

std::atomic<unsigned int> DerivationFramework::MenuTruthThinning::m_totpart

mutableprivate

Definition at line 136 of file MenuTruthThinning.h.

m_verticesKey

SG::ThinningHandleKey<xAOD::TruthVertexContainer> DerivationFramework::MenuTruthThinning::m_verticesKey { this, "VerticesKey", "TruthVertices", "TruthVertex container name" }

private

Definition at line 57 of file MenuTruthThinning.h.

m_writeAllLeptons

bool DerivationFramework::MenuTruthThinning::m_writeAllLeptons

private

Parameter: Write all leptons.

Definition at line 105 of file MenuTruthThinning.h.

m_writeAllStable

bool DerivationFramework::MenuTruthThinning::m_writeAllStable

private

Definition at line 111 of file MenuTruthThinning.h.

m_writeBHadrons

bool DerivationFramework::MenuTruthThinning::m_writeBHadrons

private

Parameter: Keep b-hadrons?

Definition at line 67 of file MenuTruthThinning.h.

m_writeBosonProducts

bool DerivationFramework::MenuTruthThinning::m_writeBosonProducts

private

Parameter: Write boson decay products.

Definition at line 96 of file MenuTruthThinning.h.

m_writeBosons

bool DerivationFramework::MenuTruthThinning::m_writeBosons

private

Parameter: Keep bosons?

Definition at line 86 of file MenuTruthThinning.h.

m_writeBSM

bool DerivationFramework::MenuTruthThinning::m_writeBSM

private

Parameter: Keep BSM particles?

Definition at line 83 of file MenuTruthThinning.h.

m_writeBSMProducts

bool DerivationFramework::MenuTruthThinning::m_writeBSMProducts

private

Parameter: Write BSM decay products.

Definition at line 99 of file MenuTruthThinning.h.

m_writeCHadrons

bool DerivationFramework::MenuTruthThinning::m_writeCHadrons

private

Parameter: Keep c-hadrons?

Definition at line 70 of file MenuTruthThinning.h.

m_writeEverything

bool DerivationFramework::MenuTruthThinning::m_writeEverything

private

Parameter: Write absolutely everything.

Definition at line 93 of file MenuTruthThinning.h.

m_writeFirstN

int DerivationFramework::MenuTruthThinning::m_writeFirstN

private

Parameter: First N particles to write.

Definition at line 120 of file MenuTruthThinning.h.

m_writeGeant

bool DerivationFramework::MenuTruthThinning::m_writeGeant

private

Parameter: Keep geant particles?

Definition at line 73 of file MenuTruthThinning.h.

m_writeHadrons

bool DerivationFramework::MenuTruthThinning::m_writeHadrons

private

Parameter: Keep hadrons?

Definition at line 64 of file MenuTruthThinning.h.

m_writeLeptonsNotFromHadrons

bool DerivationFramework::MenuTruthThinning::m_writeLeptonsNotFromHadrons

private

Parameter: Write all leptons.

Definition at line 108 of file MenuTruthThinning.h.

m_writeNotPhysical

bool DerivationFramework::MenuTruthThinning::m_writeNotPhysical

private

Parameter: Write particles with status code 3.

Definition at line 114 of file MenuTruthThinning.h.

m_writePartons

bool DerivationFramework::MenuTruthThinning::m_writePartons

private

Parameter: Keep partons?

Definition at line 61 of file MenuTruthThinning.h.

m_writeTauHad

bool DerivationFramework::MenuTruthThinning::m_writeTauHad

private

Parameter: Keep hadronic tau decays?

Definition at line 80 of file MenuTruthThinning.h.

m_writeTopAndDecays

bool DerivationFramework::MenuTruthThinning::m_writeTopAndDecays

private

Parameter: Write top and decay products.

Definition at line 102 of file MenuTruthThinning.h.

m_writettHFHadrons

bool DerivationFramework::MenuTruthThinning::m_writettHFHadrons

private

Parameter: Write particles for tt+HF classification.

Definition at line 117 of file MenuTruthThinning.h.

---

The documentation for this class was generated from the following files:

• MenuTruthThinning.h
• MenuTruthThinning.cxx