CMS 3D CMS Logo

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

#include <SingleTopTChannelLeptonDQM_miniAOD.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_
 
std::unique_ptr< StringCutObjectSelector< pat::Electron > > elecIso_
 extra isolation criterion on electron More...
 
edm::EDGetTokenT< edm::View< pat::Electron > > elecs_
 
std::unique_ptr< StringCutObjectSelector< pat::Electron > > 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_
 
std::string jetCorrector_
 jetCorrector More...
 
edm::EDGetTokenT< reco::JetIDValueMapjetIDLabel_
 jetID as an extra selection type More...
 
std::unique_ptr< StringCutObjectSelector< reco::JetID > > jetIDSelect_
 extra jetID selection on calo jets More...
 
edm::EDGetTokenT< edm::View< pat::Jet > > jets_
 input sources for monitoring More...
 
std::unique_ptr< StringCutObjectSelector< pat::Jet > > jetSelect
 
std::string jetSelect_
 
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< pat::MET > > > mets_
 considers a vector of METs More...
 
std::unique_ptr< StringCutObjectSelector< pat::Muon > > muonIso_
 extra isolation criterion on muon More...
 
edm::EDGetTokenT< edm::View< pat::Muon > > muons_
 
std::unique_ptr< StringCutObjectSelector< pat::Muon > > 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 30 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Member Enumeration Documentation

◆ Level

different verbosity levels

Enumerator
STANDARD 
VERBOSE 
DEBUG 

Definition at line 33 of file SingleTopTChannelLeptonDQM_miniAOD.h.

33 { STANDARD, VERBOSE, DEBUG };

Constructor & Destructor Documentation

◆ MonitorEnsemble()

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

default contructor

Definition at line 31 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

32  : label_(label),
33  elecIso_(nullptr),
34  elecSelect_(nullptr),
35  pvSelect_(nullptr),
36  muonIso_(nullptr),
37  muonSelect_(nullptr),
38  jetIDSelect_(nullptr),
39  jetSelect(nullptr),
40  includeBTag_(false),
41  lowerEdge_(-1.),
42  upperEdge_(-1.),
43  logged_(0) {
44  // sources have to be given; this PSet is not optional
45  edm::ParameterSet sources = cfg.getParameter<edm::ParameterSet>("sources");
46  // muons_ = iC.consumes<edm::View<reco::PFCandidate> >(
47  // sources.getParameter<edm::InputTag>("muons"));
48 
49  muons_ = iC.consumes<edm::View<pat::Muon>>(sources.getParameter<edm::InputTag>("muons"));
50 
51  elecs_ = iC.consumes<edm::View<pat::Electron>>(sources.getParameter<edm::InputTag>("elecs"));
52  pvs_ = iC.consumes<edm::View<reco::Vertex>>(sources.getParameter<edm::InputTag>("pvs"));
53  jets_ = iC.consumes<edm::View<pat::Jet>>(sources.getParameter<edm::InputTag>("jets"));
54  for (edm::InputTag const& tag : sources.getParameter<std::vector<edm::InputTag>>("mets"))
55  mets_.push_back(iC.consumes<edm::View<pat::MET>>(tag));
56  // electronExtras are optional; they may be omitted or
57  // empty
58  if (cfg.existsAs<edm::ParameterSet>("elecExtras")) {
59  edm::ParameterSet elecExtras = cfg.getParameter<edm::ParameterSet>("elecExtras");
60  // select is optional; in case it's not found no
61  // selection will be applied
62  if (elecExtras.existsAs<std::string>("select")) {
63  elecSelect_ =
64  std::make_unique<StringCutObjectSelector<pat::Electron>>(elecExtras.getParameter<std::string>("select"));
65  }
66  // isolation is optional; in case it's not found no
67  // isolation will be applied
68  if (elecExtras.existsAs<std::string>("isolation")) {
69  elecIso_ =
70  std::make_unique<StringCutObjectSelector<pat::Electron>>(elecExtras.getParameter<std::string>("isolation"));
71  }
72 
73  if (elecExtras.existsAs<std::string>("rho")) {
74  rhoTag = elecExtras.getParameter<edm::InputTag>("rho");
75  }
76  // electronId is optional; in case it's not found the
77  // InputTag will remain empty
78  if (elecExtras.existsAs<edm::ParameterSet>("electronId")) {
79  edm::ParameterSet elecId = elecExtras.getParameter<edm::ParameterSet>("electronId");
81  eidCutValue_ = elecId.getParameter<double>("cutValue");
82  }
83  }
84  // pvExtras are opetional; they may be omitted or empty
85  if (cfg.existsAs<edm::ParameterSet>("pvExtras")) {
86  edm::ParameterSet pvExtras = cfg.getParameter<edm::ParameterSet>("pvExtras");
87  // select is optional; in case it's not found no
88  // selection will be applied
89  if (pvExtras.existsAs<std::string>("select")) {
90  pvSelect_ =
91  std::make_unique<StringCutObjectSelector<reco::Vertex>>(pvExtras.getParameter<std::string>("select"));
92  }
93  }
94  // muonExtras are optional; they may be omitted or empty
95  if (cfg.existsAs<edm::ParameterSet>("muonExtras")) {
96  edm::ParameterSet muonExtras = cfg.getParameter<edm::ParameterSet>("muonExtras");
97  // select is optional; in case it's not found no
98  // selection will be applied
99  if (muonExtras.existsAs<std::string>("select")) {
100  muonSelect_ =
101  std::make_unique<StringCutObjectSelector<pat::Muon>>(muonExtras.getParameter<std::string>("select"));
102  }
103  // isolation is optional; in case it's not found no
104  // isolation will be applied
105  if (muonExtras.existsAs<std::string>("isolation")) {
106  muonIso_ =
107  std::make_unique<StringCutObjectSelector<pat::Muon>>(muonExtras.getParameter<std::string>("isolation"));
108  }
109  }
110 
111  // jetExtras are optional; they may be omitted or
112  // empty
113  if (cfg.existsAs<edm::ParameterSet>("jetExtras")) {
114  edm::ParameterSet jetExtras = cfg.getParameter<edm::ParameterSet>("jetExtras");
115  // jetCorrector is optional; in case it's not found
116  // the InputTag will remain empty
117  if (jetExtras.existsAs<std::string>("jetCorrector")) {
118  jetCorrector_ = jetExtras.getParameter<std::string>("jetCorrector");
119  }
120  // read jetID information if it exists
121  if (jetExtras.existsAs<edm::ParameterSet>("jetID")) {
122  edm::ParameterSet jetID = jetExtras.getParameter<edm::ParameterSet>("jetID");
124  jetIDSelect_ =
125  std::make_unique<StringCutObjectSelector<reco::JetID>>(jetID.getParameter<std::string>("select"));
126  }
127  // select is optional; in case it's not found no
128  // selection will be applied (only implemented for
129  // CaloJets at the moment)
130  if (jetExtras.existsAs<std::string>("select")) {
131  jetSelect_ = jetExtras.getParameter<std::string>("select");
132  jetSelect = std::make_unique<StringCutObjectSelector<pat::Jet>>(jetSelect_);
133  }
134  }
135 
136  // triggerExtras are optional; they may be omitted or empty
137  if (cfg.existsAs<edm::ParameterSet>("triggerExtras")) {
138  edm::ParameterSet triggerExtras = cfg.getParameter<edm::ParameterSet>("triggerExtras");
140  triggerPaths_ = triggerExtras.getParameter<std::vector<std::string>>("paths");
141  }
142 
143  // massExtras is optional; in case it's not found no mass
144  // window cuts are applied for the same flavor monitor
145  // histograms
146  if (cfg.existsAs<edm::ParameterSet>("massExtras")) {
147  edm::ParameterSet massExtras = cfg.getParameter<edm::ParameterSet>("massExtras");
148  lowerEdge_ = massExtras.getParameter<double>("lowerEdge");
149  upperEdge_ = massExtras.getParameter<double>("upperEdge");
150  }
151 
152  // setup the verbosity level for booking histograms;
153  // per default the verbosity level will be set to
154  // STANDARD. This will also be the chosen level in
155  // the case when the monitoring PSet is not found
157  if (cfg.existsAs<edm::ParameterSet>("monitoring")) {
158  edm::ParameterSet monitoring = cfg.getParameter<edm::ParameterSet>("monitoring");
159  if (monitoring.getParameter<std::string>("verbosity") == "DEBUG")
160  verbosity_ = DEBUG;
161  if (monitoring.getParameter<std::string>("verbosity") == "VERBOSE")
163  if (monitoring.getParameter<std::string>("verbosity") == "STANDARD")
165  }
166  // and don't forget to do the histogram booking
167  directory_ = cfg.getParameter<std::string>("directory");
168  // book(ibooker);
169  }

