MSVtxValidationAlgUtils Namespace Reference

Classes
struct	ActiveVertex
struct	VtxIso

Functions
double	getCTau (const xAOD::TruthVertex *decVtx)
double	getCalEnergyLogRatio (double EMF)
bool	comparePt (const xAOD::TruthParticle *part1, const xAOD::TruthParticle *part2)
std::vector< const xAOD::TruthParticle * >	getChildren (const xAOD::TruthParticle *mother)
std::vector< const xAOD::TruthParticle * >	getStableChildrenRecursive (const xAOD::TruthParticle *particle, bool findOnlyGenStable, std::unordered_set< const xAOD::TruthParticle * > &visited)
std::vector< const xAOD::TruthParticle * >	getStableChildren (const xAOD::TruthParticle *particle, bool findOnlyGenStable)
std::vector< const xAOD::TruthParticle * >	getDecayProducts (const xAOD::TruthVertex *vtx)
std::vector< ActiveVertex >	getActiveVertices (const xAOD::Jet *jet, const xAOD::TruthParticleContainer &truthParticles)
VtxIso	getIso (const xAOD::Vertex *vtx, const xAOD::TrackParticleContainer &Tracks, const xAOD::JetContainer &Jets, double trackIso_pT, double softTrackIso_R, double jetIso_pT, double jetIso_LogRatio)

Function Documentation

comparePt()

bool MSVtxValidationAlgUtils::comparePt	(	const xAOD::TruthParticle *	part1,
		const xAOD::TruthParticle *	part2
	)

Definition at line 33 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                              {
    return (part1->pt()>part2->pt());
}

getActiveVertices()

std::vector< ActiveVertex > MSVtxValidationAlgUtils::getActiveVertices	(	const xAOD::Jet *	jet,
		const xAOD::TruthParticleContainer &	truthParticles
	)

Definition at line 107 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                                               {
    // Finds active vertices in the ancestry tree of the Truth particles close to the jet 
    // Production vertices passing the minimal activity requirements are stored as ActiveVertex objects 
    
    std::set<const xAOD::TruthVertex*> seenVertices; // vertices already seen in the ancestry tree
    std::vector<ActiveVertex> actVertices; // holds all the active vertices associated with the jet
    std::vector<const xAOD::TruthParticle*> children; // decay products of the current vertex
    size_t nChildren{0};        
 
    for (const xAOD::TruthParticle* tp : truthParticles){
        if (!tp || !tp->isStable() || jet->p4().DeltaR(tp->p4()) > 0.4) continue; // only final state particles close to the jet axis will pass
        if (!HepMC::is_simulation_particle(tp)) continue; // only simulated particles will pass.
 
        size_t decayDepth{0};
        const xAOD::TruthParticle* current = tp;
        while (current) {
            if (decayDepth > 200) break; // safety break 
            // prevent loop from going too deep where the truth record contains information used for generator internal book keeping 
            // add check on HepMC::is_simulation_particle(current) to limit to GEANT4 layer of ancestry
            if (!MC::isPhysical(current)) break;
            const xAOD::TruthVertex* prodVtx = current->prodVtx();
            if (!prodVtx || seenVertices.count(prodVtx)) break;  // No more ancestry or already visited. 
            if (prodVtx->v4().Mag2() < 0) break; // minimal requirements on the vertex: physical spacetime interval
            seenVertices.insert(prodVtx);
            decayDepth++;
 
            children = getDecayProducts(prodVtx);
            nChildren = children.size();
            // requirements for an active vertex
            if (nChildren >= 2) {
                TLorentzVector vtx4Vec;
                float scalarPtSum{0.};
                for (const xAOD::TruthParticle* child : children) {
                    vtx4Vec += child->p4();
                    scalarPtSum += child->p4().Pt();
                }
                actVertices.push_back({prodVtx, vtx4Vec.E(), vtx4Vec.M(), vtx4Vec.Pt(), scalarPtSum, nChildren, decayDepth});
            }
            current = prodVtx->incomingParticle(0); // move up one level in the ancestor tree
        }
    }
 
    return actVertices;
}

getCalEnergyLogRatio()

double MSVtxValidationAlgUtils::getCalEnergyLogRatio

(

double

EMF

)

Definition at line 16 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                       {
    double logRatio{-99.};
    double zero{0.};
    double one{1.};
 
    if (CxxUtils::fpcompare::greater(EMF,zero)){
        if(CxxUtils::fpcompare::greater_equal(EMF,one)) logRatio = -999.;
        else logRatio = std::log10(double(1./EMF - 1.));
    } 
    else {
        logRatio = 999;
    }
 
    return logRatio;
}

getChildren()

std::vector< const xAOD::TruthParticle * > MSVtxValidationAlgUtils::getChildren

(

const xAOD::TruthParticle *

mother

)

Definition at line 38 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                {
    // returns pointers to children of a given particle
    if (!mother) return {};
 
    std::vector<const xAOD::TruthParticle*> children{};
    for (size_t i=0; i<mother->nChildren(); ++i) {
        const xAOD::TruthParticle* child = mother->child(i);
        // avoid infinite loop
        if (!child || child == mother) continue;
        children.push_back(child);
    }
 
    return children;
}

getCTau()

double MSVtxValidationAlgUtils::getCTau

(

const xAOD::TruthVertex *

decVtx

)

Definition at line 10 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                             {
    // compute the ctau for a given decay vertex via: ctau = |r|/gamma
    return decVtx->v4().Vect().Mag() / decVtx->genvecV4().Gamma();
}

