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

#include <MuonIsolationDQM.h>

Inheritance diagram for MuonIsolationDQM:
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
 
 MuonIsolationDQM (const edm::ParameterSet &)
 
 ~MuonIsolationDQM () 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
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndRuns () const final
 

Private Types

enum  { NOAXIS, XAXIS, YAXIS, ZAXIS }
 
typedef edm::Handle< reco::IsoDepositMap > MuIsoDepHandle
 
typedef edm::RefToBase< reco::Muon > MuonBaseRef
 
typedef edm::View< reco::Muon >::const_iterator MuonIterator
 

Private Member Functions

void FillHistos (int)
 
void FillNVtxHistos (int)
 
int GetNVtxBin (int)
 
TH1 * GetTH1FromMonitorElement (MonitorElement *me)
 
void InitHistos ()
 
void InitStatics ()
 
void NormalizeHistos ()
 
void RecordData (const reco::Muon &muon)
 

Private Attributes

std::vector< std::string > axis_titles
 
std::vector< std::string > axis_titles_NVtxs
 
std::string dirName
 
std::vector< MonitorElement * > h_1D
 
std::vector< MonitorElement * > h_1D_NVTX
 
std::vector< MonitorElement * > h_2D
 
MonitorElement * h_nMuons
 
std::vector< int > isContinuous
 
double L_BIN_WIDTH
 
int LOG_BINNING_ENABLED
 
double LOG_BINNING_RATIO
 
std::vector< std::string > main_titles
 
std::vector< std::string > main_titles_NVtxs
 
const typedef reco::IsoDeposit MuIsoDepRef
 
std::vector< std::string > names
 
std::vector< std::string > names_2D
 
std::vector< std::string > names_NVtxs
 
int nEvents
 
int nGLBMuons
 
int nSTAMuons
 
int nTRKMuons
 
int NUM_LOG_BINS
 
std::vector< std::vector< double > > param
 
bool requireGLBMuon
 
bool requireSTAMuon
 
bool requireTRKMuon
 
std::string rootfilename
 
double S_BIN_WIDTH
 
double theData [NUM_VARS]
 
double theData2D [NUM_VARS_2D]
 
double theDataNVtx [NUM_VARS_NVTX]
 
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
 
edm::EDGetTokenT< reco::VertexCollection > theVertexCollectionLabel_
 
std::string title_cone
 
std::string title_sam
 
std::vector< std::string > titles_2D
 

Static Private Attributes

static const int NUM_VARS = 48
 
static const int NUM_VARS_2D = 10
 
static const int NUM_VARS_NVTX = 6
 

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 >
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache 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< DQMEDAnalyzerGlobalCache > initializeGlobalCache (edm::ParameterSet const &)
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 
unsigned int streamId_
 

Detailed Description

Definition at line 59 of file MuonIsolationDQM.h.

Member Typedef Documentation

◆ MuIsoDepHandle

typedef edm::Handle<reco::IsoDepositMap> MuonIsolationDQM::MuIsoDepHandle
private

Definition at line 63 of file MuonIsolationDQM.h.

◆ MuonBaseRef

typedef edm::RefToBase<reco::Muon> MuonIsolationDQM::MuonBaseRef
private

Definition at line 62 of file MuonIsolationDQM.h.

◆ MuonIterator

typedef edm::View<reco::Muon>::const_iterator MuonIsolationDQM::MuonIterator
private

Definition at line 61 of file MuonIsolationDQM.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
private
Enumerator
NOAXIS 
XAXIS 
YAXIS 
ZAXIS 

Definition at line 151 of file MuonIsolationDQM.h.

151 { NOAXIS, XAXIS, YAXIS, ZAXIS };

Constructor & Destructor Documentation

◆ MuonIsolationDQM()

MuonIsolationDQM::MuonIsolationDQM ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 55 of file MuonIsolationDQM.cc.

55  {
56 #ifdef DEBUG
57  cout << " Initialise Constructor " << endl;
58 #endif
59  requireSTAMuon = iConfig.getUntrackedParameter<bool>("requireSTAMuon");
60  requireTRKMuon = iConfig.getUntrackedParameter<bool>("requireTRKMuon");
61  requireGLBMuon = iConfig.getUntrackedParameter<bool>("requireGLBMuon");
62  dirName = iConfig.getParameter<std::string>("directory");
63 
64  //--------Initialize tags-------
65  theMuonCollectionLabel_ =
66  consumes<edm::View<reco::Muon> >(iConfig.getUntrackedParameter<edm::InputTag>("Global_Muon_Label"));
67  theVertexCollectionLabel_ =
68  consumes<reco::VertexCollection>(iConfig.getUntrackedParameter<edm::InputTag>("vertexLabel"));
69 
70  //-------Initialize Counterse----------------
71  nEvents = 0;
72  nSTAMuons = 0;
73  nTRKMuons = 0;
74  nGLBMuons = 0;
75 
76  InitStatics();
77 
78  //------"allocate" space for the data vectors-------
79  h_1D.resize(NUM_VARS);
80  h_2D.resize(NUM_VARS_2D);
81  h_1D_NVTX.resize(NUM_VARS_NVTX);
82 }

References gather_cfg::cout, TrackerOfflineValidation_Dqm_cff::dirName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nEvents, muonIsolationDQM_cff::requireGLBMuon, muonIsolationDQM_cff::requireSTAMuon, muonIsolationDQM_cff::requireTRKMuon, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~MuonIsolationDQM()

MuonIsolationDQM::~MuonIsolationDQM ( )
override

Definition at line 87 of file MuonIsolationDQM.cc.

87  {
88 #ifdef DEBUG
89  cout << "Calling destructor" << endl;
90 #endif
91  //Deallocate memory
92 }

References gather_cfg::cout.

Member Function Documentation

◆ analyze()

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

Reimplemented from DQMEDAnalyzer.

Definition at line 595 of file MuonIsolationDQM.cc.

595  {
596  ++nEvents;
597  edm::LogInfo("Tutorial") << "\nInvestigating event #" << nEvents << "\n";
598 #ifdef DEBUG
599  cout << "[MuonIsolationDQM]: analyze()" << endl;
600 #endif
601 
602  // Get Muon Collection
603  edm::Handle<edm::View<reco::Muon> > muons;
604  iEvent.getByToken(theMuonCollectionLabel_, muons);
605 
606 #ifdef DEBUG
607  cout << "[MuonIsolationDQM]: Number of muons -> " << muons->size() << endl;
608 #endif
609 
610  int theMuonData = muons->size();
611  h_nMuons->Fill(theMuonData);
612 #ifdef DEBUG
613  cout << "[MuonIsolationDQM]: Vertex is Valid" << endl;
614 #endif
615 
616  //Get Vertex Information
617  int _numPV = 0;
618  edm::Handle<reco::VertexCollection> vertexHandle;
619  iEvent.getByToken(theVertexCollectionLabel_, vertexHandle);
620 
621  if (vertexHandle.isValid()) {
622  reco::VertexCollection vertex = *(vertexHandle.product());
623  for (reco::VertexCollection::const_iterator v = vertex.begin(); v != vertex.end(); ++v) {
624  if (v->isFake())
625  continue;
626  if (v->ndof() < 4)
627  continue;
628  if (fabs(v->z()) > 24.0)
629  continue;
630  ++_numPV;
631  }
632  }
633 
634 #ifdef DEBUG
635  cout << "[MuonIsolationDQM]: Vertex is Valid" << endl;
636 #endif
637  // Get Muon Collection
638 
639  //Fill historgams concerning muon isolation
640  for (edm::View<reco::Muon>::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
641  if (requireSTAMuon && muon->isStandAloneMuon()) {
642  ++nSTAMuons;
643  RecordData(*muon);
644  FillHistos(_numPV);
645  } else if (requireTRKMuon && muon->isTrackerMuon()) {
646  ++nTRKMuons;
647  RecordData(*muon);
648  FillHistos(_numPV);
649  } else if (requireGLBMuon && muon->isGlobalMuon()) {
650  ++nGLBMuons;
651  RecordData(*muon);
652  FillHistos(_numPV);
653  FillNVtxHistos(_numPV);
654  }
655  }
656 }

