#include <TtresdRmin.h>

Collaboration diagram for TtresdRmin:

Public Member Functions
	TtresdRmin (TString units="GeV", bool verbose=false)
	Class constructor. More...

virtual	~TtresdRmin ()
	Class destructor. More...

void	Init (double highJetMass=60.0, double minJetPt=25.0, double dRCutOffset=2.5, double dRCutMassFactor=0.015)
	Initialize the tool, giving specific parameters needed for dRmin cut Paramters are given in GeV scale, will be scaled internally to MeV if set. More...

bool	SetupEvent (vector< TLorentzVector * > &jets, TLorentzVector &lepton, TLorentzVector &neutrino)
	Set the objects needed to calculate the invariant ttbar mass. More...

vector< int >	GetJetIndices ()
	Get jet indicices identified by dRmin method, used for mtt reconstruction. More...

int	GetHighMassJetIndex ()
	Get index of jet with highest mass. More...

double	GetHighestJetMass ()
	Get highest jet mass found in jet vector. More...

void	SetHighMThreshold (double highJetMass)
	Set threshold for highMass case. More...

double	GetMttdRmin ()
	Get calculated invariant ttbar mass for dRmin method. More...

TLorentzVector	GetTlvdRmin ()
	Get the TLorentzVector for the ttbar system. More...

Private Member Functions
bool	FindJetIndices_HM ()
	Find the jet indices used for mtt reconstruction for events containing a high mass jet. More...

bool	FindJetIndices_LM ()
	Find the jet indices used for mtt reconstruction for events NOT containing a high mass jet. More...

Private Attributes
bool	m_init
	Flag of initialization status of tool. More...

bool	m_verbose
	Flag for verbosity level. More...

double	m_highM
	Jet mass threshold to identify it as a "high mass jet". More...

double	m_Units
	Unit factor used, GeV = 1, MeV = 1000. More...

double	m_minJetPt
	Minimum jet pT used for highM jet search. More...

double	m_dRCutOffset
	Offset of dRmin specific cut on jet_dR. More...

double	m_dRCutMassFactor
	Jet mass factor of dRmin specific cut on jet_dR. More...

int	m_highMassJetIndex
	Index of highest mass jet for current event. More...

double	m_highestJetMass
	Highest jet mass found for current event. More...

vector< TLorentzVector * >	m_jets
	Input vector of jet objects. More...

TLorentzVector	m_lepton
	Input lepton object. More...

TLorentzVector	m_neutrino
	Input neutrino object. More...

TLorentzVector	m_Wleptonic
	leptonic W object reconstructed from lepton and neutrino More...

vector< int >	m_jet_indices
	Indices of jets used for reconstruction of invariant ttbar mass. More...

Detailed Description

Definition at line 16 of file TtresdRmin.h.

Constructor & Destructor Documentation

◆ TtresdRmin()

TtresdRmin::TtresdRmin	(	TString	units = `"GeV"`,
		bool	verbose = `false`
	)

Class constructor.

Set unit used in four momentum vector calculations

Parameters

units : Unit used for reconstruction, MeV or GeV

Definition at line 35 of file TtresdRmin.cxx.

                                                   {
   m_init = false;
   m_verbose = verbose;
  
   if (units == "MeV") {
     m_Units = 1000.0;
   } else if (units == "GeV") {
     m_Units = 1.0;
   } else {
     cout << " ERROR in in TtresdRmin :: Unit setting '" << units << "' doest not match pattern of GeV or MeV" << endl;
   }
 }

◆ ~TtresdRmin()

virtual TtresdRmin::~TtresdRmin ( )

inlinevirtual

Class destructor.

Definition at line 31 of file TtresdRmin.h.

31 {};

Member Function Documentation

◆ FindJetIndices_HM()

bool TtresdRmin::FindJetIndices_HM ( )

private

Find the jet indices used for mtt reconstruction for events containing a high mass jet.

Returns: bool : success of function

