CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
TopSingleLepton::MonitorEnsemble Class Reference

#include <TopSingleLeptonDQM.h>

Public Types

enum  Level { STANDARD, VERBOSE, DEBUG }
 different verbosity levels More...
 

Public Member Functions

void book (DQMStore::IBooker &ibooker)
 book histograms in subdirectory directory More...
 
void fill (const edm::Event &event, const edm::EventSetup &setup)
 fill monitor histograms with electronId and jetCorrections More...
 
 MonitorEnsemble (const char *label, const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
 default contructor More...
 
 ~MonitorEnsemble ()
 default destructor More...
 

Private Member Functions

bool booked (const std::string histName) const
 check if histogram was booked More...
 
void fill (const edm::Event &event, const edm::TriggerResults &triggerTable, std::string channel, const std::vector< std::string > labels) const
 fill trigger monitoring histograms More...
 
void fill (const std::string histName, double value) const
 fill histogram if it had been booked before More...
 
void fill (const std::string histName, double xValue, double yValue) const
 fill histogram if it had been booked before (2-dim version) More...
 
void fill (const std::string histName, double xValue, double yValue, double zValue) const
 fill histogram if it had been booked before (2-dim version) More...
 
std::string monitorPath (const std::string &label) const
 
std::string selectionPath (const std::string &label) const
 
void triggerBinLabels (std::string channel, const std::vector< std::string > labels)
 set configurable labels for trigger monitoring histograms More...
 

Private Attributes

edm::EDGetTokenT< reco::JetTagCollectionbtagCSV_
 
double btagCSVWP_
 
edm::EDGetTokenT< reco::JetTagCollectionbtagEff_
 btag discriminator labels More...
 
double btagEffWP_
 btag working points More...
 
edm::EDGetTokenT< reco::JetTagCollectionbtagPur_
 
double btagPurWP_
 
edm::EDGetTokenT< reco::JetTagCollectionbtagVtx_
 
double btagVtxWP_
 
std::string directory_
 
double eidCutValue_
 
edm::EDGetTokenT< edm::View< reco::PFCandidate > > elecs_
 
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > elecSelect_
 extra selection on electrons More...
 
edm::EDGetTokenT< edm::ValueMap< float > > electronId_
 electronId label More...
 
std::map< std::string, MonitorElement * > hists_
 histogram container More...
 
bool includeBTag_
 
edm::ESGetToken< JetCorrector, JetCorrectionsRecordjetCorrector_
 jetCorrector More...
 
edm::EDGetTokenT< reco::JetIDValueMapjetIDLabel_
 jetID as an extra selection type More...
 
std::unique_ptr< StringCutObjectSelector< reco::JetID > > jetIDSelect_
 
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetlooseSelection_
 
edm::EDGetTokenT< edm::View< reco::Jet > > jets_
 input sources for monitoring More...
 
std::string jetSelect_
 extra selection on jets More...
 
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetSelection_
 
std::string label_
 instance label More...
 
int logged_
 number of logged interesting events More...
 
double lowerEdge_
 mass window upper and lower edge More...
 
std::vector< edm::EDGetTokenT< edm::View< reco::MET > > > mets_
 considers a vector of METs More...
 
edm::EDGetTokenT< reco::JetCorrectormJetCorrector
 
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > muonIso_
 extra isolation criterion on muon More...
 
edm::EDGetTokenT< edm::View< reco::PFCandidate > > muons_
 
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > muonSelect_
 extra selection on muons More...
 
edm::EDGetTokenT< edm::View< reco::Vertex > > pvs_
 
std::unique_ptr< StringCutObjectSelector< reco::Vertex > > pvSelect_
 
edm::InputTag rhoTag
 
edm::EDConsumerBase tmpConsumerBase
 
std::vector< std::string > triggerPaths_
 
edm::EDGetTokenT< edm::TriggerResultstriggerTable_
 trigger table More...
 
double upperEdge_
 
Level verbosity_
 verbosity level for booking More...
 

Detailed Description

Definition at line 53 of file TopSingleLeptonDQM.h.

Member Enumeration Documentation

◆ Level

Constructor & Destructor Documentation

◆ MonitorEnsemble()

MonitorEnsemble::MonitorEnsemble ( const char *  label,
const edm::ParameterSet cfg,
edm::ConsumesCollector &&  iC 
)

default contructor

Definition at line 24 of file TopSingleLeptonDQM.cc.

References btagCSV_, btagCSVWP_, looper::cfg, DEBUG, directory_, eidCutValue_, singleTopDQM_cfi::elecExtras, elecs_, elecSelect_, electronId_, edm::ParameterSet::getParameter(), includeBTag_, jetCorrector_, singleTopDQM_cfi::jetExtras, jetIDLabel_, jetIDSelect_, jetlooseSelection_, jets_, jetSelect_, jetSelection_, lowerEdge_, singleTopDQM_cfi::massExtras, mets_, singleTopDQM_cfi::monitoring, singleTopDQM_cfi::muonExtras, muonIso_, muons_, muonSelect_, singleTopDQM_cfi::pvExtras, pvs_, pvSelect_, rhoTag, CalibrationSummaryClient_cfi::sources, STANDARD, AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, triggerPaths_, triggerTable_, upperEdge_, VERBOSE, and verbosity_.

25  : label_(label),
26  elecSelect_(nullptr),
27  pvSelect_(nullptr),
28  muonIso_(nullptr),
29  muonSelect_(nullptr),
30  jetIDSelect_(nullptr),
31  jetlooseSelection_(nullptr),
32  jetSelection_(nullptr),
33  includeBTag_(false),
34  lowerEdge_(-1.),
35  upperEdge_(-1.),
36  logged_(0) {
37  // sources have to be given; this PSet is not optional
38  edm::ParameterSet sources = cfg.getParameter<edm::ParameterSet>("sources");
41  pvs_ = iC.consumes<edm::View<reco::Vertex>>(sources.getParameter<edm::InputTag>("pvs"));
42  jets_ = iC.consumes<edm::View<reco::Jet>>(sources.getParameter<edm::InputTag>("jets"));
43  for (edm::InputTag const& tag : sources.getParameter<std::vector<edm::InputTag>>("mets"))
44  mets_.push_back(iC.consumes<edm::View<reco::MET>>(tag));
45  // electronExtras are optional; they may be omitted or
46  // empty
47  if (cfg.existsAs<edm::ParameterSet>("elecExtras")) {
48  // rho for PF isolation with EA corrections
49  // eventrhoToken_ =
50  // iC.consumes<double>(edm::InputTag("fixedGridRhoFastjetAll"));
51 
52  edm::ParameterSet elecExtras = cfg.getParameter<edm::ParameterSet>("elecExtras");
53  // select is optional; in case it's not found no
54  // selection will be applied
55  if (elecExtras.existsAs<std::string>("select")) {
56  elecSelect_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
57  elecExtras.getParameter<std::string>("select"));
58  }
59 
60  if (elecExtras.existsAs<std::string>("rho")) {
61  rhoTag = elecExtras.getParameter<edm::InputTag>("rho");
62  }
63  // electronId is optional; in case it's not found the
64  // InputTag will remain empty
65  if (elecExtras.existsAs<edm::ParameterSet>("electronId")) {
66  edm::ParameterSet elecId = elecExtras.getParameter<edm::ParameterSet>("electronId");
68  eidCutValue_ = elecId.getParameter<double>("cutValue");
69  }
70  }
71  // pvExtras are optional; they may be omitted or empty
72  if (cfg.existsAs<edm::ParameterSet>("pvExtras")) {
73  edm::ParameterSet pvExtras = cfg.getParameter<edm::ParameterSet>("pvExtras");
74  // select is optional; in case it's not found no
75  // selection will be applied
76  if (pvExtras.existsAs<std::string>("select")) {
77  pvSelect_ =
78  std::make_unique<StringCutObjectSelector<reco::Vertex>>(pvExtras.getParameter<std::string>("select"));
79  }
80  }
81  // muonExtras are optional; they may be omitted or empty
82  if (cfg.existsAs<edm::ParameterSet>("muonExtras")) {
83  edm::ParameterSet muonExtras = cfg.getParameter<edm::ParameterSet>("muonExtras");
84  // select is optional; in case it's not found no
85  // selection will be applied
86  if (muonExtras.existsAs<std::string>("select")) {
87  muonSelect_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
88  muonExtras.getParameter<std::string>("select"));
89  }
90  // isolation is optional; in case it's not found no
91  // isolation will be applied
92  if (muonExtras.existsAs<std::string>("isolation")) {
93  muonIso_ = std::make_unique<StringCutObjectSelector<reco::PFCandidate>>(
94  muonExtras.getParameter<std::string>("isolation"));
95  }
96  }
97 
98  // jetExtras are optional; they may be omitted or
99  // empty
100  if (cfg.existsAs<edm::ParameterSet>("jetExtras")) {
101  edm::ParameterSet jetExtras = cfg.getParameter<edm::ParameterSet>("jetExtras");
102  // jetCorrector is optional; in case it's not found
103  // the InputTag will remain empty
104  if (jetExtras.existsAs<std::string>("jetCorrector")) {
105  jetCorrector_ = iC.esConsumes(edm::ESInputTag("", jetExtras.getParameter<std::string>("jetCorrector")));
106  }
107  // read jetID information if it exists
108  if (jetExtras.existsAs<edm::ParameterSet>("jetID")) {
109  edm::ParameterSet jetID = jetExtras.getParameter<edm::ParameterSet>("jetID");
111  jetIDSelect_ =
112  std::make_unique<StringCutObjectSelector<reco::JetID>>(jetID.getParameter<std::string>("select"));
113  }
114  // select is optional; in case it's not found no
115  // selection will be applied (only implemented for
116  // CaloJets at the moment)
117  if (jetExtras.existsAs<std::string>("select")) {
118  jetSelect_ = jetExtras.getParameter<std::string>("select");
119  jetSelection_ = std::make_unique<StringCutObjectSelector<reco::PFJet>>(jetSelect_);
120  jetlooseSelection_ = std::make_unique<StringCutObjectSelector<reco::PFJet>>(
121  "chargedHadronEnergyFraction()>0 && chargedMultiplicity()>0 && chargedEmEnergyFraction()<0.99 && "
122  "neutralHadronEnergyFraction()<0.99 && neutralEmEnergyFraction()<0.99 && "
123  "(chargedMultiplicity()+neutralMultiplicity())>1");
124  }
125  // jetBDiscriminators are optional; in case they are
126  // not found the InputTag will remain empty; they
127  // consist of pairs of edm::JetFlavorAssociation's &
128  // corresponding working points
129  includeBTag_ = jetExtras.existsAs<edm::ParameterSet>("jetBTaggers");
130  if (includeBTag_) {
131  edm::ParameterSet btagCSV =
132  jetExtras.getParameter<edm::ParameterSet>("jetBTaggers").getParameter<edm::ParameterSet>("cvsVertex");
134  btagCSVWP_ = btagCSV.getParameter<double>("workingPoint");
135  }
136  }
137 
138  // triggerExtras are optional; they may be omitted or empty
139  if (cfg.existsAs<edm::ParameterSet>("triggerExtras")) {
140  edm::ParameterSet triggerExtras = cfg.getParameter<edm::ParameterSet>("triggerExtras");
142  triggerPaths_ = triggerExtras.getParameter<std::vector<std::string>>("paths");
143  }
144 
145  // massExtras is optional; in case it's not found no mass
146  // window cuts are applied for the same flavor monitor
147  // histograms
148  if (cfg.existsAs<edm::ParameterSet>("massExtras")) {
149  edm::ParameterSet massExtras = cfg.getParameter<edm::ParameterSet>("massExtras");
150  lowerEdge_ = massExtras.getParameter<double>("lowerEdge");
151  upperEdge_ = massExtras.getParameter<double>("upperEdge");
152  }
153 
154  // setup the verbosity level for booking histograms;
155  // per default the verbosity level will be set to
156  // STANDARD. This will also be the chosen level in
157  // the case when the monitoring PSet is not found
159  if (cfg.existsAs<edm::ParameterSet>("monitoring")) {
160  edm::ParameterSet monitoring = cfg.getParameter<edm::ParameterSet>("monitoring");
161  if (monitoring.getParameter<std::string>("verbosity") == "DEBUG")
162  verbosity_ = DEBUG;
163  if (monitoring.getParameter<std::string>("verbosity") == "VERBOSE")
165  if (monitoring.getParameter<std::string>("verbosity") == "STANDARD")
167  }
168  // and don't forget to do the histogram booking
169  directory_ = cfg.getParameter<std::string>("directory");
170  }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
int logged_
number of logged interesting events
edm::EDGetTokenT< reco::JetIDValueMap > jetIDLabel_
jetID as an extra selection type
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::EDGetTokenT< edm::View< reco::PFCandidate > > muons_
std::vector< std::string > triggerPaths_
edm::EDGetTokenT< edm::TriggerResults > triggerTable_
trigger table
std::unique_ptr< StringCutObjectSelector< reco::JetID > > jetIDSelect_
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > muonSelect_
extra selection on muons
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > elecSelect_
extra selection on electrons
edm::EDGetTokenT< reco::JetTagCollection > btagCSV_
std::vector< edm::EDGetTokenT< edm::View< reco::MET > > > mets_
considers a vector of METs
edm::EDGetTokenT< edm::View< reco::Jet > > jets_
input sources for monitoring
Level verbosity_
verbosity level for booking
double lowerEdge_
mass window upper and lower edge
char const * label
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetSelection_
std::string jetSelect_
extra selection on jets
edm::EDGetTokenT< edm::View< reco::Vertex > > pvs_
edm::EDGetTokenT< edm::View< reco::PFCandidate > > elecs_
edm::ESGetToken< JetCorrector, JetCorrectionsRecord > jetCorrector_
jetCorrector
std::string label_
instance label
std::unique_ptr< StringCutObjectSelector< reco::Vertex > > pvSelect_
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetlooseSelection_
edm::EDGetTokenT< edm::ValueMap< float > > electronId_
electronId label
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > muonIso_
extra isolation criterion on muon