References gather_cfg::cout, iEvent, edm::HandleBase::isValid(), HLT_2018_cff::muon, PDWG_BPHSkim_cff::muons, nEvents, edm::Handle< T >::product(), muonIsolationDQM_cff::requireGLBMuon, muonIsolationDQM_cff::requireSTAMuon, muonIsolationDQM_cff::requireTRKMuon, findQualityFiles::v, and bphysicsOniaDQM_cfi::vertex.

◆ bookHistograms()

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

GetTH1FromMonitorElement(h_1D_NVTX[var])->Sumw2();

Implements DQMEDAnalyzer.

Definition at line 758 of file MuonIsolationDQM.cc.

760  {
761  ibooker.cd();
762  ibooker.setCurrentFolder(dirName);
763 
764  ibooker.cd();
765  ibooker.setCurrentFolder(dirName);
766 
767  //---initialize number of muons histogram---
768  h_nMuons = ibooker.book1D("nMuons", title_sam + "Number of Muons", 20, 0., 20.);
769  h_nMuons->setAxisTitle("Number of Muons", XAXIS);
770  h_nMuons->setAxisTitle("Fraction of Events", YAXIS);
771 
772  //---Initialize 1D Histograms---
773  for (int var = 0; var < NUM_VARS; var++) {
774  h_1D[var] = ibooker.book1D(
775  names[var], title_sam + main_titles[var] + title_cone, (int)param[var][0], param[var][1], param[var][2]);
776  h_1D[var]->setAxisTitle(axis_titles[var], XAXIS);
777  // GetTH1FromMonitorElement(h_1D[var])->Sumw2();
778  } //Finish 1D
779 
780  //----Initialize 2D Histograms
781  for (int var = 0; var < NUM_VARS_2D; var++) {
782  h_2D[var] = ibooker.bookProfile(names_2D[var] + "_VsPV", titles_2D[var] + " Vs PV", 50, 0.5, 50.5, 20, 0.0, 20.0);
783 
784  h_2D[var]->setAxisTitle("Number of PV", XAXIS);
785  h_2D[var]->setAxisTitle(titles_2D[var] + " (GeV)", YAXIS);
786  // h_2D[var]->getTH1()->Sumw2();
787  }
788 
789  //-----Initialise PU-Binned histograms
790  for (int var = 0; var < NUM_VARS_NVTX; var++) {
791  h_1D_NVTX[var] = ibooker.book1D(names_NVtxs[var], main_titles_NVtxs[var], 50, 0.0, 10.0);
792  h_1D_NVTX[var]->setAxisTitle(axis_titles_NVtxs[var], XAXIS);
794  }
795 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), dqm::implementation::NavigatorBase::cd(), TrackerOfflineValidation_Dqm_cff::dirName, names, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), and trigObjTnPSource_cfi::var.

◆ FillHistos()

void MuonIsolationDQM::FillHistos ( int  numPV)
private

Definition at line 804 of file MuonIsolationDQM.cc.

804  {
805 #ifdef DEBUG
806  cout << "FillHistos( " << numPV << " )" << endl;
807 #endif
808 
809  //----------Fill 1D histograms---------------
810  for (int var = 0; var < NUM_VARS; var++) {
811  h_1D[var]->Fill(theData[var]);
812  // cd_plots[var]->Fill(theData[var]);//right now, this is a regular PDF (just like h_1D)
813  //OFBin if (theData[var] > param[var][2]) {
814  //OFBin // fill the overflow bin
815  //OFBin overFlowBin = (int) param[var][0] + 1;
816  //OFBin overFlow = GetTH1FromMonitorElement(h_1D[var])->GetBinContent(overFlowBin);
817  //OFBin GetTH1FromMonitorElement(h_1D[var])->SetBinContent(overFlowBin, overFlow + 1);
818  //OFBin }
819  } //Finish 1D
820 
821  for (int var = 0; var < NUM_VARS_2D; var++) {
822  h_2D[var]->Fill(numPV, theData2D[var]);
823  }
824 
825 #ifdef DEBUG
826  cout << "FillHistos( " << numPV << " ): DONE" << endl;
827 #endif
828 }

References gather_cfg::cout, and trigObjTnPSource_cfi::var.

◆ FillNVtxHistos()

void MuonIsolationDQM::FillNVtxHistos ( int  PV)
private

Definition at line 829 of file MuonIsolationDQM.cc.

829  {
830  if (PV < 15) {
831  h_1D_NVTX[0]->Fill(theDataNVtx[0]);
832  h_1D_NVTX[3]->Fill(theDataNVtx[3]);
833  }
834  if (PV >= 15 && PV < 30) {
835  h_1D_NVTX[1]->Fill(theDataNVtx[1]);
836  h_1D_NVTX[4]->Fill(theDataNVtx[4]);
837  }
838  if (PV >= 30) {
839  h_1D_NVTX[2]->Fill(theDataNVtx[2]);
840  h_1D_NVTX[5]->Fill(theDataNVtx[5]);
841  }
842 }

References nanoDQM_cfi::PV.

◆ GetNVtxBin()

int MuonIsolationDQM::GetNVtxBin ( int  )
private

◆ GetTH1FromMonitorElement()

TH1 * MuonIsolationDQM::GetTH1FromMonitorElement ( MonitorElement *  me)
private

Definition at line 844 of file MuonIsolationDQM.cc.

844 { return me->getTH1(); }

References hlt_dqm_clientPB-live_cfg::me.

◆ InitHistos()

void MuonIsolationDQM::InitHistos ( )
private

◆ InitStatics()

void MuonIsolationDQM::InitStatics ( )
private

Definition at line 97 of file MuonIsolationDQM.cc.

