CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Types | Private Member Functions | Private Attributes
TtDilepLRSignalSelObservables Class Reference

#include <TtDilepLRSignalSelObservables.h>

Public Types

typedef std::pair< unsigned int, bool > IntBoolPair
 

Public Member Functions

std::vector< IntBoolPair > operator() (TtDilepEvtSolution &, const edm::Event &iEvent, bool matchOnly=false)
 
 TtDilepLRSignalSelObservables (edm::ConsumesCollector &&iC, const edm::EDGetTokenT< std::vector< pat::Jet > > &jetSourceToken)
 
 ~TtDilepLRSignalSelObservables ()
 

Private Types

typedef std::pair< unsigned int, double > IntDblPair
 

Private Member Functions

double delta (double phi1, double phi2)
 
void fillMinMax (double v1, double v2, int obsNbr, std::vector< IntDblPair > &varList, bool match1, bool match2, std::vector< IntBoolPair > &matchList)
 

Private Attributes

int count1
 
int count2
 
int count3
 
int count4
 
int count5
 
std::vector< IntBoolPair > evtselectVarMatch
 
std::vector< IntDblPair > evtselectVarVal
 
edm::EDGetTokenT< TtGenEvent > genEvtToken_
 
edm::EDGetTokenT< std::vector< pat::Jet > > jetSourceToken_
 

Detailed Description

Definition at line 16 of file TtDilepLRSignalSelObservables.h.

Member Typedef Documentation

typedef std::pair<unsigned int, bool> TtDilepLRSignalSelObservables::IntBoolPair

Definition at line 22 of file TtDilepLRSignalSelObservables.h.

typedef std::pair<unsigned int, double> TtDilepLRSignalSelObservables::IntDblPair
private

Definition at line 26 of file TtDilepLRSignalSelObservables.h.

Constructor & Destructor Documentation

TtDilepLRSignalSelObservables::TtDilepLRSignalSelObservables ( edm::ConsumesCollector &&  iC,
const edm::EDGetTokenT< std::vector< pat::Jet > > &  jetSourceToken 
)

Definition at line 11 of file TtDilepLRSignalSelObservables.cc.

References count1, count2, count3, count4, and count5.

13  : jetSourceToken_(jetSourceToken), genEvtToken_(iC.consumes<TtGenEvent>(edm::InputTag("genEvt"))) {
14  count1 = 0;
15  count2 = 0;
16  count3 = 0;
17  count4 = 0;
18  count5 = 0;
19  count3 = 0;
20 }
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:49
TtGenEvent
Class derived from the TopGenEvent for ttbar events.
Definition: TtGenEvent.h:18
TtDilepLRSignalSelObservables::genEvtToken_
edm::EDGetTokenT< TtGenEvent > genEvtToken_
Definition: TtDilepLRSignalSelObservables.h:38
TtDilepLRSignalSelObservables::jetSourceToken_
edm::EDGetTokenT< std::vector< pat::Jet > > jetSourceToken_
Definition: TtDilepLRSignalSelObservables.h:37
TtDilepLRSignalSelObservables::~TtDilepLRSignalSelObservables ( )

Definition at line 23 of file TtDilepLRSignalSelObservables.cc.

23 {}

Member Function Documentation

double TtDilepLRSignalSelObservables::delta ( double  phi1,
double  phi2 
)
private

Definition at line 219 of file TtDilepLRSignalSelObservables.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, and M_PI.

Referenced by operator()().

219  {
220  double deltaPhi = phi1 - phi2;
221  while (deltaPhi > M_PI)
222  deltaPhi -= 2 * M_PI;
223  while (deltaPhi <= -M_PI)
224  deltaPhi += 2 * M_PI;
225  return deltaPhi;
226 }
SiPixelRawToDigiRegional_cfi.deltaPhi
deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
void TtDilepLRSignalSelObservables::fillMinMax ( double  v1,
double  v2,
int  obsNbr,
std::vector< IntDblPair > &  varList,
bool  match1,
bool  match2,
std::vector< IntBoolPair > &  matchList 
)
private

Definition at line 198 of file TtDilepLRSignalSelObservables.cc.

Referenced by operator()().