getDecayProducts()

std::vector< const xAOD::TruthParticle * > MSVtxValidationAlgUtils::getDecayProducts

(

const xAOD::TruthVertex *

vtx

)

Definition at line 91 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                {
    // Returns the decay products of a given vertex. Use this instead of vtx->particles_out() to remove null pointers. 
    std::vector<const xAOD::TruthParticle*> decayProducts;
    if (!vtx) return decayProducts;
 
    for (size_t i = 0; i < vtx->nOutgoingParticles(); ++i) {
        const xAOD::TruthParticle* tp = vtx->outgoingParticle(i);
        if (!tp) continue;
        decayProducts.push_back(tp);
    }
 
    return decayProducts;
}

getIso()

VtxIso MSVtxValidationAlgUtils::getIso	(	const xAOD::Vertex *	vtx,
		const xAOD::TrackParticleContainer &	Tracks,
		const xAOD::JetContainer &	Jets,
		double	trackIso_pT,
		double	softTrackIso_R,
		double	jetIso_pT,
		double	jetIso_LogRatio
	)

Definition at line 153 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                                  {
    // compute the isolation metrics of the MS vertex: 
    // - delta R to closest hard track
    // - delta R to closest punch-through candidate jet
    // - sum of soft track pT in a cone around the vertex 
 
    VtxIso iso{};
    const Amg::Vector3D vtx_pos = vtx->position();
 
    // isolation towards tracks 
    Amg::Vector3D softTrack_pTsum{Amg::Vector3D::Zero()};
    double hardTrack_mindR{99.}; 
    for(const xAOD::TrackParticle* Track : Tracks){
        double dR = xAOD::P4Helpers::deltaR(vtx_pos.eta(), vtx_pos.phi(), Track->eta(), Track->phi());
        // hard tracks
        if(Track->pt() >= trackIso_pT && dR < hardTrack_mindR) hardTrack_mindR = dR; 
        // soft tracks
        if(Track->pt() < trackIso_pT && dR < softTrackIso_R) softTrack_pTsum += Amg::Vector3D(Track->p4()[0], Track->p4()[1], Track->p4()[2]);
    }
 
    iso.track_mindR = hardTrack_mindR != 99. ? hardTrack_mindR : -1.;
    iso.track_pTsum = softTrack_pTsum.mag() != 0. ? softTrack_pTsum.perp()/Gaudi::Units::GeV : -1.;
 
    // isolation towards jets
    double jet_mindR{99.}; 
    for(const xAOD::Jet* Jet : Jets){
        if(Jet->pt() < jetIso_pT) continue;
        // if(!Jet->getAttribute<char>("passJVT")) continue; 
        double logratio = getCalEnergyLogRatio(Jet->getAttribute<float>("EMFrac"));
        if(logratio >= jetIso_LogRatio) continue;
 
        double dR = xAOD::P4Helpers::deltaR(vtx_pos.eta(), vtx_pos.phi(), Jet->eta(), Jet->phi());
        if(dR < jet_mindR){
            jet_mindR = dR;
        } 
    }
    iso.jet_mindR = jet_mindR != 99. ? jet_mindR : -1.;
 
    return iso;
}

getStableChildren()

std::vector< const xAOD::TruthParticle * > MSVtxValidationAlgUtils::getStableChildren	(	const xAOD::TruthParticle *	particle,
		bool	findOnlyGenStable
	)

Definition at line 83 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                                                {
    // Finds the stable decay products of a given particle. Can either find all stable particles or only those that are generator stable only.
    // Interface to the recursive function that traverses the decay chain of the particle.
    std::unordered_set<const xAOD::TruthParticle*> visited; // keeps track of visited particles to avoid infinite loops in decay chains
    return getStableChildrenRecursive(particle, findOnlyGenStable, visited);
}

getStableChildrenRecursive()

std::vector< const xAOD::TruthParticle * > MSVtxValidationAlgUtils::getStableChildrenRecursive	(	const xAOD::TruthParticle *	particle,
		bool	findOnlyGenStable,
		std::unordered_set< const xAOD::TruthParticle * > &	visited
	)

Definition at line 54 of file MuonSpectrometer/MuonValidation/MuonVertexValidation/src/Utils.cxx.

                                                                                                                                                                              {
    // Can either find all stable particles or only those that are generator stable only.
    
    std::vector<const xAOD::TruthParticle*> stableChildren;
    // Base case: if the particle is null or has already been visited, return empty vector
    if (!particle || visited.count(particle)) return stableChildren;
    visited.insert(particle);
 
    // Return the particle if it is stable
    bool particleIsStable = findOnlyGenStable ? particle->isGenStable() : particle->isStable();
    if (particleIsStable) {
        stableChildren.push_back(particle);
        return stableChildren;
    }
 
    // Recursive case: get the children of the particle and traverse their decay chains
    std::vector<const xAOD::TruthParticle*> children = getChildren(particle);
    for (const xAOD::TruthParticle* child : children) {
        if (!child) continue;
        std::vector<const xAOD::TruthParticle*> grandChildren = getStableChildrenRecursive(child, findOnlyGenStable, visited);
        for (const xAOD::TruthParticle* c : grandChildren) {
            if (std::none_of(stableChildren.begin(), stableChildren.end(), [&](const auto& x) {return x==c;})) stableChildren.push_back(c);
        }
    }
 
    return stableChildren;
}