Definition at line 153 of file TtresdRmin.cxx.

                                    {
   //===== Find closest jet to high mass jet =====
   double dR_had = 999.9;
   int closest_jetIdx_had = -1;
  
   //Loop m_jets
   for (int ijet = 0; ijet < (int) m_jets.size(); ijet++) {
     if (ijet == m_highMassJetIndex) continue; //next jet
     double dR = m_jets.at(m_highMassJetIndex)->DeltaR(*(m_jets.at(ijet)));
     if (dR < dR_had) {
       closest_jetIdx_had = ijet;
       dR_had = dR;
     }
   }//end of jet loop
  
  
   //===== Find closest jet to leptonic W-boson =====
   double dR_lep = 999.9;
   int closest_jetIdx_lep = -1;
  
   //Loop jets
   for (int ijet = 0; ijet < (int) m_jets.size(); ijet++) {
     if (ijet == m_highMassJetIndex) continue; //next jet
     if (ijet == closest_jetIdx_had) continue; //next jet
     double dR = m_Wleptonic.DeltaR(*(m_jets.at(ijet)));
     if (dR < dR_lep) {
       closest_jetIdx_lep = ijet;
       dR_lep = dR;
     }
   }//end of jet loop
  
   //===== Save jet indices used for reconstruction =====
   m_jet_indices.push_back(m_highMassJetIndex);
  
   if (closest_jetIdx_had >= 0) {
     m_jet_indices.push_back(closest_jetIdx_had);
   } else {
     if (m_verbose) cout << " WARNING in TtresdRmin :: Can not find jet close to high mass jet" << endl;
   }
  
   if (closest_jetIdx_lep >= 0) {
     m_jet_indices.push_back(closest_jetIdx_lep);
   } else {
     if (m_verbose) cout << " WARNING in TtresdRmin :: Can not find jet close to leptonic W-boson" << endl;
   }
  
   return true;
 }

◆ FindJetIndices_LM()

bool TtresdRmin::FindJetIndices_LM ( )

private

Find the jet indices used for mtt reconstruction for events NOT containing a high mass jet.

Returns: bool : success of function

Definition at line 205 of file TtresdRmin.cxx.

                                    {
   //jetflag vector to flag jets rejected by dRmin method
   vector<bool> jetflags(m_jets.size(), true);
  
   //do dRmin procedure if at least 4 jets are available
   bool repeat_dRmin = ((int) m_jets.size() > 3);
  
   //Loop dRmin method
   while (repeat_dRmin) {
     repeat_dRmin = false;
  
     double dR_max = 0.0;  //Maximum dR passing dRmin cut requirement
     int index_max = -1;   //Index of jet with dR_max
  
     //Number of goodflagged processed jets (outer loop)
     int nijet = 0;
  
     //Outer jet loop
     for (int ijet = 0; ijet < (int) m_jets.size(); ijet++) {
       if (!jetflags.at(ijet)) continue; //next outer jet
  
       //Use maximum of 4 goodflagged jets
       if (nijet > 3) break;
       nijet++;
  
       double dR_min = 999.9;
  
       double dR_lep = m_jets.at(ijet)->DeltaR(m_lepton);
       dR_min = dR_lep;
  
       //Number of goodflagged processed jets (inner loop)
       int njjet = 0;
  
       //Inner jet loop
       for (int jjet = 0; jjet < (int) m_jets.size(); jjet++) {
         if (!jetflags.at(jjet)) continue; //next inner jet
  
         //Use maximum of 4 goodflagged jets
         if (njjet > 3) break;
         njjet++;
  
         if (jjet == ijet) continue; //next inner jet
  
         double dR_jet = m_jets.at(ijet)->DeltaR(*(m_jets.at(jjet)));
         if (dR_jet < dR_min) {
           dR_min = dR_jet;
         }
       } //end of inner jet loop
  
       if (dR_min > dR_max) {
         //dRmin specific cut (defined for GeV, needs m_Units correction on mass)
         double dR_cut = m_dRCutOffset - m_dRCutMassFactor * m_jets.at(ijet)->M() / m_Units;
         if (dR_min > dR_cut) {
           dR_max = dR_min;
           index_max = ijet;
         }
       }
     }//end of outer jet loop
  
     if (index_max >= 0) {
       jetflags.at(index_max) = false;
       if (count(jetflags.begin(), jetflags.end(), true) > 3) repeat_dRmin = true;
     }
   }//end of while(repeat)
  
   int nijet = 0;
   //Save selected jet indices
   for (int ijet = 0; ijet < (int) m_jets.size(); ijet++) {
     if (!jetflags.at(ijet)) continue; //next jet
  
     //Use maximum of 4 goodflagged jets
     if (nijet > 3) break;
     nijet++;
  
     m_jet_indices.push_back(ijet);
   }
  
   return true;
 }

◆ GetHighestJetMass()