97  {
98 #ifdef DEBUG
99  cout << " InitStatistics() " << endl;
100 #endif
101  //-----------Initialize primitives-----------
102  S_BIN_WIDTH = 1.0; //in GeV
103  L_BIN_WIDTH = 2.0; //in GeV
104  LOG_BINNING_ENABLED = 1;
105  NUM_LOG_BINS = 15;
106  LOG_BINNING_RATIO = 1.1;
107  //ratio by which each bin is wider than the last for log binning
108  //i.e. bin widths are (x), (r*x), (r^2*x), ..., (r^(nbins)*x)
109 
110  //-------Initialize Titles---------
111  title_sam = ""; //"[Sample b-jet events] ";
112  title_cone = ""; //" [in R=0.3 IsoDeposit Cone]";
113  //The above two pieces of info will be printed on the title of the whole page,
114  //not for each individual histogram
115  // title_cd = "C.D. of ";
116 
117  //-------"Allocate" memory for vectors
118  main_titles.resize(NUM_VARS);
119  axis_titles.resize(NUM_VARS);
120  names.resize(NUM_VARS);
121  param.resize(NUM_VARS, vector<double>(3));
122  isContinuous.resize(NUM_VARS);
123 
124  titles_2D.resize(NUM_VARS_2D);
125  names_2D.resize(NUM_VARS_2D);
126 
127  main_titles_NVtxs.resize(NUM_VARS_NVTX);
128  axis_titles_NVtxs.resize(NUM_VARS_NVTX);
129  names_NVtxs.resize(NUM_VARS_NVTX);
130 
131 #ifdef DEBUG
132  cout << "InitStatistics(): vectors resized " << endl;
133 #endif
134  //-----Titles of the plots-----------
135  main_titles[0] = "Total Tracker Momentum, #Delta R = 0.3";
136  main_titles[1] = "Total EM Cal Energy, #Delta R = 0.3";
137  main_titles[2] = "Total Had Cal Energy, #Delta R = 0.3";
138  main_titles[3] = "Total HO Cal Energy, #Delta R = 0.3";
139  main_titles[4] = "Number of Tracker Tracks, #Delta R = 0.3";
140  main_titles[5] = "Number of Jets around Muon, #Delta R = 0.3";
141  main_titles[6] = "Tracker p_{T} within veto cone, #Delta R = 0.3";
142  main_titles[7] = "EM E_{T} within veto cone, #Delta R = 0.3";
143  main_titles[8] = "Had E_{T} within veto cone, #Delta R = 0.3";
144  main_titles[9] = "HO E_{T} within veto cone, #Delta R = 0.3";
145  main_titles[10] = "Average Momentum per Track, #Delta R = 0.3";
146  main_titles[11] = "Weighted Energy, #Delta R = 0.3";
147 
148  main_titles[12] = "Total Tracker Momentum, #Delta R = 0.5";
149  main_titles[13] = "Total EM Cal Energy, #Delta R = 0.5";
150  main_titles[14] = "Total Had Cal Energy, #Delta R = 0.5";
151  main_titles[15] = "Total HO Cal Energy, #Delta R = 0.5";
152  main_titles[16] = "Number of Tracker Tracks, #Delta R = 0.5";
153  main_titles[17] = "Number of Jets around Muon, #Delta R = 0.5";
154  main_titles[18] = "Tracker p_{T} within veto cone, #Delta R = 0.5";
155  main_titles[19] = "EM E_{T} within veto cone, #Delta R = 0.5";
156  main_titles[20] = "Had E_{T} within veto cone, #Delta R = 0.5";
157  main_titles[21] = "HO E_{T} within veto cone, #Delta R = 0.5";
158  main_titles[22] = "Average Momentum per Track, #Delta R = 0.5";
159  main_titles[23] = "Weighted Energy, #Delta R = 0.5";
160 
161  main_titles[24] = "Relative Detector-Based Isolation, #Delta R = 0.3";
162  main_titles[25] = "Relative Detector-Based Isolation, #Delta R = 0.5";
163 
164  //-----Titles of the plots-----------
165  main_titles[26] = "Sum PF Charged Hadron Pt, #Delta R = 0.3";
166  main_titles[27] = "Sum PF Neutral Hadron Pt, #Delta R = 0.3";
167  main_titles[28] = "Sum PF Photon Et, #Delta R = 0.3";
168  main_titles[29] = "Sum PF Neutral Hadron Pt (Higher Pt threshold), #Delta R = 0.3";
169  main_titles[30] = "Sum PF Photon Et (Higher Pt threshold), #Delta R = 0.3";
170  main_titles[31] = "Sum PF Charged Particles Pt not from PV (for Pu corrections), #Delta R = 0.3";
171 
172  //-----Titles of the plots-----------
173  main_titles[32] = "Sum PF Charged Hadron Pt, #Delta R = 0.4";
174  main_titles[33] = "Sum PF Neutral Hadron Pt, #Delta R = 0.4";
175  main_titles[34] = "Sum PF Photon Et, #Delta R = 0.4";
176  main_titles[35] = "Sum PF Neutral Hadron Pt (Higher Pt threshold), #Delta R = 0.4";
177  main_titles[36] = "Sum PF Photon Et (Higher Pt threshold), #Delta R = 0.4";
178  main_titles[37] = "Sum PF Charged Particles Pt not from PV (for Pu corrections), #Delta R = 0.4";
179 
180  main_titles[38] = "Relative PF Isolation, #Delta R = 0.3";
181  main_titles[39] = "Relative PF Isolation, #Delta R = 0.4";
182 
183  main_titles[40] = "Relative PF Isolation (Higher Pt threshold), #Delta R = 0.3";
184  main_titles[41] = "Relative PF Isolation (Higher Pt threshold), #Delta R = 0.4";
185 
186  main_titles[42] = "Sum DR Isolation Profile for Charged Hadron, #Delta R = 0.4";
187 
188  main_titles[43] = "Sum DR Isolation Profile for Neutral Hadron, #Delta R = 0.4";
189 
190  main_titles[44] = "Sum DR Isolation Profile for Photon, #Delta R = 0.4";
191 
192  main_titles[45] = "Mean DR Isolation Profile for Charged Hadron, #Delta R = 0.4";
193 
194  main_titles[46] = "Mean DR Isolation Profile for Neutral Hadron, #Delta R = 0.4";
195 
196  main_titles[47] = "Mean DR Isolation Profile for Photon, #Delta R = 0.4";
197 
198 #ifdef DEBUG
199  cout << "InitStatistics(): main titles 1D DONE " << endl;
200 #endif
201  titles_2D[0] = "Total Tracker Momentum, #Delta R = 0.3";
202  titles_2D[1] = "Total EM Cal Energy, #Delta R = 0.3";
203  titles_2D[2] = "Total Had Cal Energy, #Delta R = 0.3";
204  titles_2D[3] = "Total HO Cal Energy, #Delta R = 0.3";
205  titles_2D[4] = "Sum PF Charged Hadron Pt, #Delta R = 0.4";
206  titles_2D[5] = "Sum PF Neutral Hadron Pt, #Delta R = 0.4";
207  titles_2D[6] = "Sum PF Photon Et, #Delta R = 0.4";
208  titles_2D[7] = "Sum PF Charged Pt Not from PV, #Delta R = 0.4";
209  titles_2D[8] = "Relative Detector-Based Isolation, #Delta R = 0.4";
210  titles_2D[9] = "Relative PF Isolation, #Delta R = 0.4";
211 
212  main_titles_NVtxs[0] = "Sum PF Neutral Hadron Pt, #DeltaR = 0.4 ( 0 < N_{Vtx} < 15)";
213  main_titles_NVtxs[1] = "Sum PF Neutral Hadron Pt, #DeltaR = 0.4 (15 < N_{Vtx} < 30)";
214  main_titles_NVtxs[2] = "Sum PF Neutral Hadron Pt, #DeltaR = 0.4 (30 < N_{Vtx})";
215  main_titles_NVtxs[3] = "Sum PF Photon Et, #DeltaR = 0.4 ( 0 < N_{Vtx} < 15)";
216  main_titles_NVtxs[4] = "Sum PF Photon Et, #DeltaR = 0.4 (15 < N_{Vtx} < 30)";
217  main_titles_NVtxs[5] = "Sum PF Photon Et, #DeltaR = 0.4 (30 < N_{Vtx})";
218 
219 #ifdef DEBUG
220  cout << "InitStatistics(): main titles 2D DONE " << endl;
221 #endif
222 
223  //------Titles on the X or Y axis------------
224  axis_titles[0] = "#Sigma p_{T} (GeV)";
225  axis_titles[1] = "#Sigma E_{T}^{EM} (GeV)";
226  axis_titles[2] = "#Sigma E_{T}^{Had} (GeV)";
227  axis_titles[3] = "#Sigma E_{T}^{HO} (GeV)";
228  axis_titles[4] = "N_{Tracks}";
229  axis_titles[5] = "N_{Jets}";
230  axis_titles[6] = "#Sigma p_{T,veto} (GeV)";
231  axis_titles[7] = "#Sigma E_{T,veto}^{EM} (GeV)";
232  axis_titles[8] = "#Sigma E_{T,veto}^{Had} (GeV)";
233  axis_titles[9] = "#Sigma E_{T,veto}^{HO} (GeV)";
234  axis_titles[10] = "#Sigma p_{T} / N_{Tracks} (GeV)";
235  axis_titles[11] = "(1.5) X #Sigma E_{T}^{EM} + #Sigma E_{T}^{Had}";
236 
237  axis_titles[12] = "#Sigma p_{T} (GeV)";
238  axis_titles[13] = "#Sigma E_{T}^{EM} (GeV)";
239  axis_titles[14] = "#Sigma E_{T}^{Had} (GeV)";
240  axis_titles[15] = "#Sigma E_{T}^{HO} (GeV)";
241  axis_titles[16] = "N_{Tracks}";
242  axis_titles[17] = "N_{Jets}";
243  axis_titles[18] = "#Sigma p_{T,veto} (GeV)";
244  axis_titles[19] = "#Sigma E_{T,veto}^{EM} (GeV)";
245  axis_titles[20] = "#Sigma E_{T,veto}^{Had} (GeV)";
246  axis_titles[21] = "#Sigma E_{T,veto}^{HO} (GeV)";
247  axis_titles[22] = "#Sigma p_{T} / N_{Tracks} (GeV)";
248  axis_titles[23] = "(1.5) X #Sigma E_{T}^{EM} + #Sigma E_{T}^{Had}";
249 
250  axis_titles[24] = "(#Sigma Tk p_{T} + #Sigma ECAL p_{T} + #Sigma HCAL p_{T})/ Mu p_{T} (GeV)";
251  axis_titles[25] = "(#Sigma Tk p_{T} + #Sigma ECAL p_{T} + #Sigma HCAL p_{T})/ Mu p_{T} (GeV)";
252 
253  axis_titles[26] = "#Sigma PFCharged p_{T}";
254  axis_titles[27] = "#Sigma PFNeutral p_{T}";
255  axis_titles[28] = "#Sigma PFPhoton p_{T}";
256  axis_titles[29] = "#Sigma PFNeutral p_{T}";
257  axis_titles[30] = "#Sigma PFPhoton p_{T}";
258  axis_titles[31] = "#Sigma PFCharged p_{T}";
259 
260  axis_titles[32] = "#Sigma PFCharged p_{T}";
261  axis_titles[33] = "#Sigma PFNeutral p_{T}";
262  axis_titles[34] = "#Sigma PFPhoton p_{T}";
263  axis_titles[35] = "#Sigma PFNeutral p_{T}";
264  axis_titles[36] = "#Sigma PFPhoton p_{T}";
265  axis_titles[37] = "#Sigma PFCharged p_{T}";
266 
267  axis_titles[38] = "(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
268  axis_titles[39] = "(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
269  axis_titles[40] = "(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
270  axis_titles[41] = "(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
271 
272  axis_titles[42] = "#Sigma DR PFCharged";
273  axis_titles[43] = "#Sigma DR PFNeutral";
274  axis_titles[44] = "#Sigma DR PFPhoton";
275 
276  axis_titles[45] = "Mean DR PFCharged";
277  axis_titles[46] = "Mean DR PFNeutral";
278  axis_titles[47] = "Mean DR PFPhoton";
279 
280  axis_titles_NVtxs[0] = "#Sigma PFNeutral p_{T}";
281  axis_titles_NVtxs[1] = "#Sigma PFNeutral p_{T}";
282  axis_titles_NVtxs[2] = "#Sigma PFNeutral p_{T}";
283  axis_titles_NVtxs[3] = "#Sigma PFPhoton p_{T}";
284  axis_titles_NVtxs[4] = "#Sigma PFPhoton p_{T}";
285  axis_titles_NVtxs[5] = "#Sigma PFPhoton p_{T}";
286 
287 #ifdef DEBUG
288  cout << "InitStatistics(): main titles 1D DONE " << endl;
289 #endif
290 
291  //-----------Names given for the root file----------
292  names[0] = "sumPt_R03";
293  names[1] = "emEt_R03";
294  names[2] = "hadEt_R03";
295  names[3] = "hoEt_R03";
296  names[4] = "nTracks_R03";
297  names[5] = "nJets_R03";
298  names[6] = "trackerVetoPt_R03";
299  names[7] = "emVetoEt_R03";
300  names[8] = "hadVetoEt_R03";
301  names[9] = "hoVetoEt_R03";
302  names[10] = "avgPt_R03";
303  names[11] = "weightedEt_R03";
304 
305  names[12] = "sumPt_R05";
306  names[13] = "emEt_R05";
307  names[14] = "hadEt_R05";
308  names[15] = "hoEt_R05";
309  names[16] = "nTracks_R05";
310  names[17] = "nJets_R05";
311  names[18] = "trackerVetoPt_R05";
312  names[19] = "emVetoEt_R05";
313  names[20] = "hadVetoEt_R05";
314  names[21] = "hoVetoEt_R05";
315  names[22] = "avgPt_R05";
316  names[23] = "weightedEt_R05";
317 
318  names[24] = "relDetIso_R03";
319  names[25] = "relDetIso_R05";
320 
321  names[26] = "pfChargedPt_R03";
322  names[27] = "pfNeutralPt_R03";
323  names[28] = "pfPhotonPt_R03";
324  names[29] = "pfNeutralPt_HT_R03";
325  names[30] = "pfPhotonPt_HT_R03";
326  names[31] = "pfChargedPt_PU_R03";
327 
328  names[32] = "pfChargedPt_R04";
329  names[33] = "pfNeutralPt_R04";
330  names[34] = "pfPhotonPt_R04";
331  names[35] = "pfNeutralPt_HT_R04";
332  names[36] = "pfPhotonPt_HT_R04";
333  names[37] = "pfChargedPt_PU_R04";
334 
335  names[38] = "relPFIso_R03";
336  names[39] = "relPFIso_R04";
337 
338  names[40] = "relPFIso_HT_R03";
339  names[41] = "relPFIso_HT_R04";
340 
341  names[42] = "SumDR_PFCharged_R04";
342  names[43] = "SumDR_PFNeutral_R04";
343  names[44] = "SumDR_PFPhoton_R04";
344 
345  names[45] = "MeanDR_PFCharged_R04";
346  names[46] = "MeanDR_PFNeutral_R04";
347  names[47] = "MeanDR_PFPhoton_R04";
348 
349 #ifdef DEBUG
350  cout << "InitStatistics(): names 1D DONE " << endl;
351 #endif
352 
353  names_2D[0] = "SumPt_R03";
354  names_2D[1] = "emEt_R03";
355  names_2D[2] = "hadEt_R03";
356  names_2D[3] = "hoEt_R03";
357  names_2D[4] = "pfChargedPt_R04";
358  names_2D[5] = "pfNeutralPt_R04";
359  names_2D[6] = "pfPhotonPt_R04";
360  names_2D[7] = "pfChargedPUPt_R04";
361  names_2D[8] = "relDetIso_R03";
362  names_2D[9] = "relPFIso_R04";
363 
364 #ifdef DEBUG
365  cout << "InitStatistics(): names 2D DONE " << endl;
366 #endif
367 
368  names_NVtxs[0] = "pfNeutralPt_R04_PV0to15";
369  names_NVtxs[1] = "pfNeutralPt_R04_PV15to30";
370  names_NVtxs[2] = "pfNeutralPt_R04_PV30toInf";
371  names_NVtxs[3] = "pfPhotonPt_R04_PV0to15";
372  names_NVtxs[4] = "pfPhotonPt_R04_PV15to30";
373  names_NVtxs[5] = "pfPhotonPt_R04_PV30toInf";
374 
375  //----------Parameters for binning of histograms---------
376  //param[var][0] is the number of bins
377  //param[var][1] is the low edge of the low bin
378  //param[var][2] is the high edge of the high bin
379  //
380  // maximum value------,
381  // |
382  // V
383  param[0][0] = (int)(20.0 / S_BIN_WIDTH);
384  param[0][1] = 0.0;
385  param[0][2] = param[0][0] * S_BIN_WIDTH;
386  param[1][0] = (int)(20.0 / S_BIN_WIDTH);
387  param[1][1] = 0.0;
388  param[1][2] = param[1][0] * S_BIN_WIDTH;
389  param[2][0] = (int)(20.0 / S_BIN_WIDTH);
390  param[2][1] = 0.0;
391  param[2][2] = param[2][0] * S_BIN_WIDTH;
392  param[3][0] = 20;
393  param[3][1] = 0.0;
394  param[3][2] = 2.0;
395  param[4][0] = 16;
396  param[4][1] = -0.5;
397  param[4][2] = param[4][0] - 0.5;
398  param[5][0] = 4;
399  param[5][1] = -0.5;
400  param[5][2] = param[5][0] - 0.5;
401  param[6][0] = (int)(40.0 / S_BIN_WIDTH);
402  param[6][1] = 0.0;
403  param[6][2] = param[6][0] * S_BIN_WIDTH;
404  param[7][0] = 20;
405  param[7][1] = 0.0;
406  param[7][2] = 10.0;
407  param[8][0] = (int)(20.0 / S_BIN_WIDTH);
408  param[8][1] = 0.0;
409  param[8][2] = param[8][0] * S_BIN_WIDTH;
410  param[9][0] = 20;
411  param[9][1] = 0.0;
412  param[9][2] = 5.0;
413  param[10][0] = (int)(15.0 / S_BIN_WIDTH);
414  param[10][1] = 0.0;
415  param[10][2] = param[10][0] * S_BIN_WIDTH;
416  param[11][0] = (int)(20.0 / S_BIN_WIDTH);
417  param[11][1] = 0.0;
418  param[11][2] = param[11][0] * S_BIN_WIDTH;
419 
420  param[12][0] = (int)(20.0 / S_BIN_WIDTH);
421  param[12][1] = 0.0;
422  param[12][2] = param[12][0] * S_BIN_WIDTH;
423  param[13][0] = (int)(20.0 / S_BIN_WIDTH);
424  param[13][1] = 0.0;
425  param[13][2] = param[13][0] * S_BIN_WIDTH;
426  param[14][0] = (int)(20.0 / S_BIN_WIDTH);
427  param[14][1] = 0.0;
428  param[14][2] = param[14][0] * S_BIN_WIDTH;
429  param[15][0] = 20;
430  param[15][1] = 0.0;
431  param[15][2] = 2.0;
432  param[16][0] = 16;
433  param[16][1] = -0.5;
434  param[16][2] = param[16][0] - 0.5;
435  param[17][0] = 4;
436  param[17][1] = -0.5;
437  param[17][2] = param[17][0] - 0.5;
438  param[18][0] = (int)(40.0 / S_BIN_WIDTH);
439  param[18][1] = 0.0;
440  param[18][2] = param[18][0] * S_BIN_WIDTH;
441  param[19][0] = 20;
442  param[19][1] = 0.0;
443  param[19][2] = 10.0;
444  param[20][0] = (int)(20.0 / S_BIN_WIDTH);
445  param[20][1] = 0.0;
446  param[20][2] = param[20][0] * S_BIN_WIDTH;
447  param[21][0] = 20;
448  param[21][1] = 0.0;
449  param[21][2] = 5.0;
450  param[22][0] = (int)(15.0 / S_BIN_WIDTH);
451  param[22][1] = 0.0;
452  param[22][2] = param[22][0] * S_BIN_WIDTH;
453  param[23][0] = (int)(20.0 / S_BIN_WIDTH);
454  param[23][1] = 0.0;
455  param[23][2] = param[23][0] * S_BIN_WIDTH;
456 
457  param[24][0] = 50;
458  param[24][1] = 0.0;
459  param[24][2] = 1.0;
460  param[25][0] = 50;
461  param[25][1] = 0.0;
462  param[25][2] = 1.0;
463 
464  param[26][0] = (int)(20.0 / S_BIN_WIDTH);
465  param[26][1] = 0.0;
466  param[26][2] = param[26][0] * S_BIN_WIDTH;
467  param[27][0] = (int)(20.0 / S_BIN_WIDTH);
468  param[27][1] = 0.0;
469  param[27][2] = param[27][0] * S_BIN_WIDTH;
470  param[28][0] = (int)(20.0 / S_BIN_WIDTH);
471  param[28][1] = 0.0;
472  param[28][2] = param[28][0] * S_BIN_WIDTH;
473  param[29][0] = (int)(20.0 / S_BIN_WIDTH);
474  param[29][1] = 0.0;
475  param[29][2] = param[29][0] * S_BIN_WIDTH;
476  param[30][0] = (int)(20.0 / S_BIN_WIDTH);
477  param[30][1] = 0.0;
478  param[30][2] = param[30][0] * S_BIN_WIDTH;
479  param[31][0] = (int)(20.0 / S_BIN_WIDTH);
480  param[31][1] = 0.0;
481  param[31][2] = param[31][0] * S_BIN_WIDTH;
482 
483  param[32][0] = (int)(20.0 / S_BIN_WIDTH);
484  param[32][1] = 0.0;
485  param[32][2] = param[32][0] * S_BIN_WIDTH;
486  param[33][0] = (int)(20.0 / S_BIN_WIDTH);
487  param[33][1] = 0.0;
488  param[33][2] = param[33][0] * S_BIN_WIDTH;
489  param[34][0] = (int)(20.0 / S_BIN_WIDTH);
490  param[34][1] = 0.0;
491  param[34][2] = param[34][0] * S_BIN_WIDTH;
492  param[35][0] = (int)(20.0 / S_BIN_WIDTH);
493  param[35][1] = 0.0;
494  param[35][2] = param[35][0] * S_BIN_WIDTH;
495  param[36][0] = (int)(20.0 / S_BIN_WIDTH);
496  param[36][1] = 0.0;
497  param[36][2] = param[36][0] * S_BIN_WIDTH;
498  param[37][0] = (int)(20.0 / S_BIN_WIDTH);
499  param[37][1] = 0.0;
500  param[37][2] = param[37][0] * S_BIN_WIDTH;
501 
502  param[38][0] = 50;
503  param[38][1] = 0.0;
504  param[38][2] = 1.0;
505  param[39][0] = 50;
506  param[39][1] = 0.0;
507  param[39][2] = 1.0;
508 
509  param[40][0] = 50;
510  param[40][1] = 0.0;
511  param[40][2] = 1.0;
512  param[41][0] = 50;
513  param[41][1] = 0.0;
514  param[41][2] = 1.0;
515 
516  param[42][0] = 50;
517  param[42][1] = 0.0;
518  param[42][2] = 5;
519  param[43][0] = 50;
520  param[43][1] = 0.0;
521  param[43][2] = 5;
522  param[44][0] = 50;
523  param[44][1] = 0.0;
524  param[44][2] = 5;
525 
526  param[45][0] = 50;
527  param[45][1] = 0.0;
528  param[45][2] = 0.4;
529  param[46][0] = 50;
530  param[46][1] = 0.0;
531  param[46][2] = 0.4;
532  param[47][0] = 50;
533  param[47][1] = 0.0;
534  param[47][2] = 0.4;
535 
536  //--------------Is the variable continuous (i.e. non-integer)?-------------
537  //---------(Log binning will only be used for continuous variables)--------
538  isContinuous[0] = 1;
539  isContinuous[1] = 1;
540  isContinuous[2] = 1;
541  isContinuous[3] = 1;
542  isContinuous[4] = 0;
543  isContinuous[5] = 0;
544  isContinuous[6] = 1;
545  isContinuous[7] = 1;
546  isContinuous[8] = 1;
547  isContinuous[9] = 1;
548  isContinuous[10] = 1;
549  isContinuous[11] = 1;
550 
551  isContinuous[12] = 1;
552  isContinuous[13] = 1;
553  isContinuous[14] = 1;
554  isContinuous[15] = 1;
555  isContinuous[16] = 0;
556  isContinuous[17] = 0;
557  isContinuous[18] = 1;
558  isContinuous[19] = 1;
559  isContinuous[20] = 1;
560  isContinuous[21] = 1;
561  isContinuous[22] = 1;
562  isContinuous[23] = 1;
563 
564  isContinuous[24] = 1;
565  isContinuous[25] = 1;
566  isContinuous[26] = 1;
567  isContinuous[27] = 1;
568  isContinuous[28] = 1;
569  isContinuous[29] = 1;
570  isContinuous[30] = 1;
571  isContinuous[31] = 1;
572  isContinuous[32] = 1;
573  isContinuous[33] = 1;
574  isContinuous[34] = 1;
575  isContinuous[35] = 1;
576  isContinuous[36] = 1;
577  isContinuous[37] = 1;
578  isContinuous[38] = 1;
579  isContinuous[39] = 1;
580  isContinuous[40] = 1;
581  isContinuous[41] = 1;
582  isContinuous[42] = 1;
583  isContinuous[43] = 1;
584  isContinuous[44] = 1;
585  isContinuous[45] = 1;
586  isContinuous[46] = 1;
587  isContinuous[47] = 1;
588 
589 #ifdef DEBUG
590  cout << "InitStatistics(): DONE " << endl;
591 #endif
592 }