◆ ~MonitorEnsemble()

TopSingleLepton::MonitorEnsemble::~MonitorEnsemble ( )
inline

default destructor

Definition at line 62 of file TopSingleLeptonDQM.h.

62 {};

Member Function Documentation

◆ book()

void MonitorEnsemble::book ( DQMStore::IBooker ibooker)

book histograms in subdirectory directory

Definition at line 172 of file TopSingleLeptonDQM.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), directory_, hists_, label_, dqm::implementation::NavigatorBase::setCurrentFolder(), STANDARD, AlCaHLTBitMon_QueryRunRegistry::string, triggerBinLabels(), triggerPaths_, VERBOSE, and verbosity_.

172  {
173  // set up the current directory path
174  std::string current(directory_);
175  current += label_;
176  ibooker.setCurrentFolder(current);
177 
178  // determine number of bins for trigger monitoring
179  // unsigned int nPaths = triggerPaths_.size();
180 
181  // --- [STANDARD] --- //
182  // Run Number
183  // hists_["RunNumb_"] = ibooker.book1D("RunNumber", "Run Nr.", 10, 0, 10);
184  // instantaneous luminosity
185  // hists_["InstLumi_"] = ibooker.book1D("InstLumi", "Inst. Lumi.", 100, 0., 1.e3);
186  // number of selected primary vertices
187  hists_["pvMult_"] = ibooker.book1D("PvMult", "N_{good pvs}", 50, 0., 50.);
188  // pt of the leading muon
189  hists_["muonPt_"] = ibooker.book1D("MuonPt", "pt(#mu TightId, TightIso)", 40, 0., 200.);
190  // muon multiplicity after tight Id
191  hists_["muonMult_"] = ibooker.book1D("MuonMult", "N_{loose}(#mu)", 10, 0., 10.);
192  // muon multiplicity after isolation
193  //hists_["muonMultIso_"] = ibooker.book1D("MuonMultIso",
194  // "N_{TightIso}(#mu)", 10, 0., 10.);
195  // muon multiplicity after tight Id and isolation
196  hists_["muonMultTight_"] = ibooker.book1D("MuonMultTight", "N_{TightIso,TightId}(#mu)", 10, 0., 10.);
197  // pt of the leading electron
198  hists_["elecPt_"] = ibooker.book1D("ElecPt", "pt(e TightId, TightIso)", 40, 0., 200.);
199  // electron multiplicity before std isolation
200  //hists_["elecMult_"] = ibooker.book1D("ElecMult", "N_{All}(e, tightId)", 10, 0., 10.);
201  // electron multiplicity after std isolation
202  //hists_["elecMultIso_"] = ibooker.book1D("ElecMultIso", "N_{Iso}(e, tightId, tightIso)", 10, 0., 10.);
203  // multiplicity of jets with pt>30
204  hists_["jetMult_"] = ibooker.book1D("JetMult", "N_{30}(jet)", 10, 0., 10.);
205  // multiplicity of loose jets with pt>30
206  hists_["jetMultLoose_"] = ibooker.book1D("JetMultLoose", "N_{30, loose}(jet)", 10, 0., 10.);
207 
208  // trigger efficiency estimates for single lepton triggers
209  // hists_["triggerEff_"] = ibooker.book1D("TriggerEff",
210  // "Eff(trigger)", nPaths, 0., nPaths);
211  // monitored trigger occupancy for single lepton triggers
212  // hists_["triggerMon_"] = ibooker.book1D("TriggerMon",
213  // "Mon(trigger)", nPaths, 0., nPaths);
214  // MET (pflow)
215  hists_["metPflow_"] = ibooker.book1D("METPflow", "MET_{Pflow}", 50, 0., 200.);
216  // W mass estimate
217  hists_["massW_"] = ibooker.book1D("MassW", "M(W)", 60, 0., 300.);
218  // Top mass estimate
219  hists_["massTop_"] = ibooker.book1D("MassTop", "M(Top)", 50, 0., 500.);
220  // b-tagged Top mass
221  hists_["massBTop_"] = ibooker.book1D("MassBTop", "M(Top, 1 b-tag)", 50, 0., 500.);
222  // set bin labels for trigger monitoring
224 
225  if (verbosity_ == STANDARD)
226  return;
227 
228  // --- [VERBOSE] --- //
229  // eta of the leading muon
230  hists_["muonEta_"] = ibooker.book1D("MuonEta", "#eta(#mu TightId, TightIso)", 30, -3., 3.);
231  // relative isolation of the candidate muon (depending on the decay channel)
232  hists_["muonRelIso_"] = ibooker.book1D("MuonRelIso", "Iso_{Rel}(#mu TightId) (#Delta#beta Corrected)", 50, 0., 1.);
233  // phi of the leading muon
234  hists_["muonPhi_"] = ibooker.book1D("MuonPhi", "#phi(#mu TightId, TightIso)", 40, -4., 4.);
235  // eta of the leading electron
236  hists_["elecEta_"] = ibooker.book1D("ElecEta", "#eta(e tightId)", 30, -3., 3.);
237  // std isolation variable of the leading electron
238  hists_["elecRelIso_"] = ibooker.book1D("ElecRelIso", "Iso_{Rel}(e TightId, TightIso)", 50, 0., 1.);
239  // phi of the leading electron
240  hists_["elecPhi_"] = ibooker.book1D("ElecPhi", "#phi(e tightId)", 40, -4., 4.);
241  // multiplicity of tight Id, tight Iso electorns
242  hists_["elecMultTight_"] = ibooker.book1D("ElecMultTight", "N_{TightIso,TightId}(e)", 10, 0., 10.);
243  // multiplicity of btagged jets (for track counting high efficiency) with
244  // pt(L2L3)>30
245  // hists_["jetMultBEff_"] = ibooker.book1D("JetMultBEff",
246  // "N_{30}(TCHE)", 10, 0., 10.);
247  // btag discriminator for track counting high efficiency for jets with
248  // pt(L2L3)>30
249  //hists_["jetBDiscEff_"] = ibooker.book1D("JetBDiscEff",
250  // "Disc_{TCHE}(jet)", 100, 0., 10.);
251  // eta of the 1. leading jet (corrected to L2+L3)
252  hists_["jet1Eta_"] = ibooker.book1D("Jet1Eta", "#eta_{30,loose}(jet1)", 60, -3., 3.);
253  // pt of the 1. leading jet (corrected to L2+L3)
254  hists_["jet1Pt_"] = ibooker.book1D("Jet1Pt", "pt_{30,loose}(jet1)", 60, 0., 300.);
255  // eta of the 2. leading jet (corrected to L2+L3)
256  hists_["jet2Eta_"] = ibooker.book1D("Jet2Eta", "#eta_{30,loose}(jet2)", 60, -3., 3.);
257  // pt of the 2. leading jet (corrected to L2+L3)
258  hists_["jet2Pt_"] = ibooker.book1D("Jet2Pt", "pt_{30,loose}(jet2)", 60, 0., 300.);
259  // eta of the 3. leading jet (corrected to L2+L3)
260  hists_["jet3Eta_"] = ibooker.book1D("Jet3Eta", "#eta_{30,loose}(jet3)", 60, -3., 3.);
261  // pt of the 3. leading jet (corrected to L2+L3)
262  hists_["jet3Pt_"] = ibooker.book1D("Jet3Pt", "pt_{30,loose}(jet3)", 60, 0., 300.);
263  // eta of the 4. leading jet (corrected to L2+L3)
264  hists_["jet4Eta_"] = ibooker.book1D("Jet4Eta", "#eta_{30,loose}(jet4)", 60, -3., 3.);
265  // pt of the 4. leading jet (corrected to L2+L3)
266  hists_["jet4Pt_"] = ibooker.book1D("Jet4Pt", "pt_{30,loose}(jet4)", 60, 0., 300.);
267  // dz for muons (to suppress cosmis)
268  hists_["muonDelZ_"] = ibooker.book1D("MuonDelZ", "d_{z}(#mu)", 50, -25., 25.);
269  // dxy for muons (to suppress cosmics)
270  hists_["muonDelXY_"] = ibooker.book2D("MuonDelXY", "d_{xy}(#mu)", 50, -0.1, 0.1, 50, -0.1, 0.1);
271 
272  // set axes titles for dxy for muons
273  hists_["muonDelXY_"]->setAxisTitle("x [cm]", 1);
274  hists_["muonDelXY_"]->setAxisTitle("y [cm]", 2);
275 
276  if (verbosity_ == VERBOSE)
277  return;
278 
279  // --- [DEBUG] --- //
280  // charged hadron isolation component of the candidate muon (depending on the
281  // decay channel)
282  hists_["muonChHadIso_"] = ibooker.book1D("MuonChHadIsoComp", "ChHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
283  // neutral hadron isolation component of the candidate muon (depending on the
284  // decay channel)
285  hists_["muonNeHadIso_"] = ibooker.book1D("MuonNeHadIsoComp", "NeHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
286  // photon isolation component of the candidate muon (depending on the decay
287  // channel)
288  hists_["muonPhIso_"] = ibooker.book1D("MuonPhIsoComp", "Photon_{IsoComponent}(#mu )", 50, 0., 5.);
289  // charged hadron isolation component of the candidate electron (depending on
290  // the decay channel)
291  hists_["elecChHadIso_"] = ibooker.book1D("ElectronChHadIsoComp", "ChHad_{IsoComponent}(e tightId)", 50, 0., 5.);
292  // neutral hadron isolation component of the candidate electron (depending on
293  // the decay channel)
294  hists_["elecNeHadIso_"] = ibooker.book1D("ElectronNeHadIsoComp", "NeHad_{IsoComponent}(e tightId)", 50, 0., 5.);
295  // photon isolation component of the candidate electron (depending on the
296  // decay channel)
297  hists_["elecPhIso_"] = ibooker.book1D("ElectronPhIsoComp", "Photon_{IsoComponent}(e tightId)", 50, 0., 5.);
298 
299  // multiplicity for combined secondary vertex
300  hists_["jetMultBCSVM_"] = ibooker.book1D("JetMultBCSVM", "N_{30}(CSVM)", 10, 0., 10.);
301  // btag discriminator for combined secondary vertex
302  hists_["jetBCSV_"] = ibooker.book1D("JetDiscCSV", "BJet Disc_{CSV}(JET)", 100, -1., 2.);
303  // pt of the 1. leading jet (uncorrected)
304  // hists_["jet1PtRaw_"] = ibooker.book1D("Jet1PtRaw", "pt_{Raw}(jet1)", 60, 0., 300.);
305  // pt of the 2. leading jet (uncorrected)
306  // hists_["jet2PtRaw_"] = ibooker.book1D("Jet2PtRaw", "pt_{Raw}(jet2)", 60, 0., 300.);
307  // pt of the 3. leading jet (uncorrected)
308  // hists_["jet3PtRaw_"] = ibooker.book1D("Jet3PtRaw", "pt_{Raw}(jet3)", 60, 0., 300.);
309  // pt of the 4. leading jet (uncorrected)
310  // hists_["jet4PtRaw_"] = ibooker.book1D("Jet4PtRaw", "pt_{Raw}(jet4)", 60, 0., 300.);
311  // selected events
312  hists_["eventLogger_"] = ibooker.book2D("EventLogger", "Logged Events", 9, 0., 9., 10, 0., 10.);
313 
314  // set axes titles for selected events
315  hists_["eventLogger_"]->getTH1()->SetOption("TEXT");
316  hists_["eventLogger_"]->setBinLabel(1, "Run", 1);
317  hists_["eventLogger_"]->setBinLabel(2, "Block", 1);
318  hists_["eventLogger_"]->setBinLabel(3, "Event", 1);
319  hists_["eventLogger_"]->setBinLabel(4, "pt_{30,loose}(jet1)", 1);
320  hists_["eventLogger_"]->setBinLabel(5, "pt_{30,loose}(jet2)", 1);
321  hists_["eventLogger_"]->setBinLabel(6, "pt_{30,loose}(jet3)", 1);
322  hists_["eventLogger_"]->setBinLabel(7, "pt_{30,loose}(jet4)", 1);
323  hists_["eventLogger_"]->setBinLabel(8, "M_{W}", 1);
324  hists_["eventLogger_"]->setBinLabel(9, "M_{Top}", 1);
325  hists_["eventLogger_"]->setAxisTitle("logged evts", 2);
326  return;
327  }
std::map< std::string, MonitorElement * > hists_
histogram container
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
std::vector< std::string > triggerPaths_
Level verbosity_
verbosity level for booking
std::string label_
instance label
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
void triggerBinLabels(std::string channel, const std::vector< std::string > labels)
set configurable labels for trigger monitoring histograms
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98

◆ booked()

bool TopSingleLepton::MonitorEnsemble::booked ( const std::string  histName) const
inlineprivate

check if histogram was booked

Definition at line 86 of file TopSingleLeptonDQM.h.

References hists_.

Referenced by fill().

86 { return hists_.find(histName) != hists_.end(); };
std::map< std::string, MonitorElement * > hists_
histogram container

◆ fill() [1/5]

void MonitorEnsemble::fill ( const edm::Event event,
const edm::EventSetup setup 
)

fill monitor histograms with electronId and jetCorrections

Definition at line 329 of file TopSingleLeptonDQM.cc.

References btagCSV_, btagCSVWP_, pfClustersFromHGC3DClusters_cfi::corrector, eidCutValue_, singleTopDQM_cfi::elecs, elecs_, elecSelect_, topSingleLeptonDQM_PU_cfi::electronId, electronId_, reco::LeafCandidate::eta(), dqmMemoryStats::float, edm::ESGetToken< ESProduct, ESRecord >::hasValidIndex(), hists_, heavyIonCSV_trainingSettings::idx, includeBTag_, edm::ESGetToken< ESProduct, ESRecord >::isInitialized(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, jetCorrector_, jetlooseSelection_, singleTopDQM_cfi::jets, jets_, jetSelection_, logged_, lowerEdge_, Calculate::massBTopQuark(), Calculate::massTopQuark(), Calculate::massWBoson(), SiStripPI::max, TopSingleLepton::MAXJETS, BTaggingMonitor_cfi::met, mets_, VarParsing::mult, PDWG_BPHSkim_cff::muons, muons_, muonSelect_, reco::Vertex::position(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), MetAnalyzer::pv(), FSQDQM_cfi::pvs, pvs_, pvSelect_, rhoTag, reco::Jet::scaleEnergy(), EgammaValidation_Wenu_cff::sel, singleTopDQM_cfi::setup, triggerPaths_, triggerTable_, funct::true, upperEdge_, trackerHitRTTI::vector, pseudoTop_cfi::wMass, and TopSingleLepton::WMASS.

Referenced by fill().

329  {
330  // fetch trigger event if configured such
332 
334  if (!event.getByToken(triggerTable_, triggerTable))
335  return;
336  }
337 
338  /*
339  ------------------------------------------------------------
340 
341  Primary Vertex Monitoring
342 
343  ------------------------------------------------------------
344  */
345  // fill monitoring plots for primary verices
347 
348  if (!event.getByToken(pvs_, pvs))
349  return;
350  const reco::Vertex& Pvertex = pvs->front();
351  unsigned int pvMult = 0;
352  for (edm::View<reco::Vertex>::const_iterator pv = pvs->begin(); pv != pvs->end(); ++pv) {
353  if (!pvSelect_ || (*pvSelect_)(*pv))
354  pvMult++;
355  }
356  fill("pvMult_", pvMult);
357 
358  /*
359  ------------------------------------------------------------
360 
361  Electron Monitoring
362 
363  ------------------------------------------------------------
364  */
365 
366  // fill monitoring plots for electrons
368  edm::Handle<double> _rhoHandle;
369  event.getByLabel(rhoTag, _rhoHandle);
370  if (!event.getByToken(elecs_, elecs))
371  return;
372 
373  // check availability of electron id
375  if (!electronId_.isUninitialized()) {
376  if (!event.getByToken(electronId_, electronId))
377  return;
378  }
379 
380  // loop electron collection
381  unsigned int eMult = 0, eMultIso = 0;
382  std::vector<const reco::PFCandidate*> isoElecs;
383  for (edm::View<reco::PFCandidate>::const_iterator elec = elecs->begin(); elec != elecs->end(); ++elec) {
384  if (elec->gsfElectronRef().isNull()) {
385  continue;
386  }
387  reco::GsfElectronRef gsf_el = elec->gsfElectronRef();
388  // restrict to electrons with good electronId
390  ? true
391  : ((double)(*electronId)[gsf_el] >=
392  eidCutValue_)) { //This Electron Id is not currently used, but we can keep this for future needs
393  if (!elecSelect_ || (*elecSelect_)(*elec)) {
394  double el_ChHadIso = gsf_el->pfIsolationVariables().sumChargedHadronPt;
395  double el_NeHadIso = gsf_el->pfIsolationVariables().sumNeutralHadronEt;
396  double el_PhIso = gsf_el->pfIsolationVariables().sumPhotonEt;
397  double absEta = std::fabs(gsf_el->superCluster()->eta());
398 
399  //Effective Area computation
400  double eA = 0;
401  if (absEta < 1.000)
402  eA = 0.1703;
403  else if (absEta < 1.479)
404  eA = 0.1715;
405  else if (absEta < 2.000)
406  eA = 0.1213;
407  else if (absEta < 2.200)
408  eA = 0.1230;
409  else if (absEta < 2.300)
410  eA = 0.1635;
411  else if (absEta < 2.400)
412  eA = 0.1937;
413  else if (absEta < 5.000)
414  eA = 0.2393;
415 
416  double rho = _rhoHandle.isValid() ? (float)(*_rhoHandle) : 0;
417  double el_pfRelIso = (el_ChHadIso + max(0., el_NeHadIso + el_PhIso - rho * eA)) / gsf_el->pt();
418 
419  //Only TightId
420  if (eMult == 0) { // Restricted to the leading tight electron
421  fill("elecRelIso_", el_pfRelIso);
422  fill("elecChHadIso_", el_ChHadIso);
423  fill("elecNeHadIso_", el_NeHadIso);
424  fill("elecPhIso_", el_PhIso);
425  }
426  ++eMult;
427 
428  if (!((el_pfRelIso < 0.0588 && absEta < 1.479) || (el_pfRelIso < 0.0571 && absEta > 1.479)))
429  continue; // PF Isolation with Effective Area Corrections according to https://twiki.cern.ch/twiki/bin/viewauth/CMS/CutBasedElectronIdentificationRun2
430 
431  // TightId and TightIso
432  if (eMultIso == 0) { //Only leading
433  fill("elecPt_", elec->pt());
434  fill("elecEta_", elec->eta());
435  fill("elecPhi_", elec->phi());
436  }
437  ++eMultIso;
438  }
439  }
440  }
441  //fill("elecMult_", eMult);
442  fill("elecMultTight_", eMultIso);
443 
444  /*
445  ------------------------------------------------------------
446 
447  Muon Monitoring
448 
449  ------------------------------------------------------------
450  */
451 
452  // fill monitoring plots for muons
453  unsigned int mMult = 0, mTight = 0, mTightId = 0;
454 
457 
458  if (!event.getByToken(muons_, muons))
459  return;
460 
461  for (edm::View<reco::PFCandidate>::const_iterator muonit = muons->begin(); muonit != muons->end(); ++muonit) {
462  if (muonit->muonRef().isNull())
463  continue;
464  reco::MuonRef muon = muonit->muonRef();
465 
466  // restrict to globalMuons
467  if (muon->isGlobalMuon()) {
468  fill("muonDelZ_", muon->innerTrack()->vz()); // CB using inner track!
469  fill("muonDelXY_", muon->innerTrack()->vx(), muon->innerTrack()->vy());
470 
471  // apply preselection
472  if ((!muonSelect_ || (*muonSelect_)(*muonit))) {
473  mMult++;
474  double chHadPt = muon->pfIsolationR04().sumChargedHadronPt;
475  double neHadEt = muon->pfIsolationR04().sumNeutralHadronEt;
476  double phoEt = muon->pfIsolationR04().sumPhotonEt;
477  double pfRelIso = (chHadPt + max(0., neHadEt + phoEt - 0.5 * muon->pfIsolationR04().sumPUPt)) /
478  muon->pt(); // CB dBeta corrected iso!
479 
480  if (!(muon->isGlobalMuon() && muon->isPFMuon() && muon->globalTrack()->normalizedChi2() < 10. &&
481  muon->globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && muon->numberOfMatchedStations() > 1 &&
482  muon->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 &&
483  muon->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
484  fabs(muon->muonBestTrack()->dxy(Pvertex.position())) < 0.2 &&
485  fabs(muon->muonBestTrack()->dz(Pvertex.position())) < 0.5))
486  continue; //Only tight muons
487 
488  if (mTightId == 0) {
489  fill("muonRelIso_", pfRelIso);
490  fill("muonChHadIso_", chHadPt);
491  fill("muonNeHadIso_", neHadEt);
492  fill("muonPhIso_", phoEt);
493  }
494  mTightId++;
495 
496  if (!(pfRelIso < 0.15))
497  continue; //Tight Iso
498 
499  if (mTight == 0) { //Leading tightId tightIso muon
500  fill("muonPt_", muon->pt());
501  fill("muonEta_", muon->eta());
502  fill("muonPhi_", muon->phi());
503  }
504  mTight++;
505  }
506  }
507  }
508  fill("muonMult_", mMult);
509  fill("muonMultTight_", mTight);
510 
511  /*
512  ------------------------------------------------------------
513 
514  Jet Monitoring
515 
516  ------------------------------------------------------------
517  */
518 
519  const JetCorrector* corrector = nullptr;
521  // check whether a jet correcto is in the event setup or not
522  if (setup.find(edm::eventsetup::EventSetupRecordKey::makeKey<JetCorrectionsRecord>())) {
523  corrector = &setup.getData(jetCorrector_);
524  } else {
525  edm::LogVerbatim("TopDiLeptonOfflineDQM") << "\n"
526  << "-----------------------------------------------------------------"
527  "-------------------- \n"
528  << " No JetCorrectionsRecord available from EventSetup: "
529  " \n"
530  << " - Jets will not be corrected. "
531  " \n"
532  << " - If you want to change this add the following lines to your "
533  "cfg file: \n"
534  << " "
535  " \n"
536  << " ## load jet corrections "
537  " \n"
538  << " "
539  "process.load(\"JetMETCorrections.Configuration."
540  "JetCorrectionServicesAllAlgos_cff\") \n"
541  << " process.prefer(\"ak5CaloL2L3\") "
542  " \n"
543  << " "
544  " \n"
545  << "-----------------------------------------------------------------"
546  "-------------------- \n";
547  }
548  }
549 
550  // check availability of the btaggers
551  edm::Handle<reco::JetTagCollection> btagEff, btagPur, btagVtx, btagCSV;
552  if (includeBTag_) {
553  if (!event.getByToken(btagCSV_, btagCSV))
554  return;
555  }
556 
557  // loop jet collection
558  std::vector<reco::Jet> correctedJets;
559  std::vector<double> JetTagValues;
560  unsigned int mult = 0, multLoose = 0, multCSV = 0;
561 
563  if (!event.getByToken(jets_, jets)) {
564  return;
565  }
566 
567  for (edm::View<reco::Jet>::const_iterator jet = jets->begin(); jet != jets->end(); ++jet) {
568  bool isLoose = false;
569  // check additional jet selection for pf jets
570  if (dynamic_cast<const reco::PFJet*>(&*jet)) {
571  reco::PFJet sel = dynamic_cast<const reco::PFJet&>(*jet);
572  if ((*jetlooseSelection_)(sel))
573  isLoose = true;
574  sel.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
575  if (!(*jetSelection_)(sel))
576  continue;
577  }
578 
579  // prepare jet to fill monitor histograms
580  reco::Jet monitorJet = *jet;
581 
582  ++mult; // determine jet (no Id) multiplicity
583  monitorJet.scaleEnergy(corrector ? corrector->correction(*jet) : 1.);
584 
585  if (isLoose) { //Loose Id
586  unsigned int idx = jet - jets->begin();
587  correctedJets.push_back(monitorJet);
588  if (includeBTag_) {
589  // fill b-discriminators
590  edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);
591  fill("jetBCSV_", (*btagCSV)[jetRef]);
592  if ((*btagCSV)[jetRef] > btagCSVWP_)
593  ++multCSV;
594  // Fill a vector with Jet b-tag WP for later M3+1tag calculation: CSV
595  // tagger
596  JetTagValues.push_back((*btagCSV)[jetRef]);
597  }
598 
599  // fill pt/eta for the leading four jets
600  if (multLoose == 0) {
601  fill("jet1Pt_", monitorJet.pt());
602  fill("jet1Eta_", monitorJet.eta());
603  };
604  if (multLoose == 1) {
605  fill("jet2Pt_", monitorJet.pt());
606  fill("jet2Eta_", monitorJet.eta());
607  }
608  if (multLoose == 2) {
609  fill("jet3Pt_", monitorJet.pt());
610  fill("jet3Eta_", monitorJet.eta());
611  }
612  if (multLoose == 3) {
613  fill("jet4Pt_", monitorJet.pt());
614  fill("jet4Eta_", monitorJet.eta());
615  }
616  multLoose++;
617  }
618  }
619  fill("jetMult_", mult);
620  fill("jetMultLoose_", multLoose);
621  fill("jetMultBCSVM_", multCSV);
622 
623  /*
624  ------------------------------------------------------------
625 
626  MET Monitoring
627 
628  ------------------------------------------------------------
629  */
630 
631  // fill monitoring histograms for met
632  for (std::vector<edm::EDGetTokenT<edm::View<reco::MET>>>::const_iterator met_ = mets_.begin(); met_ != mets_.end();
633  ++met_) {
635  if (!event.getByToken(*met_, met))
636  continue;
637  if (met->begin() != met->end()) { //If we ever have to use more than one type of met again
638  unsigned int idx = met_ - mets_.begin();
639  if (idx == 0)
640  fill("metPflow_", met->begin()->et());
641  }
642  }
643 
644  /*
645  ------------------------------------------------------------
646 
647  Event Monitoring
648 
649  ------------------------------------------------------------
650  */
651 
652  // fill W boson and top mass estimates
653 
654  Calculate eventKinematics(MAXJETS, WMASS);
655  double wMass = eventKinematics.massWBoson(correctedJets);
656  double topMass = eventKinematics.massTopQuark(correctedJets);
657  if (wMass >= 0 && topMass >= 0) {
658  fill("massW_", wMass);
659  fill("massTop_", topMass);
660  }
661 
662  // Fill M3 with Btag (CSV Tight) requirement
663 
664  if (!includeBTag_)
665  return;
666  if (correctedJets.size() != JetTagValues.size())
667  return;
668  double btopMass = eventKinematics.massBTopQuark(correctedJets, JetTagValues, btagCSVWP_);
669  if (btopMass >= 0)
670  fill("massBTop_", btopMass);
671 
672  // fill plots for trigger monitoring
673  if ((lowerEdge_ == -1. && upperEdge_ == -1.) || (lowerEdge_ < wMass && wMass < upperEdge_)) {
675  fill(event, *triggerTable, "trigger", triggerPaths_);
676  if (logged_ <= hists_.find("eventLogger_")->second->getNbinsY()) {
677  // log runnumber, lumi block, event number & some
678  // more pysics infomation for interesting events
679  fill("eventLogger_", 0.5, logged_ + 0.5, event.eventAuxiliary().run());
680  fill("eventLogger_", 1.5, logged_ + 0.5, event.eventAuxiliary().luminosityBlock());
681  fill("eventLogger_", 2.5, logged_ + 0.5, event.eventAuxiliary().event());
682  if (!correctedJets.empty())
683  fill("eventLogger_", 3.5, logged_ + 0.5, correctedJets[0].pt());
684  if (correctedJets.size() > 1)
685  fill("eventLogger_", 4.5, logged_ + 0.5, correctedJets[1].pt());
686  if (correctedJets.size() > 2)
687  fill("eventLogger_", 5.5, logged_ + 0.5, correctedJets[2].pt());
688  if (correctedJets.size() > 3)
689  fill("eventLogger_", 6.5, logged_ + 0.5, correctedJets[3].pt());
690  fill("eventLogger_", 7.5, logged_ + 0.5, wMass);
691  fill("eventLogger_", 8.5, logged_ + 0.5, topMass);
692  ++logged_;
693  }
694  }
695  }
std::map< std::string, MonitorElement * > hists_
histogram container
Log< level::Info, true > LogVerbatim
int logged_
number of logged interesting events
edm::EDGetTokenT< edm::View< reco::PFCandidate > > muons_
virtual void scaleEnergy(double fScale)
scale energy of the jet
double pt() const final
transverse momentum
std::vector< std::string > triggerPaths_
const Point & position() const
position
Definition: Vertex.h:127
Base class for all types of Jets.
Definition: Jet.h:20
edm::EDGetTokenT< edm::TriggerResults > triggerTable_
trigger table
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > muonSelect_
extra selection on muons
constexpr bool isUninitialized() const noexcept
Definition: EDGetToken.h:99
std::unique_ptr< StringCutObjectSelector< reco::PFCandidate > > elecSelect_
extra selection on electrons
edm::EDGetTokenT< reco::JetTagCollection > btagCSV_
Jets made from PFObjects.
Definition: PFJet.h:20
std::vector< edm::EDGetTokenT< edm::View< reco::MET > > > mets_
considers a vector of METs
electronId
when omitted electron plots will be filled w/o cut on electronId
edm::EDGetTokenT< edm::View< reco::Jet > > jets_
input sources for monitoring
static const double WMASS
constexpr bool hasValidIndex() const noexcept
Definition: ESGetToken.h:74
double lowerEdge_
mass window upper and lower edge
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetSelection_
Helper class for the calculation of a top and a W boson mass estime.
constexpr bool isInitialized() const noexcept
Definition: ESGetToken.h:72
def pv(vc)
Definition: MetAnalyzer.py:7
edm::EDGetTokenT< edm::View< reco::Vertex > > pvs_
edm::EDGetTokenT< edm::View< reco::PFCandidate > > elecs_
edm::ESGetToken< JetCorrector, JetCorrectionsRecord > jetCorrector_
jetCorrector
bool isValid() const
Definition: HandleBase.h:70
std::unique_ptr< StringCutObjectSelector< reco::Vertex > > pvSelect_
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
void fill(const edm::Event &event, const edm::EventSetup &setup)
fill monitor histograms with electronId and jetCorrections
std::unique_ptr< StringCutObjectSelector< reco::PFJet > > jetlooseSelection_
static const unsigned int MAXJETS
edm::EDGetTokenT< edm::ValueMap< float > > electronId_
electronId label
Definition: event.py:1
double eta() const final
momentum pseudorapidity