double TtresdRmin::GetHighestJetMass ( )

inline

Get highest jet mass found in jet vector.

Returns: double : jet mass

Definition at line 77 of file TtresdRmin.h.

77 {return m_highestJetMass;};

◆ GetHighMassJetIndex()

int TtresdRmin::GetHighMassJetIndex ( )

inline

Get index of jet with highest mass.

Returns: int : jet index

Definition at line 70 of file TtresdRmin.h.

70 {return m_highMassJetIndex;};

◆ GetJetIndices()

vector<int> TtresdRmin::GetJetIndices ( )

inline

Get jet indicices identified by dRmin method, used for mtt reconstruction.

Returns: vector<int> : Indices of identified jets in passed jet vector

Definition at line 63 of file TtresdRmin.h.

63 {return m_jet_indices;};

◆ GetMttdRmin()

double TtresdRmin::GetMttdRmin ( )

Get calculated invariant ttbar mass for dRmin method.

Returns: double : invariant ttbar mass

Definition at line 119 of file TtresdRmin.cxx.

                                {
   double mtt_dRmin = 0.0;
  
   if (m_jet_indices.size() > 0) {
     TLorentzVector dRsumJet(0, 0, 0, 0);
     for (unsigned int ijet = 0; ijet < m_jet_indices.size(); ijet++) {
       dRsumJet += *(m_jets.at(m_jet_indices.at(ijet)));
     }
     mtt_dRmin = (dRsumJet + m_Wleptonic).M();
   }
  
   return mtt_dRmin;
 }

◆ GetTlvdRmin()

TLorentzVector TtresdRmin::GetTlvdRmin ( )

Get the TLorentzVector for the ttbar system.

Definition at line 136 of file TtresdRmin.cxx.

                                        {
   TLorentzVector tt_dRmin(0, 0, 0, 0);
  
   if (m_jet_indices.size() > 0) {
     TLorentzVector dRsumJet(0, 0, 0, 0);
     for (unsigned int ijet = 0; ijet < (unsigned int) m_jet_indices.size(); ijet++) {
       dRsumJet += *(m_jets.at(m_jet_indices.at(ijet)));
     }
     tt_dRmin = dRsumJet + m_Wleptonic;
   }
  
   return tt_dRmin;
 }

◆ Init()

void TtresdRmin::Init	(	double	highJetMass = `60.0`,
		double	minJetPt = `25.0`,
		double	dRCutOffset = `2.5`,
		double	dRCutMassFactor = `0.015`
	)

Initialize the tool, giving specific parameters needed for dRmin cut Paramters are given in GeV scale, will be scaled internally to MeV if set.

Parameters

highJetMass	: jet mass for highMass jet criteria
minJetPt	: minimum pT of jets used for dRmin method
dRCutOffset	: Offset for dRmin specific cut
dRCutMassFactor	: Factor applied to jet mass in dRmin specific cut

Returns: void

Definition at line 51 of file TtresdRmin.cxx.

                                                                                                      {
   //Set parameters
   m_highM = highJetMass;
   m_minJetPt = minJetPt;
   m_dRCutOffset = dRCutOffset;
   m_dRCutMassFactor = dRCutMassFactor;
  
   //Reset
   m_jet_indices.clear();
   m_highMassJetIndex = -1;
   m_highestJetMass = 0.0;
  
   //Initialization status
   m_init = true;
 }

◆ SetHighMThreshold()

void TtresdRmin::SetHighMThreshold ( double highJetMass )

inline

Set threshold for highMass case.

Parameters

highJetMass threshold for the high jet mass case

Definition at line 84 of file TtresdRmin.h.

84 {m_highM = highJetMass;};

◆ SetupEvent()

bool TtresdRmin::SetupEvent	(	vector< TLorentzVector * > &	jets,
		TLorentzVector &	lepton,
		TLorentzVector &	neutrino
	)

Set the objects needed to calculate the invariant ttbar mass.

Parameters

jets	: Jet container
lepton	: Selected lepton
neutrino	: Reconstructed neutrino

Returns: bool : Success of setup

