TrigBjetMonitorAlgorithm.cxx File Reference

#include "TrigBjetMonitorAlgorithm.h"
#include "EventPrimitives/EventPrimitivesHelpers.h"
#include "TrigCompositeUtils/TrigCompositeUtils.h"
#include "TrigDecisionTool/TrigDecisionTool.h"
#include "xAODJet/JetContainer.h"
#include "xAODBTagging/BTaggingContainer.h"

Include dependency graph for TrigBjetMonitorAlgorithm.cxx:

Go to the source code of this file.

Functions
bool	LLR (double pu, double pc, double pb, double &w)
bool	LLRW (float pqcd, float ptop, float phbb, float &w)
bool	CalcRelPt (float muonPt, float muonEta, float muonPhi, float jetPt, float jetEta, float jetPhi, float &RelPt)
float	phiCorr (float phi)

Function Documentation

CalcRelPt()

bool CalcRelPt	(	float	muonPt,
		float	muonEta,
		float	muonPhi,
		float	jetPt,
		float	jetEta,
		float	jetPhi,
		float &	RelPt
	)

Definition at line 63 of file TrigBjetMonitorAlgorithm.cxx.

                                                                                                                   {
 
    bool r = false;
    RelPt = -20.;
 
    float muonT, muonX, muonY, muonZ, muon, jetT, jetX, jetY, jetZ, jet, scprod;
 
    muonT = 2.*atan( exp(-muonEta) );
    jetT = 2.*atan( exp(-jetEta) );
    if ( (std::abs(muonT) > 0.) && (std::abs(jetT) > 0.) ) {
        muon = muonPt/std::abs( sin(muonT) );
        muonX = muonPt*cos(muonPhi);
        muonY = muonPt*sin(muonPhi);
        muonZ = muon*cos(muonT);
        jet = jetPt/std::abs( sin(jetT) );
        jetX = jetPt*cos(jetPhi);
        jetY = jetPt*sin(jetPhi);
        jetZ = jet*cos(jetT);
        scprod = (muonX*jetX + muonY*jetY + muonZ*jetZ)/(muon*jet);
        scprod *= scprod;
        if ( (1. - scprod) > 0. ) {
            RelPt = muon * sqrt(1. - scprod);
            r = true;
        }
    }
 
    return r;
 
}

LLR()

bool LLR	(	double	pu,
		double	pc,
		double	pb,
		double &	w
	)

Definition at line 34 of file TrigBjetMonitorAlgorithm.cxx.

                                                      {
    w = -100.;
    bool ll = false;
    float cfrac(0.018); // DG 2022/07/28
    if (pb > 0.) {
        double denom = pu*(1.-cfrac)+pc*cfrac;
        if (denom > 0.) {
            w = log(pb/denom);
            ll = true;
        }
    }
    return ll;
}

LLRW()

bool LLRW	(	float	pqcd,
		float	ptop,
		float	phbb,
		float &	w
	)

Definition at line 48 of file TrigBjetMonitorAlgorithm.cxx.

                                                         { // RJ 17/02/2025
    w = -100.;
    bool ll = false;
    float topfrac(0.25);
    if (phbb > 0.) {
        double denom = pqcd*(1.-topfrac)+ptop*topfrac;
        if (denom > 0.) {
            w = log(phbb/denom);
            ll = true;
        }
    }
    return ll;
}

phiCorr()

float phiCorr

(

float

phi

)

Definition at line 93 of file TrigBjetMonitorAlgorithm.cxx.

                         {
    if (phi < -M_PI) phi += 2*M_PI;
    if (phi >  M_PI) phi -= 2*M_PI;
    return phi;
}