◆ fill() [2/5]

void MonitorEnsemble::fill ( const edm::Event event,
const edm::TriggerResults triggerTable,
std::string  channel,
const std::vector< std::string >  labels 
) const
inlineprivate

fill trigger monitoring histograms

Definition at line 196 of file TopSingleLeptonDQM.h.

References accept(), fill(), hists_, heavyIonCSV_trainingSettings::idx, SummaryClient_cfi::labels, monitorPath(), and selectionPath().

199  {
200  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
201  if (accept(event, triggerTable, monitorPath(labels[idx]))) {
202  fill(channel + "Mon_", idx + 0.5);
203  // take care to fill triggerMon_ before evts is being called
204  int evts = hists_.find(channel + "Mon_")->second->getBinContent(idx + 1);
205  double value = hists_.find(channel + "Eff_")->second->getBinContent(idx + 1);
206  fill(
207  channel + "Eff_", idx + 0.5, 1. / evts * (accept(event, triggerTable, selectionPath(labels[idx])) - value));
208  }
209  }
210  }
std::map< std::string, MonitorElement * > hists_
histogram container
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
std::string monitorPath(const std::string &label) const
std::string selectionPath(const std::string &label) const
Definition: value.py:1
void fill(const edm::Event &event, const edm::EventSetup &setup)
fill monitor histograms with electronId and jetCorrections
Definition: event.py:1

