#include <TopHLTOfflineDQMHelper.h>

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

double	masslb (reco::RecoCandidate *mu, const reco::MET &met, const reco::Jet &b)
	calculate mlb estimate More...

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

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

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

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

	~CalculateHLT ()
	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)

Private Attributes
bool	failed_
	indicate failed associations 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	mlb_
	cache of mlb estimate 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

Definition at line 72 of file TopHLTOfflineDQMHelper.h.

Constructor & Destructor Documentation

CalculateHLT::CalculateHLT	(	int	maxNJets,
		double	wMass
	)

default constructor

Definition at line 81 of file TopHLTOfflineDQMHelper.cc.

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

CalculateHLT::~CalculateHLT ( )

inline

default destructor

Definition at line 77 of file TopHLTOfflineDQMHelper.h.

References b, fwrapper::jets, masslb(), massTopQuark(), massWBoson(), RazorAnalyzer::met, RPCpg::mu, operator()(), tmassTopQuark(), and tmassWBoson().

77 {};

Member Function Documentation

double CalculateHLT::masslb	(	reco::RecoCandidate *	mu,
		const reco::MET &	met,
		const reco::Jet &	b
	)

calculate mlb estimate

Definition at line 143 of file TopHLTOfflineDQMHelper.cc.

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

Referenced by HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), and ~CalculateHLT().

 {
   if(!failed_&& mlb_<0) operator()(b,mu,met);
   return mlb_;
 }

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

calculate top quark mass estimate

Definition at line 96 of file TopHLTOfflineDQMHelper.cc.

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

Referenced by HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), and ~CalculateHLT().

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

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

calculate W boson mass estimate

Definition at line 88 of file TopHLTOfflineDQMHelper.cc.

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

Referenced by HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), and ~CalculateHLT().

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

void CalculateHLT::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 173 of file TopHLTOfflineDQMHelper.cc.

References failed_, training_settings::idx, createfilelist::int, ResonanceBuilder::mass, massTopQuark_, massWBoson_, maxNJets_, MuonErrorMatrixAnalyzer_cfi::maxPt, and wMass_.

Referenced by masslb(), massTopQuark(), massWBoson(), tmassTopQuark(), tmassWBoson(), and ~CalculateHLT().

 {
   
   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 CalculateHLT::operator()	(	const reco::Jet &	bJet,
		reco::RecoCandidate *	lepton,
		const reco::MET &	met
	)

private

Definition at line 158 of file TopHLTOfflineDQMHelper.cc.

References mlb_, 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();
   //  cout<<"in calculate:\n\t"<<bJet.pt()<<"\t"<<lepton->pt()<<"\t"<<met.pt()<<endl;
   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()));
   reco::Particle::LorentzVector lb = bJet.p4() + lepton->p4();  
   mlb_ = lb.mass();
 }

double CalculateHLT::tmassTopQuark	(	reco::RecoCandidate *	mu,
		const reco::MET &	met,
		const reco::Jet &	b
	)

calculate top quark transverse mass estimate

Definition at line 151 of file TopHLTOfflineDQMHelper.cc.

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

Referenced by HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), and ~CalculateHLT().

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

double CalculateHLT::tmassWBoson	(	reco::RecoCandidate *	mu,
		const reco::MET &	met,
		const reco::Jet &	b
	)

calculate W boson transverse mass estimate

Definition at line 135 of file TopHLTOfflineDQMHelper.cc.

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

Referenced by HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), and ~CalculateHLT().

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

Member Data Documentation

bool CalculateHLT::failed_

private

indicate failed associations

Definition at line 105 of file TopHLTOfflineDQMHelper.h.

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

double CalculateHLT::massTopQuark_

private

cache of top quark mass estimate

Definition at line 113 of file TopHLTOfflineDQMHelper.h.

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

double CalculateHLT::massWBoson_

private

cache of w boson mass estimate

Definition at line 111 of file TopHLTOfflineDQMHelper.h.

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

int CalculateHLT::maxNJets_

private

max. number of jets to be considered

Definition at line 107 of file TopHLTOfflineDQMHelper.h.

Referenced by operator()().

double CalculateHLT::mlb_

private

cache of mlb estimate

Definition at line 119 of file TopHLTOfflineDQMHelper.h.

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

double CalculateHLT::tmassTopQuark_

private

cache of top quark transverse mass estimate

Definition at line 117 of file TopHLTOfflineDQMHelper.h.

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

double CalculateHLT::tmassWBoson_

private

cache of W boson transverse mass estimate

Definition at line 115 of file TopHLTOfflineDQMHelper.h.

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

double CalculateHLT::wMass_

private

paramater of the w boson mass

Definition at line 109 of file TopHLTOfflineDQMHelper.h.

Referenced by operator()().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation