CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
DiMuonHistograms Class Reference

#include <DiMuonHistograms.h>

Inheritance diagram for DiMuonHistograms:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 DiMuonHistograms (const edm::ParameterSet &pset)
 
 ~DiMuonHistograms () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 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 EDProduceroperator= (const EDProducer &)=delete
 

Private Attributes

int etaBBin
 
int etaBin
 
double etaBMax
 
double etaBMin
 
double EtaCutMax
 
double EtaCutMin
 
int etaEBin
 
double etaECMax
 
double etaECMin
 
std::string EtaName [3]
 
int etaOvlpBin
 
std::vector< MonitorElement * > GlbGlbMuon_HM
 
std::vector< MonitorElement * > GlbGlbMuon_LM
 
double HighMassMax
 
double HighMassMin
 
std::vector< MonitorElement * > LooseLooseMuon
 
std::vector< MonitorElement * > LooseLooseMuonBadFrac
 
double LowMassMax
 
double LowMassMin
 
std::vector< MonitorElement * > MediumMediumMuon
 
std::vector< MonitorElement * > MediumMediumMuonBadFrac
 
std::string metname
 
int nGlbGlb
 
int nLooseLoose
 
int nMediumMedium
 
int nTightTight
 
edm::ParameterSet parameters
 
std::vector< MonitorElement * > SoftSoftMuon
 
std::vector< MonitorElement * > SoftSoftMuonBadFrac
 
std::vector< MonitorElement * > StaTrkMuon_HM
 
std::vector< MonitorElement * > StaTrkMuon_LM
 
MonitorElementtest
 
edm::EDGetTokenT< reco::BeamSpottheBeamSpotLabel_
 
std::string theFolder
 
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
 
MuonServiceProxytheService
 
edm::EDGetTokenT< reco::VertexCollectiontheVertexLabel_
 
std::vector< MonitorElement * > TightTightMuon
 
std::vector< MonitorElement * > TightTightMuonBadFrac
 
std::vector< MonitorElement * > TrkTrkMuon_HM
 
std::vector< MonitorElement * > TrkTrkMuon_LM
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
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
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr< DQMEDAnalyzerGlobalCacheinitializeGlobalCache (edm::ParameterSet const &)
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 
unsigned int streamId_
 

Detailed Description

Class DiMuonHistograms

DQM monitoring for dimuon mass

Author: S.Folgueras, U. Oviedo

Definition at line 31 of file DiMuonHistograms.h.

Constructor & Destructor Documentation

◆ DiMuonHistograms()

DiMuonHistograms::DiMuonHistograms ( const edm::ParameterSet pset)

Definition at line 31 of file DiMuonHistograms.cc.

31  {
32  // initialise parameters:
33  parameters = pSet;
34 
35  // counter
36  nTightTight = 0;
37  nMediumMedium = 0;
38  nLooseLoose = 0;
39  nGlbGlb = 0;
40 
41  // declare consumes:
42  theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> >(parameters.getParameter<edm::InputTag>("MuonCollection"));
43  theVertexLabel_ = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
44 
45  theBeamSpotLabel_ = mayConsume<reco::BeamSpot>(parameters.getParameter<edm::InputTag>("BeamSpotLabel"));
46 
47  etaBin = parameters.getParameter<int>("etaBin");
48  etaBBin = parameters.getParameter<int>("etaBBin");
49  etaEBin = parameters.getParameter<int>("etaEBin");
50 
51  etaBMin = parameters.getParameter<double>("etaBMin");
52  etaBMax = parameters.getParameter<double>("etaBMax");
53  etaECMin = parameters.getParameter<double>("etaECMin");
54  etaECMax = parameters.getParameter<double>("etaECMax");
55 
56  LowMassMin = parameters.getParameter<double>("LowMassMin");
57  LowMassMax = parameters.getParameter<double>("LowMassMax");
58  HighMassMin = parameters.getParameter<double>("HighMassMin");
59  HighMassMax = parameters.getParameter<double>("HighMassMax");
60 
61  theFolder = parameters.getParameter<string>("folder");
62 }