References gather_cfg::cout, createfilelist::int, and names.

◆ NormalizeHistos()

void MuonIsolationDQM::NormalizeHistos ( )
private

Definition at line 797 of file MuonIsolationDQM.cc.

797  {
798  for (int var = 0; var < NUM_VARS; var++) {
799  double entries = GetTH1FromMonitorElement(h_1D[var])->GetEntries();
800  GetTH1FromMonitorElement(h_1D[var])->Scale(1. / entries);
801  }
802 }

References trigObjTnPSource_cfi::var.

◆ RecordData()

void MuonIsolationDQM::RecordData ( const reco::Muon &  muon)
private

Definition at line 659 of file MuonIsolationDQM.cc.

659  {
660 #ifdef DEBUG
661  std::cout << "RecordData()" << endl;
662 #endif
663  float MuPt = muon.pt();
664 
665  theData[0] = muon.isolationR03().sumPt;
666  theData[1] = muon.isolationR03().emEt;
667  theData[2] = muon.isolationR03().hadEt;
668  theData[3] = muon.isolationR03().hoEt;
669 
670  theData[4] = muon.isolationR03().nTracks;
671  theData[5] = muon.isolationR03().nJets;
672  theData[6] = muon.isolationR03().trackerVetoPt;
673  theData[7] = muon.isolationR03().emVetoEt;
674  theData[8] = muon.isolationR03().hadVetoEt;
675  theData[9] = muon.isolationR03().hoVetoEt;
676 
677  // make sure nTracks != 0 before filling this one
678  if (theData[4] != 0)
679  theData[10] = (double)theData[0] / (double)theData[4];
680  else
681  theData[10] = -99;
682 
683  theData[11] = 1.5 * theData[1] + theData[2];
684 
685  theData[12] = muon.isolationR05().sumPt;
686  theData[13] = muon.isolationR05().emEt;
687  theData[14] = muon.isolationR05().hadEt;
688  theData[15] = muon.isolationR05().hoEt;
689 
690  theData[16] = muon.isolationR05().nTracks;
691  theData[17] = muon.isolationR05().nJets;
692  theData[18] = muon.isolationR05().trackerVetoPt;
693  theData[19] = muon.isolationR05().emVetoEt;
694  theData[20] = muon.isolationR05().hadVetoEt;
695  theData[21] = muon.isolationR05().hoVetoEt;
696 
697  // make sure nTracks != 0 before filling this one
698  if (theData[16] != 0)
699  theData[22] = (double)theData[12] / (double)theData[16];
700  else
701  theData[22] = -99;
702 
703  theData[23] = 1.5 * theData[13] + theData[14];
704 
705  theData[24] = (theData[0] + theData[1] + theData[2]) / MuPt;
706  theData[25] = (theData[12] + theData[13] + theData[14]) / MuPt;
707 
708  theData[26] = muon.pfIsolationR03().sumChargedHadronPt;
709  theData[27] = muon.pfIsolationR03().sumNeutralHadronEt;
710  theData[28] = muon.pfIsolationR03().sumPhotonEt;
711  theData[29] = muon.pfIsolationR03().sumNeutralHadronEtHighThreshold;
712  theData[30] = muon.pfIsolationR03().sumPhotonEtHighThreshold;
713  theData[31] = muon.pfIsolationR03().sumPUPt;
714 
715  theData[32] = muon.pfIsolationR04().sumChargedHadronPt;
716  theData[33] = muon.pfIsolationR04().sumNeutralHadronEt;
717  theData[34] = muon.pfIsolationR04().sumPhotonEt;
718  theData[35] = muon.pfIsolationR04().sumNeutralHadronEtHighThreshold;
719  theData[36] = muon.pfIsolationR04().sumPhotonEtHighThreshold;
720  theData[37] = muon.pfIsolationR04().sumPUPt;
721 
722  theData[38] = (theData[26] + theData[27] + theData[28]) / MuPt;
723  theData[39] = (theData[32] + theData[33] + theData[34]) / MuPt;
724 
725  theData[40] = (theData[26] + theData[29] + theData[30]) / MuPt;
726  theData[41] = (theData[32] + theData[35] + theData[36]) / MuPt;
727 
728  theData[42] = muon.pfSumDRIsoProfileR04().sumChargedHadronPt;
729  theData[43] = muon.pfSumDRIsoProfileR04().sumNeutralHadronEt;
730  theData[44] = muon.pfSumDRIsoProfileR04().sumPhotonEt;
731  theData[45] = muon.pfMeanDRIsoProfileR04().sumChargedHadronPt;
732  theData[46] = muon.pfMeanDRIsoProfileR04().sumNeutralHadronEt;
733  theData[47] = muon.pfMeanDRIsoProfileR04().sumPhotonEt;
734 
735  //--------------Filling the 2D Histos Data -------- //
736  theData2D[0] = muon.isolationR03().sumPt;
737  theData2D[1] = muon.isolationR03().emEt;
738  theData2D[2] = muon.isolationR03().hadEt;
739  theData2D[3] = muon.isolationR03().hoEt;
740 
741  theData2D[4] = muon.pfIsolationR04().sumChargedHadronPt;
742  theData2D[5] = muon.pfIsolationR04().sumNeutralHadronEt;
743  theData2D[6] = muon.pfIsolationR04().sumPhotonEt;
744  theData2D[7] = muon.pfIsolationR04().sumPUPt;
745 
746  theData2D[8] = theData2D[0] + theData2D[1] + theData2D[2] + theData2D[3] / MuPt; //Det RelIso;
747  theData2D[9] = theData2D[4] + theData2D[5] + theData2D[6] / MuPt; //PF RelIso;
748 
749  //-----------Filling the NVTX 1D HISTOS DATA ------------- //
750  theDataNVtx[0] = muon.pfIsolationR04().sumNeutralHadronEt;
751  theDataNVtx[1] = theDataNVtx[0];
752  theDataNVtx[2] = theDataNVtx[0];
753 
754  theDataNVtx[3] = muon.pfIsolationR04().sumPhotonEt;
755  theDataNVtx[4] = theDataNVtx[3];
756  theDataNVtx[5] = theDataNVtx[3];
757 }