◆ fill() [3/5]

void TopSingleLepton::MonitorEnsemble::fill ( const std::string  histName,
double  value 
) const
inlineprivate

fill histogram if it had been booked before

Definition at line 88 of file TopSingleLeptonDQM.h.

References booked(), and hists_.

88  {
89  if (booked(histName))
90  hists_.find(histName)->second->Fill(value);
91  };
std::map< std::string, MonitorElement * > hists_
histogram container
Definition: value.py:1
bool booked(const std::string histName) const
check if histogram was booked

◆ fill() [4/5]

void TopSingleLepton::MonitorEnsemble::fill ( const std::string  histName,
double  xValue,
double  yValue 
) const
inlineprivate

fill histogram if it had been booked before (2-dim version)

Definition at line 93 of file TopSingleLeptonDQM.h.

References booked(), and hists_.

93  {
94  if (booked(histName))
95  hists_.find(histName)->second->Fill(xValue, yValue);
96  };
std::map< std::string, MonitorElement * > hists_
histogram container
bool booked(const std::string histName) const
check if histogram was booked

◆ fill() [5/5]

void TopSingleLepton::MonitorEnsemble::fill ( const std::string  histName,
double  xValue,
double  yValue,
double  zValue 
) const
inlineprivate

fill histogram if it had been booked before (2-dim version)