References DQMOfflineHeavyIons_cff::etaBBin, etaBin(), diMuonHistograms_cfi::etaBMax, diMuonHistograms_cfi::etaBMin, DQMOfflineHeavyIons_cff::etaEBin, diMuonHistograms_cfi::etaECMax, diMuonHistograms_cfi::etaECMin, DQMOfflineHeavyIons_cff::HighMassMax, DQMOfflineHeavyIons_cff::HighMassMin, DQMOfflineHeavyIons_cff::LowMassMax, and DQMOfflineHeavyIons_cff::LowMassMin.

◆ ~DiMuonHistograms()

DiMuonHistograms::~DiMuonHistograms ( )
override

Definition at line 64 of file DiMuonHistograms.cc.

64 {}

Member Function Documentation

◆ analyze()

void DiMuonHistograms::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 142 of file DiMuonHistograms.cc.

142  {
143  LogTrace(metname) << "[DiMuonHistograms] Analyze the mu in different eta regions";
145  iEvent.getByToken(theMuonCollectionLabel_, muons);
146 
147  // =================================================================================
148  // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
149  reco::Vertex::Point posVtx;
150  reco::Vertex::Error errVtx;
151  unsigned int theIndexOfThePrimaryVertex = 999.;
152 
154  iEvent.getByToken(theVertexLabel_, vertex);
155  if (vertex.isValid()) {
156  for (unsigned int ind = 0; ind < vertex->size(); ++ind) {
157  if ((*vertex)[ind].isValid() && !((*vertex)[ind].isFake())) {
158  theIndexOfThePrimaryVertex = ind;
159  break;
160  }
161  }
162  }
163 
164  if (theIndexOfThePrimaryVertex < 100) {
165  posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
166  errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
167  } else {
168  LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
169 
170  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
171  iEvent.getByToken(theBeamSpotLabel_, recoBeamSpotHandle);
172  reco::BeamSpot bs = *recoBeamSpotHandle;
173 
174  posVtx = bs.position();
175  errVtx(0, 0) = bs.BeamWidthX();
176  errVtx(1, 1) = bs.BeamWidthY();
177  errVtx(2, 2) = bs.sigmaZ();
178  }
179 
180  const reco::Vertex vtx(posVtx, errVtx);
181 
182  if (!muons.isValid())
183  return;
184 
185  // Loop on muon collection
186  TLorentzVector Mu1, Mu2;
187  float charge = 99.;
188  float InvMass = -99.;
189 
190  //Eta regions
191  double EtaCutMin[] = {0, etaBMin, etaECMin};
192  double EtaCutMax[] = {2.4, etaBMax, etaECMax};
193 
194  for (edm::View<reco::Muon>::const_iterator muon1 = muons->begin(); muon1 != muons->end(); ++muon1) {
195  LogTrace(metname) << "[DiMuonHistograms] loop over 1st muon" << endl;
196 
197  // Loop on second muons to fill invariant mass plots
198  for (edm::View<reco::Muon>::const_iterator muon2 = muon1; muon2 != muons->end(); ++muon2) {
199  LogTrace(metname) << "[DiMuonHistograms] loop over 2nd muon" << endl;
200  if (muon1 == muon2)
201  continue;
202 
203  // Global-Global Muon
204  if (muon1->isGlobalMuon() && muon2->isGlobalMuon()) {
205  LogTrace(metname) << "[DiMuonHistograms] Glb-Glb pair" << endl;
206  reco::TrackRef recoCombinedGlbTrack1 = muon1->combinedMuon();
207  reco::TrackRef recoCombinedGlbTrack2 = muon2->combinedMuon();
208  Mu1.SetPxPyPzE(recoCombinedGlbTrack1->px(),
209  recoCombinedGlbTrack1->py(),
210  recoCombinedGlbTrack1->pz(),
211  recoCombinedGlbTrack1->p());
212  Mu2.SetPxPyPzE(recoCombinedGlbTrack2->px(),
213  recoCombinedGlbTrack2->py(),
214  recoCombinedGlbTrack2->pz(),
215  recoCombinedGlbTrack2->p());
216 
217  charge = recoCombinedGlbTrack1->charge() * recoCombinedGlbTrack2->charge();
218  if (charge < 0) {
219  InvMass = (Mu1 + Mu2).M();
220  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
221  if (fabs(recoCombinedGlbTrack1->eta()) > EtaCutMin[iEtaRegion] &&
222  fabs(recoCombinedGlbTrack1->eta()) < EtaCutMax[iEtaRegion] &&
223  fabs(recoCombinedGlbTrack2->eta()) > EtaCutMin[iEtaRegion] &&
224  fabs(recoCombinedGlbTrack2->eta()) < EtaCutMax[iEtaRegion]) {
225  if (InvMass < LowMassMax)
226  GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
227  if (InvMass > HighMassMin)
228  GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
229  }
230  }
231  }
232  // Also Tight-Tight Muon Selection
233  if (muon::isTightMuon(*muon1, vtx) && muon::isTightMuon(*muon2, vtx)) {
234  test->Fill(InvMass);
235  LogTrace(metname) << "[DiMuonHistograms] Tight-Tight pair" << endl;
236  if (charge < 0) {
237  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
238  if (fabs(recoCombinedGlbTrack1->eta()) > EtaCutMin[iEtaRegion] &&
239  fabs(recoCombinedGlbTrack1->eta()) < EtaCutMax[iEtaRegion] &&
240  fabs(recoCombinedGlbTrack2->eta()) > EtaCutMin[iEtaRegion] &&
241  fabs(recoCombinedGlbTrack2->eta()) < EtaCutMax[iEtaRegion]) {
242  if (InvMass > 55. && InvMass < 125.) {
243  TightTightMuon[iEtaRegion]->Fill(InvMass);
244  TightTightMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost() / muon1->innerTrack()->found());
245  }
246  }
247  }
248  }
249  }
250  // Also Medium-Medium Muon Selection
251  if (muon::isMediumMuon(*muon1) && muon::isMediumMuon(*muon2)) {
252  test->Fill(InvMass);
253  LogTrace(metname) << "[DiMuonHistograms] Medium-Medium pair" << endl;
254  if (charge < 0) {
255  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
256  if (fabs(recoCombinedGlbTrack1->eta()) > EtaCutMin[iEtaRegion] &&
257  fabs(recoCombinedGlbTrack1->eta()) < EtaCutMax[iEtaRegion] &&
258  fabs(recoCombinedGlbTrack2->eta()) > EtaCutMin[iEtaRegion] &&
259  fabs(recoCombinedGlbTrack2->eta()) < EtaCutMax[iEtaRegion]) {
260  if (InvMass > 55. && InvMass < 125.) {
261  MediumMediumMuon[iEtaRegion]->Fill(InvMass);
262  MediumMediumMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost() / muon1->innerTrack()->found());
263  }
264  }
265  }
266  }
267  }
268  // Also Loose-Loose Muon Selection
269  if (muon::isLooseMuon(*muon1) && muon::isLooseMuon(*muon2)) {
270  test->Fill(InvMass);
271  LogTrace(metname) << "[DiMuonHistograms] Loose-Loose pair" << endl;
272  if (charge < 0) {
273  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
274  if (fabs(recoCombinedGlbTrack1->eta()) > EtaCutMin[iEtaRegion] &&
275  fabs(recoCombinedGlbTrack1->eta()) < EtaCutMax[iEtaRegion] &&
276  fabs(recoCombinedGlbTrack2->eta()) > EtaCutMin[iEtaRegion] &&
277  fabs(recoCombinedGlbTrack2->eta()) < EtaCutMax[iEtaRegion]) {
278  if (InvMass > 55. && InvMass < 125.) {
279  LooseLooseMuon[iEtaRegion]->Fill(InvMass);
280  LooseLooseMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost() / muon1->innerTrack()->found());
281  }
282  }
283  }
284  }
285  }
286  }
287 
288  // Now check for STA-TRK
289  if (muon2->isStandAloneMuon() && muon1->isTrackerMuon()) {
290  LogTrace(metname) << "[DiMuonHistograms] STA-Trk pair" << endl;
291  reco::TrackRef recoStaTrack = muon2->standAloneMuon();
292  reco::TrackRef recoTrack = muon1->track();
293  Mu2.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(), recoStaTrack->pz(), recoStaTrack->p());
294  Mu1.SetPxPyPzE(recoTrack->px(), recoTrack->py(), recoTrack->pz(), recoTrack->p());
295 
296  charge = recoStaTrack->charge() * recoTrack->charge();
297  if (charge < 0) {
298  InvMass = (Mu1 + Mu2).M();
299  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
300  if (fabs(recoStaTrack->eta()) > EtaCutMin[iEtaRegion] &&
301  fabs(recoStaTrack->eta()) < EtaCutMax[iEtaRegion] && fabs(recoTrack->eta()) > EtaCutMin[iEtaRegion] &&
302  fabs(recoTrack->eta()) < EtaCutMax[iEtaRegion]) {
303  if (InvMass < LowMassMax)
304  StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
305  if (InvMass > HighMassMin)
306  StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
307  }
308  }
309  }
310  }
311  if (muon1->isStandAloneMuon() && muon2->isTrackerMuon()) {
312  LogTrace(metname) << "[DiMuonHistograms] STA-Trk pair" << endl;
313  reco::TrackRef recoStaTrack = muon1->standAloneMuon();
314  reco::TrackRef recoTrack = muon2->track();
315  Mu1.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(), recoStaTrack->pz(), recoStaTrack->p());
316  Mu2.SetPxPyPzE(recoTrack->px(), recoTrack->py(), recoTrack->pz(), recoTrack->p());
317 
318  charge = recoStaTrack->charge() * recoTrack->charge();
319  if (charge < 0) {
320  InvMass = (Mu1 + Mu2).M();
321  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
322  if (fabs(recoStaTrack->eta()) > EtaCutMin[iEtaRegion] &&
323  fabs(recoStaTrack->eta()) < EtaCutMax[iEtaRegion] && fabs(recoTrack->eta()) > EtaCutMin[iEtaRegion] &&
324  fabs(recoTrack->eta()) < EtaCutMax[iEtaRegion]) {
325  if (InvMass < LowMassMax)
326  StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
327  if (InvMass > HighMassMin)
328  StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
329  }
330  }
331  }
332  }
333 
334  // TRK-TRK dimuon
335  if (muon1->isTrackerMuon() && muon2->isTrackerMuon()) {
336  LogTrace(metname) << "[DiMuonHistograms] Trk-Trk dimuon pair" << endl;
337  reco::TrackRef recoTrack2 = muon2->track();
338  reco::TrackRef recoTrack1 = muon1->track();
339  Mu2.SetPxPyPzE(recoTrack2->px(), recoTrack2->py(), recoTrack2->pz(), recoTrack2->p());
340  Mu1.SetPxPyPzE(recoTrack1->px(), recoTrack1->py(), recoTrack1->pz(), recoTrack1->p());
341 
342  charge = recoTrack1->charge() * recoTrack2->charge();
343  if (charge < 0) {
344  InvMass = (Mu1 + Mu2).M();
345  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
346  if (fabs(recoTrack1->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta()) < EtaCutMax[iEtaRegion] &&
347  fabs(recoTrack2->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta()) < EtaCutMax[iEtaRegion]) {
348  if (InvMass < LowMassMax)
349  TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
350  if (InvMass > HighMassMin)
351  TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
352  }
353  }
354  }
355 
356  LogTrace(metname) << "[DiMuonHistograms] Soft-Soft pair" << endl;
357 
358  if (muon::isSoftMuon(*muon1, vtx) && muon::isSoftMuon(*muon2, vtx)) {
359  if (charge < 0) {
360  InvMass = (Mu1 + Mu2).M();
361  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
362  if (fabs(recoTrack1->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta()) < EtaCutMax[iEtaRegion] &&
363  fabs(recoTrack2->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta()) < EtaCutMax[iEtaRegion]) {
364  SoftSoftMuon[iEtaRegion]->Fill(InvMass);
365  SoftSoftMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost() / muon1->innerTrack()->found());
366  }
367  }
368  }
369  }
370  }
371  } //muon2
372  } //Muon1
373 }