References looper::cfg, DEBUG, directory_, eidCutValue_, singleTopDQM_cfi::elecExtras, elecIso_, elecs_, elecSelect_, electronId_, edm::ParameterSet::getParameter(), jetCorrector_, singleTopDQM_cfi::jetExtras, jetIDLabel_, jetIDSelect_, jets_, jetSelect, jetSelect_, 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, GlobalPosition_Frontier_DevDB_cff::tag, triggerPaths_, triggerTable_, upperEdge_, VERBOSE, and verbosity_.

◆ ~MonitorEnsemble()

SingleTopTChannelLepton_miniAOD::MonitorEnsemble::~MonitorEnsemble ( )
inline

default destructor

Definition at line 39 of file SingleTopTChannelLeptonDQM_miniAOD.h.

39 {};

Member Function Documentation

◆ book()

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

book histograms in subdirectory directory

Definition at line 171 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

171  {
172  // set up the current directory path
173  std::string current(directory_);
174  current += label_;
175  ibooker.setCurrentFolder(current);
176 
177  // determine number of bins for trigger monitoring
178  //unsigned int nPaths = triggerPaths_.size();
179 
180  // --- [STANDARD] --- //
181  // Run Number
182  //hists_["RunNumb_"] = ibooker.book1D("RunNumber", "Run Nr.", 1.e4, 1.5e5, 3.e5);
183  // instantaneous luminosity
184  //hists_["InstLumi_"] = ibooker.book1D("InstLumi", "Inst. Lumi.", 100, 0., 1.e3);
185  // number of selected primary vertices
186  hists_["pvMult_"] = ibooker.book1D("PvMult", "N_{good pvs}", 50, 0., 50.);
187  // pt of the leading muon
188  hists_["muonPt_"] = ibooker.book1D("MuonPt", "pt(#mu TightId, TightIso)", 40, 0., 200.);
189  // muon multiplicity before std isolation
190  hists_["muonMult_"] = ibooker.book1D("MuonMult", "N_{loose}(#mu)", 10, 0., 10.);
191  // muon multiplicity after std isolation
192  //hists_["muonMultIso_"] = ibooker.book1D("MuonMultIso",
193  // "N_{TightIso}(#mu)", 10, 0., 10.);
194 
195  hists_["muonMultTight_"] = ibooker.book1D("MuonMultTight", "N_{TightIso,TightId}(#mu)", 10, 0., 10.);
196 
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_{looseId}(e)", 10, 0., 10.);
201  // electron multiplicity after std isolation
202  //hists_["elecMultIso_"] = ibooker.book1D("ElecMultIso", "N_{Iso}(e)", 10, 0., 10.);
203  // multiplicity of jets with pt>20 (corrected to L2+L3)
204  hists_["jetMult_"] = ibooker.book1D("JetMult", "N_{30}(jet)", 10, 0., 10.);
205  hists_["jetLooseMult_"] = ibooker.book1D("JetLooseMult", "N_{30,loose}(jet)", 10, 0., 10.);
206 
207  // trigger efficiency estimates for single lepton triggers
208  //hists_["triggerEff_"] = ibooker.book1D("TriggerEff",
209  // "Eff(trigger)", nPaths, 0., nPaths);
210  // monitored trigger occupancy for single lepton triggers
211  //hists_["triggerMon_"] = ibooker.book1D("TriggerMon",
212  // "Mon(trigger)", nPaths, 0., nPaths);
213  // MET (calo)
214  hists_["slimmedMETs_"] = ibooker.book1D("slimmedMETs", "MET_{slimmed}", 40, 0., 200.);
215  // W mass estimate
216  hists_["massW_"] = ibooker.book1D("MassW", "M(W)", 60, 0., 300.);
217  // Top mass estimate
218  hists_["massTop_"] = ibooker.book1D("MassTop", "M(Top)", 50, 0., 500.);
219  // b-tagged Top mass
220  //hists_["massBTop_"] = ibooker.book1D("MassBTop", "M(Top, 1 b-tag)", 50, 0., 500.);
221 
222  // W mass transverse estimate mu
223  hists_["MTWm_"] = ibooker.book1D("MTWm", "M_{T}^{W}(#mu)", 60, 0., 300.);
224  // Top mass transverse estimate mu
225  hists_["mMTT_"] = ibooker.book1D("mMTT", "M_{T}^{t}(#mu)", 50, 0., 500.);
226 
227  // W mass transverse estimate e
228  hists_["MTWe_"] = ibooker.book1D("MTWe", "M_{T}^{W}(e)", 60, 0., 300.);
229  // Top mass transverse estimate e
230  hists_["eMTT_"] = ibooker.book1D("eMTT", "M_{T}^{t}(e)", 50, 0., 500.);
231 
232  // set bin labels for trigger monitoring
234 
235  if (verbosity_ == STANDARD)
236  return;
237 
238  // --- [VERBOSE] --- //
239  // eta of the leading muon
240  hists_["muonEta_"] = ibooker.book1D("MuonEta", "#eta(#mu TightId,TightIso)", 30, -3., 3.);
241  // relative isolation of the candidate muon (depending on the decay channel)
242  hists_["muonPhi_"] = ibooker.book1D("MuonPhi", "#phi(#mu TightId,TightIso)", 40, -4., 4.);
243  hists_["muonRelIso_"] = ibooker.book1D("MuonRelIso", "Iso_{Rel}(#mu TightId) (#Delta#beta Corrected)", 50, 0., 1.);
244 
245  // eta of the leading electron
246  hists_["elecEta_"] = ibooker.book1D("ElecEta", "#eta(e TightId, TightIso)", 30, -3., 3.);
247  hists_["elecPhi_"] = ibooker.book1D("ElecPhi", "#phi(e TightId, TightIso)", 40, -4., 4.);
248  // std isolation variable of the leading electron
249  hists_["elecRelIso_"] = ibooker.book1D("ElecRelIso", "Iso_{Rel}(e TightId)", 50, 0., 1.);
250 
251  hists_["elecMultTight_"] = ibooker.book1D("ElecMultTight", "N_{TightIso,TightId}(e)", 10, 0., 10.);
252 
253  // multiplicity of btagged jets (for track counting high efficiency) with
254  // pt(L2L3)>30
255  //hists_["jetMultBEff_"] = ibooker.book1D("JetMultBEff",
256  // "N_{30}(TCHE)", 10, 0., 10.);
257  // btag discriminator for track counting high efficiency for jets with
258  // pt(L2L3)>30
259  //hists_["jetBDiscEff_"] = ibooker.book1D("JetBDiscEff",
260  // "Disc_{TCHE}(jet)", 100, 0., 10.);
261  // eta of the 1. leading jet (corrected to L2+L3)
262  hists_["jet1Eta_"] = ibooker.book1D("Jet1Eta", "#eta_{30,loose}(jet1)", 60, -3., 3.);
263  // pt of the 1. leading jet (corrected to L2+L3)
264  hists_["jet1Pt_"] = ibooker.book1D("Jet1Pt", "pt_{30,loose}(jet1)", 60, 0., 300.);
265  // eta of the 2. leading jet (corrected to L2+L3)
266  hists_["jet2Eta_"] = ibooker.book1D("Jet2Eta", "#eta_{30,loose}(jet2)", 60, -3., 3.);
267  // pt of the 2. leading jet (corrected to L2+L3)
268  hists_["jet2Pt_"] = ibooker.book1D("Jet2Pt", "pt_{30,loose}(jet2)", 60, 0., 300.);
269  // eta of the 3. leading jet (corrected to L2+L3)
270  //hists_["jet3Eta_"] = ibooker.book1D("Jet3Eta", "#eta_{30,loose}(jet3)", 60, -3., 3.);
271  // pt of the 3. leading jet (corrected to L2+L3)
272  //hists_["jet3Pt_"] = ibooker.book1D("Jet3Pt", "pt_{30,loose}(jet3)", 60, 0., 300.);
273  // eta of the 4. leading jet (corrected to L2+L3)
274  //hists_["jet4Eta_"] = ibooker.book1D("Jet4Eta", "#eta_{30,loose}(jet4)", 60, -3., 3.);
275  // pt of the 4. leading jet (corrected to L2+L3)
276  //hists_["jet4Pt_"] = ibooker.book1D("Jet4Pt", "pt_{30,loose}(jet4)", 60, 0., 300.);
277  // MET (tc)
278  hists_["slimmedMETsNoHF_"] = ibooker.book1D("slimmedMETsNoHF", "MET_{slimmedNoHF}", 40, 0., 200.);
279  // MET (pflow)
280  hists_["slimmedMETsPuppi_"] = ibooker.book1D("slimmedMETsPuppi", "MET_{slimmedPuppi}", 40, 0., 200.);
281  // dz for muons (to suppress cosmis)
282  hists_["muonDelZ_"] = ibooker.book1D("MuonDelZ", "d_{z}(#mu)", 50, -25., 25.);
283  // dxy for muons (to suppress cosmics)
284  hists_["muonDelXY_"] = ibooker.book2D("MuonDelXY", "d_{xy}(#mu)", 50, -0.1, 0.1, 50, -0.1, 0.1);
285 
286  // set axes titles for dxy for muons
287  hists_["muonDelXY_"]->setAxisTitle("x [cm]", 1);
288  hists_["muonDelXY_"]->setAxisTitle("y [cm]", 2);
289 
290  if (verbosity_ == VERBOSE)
291  return;
292 
293  // --- [DEBUG] --- //
294  // charged hadron isolation component of the candidate muon (depending on the
295  // decay channel)
296  hists_["muonChHadIso_"] = ibooker.book1D("MuonChHadIsoComp", "ChHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
297  // neutral hadron isolation component of the candidate muon (depending on the
298  // decay channel)
299  hists_["muonNeHadIso_"] = ibooker.book1D("MuonNeHadIsoComp", "NeHad_{IsoComponent}(#mu TightId)", 50, 0., 5.);
300  // photon isolation component of the candidate muon (depending on the decay
301  // channel)
302  hists_["muonPhIso_"] = ibooker.book1D("MuonPhIsoComp", "Photon_{IsoComponent}(#mu TightId)", 50, 0., 5.);
303  // charged hadron isolation component of the candidate electron (depending on
304  // the decay channel)
305  hists_["elecChHadIso_"] = ibooker.book1D("ElectronChHadIsoComp", "ChHad_{IsoComponent}(e TightId)", 50, 0., 5.);
306  // neutral hadron isolation component of the candidate electron (depending on
307  // the decay channel)
308  hists_["elecNeHadIso_"] = ibooker.book1D("ElectronNeHadIsoComp", "NeHad_{IsoComponent}(e TightId)", 50, 0., 5.);
309  // photon isolation component of the candidate electron (depending on the
310  // decay channel)
311  hists_["elecPhIso_"] = ibooker.book1D("ElectronPhIsoComp", "Photon_{IsoComponent}(e TightId)", 50, 0., 5.);
312  // multiplicity of btagged jets (for track counting high purity) with
313  // pt(L2L3)>30
314  //hists_["jetMultBPur_"] = ibooker.book1D("JetMultBPur",
315  // "N_{30}(TCHP)", 10, 0., 10.);
316  // btag discriminator for track counting high purity
317  //hists_["jetBDiscPur_"] = ibooker.book1D("JetBDiscPur",
318  // "Disc_{TCHP}(Jet)", 100, 0., 10.);
319  // multiplicity of btagged jets (for simple secondary vertex) with pt(L2L3)>30
320  //hists_["jetMultBVtx_"] = ibooker.book1D("JetMultBVtx",
321  // "N_{30}(SSVHE)", 10, 0., 10.);
322  // btag discriminator for simple secondary vertex
323  //hists_["jetBDiscVtx_"] = ibooker.book1D("JetBDiscVtx",
324  // "Disc_{SSVHE}(Jet)", 35, -1., 6.);
325  // multiplicity for combined secondary vertex
326  hists_["jetMultBCSVM_"] = ibooker.book1D("JetMultBCSVM", "N_{30}(CSVM)", 10, 0., 10.);
327  // btag discriminator for combined secondary vertex
328  hists_["jetBCSV_"] = ibooker.book1D("JetDiscCSV", "BJet Disc_{CSV}(JET)", 100, -1., 2.);
329  // pt of the 1. leading jet (uncorrected)
330  //hists_["jet1PtRaw_"] = ibooker.book1D("Jet1PtRaw", "pt_{Raw}(jet1)", 60, 0., 300.);
331  // pt of the 2. leading jet (uncorrected)
332  //hists_["jet2PtRaw_"] = ibooker.book1D("Jet2PtRaw", "pt_{Raw}(jet2)", 60, 0., 300.);
333  // pt of the 3. leading jet (uncorrected)
334  //hists_["jet3PtRaw_"] = ibooker.book1D("Jet3PtRaw", "pt_{Raw}(jet3)", 60, 0., 300.);
335  // pt of the 4. leading jet (uncorrected)
336  //hists_["jet4PtRaw_"] = ibooker.book1D("Jet4PtRaw", "pt_{Raw}(jet4)", 60, 0., 300.);
337  // selected events
338  hists_["eventLogger_"] = ibooker.book2D("EventLogger", "Logged Events", 9, 0., 9., 10, 0., 10.);
339 
340  // set axes titles for selected events
341  hists_["eventLogger_"]->getTH1()->SetOption("TEXT");
342  hists_["eventLogger_"]->setBinLabel(1, "Run", 1);
343  hists_["eventLogger_"]->setBinLabel(2, "Block", 1);
344  hists_["eventLogger_"]->setBinLabel(3, "Event", 1);
345  hists_["eventLogger_"]->setBinLabel(4, "pt_{L2L3}(jet1)", 1);
346  hists_["eventLogger_"]->setBinLabel(5, "pt_{L2L3}(jet2)", 1);
347  hists_["eventLogger_"]->setBinLabel(6, "pt_{L2L3}(jet3)", 1);
348  hists_["eventLogger_"]->setBinLabel(7, "pt_{L2L3}(jet4)", 1);
349  hists_["eventLogger_"]->setBinLabel(8, "M_{W}", 1);
350  hists_["eventLogger_"]->setBinLabel(9, "M_{Top}", 1);
351  hists_["eventLogger_"]->setAxisTitle("logged evts", 2);
352  return;
353  }

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_.

◆ booked()

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

check if histogram was booked

Definition at line 63 of file SingleTopTChannelLeptonDQM_miniAOD.h.

63 { return hists_.find(histName) != hists_.end(); };

References hists_.

Referenced by fill().

◆ fill() [1/5]

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

fill monitor histograms with electronId and jetCorrections

Definition at line 355 of file SingleTopTChannelLeptonDQM_miniAOD.cc.

355  {
356  // fetch trigger event if configured such
358 
360  if (!event.getByToken(triggerTable_, triggerTable))
361  return;
362  }
363 
364  /*
365  ------------------------------------------------------------
366 
367  Primary Vertex Monitoring
368 
369  ------------------------------------------------------------
370  */
371  // fill monitoring plots for primary verices
373  if (!event.getByToken(pvs_, pvs))
374  return;
375  const reco::Vertex& pver = pvs->front();
376 
377  unsigned int pvMult = 0;
378  if (pvs.isValid()) {
379  for (edm::View<reco::Vertex>::const_iterator pv = pvs->begin(); pv != pvs->end(); ++pv) {
380  bool isGood =
381  (!(pv->isFake()) && (pv->ndof() > 4.0) && (abs(pv->z()) < 24.0) && (abs(pv->position().Rho()) < 2.0));
382  if (!isGood)
383  continue;
384  pvMult++;
385  }
386  //std::cout<<" npv "<<testn<<endl;
387  }
388 
389  fill("pvMult_", pvMult);
390 
391  /*
392  ------------------------------------------------------------
393 
394  Run and Inst. Luminosity information (Inst. Lumi. filled now with a dummy
395  value=5.0)
396 
397  ------------------------------------------------------------
398  */
399 
400  //if (!event.eventAuxiliary().run()) return;
401 
402  //fill("RunNumb_", event.eventAuxiliary().run());
403 
404  //double dummy = 5.;
405  //fill("InstLumi_", dummy);
406 
407  /*
408  ------------------------------------------------------------
409 
410  Electron Monitoring
411 
412  ------------------------------------------------------------
413  */
414 
415  // fill monitoring plots for electrons
417  if (!event.getByToken(elecs_, elecs))
418  return;
419 
420  edm::Handle<double> _rhoHandle;
421  event.getByLabel(rhoTag, _rhoHandle);
422  //if (!event.getByToken(elecs_, elecs)) return;
423 
424  // check availability of electron id
426  if (!electronId_.isUninitialized()) {
427  if (!event.getByToken(electronId_, electronId))
428  return;
429  }
430 
431  // loop electron collection
432  unsigned int eMultIso = 0, eMult = 0;
433  std::vector<const pat::Electron*> isoElecs;
434 
436 
437  for (edm::View<pat::Electron>::const_iterator elec = elecs->begin(); elec != elecs->end(); ++elec) {
438  if (true) { //loose id
439  if (!elecSelect_ || (*elecSelect_)(*elec)) {
440  double el_ChHadIso = elec->pfIsolationVariables().sumChargedHadronPt;
441  double el_NeHadIso = elec->pfIsolationVariables().sumNeutralHadronEt;
442  double el_PhIso = elec->pfIsolationVariables().sumPhotonEt;
443 
444  double rho = _rhoHandle.isValid() ? (float)(*_rhoHandle) : 0;
445  double absEta = abs(elec->superCluster()->eta());
446  double eA = 0;
447  if (absEta < 1.000)
448  eA = 0.1703;
449  else if (absEta < 1.479)
450  eA = 0.1715;
451  else if (absEta < 2.000)
452  eA = 0.1213;
453  else if (absEta < 2.200)
454  eA = 0.1230;
455  else if (absEta < 2.300)
456  eA = 0.1635;
457  else if (absEta < 2.400)
458  eA = 0.1937;
459  else if (absEta < 5.000)
460  eA = 0.2393;
461 
462  double el_pfRelIso = (el_ChHadIso + max(0., el_NeHadIso + el_PhIso - rho * eA)) / elec->pt();
463 
464  ++eMult;
465 
466  if (eMult == 1) {
467  fill("elecRelIso_", el_pfRelIso);
468  fill("elecChHadIso_", el_ChHadIso);
469  fill("elecNeHadIso_", el_NeHadIso);
470  fill("elecPhIso_", el_PhIso);
471  }
472  //loose Iso
473  //if(!((el_pfRelIso<0.0994 && absEta<1.479)||(el_pfRelIso<0.107 && absEta>1.479)))continue;
474 
475  //tight Iso
476  if (!((el_pfRelIso < 0.0588 && absEta < 1.479) || (el_pfRelIso < 0.0571 && absEta > 1.479)))
477  continue;
478  ++eMultIso;
479 
480  if (eMultIso == 1) {
481  // restrict to the leading electron
482  e = *elec;
483 
484  fill("elecPt_", elec->pt());
485  fill("elecEta_", elec->eta());
486  fill("elecPhi_", elec->phi());
487  }
488  }
489  }
490  }
491  //fill("elecMult_", eMult);
492  fill("elecMultTight_", eMultIso);
493 
494  /*
495  ------------------------------------------------------------
496 
497  Muon Monitoring
498 
499  ------------------------------------------------------------
500  */
501 
502  // fill monitoring plots for muons
503  unsigned int mMult = 0, mTight = 0;
504 
506 
507  pat::Muon mu;
508 
510 
511  if (!event.getByToken(muons_, muons))
512  return;
513 
514  for (edm::View<pat::Muon>::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
515  // restrict to globalMuons
516  if (muon->isGlobalMuon()) {
517  fill("muonDelZ_", muon->innerTrack()->vz()); // CB using inner track!
518  fill("muonDelXY_", muon->innerTrack()->vx(), muon->innerTrack()->vy());
519 
520  // apply preselection loose muon
521  if (!muonSelect_ || (*muonSelect_)(*muon)) {
522  //loose muon count
523  ++mMult;
524 
525  double chHadPt = muon->pfIsolationR04().sumChargedHadronPt;
526  double neHadEt = muon->pfIsolationR04().sumNeutralHadronEt;
527  double phoEt = muon->pfIsolationR04().sumPhotonEt;
528 
529  double pfRelIso = (chHadPt + max(0., neHadEt + phoEt - 0.5 * muon->pfIsolationR04().sumPUPt)) /
530  muon->pt(); // CB dBeta corrected iso!
531 
532  if (!(muon->isGlobalMuon() && muon->isPFMuon() && muon->globalTrack()->normalizedChi2() < 10. &&
533  muon->globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && muon->numberOfMatchedStations() > 1 &&
534  fabs(muon->muonBestTrack()->dxy(pver.position())) < 0.2 &&
535  fabs(muon->muonBestTrack()->dz(pver.position())) < 0.5 &&
536  muon->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 &&
537  muon->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5))
538  continue;
539 
540  if (mMult == 1) {
541  // restrict to leading muon
542  fill("muonRelIso_", pfRelIso);
543  fill("muonChHadIso_", chHadPt);
544  fill("muonNeHadIso_", neHadEt);
545  fill("muonPhIso_", phoEt);
546  fill("muonRelIso_", pfRelIso);
547  }
548 
549  if (!(pfRelIso < 0.15))
550  continue;
551 
552  ++mTight;
553 
554  //tight id
555  if (mTight == 1) {
556  // restrict to leading muon
557  mu = *muon;
558  fill("muonPt_", muon->pt());
559  fill("muonEta_", muon->eta());
560  fill("muonPhi_", muon->phi());
561  }
562  }
563  }
564  }
565  fill("muonMult_", mMult); //loose
566  fill("muonMultTight_", mTight); //tight id & iso
567 
568  /*
569  ------------------------------------------------------------
570 
571  Jet Monitoring
572 
573  ------------------------------------------------------------
574  */
575 
576  // loop jet collection
577  std::vector<pat::Jet> correctedJets;
578  std::vector<double> JetTagValues;
579  pat::Jet TaggedJetCand;
580  vector<double> bJetDiscVal;
581 
582  unsigned int mult = 0, loosemult = 0, multBCSVM = 0;
583 
585  if (!event.getByToken(jets_, jets)) {
586  return;
587  }
588 
589  for (edm::View<pat::Jet>::const_iterator jet = jets->begin(); jet != jets->end(); ++jet) {
590  // check jetID for calo jets
591  //unsigned int idx = jet - jets->begin();
592 
593  pat::Jet sel = *jet;
594 
595  if (jetSelect == nullptr)
596  jetSelect = std::make_unique<StringCutObjectSelector<pat::Jet>>(jetSelect_);
597 
598  if (!(*jetSelect)(sel))
599  continue;
600  // if (!jetSelect(sel)) continue;
601 
602  // prepare jet to fill monitor histograms
603  pat::Jet monitorJet = *jet;
604 
605  ++mult;
606 
607  if (monitorJet.chargedHadronEnergyFraction() > 0 && monitorJet.chargedMultiplicity() > 0 &&
608  monitorJet.chargedEmEnergyFraction() < 0.99 && monitorJet.neutralHadronEnergyFraction() < 0.99 &&
609  monitorJet.neutralEmEnergyFraction() < 0.99 &&
610  (monitorJet.chargedMultiplicity() + monitorJet.neutralMultiplicity()) > 1) {
611  correctedJets.push_back(monitorJet);
612  ++loosemult; // determine jet multiplicity
613 
614  fill("jetBCSV_",
615  monitorJet.bDiscriminator(
616  "pfCombinedInclusiveSecondaryVertexV2BJetTags")); //hard coded discriminator and value right now.
617  if (monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags") > 0.89) {
618  if (multBCSVM == 0) {
619  TaggedJetCand = monitorJet;
620  bJetDiscVal.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
621  } else if (multBCSVM == 1) {
622  bJetDiscVal.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
623  if (bJetDiscVal[1] > bJetDiscVal[0])
624  TaggedJetCand = monitorJet;
625  }
626 
627  ++multBCSVM;
628  }
629 
630  // Fill a vector with Jet b-tag WP for later M3+1tag calculation: CSV
631  // tagger
632  JetTagValues.push_back(monitorJet.bDiscriminator("pfCombinedInclusiveSecondaryVertexV2BJetTags"));
633  // }
634  // fill pt (raw or L2L3) for the leading four jets
635  if (loosemult == 1) {
636  //cout<<" jet id= "<<monitorJet.chargedHadronEnergyFraction()<<endl;
637 
638  fill("jet1Pt_", monitorJet.pt());
639  //fill("jet1PtRaw_", jet->pt());
640  fill("jet1Eta_", monitorJet.eta());
641  };
642  if (loosemult == 2) {
643  fill("jet2Pt_", monitorJet.pt());
644  //fill("jet2PtRaw_", jet->pt());
645  fill("jet2Eta_", monitorJet.eta());
646  }
647  }
648  }
649  fill("jetMult_", mult);
650  fill("jetLooseMult_", loosemult);
651  fill("jetMultBCSVM_", multBCSVM);
652 
653  /*
654  ------------------------------------------------------------
655 
656  MET Monitoring
657 
658  ------------------------------------------------------------
659  */
660 
661  // fill monitoring histograms for met
662 
663  pat::MET mET;
664 
665  for (std::vector<edm::EDGetTokenT<edm::View<pat::MET>>>::const_iterator met_ = mets_.begin(); met_ != mets_.end();
666  ++met_) {
668  if (!event.getByToken(*met_, met))
669  continue;
670  if (met->begin() != met->end()) {
671  unsigned int idx = met_ - mets_.begin();
672  if (idx == 0)
673  fill("slimmedMETs_", met->begin()->et());
674  if (idx == 1)
675  fill("slimmedMETsNoHF_", met->begin()->et());
676  if (idx == 2)
677  fill("slimmedMETsPuppi_", met->begin()->et());
678  }
679  }
680 
681  /*
682  ------------------------------------------------------------
683 
684  Event Monitoring
685 
686  ------------------------------------------------------------
687  */
688 
689  // fill W boson and top mass estimates
690 
691  Calculate_miniAOD eventKinematics(MAXJETS, WMASS);
692  double wMass = eventKinematics.massWBoson(correctedJets);
693  double topMass = eventKinematics.massTopQuark(correctedJets);
694  if (wMass >= 0 && topMass >= 0) {
695  fill("massW_", wMass);
696  fill("massTop_", topMass);
697  }
698 
699  // Fill M3 with Btag (CSV Tight) requirement
700 
701  // if (!includeBTag_) return;
702  //if (correctedJets.size() != JetTagValues.size()) return;
703  //double btopMass =
704  // eventKinematics.massBTopQuark(correctedJets, JetTagValues, 0.89); //hard coded CSVv2 value
705 
706  //if (btopMass >= 0) fill("massBTop_", btopMass);
707 
708  // fill plots for trigger monitoring
709  if ((lowerEdge_ == -1. && upperEdge_ == -1.) || (lowerEdge_ < wMass && wMass < upperEdge_)) {
711  fill(event, *triggerTable, "trigger", triggerPaths_);
712  if (logged_ <= hists_.find("eventLogger_")->second->getNbinsY()) {
713  // log runnumber, lumi block, event number & some
714  // more pysics infomation for interesting events
715  fill("eventLogger_", 0.5, logged_ + 0.5, event.eventAuxiliary().run());
716  fill("eventLogger_", 1.5, logged_ + 0.5, event.eventAuxiliary().luminosityBlock());
717  fill("eventLogger_", 2.5, logged_ + 0.5, event.eventAuxiliary().event());
718  if (!correctedJets.empty())
719  fill("eventLogger_", 3.5, logged_ + 0.5, correctedJets[0].pt());
720  if (correctedJets.size() > 1)
721  fill("eventLogger_", 4.5, logged_ + 0.5, correctedJets[1].pt());
722  if (correctedJets.size() > 2)
723  fill("eventLogger_", 5.5, logged_ + 0.5, correctedJets[2].pt());
724  if (correctedJets.size() > 3)
725  fill("eventLogger_", 6.5, logged_ + 0.5, correctedJets[3].pt());
726 
727  fill("eventLogger_", 7.5, logged_ + 0.5, wMass);
728  fill("eventLogger_", 8.5, logged_ + 0.5, topMass);
729  ++logged_;
730  }
731  }
732 
733  if (multBCSVM != 0 && mTight == 1) {
734  double mtW = eventKinematics.tmassWBoson(&mu, mET, TaggedJetCand);
735  fill("MTWm_", mtW);
736  double MTT = eventKinematics.tmassTopQuark(&mu, mET, TaggedJetCand);
737  fill("mMTT_", MTT);
738  }
739 
740  if (multBCSVM != 0 && eMultIso == 1) {
741  double mtW = eventKinematics.tmassWBoson(&e, mET, TaggedJetCand);
742  fill("MTWe_", mtW);
743  double MTT = eventKinematics.tmassTopQuark(&e, mET, TaggedJetCand);
744  fill("eMTT_", MTT);
745  }
746  }