Definition at line 98 of file TopSingleLeptonDQM.h.

References booked(), and hists_.

98  {
99  if (booked(histName))
100  hists_.find(histName)->second->Fill(xValue, yValue, zValue);
101  };
std::map< std::string, MonitorElement * > hists_
histogram container
bool booked(const std::string histName) const
check if histogram was booked

◆ monitorPath()

std::string TopSingleLepton::MonitorEnsemble::monitorPath ( const std::string &  label) const
inlineprivate

deduce monitorPath from label, the label is expected to be of type 'selectionPath:monitorPath'

Definition at line 72 of file TopSingleLeptonDQM.h.

References label.

Referenced by fill(), and triggerBinLabels().

72 { return label.substr(label.find(':') + 1); };
char const * label

◆ selectionPath()

std::string TopSingleLepton::MonitorEnsemble::selectionPath ( const std::string &  label) const
inlineprivate

deduce selectionPath from label, the label is expected to be of type 'selectionPath:monitorPath'

Definition at line 75 of file TopSingleLeptonDQM.h.

References label.

Referenced by fill(), and triggerBinLabels().

75 { return label.substr(0, label.find(':')); };
char const * label

◆ triggerBinLabels()

void MonitorEnsemble::triggerBinLabels ( std::string  channel,
const std::vector< std::string >  labels 
)
inlineprivate