References cms::cuda::bs, ALCARECOTkAlJpsiMuMu_cff::charge, relativeConstraints::error, diMuonHistograms_cfi::etaBMax, diMuonHistograms_cfi::etaBMin, diMuonHistograms_cfi::etaECMax, diMuonHistograms_cfi::etaECMin, DQMOfflineHeavyIons_cff::HighMassMin, iEvent, muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), muon::isTightMuon(), LogTrace, DQMOfflineHeavyIons_cff::LowMassMax, metname, PDWG_BPHSkim_cff::muons, position, bphysicsOniaDQM_cfi::vertex, and extraflags_cff::vtx.

◆ bookHistograms()

void DiMuonHistograms::bookHistograms ( DQMStore::IBooker ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 66 of file DiMuonHistograms.cc.

68  {
69  ibooker.cd();
70  ibooker.setCurrentFolder(theFolder);
71 
72  int nBin[3] = {etaBin, etaBBin, etaEBin};
73  EtaName[0] = "";
74  EtaName[1] = "_Barrel";
75  EtaName[2] = "_EndCap";
76  test = ibooker.book1D("test", "InvMass_{Tight,Tight}", 100, 0., 200.);
77  for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
78  GlbGlbMuon_LM.push_back(ibooker.book1D("GlbGlbMuon_LM" + EtaName[iEtaRegion],
79  "InvMass_{GLB,GLB}" + EtaName[iEtaRegion],
80  nBin[iEtaRegion],
81  LowMassMin,
82  LowMassMax));
83  TrkTrkMuon_LM.push_back(ibooker.book1D("TrkTrkMuon_LM" + EtaName[iEtaRegion],
84  "InvMass_{TRK,TRK}" + EtaName[iEtaRegion],
85  nBin[iEtaRegion],
86  LowMassMin,
87  LowMassMax));
88  StaTrkMuon_LM.push_back(ibooker.book1D("StaTrkMuon_LM" + EtaName[iEtaRegion],
89  "InvMass_{STA,TRK}" + EtaName[iEtaRegion],
90  nBin[iEtaRegion],
91  LowMassMin,
92  LowMassMax));
93 
94  GlbGlbMuon_HM.push_back(ibooker.book1D("GlbGlbMuon_HM" + EtaName[iEtaRegion],
95  "InvMass_{GLB,GLB}" + EtaName[iEtaRegion],
96  nBin[iEtaRegion],
98  HighMassMax));
99  TrkTrkMuon_HM.push_back(ibooker.book1D("TrkTrkMuon_HM" + EtaName[iEtaRegion],
100  "InvMass_{TRK,TRK}" + EtaName[iEtaRegion],
101  nBin[iEtaRegion],
102  HighMassMin,
103  HighMassMax));
104  StaTrkMuon_HM.push_back(ibooker.book1D("StaTrkMuon_HM" + EtaName[iEtaRegion],
105  "InvMass_{STA,TRK}" + EtaName[iEtaRegion],
106  nBin[iEtaRegion],
107  HighMassMin,
108  HighMassMax));
109 
110  // arround the Z peak
111  TightTightMuon.push_back(ibooker.book1D("TightTightMuon" + EtaName[iEtaRegion],
112  "InvMass_{Tight,Tight}" + EtaName[iEtaRegion],
113  nBin[iEtaRegion],
114  LowMassMin,
115  LowMassMax));
116  MediumMediumMuon.push_back(ibooker.book1D("MediumMediumMuon" + EtaName[iEtaRegion],
117  "InvMass_{Medium,Medium}" + EtaName[iEtaRegion],
118  nBin[iEtaRegion],
119  LowMassMin,
120  LowMassMax));
121  LooseLooseMuon.push_back(ibooker.book1D("LooseLooseMuon" + EtaName[iEtaRegion],
122  "InvMass_{Loose,Loose}" + EtaName[iEtaRegion],
123  nBin[iEtaRegion],
124  LowMassMin,
125  LowMassMax));
126  //Fraction of bad hits in the tracker track to the total
127  TightTightMuonBadFrac.push_back(ibooker.book1D(
128  "TightTightMuonBadFrac" + EtaName[iEtaRegion], "BadFrac_{Tight,Tight}" + EtaName[iEtaRegion], 10, 0, 0.4));
129  MediumMediumMuonBadFrac.push_back(ibooker.book1D(
130  "MediumMediumMuonBadFrac" + EtaName[iEtaRegion], "BadFrac_{Medium,Medium}" + EtaName[iEtaRegion], 10, 0, 0.4));
131  LooseLooseMuonBadFrac.push_back(ibooker.book1D(
132  "LooseLooseMuonBadFrac" + EtaName[iEtaRegion], "BadFrac_{Loose,Loose}" + EtaName[iEtaRegion], 10, 0, 0.4));
133 
134  // low-mass resonances
135  SoftSoftMuon.push_back(ibooker.book1D(
136  "SoftSoftMuon" + EtaName[iEtaRegion], "InvMass_{Soft,Soft}" + EtaName[iEtaRegion], nBin[iEtaRegion], 0.0, 55.0));
137  SoftSoftMuonBadFrac.push_back(ibooker.book1D(
138  "SoftSoftMuonBadFrac" + EtaName[iEtaRegion], "BadFrac_{Soft,Soft}" + EtaName[iEtaRegion], 10, 0, 0.4));
139  }
140 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), DQMOfflineHeavyIons_cff::etaBBin, etaBin(), DQMOfflineHeavyIons_cff::etaEBin, DQMOfflineHeavyIons_cff::HighMassMax, DQMOfflineHeavyIons_cff::HighMassMin, DQMOfflineHeavyIons_cff::LowMassMax, DQMOfflineHeavyIons_cff::LowMassMin, PFElectronDQMAnalyzer_cfi::nBin, and dqm::implementation::NavigatorBase::setCurrentFolder().