References funct::abs(), pat::Jet::bDiscriminator(), pat::Jet::chargedEmEnergyFraction(), pat::Jet::chargedHadronEnergyFraction(), pat::Jet::chargedMultiplicity(), MillePedeFileConverter_cfg::e, singleTopDQM_cfi::elecs, elecs_, elecSelect_, topSingleLeptonDQM_PU_cfi::electronId, electronId_, reco::LeafCandidate::eta(), dqmMemoryStats::float, hists_, heavyIonCSV_trainingSettings::idx, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, singleTopDQM_cfi::jets, jets_, jetSelect, jetSelect_, logged_, lowerEdge_, Calculate_miniAOD::massTopQuark(), Calculate_miniAOD::massWBoson(), SiStripPI::max, SingleTopTChannelLepton_miniAOD::MAXJETS, BTaggingMonitor_cfi::met, mets_, amptDefaultParameters_cff::mu, VarParsing::mult, HLT_FULL_cff::muon, PDWG_BPHSkim_cff::muons, muons_, muonSelect_, pat::Jet::neutralEmEnergyFraction(), pat::Jet::neutralHadronEnergyFraction(), pat::Jet::neutralMultiplicity(), reco::Vertex::position(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), MetAnalyzer::pv(), FSQDQM_cfi::pvs, pvs_, rhoTag, EgammaValidation_Wenu_cff::sel, Calculate_miniAOD::tmassTopQuark(), Calculate_miniAOD::tmassWBoson(), triggerPaths_, triggerTable_, upperEdge_, trackerHitRTTI::vector, pseudoTop_cfi::wMass, and SingleTopTChannelLepton_miniAOD::WMASS.

