Helper class for the calculation of a top and a W boson mass estime. More...

#include "DQM/Physics/interface/TopDQMHelpers.h"

Public Member Functions
	Calculate (int maxNJets, double wMass)
	default constructor More...

double	massBTopQuark (const std::vector< reco::Jet > &jets, std::vector< double > VbtagWP, double btagWP_)
	calculate b-tagged t-quark mass estimate More...

double	massTopQuark (const std::vector< reco::Jet > &jets)
	calculate t-quark mass estimate More...

double	massWBoson (const std::vector< reco::Jet > &jets)
	calculate W boson mass estimate More...

double	tmassTopQuark (reco::RecoCandidate *lep, const reco::MET &met, const reco::Jet &b)
	calculate top quark transverse mass estimate More...

double	tmassWBoson (reco::RecoCandidate *lep, const reco::MET &met, const reco::Jet &b)
	calculate W boson transverse mass estimate More...

	~Calculate ()
	default destructor More...

Private Member Functions
void	operator() (const reco::Jet &bJet, reco::RecoCandidate *lepton, const reco::MET &met)
	do the calculation of the transverse top and W masses More...

void	operator() (const std::vector< reco::Jet > &jets)

void	operator2 (const std::vector< reco::Jet > &, std::vector< double >, double)
	do the calculation of the t-quark mass with one b-jet More...

Private Attributes
bool	failed_
	indicate failed associations More...

double	massBTopQuark_
	cache of b-tagged top quark mass estimate More...

double	massTopQuark_
	cache of top quark mass estimate More...

double	massWBoson_
	cache of w boson mass estimate More...

int	maxNJets_
	max. number of jets to be considered More...

double	tmassTopQuark_
	cache of top quark transverse mass estimate More...

double	tmassWBoson_
	cache of W boson transverse mass estimate More...

double	wMass_
	paramater of the w boson mass More...

Detailed Description

Helper class for the calculation of a top and a W boson mass estime.

Helper class for the calculation of a top and a W boson mass estimate. The core implementation originates from the plugin TtSemiLepHypMaxSumPtWMass in TopQuarkAnalysis/TopJetCombination package. It may be extended to include b tag information.

Definition at line 127 of file TopDQMHelpers.h.

Constructor & Destructor Documentation

◆ Calculate()

Calculate::Calculate	(	int	maxNJets,
		double	wMass
	)

default constructor

Definition at line 3 of file TopDQMHelpers.cc.

     : failed_(false),
       maxNJets_(maxNJets),
       wMass_(wMass),
       massWBoson_(-1.),
       massTopQuark_(-1.),
       massBTopQuark_(-1.),
       tmassWBoson_(-1),
       tmassTopQuark_(-1) {}

◆ ~Calculate()

Calculate::~Calculate ( )

inline

default destructor

Definition at line 132 of file TopDQMHelpers.h.

132 {};

Member Function Documentation

◆ massBTopQuark()

double Calculate::massBTopQuark	(	const std::vector< reco::Jet > &	jets,
		std::vector< double >	VbtagWP,
		double	btagWP_
	)

calculate b-tagged t-quark mass estimate

Definition at line 25 of file TopDQMHelpers.cc.

                                                                                                          {
   if (!failed_ && massBTopQuark_ < 0)
     operator2(jets, VbtagWP, btagWP_);
   return massBTopQuark_;
 }

References failed_, singleTopDQM_cfi::jets, massBTopQuark_, and operator2().

Referenced by TopSingleLepton::MonitorEnsemble::fill().

◆ massTopQuark()

double Calculate::massTopQuark ( const std::vector< reco::Jet > & jets )

calculate t-quark mass estimate

Definition at line 19 of file TopDQMHelpers.cc.

                                                              {
   if (!failed_ && massTopQuark_ < 0)
     operator()(jets);
   return massTopQuark_;
 }

References failed_, singleTopDQM_cfi::jets, massTopQuark_, and operator()().

Referenced by TopSingleLepton::MonitorEnsemble::fill(), and SingleTopTChannelLepton::MonitorEnsemble::fill().

◆ massWBoson()

double Calculate::massWBoson ( const std::vector< reco::Jet > & jets )

calculate W boson mass estimate

Definition at line 13 of file TopDQMHelpers.cc.

                                                            {
   if (!failed_ && massWBoson_ < 0)
     operator()(jets);
   return massWBoson_;
 }

References failed_, singleTopDQM_cfi::jets, massWBoson_, and operator()().

Referenced by TopSingleLepton::MonitorEnsemble::fill(), and SingleTopTChannelLepton::MonitorEnsemble::fill().

◆ operator()() [1/2]

void Calculate::operator()	(	const reco::Jet &	bJet,
		reco::RecoCandidate *	lepton,
		const reco::MET &	met
	)

private

do the calculation of the transverse top and W masses