Member Data Documentation

◆ etaBBin

int DiMuonHistograms::etaBBin
private

Definition at line 52 of file DiMuonHistograms.h.

◆ etaBin

int DiMuonHistograms::etaBin
private

Definition at line 51 of file DiMuonHistograms.h.

◆ etaBMax

double DiMuonHistograms::etaBMax
private

Definition at line 62 of file DiMuonHistograms.h.

◆ etaBMin

double DiMuonHistograms::etaBMin
private

Definition at line 61 of file DiMuonHistograms.h.

◆ EtaCutMax

double DiMuonHistograms::EtaCutMax
private

Definition at line 60 of file DiMuonHistograms.h.

◆ EtaCutMin

double DiMuonHistograms::EtaCutMin
private

Definition at line 59 of file DiMuonHistograms.h.

◆ etaEBin

int DiMuonHistograms::etaEBin
private

Definition at line 53 of file DiMuonHistograms.h.

◆ etaECMax

double DiMuonHistograms::etaECMax
private

Definition at line 64 of file DiMuonHistograms.h.

◆ etaECMin

double DiMuonHistograms::etaECMin
private

Definition at line 63 of file DiMuonHistograms.h.

◆ EtaName

std::string DiMuonHistograms::EtaName[3]
private

Definition at line 57 of file DiMuonHistograms.h.