set configurable labels for trigger monitoring histograms

Definition at line 188 of file TopSingleLeptonDQM.h.

References hists_, heavyIonCSV_trainingSettings::idx, SummaryClient_cfi::labels, monitorPath(), and selectionPath().

Referenced by book().

188  {
189  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
190  hists_[channel + "Mon_"]->setBinLabel(idx + 1, "[" + monitorPath(labels[idx]) + "]", 1);
191  hists_[channel + "Eff_"]->setBinLabel(
192  idx + 1, "[" + selectionPath(labels[idx]) + "]|[" + monitorPath(labels[idx]) + "]", 1);
193  }
194  }
std::map< std::string, MonitorElement * > hists_
histogram container
std::string monitorPath(const std::string &label) const
std::string selectionPath(const std::string &label) const

Member Data Documentation

◆ btagCSV_

edm::EDGetTokenT<reco::JetTagCollection> TopSingleLepton::MonitorEnsemble::btagCSV_
private

Definition at line 172 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ btagCSVWP_

double TopSingleLepton::MonitorEnsemble::btagCSVWP_
private

Definition at line 174 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ btagEff_

edm::EDGetTokenT<reco::JetTagCollection> TopSingleLepton::MonitorEnsemble::btagEff_
private

btag discriminator labels

Definition at line 172 of file TopSingleLeptonDQM.h.