Referenced by fill().

◆ 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 156 of file SingleTopTChannelLeptonDQM_miniAOD.h.

159  {
160  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
161  if (accept(event, triggerTable, monitorPath(labels[idx]))) {
162  fill(channel + "Mon_", idx + 0.5);
163  // take care to fill triggerMon_ before evts is being called
164  int evts = hists_.find(channel + "Mon_")->second->getBinContent(idx + 1);
165  double value = hists_.find(channel + "Eff_")->second->getBinContent(idx + 1);
166  fill(
167  channel + "Eff_", idx + 0.5, 1. / evts * (accept(event, triggerTable, selectionPath(labels[idx])) - value));
168  }
169  }
170  }

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

◆ fill() [3/5]

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

fill histogram if it had been booked before

Definition at line 65 of file SingleTopTChannelLeptonDQM_miniAOD.h.

65  {
66  if (booked(histName))
67  hists_.find(histName)->second->Fill(value);
68  };

References booked(), and hists_.

◆ fill() [4/5]

void SingleTopTChannelLepton_miniAOD::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 70 of file SingleTopTChannelLeptonDQM_miniAOD.h.

70  {
71  if (booked(histName))
72  hists_.find(histName)->second->Fill(xValue, yValue);
73  };

References booked(), and hists_.

