TtSemiLepJetComb Class Reference

Common calculator class to keep multivariate analysis variables for jet combinations in semi-leptonic ttbar decays. More...

#include "TopQuarkAnalysis/TopJetCombination/interface/TtSemiLepJetComb.h"

Public Member Functions
double	addUserVar (std::string key, double value)
	add an arbitary user defined variable with given key and value More...
double	bOverLightQPt () const
double	bQuarkVar (JetComb::DecayType decay, JetComb::VarType var) const
	b quark candidate variable More...
double	bTag (JetComb::DecayType decay, JetComb::BTagAlgo algo) const
	b-tag value of a b candidate More...
double	combinedBTags (JetComb::BTagAlgo algo, JetComb::Operator op) const
	combined b-tag values of the two b candidates More...
double	combinedBTagsForLightQuarks (JetComb::BTagAlgo algo, JetComb::Operator op) const
	combined b-tag values of the two light quark candidates More...
double	compareBLepton (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the b and the lepton candidate More...
double	compareBNeutrino (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the b and the neutrino candidate More...
double	compareHadBLepB (JetComb::CompType comp) const
	comparison between the two b candidates More...
double	compareHadTopLepTop (JetComb::CompType comp) const
	comparison between the two top candidates More...
double	compareHadWLepW (JetComb::CompType comp) const
	comparison between the two W candidates More...
double	compareLeptonNeutrino (JetComb::CompType comp) const
	comparison between the lepton and the neutrino candidate More...
double	compareLightQuarks (JetComb::CompType comp) const
	comparison between the lightQ and the lightQBar candidate More...
double	compareTopB (JetComb::DecayType dec1, JetComb::DecayType dec2, JetComb::CompType comp) const
	comparison between the top and the b candidate More...
double	compareTopLepton (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the top and the lepton candidate More...
double	compareTopNeutrino (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the top and the neutrino candidate More...
double	compareTopW (JetComb::DecayType dec1, JetComb::DecayType dec2, JetComb::CompType comp) const
	comparison between the top and the W candidate More...
double	compareWB (JetComb::DecayType dec1, JetComb::DecayType dec2, JetComb::CompType comp) const
	comparison between the W and the b candidate More...
double	compareWLepton (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the W and the lepton candidate More...
double	compareWNeutrino (JetComb::DecayType decay, JetComb::CompType comp) const
	comparison between the W and the neutrino candidate More...
double	leptonVar (JetComb::VarType var) const
	lepton candidate variable More...
double	lightQVar (bool qbar, JetComb::VarType var) const
	light quark candidate variable More...
double	neutrinoVar (JetComb::VarType var) const
	neutrino candidate variable More...
double	relativePtHadronicTop () const
double	topVar (JetComb::DecayType decay, JetComb::VarType var) const
	top quark candidate variable More...
	TtSemiLepJetComb ()
	emtpy constructor More...
	TtSemiLepJetComb (const std::vector< pat::Jet > &, const std::vector< int > &, const math::XYZTLorentzVector &, const pat::MET &)
	default constructor More...
double	userVar (const std::string &key) const
	receive user defined variable value with given key More...
double	wBosonVar (JetComb::DecayType decay, JetComb::VarType var) const
	W boson candidate variable. More...
	~TtSemiLepJetComb ()
	default destructor More...

Private Member Functions
const pat::Jet &	bQuark (JetComb::DecayType decay) const
	return leptonic or hadronic b candidate depending on argument More...
double	bTag (const pat::Jet &jet, JetComb::BTagAlgo algo) const
	b-tag value of one of the 4 jets More...
void	deduceMothers ()
	reconstruct mother candidates from final state candidates More...
const pat::Jet &	lightQ (bool qbar=false) const
	return lightQ or lightQBar candidate depending on argument More...
double	lightQVar (JetComb::VarType var) const
	light quark candidate variable with default on q and not qbar More...
const math::XYZTLorentzVector &	top (JetComb::DecayType decay) const
	return leptonic or hadronic top candidate depending on argument More...
const math::XYZTLorentzVector &	wBoson (JetComb::DecayType decay) const
	return leptonic or hadronic W candidate depending on argument More...

Private Attributes
pat::Jet	hadBJet_
	hadronic b jet More...
pat::Jet	hadQBarJet_
	lightQBar jet More...
pat::Jet	hadQJet_
	lightQ jet More...
math::XYZTLorentzVector	hadTop_
	hadronic top 4-vector More...
math::XYZTLorentzVector	hadW_
	hadronic W 4-vector More...
pat::Jet	lepBJet_
	leptonic b jet More...
math::XYZTLorentzVector	lepton_
	lepton 4-vector More...
math::XYZTLorentzVector	lepTop_
	leptonic top 4-vector More...
math::XYZTLorentzVector	lepW_
	leptonic W 4-vector More...
pat::MET	neutrino_
	neutrino candidate More...
std::map< std::string, double >	userVariables_
	map for user defined variables More...

Detailed Description

Common calculator class to keep multivariate analysis variables for jet combinations in semi-leptonic ttbar decays.

This class is currently used by TtSemiLepJetCombEval.h, which is included in both the TtSemiLepJetCombMVAComputer and the TtSemiLepJetCombMVATrainer.

Definition at line 37 of file TtSemiLepJetComb.h.

Constructor & Destructor Documentation

TtSemiLepJetComb::TtSemiLepJetComb

(

)

emtpy constructor

Definition at line 6 of file TtSemiLepJetComb.cc.

{
}

TtSemiLepJetComb::TtSemiLepJetComb	(	const std::vector< pat::Jet > &	jets,
		const std::vector< int > &	combination,
		const math::XYZTLorentzVector &	lepton,
		const pat::MET &	neutrino
	)

default constructor

Definition at line 10 of file TtSemiLepJetComb.cc.

References deduceMothers(), TtSemiLepEvtPartons::HadB, hadBJet_, hadQBarJet_, hadQJet_, TtSemiLepEvtPartons::LepB, lepBJet_, lepton_, TtSemiLepEvtPartons::LightQ, TtSemiLepEvtPartons::LightQBar, and neutrino_.

{ 
  // receive right jet association
  // from jet-parton matching 
  hadQJet_    = jets[ combination[TtSemiLepEvtPartons::LightQ   ] ];
  hadQBarJet_ = jets[ combination[TtSemiLepEvtPartons::LightQBar] ];
  hadBJet_    = jets[ combination[TtSemiLepEvtPartons::HadB     ] ];
  lepBJet_    = jets[ combination[TtSemiLepEvtPartons::LepB     ] ]; 
  lepton_     = lepton;
  neutrino_   = neutrino;
  // create mother candidates from 
  // final-state candidates
  deduceMothers();
}

TtSemiLepJetComb::~TtSemiLepJetComb

(

)

default destructor

Definition at line 26 of file TtSemiLepJetComb.cc.

{
}

Member Function Documentation

double TtSemiLepJetComb::addUserVar ( std::string  key, double  value  )

inline

add an arbitary user defined variable with given key and value

Definition at line 106 of file TtSemiLepJetComb.h.

References crabWrapper::key, and relativeConstraints::value.

{ return userVariables_[key]=value;};

double TtSemiLepJetComb::bOverLightQPt

(

)

const

scalar sum of the pt of the b candidates relative to scalar sum of the pt of the light quark candidates

Definition at line 265 of file TtSemiLepJetComb.cc.

References bQuark(), JetComb::kHad, JetComb::kLep, lightQ(), p4, reco::LeafCandidate::p4(), and EnergyCorrector::pt.

Referenced by evaluateTtSemiLepJetComb().

{
  return (bQuark(JetComb::kHad).p4().pt()+bQuark(JetComb::kLep).p4().pt())/(lightQ().p4().pt()+lightQ(true).p4().pt());  
}

const pat::Jet& TtSemiLepJetComb::bQuark ( JetComb::DecayType  decay ) const

inlineprivate

return leptonic or hadronic b candidate depending on argument

Definition at line 122 of file TtSemiLepJetComb.h.

References JetComb::kHad.

Referenced by bOverLightQPt(), bQuarkVar(), compareBLepton(), compareBNeutrino(), compareHadBLepB(), compareTopB(), compareWB(), and relativePtHadronicTop().

{ return (decay==JetComb::kHad ? hadBJet_ : lepBJet_); }

double TtSemiLepJetComb::bQuarkVar	(	JetComb::DecayType	decay,
		JetComb::VarType	var
	)		const

b quark candidate variable

Definition at line 54 of file TtSemiLepJetComb.cc.

References bQuark(), JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, and reco::LeafCandidate::p4().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(var){
  case JetComb::kMass  : return -9999.;
  case JetComb::kPt    : return bQuark(decay).p4().pt();
  case JetComb::kEta   : return bQuark(decay).p4().eta();
  case JetComb::kPhi   : return bQuark(decay).p4().phi();
  case JetComb::kTheta : return bQuark(decay).p4().theta();
  };
  return -9999.;
}

double TtSemiLepJetComb::bTag ( JetComb::DecayType  decay, JetComb::BTagAlgo  algo  ) const

inline

b-tag value of a b candidate

Definition at line 99 of file TtSemiLepJetComb.h.

References patPFMETCorrections_cff::algo, and bTag().

Referenced by bTag(), combinedBTags(), combinedBTagsForLightQuarks(), and evaluateTtSemiLepJetComb().

{ return bTag(bQuark(decay), algo); };

double TtSemiLepJetComb::bTag ( const pat::Jet &  jet, JetComb::BTagAlgo  algo  ) const

private

b-tag value of one of the 4 jets

Definition at line 301 of file TtSemiLepJetComb.cc.

References pat::Jet::bDiscriminator(), JetComb::kBProbability, JetComb::kCombSecondVtx, JetComb::kCombSecondVtxMVA, JetComb::kProbability, JetComb::kSimpleSecondVtx, JetComb::kSofMuonByIP3d, JetComb::kSoftElec, JetComb::kSoftMuon, JetComb::kSoftMuonByPt, JetComb::kTrackCountHighEff, and JetComb::kTrackCountHighPur.

{
  switch(algo){
  case JetComb::kTrackCountHighEff : return jet.bDiscriminator("trackCountingHighEffBJetTags"      );
  case JetComb::kTrackCountHighPur : return jet.bDiscriminator("trackCountingHighPurBJetTags"      );
  case JetComb::kSoftMuon          : return jet.bDiscriminator("softMuonBJetTags"                  );
  case JetComb::kSoftMuonByPt      : return jet.bDiscriminator("softMuonByPtBJetTags"              );
  case JetComb::kSofMuonByIP3d     : return jet.bDiscriminator("softMuonByIP3dBJetTags"            );
  case JetComb::kSoftElec          : return jet.bDiscriminator("softElectronBJetTags"              );
  case JetComb::kProbability       : return jet.bDiscriminator("jetProbabilityBJetTags"            );
  case JetComb::kBProbability      : return jet.bDiscriminator("jetBProbabilityBJetTags"           );
  case JetComb::kSimpleSecondVtx   : return jet.bDiscriminator("simpleSecondaryVertexBJetTags"     );
  case JetComb::kCombSecondVtx     : return jet.bDiscriminator("combinedSecondaryVertexBJetTags"   );
  case JetComb::kCombSecondVtxMVA  : return jet.bDiscriminator("combinedSecondaryVertexMVABJetTags");
  };
  return -9999.;
}

double TtSemiLepJetComb::combinedBTags	(	JetComb::BTagAlgo	algo,
		JetComb::Operator	op
	)		const

combined b-tag values of the two b candidates

Definition at line 270 of file TtSemiLepJetComb.cc.

References bTag(), JetComb::kAdd, JetComb::kHad, JetComb::kLep, and JetComb::kMult.

Referenced by evaluateTtSemiLepJetComb().

{
  switch(op){
  case JetComb::kAdd  : return bTag(JetComb::kHad, algo) + bTag(JetComb::kLep, algo);
  case JetComb::kMult : return bTag(JetComb::kHad, algo) * bTag(JetComb::kLep, algo);
  };
  return -9999.;
}

double TtSemiLepJetComb::combinedBTagsForLightQuarks	(	JetComb::BTagAlgo	algo,
		JetComb::Operator	op
	)		const

combined b-tag values of the two light quark candidates

Definition at line 279 of file TtSemiLepJetComb.cc.

References bTag(), JetComb::kAdd, JetComb::kMult, and lightQ().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(op){
  case JetComb::kAdd  : return bTag(lightQ(), algo) + bTag(lightQ(true), algo);
  case JetComb::kMult : return bTag(lightQ(), algo) * bTag(lightQ(true), algo);
  };
  return -9999.;
}

double TtSemiLepJetComb::compareBLepton	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the b and the lepton candidate

Definition at line 234 of file TtSemiLepJetComb.cc.

References bQuark(), allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, lepton_, and p4.

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (bQuark(decay).p4(), lepton_);
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(bQuark(decay).p4(), lepton_));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (bQuark(decay).p4(), lepton_);
  };
  return -9999.;
}