References gather_cfg::cout.

Member Data Documentation

◆ axis_titles

std::vector<std::string> MuonIsolationDQM::axis_titles
private

Definition at line 118 of file MuonIsolationDQM.h.

◆ axis_titles_NVtxs

std::vector<std::string> MuonIsolationDQM::axis_titles_NVtxs
private

Definition at line 128 of file MuonIsolationDQM.h.

◆ dirName

std::string MuonIsolationDQM::dirName
private

Definition at line 97 of file MuonIsolationDQM.h.

Referenced by addOnTests.testit::run(), analyzer.Analyzer::write(), and core.JSONAnalyzer.JSONAnalyzer::write().

◆ h_1D

std::vector<MonitorElement*> MuonIsolationDQM::h_1D
private

Definition at line 138 of file MuonIsolationDQM.h.

◆ h_1D_NVTX

std::vector<MonitorElement*> MuonIsolationDQM::h_1D_NVTX
private

Definition at line 140 of file MuonIsolationDQM.h.

◆ h_2D

std::vector<MonitorElement*> MuonIsolationDQM::h_2D
private

Definition at line 139 of file MuonIsolationDQM.h.

◆ h_nMuons

MonitorElement* MuonIsolationDQM::h_nMuons
private

Definition at line 137 of file MuonIsolationDQM.h.