Definition at line 70 of file TtresdRmin.cxx.

                                                                                                            {
   if (!m_init) {
     cout << " ERROR in TtresdRmin :: Tool not initialized! Initialize tool before usage!!!" << endl;
     return false;
   }
  
   //Set objects
   m_jets = jets;
   m_lepton = lepton;
   m_neutrino = neutrino;
   m_Wleptonic = lepton + neutrino;
  
   //Reset
   m_jet_indices.clear();
   m_highMassJetIndex = -1;
   m_highestJetMass = 0.0;
  
  
   //===== Search for high mass jet =====
   int highMJetIdx = -1;                  //Index of jet with highest mass above threshold
   double highM = m_highM * m_Units;     //Threshold for high mass jet
   double minJetPt = m_minJetPt * m_Units;  //Threshold for jet pT
  
   //Loop jets
   for (int ijet = 0; ijet < (int) jets.size(); ijet++) {
     //pT requirmenet
     if (jets.at(ijet)->Pt() <= minJetPt) continue; //next jet
  
     //high mass requirement
     if (jets.at(ijet)->M() <= highM) continue; //next jet
  
     highMJetIdx = ijet;
     highM = jets.at(ijet)->M();
   }//end of jet loop
  
   if (highMJetIdx >= 0) {
     m_highMassJetIndex = highMJetIdx;
     m_highestJetMass = highM;
     FindJetIndices_HM();
   } else {
     FindJetIndices_LM();
   }
  
   return true;
 }

Member Data Documentation

◆ m_dRCutMassFactor

double TtresdRmin::m_dRCutMassFactor

private

Jet mass factor of dRmin specific cut on jet_dR.

Definition at line 122 of file TtresdRmin.h.

◆ m_dRCutOffset

double TtresdRmin::m_dRCutOffset

private

Offset of dRmin specific cut on jet_dR.

Definition at line 121 of file TtresdRmin.h.

◆ m_highestJetMass

double TtresdRmin::m_highestJetMass

private

Highest jet mass found for current event.

Definition at line 125 of file TtresdRmin.h.

◆ m_highM

double TtresdRmin::m_highM

private

Jet mass threshold to identify it as a "high mass jet".

Definition at line 118 of file TtresdRmin.h.

◆ m_highMassJetIndex

int TtresdRmin::m_highMassJetIndex

private

Index of highest mass jet for current event.

Definition at line 124 of file TtresdRmin.h.

◆ m_init

bool TtresdRmin::m_init

private

Flag of initialization status of tool.

Definition at line 116 of file TtresdRmin.h.

◆ m_jet_indices

vector<int> TtresdRmin::m_jet_indices

private

Indices of jets used for reconstruction of invariant ttbar mass.

Definition at line 132 of file TtresdRmin.h.

◆ m_jets

vector<TLorentzVector*> TtresdRmin::m_jets

private

Input vector of jet objects.

Definition at line 127 of file TtresdRmin.h.

◆ m_lepton

TLorentzVector TtresdRmin::m_lepton

private

Input lepton object.

Definition at line 128 of file TtresdRmin.h.

◆ m_minJetPt

double TtresdRmin::m_minJetPt

private

Minimum jet pT used for highM jet search.

Definition at line 120 of file TtresdRmin.h.

◆ m_neutrino

TLorentzVector TtresdRmin::m_neutrino

private

Input neutrino object.

Definition at line 129 of file TtresdRmin.h.

◆ m_Units

double TtresdRmin::m_Units

private

Unit factor used, GeV = 1, MeV = 1000.

Definition at line 119 of file TtresdRmin.h.

◆ m_verbose

bool TtresdRmin::m_verbose

private

Flag for verbosity level.

Definition at line 117 of file TtresdRmin.h.

◆ m_Wleptonic

TLorentzVector TtresdRmin::m_Wleptonic

private

leptonic W object reconstructed from lepton and neutrino

Definition at line 130 of file TtresdRmin.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TtresdRmin()

◆ ~TtresdRmin()

Member Function Documentation

◆ FindJetIndices_HM()

◆ FindJetIndices_LM()

◆ GetHighestJetMass()

◆ GetHighMassJetIndex()

◆ GetJetIndices()

◆ GetMttdRmin()

◆ GetTlvdRmin()

◆ Init()

◆ SetHighMThreshold()

◆ SetupEvent()

Member Data Documentation

◆ m_dRCutMassFactor

◆ m_dRCutOffset

◆ m_highestJetMass

◆ m_highM

◆ m_highMassJetIndex

◆ m_init

◆ m_jet_indices

◆ m_jets

◆ m_lepton

◆ m_minJetPt

◆ m_neutrino

◆ m_Units

◆ m_verbose

◆ m_Wleptonic