204  {
205  if (v1 < v2) {
206  varList.push_back(IntDblPair(obsNbr, v1));
207  varList.push_back(IntDblPair(obsNbr + 1, v2));
208  matchList.push_back(IntBoolPair(obsNbr, match1));
209  matchList.push_back(IntBoolPair(obsNbr + 1, match2));
210 
211  } else {
212  varList.push_back(IntDblPair(obsNbr, v2));
213  varList.push_back(IntDblPair(obsNbr + 1, v1));
214  matchList.push_back(IntBoolPair(obsNbr, match2));
215  matchList.push_back(IntBoolPair(obsNbr + 1, match1));
216  }
217 }
TtDilepLRSignalSelObservables::IntBoolPair
std::pair< unsigned int, bool > IntBoolPair
Definition: TtDilepLRSignalSelObservables.h:22
TtDilepLRSignalSelObservables::IntDblPair
std::pair< unsigned int, double > IntDblPair
Definition: TtDilepLRSignalSelObservables.h:26
std::vector< TtDilepLRSignalSelObservables::IntBoolPair > TtDilepLRSignalSelObservables::operator() ( TtDilepEvtSolution &  solution,
const edm::Event &  iEvent,
bool  matchOnly = false 
)

Definition at line 25 of file TtDilepLRSignalSelObservables.cc.

References funct::abs(), TopGenEvent::b(), TopGenEvent::bBar(), count1, count2, count3, count4, count5, delta(), SiPixelRawToDigiRegional_cfi::deltaPhi, flavorHistoryFilter_cfi::dr, EgHLTOffHistBins_cfi::et, reco::LeafCandidate::et(), evtselectVarMatch, evtselectVarVal, edm::HandleBase::failedToGet(), fillMinMax(), nano_cff::genEvent, genEvtToken_, edm::Event::getByToken(), TtDilepEvtSolution::getCalJetB(), TtDilepEvtSolution::getCalJetBbar(), TtDilepEvtSolution::getGenLepm(), TtDilepEvtSolution::getGenLepp(), TtDilepEvtSolution::getJetB(), TtDilepEvtSolution::getJetBbar(), TtDilepEvtSolution::getLeptNeg(), TtDilepEvtSolution::getLeptPos(), TtDilepEvtSolution::getWmDecay(), TtDilepEvtSolution::getWpDecay(), mps_fire::i, createfilelist::int, TtGenEvent::isFullLeptonic(), TtGenEvent::isSemiLeptonic(), TtGenEvent::isTtBar(), edm::HandleBase::isValid(), singleTopDQM_cfi::jets, jetSourceToken_, M_PI, EgHLTOffHistBins_cfi::mass, TopGenEvent::numberOfBQuarks(), reco::Particle::p4(), reco::LeafCandidate::p4(), pat::Jet::partonFlavour(), reco::LeafCandidate::pdgId(), createTree::pp, HLT_2018_cff::pt1, HLT_2018_cff::pt2, TtDilepEvtSolution::setLRSignalEvtObservables(), and TtGenEvent::singleLepton().