double TtSemiLepJetComb::compareBNeutrino	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the b and the neutrino candidate

Definition at line 245 of file TtSemiLepJetComb.cc.

References bQuark(), allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, neutrino_, p4, and reco::LeafCandidate::p4().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (bQuark(decay).p4(), neutrino_.p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(bQuark(decay).p4(), neutrino_.p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (bQuark(decay).p4(), neutrino_.p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareHadBLepB

(

JetComb::CompType

comp

)

const

comparison between the two b candidates

Definition at line 124 of file TtSemiLepJetComb.cc.

References bQuark(), allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, JetComb::kHad, JetComb::kLep, and p4.

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (bQuark(JetComb::kHad).p4(), bQuark(JetComb::kLep).p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(bQuark(JetComb::kHad).p4(), bQuark(JetComb::kLep).p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (bQuark(JetComb::kHad).p4(), bQuark(JetComb::kLep).p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareHadTopLepTop

(

JetComb::CompType

comp

)

const

comparison between the two top candidates

Definition at line 102 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, JetComb::kHad, JetComb::kLep, and top().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return top(JetComb::kHad).mass() - top(JetComb::kLep).mass();  
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (top(JetComb::kHad), top(JetComb::kLep));
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(top(JetComb::kHad), top(JetComb::kLep)));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (top(JetComb::kHad), top(JetComb::kLep));
  };
  return -9999.;
}

double TtSemiLepJetComb::compareHadWLepW

(

JetComb::CompType

comp

)

const

comparison between the two W candidates

Definition at line 113 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, JetComb::kHad, JetComb::kLep, and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return wBoson(JetComb::kHad).mass() - wBoson(JetComb::kLep).mass();
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (wBoson(JetComb::kHad), wBoson(JetComb::kLep));
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(wBoson(JetComb::kHad), wBoson(JetComb::kLep)));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (wBoson(JetComb::kHad), wBoson(JetComb::kLep));
  };
  return -9999.;
}

double TtSemiLepJetComb::compareLeptonNeutrino

(

JetComb::CompType

comp

)

const

comparison between the lepton and the neutrino candidate

Definition at line 146 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, lepton_, neutrino_, and reco::LeafCandidate::p4().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (lepton_, neutrino_.p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(lepton_, neutrino_.p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (lepton_, neutrino_.p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareLightQuarks

(

JetComb::CompType

comp

)

const

comparison between the lightQ and the lightQBar candidate

Definition at line 135 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, lightQ(), and p4.

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (lightQ().p4(), lightQ(true).p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(lightQ().p4(), lightQ(true).p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (lightQ().p4(), lightQ(true).p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareTopB	(	JetComb::DecayType	dec1,
		JetComb::DecayType	dec2,
		JetComb::CompType	comp
	)		const

comparison between the top and the b candidate

Definition at line 168 of file TtSemiLepJetComb.cc.

References bQuark(), allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, p4, and top().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (top(dec1), bQuark(dec2).p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(top(dec1), bQuark(dec2).p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (top(dec1), bQuark(dec2).p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareTopLepton	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the top and the lepton candidate

Definition at line 190 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, lepton_, and top().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (top(decay), lepton_);
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(top(decay), lepton_));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (top(decay), lepton_);
  };
  return -9999.;
}

double TtSemiLepJetComb::compareTopNeutrino	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the top and the neutrino candidate

Definition at line 201 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, neutrino_, reco::LeafCandidate::p4(), and top().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (top(decay), neutrino_.p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(top(decay), neutrino_.p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (top(decay), neutrino_.p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareTopW	(	JetComb::DecayType	dec1,
		JetComb::DecayType	dec2,
		JetComb::CompType	comp
	)		const

comparison between the top and the W candidate

Definition at line 157 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, top(), and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return top(dec1).mass() - wBoson(dec2).mass();
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (top(dec1), wBoson(dec2));
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(top(dec1), wBoson(dec2)));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (top(dec1), wBoson(dec2));
  };
  return -9999.;
}

double TtSemiLepJetComb::compareWB	(	JetComb::DecayType	dec1,
		JetComb::DecayType	dec2,
		JetComb::CompType	comp
	)		const

comparison between the W and the b candidate

Definition at line 179 of file TtSemiLepJetComb.cc.

References bQuark(), allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, p4, and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (wBoson(dec1), bQuark(dec2).p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(wBoson(dec1), bQuark(dec2).p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (wBoson(dec1), bQuark(dec2).p4());
  };
  return -9999.;
}

double TtSemiLepJetComb::compareWLepton	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the W and the lepton candidate

Definition at line 212 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, lepton_, and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (wBoson(decay), lepton_);
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(wBoson(decay), lepton_));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (wBoson(decay), lepton_);
  };
  return -9999.;
}

double TtSemiLepJetComb::compareWNeutrino	(	JetComb::DecayType	decay,
		JetComb::CompType	comp
	)		const

comparison between the W and the neutrino candidate

Definition at line 223 of file TtSemiLepJetComb.cc.

References allConversions_cfi::DeltaPhi, HiRegitMuonDetachedTripletStep_cff::DeltaR, JetComb::kDeltaM, JetComb::kDeltaPhi, JetComb::kDeltaR, JetComb::kDeltaTheta, neutrino_, reco::LeafCandidate::p4(), and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(comp){
  case JetComb::kDeltaM     : return -9999.;
  case JetComb::kDeltaR     : return      ROOT::Math::VectorUtil::DeltaR  (wBoson(decay), neutrino_.p4());
  case JetComb::kDeltaPhi   : return fabs(ROOT::Math::VectorUtil::DeltaPhi(wBoson(decay), neutrino_.p4()));
  case JetComb::kDeltaTheta : return      ROOT::Math::VectorUtil::Angle   (wBoson(decay), neutrino_.p4());
  };
  return -9999.;
}

void TtSemiLepJetComb::deduceMothers ( )

private

reconstruct mother candidates from final state candidates

Definition at line 293 of file TtSemiLepJetComb.cc.

References hadBJet_, hadQBarJet_, hadQJet_, hadTop_, hadW_, lepBJet_, lepton_, lepTop_, lepW_, neutrino_, and reco::LeafCandidate::p4().

Referenced by TtSemiLepJetComb().

{
  hadW_   = hadQJet_.p4() + hadQBarJet_.p4();
  lepW_   = lepton_       + neutrino_  .p4();
  hadTop_ = hadW_         + hadBJet_   .p4();
  lepTop_ = lepW_         + lepBJet_   .p4();
}

double TtSemiLepJetComb::leptonVar

(

JetComb::VarType

var

)

const

lepton candidate variable

Definition at line 78 of file TtSemiLepJetComb.cc.

References JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, and lepton_.

{
  switch(var){
  case JetComb::kMass  : return -9999.;
  case JetComb::kPt    : return lepton_.pt();
  case JetComb::kEta   : return lepton_.eta();
  case JetComb::kPhi   : return lepton_.phi();
  case JetComb::kTheta : return lepton_.theta();
  };
  return -9999.;
}

const pat::Jet& TtSemiLepJetComb::lightQ ( bool  qbar = false ) const

inlineprivate

return lightQ or lightQBar candidate depending on argument

Definition at line 120 of file TtSemiLepJetComb.h.

Referenced by bOverLightQPt(), combinedBTagsForLightQuarks(), compareLightQuarks(), lightQVar(), and relativePtHadronicTop().

{ return (qbar ? hadQBarJet_ : hadQJet_); }

double TtSemiLepJetComb::lightQVar	(	bool	qbar,
		JetComb::VarType	var
	)		const

light quark candidate variable

Definition at line 66 of file TtSemiLepJetComb.cc.

References JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, lightQ(), and reco::LeafCandidate::p4().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(var){
  case JetComb::kMass  : return -9999.;
  case JetComb::kPt    : return lightQ(qbar).p4().pt();
  case JetComb::kEta   : return lightQ(qbar).p4().eta();
  case JetComb::kPhi   : return lightQ(qbar).p4().phi();
  case JetComb::kTheta : return lightQ(qbar).p4().theta();
  };
  return -9999.;
}

double TtSemiLepJetComb::lightQVar ( JetComb::VarType  var ) const

inlineprivate

light quark candidate variable with default on q and not qbar

Definition at line 118 of file TtSemiLepJetComb.h.

References lightQVar().

Referenced by lightQVar().

{ return lightQVar(false, var); };

double TtSemiLepJetComb::neutrinoVar

(

JetComb::VarType

var

)

const

neutrino candidate variable

Definition at line 90 of file TtSemiLepJetComb.cc.

References JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, neutrino_, and reco::LeafCandidate::p4().

{
  switch(var){
  case JetComb::kMass  : return -9999.;
  case JetComb::kPt    : return neutrino_.p4().pt();
  case JetComb::kEta   : return neutrino_.p4().eta();
  case JetComb::kPhi   : return neutrino_.p4().phi();
  case JetComb::kTheta : return neutrino_.p4().theta();
  };
  return -9999.;
}

double TtSemiLepJetComb::relativePtHadronicTop

(

)

const

pt of the hadronic top candidate relative to pt of the sum of all other reconstruction possibilities (returns values between 0 and 1)

Definition at line 256 of file TtSemiLepJetComb.cc.

References bQuark(), JetComb::kHad, JetComb::kLep, lightQ(), reco::LeafCandidate::p4(), EnergyCorrector::pt, and top().

Referenced by evaluateTtSemiLepJetComb().

{
  return top(JetComb::kHad).pt()/(top(JetComb::kHad).pt()                                               +
             (lightQ()    .p4()+lightQ(true).p4()+bQuark(JetComb::kLep).p4()).pt()          +
             (lightQ()    .p4()+bQuark(JetComb::kHad).p4()+bQuark(JetComb::kLep).p4()).pt() +
             (bQuark(JetComb::kHad).p4()+lightQ(true).p4()+bQuark(JetComb::kLep).p4()).pt()
             ); 
}

const math::XYZTLorentzVector& TtSemiLepJetComb::top ( JetComb::DecayType  decay ) const

inlineprivate

return leptonic or hadronic top candidate depending on argument

Definition at line 126 of file TtSemiLepJetComb.h.

References JetComb::kHad.

Referenced by compareHadTopLepTop(), compareTopB(), compareTopLepton(), compareTopNeutrino(), compareTopW(), relativePtHadronicTop(), and topVar().

{ return (decay==JetComb::kHad ? hadTop_ : lepTop_); }  

double TtSemiLepJetComb::topVar	(	JetComb::DecayType	decay,
		JetComb::VarType	var
	)		const

top quark candidate variable

Definition at line 30 of file TtSemiLepJetComb.cc.

References JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, and top().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(var){
  case JetComb::kMass  : return top(decay).mass();
  case JetComb::kPt    : return top(decay).pt();
  case JetComb::kEta   : return top(decay).eta();
  case JetComb::kPhi   : return top(decay).phi();
  case JetComb::kTheta : return top(decay).theta();
  };
  return -9999.;
}

