Calculate Class Reference

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 std::vector< reco::Jet > &jets)
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	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 69 of file TopDQMHelpers.h.

Constructor & Destructor Documentation

Calculate::Calculate	(	int	maxNJets,
		double	wMass
	)

default constructor

Definition at line 4 of file TopDQMHelpers.cc.

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

Calculate::~Calculate ( )

inline

default destructor

Definition at line 74 of file TopDQMHelpers.h.

{};

Member Function Documentation

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 23 of file TopDQMHelpers.cc.

References failed_, massBTopQuark_, and operator2().

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

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

double Calculate::massTopQuark

(

const std::vector< reco::Jet > &

jets

)

calculate t-quark mass estimate

Definition at line 16 of file TopDQMHelpers.cc.

References failed_, massTopQuark_, and operator()().

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

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

double Calculate::massWBoson

(

const std::vector< reco::Jet > &

jets

)

calculate W boson mass estimate

Definition at line 10 of file TopDQMHelpers.cc.

References failed_, massWBoson_, and operator()().

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

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

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 55 of file TopDQMHelpers.cc.

References failed_, customizeTrackingMonitorSeedNumber::idx, massTopQuark_, massWBoson_, maxNJets_, and wMass_.

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

{
  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=0; 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();
}

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 42 of file TopDQMHelpers.cc.

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

                                                                                             {
  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()));
}

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.

References failed_, customizeTrackingMonitorSeedNumber::idx, massBTopQuark_, and maxNJets_.

Referenced by massBTopQuark().

{
  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=0; 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();
}

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

calculate top quark transverse mass estimate

Definition at line 36 of file TopDQMHelpers.cc.

References operator()(), and tmassTopQuark_.

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

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

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

calculate W boson transverse mass estimate

Definition at line 29 of file TopDQMHelpers.cc.

References operator()(), and tmassWBoson_.

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

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

Member Data Documentation

bool Calculate::failed_

private

indicate failed associations

Definition at line 102 of file TopDQMHelpers.h.

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

double Calculate::massBTopQuark_

private

cache of b-tagged top quark mass estimate

Definition at line 112 of file TopDQMHelpers.h.

Referenced by massBTopQuark(), and operator2().

double Calculate::massTopQuark_

private

cache of top quark mass estimate

Definition at line 110 of file TopDQMHelpers.h.

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

double Calculate::massWBoson_

private

cache of w boson mass estimate

Definition at line 108 of file TopDQMHelpers.h.

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

int Calculate::maxNJets_

private

max. number of jets to be considered

Definition at line 104 of file TopDQMHelpers.h.

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

double Calculate::tmassTopQuark_

private

cache of top quark transverse mass estimate

Definition at line 116 of file TopDQMHelpers.h.

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

double Calculate::tmassWBoson_

private

cache of W boson transverse mass estimate

Definition at line 114 of file TopDQMHelpers.h.

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

double Calculate::wMass_

private

paramater of the w boson mass

Definition at line 106 of file TopDQMHelpers.h.

Referenced by operator()().