26  {
27  evtselectVarVal.clear();
28  evtselectVarMatch.clear();
29 
30  // Check whether the objects are matched:
31  bool matchB1 = false;
32  bool matchB2 = false;
33  bool matchB = false;
34  bool matchBbar = false;
35  bool matchLeptPos = false;
36  bool matchLeptNeg = false;
37 
38  edm::Handle<TtGenEvent> genEvent;
39  iEvent.getByToken(genEvtToken_, genEvent);
40 
41  if (!genEvent.failedToGet() && genEvent.isValid()) {
42  // std::cout <<std::endl;
43  double dr, dr1, dr2;
44 
45  if (genEvent->isFullLeptonic()) {
46  // Match the leptons, by type and deltaR
47  dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptPos(), *(solution.getGenLepp()));
48  matchLeptPos = ((((solution.getWpDecay() == "electron") && (std::abs(solution.getGenLepp()->pdgId()) == 11)) ||
49  ((solution.getWpDecay() == "muon") && (std::abs(solution.getGenLepp()->pdgId()) == 13))) &&
50  (dr < 0.1));
51 
52  dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptNeg(), *(solution.getGenLepm()));
53  matchLeptNeg = ((((solution.getWmDecay() == "electron") && (std::abs(solution.getGenLepm()->pdgId()) == 11)) ||
54  ((solution.getWmDecay() == "muon") && (std::abs(solution.getGenLepm()->pdgId()) == 13))) &&
55  (dr < 0.1));
56  }
57 
58  if (genEvent->isSemiLeptonic()) {
59  int id = genEvent->singleLepton()->pdgId();
60 
61  if (id > 0) {
62  dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptNeg(), *(genEvent->singleLepton()));
63  matchLeptNeg = ((((solution.getWmDecay() == "electron") && (id == 11)) ||
64  ((solution.getWmDecay() == "muon") && (id == 13))) &&
65  (dr < 0.1));
66  } else {
67  dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptPos(), *(genEvent->singleLepton()));
68  matchLeptPos = ((((solution.getWpDecay() == "electron") && (id == -11)) ||
69  ((solution.getWpDecay() == "muon") && (id == -13))) &&
70  (dr < 0.1));
71  }
72  }
73 
74  if (genEvent->isTtBar() && genEvent->numberOfBQuarks() > 1) {
75  if (solution.getJetB().partonFlavour() == 5)
76  ++count1;
77  if (solution.getJetBbar().partonFlavour() == 5)
78  ++count1;
79 
80  dr1 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetB(), *(genEvent->b()));
81  dr2 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetB(), *(genEvent->bBar()));
82 
83  matchB1 = ((dr1 < 0.4) || (dr2 < 0.4));
84  matchB = ((solution.getJetB().partonFlavour() == 5) && (dr1 < 0.4));
85  if (matchB)
86  ++count3;
87  matchB = ((dr1 < 0.4));
88  if (dr1 < 0.5)
89  ++count2;
90  if (dr1 < 0.4)
91  ++count4;
92  if (dr1 < 0.3)
93  ++count5;
94 
95  dr1 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetBbar(), *(genEvent->b()));
96  dr2 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetBbar(), *(genEvent->bBar()));
97 
98  matchBbar = ((solution.getJetBbar().partonFlavour() == 5) && (dr2 < 0.4));
99  if (matchBbar)
100  ++count3;
101  matchBbar = ((dr2 < 0.4));
102  matchB2 = ((dr1 < 0.4) || (dr2 < 0.4));
103  if (dr2 < 0.5)
104  ++count2;
105  if (dr2 < 0.4)
106  ++count4;
107  if (dr2 < 0.3)
108  ++count5;
109  }
110  }
111 
112  edm::Handle<std::vector<pat::Jet> > jets;
113  iEvent.getByToken(jetSourceToken_, jets);
114 
115  // Lower / Higher of both jet angles
116 
117  double v1 = std::abs(solution.getJetB().p4().theta() - M_PI / 2);
118  double v2 = std::abs(solution.getJetBbar().p4().theta() - M_PI / 2);
119  fillMinMax(v1, v2, 1, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);
120 
121  // Lower / Higher of both jet pT
122 
123  double pt1 = solution.getJetB().p4().pt();
124  double pt2 = solution.getJetBbar().p4().pt();
125  fillMinMax(pt1, pt2, 3, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);
126 
127  // Lower / Higher of both lepton pT
128 
129  pt1 = solution.getLeptPos().p4().pt();
130  pt2 = solution.getLeptNeg().p4().pt();
131  fillMinMax(pt1, pt2, 5, evtselectVarVal, matchLeptPos, matchLeptNeg, evtselectVarMatch);
132 
133  // delta theta btw the b-jets
134 
135  double deltaPhi = std::abs(delta(solution.getJetB().p4().phi(), solution.getJetBbar().p4().phi()));
136 
137  evtselectVarVal.push_back(IntDblPair(7, deltaPhi));
138  evtselectVarMatch.push_back(IntBoolPair(7, matchB1 && matchB2));
139 
140  // delta phi btw the b-jets
141 
142  double deltaTheta = std::abs(delta(solution.getJetBbar().p4().theta(), solution.getJetB().p4().theta()));
143 
144  evtselectVarVal.push_back(IntDblPair(8, deltaTheta));
145  evtselectVarMatch.push_back(IntBoolPair(8, matchB1 && matchB2));
146 
147  // Lower / Higher of phi difference between the b and associated lepton
148 
149  double deltaPhi1 = std::abs(delta(solution.getJetB().p4().phi(), solution.getLeptPos().p4().phi()));
150  double deltaPhi2 = std::abs(delta(solution.getJetBbar().p4().phi(), solution.getLeptNeg().p4().phi()));
151 
152  fillMinMax(
153  deltaPhi1, deltaPhi2, 9, evtselectVarVal, matchB && matchLeptPos, matchBbar && matchLeptNeg, evtselectVarMatch);
154 
155  // Lower / Higher of theta difference between the b and associated lepton
156 
157  double deltaTheta1 = std::abs(solution.getJetB().p4().theta() - solution.getLeptPos().p4().theta());
158  double deltaTheta2 = std::abs(solution.getJetBbar().p4().theta() - solution.getLeptNeg().p4().theta());
159  fillMinMax(deltaTheta1,
160  deltaTheta2,
161  11,
162  evtselectVarVal,
163  matchB && matchLeptPos,
164  matchBbar && matchLeptNeg,
165  evtselectVarMatch);
166 
167  // Invariant Mass of lepton pair
168 
169  math::XYZTLorentzVector pp = solution.getLeptPos().p4() + solution.getLeptNeg().p4();
170  double mass = pp.mass();
171  evtselectVarVal.push_back(IntDblPair(13, mass));
172  evtselectVarMatch.push_back(IntBoolPair(13, matchLeptNeg && matchLeptPos));
173 
174  evtselectVarVal.push_back(IntDblPair(13, mass));
175  evtselectVarMatch.push_back(IntBoolPair(13, matchLeptNeg && matchLeptPos));
176 
177  std::vector<pat::Jet> jet3;
178  for (int i = 0; i < (int)jets->size(); ++i) {
179  if (((*jets)[i].et() < solution.getJetB().et()) && ((*jets)[i].et() < solution.getJetBbar().et())) {
180  jet3.push_back((*jets)[i]);
181  }
182  }
183  double jet1Ratio = 0., jet2Ratio = 0.;
184  if (!jet3.empty()) {
185  jet1Ratio = jet3[0].et() / solution.getJetB().et();
186  jet2Ratio = jet3[0].et() / solution.getJetBbar().et();
187  }
188  fillMinMax(jet1Ratio, jet2Ratio, 14, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);
189 
190  evtselectVarVal.push_back(IntDblPair(16, jets->size()));
191  evtselectVarMatch.push_back(IntBoolPair(16, matchB && matchBbar));
192 
193  if (!matchOnly)
194  solution.setLRSignalEvtObservables(evtselectVarVal);
195  return evtselectVarMatch;
196 }
TtDilepEvtSolution::getCalJetBbar
pat::Jet getCalJetBbar() const
Definition: TtDilepEvtSolution.h:119
TtDilepEvtSolution::getGenLepm
const reco::GenParticle * getGenLepm() const
Definition: TtDilepEvtSolution.h:100
TtDilepLRSignalSelObservables::delta
double delta(double phi1, double phi2)
Definition: TtDilepLRSignalSelObservables.cc:219
TopGenEvent::b
const reco::GenParticle * b() const
return b quark if available; 0 else
Definition: TopGenEvent.h:97
reco::LeafCandidate::pdgId
int pdgId() const final
PDG identifier.
Definition: LeafCandidate.h:176
EgHLTOffHistBins_cfi.et
et
Definition: EgHLTOffHistBins_cfi.py:8
mps_fire.i
i
Definition: mps_fire.py:341
TtGenEvent::isSemiLeptonic
bool isSemiLeptonic(bool excludeTauLeptons=false) const
check if the event can be classified as semi-laptonic
Definition: TtGenEvent.h:38
TtDilepLRSignalSelObservables::fillMinMax
void fillMinMax(double v1, double v2, int obsNbr, std::vector< IntDblPair > &varList, bool match1, bool match2, std::vector< IntBoolPair > &matchList)
Definition: TtDilepLRSignalSelObservables.cc:198
DeltaR
Definition: deltaR.h:58
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
nano_cff.genEvent
genEvent
Definition: nano_cff.py:97
HLT_2018_cff.pt2
pt2
Definition: HLT_2018_cff.py:8539
SiPixelRawToDigiRegional_cfi.deltaPhi
deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
TopGenEvent::bBar
const reco::GenParticle * bBar() const
return anti-b quark if available; 0 else
Definition: TopGenEvent.h:99
reco::Particle::p4
const LorentzVector & p4() const
four-momentum Lorentz vector
Definition: Particle.h:88
TtDilepEvtSolution::getWmDecay
std::string getWmDecay() const
Definition: TtDilepEvtSolution.h:125
TtDilepEvtSolution::getWpDecay
std::string getWpDecay() const
Definition: TtDilepEvtSolution.h:124
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
TtGenEvent::isFullLeptonic
bool isFullLeptonic(bool excludeTauLeptons=false) const
check if the event can be classified as full leptonic
Definition: TtGenEvent.h:42
TtDilepLRSignalSelObservables::IntBoolPair
std::pair< unsigned int, bool > IntBoolPair
Definition: TtDilepLRSignalSelObservables.h:22
TtDilepEvtSolution::setLRSignalEvtObservables
void setLRSignalEvtObservables(const std::vector< std::pair< unsigned int, double > > &)
Definition: TtDilepEvtSolution.cc:148
reco::LeafCandidate::et
double et() const final
transverse energy
Definition: LeafCandidate.h:127
TtDilepLRSignalSelObservables::IntDblPair
std::pair< unsigned int, double > IntDblPair
Definition: TtDilepLRSignalSelObservables.h:26
TtDilepEvtSolution::getLeptNeg
reco::Particle getLeptNeg() const
Definition: TtDilepEvtSolution.cc:104
TtDilepEvtSolution::getJetBbar
pat::Jet getJetBbar() const
Definition: TtDilepEvtSolution.cc:35
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
TtDilepLRSignalSelObservables::genEvtToken_
edm::EDGetTokenT< TtGenEvent > genEvtToken_
Definition: TtDilepLRSignalSelObservables.h:38
reco::LeafCandidate::p4
const LorentzVector & p4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:114
pat::Jet::partonFlavour
int partonFlavour() const
return the parton-based flavour of the jet
TtDilepEvtSolution::getLeptPos
reco::Particle getLeptPos() const
Definition: TtDilepEvtSolution.cc:52
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
edm::HandleBase::failedToGet
bool failedToGet() const
Definition: HandleBase.h:72
TopGenEvent::numberOfBQuarks
int numberOfBQuarks(bool fromTopQuark=true) const
return number of b quarks in the decay chain
Definition: TopGenEvent.cc:99
TtDilepEvtSolution::getJetB
pat::Jet getJetB() const
Definition: TtDilepEvtSolution.cc:22
TtDilepLRSignalSelObservables::jetSourceToken_
edm::EDGetTokenT< std::vector< pat::Jet > > jetSourceToken_
Definition: TtDilepLRSignalSelObservables.h:37
TtDilepEvtSolution::getCalJetB
pat::Jet getCalJetB() const
Definition: TtDilepEvtSolution.h:117
TtDilepEvtSolution::getGenLepp
const reco::GenParticle * getGenLepp() const
Definition: TtDilepEvtSolution.h:70
HLT_2018_cff.pt1
pt1
Definition: HLT_2018_cff.py:8537
singleTopDQM_cfi.jets
jets
Definition: singleTopDQM_cfi.py:45
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
TtDilepLRSignalSelObservables::evtselectVarMatch
std::vector< IntBoolPair > evtselectVarMatch
Definition: TtDilepLRSignalSelObservables.h:41
TtGenEvent::isTtBar
bool isTtBar() const
check if the event can be classified as ttbar
Definition: TtGenEvent.h:28
TtGenEvent::singleLepton
const reco::GenParticle * singleLepton(bool excludeTauLeptons=false) const
return single lepton if available; 0 else
Definition: TtGenEvent.cc:97

Member Data Documentation

int TtDilepLRSignalSelObservables::count1
private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

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

int TtDilepLRSignalSelObservables::count2
private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

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

int TtDilepLRSignalSelObservables::count3
private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

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

int TtDilepLRSignalSelObservables::count4
private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

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

int TtDilepLRSignalSelObservables::count5
private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

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

std::vector<IntBoolPair> TtDilepLRSignalSelObservables::evtselectVarMatch
private

Definition at line 41 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

std::vector<IntDblPair> TtDilepLRSignalSelObservables::evtselectVarVal
private

Definition at line 40 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

edm::EDGetTokenT<TtGenEvent> TtDilepLRSignalSelObservables::genEvtToken_
private

Definition at line 38 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

edm::EDGetTokenT<std::vector<pat::Jet> > TtDilepLRSignalSelObservables::jetSourceToken_
private

Definition at line 37 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

Generated for CMSSW Reference Manual by doxygen 1.8.11