double TtSemiLepJetComb::userVar ( const std::string &  key ) const

inline

receive user defined variable value with given key

Definition at line 108 of file TtSemiLepJetComb.h.

References TtFullHadHypKinFit_cfi::bTag, and metsig::jet.

                                             {
    return (userVariables_.find(key)!=userVariables_.end() ? userVariables_.find(key)->second : -9999.);};

const math::XYZTLorentzVector& TtSemiLepJetComb::wBoson ( JetComb::DecayType  decay ) const

inlineprivate

return leptonic or hadronic W candidate depending on argument

Definition at line 124 of file TtSemiLepJetComb.h.

References JetComb::kHad.

Referenced by compareHadWLepW(), compareTopW(), compareWB(), compareWLepton(), compareWNeutrino(), and wBosonVar().

{ return (decay==JetComb::kHad ? hadW_ : lepW_); }

double TtSemiLepJetComb::wBosonVar	(	JetComb::DecayType	decay,
		JetComb::VarType	var
	)		const

W boson candidate variable.

Definition at line 42 of file TtSemiLepJetComb.cc.

References JetComb::kEta, JetComb::kMass, JetComb::kPhi, JetComb::kPt, JetComb::kTheta, and wBoson().

Referenced by evaluateTtSemiLepJetComb().