◆ isContinuous

std::vector<int> MuonIsolationDQM::isContinuous
private

Definition at line 121 of file MuonIsolationDQM.h.

◆ L_BIN_WIDTH

double MuonIsolationDQM::L_BIN_WIDTH
private

Definition at line 104 of file MuonIsolationDQM.h.

◆ LOG_BINNING_ENABLED

int MuonIsolationDQM::LOG_BINNING_ENABLED
private

Definition at line 106 of file MuonIsolationDQM.h.

◆ LOG_BINNING_RATIO

double MuonIsolationDQM::LOG_BINNING_RATIO
private

Definition at line 108 of file MuonIsolationDQM.h.

◆ main_titles

std::vector<std::string> MuonIsolationDQM::main_titles
private

Definition at line 117 of file MuonIsolationDQM.h.

◆ main_titles_NVtxs

std::vector<std::string> MuonIsolationDQM::main_titles_NVtxs
private

Definition at line 126 of file MuonIsolationDQM.h.

◆ MuIsoDepRef

const typedef reco::IsoDeposit MuonIsolationDQM::MuIsoDepRef
private

Definition at line 64 of file MuonIsolationDQM.h.

◆ names

std::vector<std::string> MuonIsolationDQM::names
private

Definition at line 119 of file MuonIsolationDQM.h.