◆ fill() [5/5]

void SingleTopTChannelLepton_miniAOD::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 75 of file SingleTopTChannelLeptonDQM_miniAOD.h.

75  {
76  if (booked(histName))
77  hists_.find(histName)->second->Fill(xValue, yValue, zValue);
78  };

References booked(), and hists_.

◆ monitorPath()

std::string SingleTopTChannelLepton_miniAOD::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 49 of file SingleTopTChannelLeptonDQM_miniAOD.h.

49 { return label.substr(label.find(':') + 1); };

References label.

Referenced by fill(), and triggerBinLabels().

◆ selectionPath()

std::string SingleTopTChannelLepton_miniAOD::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 52 of file SingleTopTChannelLeptonDQM_miniAOD.h.

52 { return label.substr(0, label.find(':')); };

References label.

Referenced by fill(), and triggerBinLabels().

◆ triggerBinLabels()

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

set configurable labels for trigger monitoring histograms

Definition at line 148 of file SingleTopTChannelLeptonDQM_miniAOD.h.

148  {
149  for (unsigned int idx = 0; idx < labels.size(); ++idx) {
150  hists_[channel + "Mon_"]->setBinLabel(idx + 1, "[" + monitorPath(labels[idx]) + "]", 1);
151  hists_[channel + "Eff_"]->setBinLabel(
152  idx + 1, "[" + selectionPath(labels[idx]) + "]|[" + monitorPath(labels[idx]) + "]", 1);
153  }
154  }

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

