CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
TtFullLepKinSolutionProducer Class Reference
Inheritance diagram for TtFullLepKinSolutionProducer:
edm::stream::EDProducer<>

Classes

struct  Compare
 

Public Member Functions

void produce (edm::Event &evt, const edm::EventSetup &iSetup) override
 
 TtFullLepKinSolutionProducer (const edm::ParameterSet &iConfig)
 
 ~TtFullLepKinSolutionProducer () override=default
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Private Member Functions

bool HasPositiveCharge (const reco::Candidate *) const
 
bool LepDiffCharge (const reco::Candidate *, const reco::Candidate *) const
 
bool PTComp (const reco::Candidate *, const reco::Candidate *) const
 

Private Attributes

bool eeChannel_
 
edm::EDGetTokenT< std::vector< pat::Electron > > electronsToken_
 
bool emuChannel_
 
std::string jetCorrLevel_
 
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
 
int maxNComb_
 
int maxNJets_
 
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
 
bool mumuChannel_
 
edm::EDGetTokenT< std::vector< pat::Muon > > muonsToken_
 
std::vector< double > nupars_
 
bool searchWrongCharge_
 
std::unique_ptr< TtFullLepKinSolver > solver
 
double tmassbegin_
 
double tmassend_
 
double tmassstep_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
using CacheTypes = CacheContexts< T... >
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T... >
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 

Detailed Description

Definition at line 14 of file TtFullLepKinSolutionProducer.cc.

Constructor & Destructor Documentation

◆ TtFullLepKinSolutionProducer()

TtFullLepKinSolutionProducer::TtFullLepKinSolutionProducer ( const edm::ParameterSet &  iConfig)
inlineexplicit

Definition at line 57 of file TtFullLepKinSolutionProducer.cc.

References eeChannel_, electronsToken_, emuChannel_, edm::ParameterSet::getParameter(), jetCorrLevel_, jetsToken_, maxNComb_, maxNJets_, metsToken_, mumuChannel_, muonsToken_, nupars_, searchWrongCharge_, solver, AlCaHLTBitMon_QueryRunRegistry::string, tmassbegin_, tmassend_, and tmassstep_.

57  {
58  // configurables
59  jetsToken_ = consumes<std::vector<pat::Jet> >(iConfig.getParameter<edm::InputTag>("jets"));
60  electronsToken_ = consumes<std::vector<pat::Electron> >(iConfig.getParameter<edm::InputTag>("electrons"));
61  muonsToken_ = consumes<std::vector<pat::Muon> >(iConfig.getParameter<edm::InputTag>("muons"));
62  metsToken_ = consumes<std::vector<pat::MET> >(iConfig.getParameter<edm::InputTag>("mets"));
63  jetCorrLevel_ = iConfig.getParameter<std::string>("jetCorrectionLevel");
64  maxNJets_ = iConfig.getParameter<int>("maxNJets");
65  maxNComb_ = iConfig.getParameter<int>("maxNComb");
66  eeChannel_ = iConfig.getParameter<bool>("eeChannel");
67  emuChannel_ = iConfig.getParameter<bool>("emuChannel");
68  mumuChannel_ = iConfig.getParameter<bool>("mumuChannel");
69  searchWrongCharge_ = iConfig.getParameter<bool>("searchWrongCharge");
70  tmassbegin_ = iConfig.getParameter<double>("tmassbegin");
71  tmassend_ = iConfig.getParameter<double>("tmassend");
72  tmassstep_ = iConfig.getParameter<double>("tmassstep");
73  nupars_ = iConfig.getParameter<std::vector<double> >("neutrino_parameters");
74 
75  // define what will be produced
76  produces<std::vector<std::vector<int> > >(); // vector of the particle inices (b, bbar, e1, e2, mu1, mu2)
77  produces<std::vector<reco::LeafCandidate> >("fullLepNeutrinos");
78  produces<std::vector<reco::LeafCandidate> >("fullLepNeutrinoBars");
79  produces<std::vector<double> >("solWeight"); //weight for a specific kin solution
80  produces<bool>("isWrongCharge"); //true if leptons have the same charge
81 
82  solver = std::make_unique<TtFullLepKinSolver>(tmassbegin_, tmassend_, tmassstep_, nupars_);
83 }
TtFullLepKinSolutionProducer::muonsToken_
edm::EDGetTokenT< std::vector< pat::Muon > > muonsToken_
Definition: TtFullLepKinSolutionProducer.cc:33
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
TtFullLepKinSolutionProducer::solver
std::unique_ptr< TtFullLepKinSolver > solver
Definition: TtFullLepKinSolutionProducer.cc:42
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TtFullLepKinSolutionProducer::metsToken_
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
Definition: TtFullLepKinSolutionProducer.cc:34
TtFullLepKinSolutionProducer::jetsToken_
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
Definition: TtFullLepKinSolutionProducer.cc:31
TtFullLepKinSolutionProducer::electronsToken_
edm::EDGetTokenT< std::vector< pat::Electron > > electronsToken_
Definition: TtFullLepKinSolutionProducer.cc:32