◆ names_2D

std::vector<std::string> MuonIsolationDQM::names_2D
private

Definition at line 124 of file MuonIsolationDQM.h.

◆ names_NVtxs

std::vector<std::string> MuonIsolationDQM::names_NVtxs
private

Definition at line 127 of file MuonIsolationDQM.h.

◆ nEvents

int MuonIsolationDQM::nEvents
private

Definition at line 145 of file MuonIsolationDQM.h.

Referenced by looper.Looper::loop().

◆ nGLBMuons

int MuonIsolationDQM::nGLBMuons
private

Definition at line 147 of file MuonIsolationDQM.h.

◆ nSTAMuons

int MuonIsolationDQM::nSTAMuons
private

Definition at line 146 of file MuonIsolationDQM.h.

◆ nTRKMuons

int MuonIsolationDQM::nTRKMuons
private

Definition at line 148 of file MuonIsolationDQM.h.

◆ NUM_LOG_BINS

int MuonIsolationDQM::NUM_LOG_BINS
private

Definition at line 107 of file MuonIsolationDQM.h.

◆ NUM_VARS

const int MuonIsolationDQM::NUM_VARS = 48
staticprivate

Definition at line 100 of file MuonIsolationDQM.h.

◆ NUM_VARS_2D

const int MuonIsolationDQM::NUM_VARS_2D = 10
staticprivate

Definition at line 101 of file MuonIsolationDQM.h.

◆ NUM_VARS_NVTX

const int MuonIsolationDQM::NUM_VARS_NVTX = 6
staticprivate

Definition at line 102 of file MuonIsolationDQM.h.

◆ param

std::vector<std::vector<double> > MuonIsolationDQM::param
private

Definition at line 120 of file MuonIsolationDQM.h.

◆ requireGLBMuon

bool MuonIsolationDQM::requireGLBMuon
private

Definition at line 109 of file MuonIsolationDQM.h.

◆ requireSTAMuon

bool MuonIsolationDQM::requireSTAMuon
private

Definition at line 110 of file MuonIsolationDQM.h.

◆ requireTRKMuon

bool MuonIsolationDQM::requireTRKMuon
private

Definition at line 111 of file MuonIsolationDQM.h.

◆ rootfilename

std::string MuonIsolationDQM::rootfilename
private

Definition at line 95 of file MuonIsolationDQM.h.

◆ S_BIN_WIDTH

double MuonIsolationDQM::S_BIN_WIDTH
private

Definition at line 105 of file MuonIsolationDQM.h.

◆ theData

double MuonIsolationDQM::theData[NUM_VARS]
private

Definition at line 132 of file MuonIsolationDQM.h.

◆ theData2D

double MuonIsolationDQM::theData2D[NUM_VARS_2D]
private

Definition at line 133 of file MuonIsolationDQM.h.

◆ theDataNVtx

double MuonIsolationDQM::theDataNVtx[NUM_VARS_NVTX]
private

Definition at line 134 of file MuonIsolationDQM.h.

◆ theMuonCollectionLabel_

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

Definition at line 92 of file MuonIsolationDQM.h.

◆ theVertexCollectionLabel_

edm::EDGetTokenT<reco::VertexCollection> MuonIsolationDQM::theVertexCollectionLabel_
private

Definition at line 91 of file MuonIsolationDQM.h.