◆ btagEffWP_

double TopSingleLepton::MonitorEnsemble::btagEffWP_
private

btag working points

Definition at line 174 of file TopSingleLeptonDQM.h.

◆ btagPur_

edm::EDGetTokenT<reco::JetTagCollection> TopSingleLepton::MonitorEnsemble::btagPur_
private

Definition at line 172 of file TopSingleLeptonDQM.h.

◆ btagPurWP_

double TopSingleLepton::MonitorEnsemble::btagPurWP_
private

Definition at line 174 of file TopSingleLeptonDQM.h.

◆ btagVtx_

edm::EDGetTokenT<reco::JetTagCollection> TopSingleLepton::MonitorEnsemble::btagVtx_
private

Definition at line 172 of file TopSingleLeptonDQM.h.

◆ btagVtxWP_

double TopSingleLepton::MonitorEnsemble::btagVtxWP_
private

Definition at line 174 of file TopSingleLeptonDQM.h.

◆ directory_

std::string TopSingleLepton::MonitorEnsemble::directory_
private

Definition at line 185 of file TopSingleLeptonDQM.h.

Referenced by book(), and MonitorEnsemble().

◆ eidCutValue_

double TopSingleLepton::MonitorEnsemble::eidCutValue_
private

electronId pattern we expect the following pattern: 0: fails 1: passes electron ID only 2: passes electron Isolation only 3: passes electron ID and Isolation only 4: passes conversion rejection 5: passes conversion rejection and ID 6: passes conversion rejection and Isolation 7: passes the whole selection As described on https://twiki.cern.ch/twiki/bin/view/CMS/SimpleCutBasedEleID