◆ ~TtFullLepKinSolutionProducer()

TtFullLepKinSolutionProducer::~TtFullLepKinSolutionProducer ( )
overridedefault

Member Function Documentation

◆ HasPositiveCharge()

bool TtFullLepKinSolutionProducer::HasPositiveCharge ( const reco::Candidate *  l) const
inlineprivate

Definition at line 53 of file TtFullLepKinSolutionProducer.cc.

References cmsLHEtoEOSManager::l.

Referenced by produce().

53  {
54  return (l->charge() > 0);
55 }

◆ LepDiffCharge()

bool TtFullLepKinSolutionProducer::LepDiffCharge ( const reco::Candidate *  l1,
const reco::Candidate *  l2 
) const
inlineprivate

Definition at line 49 of file TtFullLepKinSolutionProducer.cc.

References reco::Candidate::charge().

Referenced by produce().

49  {
50  return (l1->charge() != l2->charge());
51 }
reco::Candidate::charge
virtual int charge() const =0
electric charge

◆ produce()

void TtFullLepKinSolutionProducer::produce ( edm::Event &  evt,
const edm::EventSetup &  iSetup 
)
inlineoverride

Definition at line 85 of file TtFullLepKinSolutionProducer.cc.

References eeChannel_, pwdgSkimBPark_cfi::electrons, electronsToken_, emuChannel_, HCALHighEnergyHPDFilter_cfi::energy, dqmdumpme::first, edm::Event::getHandle(), HasPositiveCharge(), mps_fire::i, cuy::ib, createfilelist::int, jetCorrLevel_, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, LepDiffCharge(), maxNComb_, maxNJets_, singleTopDQM_cfi::mets, metsToken_, eostools::move(), mumuChannel_, PDWG_BPHSkim_cff::muons, muonsToken_, TtFullLepKinSolver::NeutrinoSolution::neutrino, TtFullLepKinSolver::NeutrinoSolution::neutrinoBar, PTComp(), edm::Event::put(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, searchWrongCharge_, edm::second(), reco::Candidate::size, solver, jetUpdater_cfi::sort, trackerHitRTTI::vector, and TtFullLepKinSolver::NeutrinoSolution::weight.

85  {
86  //create vectors fo runsorted output
87  std::vector<std::vector<int> > idcsV;
88  std::vector<reco::LeafCandidate> nusV;
89  std::vector<reco::LeafCandidate> nuBarsV;
90  std::vector<std::pair<double, int> > weightsV;
91 
92  //create pointer for products
93  std::unique_ptr<std::vector<std::vector<int> > > pIdcs(new std::vector<std::vector<int> >);
94  std::unique_ptr<std::vector<reco::LeafCandidate> > pNus(new std::vector<reco::LeafCandidate>);
95  std::unique_ptr<std::vector<reco::LeafCandidate> > pNuBars(new std::vector<reco::LeafCandidate>);
96  std::unique_ptr<std::vector<double> > pWeight(new std::vector<double>);
97  std::unique_ptr<bool> pWrongCharge(new bool);
98 
99  const edm::Handle<std::vector<pat::Jet> >& jets = evt.getHandle(jetsToken_);
100  const edm::Handle<std::vector<pat::Electron> >& electrons = evt.getHandle(electronsToken_);
101  const edm::Handle<std::vector<pat::Muon> >& muons = evt.getHandle(muonsToken_);
102  const edm::Handle<std::vector<pat::MET> >& mets = evt.getHandle(metsToken_);
103 
104  int selMuon1 = -1, selMuon2 = -1;
105  int selElectron1 = -1, selElectron2 = -1;
106  bool ee = false;
107  bool emu = false;
108  bool mumu = false;
109  bool isWrongCharge = false;
110  bool jetsFound = false;
111  bool METFound = false;
112  bool electronsFound = false;
113  bool electronMuonFound = false;
114  bool muonsFound = false;
115 
116  //select Jets (TopJet vector is sorted on ET)
117  if (jets->size() >= 2) {
118  jetsFound = true;
119  }
120 
121  //select MET (TopMET vector is sorted on ET)
122  if (!mets->empty()) {
123  METFound = true;
124  }
125 
126  // If we have electrons and muons available,
127  // build a solutions with electrons and muons.
128  if (muons->size() + electrons->size() >= 2) {
129  // select leptons
130  if (electrons->empty())
131  mumu = true;
132  else if (muons->empty())
133  ee = true;
134  else if (electrons->size() == 1) {
135  if (muons->size() == 1)
136  emu = true;
137  else if (PTComp(&(*electrons)[0], &(*muons)[1]))
138  emu = true;
139  else
140  mumu = true;
141  } else if (electrons->size() > 1) {
142  if (PTComp(&(*electrons)[1], &(*muons)[0]))
143  ee = true;
144  else if (muons->size() == 1)
145  emu = true;
146  else if (PTComp(&(*electrons)[0], &(*muons)[1]))
147  emu = true;
148  else
149  mumu = true;
150  }
151  if (ee && eeChannel_) {
152  if (LepDiffCharge(&(*electrons)[0], &(*electrons)[1]) || searchWrongCharge_) {
153  if (HasPositiveCharge(&(*electrons)[0]) || !LepDiffCharge(&(*electrons)[0], &(*electrons)[1])) {
154  selElectron1 = 0;
155  selElectron2 = 1;
156  } else {
157  selElectron1 = 1;
158  selElectron2 = 0;
159  }
160  electronsFound = true;
161  if (!LepDiffCharge(&(*electrons)[0], &(*electrons)[1]))
162  isWrongCharge = true;
163  }
164  } else if (emu && emuChannel_) {
165  if (LepDiffCharge(&(*electrons)[0], &(*muons)[0]) || searchWrongCharge_) {
166  selElectron1 = 0;
167  selMuon1 = 0;
168  electronMuonFound = true;
169  if (!LepDiffCharge(&(*electrons)[0], &(*muons)[0]))
170  isWrongCharge = true;
171  }
172  } else if (mumu && mumuChannel_) {
173  if (LepDiffCharge(&(*muons)[0], &(*muons)[1]) || searchWrongCharge_) {
174  if (HasPositiveCharge(&(*muons)[0]) || !LepDiffCharge(&(*muons)[0], &(*muons)[1])) {
175  selMuon1 = 0;
176  selMuon2 = 1;
177  } else {
178  selMuon1 = 1;
179  selMuon2 = 0;
180  }
181  muonsFound = true;
182  if (!LepDiffCharge(&(*muons)[0], &(*muons)[1]))
183  isWrongCharge = true;
184  }
185  }
186  }
187 
188  *pWrongCharge = isWrongCharge;
189 
190  // Check that the above work makes sense
191  if (int(ee) + int(emu) + int(mumu) > 1) {
192  edm::LogWarning("TtFullLepKinSolutionProducer") << "Lepton selection criteria uncorrectly defined";
193  }
194 
195  // Check if the leptons for the required Channel are available
196  bool correctLeptons =
197  ((electronsFound && eeChannel_) || (muonsFound && mumuChannel_) || (electronMuonFound && emuChannel_));
198  // Check for equally charged leptons if for wrong charge combinations is searched
199  if (isWrongCharge) {
200  correctLeptons = correctLeptons && searchWrongCharge_;
201  }
202 
203  if (correctLeptons && METFound && jetsFound) {
204  // run over all jets if input parameter maxNJets is -1 or
205  // adapt maxNJets if number of present jets is smaller than selected
206  // number of jets
207  int stop = maxNJets_;
208  if (jets->size() < static_cast<unsigned int>(stop) || stop < 0)
209  stop = jets->size();
210 
211  // counter for number of found kinematic solutions
212  int nSol = 0;
213 
214  // consider all permutations
215  for (int ib = 0; ib < stop; ib++) {
216  // second loop of the permutations
217  for (int ibbar = 0; ibbar < stop; ibbar++) {
218  // avoid the diagonal: b and bbar must be distinct jets
219  if (ib == ibbar)
220  continue;
221 
222  std::vector<int> idcs;
223 
224  // push back the indices of the jets
225  idcs.push_back(ib);
226  idcs.push_back(ibbar);
227 
228  TLorentzVector LV_l1;
229  TLorentzVector LV_l2;
230  TLorentzVector LV_b = TLorentzVector((*jets)[ib].correctedJet(jetCorrLevel_, "bottom").px(),
231  (*jets)[ib].correctedJet(jetCorrLevel_, "bottom").py(),
232  (*jets)[ib].correctedJet(jetCorrLevel_, "bottom").pz(),
233  (*jets)[ib].correctedJet(jetCorrLevel_, "bottom").energy());
234  TLorentzVector LV_bbar = TLorentzVector((*jets)[ibbar].correctedJet(jetCorrLevel_, "bottom").px(),
235  (*jets)[ibbar].correctedJet(jetCorrLevel_, "bottom").py(),
236  (*jets)[ibbar].correctedJet(jetCorrLevel_, "bottom").pz(),
237  (*jets)[ibbar].correctedJet(jetCorrLevel_, "bottom").energy());
238 
239  double xconstraint = 0, yconstraint = 0;
240 
241  if (ee) {
242  idcs.push_back(selElectron1);
243  LV_l1.SetXYZT((*electrons)[selElectron1].px(),
244  (*electrons)[selElectron1].py(),
245  (*electrons)[selElectron1].pz(),
246  (*electrons)[selElectron1].energy());
247  xconstraint += (*electrons)[selElectron1].px();
248  yconstraint += (*electrons)[selElectron1].py();
249 
250  idcs.push_back(selElectron2);
251  LV_l2.SetXYZT((*electrons)[selElectron2].px(),
252  (*electrons)[selElectron2].py(),
253  (*electrons)[selElectron2].pz(),
254  (*electrons)[selElectron2].energy());
255  xconstraint += (*electrons)[selElectron2].px();
256  yconstraint += (*electrons)[selElectron2].py();
257 
258  idcs.push_back(-1);
259  idcs.push_back(-1);
260  }
261 
262  else if (emu) {
263  if (!isWrongCharge) {
264  if (HasPositiveCharge(&(*electrons)[selElectron1])) {
265  idcs.push_back(selElectron1);
266  LV_l1.SetXYZT((*electrons)[selElectron1].px(),
267  (*electrons)[selElectron1].py(),
268  (*electrons)[selElectron1].pz(),
269  (*electrons)[selElectron1].energy());
270  xconstraint += (*electrons)[selElectron1].px();
271  yconstraint += (*electrons)[selElectron1].py();
272 
273  idcs.push_back(-1);
274  idcs.push_back(-1);
275 
276  idcs.push_back(selMuon1);
277  LV_l2.SetXYZT((*muons)[selMuon1].px(),
278  (*muons)[selMuon1].py(),
279  (*muons)[selMuon1].pz(),
280  (*muons)[selMuon1].energy());
281  xconstraint += (*muons)[selMuon1].px();
282  yconstraint += (*muons)[selMuon1].py();
283  } else {
284  idcs.push_back(-1);
285 
286  idcs.push_back(selMuon1);
287  LV_l1.SetXYZT((*muons)[selMuon1].px(),
288  (*muons)[selMuon1].py(),
289  (*muons)[selMuon1].pz(),
290  (*muons)[selMuon1].energy());
291  xconstraint += (*muons)[selMuon1].px();
292  yconstraint += (*muons)[selMuon1].py();
293 
294  idcs.push_back(selElectron1);
295  LV_l2.SetXYZT((*electrons)[selElectron1].px(),
296  (*electrons)[selElectron1].py(),
297  (*electrons)[selElectron1].pz(),
298  (*electrons)[selElectron1].energy());
299  xconstraint += (*electrons)[selElectron1].px();
300  yconstraint += (*electrons)[selElectron1].py();
301 
302  idcs.push_back(-1);
303  }
304  } else { // means "if wrong charge"
305  if (HasPositiveCharge(&(*electrons)[selElectron1])) { // both leps positive
306  idcs.push_back(selElectron1);
307  LV_l1.SetXYZT((*electrons)[selElectron1].px(),
308  (*electrons)[selElectron1].py(),
309  (*electrons)[selElectron1].pz(),
310  (*electrons)[selElectron1].energy());
311  xconstraint += (*electrons)[selElectron1].px();
312  yconstraint += (*electrons)[selElectron1].py();
313 
314  idcs.push_back(-1);
315 
316  idcs.push_back(selMuon1);
317  LV_l2.SetXYZT((*muons)[selMuon1].px(),
318  (*muons)[selMuon1].py(),
319  (*muons)[selMuon1].pz(),
320  (*muons)[selMuon1].energy());
321  xconstraint += (*muons)[selMuon1].px();
322  yconstraint += (*muons)[selMuon1].py();
323 
324  idcs.push_back(-1);
325  } else { // both leps negative
326  idcs.push_back(-1);
327 
328  idcs.push_back(selElectron1);
329  LV_l2.SetXYZT((*electrons)[selElectron1].px(),
330  (*electrons)[selElectron1].py(),
331  (*electrons)[selElectron1].pz(),
332  (*electrons)[selElectron1].energy());
333  xconstraint += (*electrons)[selElectron1].px();
334  yconstraint += (*electrons)[selElectron1].py();
335 
336  idcs.push_back(-1);
337 
338  idcs.push_back(selMuon1);
339  LV_l1.SetXYZT((*muons)[selMuon1].px(),
340  (*muons)[selMuon1].py(),
341  (*muons)[selMuon1].pz(),
342  (*muons)[selMuon1].energy());
343  xconstraint += (*muons)[selMuon1].px();
344  yconstraint += (*muons)[selMuon1].py();
345  }
346  }
347  }
348 
349  else if (mumu) {
350  idcs.push_back(-1);
351  idcs.push_back(-1);
352 
353  idcs.push_back(selMuon1);
354  LV_l1.SetXYZT(
355  (*muons)[selMuon1].px(), (*muons)[selMuon1].py(), (*muons)[selMuon1].pz(), (*muons)[selMuon1].energy());
356  xconstraint += (*muons)[selMuon1].px();
357  yconstraint += (*muons)[selMuon1].py();
358 
359  idcs.push_back(selMuon2);
360  LV_l2.SetXYZT(
361  (*muons)[selMuon2].px(), (*muons)[selMuon2].py(), (*muons)[selMuon2].pz(), (*muons)[selMuon2].energy());
362  xconstraint += (*muons)[selMuon2].px();
363  yconstraint += (*muons)[selMuon2].py();
364  }
365 
366  xconstraint += (*jets)[ib].px() + (*jets)[ibbar].px() + (*mets)[0].px();
367  yconstraint += (*jets)[ib].py() + (*jets)[ibbar].py() + (*mets)[0].py();
368 
369  // calculate neutrino momenta and eventweight
370  solver->SetConstraints(xconstraint, yconstraint);
371  TtFullLepKinSolver::NeutrinoSolution nuSol = solver->getNuSolution(LV_l1, LV_l2, LV_b, LV_bbar);
372 
373  // add solution to the vectors of solutions if solution exists
374  if (nuSol.weight > 0) {
375  // add the leptons and jets indices to the vector of combinations
376  idcsV.push_back(idcs);
377 
378  // add the neutrinos
379  nusV.push_back(nuSol.neutrino);
380  nuBarsV.push_back(nuSol.neutrinoBar);
381 
382  // add the solution weight
383  weightsV.push_back(std::make_pair(nuSol.weight, nSol));
384 
385  nSol++;
386  }
387  }
388  }
389  }
390 
391  if (weightsV.empty()) {
392  //create dmummy vector
393  std::vector<int> idcs;
394  idcs.reserve(6);
395  for (int i = 0; i < 6; ++i)
396  idcs.push_back(-1);
397 
398  idcsV.push_back(idcs);
399  weightsV.push_back(std::make_pair(-1, 0));
400  reco::LeafCandidate nu;
401  nusV.push_back(nu);
402  reco::LeafCandidate nuBar;
403  nuBarsV.push_back(nuBar);
404  }
405 
406  // check if all vectors have correct length
407  int weightL = weightsV.size();
408  int nuL = nusV.size();
409  int nuBarL = nuBarsV.size();
410  int idxL = idcsV.size();
411 
412  if (weightL != nuL || weightL != nuBarL || weightL != idxL) {
413  edm::LogWarning("TtFullLepKinSolutionProducer")
414  << "Output vectors are of different length:"
415  << "\n weight: " << weightL << "\n nu: " << nuL << "\n nubar: " << nuBarL << "\n idcs: " << idxL;
416  }
417 
418  // sort vectors by weight in decreasing order
419  if (weightsV.size() > 1) {
420  sort(weightsV.begin(), weightsV.end(), Compare());
421  }
422 
423  // determine the number of solutions which is written in the event
424  int stop = weightL;
425  if (maxNComb_ > 0 && maxNComb_ < stop)
426  stop = maxNComb_;
427 
428  for (int i = 0; i < stop; ++i) {
429  pWeight->push_back(weightsV[i].first);
430  pNus->push_back(nusV[weightsV[i].second]);
431  pNuBars->push_back(nuBarsV[weightsV[i].second]);
432  pIdcs->push_back(idcsV[weightsV[i].second]);
433  }
434 
435  // put the results in the event
436  evt.put(std::move(pIdcs));
437  evt.put(std::move(pNus), "fullLepNeutrinos");
438  evt.put(std::move(pNuBars), "fullLepNeutrinoBars");
439  evt.put(std::move(pWeight), "solWeight");
440  evt.put(std::move(pWrongCharge), "isWrongCharge");
441 }
TtFullLepKinSolutionProducer::muonsToken_
edm::EDGetTokenT< std::vector< pat::Muon > > muonsToken_
Definition: TtFullLepKinSolutionProducer.cc:33
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
mps_fire.i
i
Definition: mps_fire.py:429
TtFullLepKinSolutionProducer::solver
std::unique_ptr< TtFullLepKinSolver > solver
Definition: TtFullLepKinSolutionProducer.cc:42
PDWG_EXODelayedJetMET_cff.jets
jets
Definition: PDWG_EXODelayedJetMET_cff.py:14
singleTopDQM_cfi.mets
mets
Definition: singleTopDQM_cfi.py:43
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
TtFullLepKinSolutionProducer::LepDiffCharge
bool LepDiffCharge(const reco::Candidate *, const reco::Candidate *) const
Definition: TtFullLepKinSolutionProducer.cc:49
TtFullLepKinSolutionProducer::metsToken_
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
Definition: TtFullLepKinSolutionProducer.cc:34
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
TtFullLepKinSolutionProducer::jetsToken_
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
Definition: TtFullLepKinSolutionProducer.cc:31
TtFullLepKinSolutionProducer::HasPositiveCharge
bool HasPositiveCharge(const reco::Candidate *) const
Definition: TtFullLepKinSolutionProducer.cc:53
pwdgSkimBPark_cfi.electrons
electrons
Definition: pwdgSkimBPark_cfi.py:6
multPhiCorr_741_25nsDY_cfi.py
py
Definition: multPhiCorr_741_25nsDY_cfi.py:12
TtFullLepKinSolutionProducer::electronsToken_
edm::EDGetTokenT< std::vector< pat::Electron > > electronsToken_
Definition: TtFullLepKinSolutionProducer.cc:32
multPhiCorr_741_25nsDY_cfi.px
px
Definition: multPhiCorr_741_25nsDY_cfi.py:10
edm::Event::getHandle
Handle< PROD > getHandle(EDGetTokenT< PROD > token) const
Definition: Event.h:564
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
TtFullLepKinSolutionProducer::PTComp
bool PTComp(const reco::Candidate *, const reco::Candidate *) const
Definition: TtFullLepKinSolutionProducer.cc:45
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
eostools.move
def move(src, dest)
Definition: eostools.py:511
cuy.ib
ib
Definition: cuy.py:661

◆ PTComp()

bool TtFullLepKinSolutionProducer::PTComp ( const reco::Candidate *  l1,
const reco::Candidate *  l2 
) const
inlineprivate

Definition at line 45 of file TtFullLepKinSolutionProducer.cc.

References reco::Candidate::pt().

Referenced by produce().

45  {
46  return (l1->pt() > l2->pt());
47 }
reco::Candidate::pt
virtual double pt() const =0
transverse momentum

Member Data Documentation

◆ eeChannel_

bool TtFullLepKinSolutionProducer::eeChannel_
private

Definition at line 38 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ electronsToken_

edm::EDGetTokenT<std::vector<pat::Electron> > TtFullLepKinSolutionProducer::electronsToken_
private

Definition at line 32 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ emuChannel_

bool TtFullLepKinSolutionProducer::emuChannel_
private

Definition at line 38 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ jetCorrLevel_

std::string TtFullLepKinSolutionProducer::jetCorrLevel_
private

Definition at line 36 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ jetsToken_

edm::EDGetTokenT<std::vector<pat::Jet> > TtFullLepKinSolutionProducer::jetsToken_
private

Definition at line 31 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ maxNComb_

int TtFullLepKinSolutionProducer::maxNComb_
private

Definition at line 37 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ maxNJets_

int TtFullLepKinSolutionProducer::maxNJets_
private

Definition at line 37 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ metsToken_

edm::EDGetTokenT<std::vector<pat::MET> > TtFullLepKinSolutionProducer::metsToken_
private

Definition at line 34 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ mumuChannel_

bool TtFullLepKinSolutionProducer::mumuChannel_
private

Definition at line 38 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ muonsToken_

edm::EDGetTokenT<std::vector<pat::Muon> > TtFullLepKinSolutionProducer::muonsToken_
private

Definition at line 33 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ nupars_

std::vector<double> TtFullLepKinSolutionProducer::nupars_
private

Definition at line 40 of file TtFullLepKinSolutionProducer.cc.

Referenced by TtFullLepKinSolutionProducer().

◆ searchWrongCharge_

bool TtFullLepKinSolutionProducer::searchWrongCharge_
private

Definition at line 38 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ solver

std::unique_ptr<TtFullLepKinSolver> TtFullLepKinSolutionProducer::solver
private

Definition at line 42 of file TtFullLepKinSolutionProducer.cc.

Referenced by produce(), and TtFullLepKinSolutionProducer().

◆ tmassbegin_

double TtFullLepKinSolutionProducer::tmassbegin_
private

Definition at line 39 of file TtFullLepKinSolutionProducer.cc.

Referenced by TtFullLepKinSolutionProducer().

◆ tmassend_

double TtFullLepKinSolutionProducer::tmassend_
private

Definition at line 39 of file TtFullLepKinSolutionProducer.cc.

Referenced by TtFullLepKinSolutionProducer().

◆ tmassstep_

double TtFullLepKinSolutionProducer::tmassstep_
private

Definition at line 39 of file TtFullLepKinSolutionProducer.cc.

Referenced by TtFullLepKinSolutionProducer().

Generated for CMSSW Reference Manual by doxygen 1.8.14