◆ title_cone

std::string MuonIsolationDQM::title_cone
private

Definition at line 114 of file MuonIsolationDQM.h.

◆ title_sam

std::string MuonIsolationDQM::title_sam
private

Definition at line 113 of file MuonIsolationDQM.h.

◆ titles_2D

std::vector<std::string> MuonIsolationDQM::titles_2D
private

Definition at line 123 of file MuonIsolationDQM.h.

PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
MuonIsolationDQM::S_BIN_WIDTH
double S_BIN_WIDTH
Definition: MuonIsolationDQM.h:105
MuonIsolationDQM::title_sam
std::string title_sam
Definition: MuonIsolationDQM.h:113
edm::Handle::product
T const * product() const
Definition: Handle.h:70
MuonIsolationDQM::title_cone
std::string title_cone
Definition: MuonIsolationDQM.h:114
MuonIsolationDQM::titles_2D
std::vector< std::string > titles_2D
Definition: MuonIsolationDQM.h:123
muon
Definition: MuonCocktails.h:17
MuonIsolationDQM::isContinuous
std::vector< int > isContinuous
Definition: MuonIsolationDQM.h:121
MuonIsolationDQM::FillHistos
void FillHistos(int)
Definition: MuonIsolationDQM.cc:804
MuonIsolationDQM::theDataNVtx
double theDataNVtx[NUM_VARS_NVTX]
Definition: MuonIsolationDQM.h:134
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
gather_cfg.cout
cout
Definition: gather_cfg.py:144
edm::LogInfo
Definition: MessageLogger.h:254
MuonIsolationDQM::param
std::vector< std::vector< double > > param
Definition: MuonIsolationDQM.h:120
MuonIsolationDQM::axis_titles
std::vector< std::string > axis_titles
Definition: MuonIsolationDQM.h:118
MuonIsolationDQM::axis_titles_NVtxs
std::vector< std::string > axis_titles_NVtxs
Definition: MuonIsolationDQM.h:128
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
findQualityFiles.v
v
Definition: findQualityFiles.py:179
HLT_2018_cff.muon
muon
Definition: HLT_2018_cff.py:10349
edm::Handle
Definition: AssociativeIterator.h:50
nanoDQM_cfi.PV
PV
Definition: nanoDQM_cfi.py:440
trigObjTnPSource_cfi.var
var
Definition: trigObjTnPSource_cfi.py:21
MuonIsolationDQM::NUM_VARS
static const int NUM_VARS
Definition: MuonIsolationDQM.h:100
MuonIsolationDQM::requireGLBMuon
bool requireGLBMuon
Definition: MuonIsolationDQM.h:109
MuonIsolationDQM::theData2D
double theData2D[NUM_VARS_2D]
Definition: MuonIsolationDQM.h:133
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
MuonIsolationDQM::theMuonCollectionLabel_
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
Definition: MuonIsolationDQM.h:92
MuonIsolationDQM::theData
double theData[NUM_VARS]
Definition: MuonIsolationDQM.h:132
MuonIsolationDQM::nTRKMuons
int nTRKMuons
Definition: MuonIsolationDQM.h:148
MuonIsolationDQM::theVertexCollectionLabel_
edm::EDGetTokenT< reco::VertexCollection > theVertexCollectionLabel_
Definition: MuonIsolationDQM.h:91
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MuonIsolationDQM::NUM_VARS_NVTX
static const int NUM_VARS_NVTX
Definition: MuonIsolationDQM.h:102
MuonIsolationDQM::L_BIN_WIDTH
double L_BIN_WIDTH
Definition: MuonIsolationDQM.h:104
MuonIsolationDQM::main_titles_NVtxs
std::vector< std::string > main_titles_NVtxs
Definition: MuonIsolationDQM.h:126
MuonIsolationDQM::XAXIS
Definition: MuonIsolationDQM.h:151
MuonIsolationDQM::GetTH1FromMonitorElement
TH1 * GetTH1FromMonitorElement(MonitorElement *me)
Definition: MuonIsolationDQM.cc:844
bphysicsOniaDQM_cfi.vertex
vertex
Definition: bphysicsOniaDQM_cfi.py:7
MuonIsolationDQM::h_2D
std::vector< MonitorElement * > h_2D
Definition: MuonIsolationDQM.h:139
createfilelist.int
int
Definition: createfilelist.py:10
iEvent
int iEvent
Definition: GenABIO.cc:224
MuonIsolationDQM::nGLBMuons
int nGLBMuons
Definition: MuonIsolationDQM.h:147
MuonIsolationDQM::main_titles
std::vector< std::string > main_titles
Definition: MuonIsolationDQM.h:117
MuonIsolationDQM::dirName
std::string dirName
Definition: MuonIsolationDQM.h:97
MuonIsolationDQM::h_nMuons
MonitorElement * h_nMuons
Definition: MuonIsolationDQM.h:137
MuonIsolationDQM::YAXIS
Definition: MuonIsolationDQM.h:151
MuonIsolationDQM::FillNVtxHistos
void FillNVtxHistos(int)
Definition: MuonIsolationDQM.cc:829
MuonIsolationDQM::LOG_BINNING_RATIO
double LOG_BINNING_RATIO
Definition: MuonIsolationDQM.h:108
MuonIsolationDQM::nSTAMuons
int nSTAMuons
Definition: MuonIsolationDQM.h:146
MuonIsolationDQM::InitStatics
void InitStatics()
Definition: MuonIsolationDQM.cc:97
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
MuonIsolationDQM::h_1D
std::vector< MonitorElement * > h_1D
Definition: MuonIsolationDQM.h:138
MuonIsolationDQM::names
std::vector< std::string > names
Definition: MuonIsolationDQM.h:119
MuonIsolationDQM::requireSTAMuon
bool requireSTAMuon
Definition: MuonIsolationDQM.h:110
MuonIsolationDQM::NOAXIS
Definition: MuonIsolationDQM.h:151
MuonIsolationDQM::ZAXIS
Definition: MuonIsolationDQM.h:151
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
MuonIsolationDQM::h_1D_NVTX
std::vector< MonitorElement * > h_1D_NVTX
Definition: MuonIsolationDQM.h:140
MuonIsolationDQM::NUM_LOG_BINS
int NUM_LOG_BINS
Definition: MuonIsolationDQM.h:107
MuonIsolationDQM::names_2D
std::vector< std::string > names_2D
Definition: MuonIsolationDQM.h:124
MuonIsolationDQM::requireTRKMuon
bool requireTRKMuon
Definition: MuonIsolationDQM.h:111
hlt_dqm_clientPB-live_cfg.me
me
Definition: hlt_dqm_clientPB-live_cfg.py:56
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
MuonIsolationDQM::RecordData
void RecordData(const reco::Muon &muon)
Definition: MuonIsolationDQM.cc:659
MuonIsolationDQM::LOG_BINNING_ENABLED
int LOG_BINNING_ENABLED
Definition: MuonIsolationDQM.h:106
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:800
edm::InputTag
Definition: InputTag.h:15
MuonIsolationDQM::NUM_VARS_2D
static const int NUM_VARS_2D
Definition: MuonIsolationDQM.h:101
MuonIsolationDQM::nEvents
int nEvents
Definition: MuonIsolationDQM.h:145
MuonIsolationDQM::names_NVtxs
std::vector< std::string > names_NVtxs
Definition: MuonIsolationDQM.h:127
Generated for CMSSW Reference Manual by doxygen 1.8.16