Referenced by book().

Member Data Documentation

◆ btagCSV_

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

Definition at line 133 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagCSVWP_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagCSVWP_
private

Definition at line 135 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagEff_

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

btag discriminator labels

Definition at line 133 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagEffWP_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagEffWP_
private

btag working points

Definition at line 135 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagPur_

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

Definition at line 133 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagPurWP_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagPurWP_
private

Definition at line 135 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagVtx_

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

Definition at line 133 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ btagVtxWP_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::btagVtxWP_
private

Definition at line 135 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ directory_

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::directory_
private

Definition at line 145 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by book(), and MonitorEnsemble().

◆ eidCutValue_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::eidCutValue_
private

Definition at line 103 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ elecIso_

std::unique_ptr<StringCutObjectSelector<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecIso_
private

extra isolation criterion on electron

Definition at line 105 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ elecs_

edm::EDGetTokenT<edm::View<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecs_
private

Definition at line 90 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ elecSelect_

std::unique_ptr<StringCutObjectSelector<pat::Electron> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecSelect_
private

extra selection on electrons

Definition at line 107 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ electronId_

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

electronId label

Definition at line 101 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ hists_

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

histogram container

Definition at line 143 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

◆ includeBTag_

bool SingleTopTChannelLepton_miniAOD::MonitorEnsemble::includeBTag_
private