◆ etaOvlpBin

int DiMuonHistograms::etaOvlpBin
private

Definition at line 54 of file DiMuonHistograms.h.

◆ GlbGlbMuon_HM

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_HM
private

Definition at line 73 of file DiMuonHistograms.h.

◆ GlbGlbMuon_LM

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_LM
private

Definition at line 72 of file DiMuonHistograms.h.

◆ HighMassMax

double DiMuonHistograms::HighMassMax
private

Definition at line 70 of file DiMuonHistograms.h.

◆ HighMassMin

double DiMuonHistograms::HighMassMin
private

Definition at line 69 of file DiMuonHistograms.h.

◆ LooseLooseMuon

std::vector<MonitorElement*> DiMuonHistograms::LooseLooseMuon
private

Definition at line 79 of file DiMuonHistograms.h.

◆ LooseLooseMuonBadFrac

std::vector<MonitorElement*> DiMuonHistograms::LooseLooseMuonBadFrac
private

Definition at line 83 of file DiMuonHistograms.h.

◆ LowMassMax

double DiMuonHistograms::LowMassMax
private

Definition at line 68 of file DiMuonHistograms.h.

◆ LowMassMin

double DiMuonHistograms::LowMassMin
private

Definition at line 67 of file DiMuonHistograms.h.