Definition at line 138 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ elecs_

edm::EDGetTokenT<edm::View<reco::PFCandidate> > TopSingleLepton::MonitorEnsemble::elecs_
private

Definition at line 113 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ elecSelect_

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > TopSingleLepton::MonitorEnsemble::elecSelect_
private

extra selection on electrons

Definition at line 145 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ electronId_

edm::EDGetTokenT<edm::ValueMap<float> > TopSingleLepton::MonitorEnsemble::electronId_
private

electronId label

Definition at line 122 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ hists_

std::map<std::string, MonitorElement*> TopSingleLepton::MonitorEnsemble::hists_
private

histogram container

Definition at line 182 of file TopSingleLeptonDQM.h.

Referenced by book(), booked(), fill(), and triggerBinLabels().

◆ includeBTag_

bool TopSingleLepton::MonitorEnsemble::includeBTag_
private

include btag information or not to be determined from the cfg

Definition at line 170 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ jetCorrector_

edm::ESGetToken<JetCorrector, JetCorrectionsRecord> TopSingleLepton::MonitorEnsemble::jetCorrector_
private

jetCorrector

Definition at line 158 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ jetIDLabel_

edm::EDGetTokenT<reco::JetIDValueMap> TopSingleLepton::MonitorEnsemble::jetIDLabel_
private

jetID as an extra selection type

Definition at line 161 of file TopSingleLeptonDQM.h.

Referenced by MonitorEnsemble().

◆ jetIDSelect_

std::unique_ptr<StringCutObjectSelector<reco::JetID> > TopSingleLepton::MonitorEnsemble::jetIDSelect_
private

Definition at line 163 of file TopSingleLeptonDQM.h.

Referenced by MonitorEnsemble().

◆ jetlooseSelection_

std::unique_ptr<StringCutObjectSelector<reco::PFJet> > TopSingleLepton::MonitorEnsemble::jetlooseSelection_
private

Definition at line 166 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ jets_

edm::EDGetTokenT<edm::View<reco::Jet> > TopSingleLepton::MonitorEnsemble::jets_
private

input sources for monitoring

Definition at line 111 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ jetSelect_

std::string TopSingleLepton::MonitorEnsemble::jetSelect_
private

extra selection on jets

Definition at line 165 of file TopSingleLeptonDQM.h.

Referenced by MonitorEnsemble().

◆ jetSelection_

std::unique_ptr<StringCutObjectSelector<reco::PFJet> > TopSingleLepton::MonitorEnsemble::jetSelection_
private

Definition at line 167 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ label_

std::string TopSingleLepton::MonitorEnsemble::label_
private

◆ logged_

int TopSingleLepton::MonitorEnsemble::logged_
private

number of logged interesting events

Definition at line 179 of file TopSingleLeptonDQM.h.

Referenced by fill().

◆ lowerEdge_

double TopSingleLepton::MonitorEnsemble::lowerEdge_
private

mass window upper and lower edge

Definition at line 176 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ mets_

std::vector<edm::EDGetTokenT<edm::View<reco::MET> > > TopSingleLepton::MonitorEnsemble::mets_
private

considers a vector of METs

Definition at line 109 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ mJetCorrector

edm::EDGetTokenT<reco::JetCorrector> TopSingleLepton::MonitorEnsemble::mJetCorrector
private

Definition at line 123 of file TopSingleLeptonDQM.h.

◆ muonIso_

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > TopSingleLepton::MonitorEnsemble::muonIso_
private

extra isolation criterion on muon

Definition at line 152 of file TopSingleLeptonDQM.h.

Referenced by MonitorEnsemble().

◆ muons_

edm::EDGetTokenT<edm::View<reco::PFCandidate> > TopSingleLepton::MonitorEnsemble::muons_
private

Definition at line 112 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ muonSelect_

std::unique_ptr<StringCutObjectSelector<reco::PFCandidate> > TopSingleLepton::MonitorEnsemble::muonSelect_
private

extra selection on muons

Definition at line 155 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ pvs_

edm::EDGetTokenT<edm::View<reco::Vertex> > TopSingleLepton::MonitorEnsemble::pvs_
private

Definition at line 114 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ pvSelect_

std::unique_ptr<StringCutObjectSelector<reco::Vertex> > TopSingleLepton::MonitorEnsemble::pvSelect_
private

extra selection on primary vertices; meant to investigate the pile-up effect

Definition at line 149 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ rhoTag

edm::InputTag TopSingleLepton::MonitorEnsemble::rhoTag
private

Definition at line 141 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ tmpConsumerBase

edm::EDConsumerBase TopSingleLepton::MonitorEnsemble::tmpConsumerBase
private

Definition at line 183 of file TopSingleLeptonDQM.h.

◆ triggerPaths_

std::vector<std::string> TopSingleLepton::MonitorEnsemble::triggerPaths_
private

trigger paths for monitoring, expected to be of form signalPath:MonitorPath

Definition at line 119 of file TopSingleLeptonDQM.h.

Referenced by book(), fill(), and MonitorEnsemble().

◆ triggerTable_

edm::EDGetTokenT<edm::TriggerResults> TopSingleLepton::MonitorEnsemble::triggerTable_
private

trigger table

Definition at line 116 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ upperEdge_

double TopSingleLepton::MonitorEnsemble::upperEdge_
private

Definition at line 176 of file TopSingleLeptonDQM.h.

Referenced by fill(), and MonitorEnsemble().

◆ verbosity_

Level TopSingleLepton::MonitorEnsemble::verbosity_
private

verbosity level for booking

Definition at line 101 of file TopSingleLeptonDQM.h.

Referenced by book(), and MonitorEnsemble().