{
  switch(var){
  case JetComb::kMass  : return wBoson(decay).mass();
  case JetComb::kPt    : return wBoson(decay).pt();
  case JetComb::kEta   : return wBoson(decay).eta();
  case JetComb::kPhi   : return wBoson(decay).phi();
  case JetComb::kTheta : return wBoson(decay).theta();
  };
  return -9999.;
}

Member Data Documentation

pat::Jet TtSemiLepJetComb::hadBJet_

private

hadronic b jet

Definition at line 135 of file TtSemiLepJetComb.h.

Referenced by deduceMothers(), and TtSemiLepJetComb().

pat::Jet TtSemiLepJetComb::hadQBarJet_

private

lightQBar jet

Definition at line 133 of file TtSemiLepJetComb.h.

Referenced by deduceMothers(), and TtSemiLepJetComb().

pat::Jet TtSemiLepJetComb::hadQJet_

private

lightQ jet

Definition at line 131 of file TtSemiLepJetComb.h.

Referenced by deduceMothers(), and TtSemiLepJetComb().

math::XYZTLorentzVector TtSemiLepJetComb::hadTop_

private

hadronic top 4-vector

Definition at line 143 of file TtSemiLepJetComb.h.

Referenced by deduceMothers().

math::XYZTLorentzVector TtSemiLepJetComb::hadW_

