CMS 3D CMS Logo
ZMuMuEfficiency.cc
Go to the documentation of this file.
1 /* \class ZMuMuEfficiency
2  *
3  * author: Pasquale Noli
4  * revised by Salvatore di Guida
5  * revised for CSA08 by Davide Piccolo
6  *
7  * Efficiency of reconstruction tracker and muon Chamber
8  *
9  */
10 
11 #include "DataFormats/Common/interface/AssociationVector.h"
12 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
13 #include "DataFormats/Candidate/interface/CandMatchMap.h"
14 #include "FWCore/Framework/interface/EDAnalyzer.h"
15 #include "DataFormats/Candidate/interface/Candidate.h"
16 #include "DataFormats/Candidate/interface/CandidateFwd.h"
17 #include "FWCore/ParameterSet/interface/ParameterSet.h"
18 #include "FWCore/Framework/interface/Event.h"
19 #include "FWCore/Framework/interface/EventSetup.h"
20 #include "FWCore/Utilities/interface/InputTag.h"
21 #include "DataFormats/Candidate/interface/OverlapChecker.h"
22 #include "DataFormats/Candidate/interface/Particle.h"
23 #include "DataFormats/Candidate/interface/Candidate.h"
24 #include "DataFormats/Common/interface/ValueMap.h"
25 #include "DataFormats/Candidate/interface/CandAssociation.h"
26 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
27 #include "TH1.h"
28 #include <vector>
29 
30 using namespace edm;
31 using namespace std;
32 using namespace reco;
33 
34 typedef ValueMap<float> IsolationCollection;
35 
36 class ZMuMuEfficiency : public edm::EDAnalyzer {
37 public:
38  ZMuMuEfficiency(const edm::ParameterSet &pset);
39 
40 private:
41  void analyze(const edm::Event &event, const edm::EventSetup &setup) override;
42  bool check_ifZmumu(const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2);
43  float getParticlePt(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2);
44  float getParticleEta(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2);
45  float getParticlePhi(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2);
46  Particle::LorentzVector getParticleP4(const int ipart,
47  const Candidate *dauGen0,
48  const Candidate *dauGen1,
49  const Candidate *dauGen2);
50  void endJob() override;
51 
52  EDGetTokenT<CandidateView> zMuMuToken_;
53  EDGetTokenT<GenParticleMatch> zMuMuMatchMapToken_;
54  EDGetTokenT<CandidateView> zMuTrackToken_;
55  EDGetTokenT<GenParticleMatch> zMuTrackMatchMapToken_;
56  EDGetTokenT<CandidateView> zMuStandAloneToken_;
57  EDGetTokenT<GenParticleMatch> zMuStandAloneMatchMapToken_;
58  EDGetTokenT<CandidateView> muonsToken_;
59  EDGetTokenT<GenParticleMatch> muonMatchMapToken_;
60  EDGetTokenT<IsolationCollection> muonIsoToken_;
61  EDGetTokenT<CandidateView> tracksToken_;
62  EDGetTokenT<IsolationCollection> trackIsoToken_;
63  EDGetTokenT<CandidateView> standAloneToken_;
64  EDGetTokenT<IsolationCollection> standAloneIsoToken_;
65  EDGetTokenT<GenParticleCollection> genParticlesToken_;
66 
67  double zMassMin_, zMassMax_, ptmin_, etamax_, isomax_;
68  unsigned int nbinsPt_, nbinsEta_;
69  reco::CandidateBaseRef globalMuonCandRef_, trackMuonCandRef_, standAloneMuonCandRef_;
70  OverlapChecker overlap_;
71 
72  //histograms for measuring tracker efficiency
73  TH1D *h_etaStandAlone_, *h_etaMuonOverlappedToStandAlone_;
74  TH1D *h_ptStandAlone_, *h_ptMuonOverlappedToStandAlone_;
75 
76  //histograms for measuring standalone efficiency
77  TH1D *h_etaTrack_, *h_etaMuonOverlappedToTrack_;
78  TH1D *h_ptTrack_, *h_ptMuonOverlappedToTrack_;
79 
80  //histograms for MC acceptance
81  TH1D *h_nZMCfound_;
82  TH1D *h_ZetaGen_, *h_ZptGen_, *h_ZmassGen_;
83  TH1D *h_muetaGen_, *h_muptGen_, *h_muIsoGen_;
84  TH1D *h_dimuonPtGen_, *h_dimuonMassGen_, *h_dimuonEtaGen_;
85  TH1D *h_ZetaGenPassed_, *h_ZptGenPassed_, *h_ZmassGenPassed_;
86  TH1D *h_muetaGenPassed_, *h_muptGenPassed_, *h_muIsoGenPassed_;
87  TH1D *h_dimuonPtGenPassed_, *h_dimuonMassGenPassed_, *h_dimuonEtaGenPassed_;
88  //histograms for invarian mass resolution
89  TH1D *h_DELTA_ZMuMuMassReco_dimuonMassGen_, *h_DELTA_ZMuStaMassReco_dimuonMassGen_,
90  *h_DELTA_ZMuTrackMassReco_dimuonMassGen_;
91 
92  int numberOfEventsWithZMuMufound, numberOfEventsWithZMuStafound;
93  int numberOfMatchedZMuSta_, numberOfMatchedSelectedZMuSta_;
94  int numberOfMatchedZMuMu_, numberOfMatchedSelectedZMuMu_;
95  int numberOfOverlappedStandAlone_, numberOfOverlappedTracks_, numberOfMatchedZMuTrack_notOverlapped;
96  int numberOfMatchedZMuTrack_exclusive, numberOfMatchedSelectedZMuTrack_exclusive;
97  int numberOfMatchedZMuTrack_matchedZMuMu, numberOfMatchedZMuTrack_matchedSelectedZMuMu;
98  int totalNumberOfevents, totalNumberOfZfound, totalNumberOfZPassed;
99  int noMCmatching, ZMuTrack_exclusive_1match, ZMuTrack_exclusive_morematch;
100  int ZMuTrackselected_exclusive_1match, ZMuTrackselected_exclusive_morematch;
101  int ZMuTrack_ZMuMu_1match, ZMuTrack_ZMuMu_2match, ZMuTrack_ZMuMu_morematch;
102 
103  int n_zMuMufound_genZsele, n_zMuStafound_genZsele, n_zMuTrkfound_genZsele;
104 };
105 
106 #include "FWCore/ServiceRegistry/interface/Service.h"
107 #include "CommonTools/UtilAlgos/interface/TFileService.h"
108 #include "DataFormats/Common/interface/Handle.h"
109 #include <iostream>
110 #include <iterator>
111 #include <cmath>
112 
113 ZMuMuEfficiency::ZMuMuEfficiency(const ParameterSet &pset)
114  : zMuMuToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuMu"))),
115  zMuMuMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuMuMatchMap"))),
116  zMuTrackToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuTrack"))),
117  zMuTrackMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuTrackMatchMap"))),
118  zMuStandAloneToken_(consumes<CandidateView>(pset.getParameter<InputTag>("zMuStandAlone"))),
119  zMuStandAloneMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("zMuStandAloneMatchMap"))),
120  muonsToken_(consumes<CandidateView>(pset.getParameter<InputTag>("muons"))),
121  muonMatchMapToken_(mayConsume<GenParticleMatch>(pset.getParameter<InputTag>("muonMatchMap"))),
122  muonIsoToken_(mayConsume<IsolationCollection>(pset.getParameter<InputTag>("muonIso"))),
123  tracksToken_(consumes<CandidateView>(pset.getParameter<InputTag>("tracks"))),
124  trackIsoToken_(mayConsume<IsolationCollection>(pset.getParameter<InputTag>("trackIso"))),
125  standAloneToken_(consumes<CandidateView>(pset.getParameter<InputTag>("standAlone"))),
126  standAloneIsoToken_(mayConsume<IsolationCollection>(pset.getParameter<InputTag>("standAloneIso"))),
127  genParticlesToken_(consumes<GenParticleCollection>(pset.getParameter<InputTag>("genParticles"))),
128 
129  zMassMin_(pset.getUntrackedParameter<double>("zMassMin")),
130  zMassMax_(pset.getUntrackedParameter<double>("zMassMax")),
131  ptmin_(pset.getUntrackedParameter<double>("ptmin")),
132  etamax_(pset.getUntrackedParameter<double>("etamax")),
133  isomax_(pset.getUntrackedParameter<double>("isomax")),
134  nbinsPt_(pset.getUntrackedParameter<unsigned int>("nbinsPt")),
135  nbinsEta_(pset.getUntrackedParameter<unsigned int>("nbinsEta")) {
136  Service<TFileService> fs;
137  TFileDirectory trackEffDir = fs->mkdir("TrackEfficiency");
138 
139  // tracker efficiency distributions
140  h_etaStandAlone_ = trackEffDir.make<TH1D>(
141  "StandAloneMuonEta", "StandAlone #eta for Z -> #mu + standalone", nbinsEta_, -etamax_, etamax_);
142  h_etaMuonOverlappedToStandAlone_ =
143  trackEffDir.make<TH1D>("MuonOverlappedToStandAloneEta",
144  "Global muon overlapped to standAlone #eta for Z -> #mu + sa",
145  nbinsEta_,
146  -etamax_,
147  etamax_);
148  h_ptStandAlone_ =
149  trackEffDir.make<TH1D>("StandAloneMuonPt", "StandAlone p_{t} for Z -> #mu + standalone", nbinsPt_, ptmin_, 100);
150  h_ptMuonOverlappedToStandAlone_ =
151  trackEffDir.make<TH1D>("MuonOverlappedToStandAlonePt",
152  "Global muon overlapped to standAlone p_{t} for Z -> #mu + sa",
153  nbinsPt_,
154  ptmin_,
155  100);
156 
157  // StandAlone efficiency distributions
158  TFileDirectory standaloneEffDir = fs->mkdir("StandaloneEfficiency");
159  h_etaTrack_ =
160  standaloneEffDir.make<TH1D>("TrackMuonEta", "Track #eta for Z -> #mu + track", nbinsEta_, -etamax_, etamax_);
161  h_etaMuonOverlappedToTrack_ = standaloneEffDir.make<TH1D>("MuonOverlappedToTrackEta",
162  "Global muon overlapped to track #eta for Z -> #mu + tk",
163  nbinsEta_,
164  -etamax_,
165  etamax_);
166  h_ptTrack_ = standaloneEffDir.make<TH1D>("TrackMuonPt", "Track p_{t} for Z -> #mu + track", nbinsPt_, ptmin_, 100);
167  h_ptMuonOverlappedToTrack_ = standaloneEffDir.make<TH1D>(
168  "MuonOverlappedToTrackPt", "Global muon overlapped to track p_{t} for Z -> #mu + tk", nbinsPt_, ptmin_, 100);
169 
170  // inv. mass resolution studies
171  TFileDirectory invMassResolutionDir = fs->mkdir("invriantMassResolution");
172  h_DELTA_ZMuMuMassReco_dimuonMassGen_ =
173  invMassResolutionDir.make<TH1D>("zMuMu_invMassResolution", "zMuMu invariant Mass Resolution", 50, -25, 25);
174  h_DELTA_ZMuStaMassReco_dimuonMassGen_ =
175  invMassResolutionDir.make<TH1D>("zMuSta_invMassResolution", "zMuSta invariant Mass Resolution", 50, -25, 25);
176  h_DELTA_ZMuTrackMassReco_dimuonMassGen_ =
177  invMassResolutionDir.make<TH1D>("zMuTrack_invMassResolution", "zMuTrack invariant Mass Resolution", 50, -25, 25);
178 
179  // generator level histograms
180  TFileDirectory genParticleDir = fs->mkdir("genParticle");
181  h_nZMCfound_ = genParticleDir.make<TH1D>("NumberOfgeneratedZeta", "n. of generated Z per event", 4, -.5, 3.5);
182  h_ZetaGen_ = genParticleDir.make<TH1D>("generatedZeta", "#eta of generated Z", 100, -5., 5.);
183  h_ZptGen_ = genParticleDir.make<TH1D>("generatedZpt", "pt of generated Z", 100, 0., 200.);
184  h_ZmassGen_ = genParticleDir.make<TH1D>("generatedZmass", "mass of generated Z", 100, 0., 200.);
185  h_muetaGen_ = genParticleDir.make<TH1D>("generatedMuonEta", "#eta of generated muons from Z decay", 100, -5., 5.);
186  h_muptGen_ = genParticleDir.make<TH1D>("generatedMuonpt", "pt of generated muons from Z decay", 100, 0., 200.);
187  h_dimuonEtaGen_ = genParticleDir.make<TH1D>("generatedDimuonEta", "#eta of generated dimuon", 100, -5., 5.);
188  h_dimuonPtGen_ = genParticleDir.make<TH1D>("generatedDimuonPt", "pt of generated dimuon", 100, 0., 200.);
189  h_dimuonMassGen_ = genParticleDir.make<TH1D>("generatedDimuonMass", "mass of generated dimuon", 100, 0., 200.);
190  h_ZetaGenPassed_ = genParticleDir.make<TH1D>("generatedZeta_passed", "#eta of generated Z after cuts", 100, -5., 5.);
191  h_ZptGenPassed_ = genParticleDir.make<TH1D>("generatedZpt_passed", "pt of generated Z after cuts", 100, 0., 200.);
192  h_ZmassGenPassed_ =
193  genParticleDir.make<TH1D>("generatedZmass_passed", "mass of generated Z after cuts", 100, 0., 200.);
194  h_muetaGenPassed_ = genParticleDir.make<TH1D>(
195  "generatedMuonEta_passed", "#eta of generated muons from Z decay after cuts", 100, -5., 5.);
196  h_muptGenPassed_ = genParticleDir.make<TH1D>(
197  "generatedMuonpt_passed", "pt of generated muons from Z decay after cuts", 100, 0., 200.);
198  h_dimuonEtaGenPassed_ =
199  genParticleDir.make<TH1D>("generatedDimuonEta_passed", "#eta of generated dimuon after cuts", 100, -5., 5.);
200  h_dimuonPtGenPassed_ =
201  genParticleDir.make<TH1D>("generatedDimuonPt_passed", "pt of generated dimuon after cuts", 100, 0., 200.);
202  h_dimuonMassGenPassed_ =
203  genParticleDir.make<TH1D>("generatedDimuonMass_passed", "mass of generated dimuon after cuts", 100, 0., 200.);
204  // to insert isolation histograms ..............
205 
206  numberOfEventsWithZMuMufound = 0;
207  numberOfEventsWithZMuStafound = 0;
208  numberOfMatchedZMuMu_ = 0;
209  numberOfMatchedSelectedZMuMu_ = 0;
210  numberOfMatchedZMuSta_ = 0;
211  numberOfMatchedSelectedZMuSta_ = 0;
212  numberOfMatchedZMuTrack_matchedZMuMu = 0;
213  numberOfMatchedZMuTrack_matchedSelectedZMuMu = 0;
214  numberOfMatchedZMuTrack_exclusive = 0;
215  numberOfMatchedSelectedZMuTrack_exclusive = 0;
216  numberOfOverlappedStandAlone_ = 0;
217  numberOfOverlappedTracks_ = 0;
218  numberOfMatchedZMuTrack_notOverlapped = 0;
219  noMCmatching = 0;
220  ZMuTrack_exclusive_1match = 0;
221  ZMuTrack_exclusive_morematch = 0;
222  ZMuTrackselected_exclusive_1match = 0;
223  ZMuTrackselected_exclusive_morematch = 0;
224  ZMuTrack_ZMuMu_1match = 0;
225  ZMuTrack_ZMuMu_2match = 0;
226  ZMuTrack_ZMuMu_morematch = 0;
227 
228  n_zMuMufound_genZsele = 0;
229  n_zMuStafound_genZsele = 0;
230  n_zMuTrkfound_genZsele = 0;
231 
232  // generator counters
233  totalNumberOfevents = 0;
234  totalNumberOfZfound = 0;
235  totalNumberOfZPassed = 0;
236 }
237 
238 void ZMuMuEfficiency::analyze(const Event &event, const EventSetup &setup) {
239  Handle<CandidateView> zMuMu;
240  Handle<GenParticleMatch> zMuMuMatchMap; //Map of Z made by Mu global + Mu global
241  Handle<CandidateView> zMuTrack;
242  Handle<GenParticleMatch> zMuTrackMatchMap; //Map of Z made by Mu + Track
243  Handle<CandidateView> zMuStandAlone;
244  Handle<GenParticleMatch> zMuStandAloneMatchMap; //Map of Z made by Mu + StandAlone
245  Handle<CandidateView> muons; //Collection of Muons
246  Handle<GenParticleMatch> muonMatchMap;
247  Handle<IsolationCollection> muonIso;
248  Handle<CandidateView> tracks; //Collection of Tracks
249  Handle<IsolationCollection> trackIso;
250  Handle<CandidateView> standAlone; //Collection of StandAlone
251  Handle<IsolationCollection> standAloneIso;
252  Handle<GenParticleCollection> genParticles; // Collection of Generatd Particles
253 
254  event.getByToken(zMuMuToken_, zMuMu);
255  event.getByToken(zMuTrackToken_, zMuTrack);
256  event.getByToken(zMuStandAloneToken_, zMuStandAlone);
257  event.getByToken(muonsToken_, muons);
258  event.getByToken(tracksToken_, tracks);
259  event.getByToken(standAloneToken_, standAlone);
260  event.getByToken(genParticlesToken_, genParticles);
261 
262  cout << "********* zMuMu size : " << zMuMu->size() << endl;
263  cout << "********* zMuStandAlone size : " << zMuStandAlone->size() << endl;
264  cout << "********* zMuTrack size : " << zMuTrack->size() << endl;
265  cout << "********* muons size : " << muons->size() << endl;
266  cout << "********* standAlone size : " << standAlone->size() << endl;
267  cout << "********* tracks size : " << tracks->size() << endl;
268  cout << "********* generated size : " << genParticles->size() << endl;
269 
270  // generator level distributions
271 
272  int nZMCfound = 0;
273  totalNumberOfevents++;
274  int ngen = genParticles->size();
275  bool ZMuMuMatchedfound = false;
276  bool ZMuMuMatchedSelectedfound = false;
277  bool ZMuStaMatchedfound = false;
278  //bool ZMuStaMatchedSelectedfound = false;
279  int ZMuTrackMatchedfound = 0;
280  int ZMuTrackMatchedSelected_exclusivefound = 0;
281 
282  double dimuonMassGen = 0;
283 
284  for (int i = 0; i < ngen; i++) {
285  const Candidate &genCand = (*genParticles)[i];
286 
287  // if((genCand.pdgId() == 23) && (genCand.status() == 2)) //this is an intermediate Z0
288  // cout << ">>> intermediate Z0 found, with " << genCand.numberOfDaughters() << " daughters" << endl;
289  if ((genCand.pdgId() == 23) && (genCand.status() == 3)) { //this is a Z0
290  if (genCand.numberOfDaughters() == 3) { // possible Z0 decays in mu+ mu-, the 3rd daughter is the same Z0
291  const Candidate *dauGen0 = genCand.daughter(0);
292  const Candidate *dauGen1 = genCand.daughter(1);
293  const Candidate *dauGen2 = genCand.daughter(2);
294  if (check_ifZmumu(dauGen0, dauGen1, dauGen2)) { // Z0 in mu+ mu-
295  totalNumberOfZfound++;
296  nZMCfound++;
297  bool checkpt = false;
298  bool checketa = false;
299  bool checkmass = false;
300  float mupluspt, muminuspt, mupluseta, muminuseta;
301  mupluspt = getParticlePt(-13, dauGen0, dauGen1, dauGen2);
302  muminuspt = getParticlePt(13, dauGen0, dauGen1, dauGen2);
303  mupluseta = getParticleEta(-13, dauGen0, dauGen1, dauGen2);
304  muminuseta = getParticleEta(13, dauGen0, dauGen1, dauGen2);
305  //float muplusphi = getParticlePhi(-13,dauGen0,dauGen1,dauGen2);
306  //float muminusphi = getParticlePhi(13,dauGen0,dauGen1,dauGen2);
307  Particle::LorentzVector pZ(0, 0, 0, 0);
308  Particle::LorentzVector muplusp4 = getParticleP4(-13, dauGen0, dauGen1, dauGen2);
309  Particle::LorentzVector muminusp4 = getParticleP4(13, dauGen0, dauGen1, dauGen2);
310  pZ = muplusp4 + muminusp4;
311  double dimuon_pt = sqrt(pZ.x() * pZ.x() + pZ.y() * pZ.y());
312  double tan_theta_half = tan(atan(dimuon_pt / pZ.z()) / 2.);
313  double dimuon_eta = 0.;
314  if (tan_theta_half > 0)
315  dimuon_eta = -log(tan(tan_theta_half));
316  if (tan_theta_half <= 0)
317  dimuon_eta = log(tan(-tan_theta_half));
318 
319  dimuonMassGen = pZ.mass(); // dimuon invariant Mass at Generator Level
320 
321  h_ZmassGen_->Fill(genCand.mass());
322  h_ZetaGen_->Fill(genCand.eta());
323  h_ZptGen_->Fill(genCand.pt());
324  h_dimuonMassGen_->Fill(pZ.mass());
325  h_dimuonEtaGen_->Fill(dimuon_eta);
326  h_dimuonPtGen_->Fill(dimuon_pt);
327  h_muetaGen_->Fill(mupluseta);
328  h_muetaGen_->Fill(muminuseta);
329  h_muptGen_->Fill(mupluspt);
330  h_muptGen_->Fill(muminuspt);
331  // dimuon 4-momentum
332  // h_mDimuonMC->Fill(pZ.mass());
333  // h_ZminusDimuonMassMC->Fill(genCand.mass()-pZ.mass());
334  // h_DeltaPhiMC->Fill(deltaPhi(muplusphi,muminusphi));
335  // if (dauGen2==23) float z_eta = dauGen2->eta();
336  // if (dauGen2==23) float Zpt = dauGen2->pt();
337  // h_DeltaPhivsZPtMC->Fill(DeltaPhi(muplusphi,muminusphi),ZPt);
338 
339  if (mupluspt > ptmin_ && muminuspt > ptmin_)
340  checkpt = true;
341  if (mupluseta < etamax_ && muminuseta < etamax_)
342  checketa = true;
343  if (genCand.mass() > zMassMin_ && genCand.mass() < zMassMax_)
344  checkmass = true;
345  if (checkpt && checketa && checkmass) {
346  totalNumberOfZPassed++;
347  h_ZmassGenPassed_->Fill(genCand.mass());
348  h_ZetaGenPassed_->Fill(genCand.eta());
349  h_ZptGenPassed_->Fill(genCand.pt());
350  h_dimuonMassGenPassed_->Fill(pZ.mass());
351  h_dimuonEtaGenPassed_->Fill(dimuon_eta);
352  h_dimuonPtGenPassed_->Fill(dimuon_pt);
353  h_muetaGenPassed_->Fill(mupluseta);
354  h_muetaGenPassed_->Fill(muminuseta);
355  h_muptGenPassed_->Fill(mupluspt);
356  h_muptGenPassed_->Fill(muminuspt);
357 
358  if (!zMuMu->empty()) {
359  n_zMuMufound_genZsele++;
360  } else if (!zMuStandAlone->empty()) {
361  n_zMuStafound_genZsele++;
362  } else {
363  n_zMuTrkfound_genZsele++;
364  }
365  }
366  }
367  }
368  }
369  }
370  h_nZMCfound_->Fill(nZMCfound); // number of Z found in the event at generator level
371 
372  //TRACK efficiency (conto numero di eventi Zmumu global e ZmuSta (ricorda che sono due campioni esclusivi)
373 
374  if (!zMuMu->empty()) {
375  numberOfEventsWithZMuMufound++;
376  event.getByToken(zMuMuMatchMapToken_, zMuMuMatchMap);
377  event.getByToken(muonIsoToken_, muonIso);
378  event.getByToken(standAloneIsoToken_, standAloneIso);
379  event.getByToken(muonMatchMapToken_, muonMatchMap);
380  for (unsigned int i = 0; i < zMuMu->size(); ++i) { //loop on candidates
381  const Candidate &zMuMuCand = (*zMuMu)[i]; //the candidate
382  CandidateBaseRef zMuMuCandRef = zMuMu->refAt(i);
383  bool isMatched = false;
384  GenParticleRef zMuMuMatch = (*zMuMuMatchMap)[zMuMuCandRef];
385 
386  if (zMuMuMatch.isNonnull()) { // ZMuMu matched
387  isMatched = true;
388  numberOfMatchedZMuMu_++;
389  }
390  CandidateBaseRef dau0 = zMuMuCand.daughter(0)->masterClone();
391  CandidateBaseRef dau1 = zMuMuCand.daughter(1)->masterClone();
392  if (isMatched)
393  ZMuMuMatchedfound = true;
394 
395  // Cuts
396  if ((dau0->pt() > ptmin_) && (dau1->pt() > ptmin_) && (fabs(dau0->eta()) < etamax_) &&
397  (fabs(dau1->eta()) < etamax_) && (zMuMuCand.mass() > zMassMin_) && (zMuMuCand.mass() < zMassMax_) &&
398  (isMatched)) {
399  //The Z daughters are already matched!
400  const double globalMuonIsolation0 = (*muonIso)[dau0];
401  const double globalMuonIsolation1 = (*muonIso)[dau1];
402  if ((globalMuonIsolation0 < isomax_) &&
403  (globalMuonIsolation1 < isomax_)) { // ZMuMu matched and selected by cuts
404  ZMuMuMatchedSelectedfound = true;
405  numberOfMatchedSelectedZMuMu_++;
406  h_etaStandAlone_->Fill(dau0->eta()); // StandAlone found dau0, eta
407  h_etaStandAlone_->Fill(dau1->eta()); // StandAlone found dau1, eta
408  h_etaMuonOverlappedToStandAlone_->Fill(dau0->eta()); // is global muon so dau0 is also found as a track, eta
409  h_etaMuonOverlappedToStandAlone_->Fill(dau1->eta()); // is global muon so dau1 is also found as a track, eta
410  h_ptStandAlone_->Fill(dau0->pt()); // StandAlone found dau0, pt
411  h_ptStandAlone_->Fill(dau1->pt()); // StandAlone found dau1, pt
412  h_ptMuonOverlappedToStandAlone_->Fill(dau0->pt()); // is global muon so dau0 is also found as a track, pt
413  h_ptMuonOverlappedToStandAlone_->Fill(dau1->pt()); // is global muon so dau1 is also found as a track, pt
414 
415  h_etaTrack_->Fill(dau0->eta()); // Track found dau0, eta
416  h_etaTrack_->Fill(dau1->eta()); // Track found dau1, eta
417  h_etaMuonOverlappedToTrack_->Fill(dau0->eta()); // is global muon so dau0 is also found as a StandAlone, eta
418  h_etaMuonOverlappedToTrack_->Fill(dau1->eta()); // is global muon so dau1 is also found as a StandAlone, eta
419  h_ptTrack_->Fill(dau0->pt()); // Track found dau0, pt
420  h_ptTrack_->Fill(dau1->pt()); // Track found dau1, pt
421  h_ptMuonOverlappedToTrack_->Fill(dau0->pt()); // is global muon so dau0 is also found as a StandAlone, pt
422  h_ptMuonOverlappedToTrack_->Fill(dau1->pt()); // is global muon so dau1 is also found as a StandAlone, pt
423 
424  h_DELTA_ZMuMuMassReco_dimuonMassGen_->Fill(zMuMuCand.mass() - dimuonMassGen);
425  // check that the two muons are matched . .per ora è solo un mio controllo
426  for (unsigned int j = 0; j < muons->size(); ++j) {
427  CandidateBaseRef muCandRef = muons->refAt(j);
428  GenParticleRef muonMatch = (*muonMatchMap)[muCandRef];
429  // if (muonMatch.isNonnull()) cout << "mu match n. " << j << endl;
430  }
431  }
432  }
433  }
434  }
435 
436  if (!zMuStandAlone->empty()) {
437  numberOfEventsWithZMuStafound++;
438  event.getByToken(zMuStandAloneMatchMapToken_, zMuStandAloneMatchMap);
439  event.getByToken(muonIsoToken_, muonIso);
440  event.getByToken(standAloneIsoToken_, standAloneIso);
441  event.getByToken(muonMatchMapToken_, muonMatchMap);
442  for (unsigned int i = 0; i < zMuStandAlone->size(); ++i) { //loop on candidates
443  const Candidate &zMuStaCand = (*zMuStandAlone)[i]; //the candidate
444  CandidateBaseRef zMuStaCandRef = zMuStandAlone->refAt(i);
445  bool isMatched = false;
446  GenParticleRef zMuStaMatch = (*zMuStandAloneMatchMap)[zMuStaCandRef];
447  if (zMuStaMatch.isNonnull()) { // ZMuSta Macthed
448  isMatched = true;
449  ZMuStaMatchedfound = true;
450  numberOfMatchedZMuSta_++;
451  }
452  CandidateBaseRef dau0 = zMuStaCand.daughter(0)->masterClone();
453  CandidateBaseRef dau1 = zMuStaCand.daughter(1)->masterClone();
454 
455  // Cuts
456  if ((dau0->pt() > ptmin_) && (dau1->pt() > ptmin_) && (fabs(dau0->eta()) < etamax_) &&
457  (fabs(dau1->eta()) < etamax_) && (zMuStaCand.mass() > zMassMin_) && (zMuStaCand.mass() < zMassMax_) &&
458  (isMatched)) {
459  CandidateBaseRef standAloneMuonCandRef_, globalMuonCandRef_;
460  if (dau0->isGlobalMuon()) {
461  standAloneMuonCandRef_ = dau1;
462  globalMuonCandRef_ = dau0;
463  }
464  if (dau1->isGlobalMuon()) {
465  standAloneMuonCandRef_ = dau0;
466  globalMuonCandRef_ = dau1;
467  }
468 
469  const double globalMuonIsolation = (*muonIso)[globalMuonCandRef_];
470  const double standAloneMuonIsolation = (*standAloneIso)[standAloneMuonCandRef_];
471 
472  if ((globalMuonIsolation < isomax_) && (standAloneMuonIsolation < isomax_)) { // ZMuSta matched and selected
473  //ZMuStaMatchedSelectedfound = true;
474  numberOfMatchedSelectedZMuSta_++;
475  h_etaStandAlone_->Fill(standAloneMuonCandRef_->eta()); //Denominator eta for measuring track efficiency
476  h_ptStandAlone_->Fill(standAloneMuonCandRef_->pt()); //Denominator pt for measuring track eff
477  h_DELTA_ZMuStaMassReco_dimuonMassGen_->Fill(
478  zMuStaCand.mass() - dimuonMassGen); // differnce between ZMuSta reco and dimuon mass gen
479  }
480  }
481  }
482  } //end loop on Candidate
483 
484  //STANDALONE efficiency
485 
486  if (!zMuTrack->empty()) {
487  event.getByToken(zMuTrackMatchMapToken_, zMuTrackMatchMap);
488  event.getByToken(muonIsoToken_, muonIso);
489  event.getByToken(trackIsoToken_, trackIso);
490  event.getByToken(muonMatchMapToken_, muonMatchMap);
491  for (unsigned int i = 0; i < zMuTrack->size(); ++i) { //loop on candidates
492  const Candidate &zMuTrkCand = (*zMuTrack)[i]; //the candidate
493  CandidateBaseRef zMuTrkCandRef = zMuTrack->refAt(i);
494  bool isMatched = false;
495  GenParticleRef zMuTrkMatch = (*zMuTrackMatchMap)[zMuTrkCandRef];
496  if (zMuTrkMatch.isNonnull()) {
497  isMatched = true;
498  }
499  CandidateBaseRef dau0 = zMuTrkCand.daughter(0)->masterClone();
500  CandidateBaseRef dau1 = zMuTrkCand.daughter(1)->masterClone();
501 
502  if (isMatched) {
503  ZMuTrackMatchedfound++;
504  if (ZMuMuMatchedfound)
505  numberOfMatchedZMuTrack_matchedZMuMu++;
506  if (ZMuMuMatchedSelectedfound)
507  numberOfMatchedZMuTrack_matchedSelectedZMuMu++;
508  if (!ZMuMuMatchedfound)
509  numberOfMatchedZMuTrack_exclusive++;
510  }
511  // Cuts
512  if ((dau0->pt() > ptmin_) && (dau1->pt() > ptmin_) && (fabs(dau0->eta()) < etamax_) &&
513  (fabs(dau1->eta()) < etamax_) && (zMuTrkCand.mass() > zMassMin_) && (zMuTrkCand.mass() < zMassMax_) &&
514  (isMatched) && !ZMuMuMatchedfound && !ZMuStaMatchedfound) {
515  // dau0 is always the global muon, dau1 is the track for ZMuTrack collection
516  const double globalMuonIsolation = (*muonIso)[dau0];
517  const double trackMuonIsolation = (*trackIso)[dau1];
518  if ((globalMuonIsolation < isomax_) &&
519  (trackMuonIsolation < isomax_)) { // ZMuTRack matched - selected without ZMuMu found (exclusive)
520  numberOfMatchedSelectedZMuTrack_exclusive++;
521  ZMuTrackMatchedSelected_exclusivefound++;
522  h_etaTrack_->Fill(dau1->eta()); //Denominator eta Sta
523  h_ptTrack_->Fill(dau1->pt()); //Denominator pt Sta
524  h_DELTA_ZMuTrackMassReco_dimuonMassGen_->Fill(zMuTrkCand.mass() - dimuonMassGen);
525  }
526  }
527  }
528  } //end loop on Candidate
529 
530  if (!ZMuMuMatchedfound && !ZMuStaMatchedfound && ZMuTrackMatchedfound == 0)
531  noMCmatching++;
532  if (!ZMuMuMatchedfound && ZMuTrackMatchedfound == 1)
533  ZMuTrack_exclusive_1match++;
534  if (!ZMuMuMatchedfound && ZMuTrackMatchedfound > 1)
535  ZMuTrack_exclusive_morematch++;
536  if (!ZMuMuMatchedfound && ZMuTrackMatchedSelected_exclusivefound == 1)
537  ZMuTrackselected_exclusive_1match++;
538  if (!ZMuMuMatchedfound && ZMuTrackMatchedSelected_exclusivefound > 1)
539  ZMuTrackselected_exclusive_morematch++;
540  if (ZMuMuMatchedfound && ZMuTrackMatchedfound == 1)
541  ZMuTrack_ZMuMu_1match++;
542  if (ZMuMuMatchedfound && ZMuTrackMatchedfound == 2)
543  ZMuTrack_ZMuMu_2match++;
544  if (ZMuMuMatchedfound && ZMuTrackMatchedfound > 2)
545  ZMuTrack_ZMuMu_morematch++;
546 }
547 
548 bool ZMuMuEfficiency::check_ifZmumu(const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2) {
549  int partId0 = dauGen0->pdgId();
550  int partId1 = dauGen1->pdgId();
551  int partId2 = dauGen2->pdgId();
552  bool muplusFound = false;
553  bool muminusFound = false;
554  bool ZFound = false;
555  if (partId0 == 13 || partId1 == 13 || partId2 == 13)
556  muminusFound = true;
557  if (partId0 == -13 || partId1 == -13 || partId2 == -13)
558  muplusFound = true;
559  if (partId0 == 23 || partId1 == 23 || partId2 == 23)
560  ZFound = true;
561  return (muplusFound && muminusFound && ZFound);
562 }
563 
564 float ZMuMuEfficiency::getParticlePt(const int ipart,
565  const Candidate *dauGen0,
566  const Candidate *dauGen1,
567  const Candidate *dauGen2) {
568  int partId0 = dauGen0->pdgId();
569  int partId1 = dauGen1->pdgId();
570  int partId2 = dauGen2->pdgId();
571  float ptpart = 0.;
572  if (partId0 == ipart) {
573  for (unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
574  const Candidate *dauMuGen = dauGen0->daughter(k);
575  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
576  ptpart = dauMuGen->pt();
577  }
578  }
579  }
580  if (partId1 == ipart) {
581  for (unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
582  const Candidate *dauMuGen = dauGen1->daughter(k);
583  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
584  ptpart = dauMuGen->pt();
585  }
586  }
587  }
588  if (partId2 == ipart) {
589  for (unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
590  const Candidate *dauMuGen = dauGen2->daughter(k);
591  if (abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() == 1) {
592  ptpart = dauMuGen->pt();
593  }
594  }
595  }
596  return ptpart;
597 }
598 
599 float ZMuMuEfficiency::getParticleEta(const int ipart,
600  const Candidate *dauGen0,
601  const Candidate *dauGen1,
602  const Candidate *dauGen2) {
603  int partId0 = dauGen0->pdgId();
604  int partId1 = dauGen1->pdgId();
605  int partId2 = dauGen2->pdgId();
606  float etapart = 0.;
607  if (partId0 == ipart) {
608  for (unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
609  const Candidate *dauMuGen = dauGen0->daughter(k);
610  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
611  etapart = dauMuGen->eta();
612  }
613  }
614  }
615  if (partId1 == ipart) {
616  for (unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
617  const Candidate *dauMuGen = dauGen1->daughter(k);
618  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
619  etapart = dauMuGen->eta();
620  }
621  }
622  }
623  if (partId2 == ipart) {
624  for (unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
625  const Candidate *dauMuGen = dauGen2->daughter(k);
626  if (abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() == 1) {
627  etapart = dauMuGen->eta();
628  }
629  }
630  }
631  return etapart;
632 }
633 
634 float ZMuMuEfficiency::getParticlePhi(const int ipart,
635  const Candidate *dauGen0,
636  const Candidate *dauGen1,
637  const Candidate *dauGen2) {
638  int partId0 = dauGen0->pdgId();
639  int partId1 = dauGen1->pdgId();
640  int partId2 = dauGen2->pdgId();
641  float phipart = 0.;
642  if (partId0 == ipart) {
643  for (unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
644  const Candidate *dauMuGen = dauGen0->daughter(k);
645  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
646  phipart = dauMuGen->phi();
647  }
648  }
649  }
650  if (partId1 == ipart) {
651  for (unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
652  const Candidate *dauMuGen = dauGen1->daughter(k);
653  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
654  phipart = dauMuGen->phi();
655  }
656  }
657  }
658  if (partId2 == ipart) {
659  for (unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
660  const Candidate *dauMuGen = dauGen2->daughter(k);
661  if (abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() == 1) {
662  phipart = dauMuGen->phi();
663  }
664  }
665  }
666  return phipart;
667 }
668 
669 Particle::LorentzVector ZMuMuEfficiency::getParticleP4(const int ipart,
670  const Candidate *dauGen0,
671  const Candidate *dauGen1,
672  const Candidate *dauGen2) {
673  int partId0 = dauGen0->pdgId();
674  int partId1 = dauGen1->pdgId();
675  int partId2 = dauGen2->pdgId();
676  Particle::LorentzVector p4part(0., 0., 0., 0.);
677  if (partId0 == ipart) {
678  for (unsigned int k = 0; k < dauGen0->numberOfDaughters(); ++k) {
679  const Candidate *dauMuGen = dauGen0->daughter(k);
680  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
681  p4part = dauMuGen->p4();
682  }
683  }
684  }
685  if (partId1 == ipart) {
686  for (unsigned int k = 0; k < dauGen1->numberOfDaughters(); ++k) {
687  const Candidate *dauMuGen = dauGen1->daughter(k);
688  if (dauMuGen->pdgId() == ipart && dauMuGen->status() == 1) {
689  p4part = dauMuGen->p4();
690  }
691  }
692  }
693  if (partId2 == ipart) {
694  for (unsigned int k = 0; k < dauGen2->numberOfDaughters(); ++k) {
695  const Candidate *dauMuGen = dauGen2->daughter(k);
696  if (abs(dauMuGen->pdgId()) == ipart && dauMuGen->status() == 1) {
697  p4part = dauMuGen->p4();
698  }
699  }
700  }
701  return p4part;
702 }
703 
704 void ZMuMuEfficiency::endJob() {
705  // double efficiencySTA =(double)numberOfOverlappedStandAlone_/(double)numberOfMatchedZMuTrack_;
706  // double errorEff_STA = sqrt( efficiencySTA*(1 - efficiencySTA)/(double)numberOfMatchedZMuTrack_);
707 
708  double myTrackEff = 2. * numberOfMatchedSelectedZMuMu_ /
709  (2. * numberOfMatchedSelectedZMuMu_ + (double)numberOfMatchedSelectedZMuSta_);
710  double myErrTrackEff = sqrt(myTrackEff * (1 - myTrackEff) /
711  (2. * numberOfMatchedSelectedZMuMu_ + (double)numberOfMatchedSelectedZMuSta_));
712 
713  double myStaEff = 2. * numberOfMatchedSelectedZMuMu_ /
714  (2. * numberOfMatchedSelectedZMuMu_ + (double)numberOfMatchedSelectedZMuTrack_exclusive);
715  double myErrStaEff = sqrt(myTrackEff * (1 - myTrackEff) /
716  (2. * numberOfMatchedSelectedZMuMu_ + (double)numberOfMatchedSelectedZMuTrack_exclusive));
717 
718  // double efficiencyTRACK =(double)numberOfOverlappedTracks_/(double)numberOfMatchedZMuSta_;
719  // double errorEff_TRACK = sqrt( efficiencyTRACK*(1 - efficiencyTRACK)/(double)numberOfMatchedZMuSta_);
720 
721  cout << "------------------------------------ Counters for MC acceptance --------------------------------" << endl;
722  cout << "totalNumberOfevents = " << totalNumberOfevents << endl;
723  cout << "totalNumberOfZfound = " << totalNumberOfZfound << endl;
724  cout << "totalNumberOfZpassed = " << totalNumberOfZPassed << endl;
725  cout << "n. of events zMuMu found (gen level selected)" << n_zMuMufound_genZsele << endl;
726  cout << "n. of events zMuSta found (gen level selected)" << n_zMuStafound_genZsele << endl;
727  cout << "n. of events zMuTrk found (gen level selected)" << n_zMuTrkfound_genZsele << endl;
728 
729  cout << "---------------------------- Counter for MC truth efficiency calculation--------------------- " << endl;
730 
731  cout << "number of events with ZMuMu found = " << numberOfEventsWithZMuMufound << endl;
732  cout << "number of events with ZMuSta found = " << numberOfEventsWithZMuStafound << endl;
733  cout << "-------------------------------------------------------------------------------------- " << endl;
734 
735  cout << "number of events without MC maching = " << noMCmatching << endl;
736  cout << "number of ZMuTrack exclsive 1 match = " << ZMuTrack_exclusive_1match << endl;
737  cout << "number of ZMuTrack exclsive more match = " << ZMuTrack_exclusive_morematch << endl;
738  cout << "number of ZMuTrack selected exclusive 1 match = " << ZMuTrackselected_exclusive_1match << endl;
739  cout << "number of ZMuTrack selected exclusive more match = " << ZMuTrackselected_exclusive_morematch << endl;
740  cout << "number of ZMuTrack ZMuMu 1 match = " << ZMuTrack_ZMuMu_1match << endl;
741  cout << "number of ZMuTrack ZMuMu 2 match = " << ZMuTrack_ZMuMu_2match << endl;
742  cout << "number of ZMuTrack ZMuMu more match = " << ZMuTrack_ZMuMu_morematch << endl;
743  cout << "numberOfMatchedZMuMu = " << numberOfMatchedZMuMu_ << endl;
744  cout << "numberOfMatchedSelectdZMuMu = " << numberOfMatchedSelectedZMuMu_ << endl;
745  cout << "numberOfMatchedZMuSta = " << numberOfMatchedZMuSta_ << endl;
746  cout << "numberOfMatchedSelectedZMuSta = " << numberOfMatchedSelectedZMuSta_ << endl;
747  cout << "numberOfMatchedZMuTrack_matchedZMuMu = " << numberOfMatchedZMuTrack_matchedZMuMu << endl;
748  cout << "numberOfMatchedZMuTrack_matchedSelectedZMuMu = " << numberOfMatchedZMuTrack_matchedSelectedZMuMu << endl;
749  cout << "numberOfMatchedZMuTrack exclusive = " << numberOfMatchedZMuTrack_exclusive << endl;
750  cout << "numberOfMatchedSelectedZMuTrack exclusive = " << numberOfMatchedSelectedZMuTrack_exclusive << endl;
751  cout << " ----------------------------- Efficiency --------------------------------- " << endl;
752  cout << "Efficiency StandAlone = " << myStaEff << " +/- " << myErrStaEff << endl;
753  cout << "Efficiency Track = " << myTrackEff << " +/- " << myErrTrackEff << endl;
754 }
755 
756 #include "FWCore/Framework/interface/MakerMacros.h"
757 
758 DEFINE_FWK_MODULE(ZMuMuEfficiency);
ZMuMuEfficiency::ZMuTrack_exclusive_morematch
int ZMuTrack_exclusive_morematch
Definition: ZMuMuEfficiency.cc:99
ZMuMuEfficiency::totalNumberOfZPassed
int totalNumberOfZPassed
Definition: ZMuMuEfficiency.cc:98
ZMuMuEfficiency::numberOfEventsWithZMuMufound
int numberOfEventsWithZMuMufound
Definition: ZMuMuEfficiency.cc:92
ZMuMuEfficiency::h_DELTA_ZMuMuMassReco_dimuonMassGen_
TH1D * h_DELTA_ZMuMuMassReco_dimuonMassGen_
Definition: ZMuMuEfficiency.cc:89
reco::Candidate::daughter
virtual const Candidate * daughter(size_type i) const =0
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode)
ZMuMuEfficiency::nbinsEta_
unsigned int nbinsEta_
Definition: ZMuMuEfficiency.cc:68
ZMuMuEfficiency::numberOfMatchedZMuTrack_matchedZMuMu
int numberOfMatchedZMuTrack_matchedZMuMu
Definition: ZMuMuEfficiency.cc:97
ZMuMuEfficiency::numberOfMatchedZMuTrack_notOverlapped
int numberOfMatchedZMuTrack_notOverlapped
Definition: ZMuMuEfficiency.cc:95
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
OverlapChecker
Definition: OverlapChecker.h:17
Handle.h
ZMuMuEfficiency::numberOfMatchedSelectedZMuMu_
int numberOfMatchedSelectedZMuMu_
Definition: ZMuMuEfficiency.cc:94
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
ZMuMuEfficiency::trackIsoToken_
EDGetTokenT< IsolationCollection > trackIsoToken_
Definition: ZMuMuEfficiency.cc:62
mps_fire.i
i
Definition: mps_fire.py:355
ZMuMuEfficiency::h_dimuonEtaGen_
TH1D * h_dimuonEtaGen_
Definition: ZMuMuEfficiency.cc:84
egammaDQMOffline_cff.standAlone
standAlone
Definition: egammaDQMOffline_cff.py:13
ZMuMuEfficiency::muonsToken_
EDGetTokenT< CandidateView > muonsToken_
Definition: ZMuMuEfficiency.cc:58
genParticles2HepMC_cfi.genParticles
genParticles
Definition: genParticles2HepMC_cfi.py:4
ZMuMuAnalysisNtupler_cff.zMuMu
zMuMu
zMuMu vector of PSet is common to all categories except zMuTrk category
Definition: ZMuMuAnalysisNtupler_cff.py:340
edm::View::empty
bool empty() const
ZMuMuEfficiency::numberOfMatchedZMuTrack_exclusive
int numberOfMatchedZMuTrack_exclusive
Definition: ZMuMuEfficiency.cc:96
ZMuMuEfficiency::ptmin_
double ptmin_
Definition: ZMuMuEfficiency.cc:67
reco::Candidate::mass
virtual double mass() const =0
mass
ZMuMuEfficiency::h_dimuonEtaGenPassed_
TH1D * h_dimuonEtaGenPassed_
Definition: ZMuMuEfficiency.cc:87
edm::View::refAt
RefToBase< value_type > refAt(size_type i) const
reco::Candidate::eta
virtual double eta() const =0
momentum pseudorapidity
edm::EDGetTokenT< CandidateView >
TFileDirectory::make
T * make(const Args &... args) const
make new ROOT object
Definition: TFileDirectory.h:53
edm
HLT enums.
Definition: AlignableModifier.h:19
gather_cfg.cout
cout
Definition: gather_cfg.py:144
reco::GenParticleCollection
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
Definition: GenParticleFwd.h:13
ZMuMuEfficiency::h_muptGenPassed_
TH1D * h_muptGenPassed_
Definition: ZMuMuEfficiency.cc:86
ZMuMuEfficiency::standAloneMuonCandRef_
reco::CandidateBaseRef standAloneMuonCandRef_
Definition: ZMuMuEfficiency.cc:69
ZMuMuEfficiency::endJob
void endJob() override
Definition: ZMuMuEfficiency.cc:704
ZMuMuEfficiency::h_ZmassGenPassed_
TH1D * h_ZmassGenPassed_
Definition: ZMuMuEfficiency.cc:85
ZMuMuEfficiency::numberOfOverlappedTracks_
int numberOfOverlappedTracks_
Definition: ZMuMuEfficiency.cc:95
ZMuMuEfficiency::n_zMuMufound_genZsele
int n_zMuMufound_genZsele
Definition: ZMuMuEfficiency.cc:103
ZMuMuEfficiency::ZMuMuEfficiency
ZMuMuEfficiency(const edm::ParameterSet &pset)
Definition: ZMuMuEfficiency.cc:113
reco::Candidate::pt
virtual double pt() const =0
transverse momentum
ZMuMuEfficiency::totalNumberOfZfound
int totalNumberOfZfound
Definition: ZMuMuEfficiency.cc:98
ZMuMuEfficiency::ZMuTrackselected_exclusive_morematch
int ZMuTrackselected_exclusive_morematch
Definition: ZMuMuEfficiency.cc:100
ZMuMuEfficiency::check_ifZmumu
bool check_ifZmumu(const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Definition: ZMuMuEfficiency.cc:548
ZMuMuEfficiency::numberOfMatchedSelectedZMuTrack_exclusive
int numberOfMatchedSelectedZMuTrack_exclusive
Definition: ZMuMuEfficiency.cc:96
ZMuMuEfficiency::ZMuTrack_ZMuMu_2match
int ZMuTrack_ZMuMu_2match
Definition: ZMuMuEfficiency.cc:101
reco::Candidate::status
virtual int status() const =0
status word
ZMuMuEfficiency::numberOfMatchedSelectedZMuSta_
int numberOfMatchedSelectedZMuSta_
Definition: ZMuMuEfficiency.cc:93
ZMuMuEfficiency::ZMuTrack_ZMuMu_1match
int ZMuTrack_ZMuMu_1match
Definition: ZMuMuEfficiency.cc:101
ZMuMuEfficiency::zMassMin_
double zMassMin_
Definition: ZMuMuEfficiency.cc:67
EDAnalyzer.h
TFileDirectory
Definition: TFileDirectory.h:24
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
ZMuMuEfficiency::ZMuTrack_ZMuMu_morematch
int ZMuTrack_ZMuMu_morematch
Definition: ZMuMuEfficiency.cc:101
edm::Handle
Definition: AssociativeIterator.h:50
ZMuMuEfficiency::isomax_
double isomax_
Definition: ZMuMuEfficiency.cc:67
ZMuMuEfficiency::h_etaMuonOverlappedToStandAlone_
TH1D * h_etaMuonOverlappedToStandAlone_
Definition: ZMuMuEfficiency.cc:73
singleTopDQM_cfi.setup
setup
Definition: singleTopDQM_cfi.py:37
ZMuMuEfficiency::globalMuonCandRef_
reco::CandidateBaseRef globalMuonCandRef_
Definition: ZMuMuEfficiency.cc:69
ZMuMuEfficiency::overlap_
OverlapChecker overlap_
Definition: ZMuMuEfficiency.cc:70
CandMatchMap.h
edm::Ref< GenParticleCollection >
edm::EDAnalyzer
Definition: EDAnalyzer.h:29
muonMatch_cfi.muonMatch
muonMatch
Definition: muonMatch_cfi.py:8
ZMuMuEfficiency::h_DELTA_ZMuStaMassReco_dimuonMassGen_
TH1D * h_DELTA_ZMuStaMassReco_dimuonMassGen_
Definition: ZMuMuEfficiency.cc:89
GenParticle.h
CandidateFwd.h
MakerMacros.h
ZMuMuEfficiency::h_nZMCfound_
TH1D * h_nZMCfound_
Definition: ZMuMuEfficiency.cc:81
ZMuMuEfficiency::h_dimuonPtGen_
TH1D * h_dimuonPtGen_
Definition: ZMuMuEfficiency.cc:84
ZMuMuEfficiency::numberOfMatchedZMuMu_
int numberOfMatchedZMuMu_
Definition: ZMuMuEfficiency.cc:94
ZMuMuEfficiency::etamax_
double etamax_
Definition: ZMuMuEfficiency.cc:67
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
ZMuMuEfficiency::noMCmatching
int noMCmatching
Definition: ZMuMuEfficiency.cc:99
ZMuMuEfficiency::h_ptMuonOverlappedToTrack_
TH1D * h_ptMuonOverlappedToTrack_
Definition: ZMuMuEfficiency.cc:78
Service.h
ZMuMuEfficiency::nbinsPt_
unsigned int nbinsPt_
Definition: ZMuMuEfficiency.cc:68
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
ZMuMuEfficiency::numberOfMatchedZMuTrack_matchedSelectedZMuMu
int numberOfMatchedZMuTrack_matchedSelectedZMuMu
Definition: ZMuMuEfficiency.cc:97
ZMuMuEfficiency::zMuStandAloneToken_
EDGetTokenT< CandidateView > zMuStandAloneToken_
Definition: ZMuMuEfficiency.cc:56
ZMuMuEfficiency::h_dimuonMassGenPassed_
TH1D * h_dimuonMassGenPassed_
Definition: ZMuMuEfficiency.cc:87
Particle.h
reco::Candidate::numberOfDaughters
virtual size_type numberOfDaughters() const =0
number of daughters
dqmdumpme.k
k
Definition: dqmdumpme.py:60
ZMuMuEfficiency::h_etaTrack_
TH1D * h_etaTrack_
Definition: ZMuMuEfficiency.cc:77
ZMuMuEfficiency::getParticlePhi
float getParticlePhi(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Definition: ZMuMuEfficiency.cc:634
edm::View::size
size_type size() const
TFileService.h
ZMuMuEfficiency::zMassMax_
double zMassMax_
Definition: ZMuMuEfficiency.cc:67
IsolationCollection
ValueMap< float > IsolationCollection
Definition: ZMuMuEfficiency.cc:34
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
edm::View
Definition: CaloClusterFwd.h:14
trackerHitRTTI::isMatched
bool isMatched(TrackingRecHit const &hit)
Definition: trackerHitRTTI.h:32
ZMuMuEfficiency::zMuTrackMatchMapToken_
EDGetTokenT< GenParticleMatch > zMuTrackMatchMapToken_
Definition: ZMuMuEfficiency.cc:55
ZMuMuEfficiency::totalNumberOfevents
int totalNumberOfevents
Definition: ZMuMuEfficiency.cc:98
edm::ParameterSet
Definition: ParameterSet.h:36
CandAssociation.h
Event.h
ZMuMuEfficiency::numberOfEventsWithZMuStafound
int numberOfEventsWithZMuStafound
Definition: ZMuMuEfficiency.cc:92
ZMuMuEfficiency::h_ZetaGenPassed_
TH1D * h_ZetaGenPassed_
Definition: ZMuMuEfficiency.cc:85
ZMuMuEfficiency::analyze
void analyze(const edm::Event &event, const edm::EventSetup &setup) override
Definition: ZMuMuEfficiency.cc:238
ZMuMuEfficiency::h_etaMuonOverlappedToTrack_
TH1D * h_etaMuonOverlappedToTrack_
Definition: ZMuMuEfficiency.cc:77
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
ZMuMuEfficiency::n_zMuStafound_genZsele
int n_zMuStafound_genZsele
Definition: ZMuMuEfficiency.cc:103
ZMuMuEfficiency::muonMatchMapToken_
EDGetTokenT< GenParticleMatch > muonMatchMapToken_
Definition: ZMuMuEfficiency.cc:59
ZMuMuEfficiency::h_muetaGen_
TH1D * h_muetaGen_
Definition: ZMuMuEfficiency.cc:83
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
edm::Service< TFileService >
createfilelist.int
int
Definition: createfilelist.py:10
ZMuMuEfficiency::standAloneToken_
EDGetTokenT< CandidateView > standAloneToken_
Definition: ZMuMuEfficiency.cc:63
analyze
example_stream void analyze(const edm::Event &, const edm::EventSetup &) override
ZMuMuEfficiency::h_ptStandAlone_
TH1D * h_ptStandAlone_
Definition: ZMuMuEfficiency.cc:74
edm::Association
Definition: Association.h:18
edm::EventSetup
Definition: EventSetup.h:57
ZMuMuEfficiency::getParticleEta
float getParticleEta(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Definition: ZMuMuEfficiency.cc:599
ZMuMuEfficiency::zMuStandAloneMatchMapToken_
EDGetTokenT< GenParticleMatch > zMuStandAloneMatchMapToken_
Definition: ZMuMuEfficiency.cc:57
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
ZMuMuEfficiency::tracksToken_
EDGetTokenT< CandidateView > tracksToken_
Definition: ZMuMuEfficiency.cc:61
InputTag.h
ZMuMuEfficiency::h_muptGen_
TH1D * h_muptGen_
Definition: ZMuMuEfficiency.cc:83
ZMuMuEfficiency::trackMuonCandRef_
reco::CandidateBaseRef trackMuonCandRef_
Definition: ZMuMuEfficiency.cc:69
reco::Candidate
Definition: Candidate.h:27
ZMuMuEfficiency::h_etaStandAlone_
TH1D * h_etaStandAlone_
Definition: ZMuMuEfficiency.cc:73
ValueMap.h
ZMuMuEfficiency::zMuMuToken_
EDGetTokenT< CandidateView > zMuMuToken_
Definition: ZMuMuEfficiency.cc:52
reco::JetExtendedAssociation::LorentzVector
math::PtEtaPhiELorentzVectorF LorentzVector
Definition: JetExtendedAssociation.h:25
ZMuMuEfficiency::genParticlesToken_
EDGetTokenT< GenParticleCollection > genParticlesToken_
Definition: ZMuMuEfficiency.cc:65
ZMuMuEfficiency::getParticlePt
float getParticlePt(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Definition: ZMuMuEfficiency.cc:564
ZMuMuEfficiency
Definition: ZMuMuEfficiency.cc:36
std
Definition: JetResolutionObject.h:76
ZMuMuEfficiency::h_ptTrack_
TH1D * h_ptTrack_
Definition: ZMuMuEfficiency.cc:78
ZMuMuEfficiency::h_muetaGenPassed_
TH1D * h_muetaGenPassed_
Definition: ZMuMuEfficiency.cc:86
ZMuMuEfficiency::n_zMuTrkfound_genZsele
int n_zMuTrkfound_genZsele
Definition: ZMuMuEfficiency.cc:103
ZMuMuEfficiency::h_dimuonPtGenPassed_
TH1D * h_dimuonPtGenPassed_
Definition: ZMuMuEfficiency.cc:87
ZMuMuEfficiency::h_ZetaGen_
TH1D * h_ZetaGen_
Definition: ZMuMuEfficiency.cc:82
edm::ValueMap< float >
ZMuMuEfficiency::h_ZmassGen_
TH1D * h_ZmassGen_
Definition: ZMuMuEfficiency.cc:82
ZMuMuEfficiency::ZMuTrack_exclusive_1match
int ZMuTrack_exclusive_1match
Definition: ZMuMuEfficiency.cc:99
ZMuMuEfficiency::muonIsoToken_
EDGetTokenT< IsolationCollection > muonIsoToken_
Definition: ZMuMuEfficiency.cc:60
EventSetup.h
edm::RefToBase< Candidate >
ZMuMuEfficiency::getParticleP4
Particle::LorentzVector getParticleP4(const int ipart, const Candidate *dauGen0, const Candidate *dauGen1, const Candidate *dauGen2)
Definition: ZMuMuEfficiency.cc:669
reco::Candidate::masterClone
virtual const CandidateBaseRef & masterClone() const =0
reco::Candidate::p4
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
ZMuMuEfficiency::numberOfMatchedZMuSta_
int numberOfMatchedZMuSta_
Definition: ZMuMuEfficiency.cc:93
ZMuMuEfficiency::h_ptMuonOverlappedToStandAlone_
TH1D * h_ptMuonOverlappedToStandAlone_
Definition: ZMuMuEfficiency.cc:74
AssociationVector.h
ZMuMuEfficiency::ZMuTrackselected_exclusive_1match
int ZMuTrackselected_exclusive_1match
Definition: ZMuMuEfficiency.cc:100
OverlapChecker.h
Candidate.h
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ParameterSet.h
ZMuMuEfficiency::h_ZptGen_
TH1D * h_ZptGen_
Definition: ZMuMuEfficiency.cc:82
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
event
Definition: event.py:1
edm::Event
Definition: Event.h:73
ZMuMuEfficiency::h_dimuonMassGen_
TH1D * h_dimuonMassGen_
Definition: ZMuMuEfficiency.cc:84
reco::Candidate::phi
virtual double phi() const =0
momentum azimuthal angle
ZMuMuEfficiency::zMuMuMatchMapToken_
EDGetTokenT< GenParticleMatch > zMuMuMatchMapToken_
Definition: ZMuMuEfficiency.cc:53
ZMuMuEfficiency::h_DELTA_ZMuTrackMassReco_dimuonMassGen_
TH1D * h_DELTA_ZMuTrackMassReco_dimuonMassGen_
Definition: ZMuMuEfficiency.cc:89
ZMuMuEfficiency::zMuTrackToken_
EDGetTokenT< CandidateView > zMuTrackToken_
Definition: ZMuMuEfficiency.cc:54
reco::Candidate::isGlobalMuon
virtual bool isGlobalMuon() const =0
ZMuMuEfficiency::standAloneIsoToken_
EDGetTokenT< IsolationCollection > standAloneIsoToken_
Definition: ZMuMuEfficiency.cc:64
edm::InputTag
Definition: InputTag.h:15
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
ZMuMuEfficiency::h_ZptGenPassed_
TH1D * h_ZptGenPassed_
Definition: ZMuMuEfficiency.cc:85
ZMuMuEfficiency::numberOfOverlappedStandAlone_
int numberOfOverlappedStandAlone_
Definition: ZMuMuEfficiency.cc:95
Generated for CMSSW Reference Manual by doxygen 1.8.16