◆ MediumMediumMuon

std::vector<MonitorElement*> DiMuonHistograms::MediumMediumMuon
private

Definition at line 80 of file DiMuonHistograms.h.

◆ MediumMediumMuonBadFrac

std::vector<MonitorElement*> DiMuonHistograms::MediumMediumMuonBadFrac
private

Definition at line 84 of file DiMuonHistograms.h.

◆ metname

std::string DiMuonHistograms::metname
private

Definition at line 48 of file DiMuonHistograms.h.

◆ nGlbGlb

int DiMuonHistograms::nGlbGlb
private

Definition at line 100 of file DiMuonHistograms.h.

◆ nLooseLoose

int DiMuonHistograms::nLooseLoose
private

Definition at line 99 of file DiMuonHistograms.h.

◆ nMediumMedium

int DiMuonHistograms::nMediumMedium
private

Definition at line 98 of file DiMuonHistograms.h.

◆ nTightTight

int DiMuonHistograms::nTightTight
private

Definition at line 97 of file DiMuonHistograms.h.

◆ parameters

edm::ParameterSet DiMuonHistograms::parameters
private

Definition at line 45 of file DiMuonHistograms.h.

◆ SoftSoftMuon

std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuon
private

Definition at line 82 of file DiMuonHistograms.h.