Definition at line 43 of file TopDQMHelpers.cc.

                                                                                              {
   double metT = sqrt(pow(met.px(), 2) + pow(met.py(), 2));
   double lepT = sqrt(pow(lepton->px(), 2) + pow(lepton->py(), 2));
   double bT = sqrt(pow(bJet.px(), 2) + pow(bJet.py(), 2));
   reco::Particle::LorentzVector WT = lepton->p4() + met.p4();
   tmassWBoson_ = sqrt(pow(metT + lepT, 2) - (WT.px() * WT.px()) - (WT.py() * WT.py()));
   reco::Particle::LorentzVector topT = WT + bJet.p4();
   tmassTopQuark_ = sqrt(pow((metT + lepT + bT), 2) - (topT.px() * topT.px()) - (topT.py() * topT.py()));
 }

References BTaggingMonitor_cfi::met, reco::LeafCandidate::p4(), funct::pow(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), tmassTopQuark_, and tmassWBoson_.

◆ operator()() [2/2]

void Calculate::operator() ( const std::vector< reco::Jet > & jets )

private

do the calculation; this is called only once per event by the first function call to return a mass estimate. The once calculated values are cached afterwards

Definition at line 53 of file TopDQMHelpers.cc.

                                                          {
   if (maxNJets_ < 0)
     maxNJets_ = jets.size();
   failed_ = jets.size() < (unsigned int)maxNJets_;
   if (failed_) {
     return;
   }
  
   // associate those jets with maximum pt of the vectorial
   // sum to the hadronic decay chain
   double maxPt = -1.;
   std::vector<int> maxPtIndices;
   maxPtIndices.push_back(-1);
   maxPtIndices.push_back(-1);
   maxPtIndices.push_back(-1);
  
   for (int idx = 0; idx < maxNJets_; ++idx) {
     for (int jdx = idx + 1; jdx < maxNJets_; ++jdx) {
       //if (jdx <= idx) continue;
       for (int kdx = 0; kdx < maxNJets_; ++kdx) {
         if (kdx == idx || kdx == jdx)
           continue;
         reco::Particle::LorentzVector sum = jets[idx].p4() + jets[jdx].p4() + jets[kdx].p4();
         if (maxPt < 0. || maxPt < sum.pt()) {
           maxPt = sum.pt();
           maxPtIndices.clear();
           maxPtIndices.push_back(idx);
           maxPtIndices.push_back(jdx);
           maxPtIndices.push_back(kdx);
         }
       }
     }
   }
   massTopQuark_ = (jets[maxPtIndices[0]].p4() + jets[maxPtIndices[1]].p4() + jets[maxPtIndices[2]].p4()).mass();
  
   // associate those jets that get closest to the W mass
   // with their invariant mass to the W boson
   double wDist = -1.;
   std::vector<int> wMassIndices;
   wMassIndices.push_back(-1);
   wMassIndices.push_back(-1);
   for (unsigned idx = 0; idx < maxPtIndices.size(); ++idx) {
     for (unsigned jdx = 0; jdx < maxPtIndices.size(); ++jdx) {
       if (jdx == idx || maxPtIndices[idx] > maxPtIndices[jdx])
         continue;
       reco::Particle::LorentzVector sum = jets[maxPtIndices[idx]].p4() + jets[maxPtIndices[jdx]].p4();
       if (wDist < 0. || wDist > fabs(sum.mass() - wMass_)) {
         wDist = fabs(sum.mass() - wMass_);
         wMassIndices.clear();
         wMassIndices.push_back(maxPtIndices[idx]);
         wMassIndices.push_back(maxPtIndices[jdx]);
       }
     }
   }
   massWBoson_ = (jets[wMassIndices[0]].p4() + jets[wMassIndices[1]].p4()).mass();
 }

References failed_, heavyIonCSV_trainingSettings::idx, createfilelist::int, singleTopDQM_cfi::jets, EgHLTOffHistBins_cfi::mass, massTopQuark_, massWBoson_, maxNJets_, MuonErrorMatrixAnalyzer_cfi::maxPt, and wMass_.

Referenced by massTopQuark(), massWBoson(), tmassTopQuark(), and tmassWBoson().

◆ operator2()

void Calculate::operator2	(	const std::vector< reco::Jet > &	jets,
		std::vector< double >	bjet,
		double	btagWP
	)

private

do the calculation of the t-quark mass with one b-jet

