TrigBjetMonitorAlgorithm.cxx File Reference

#include "TrigBjetMonitorAlgorithm.h"
#include "EventPrimitives/EventPrimitivesHelpers.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	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 45 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 29 of file TrigBjetMonitorAlgorithm.cxx.

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

phiCorr()

float phiCorr

(

float

phi

)

Definition at line 75 of file TrigBjetMonitorAlgorithm.cxx.

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