private

hadronic W 4-vector

Definition at line 145 of file TtSemiLepJetComb.h.

Referenced by deduceMothers().

pat::Jet TtSemiLepJetComb::lepBJet_

private

leptonic b jet

Definition at line 137 of file TtSemiLepJetComb.h.

Referenced by deduceMothers(), and TtSemiLepJetComb().

math::XYZTLorentzVector TtSemiLepJetComb::lepton_

private

lepton 4-vector

Definition at line 141 of file TtSemiLepJetComb.h.

Referenced by compareBLepton(), compareLeptonNeutrino(), compareTopLepton(), compareWLepton(), deduceMothers(), leptonVar(), and TtSemiLepJetComb().

math::XYZTLorentzVector TtSemiLepJetComb::lepTop_

private

leptonic top 4-vector

Definition at line 147 of file TtSemiLepJetComb.h.

Referenced by deduceMothers().

math::XYZTLorentzVector TtSemiLepJetComb::lepW_

private

leptonic W 4-vector

Definition at line 149 of file TtSemiLepJetComb.h.

Referenced by deduceMothers().

pat::MET TtSemiLepJetComb::neutrino_

private

neutrino candidate

Definition at line 139 of file TtSemiLepJetComb.h.

Referenced by compareBNeutrino(), compareLeptonNeutrino(), compareTopNeutrino(), compareWNeutrino(), deduceMothers(), neutrinoVar(), and TtSemiLepJetComb().

std::map<std::string, double> TtSemiLepJetComb::userVariables_

private

map for user defined variables

Definition at line 151 of file TtSemiLepJetComb.h.