|
|
#include <QcdUeDQM.h>
|
void | book1D (DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, bool sumw2=true, bool sbox=true) |
|
void | bookProfile (DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, double y1, double y2, bool sumw2=true, bool sbox=true) |
|
void | fill1D (std::vector< MonitorElement * > &mes, double val, double w=1.) |
|
void | fill1D (std::vector< TH1F * > &hs, double val, double w=1.) |
|
void | fill2D (std::vector< MonitorElement * > &mes, double valx, double valy, double w=1.) |
|
void | fill2D (std::vector< TH2F * > &hs, double valx, double valy, double w=1.) |
|
void | fill3D (std::vector< TH3F * > &hs, int gbin, double w=1.) |
|
void | fillChargedJetSpectra (const edm::Handle< reco::TrackJetCollection > trackJets) |
|
void | fillHltBits (const edm::Event &iEvent, const edm::EventSetup &iSetup) |
|
void | fillProfile (std::vector< MonitorElement * > &mes, double valx, double valy, double w=1.) |
|
void | fillProfile (std::vector< TProfile * > &hs, double valx, double valy, double w=1.) |
|
void | fillpTMaxRelated (const std::vector< const reco::Track * > &track) |
|
void | fillUE_with_ChargedJets (const std::vector< const reco::Track * > &track, const edm::Handle< reco::TrackJetCollection > &trackJets) |
|
void | fillUE_with_MaxpTtrack (const std::vector< const reco::Track * > &track) |
|
bool | fillVtxPlots (const reco::BeamSpot *bs, const edm::Handle< reco::VertexCollection > vtxColl) |
|
template<typename TYPE > |
void | getProduct (const std::string name, edm::Handle< TYPE > &prod, const edm::Event &event) const |
|
template<typename TYPE > |
bool | getProductSafe (const std::string name, edm::Handle< TYPE > &prod, const edm::Event &event) const |
|
void | setLabel1D (std::vector< MonitorElement * > &mes) |
|
bool | trackSelection (const reco::Track &trk, const reco::BeamSpot *bs, const reco::Vertex &vtx, int sizevtx) |
|
Definition at line 53 of file QcdUeDQM.h.
◆ container
◆ QcdUeDQM()
◆ ~QcdUeDQM()
◆ analyze()
Reimplemented from DQMEDAnalyzer.
Definition at line 636 of file QcdUeDQM.cc.
668 std::sort(tracks_sort.begin(), tracks_sort.end(),
PtSorter());
681 for (reco::TrackCollection::const_iterator Trk = tracks_sort.begin(); Trk != tracks_sort.end(); ++Trk) {
690 if (!trkJets->empty())
References pwdgSkimBPark_cfi::beamSpot, bsLabel_, caloJetLabel_, chargedJetLabel_, fill1D(), fillChargedJetSpectra(), fillHltBits(), fillpTMaxRelated(), fillUE_with_ChargedJets(), fillUE_with_MaxpTtrack(), fillVtxPlots(), hNevts_, iEvent, isHltConfigSuccessful_, selected_, trackLabel_, PDWG_EXOHSCP_cff::tracks, trackSelection(), vtx1, and vtxLabel_.
◆ book1D()
void QcdUeDQM::book1D |
( |
DQMStore::IBooker & |
, |
|
|
std::vector< MonitorElement * > & |
mes, |
|
|
const std::string & |
name, |
|
|
const std::string & |
title, |
|
|
int |
nx, |
|
|
double |
x1, |
|
|
double |
x2, |
|
|
bool |
sumw2 = true , |
|
|
bool |
sbox = true |
|
) |
| |
|
private |
Definition at line 696 of file QcdUeDQM.cc.
714 TH1 *h1 =
e->getTH1();
716 if (0 == h1->GetSumw2N()) {
References dqm::implementation::IBooker::book1D(), MillePedeFileConverter_cfg::e, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, hltTrgUsedNames_, mps_fire::i, Skims_PA_cff::name, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, overlapproblemtsosanalyzer_cfi::title, testProducerWithPsetDescEmpty_cfi::x1, and testProducerWithPsetDescEmpty_cfi::x2.
Referenced by bookHistograms().
◆ bookHistograms()
Implements DQMEDAnalyzer.
Definition at line 164 of file QcdUeDQM.cc.
166 iBooker.setCurrentFolder(
"Physics/QcdUe");
170 const double x1 = -0.5;
171 const double x2 =
Nx - 0.5;
181 book1D(iBooker,
hNevts_,
"hNevts",
"number of events", 2, 0, 2);
182 book1D(iBooker,
hNtrackerLayer_,
"hNtrackerLayer",
"number of tracker layers;multiplicity", 20, -0.5, 19.5);
186 "hNtrackerStripPixelLayer",
187 "number of strip + pixel layers;multiplicity",
192 book1D(iBooker,
hTrkPt_,
"hTrkPt",
"pT of all tracks", 50, 0., 50.);
193 book1D(iBooker,
hTrkEta_,
"hTrkEta",
"eta of all tracks", 40, -4., 4.);
194 book1D(iBooker,
hTrkPhi_,
"hTrkPhi",
"phi of all tracks", 40, -4., 4.);
197 "hRatioDxySigmaDxyBS",
198 "ratio of transverse impact parameter and its significance wrt beam spot",
204 "hRatioDxySigmaDxyPV",
205 "ratio of transverse impact parameter and its significance wrt PV",
212 "ratio of longitudinal impact parameter and its significance wrt beam "
220 "ratio of longitudinal impact parameter and its significance wrt PV",
227 book1D(iBooker,
hNgoodTrk_,
"hNgoodTrk",
"number of good tracks", 50, -0.5, 49.5);
229 book1D(iBooker,
hGoodTrkPt500_,
"hGoodTrkPt500",
"pT of all good tracks with pT > 500 MeV", 50, 0., 50.);
230 book1D(iBooker,
hGoodTrkEta500_,
"hGoodTrkEta500",
"eta of all good tracks pT > 500 MeV", 40, -4., 4.);
231 book1D(iBooker,
hGoodTrkPhi500_,
"hGoodTrkPhi500",
"phi of all good tracks pT > 500 MeV", 40, -4., 4.);
233 book1D(iBooker,
hGoodTrkPt900_,
"hGoodTrkPt900",
"pT of all good tracks with pT > 900 MeV", 50, 0., 50.);
234 book1D(iBooker,
hGoodTrkEta900_,
"hGoodTrkEta900",
"eta of all good tracks pT > 900 MeV", 40, -4., 4.);
235 book1D(iBooker,
hGoodTrkPhi900_,
"hGoodTrkPhi900",
"phi of all good tracks pT > 900 MeV", 40, -4., 4.);
238 book1D(iBooker,
hVertex_z_,
"hVertex_z",
"z position of vertex; z[cm]", 200, -50, 50);
239 book1D(iBooker,
hVertex_y_,
"hVertex_y",
"y position of vertex; y[cm]", 100, -5, 5);
240 book1D(iBooker,
hVertex_x_,
"hVertex_x",
"x position of vertex; x[cm]", 100, -5, 5);
243 book1D(iBooker,
hVertex_z_bs_,
"hVertex_z_bs",
"z position of vertex from beamspot; z[cm]", 200, -50, 50);
245 book1D(iBooker,
hBeamSpot_z_,
"hBeamSpot_z",
"z position of beamspot; z[cm]", 100, -20, 20);
246 book1D(iBooker,
hBeamSpot_y_,
"hBeamSpot_y",
"y position of beamspot; y[cm]", 50, -10, 10);
247 book1D(iBooker,
hBeamSpot_x_,
"hBeamSpot_x",
"x position of beamspot; x[cm]", 50, -10, 10);
251 const double x1 = 0.0;
252 const double x2 = 50.0;
255 "hLeadingTrack_pTSpectrum",
256 "pT spectrum of leading track;pT(GeV/c)",
262 "hLeadingChargedJet_pTSpectrum",
263 "pT spectrum of leading track jet;pT(GeV/c)",
271 const double x1 = -4.;
272 const double x2 = 4.;
275 "hLeadingTrack_phiSpectrum",
276 "#phi spectrum of leading track;#phi",
282 "hLeadingChargedJet_phiSpectrum",
283 "#phi spectrum of leading track jet;#phi",
291 const double x1 = -4.;
292 const double x2 = 4.;
295 "hLeadingTrack_etaSpectrum",
296 "#eta spectrum of leading track;#eta",
302 "hLeadingChargedJet_etaSpectrum",
303 "#eta spectrum of leading track jet;#eta",
311 const double x1 = 0.0;
312 const double x2 = 75.0;
313 const double y1 = 0.;
314 const double y2 = 10.;
317 "hdNdEtadPhi_pTMax_Toward500",
318 "Average number of tracks (pT > 500 MeV) in toward region vs "
319 "leading track pT;pT(GeV/c);dN/d#eta d#phi",
329 "hdNdEtadPhi_pTMax_Transverse500",
330 "Average number of tracks (pT > 500 MeV) in transverse region "
331 "vs leading track pT;pT(GeV/c);dN/d#eta d#phi",
341 "hdNdEtadPhi_pTMax_Away500",
342 "Average number of tracks (pT > 500 MeV) in away region vs "
343 "leading track pT;pT(GeV/c);dN/d#eta d#phi",
353 "hdNdEtadPhi_trackJet_Toward500",
354 "Average number of tracks (pT > 500 MeV) in toward region vs "
355 "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
363 "hdNdEtadPhi_trackJet_Transverse500",
364 "Average number of tracks (pT > 500 MeV) in transverse region "
365 "vs leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
375 "hdNdEtadPhi_trackJet_Away500",
376 "Average number of tracks (pT > 500 MeV) in away region vs "
377 "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
388 "hpTSumdEtadPhi_pTMax_Toward500",
389 "Average number of tracks (pT > 500 MeV) in toward region vs "
390 "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
400 "hpTSumdEtadPhi_pTMax_Transverse500",
401 "Average number of tracks (pT > 500 MeV) in transverse region "
402 "vs leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
412 "hpTSumdEtadPhi_pTMax_Away500",
413 "Average number of tracks (pT > 500 MeV) in away region vs "
414 "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
424 "hpTSumdEtadPhi_trackJet_Toward500",
425 "Average number of tracks (pT > 500 MeV) in toward region vs "
426 "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
436 "hpTSumdEtadPhi_trackJet_Transverse500",
437 "Average number of tracks (pT > 500 MeV) in transverse region "
438 "vs leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
448 "hpTSumdEtadPhi_trackJet_Away500",
449 "Average number of tracks (pT > 500 MeV) in away region vs "
450 "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
461 "hdNdEtadPhi_pTMax_Toward900",
462 "Average number of tracks (pT > 900 MeV) in toward region vs "
463 "leading track pT;pT(GeV/c);dN/d#eta d#phi",
473 "hdNdEtadPhi_pTMax_Transverse900",
474 "Average number of tracks (pT > 900 MeV) in transverse region "
475 "vs leading track pT;pT(GeV/c);dN/d#eta d#phi",
485 "hdNdEtadPhi_pTMax_Away900",
486 "Average number of tracks (pT > 900 MeV) in away region vs "
487 "leading track pT;pT(GeV/c);dN/d#eta d#phi",
497 "hdNdEtadPhi_trackJet_Toward900",
498 "Average number of tracks (pT > 900 MeV) in toward region vs "
499 "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
507 "hdNdEtadPhi_trackJet_Transverse900",
508 "Average number of tracks (pT > 900 MeV) in transverse region "
509 "vs leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
519 "hdNdEtadPhi_trackJet_Away900",
520 "Average number of tracks (pT > 900 MeV) in away region vs "
521 "leading track jet pT;pT(GeV/c);dN/d#eta d#phi",
532 "hpTSumdEtadPhi_pTMax_Toward900",
533 "Average number of tracks (pT > 900 MeV) in toward region vs "
534 "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
544 "hpTSumdEtadPhi_pTMax_Transverse900",
545 "Average number of tracks (pT > 900 MeV) in transverse region "
546 "vs leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
556 "hpTSumdEtadPhi_pTMax_Away900",
557 "Average number of tracks (pT > 900 MeV) in away region vs "
558 "leading track pT;pT(GeV/c);dpTSum/d#eta d#phi",
568 "hpTSumdEtadPhi_trackJet_Toward900",
569 "Average number of tracks (pT > 900 MeV) in toward region vs "
570 "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
580 "hpTSumdEtadPhi_trackJet_Transverse900",
581 "Average number of tracks (pT > 900 MeV) in transverse region "
582 "vs leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
592 "hpTSumdEtadPhi_trackJet_Away900",
593 "Average number of tracks (pT > 900 MeV) in away region vs "
594 "leading track jet pT;pT(GeV/c);dpTSum/d#eta d#phi",
606 const double x1 = 0.0;
607 const double x2 = 20.0;
614 const double x1 = -180.0;
615 const double x2 = 180.0;
619 "hdPhi_maxpTTrack_tracks",
620 "delta phi between leading tracks and other "
621 "tracks;#Delta#phi(leading track-track)",
627 "hdPhi_chargedJet_tracks",
628 "delta phi between leading charged jet and "
629 "tracks;#Delta#phi(leading charged jet-track)",
References book1D(), dqm::implementation::IBooker::book2D(), bookProfile(), dqm::impl::MonitorElement::getTH1(), h2TrigCorr_, hBeamSpot_x_, hBeamSpot_y_, hBeamSpot_z_, hChargedJetMulti_, hdNdEtadPhi_pTMax_Away500_, hdNdEtadPhi_pTMax_Away900_, hdNdEtadPhi_pTMax_Toward500_, hdNdEtadPhi_pTMax_Toward900_, hdNdEtadPhi_pTMax_Transverse500_, hdNdEtadPhi_pTMax_Transverse900_, hdNdEtadPhi_trackJet_Away500_, hdNdEtadPhi_trackJet_Away900_, hdNdEtadPhi_trackJet_Toward500_, hdNdEtadPhi_trackJet_Toward900_, hdNdEtadPhi_trackJet_Transverse500_, hdNdEtadPhi_trackJet_Transverse900_, hdPhi_chargedJet_tracks_, hdPhi_maxpTTrack_tracks_, hGoodTrkEta500_, hGoodTrkEta900_, hGoodTrkPhi500_, hGoodTrkPhi900_, hGoodTrkPt500_, hGoodTrkPt900_, hLeadingChargedJet_etaSpectrum_, hLeadingChargedJet_phiSpectrum_, hLeadingChargedJet_pTSpectrum_, hLeadingTrack_etaSpectrum_, hLeadingTrack_phiSpectrum_, hLeadingTrack_pTSpectrum_, hltTrgUsedNames_, hNevts_, hNgoodTrk_, hNtrackerLayer_, hNtrackerPixelLayer_, hNtrackerStripPixelLayer_, hNvertices_, hpTSumdEtadPhi_pTMax_Away500_, hpTSumdEtadPhi_pTMax_Away900_, hpTSumdEtadPhi_pTMax_Toward500_, hpTSumdEtadPhi_pTMax_Toward900_, hpTSumdEtadPhi_pTMax_Transverse500_, hpTSumdEtadPhi_pTMax_Transverse900_, hpTSumdEtadPhi_trackJet_Away500_, hpTSumdEtadPhi_trackJet_Away900_, hpTSumdEtadPhi_trackJet_Toward500_, hpTSumdEtadPhi_trackJet_Toward900_, hpTSumdEtadPhi_trackJet_Transverse500_, hpTSumdEtadPhi_trackJet_Transverse900_, hRatioDxySigmaDxyBS_, hRatioDxySigmaDxyPV_, hRatioDzSigmaDzBS_, hRatioDzSigmaDzPV_, hRatioPtErrorPt_, hTrkChi2_, hTrkEta_, hTrkNdof_, hTrkPhi_, hTrkPt_, hVertex_ndof_, hVertex_rho_, hVertex_x_, hVertex_y_, hVertex_z_, hVertex_z_bs_, mps_fire::i, PixelTestBeamValidation_cfi::Nx, dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.
◆ bookProfile()
void QcdUeDQM::bookProfile |
( |
DQMStore::IBooker & |
, |
|
|
std::vector< MonitorElement * > & |
mes, |
|
|
const std::string & |
name, |
|
|
const std::string & |
title, |
|
|
int |
nx, |
|
|
double |
x1, |
|
|
double |
x2, |
|
|
double |
y1, |
|
|
double |
y2, |
|
|
bool |
sumw2 = true , |
|
|
bool |
sbox = true |
|
) |
| |
|
private |
Definition at line 726 of file QcdUeDQM.cc.
References dqm::implementation::IBooker::bookProfile(), MillePedeFileConverter_cfg::e, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, hltTrgUsedNames_, mps_fire::i, Skims_PA_cff::name, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, overlapproblemtsosanalyzer_cfi::title, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.
Referenced by bookHistograms().
◆ dqmBeginRun()
Reimplemented from DQMEDAnalyzer.
Definition at line 106 of file QcdUeDQM.cc.
109 bool isHltCfgChange =
false;
References gather_cfg::cout, CP, newFWLiteAna::found, hltConfig, hltProcNames_, hltResName_, hltTrgBits_, hltTrgDeci_, hltTrgNames_, hltTrgUsedNames_, hltUsedResName_, mps_fire::i, HLTConfigProvider::init(), isHltConfigSuccessful_, dqmiolumiharvest::j, writedatasetfile::run, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().
◆ fill1D() [1/2]
void QcdUeDQM::fill1D |
( |
std::vector< MonitorElement * > & |
mes, |
|
|
double |
val, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fill1D() [2/2]
void QcdUeDQM::fill1D |
( |
std::vector< TH1F * > & |
hs, |
|
|
double |
val, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
Definition at line 754 of file QcdUeDQM.cc.
757 for (
size_t i = 0;
i < hs.size(); ++
i) {
760 hs.at(
i)->Fill(
val,
w);
References hltTrgDeci_, mps_fire::i, heppy_batch::val, and w.
Referenced by analyze(), fillChargedJetSpectra(), fillpTMaxRelated(), fillUE_with_ChargedJets(), fillUE_with_MaxpTtrack(), fillVtxPlots(), and trackSelection().
◆ fill2D() [1/2]
void QcdUeDQM::fill2D |
( |
std::vector< MonitorElement * > & |
mes, |
|
|
double |
valx, |
|
|
double |
valy, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fill2D() [2/2]
void QcdUeDQM::fill2D |
( |
std::vector< TH2F * > & |
hs, |
|
|
double |
valx, |
|
|
double |
valy, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fill3D()
void QcdUeDQM::fill3D |
( |
std::vector< TH3F * > & |
hs, |
|
|
int |
gbin, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fillChargedJetSpectra()
◆ fillHltBits()
Definition at line 1251 of file QcdUeDQM.cc.
1286 if (tbit < triggerResultsHLT->
size()) {
References edm::HLTGlobalStatus::accept(), dqm::impl::MonitorElement::Fill(), getProduct(), h2TrigCorr_, hltTrgBits_, hltTrgDeci_, hltUsedResName_, mps_fire::i, iEvent, dqmiolumiharvest::j, and findQualityFiles::size.
Referenced by analyze().
◆ fillProfile() [1/2]
void QcdUeDQM::fillProfile |
( |
std::vector< MonitorElement * > & |
mes, |
|
|
double |
valx, |
|
|
double |
valy, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fillProfile() [2/2]
void QcdUeDQM::fillProfile |
( |
std::vector< TProfile * > & |
hs, |
|
|
double |
valx, |
|
|
double |
valy, |
|
|
double |
w = 1. |
|
) |
| |
|
private |
◆ fillpTMaxRelated()
void QcdUeDQM::fillpTMaxRelated |
( |
const std::vector< const reco::Track * > & |
track | ) |
|
|
private |
Definition at line 957 of file QcdUeDQM.cc.
959 if (!
track.empty()) {
964 for (
size_t i = 0;
i <
track.size();
i++) {
References PVValHelper::eta, fill1D(), hGoodTrkEta500_, hGoodTrkEta900_, hGoodTrkPhi500_, hGoodTrkPhi900_, hGoodTrkPt500_, hGoodTrkPt900_, hLeadingTrack_etaSpectrum_, hLeadingTrack_phiSpectrum_, hLeadingTrack_pTSpectrum_, hNgoodTrk_, mps_fire::i, phi, DiDispStaMuonMonitor_cfi::pt, and HLT_2018_cff::track.
Referenced by analyze().
◆ fillUE_with_ChargedJets()
Definition at line 1082 of file QcdUeDQM.cc.
1085 double nTrk500_TransReg = 0;
1086 double nTrk500_AwayReg = 0;
1087 double nTrk500_TowardReg = 0;
1089 double pTSum500_TransReg = 0;
1090 double pTSum500_AwayReg = 0;
1091 double pTSum500_TowardReg = 0;
1093 double nTrk900_TransReg = 0;
1094 double nTrk900_AwayReg = 0;
1095 double nTrk900_TowardReg = 0;
1097 double pTSum900_TransReg = 0;
1098 double pTSum900_AwayReg = 0;
1099 double pTSum900_TowardReg = 0;
1101 if (!(trackJets->empty()) && (trackJets->begin())->
pt() > 1.) {
1102 double jetPhi = (trackJets->begin())->
phi();
1103 for (
size_t i = 0;
i <
track.size();
i++) {
1106 if (fabs(dphi) > 60. && fabs(dphi) < 120.) {
1107 pTSum500_TransReg = pTSum500_TransReg +
track[
i]->pt();
1110 pTSum900_TransReg = pTSum900_TransReg +
track[
i]->pt();
1115 if (fabs(dphi) > 120. && fabs(dphi) < 180.) {
1116 pTSum500_AwayReg = pTSum500_AwayReg +
track[
i]->pt();
1119 pTSum900_AwayReg = pTSum900_AwayReg +
track[
i]->pt();
1123 if (fabs(dphi) < 60.) {
1124 pTSum500_TowardReg = pTSum500_TowardReg +
track[
i]->pt();
1125 nTrk500_TowardReg++;
1127 pTSum900_TowardReg = pTSum900_TowardReg +
track[
i]->pt();
1128 nTrk900_TowardReg++;
References SiPixelRawToDigiRegional_cfi::deltaPhi, fill1D(), fillProfile(), hdNdEtadPhi_trackJet_Away500_, hdNdEtadPhi_trackJet_Away900_, hdNdEtadPhi_trackJet_Toward500_, hdNdEtadPhi_trackJet_Toward900_, hdNdEtadPhi_trackJet_Transverse500_, hdNdEtadPhi_trackJet_Transverse900_, hdPhi_chargedJet_tracks_, hpTSumdEtadPhi_trackJet_Away500_, hpTSumdEtadPhi_trackJet_Away900_, hpTSumdEtadPhi_trackJet_Toward500_, hpTSumdEtadPhi_trackJet_Toward900_, hpTSumdEtadPhi_trackJet_Transverse500_, hpTSumdEtadPhi_trackJet_Transverse900_, mps_fire::i, reco::btau::jetPhi, phi, PI, DiDispStaMuonMonitor_cfi::pt, HLT_2018_cff::track, and w.
Referenced by analyze().
◆ fillUE_with_MaxpTtrack()
void QcdUeDQM::fillUE_with_MaxpTtrack |
( |
const std::vector< const reco::Track * > & |
track | ) |
|
|
private |
Definition at line 1012 of file QcdUeDQM.cc.
1015 double nTrk500_TransReg = 0;
1016 double nTrk500_AwayReg = 0;
1017 double nTrk500_TowardReg = 0;
1019 double pTSum500_TransReg = 0;
1020 double pTSum500_AwayReg = 0;
1021 double pTSum500_TowardReg = 0;
1023 double nTrk900_TransReg = 0;
1024 double nTrk900_AwayReg = 0;
1025 double nTrk900_TowardReg = 0;
1027 double pTSum900_TransReg = 0;
1028 double pTSum900_AwayReg = 0;
1029 double pTSum900_TowardReg = 0;
1030 if (!
track.empty()) {
1032 for (
size_t i = 1;
i <
track.size();
i++) {
1035 if (fabs(dphi) > 60. && fabs(dphi) < 120.) {
1036 pTSum500_TransReg = pTSum500_TransReg +
track[
i]->pt();
1039 pTSum900_TransReg = pTSum900_TransReg +
track[
i]->pt();
1044 if (fabs(dphi) > 120. && fabs(dphi) < 180.) {
1045 pTSum500_AwayReg = pTSum500_AwayReg +
track[
i]->pt();
1048 pTSum900_AwayReg = pTSum900_AwayReg +
track[
i]->pt();
1053 if (fabs(dphi) < 60.) {
1054 pTSum500_TowardReg = pTSum500_TowardReg +
track[
i]->pt();
1055 nTrk500_TowardReg++;
1057 pTSum900_TowardReg = pTSum900_TowardReg +
track[
i]->pt();
1058 nTrk900_TowardReg++;
References SiPixelRawToDigiRegional_cfi::deltaPhi, fill1D(), fillProfile(), hdNdEtadPhi_pTMax_Away500_, hdNdEtadPhi_pTMax_Away900_, hdNdEtadPhi_pTMax_Toward500_, hdNdEtadPhi_pTMax_Toward900_, hdNdEtadPhi_pTMax_Transverse500_, hdNdEtadPhi_pTMax_Transverse900_, hdPhi_maxpTTrack_tracks_, hpTSumdEtadPhi_pTMax_Away500_, hpTSumdEtadPhi_pTMax_Away900_, hpTSumdEtadPhi_pTMax_Toward500_, hpTSumdEtadPhi_pTMax_Toward900_, hpTSumdEtadPhi_pTMax_Transverse500_, hpTSumdEtadPhi_pTMax_Transverse900_, mps_fire::i, phi, PI, DiDispStaMuonMonitor_cfi::pt, HLT_2018_cff::track, and w.
Referenced by analyze().
◆ fillVtxPlots()
Definition at line 934 of file QcdUeDQM.cc.
936 bool goodVtx =
false;
938 for (reco::VertexCollection::const_iterator vertexIt = theVertices.begin(); vertexIt != theVertices.end();
947 if (fabs(vertexIt->z() -
bs->z0()) <
diffvtxbs_ && vertexIt->ndof() >= 4 && vertexIt->position().rho() <= 2.0) {
References cms::cuda::bs, diffvtxbs_, fill1D(), hNvertices_, hVertex_ndof_, hVertex_rho_, hVertex_x_, hVertex_y_, hVertex_z_, hVertex_z_bs_, edm::Handle< T >::product(), and vtx1.
Referenced by analyze().
◆ getProduct()
template<typename TYPE >
void QcdUeDQM::getProduct |
( |
const std::string |
name, |
|
|
edm::Handle< TYPE > & |
prod, |
|
|
const edm::Event & |
event |
|
) |
| const |
|
inlineprivate |
◆ getProductSafe()
template<typename TYPE >
bool QcdUeDQM::getProductSafe |
( |
const std::string |
name, |
|
|
edm::Handle< TYPE > & |
prod, |
|
|
const edm::Event & |
event |
|
) |
| const |
|
inlineprivate |
◆ setLabel1D()
◆ trackSelection()
Definition at line 832 of file QcdUeDQM.cc.
833 bool goodTrk =
false;
870 bool hasPIX1 =
false;
891 bool quality_ok =
true;
894 for (
unsigned int i = 0;
i <
quality_.size(); ++
i) {
909 if (hasPIX1 && pxlLayerMinCutbool && layerMinCutbool &&
920 if (hasPIX1 && pxlLayerMinCutbool && layerMinCutbool &&
References reco::TrackBase::algo(), algorithm_, allowTriplets_, cms::cuda::bs, bsuse_, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::TrackBase::eta(), fill1D(), spr::find(), reco::HitPattern::getLayer(), hBeamSpot_x_, hBeamSpot_y_, hBeamSpot_z_, reco::TrackBase::hitPattern(), hNtrackerLayer_, hNtrackerPixelLayer_, hNtrackerStripPixelLayer_, hRatioDxySigmaDxyBS_, hRatioDxySigmaDxyPV_, hRatioDzSigmaDzBS_, hRatioDzSigmaDzPV_, hRatioPtErrorPt_, hTrkChi2_, hTrkEta_, hTrkNdof_, hTrkPhi_, hTrkPt_, mps_fire::i, lip_, maxChi2_, maxRapidity_, min3DHit_, minHit_, minRapidity_, reco::TrackBase::ndof(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidStripLayersWithMonoAndStereo(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::phi(), reco::HitPattern::pixelHitFilter(), reco::HitPattern::pixelLayersWithMeasurement(), reco::TrackBase::pt(), ptErr_pt_, reco::TrackBase::ptError(), ptMin_, pxlLayerMinCut_, reco::TrackBase::quality(), quality_, requirePIX1_, tip_, reco::HitPattern::TRACK_HITS, reco::HitPattern::trackerLayersWithMeasurement(), reco::HitPattern::validHitFilter(), and badGlobalMuonTaggersAOD_cff::vtx.
Referenced by analyze().
◆ algorithm_
◆ allowTriplets_
bool QcdUeDQM::allowTriplets_ |
|
private |
◆ bsLabel_
◆ bsPos_
◆ bsuse_
◆ caloJetLabel_
◆ chargedJetLabel_
◆ diffvtxbs_
double QcdUeDQM::diffvtxbs_ |
|
private |
◆ h2TrigCorr_
◆ hBeamSpot_x_
◆ hBeamSpot_y_
◆ hBeamSpot_z_
◆ hChargedJetConstituent_
◆ hChargedJetMulti_
◆ hdNdEtadPhi_pTMax_Away500_
◆ hdNdEtadPhi_pTMax_Away900_
◆ hdNdEtadPhi_pTMax_Toward500_
◆ hdNdEtadPhi_pTMax_Toward900_
◆ hdNdEtadPhi_pTMax_Transverse500_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Transverse500_ |
|
private |
◆ hdNdEtadPhi_pTMax_Transverse900_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_pTMax_Transverse900_ |
|
private |
◆ hdNdEtadPhi_trackJet_Away500_
◆ hdNdEtadPhi_trackJet_Away900_
◆ hdNdEtadPhi_trackJet_Toward500_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Toward500_ |
|
private |
◆ hdNdEtadPhi_trackJet_Toward900_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Toward900_ |
|
private |
◆ hdNdEtadPhi_trackJet_Transverse500_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Transverse500_ |
|
private |
◆ hdNdEtadPhi_trackJet_Transverse900_
std::vector<MonitorElement *> QcdUeDQM::hdNdEtadPhi_trackJet_Transverse900_ |
|
private |
◆ hdPhi_chargedJet_tracks_
◆ hdPhi_maxpTTrack_tracks_
◆ hGoodTrkEta500_
◆ hGoodTrkEta900_
◆ hGoodTrkPhi500_
◆ hGoodTrkPhi900_
◆ hGoodTrkPt500_
◆ hGoodTrkPt900_
◆ hLeadingChargedJet_etaSpectrum_
std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_etaSpectrum_ |
|
private |
◆ hLeadingChargedJet_phiSpectrum_
std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_phiSpectrum_ |
|
private |
◆ hLeadingChargedJet_pTSpectrum_
std::vector<MonitorElement *> QcdUeDQM::hLeadingChargedJet_pTSpectrum_ |
|
private |
◆ hLeadingTrack_etaSpectrum_
◆ hLeadingTrack_phiSpectrum_
◆ hLeadingTrack_pTSpectrum_
◆ hltConfig
◆ hltProcNames_
std::vector<std::string> QcdUeDQM::hltProcNames_ |
|
private |
◆ hltResName_
std::string QcdUeDQM::hltResName_ |
|
private |
◆ hltTrgBits_
std::vector<int> QcdUeDQM::hltTrgBits_ |
|
private |
◆ hltTrgDeci_
std::vector<bool> QcdUeDQM::hltTrgDeci_ |
|
private |
◆ hltTrgNames_
std::vector<std::string> QcdUeDQM::hltTrgNames_ |
|
private |
◆ hltTrgUsedNames_
std::vector<std::string> QcdUeDQM::hltTrgUsedNames_ |
|
private |
◆ hltUsedResName_
std::string QcdUeDQM::hltUsedResName_ |
|
private |
◆ hNevts_
◆ hNgoodTrk_
◆ hNtrackerLayer_
◆ hNtrackerPixelLayer_
◆ hNtrackerStripPixelLayer_
◆ hNvertices_
◆ hpTSumdEtadPhi_pTMax_Away500_
◆ hpTSumdEtadPhi_pTMax_Away900_
◆ hpTSumdEtadPhi_pTMax_Toward500_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Toward500_ |
|
private |
◆ hpTSumdEtadPhi_pTMax_Toward900_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Toward900_ |
|
private |
◆ hpTSumdEtadPhi_pTMax_Transverse500_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Transverse500_ |
|
private |
◆ hpTSumdEtadPhi_pTMax_Transverse900_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_pTMax_Transverse900_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Away500_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Away500_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Away900_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Away900_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Toward500_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Toward500_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Toward900_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Toward900_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Transverse500_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Transverse500_ |
|
private |
◆ hpTSumdEtadPhi_trackJet_Transverse900_
std::vector<MonitorElement *> QcdUeDQM::hpTSumdEtadPhi_trackJet_Transverse900_ |
|
private |
◆ hRatioDxySigmaDxyBS_
◆ hRatioDxySigmaDxyPV_
◆ hRatioDzSigmaDzBS_
◆ hRatioDzSigmaDzPV_
◆ hRatioPtErrorPt_
◆ hTrkChi2_
◆ hTrkEta_
◆ hTrkNdof_
◆ hTrkPhi_
◆ hTrkPt_
◆ hVertex_ndof_
◆ hVertex_rho_
◆ hVertex_x_
◆ hVertex_y_
◆ hVertex_z_
◆ hVertex_z_bs_
◆ isHltConfigSuccessful_
bool QcdUeDQM::isHltConfigSuccessful_ |
|
private |
◆ lip_
◆ maxChi2_
double QcdUeDQM::maxChi2_ |
|
private |
◆ maxRapidity_
double QcdUeDQM::maxRapidity_ |
|
private |
◆ min3DHit_
◆ minHit_
◆ minRapidity_
double QcdUeDQM::minRapidity_ |
|
private |
◆ ptErr_pt_
double QcdUeDQM::ptErr_pt_ |
|
private |
◆ ptMin_
◆ pxlLayerMinCut_
double QcdUeDQM::pxlLayerMinCut_ |
|
private |
◆ quality_
◆ repSumMap_
◆ repSummary_
◆ requirePIX1_
bool QcdUeDQM::requirePIX1_ |
|
private |
◆ selected_
◆ tgeo_
◆ tip_
◆ trackLabel_
◆ verbose_
◆ vtx1
◆ vtxLabel_
◆ vtxntk_
edm::EDGetTokenT< reco::TrackJetCollection > chargedJetLabel_
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Away900_
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Toward500_
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Transverse900_
std::vector< std::string > hltProcNames_
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
const TrackerGeometry * tgeo_
std::vector< MonitorElement * > hVertex_ndof_
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Away500_
T const * product() const
edm::EDGetTokenT< reco::VertexCollection > vtxLabel_
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Toward500_
std::vector< MonitorElement * > hTrkPhi_
void bookProfile(DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, double y1, double y2, bool sumw2=true, bool sbox=true)
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Transverse900_
static bool validHitFilter(uint16_t pattern)
std::vector< Vertex > VertexCollection
collection of Vertex objects
void fillUE_with_ChargedJets(const std::vector< const reco::Track * > &track, const edm::Handle< reco::TrackJetCollection > &trackJets)
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Transverse500_
std::vector< MonitorElement * > hdPhi_chargedJet_tracks_
std::vector< MonitorElement * > hNtrackerStripPixelLayer_
bool trackSelection(const reco::Track &trk, const reco::BeamSpot *bs, const reco::Vertex &vtx, int sizevtx)
MonitorElement * h2TrigCorr_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Toward900_
std::vector< reco::TrackBase::TrackQuality > quality_
double ndof() const
number of degrees of freedom of the fit
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Toward900_
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
std::vector< MonitorElement * > hGoodTrkPt500_
void fillpTMaxRelated(const std::vector< const reco::Track * > &track)
double dxyError() const
error on dxy
std::vector< MonitorElement * > hTrkPt_
void fillChargedJetSpectra(const edm::Handle< reco::TrackJetCollection > trackJets)
int pixelLayersWithMeasurement() const
std::vector< MonitorElement * > hRatioDxySigmaDxyBS_
double pt() const
track transverse momentum
std::vector< MonitorElement * > hRatioDxySigmaDxyPV_
std::vector< MonitorElement * > hNgoodTrk_
std::vector< MonitorElement * > hBeamSpot_z_
std::vector< MonitorElement * > hRatioPtErrorPt_
void fillUE_with_MaxpTtrack(const std::vector< const reco::Track * > &track)
std::vector< MonitorElement * > hVertex_x_
std::vector< std::string > hltTrgUsedNames_
int trackerLayersWithMeasurement() const
std::vector< MonitorElement * > hLeadingChargedJet_pTSpectrum_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Away500_
HLTConfigProvider hltConfig
std::vector< bool > hltTrgDeci_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Away900_
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
edm::EDGetTokenT< reco::CaloJetCollection > caloJetLabel_
double phi() const
azimuthal angle of momentum vector
std::vector< MonitorElement * > hLeadingChargedJet_etaSpectrum_
unsigned int size() const
number of trigger paths in trigger table
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Away500_
static bool pixelHitFilter(uint16_t pattern)
static uint32_t getLayer(uint16_t pattern)
std::vector< MonitorElement * > hBeamSpot_y_
void book1D(DQMStore::IBooker &, std::vector< MonitorElement * > &mes, const std::string &name, const std::string &title, int nx, double x1, double x2, bool sumw2=true, bool sbox=true)
std::vector< MonitorElement * > hVertex_rho_
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Toward500_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Toward500_
void fillHltBits(const edm::Event &iEvent, const edm::EventSetup &iSetup)
double eta() const
pseudorapidity of momentum vector
void fill1D(std::vector< TH1F * > &hs, double val, double w=1.)
std::vector< MonitorElement * > hRatioDzSigmaDzPV_
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Away900_
MonitorElement * repSummary_
double dzError() const
error on dz
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Away900_
std::vector< MonitorElement * > hLeadingTrack_etaSpectrum_
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
edm::EDGetTokenT< reco::TrackCollection > trackLabel_
std::vector< MonitorElement * > hNtrackerPixelLayer_
std::vector< MonitorElement * > hGoodTrkPt900_
std::vector< MonitorElement * > hVertex_z_
std::vector< int > hltTrgBits_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Transverse900_
bool isHltConfigSuccessful_
std::vector< MonitorElement * > hTrkNdof_
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Transverse500_
std::vector< MonitorElement * > hLeadingTrack_pTSpectrum_
std::vector< MonitorElement * > hLeadingChargedJet_phiSpectrum_
std::vector< MonitorElement * > hGoodTrkEta500_
void fillProfile(std::vector< TProfile * > &hs, double valx, double valy, double w=1.)
std::vector< MonitorElement * > hLeadingTrack_phiSpectrum_
std::vector< MonitorElement * > hTrkEta_
TrackAlgorithm algo() const
static TrackQuality qualityByName(const std::string &name)
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
std::vector< MonitorElement * > hGoodTrkPhi900_
std::string hltUsedResName_
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
std::vector< MonitorElement * > hBeamSpot_x_
std::vector< MonitorElement * > hGoodTrkEta900_
std::vector< MonitorElement * > hdNdEtadPhi_trackJet_Transverse500_
std::vector< MonitorElement * > hdPhi_maxpTTrack_tracks_
bool fillVtxPlots(const reco::BeamSpot *bs, const edm::Handle< reco::VertexCollection > vtxColl)
std::vector< MonitorElement * > hGoodTrkPhi500_
std::vector< MonitorElement * > hRatioDzSigmaDzBS_
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Away500_
std::vector< MonitorElement * > hChargedJetMulti_
bool accept() const
Has at least one path accepted the event?
std::vector< reco::TrackBase::TrackAlgorithm > algorithm_
const std::string & triggerName(unsigned int triggerIndex) const
std::vector< MonitorElement * > hTrkChi2_
void getProduct(const std::string name, edm::Handle< TYPE > &prod, const edm::Event &event) const
static TrackAlgorithm algoByName(const std::string &name)
std::vector< MonitorElement * > hdNdEtadPhi_pTMax_Transverse900_
std::vector< MonitorElement * > hVertex_y_
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Transverse500_
std::vector< MonitorElement * > hNvertices_
std::vector< MonitorElement * > hpTSumdEtadPhi_trackJet_Toward900_
std::vector< MonitorElement * > hVertex_z_bs_
MonitorElement * repSumMap_
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
std::vector< MonitorElement * > hNevts_
bool quality(const TrackQuality) const
Track quality.
std::vector< Track > TrackCollection
collection of Tracks
std::vector< MonitorElement * > hNtrackerLayer_
int numberOfValidStripLayersWithMonoAndStereo(uint16_t stripdet, uint16_t layer) const
std::vector< MonitorElement * > hpTSumdEtadPhi_pTMax_Toward900_
edm::EDGetTokenT< reco::BeamSpot > bsLabel_
std::vector< std::string > hltTrgNames_