Definition at line 110 of file TopDQMHelpers.cc.

                                                                                                {
   if (maxNJets_ < 0)
     maxNJets_ = jets.size();
   failed_ = jets.size() < (unsigned int)maxNJets_;
   if (failed_) {
     return;
   }
   if (jets.size() != bjet.size()) {
     return;
   }
  
   // associate those jets with maximum pt of the vectorial
   // sum to the hadronic decay chain. Require ONLY 1 btagged jet
   double maxBPt = -1.;
   std::vector<int> maxBPtIndices;
   maxBPtIndices.push_back(-1);
   maxBPtIndices.push_back(-1);
   maxBPtIndices.push_back(-1);
   for (int idx = 0; idx < maxNJets_; ++idx) {
     for (int jdx = idx + 1; jdx < maxNJets_; ++jdx) {
       //if (jdx <= idx) continue;
       for (int kdx = 0; kdx < maxNJets_; ++kdx) {
         if (kdx == idx || kdx == jdx)
           continue;
         // require only 1b-jet
         if ((bjet[idx] > btagWP && bjet[jdx] <= btagWP && bjet[kdx] <= btagWP) ||
             (bjet[idx] <= btagWP && bjet[jdx] > btagWP && bjet[kdx] <= btagWP) ||
             (bjet[idx] <= btagWP && bjet[jdx] <= btagWP && bjet[kdx] > btagWP)) {
           reco::Particle::LorentzVector sum = jets[idx].p4() + jets[jdx].p4() + jets[kdx].p4();
           if (maxBPt < 0. || maxBPt < sum.pt()) {
             maxBPt = sum.pt();
             maxBPtIndices.clear();
             maxBPtIndices.push_back(idx);
             maxBPtIndices.push_back(jdx);
             maxBPtIndices.push_back(kdx);
           }
         }
       }
     }
   }
   if (maxBPtIndices[0] < 0 || maxBPtIndices[1] < 0 || maxBPtIndices[2] < 0)
     return;
   massBTopQuark_ = (jets[maxBPtIndices[0]].p4() + jets[maxBPtIndices[1]].p4() + jets[maxBPtIndices[2]].p4()).mass();
 }

References failed_, heavyIonCSV_trainingSettings::idx, createfilelist::int, singleTopDQM_cfi::jets, EgHLTOffHistBins_cfi::mass, massBTopQuark_, and maxNJets_.

Referenced by massBTopQuark().

◆ tmassTopQuark()

double Calculate::tmassTopQuark	(	reco::RecoCandidate *	lep,
		const reco::MET &	met,
		const reco::Jet &	b
	)

calculate top quark transverse mass estimate

Definition at line 37 of file TopDQMHelpers.cc.

                                                                                                {
   if (tmassTopQuark_ < 0)
     operator()(b, lepton, met);
   return tmassTopQuark_;
 }

References b, BTaggingMonitor_cfi::met, operator()(), and tmassTopQuark_.

Referenced by SingleTopTChannelLepton::MonitorEnsemble::fill().

◆ tmassWBoson()

double Calculate::tmassWBoson	(	reco::RecoCandidate *	lep,
		const reco::MET &	met,
		const reco::Jet &	b
	)

calculate W boson transverse mass estimate

Definition at line 31 of file TopDQMHelpers.cc.

                                                                                          {
   if (tmassWBoson_ < 0)
     operator()(b, mu, met);
   return tmassWBoson_;
 }

References b, BTaggingMonitor_cfi::met, amptDefaultParameters_cff::mu, operator()(), and tmassWBoson_.

Referenced by SingleTopTChannelLepton::MonitorEnsemble::fill().

Member Data Documentation

◆ failed_

bool Calculate::failed_

private

indicate failed associations

Definition at line 160 of file TopDQMHelpers.h.

Referenced by massBTopQuark(), massTopQuark(), massWBoson(), operator()(), and operator2().

◆ massBTopQuark_

double Calculate::massBTopQuark_

private

cache of b-tagged top quark mass estimate

Definition at line 170 of file TopDQMHelpers.h.

Referenced by massBTopQuark(), and operator2().

◆ massTopQuark_

double Calculate::massTopQuark_

private

cache of top quark mass estimate

Definition at line 168 of file TopDQMHelpers.h.

Referenced by massTopQuark(), and operator()().

◆ massWBoson_

double Calculate::massWBoson_

private

cache of w boson mass estimate

Definition at line 166 of file TopDQMHelpers.h.

Referenced by massWBoson(), and operator()().

◆ maxNJets_

int Calculate::maxNJets_

private

max. number of jets to be considered

Definition at line 162 of file TopDQMHelpers.h.

Referenced by operator()(), and operator2().

◆ tmassTopQuark_

double Calculate::tmassTopQuark_

private

cache of top quark transverse mass estimate

Definition at line 174 of file TopDQMHelpers.h.

Referenced by operator()(), and tmassTopQuark().

◆ tmassWBoson_

double Calculate::tmassWBoson_

private

cache of W boson transverse mass estimate

Definition at line 172 of file TopDQMHelpers.h.

Referenced by operator()(), and tmassWBoson().

◆ wMass_

double Calculate::wMass_

private

paramater of the w boson mass

Definition at line 164 of file TopDQMHelpers.h.

Referenced by operator()().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Calculate()

◆ ~Calculate()

Member Function Documentation

◆ massBTopQuark()

◆ massTopQuark()

◆ massWBoson()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ operator2()

◆ tmassTopQuark()

◆ tmassWBoson()

Member Data Documentation

◆ failed_

◆ massBTopQuark_

◆ massTopQuark_

◆ massWBoson_

◆ maxNJets_

◆ tmassTopQuark_

◆ tmassWBoson_

◆ wMass_