include btag information or not to be determined from the cfg

Definition at line 131 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ jetCorrector_

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetCorrector_
private

jetCorrector

Definition at line 120 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ jetIDLabel_

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

jetID as an extra selection type

Definition at line 122 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ jetIDSelect_

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

extra jetID selection on calo jets

Definition at line 124 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ jets_

edm::EDGetTokenT<edm::View<pat::Jet> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jets_
private

input sources for monitoring

Definition at line 88 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ jetSelect

std::unique_ptr<StringCutObjectSelector<pat::Jet> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetSelect
private

Definition at line 128 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ jetSelect_

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetSelect_
private

extra selection on jets (here given as std::string as it depends on the the jet type, which selections are valid and which not)

Definition at line 127 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ label_

std::string SingleTopTChannelLepton_miniAOD::MonitorEnsemble::label_
private

◆ logged_

int SingleTopTChannelLepton_miniAOD::MonitorEnsemble::logged_
private

number of logged interesting events

Definition at line 140 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill().

◆ lowerEdge_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::lowerEdge_
private

mass window upper and lower edge

Definition at line 137 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ mets_

std::vector<edm::EDGetTokenT<edm::View<pat::MET> > > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::mets_
private

considers a vector of METs

Definition at line 86 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ muonIso_

std::unique_ptr<StringCutObjectSelector<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonIso_
private

extra isolation criterion on muon

Definition at line 114 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ muons_

edm::EDGetTokenT<edm::View<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muons_
private

Definition at line 89 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ muonSelect_

std::unique_ptr<StringCutObjectSelector<pat::Muon> > SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonSelect_
private

extra selection on muons

Definition at line 117 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ pvs_

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

Definition at line 91 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ pvSelect_

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

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

Definition at line 111 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by MonitorEnsemble().

◆ rhoTag

edm::InputTag SingleTopTChannelLepton_miniAOD::MonitorEnsemble::rhoTag
private

Definition at line 98 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ tmpConsumerBase

edm::EDConsumerBase SingleTopTChannelLepton_miniAOD::MonitorEnsemble::tmpConsumerBase
private

Definition at line 144 of file SingleTopTChannelLeptonDQM_miniAOD.h.

◆ triggerPaths_

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

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

Definition at line 96 of file SingleTopTChannelLeptonDQM_miniAOD.h.

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

◆ triggerTable_

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

trigger table

Definition at line 93 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ upperEdge_

double SingleTopTChannelLepton_miniAOD::MonitorEnsemble::upperEdge_
private

Definition at line 137 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by fill(), and MonitorEnsemble().

◆ verbosity_

Level SingleTopTChannelLepton_miniAOD::MonitorEnsemble::verbosity_
private

verbosity level for booking

Definition at line 78 of file SingleTopTChannelLeptonDQM_miniAOD.h.

Referenced by book(), and MonitorEnsemble().

SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetSelect
std::unique_ptr< StringCutObjectSelector< pat::Jet > > jetSelect
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:128
SummaryClient_cfi.labels
labels
Definition: SummaryClient_cfi.py:61
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecSelect_
std::unique_ptr< StringCutObjectSelector< pat::Electron > > elecSelect_
extra selection on electrons
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:107
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::pvs_
edm::EDGetTokenT< edm::View< reco::Vertex > > pvs_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:91
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
pat::Jet::bDiscriminator
float bDiscriminator(const std::string &theLabel) const
-— methods for accessing b-tagging info -—
muon
Definition: MuonCocktails.h:17
SingleTopTChannelLepton_miniAOD::WMASS
static const double WMASS
Definition: SingleTopTChannelLeptonDQM_miniAOD.cc:29
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
amptDefaultParameters_cff.mu
mu
Definition: amptDefaultParameters_cff.py:16
edm::EDGetTokenT
Definition: EDGetToken.h:33
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::triggerBinLabels
void triggerBinLabels(std::string channel, const std::vector< std::string > labels)
set configurable labels for trigger monitoring histograms
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:148
ZElectronSkim_cff.rho
rho
Definition: ZElectronSkim_cff.py:38
singleTopDQM_cfi.elecExtras
elecExtras
Definition: singleTopDQM_cfi.py:52
pat::Jet::neutralHadronEnergyFraction
float neutralHadronEnergyFraction() const
neutralHadronEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:404
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::eidCutValue_
double eidCutValue_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:103
edm::EDGetTokenT::isUninitialized
constexpr bool isUninitialized() const noexcept
Definition: EDGetToken.h:99
CalibrationSummaryClient_cfi.sources
sources
Definition: CalibrationSummaryClient_cfi.py:23
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:114
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::electronId_
edm::EDGetTokenT< edm::ValueMap< float > > electronId_
electronId label
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:101
singleTopDQM_cfi.jets
jets
Definition: singleTopDQM_cfi.py:42
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
pat::Muon
Analysis-level muon class.
Definition: Muon.h:51
topSingleLeptonDQM_PU_cfi.electronId
electronId
when omitted electron plots will be filled w/o cut on electronId
Definition: topSingleLeptonDQM_PU_cfi.py:39
edm::Handle< edm::TriggerResults >
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::triggerTable_
edm::EDGetTokenT< edm::TriggerResults > triggerTable_
trigger table
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:93
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::hists_
std::map< std::string, MonitorElement * > hists_
histogram container
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:143
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jets_
edm::EDGetTokenT< edm::View< pat::Jet > > jets_
input sources for monitoring
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:88
BTaggingMonitor_cfi.met
met
Definition: BTaggingMonitor_cfi.py:84
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::includeBTag_
bool includeBTag_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:131
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::booked
bool booked(const std::string histName) const
check if histogram was booked
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:63
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:30
pat::Jet
Analysis-level calorimeter jet class.
Definition: Jet.h:77
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::VERBOSE
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:33
singleTopDQM_cfi.pvExtras
pvExtras
Definition: singleTopDQM_cfi.py:49
GlobalPosition_Frontier_DevDB_cff.tag
tag
Definition: GlobalPosition_Frontier_DevDB_cff.py:11
HLT_FULL_cff.muon
muon
Definition: HLT_FULL_cff.py:11773
singleTopDQM_cfi.muonExtras
muonExtras
Definition: singleTopDQM_cfi.py:56
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonSelect_
std::unique_ptr< StringCutObjectSelector< pat::Muon > > muonSelect_
extra selection on muons
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:117
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:55
singleTopDQM_cfi.elecs
elecs
Definition: singleTopDQM_cfi.py:41
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetSelect_
std::string jetSelect_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:127
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::selectionPath
std::string selectionPath(const std::string &label) const
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:52
pat::Jet::chargedMultiplicity
int chargedMultiplicity() const
chargedMultiplicity
Definition: Jet.h:746
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::STANDARD
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:33
singleTopDQM_cfi.massExtras
massExtras
Definition: singleTopDQM_cfi.py:70
edm::View
Definition: CaloClusterFwd.h:14
edm::ParameterSet
Definition: ParameterSet.h:47
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
pat::Jet::neutralEmEnergyFraction
float neutralEmEnergyFraction() const
neutralEmEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:412
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::lowerEdge_
double lowerEdge_
mass window upper and lower edge
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:137
Calculate_miniAOD
Definition: TopDQMHelpers.h:70
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::directory_
std::string directory_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:145
pat::Jet::neutralMultiplicity
int neutralMultiplicity() const
neutralMultiplicity
Definition: Jet.h:466
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetCorrector_
std::string jetCorrector_
jetCorrector
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:120
value
Definition: value.py:1
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::label_
std::string label_
instance label
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:84
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
pat::MET
Analysis-level MET class.
Definition: MET.h:40
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecIso_
std::unique_ptr< StringCutObjectSelector< pat::Electron > > elecIso_
extra isolation criterion on electron
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:105
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetIDLabel_
edm::EDGetTokenT< reco::JetIDValueMap > jetIDLabel_
jetID as an extra selection type
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:122
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::logged_
int logged_
number of logged interesting events
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:140
FSQDQM_cfi.pvs
pvs
Definition: FSQDQM_cfi.py:12
looper.cfg
cfg
Definition: looper.py:297
pat::Jet::chargedEmEnergyFraction
float chargedEmEnergyFraction() const
chargedEmEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:408
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::DEBUG
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:33
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muonIso_
std::unique_ptr< StringCutObjectSelector< pat::Muon > > muonIso_
extra isolation criterion on muon
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:114
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::rhoTag
edm::InputTag rhoTag
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:98
singleTopDQM_cfi.jetExtras
jetExtras
Definition: singleTopDQM_cfi.py:60
pat::Jet::chargedHadronEnergyFraction
float chargedHadronEnergyFraction() const
chargedHadronEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:400
singleTopDQM_cfi.monitoring
monitoring
Definition: singleTopDQM_cfi.py:46
metsig::jet
Definition: SignAlgoResolutions.h:47
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::monitorPath
std::string monitorPath(const std::string &label) const
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:49
edm::ValueMap< float >
dqm::implementation::IBooker::book2D
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:177
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::verbosity_
Level verbosity_
verbosity level for booking
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:78
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::muons_
edm::EDGetTokenT< edm::View< pat::Muon > > muons_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:89
pseudoTop_cfi.wMass
wMass
Definition: pseudoTop_cfi.py:12
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::upperEdge_
double upperEdge_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:137
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::jetIDSelect_
std::unique_ptr< StringCutObjectSelector< reco::JetID > > jetIDSelect_
extra jetID selection on calo jets
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:124
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::elecs_
edm::EDGetTokenT< edm::View< pat::Electron > > elecs_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:90
VarParsing.mult
mult
Definition: VarParsing.py:659
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::mets_
std::vector< edm::EDGetTokenT< edm::View< pat::MET > > > mets_
considers a vector of METs
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:86
pat::Electron
Analysis-level electron class.
Definition: Electron.h:51
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::triggerPaths_
std::vector< std::string > triggerPaths_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:96
EgammaValidation_Wenu_cff.sel
sel
Definition: EgammaValidation_Wenu_cff.py:33
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::pvSelect_
std::unique_ptr< StringCutObjectSelector< reco::Vertex > > pvSelect_
Definition: SingleTopTChannelLeptonDQM_miniAOD.h:111
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
event
Definition: event.py:1
edm::InputTag
Definition: InputTag.h:15
SingleTopTChannelLepton_miniAOD::MAXJETS
static const unsigned int MAXJETS
Definition: SingleTopTChannelLeptonDQM_miniAOD.cc:26
label
const char * label
Definition: PFTauDecayModeTools.cc:11
edm::TriggerResults
Definition: TriggerResults.h:35
reco::Vertex
Definition: Vertex.h:35
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
SingleTopTChannelLepton_miniAOD::MonitorEnsemble::fill
void fill(const edm::Event &event, const edm::EventSetup &setup)
fill monitor histograms with electronId and jetCorrections
Definition: SingleTopTChannelLeptonDQM_miniAOD.cc:355