◆ SoftSoftMuonBadFrac

std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuonBadFrac
private

Definition at line 86 of file DiMuonHistograms.h.

◆ StaTrkMuon_HM

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_HM
private

Definition at line 75 of file DiMuonHistograms.h.

◆ StaTrkMuon_LM

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_LM
private

Definition at line 74 of file DiMuonHistograms.h.

◆ test

MonitorElement* DiMuonHistograms::test
private

◆ theBeamSpotLabel_

edm::EDGetTokenT<reco::BeamSpot> DiMuonHistograms::theBeamSpotLabel_
private

Definition at line 93 of file DiMuonHistograms.h.

◆ theFolder

std::string DiMuonHistograms::theFolder
private

Definition at line 95 of file DiMuonHistograms.h.

◆ theMuonCollectionLabel_

edm::EDGetTokenT<edm::View<reco::Muon> > DiMuonHistograms::theMuonCollectionLabel_
private

Definition at line 91 of file DiMuonHistograms.h.

◆ theService

MuonServiceProxy* DiMuonHistograms::theService
private

Definition at line 44 of file DiMuonHistograms.h.

◆ theVertexLabel_

edm::EDGetTokenT<reco::VertexCollection> DiMuonHistograms::theVertexLabel_
private

Definition at line 92 of file DiMuonHistograms.h.

◆ TightTightMuon

std::vector<MonitorElement*> DiMuonHistograms::TightTightMuon
private

Definition at line 81 of file DiMuonHistograms.h.

◆ TightTightMuonBadFrac

std::vector<MonitorElement*> DiMuonHistograms::TightTightMuonBadFrac
private

Definition at line 85 of file DiMuonHistograms.h.

◆ TrkTrkMuon_HM

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_HM
private

Definition at line 77 of file DiMuonHistograms.h.

◆ TrkTrkMuon_LM

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_LM
private

Definition at line 76 of file DiMuonHistograms.h.

DiMuonHistograms::TrkTrkMuon_HM
std::vector< MonitorElement * > TrkTrkMuon_HM
Definition: DiMuonHistograms.h:77
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
DiMuonHistograms::nMediumMedium
int nMediumMedium
Definition: DiMuonHistograms.h:98
DiMuonHistograms::GlbGlbMuon_HM
std::vector< MonitorElement * > GlbGlbMuon_HM
Definition: DiMuonHistograms.h:73
muon::isSoftMuon
bool isSoftMuon(const reco::Muon &, const reco::Vertex &, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:916
DiMuonHistograms::TrkTrkMuon_LM
std::vector< MonitorElement * > TrkTrkMuon_LM
Definition: DiMuonHistograms.h:76
DiMuonHistograms::nGlbGlb
int nGlbGlb
Definition: DiMuonHistograms.h:100
DiMuonHistograms::HighMassMax
double HighMassMax
Definition: DiMuonHistograms.h:70
reco::Vertex::Error
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:44
DiMuonHistograms::nLooseLoose
int nLooseLoose
Definition: DiMuonHistograms.h:99
DiMuonHistograms::LowMassMax
double LowMassMax
Definition: DiMuonHistograms.h:68
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
DiMuonHistograms::theMuonCollectionLabel_
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
Definition: DiMuonHistograms.h:91
DiMuonHistograms::TightTightMuon
std::vector< MonitorElement * > TightTightMuon
Definition: DiMuonHistograms.h:81
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
edm::Handle
Definition: AssociativeIterator.h:50
relativeConstraints.error
error
Definition: relativeConstraints.py:53
PFElectronDQMAnalyzer_cfi.nBin
nBin
Definition: PFElectronDQMAnalyzer_cfi.py:25
edm::Ref< TrackCollection >
DiMuonHistograms::etaBMin
double etaBMin
Definition: DiMuonHistograms.h:61
cms::cuda::bs
bs
Definition: HistoContainer.h:76
DiMuonHistograms::etaBin
int etaBin
Definition: DiMuonHistograms.h:51
test
Definition: SmallWORMDict.h:13
DiMuonHistograms::HighMassMin
double HighMassMin
Definition: DiMuonHistograms.h:69
muon::isLooseMuon
bool isLooseMuon(const reco::Muon &)
Definition: MuonSelectors.cc:895
reco::BeamSpot
Definition: BeamSpot.h:21
dqm::implementation::NavigatorBase::cd
virtual void cd()
Definition: DQMStore.cc:29
DiMuonHistograms::MediumMediumMuonBadFrac
std::vector< MonitorElement * > MediumMediumMuonBadFrac
Definition: DiMuonHistograms.h:84
DiMuonHistograms::etaECMax
double etaECMax
Definition: DiMuonHistograms.h:64
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
DiMuonHistograms::etaBBin
int etaBBin
Definition: DiMuonHistograms.h:52
DiMuonHistograms::StaTrkMuon_LM
std::vector< MonitorElement * > StaTrkMuon_LM
Definition: DiMuonHistograms.h:74
DiMuonHistograms::LooseLooseMuonBadFrac
std::vector< MonitorElement * > LooseLooseMuonBadFrac
Definition: DiMuonHistograms.h:83
bphysicsOniaDQM_cfi.vertex
vertex
Definition: bphysicsOniaDQM_cfi.py:7
DiMuonHistograms::nTightTight
int nTightTight
Definition: DiMuonHistograms.h:97
DiMuonHistograms::EtaName
std::string EtaName[3]
Definition: DiMuonHistograms.h:57
DiMuonHistograms::TightTightMuonBadFrac
std::vector< MonitorElement * > TightTightMuonBadFrac
Definition: DiMuonHistograms.h:85
DiMuonHistograms::theFolder
std::string theFolder
Definition: DiMuonHistograms.h:95
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:880
position
static int position[264][3]
Definition: ReadPGInfo.cc:289
iEvent
int iEvent
Definition: GenABIO.cc:224
DiMuonHistograms::GlbGlbMuon_LM
std::vector< MonitorElement * > GlbGlbMuon_LM
Definition: DiMuonHistograms.h:72
DiMuonHistograms::theBeamSpotLabel_
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
Definition: DiMuonHistograms.h:93
DiMuonHistograms::SoftSoftMuonBadFrac
std::vector< MonitorElement * > SoftSoftMuonBadFrac
Definition: DiMuonHistograms.h:86
DiMuonHistograms::theVertexLabel_
edm::EDGetTokenT< reco::VertexCollection > theVertexLabel_
Definition: DiMuonHistograms.h:92
muon::isMediumMuon
bool isMediumMuon(const reco::Muon &, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:899
DiMuonHistograms::LooseLooseMuon
std::vector< MonitorElement * > LooseLooseMuon
Definition: DiMuonHistograms.h:79
DiMuonHistograms::SoftSoftMuon
std::vector< MonitorElement * > SoftSoftMuon
Definition: DiMuonHistograms.h:82
DiMuonHistograms::etaECMin
double etaECMin
Definition: DiMuonHistograms.h:63
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:40
DiMuonHistograms::EtaCutMax
double EtaCutMax
Definition: DiMuonHistograms.h:60
DiMuonHistograms::StaTrkMuon_HM
std::vector< MonitorElement * > StaTrkMuon_HM
Definition: DiMuonHistograms.h:75
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:19
DiMuonHistograms::etaBMax
double etaBMax
Definition: DiMuonHistograms.h:62
DiMuonHistograms::metname
std::string metname
Definition: DiMuonHistograms.h:48
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
DiMuonHistograms::etaEBin
int etaEBin
Definition: DiMuonHistograms.h:53
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
DiMuonHistograms::EtaCutMin
double EtaCutMin
Definition: DiMuonHistograms.h:59
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:234
DiMuonHistograms::parameters
edm::ParameterSet parameters
Definition: DiMuonHistograms.h:45
edm::InputTag
Definition: InputTag.h:15
DiMuonHistograms::LowMassMin
double LowMassMin
Definition: DiMuonHistograms.h:67
reco::Vertex
Definition: Vertex.h:35
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
DiMuonHistograms::MediumMediumMuon
std::vector< MonitorElement * > MediumMediumMuon
Definition: DiMuonHistograms.h:80