4 : monitorDir(ps.getUntrackedParameter<
std::
string>(
"monitorDir",
"")),
6 ps.getParameter<
edm::InputTag>(
7 "calol2JetCollectionData"))),
9 ps.getParameter<
edm::InputTag>(
10 "calol2JetCollectionEmul"))),
12 ps.getParameter<
edm::InputTag>(
13 "calol2EGammaCollectionData"))),
15 ps.getParameter<
edm::InputTag>(
16 "calol2EGammaCollectionEmul"))),
18 ps.getParameter<
edm::InputTag>(
19 "calol2TauCollectionData"))),
21 ps.getParameter<
edm::InputTag>(
22 "calol2TauCollectionEmul"))),
24 ps.getParameter<
edm::InputTag>(
25 "calol2EtSumCollectionData"))),
27 ps.getParameter<
edm::InputTag>(
28 "calol2EtSumCollectionEmul"))),
30 enable2DComp(ps.getUntrackedParameter<
bool> (
"enable2DComp",
false))
56 50,-5.,5.,25,-3.2,3.2);
61 50,-5.,5.,25,-3.2,3.2);
84 30,-3.,3.,25,-3.2,3.2);
89 30,-3.,3.,25,-3.2,3.2);
94 "Iso Eg iE_{T}", 1400, 0, 1399);
96 "Iso Eg i#eta", 227, -113.5, 113.5);
98 "Iso Eg i#phi", 288, -0.5, 143.5);
100 "Iso Eg iE_{T}", 1400, 0, 1399);
102 "Iso Eg i#eta", 227, -113.5, 113.5);
104 "Iso Eg i#phi", 288, -0.5, 143.5);
108 30,-3.,3.,25,-3.2,3.2);
113 30,-3.,3.,25,-3.2,3.2);
135 30,-3.,3.,25,-3.2,3.2);
140 30,-3.,3.,25,-3.2,3.2);
145 "Iso Tau iE_{T}", 1400, 0, 1399);
147 "Iso Tau i#eta", 227, -113.5, 113.5);
149 "Iso Tau i#phi", 288, -0.5, 143.5);
151 "Iso Tau iE_{T}", 1400, 0, 1399);
153 "Iso Tau i#eta", 227, -113.5, 113.5);
155 "Iso Tau i#phi", 288, -0.5, 143.5);
159 30,-3.,3.,25,-3.2,3.2);
164 30,-3.,3.,25,-3.2,3.2);
172 ettData = ibooker.
book1D(
"Problematic ETT Sum - Data",
"ETT iE_{T}",
174 ettEmul = ibooker.
book1D(
"Problematic ETT Sum - Emulator",
"ETT iE_{T}",
176 ettHFData = ibooker.
book1D(
"Problematic ETTHF Sum - Data",
"ETTHF iE_{T}",
178 ettHFEmul = ibooker.
book1D(
"Problematic ETTHF Sum - Emulator",
"ETTHF iE_{T}",
180 ettEmData = ibooker.
book1D(
"Problematic ETTEM Sum - Data",
"ETTEM iE_{T}",
182 ettEmEmul = ibooker.
book1D(
"Problematic ETTEM Sum - Emulator",
"ETTEM iE_{T}",
186 httData = ibooker.
book1D(
"Problematic HTT Sum - Data",
"HTT iE_{T}",
188 httEmul = ibooker.
book1D(
"Problematic HTT Sum - Emulator",
"HTT iE_{T}",
190 httHFData = ibooker.
book1D(
"Problematic HTTHF Sum - Data",
"HTTHF iE_{T}",
192 httHFEmul = ibooker.
book1D(
"Problematic HTTHF Sum - Emulator",
"HTTHF iE_{T}",
196 metEtData = ibooker.
book1D(
"Problematic MET Sum Et - Data",
"MET iE_{T}",
198 metEtEmul = ibooker.
book1D(
"Problematic MET Sum Et - Emulator",
"MET iE_{T}",
202 metPhiEmul = ibooker.
book1D(
"Problematic MET Sum phi - Emulator",
"MET i#phi",
206 "METHF iE_{T}", 7000, -0.5, 6999.5);
208 "METHF iE_{T}", 7000, -0.5, 6999.5);
210 "METHF i#phi", 1008, -0.5, 1007.5);
212 "METHF i#phi", 1008, -0.5, 1007.5);
215 mhtEtData = ibooker.
book1D(
"Problematic MHT Sum Et - Data",
"MHT iE_{T}",
217 mhtEtEmul = ibooker.
book1D(
"Problematic MHT Sum Et - Emulator",
"MHT iE_{T}",
221 mhtPhiEmul = ibooker.
book1D(
"Problematic MHT Sum phi - Emulator",
"MHT i#phi",
225 "MHTHF iE_{T}", 7000, -0.5, 6999.5);
227 "MHTHF iE_{T}", 7000, -0.5, 6999.5);
229 "MHTHF i#phi", 1008, -0.5, 1007.5);
231 "MHTHF i#phi", 1008, -0.5, 1007.5);
253 "", 5904, -0.5, 5903.5);
255 "", 5904, -0.5, 5903.5);
262 "", 256, -0.5, 255.5);
264 "", 256, -0.5, 255.5);
288 "CaloL2 Object Agreement Summary",
289 "CaloL2 event-by-event object agreement fractions", 10, 1, 11);
303 "Jet Agreement Summary",
"Jet Agreement Summary", 4, 1, 5);
310 "EG Agreement Summary",
"EG Agreement Summary", 8, 1, 9);
321 "Tau Agreement Summary",
"Tau Agreement Summary", 8, 1, 9);
332 "Energy Sum Agreement Summary",
"Sum Agreement Summary", 18, 1, 19);
357 "Problem Summary",
"Problematic Event Summary", 8, 1, 9);
372 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"L1TdeStage2CaloLayer2: analyse " 395 bool eventGood =
true;
421 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"l1t calol2 dqm: jet problem " 429 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"l1t calol2 dqm: eg problem " 437 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"l1t calol2 dqm: tau problem " 445 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"l1t calol2 dqm: sum problem " 470 bool eventGood =
true;
529 if (dataIt->hwPt() != emulIt->hwPt()) {
535 if (dataIt->hwPhi() != emulIt->hwPhi()){
542 if (dataIt->hwEta() != emulIt->hwEta()) {
548 if (etGood && posGood) {
563 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"--- jet ---"<< std::endl;
564 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data jet Et = " 565 << dataIt->hwPt() << std::endl;
566 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul jet Et = " 567 << emulIt->hwPt() << std::endl;
568 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data jet phi = " 569 << dataIt->hwPhi() << std::endl;
570 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul jet phi = " 571 << emulIt->hwPhi() << std::endl;
572 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data jet eta = " 573 << dataIt->hwEta() << std::endl;
574 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul jet eta = " 575 << emulIt->hwEta() << std::endl;
576 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"---"<< std::endl;
581 if (posGood && !etGood) {
586 if (!posGood && etGood) {
617 bool eventGood =
true;
632 if (dataIt->hwIso()) {
655 if(emulIt->hwIso()) {
686 bool iso = dataIt->hwIso();
689 if (dataIt->hwPt() != emulIt->hwPt()) {
695 if (dataIt->hwPhi() != emulIt->hwPhi()) {
701 if (dataIt->hwEta() != emulIt->hwEta()) {
707 if (posGood && etGood) {
741 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"--- eg ---"<< std::endl;
742 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data eg Et = " 743 << dataIt->hwPt() << std::endl;
744 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul eg Et = " 745 << emulIt->hwPt() << std::endl;
746 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data eg phi = " 747 << dataIt->hwPhi() << std::endl;
748 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul eg phi = " 749 << emulIt->hwPhi() << std::endl;
750 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data eg eta = " 751 << dataIt->hwEta() << std::endl;
752 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul eg eta = " 753 << emulIt->hwEta() << std::endl;
754 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"---"<< std::endl;
759 if (posGood && !etGood) {
768 if (!posGood && etGood) {
807 bool eventGood =
true;
825 if (dataIt->hwIso()) {
851 if(emulIt->hwIso()) {
883 bool iso = dataIt->hwIso();
886 if (dataIt->hwPt() != emulIt->hwPt()) {
892 if (dataIt->hwPhi() != emulIt->hwPhi()) {
898 if (dataIt->hwEta() != emulIt->hwEta()) {
904 if (posGood && etGood) {
938 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"--- tau ---"<< std::endl;
939 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data tau Et = " 940 << dataIt->hwPt() << std::endl;
941 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul tau Et = " 942 << emulIt->hwPt() << std::endl;
943 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data tau phi = " 944 << dataIt->hwPhi() << std::endl;
945 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul tau phi = " 946 << emulIt->hwPhi() << std::endl;
947 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"data tau eta = " 948 << dataIt->hwEta() << std::endl;
949 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"emul tau eta = " 950 << emulIt->hwEta() << std::endl;
951 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"---"<< std::endl;
956 if (posGood && !etGood) {
965 if (!posGood && etGood) {
1005 bool eventGood =
true;
1008 bool phiGood =
true;
1028 if (dataIt->getType()==emulIt->getType()){
1035 dataEt = dataIt->hwPt();
1036 emulEt = emulIt->hwPt();
1038 if (dataEt != emulEt) {
1051 if (dataEt != emulEt)
1055 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1068 dataEt = dataIt->hwPt();
1069 emulEt = emulIt->hwPt();
1071 if (dataEt != emulEt) {
1084 if (dataEt != emulEt)
1088 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1101 dataEt = dataIt->hwPt();
1102 emulEt = emulIt->hwPt();
1104 if (dataEt != emulEt) {
1116 if (dataEt != emulEt)
1120 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1132 && dataIt->hwPt() != 0) {
1134 dataEt = dataIt->hwPt();
1135 emulEt = emulIt->hwPt();
1137 dataPhi = dataIt->hwPhi();
1138 emulPhi = emulIt->hwPhi();
1140 if (dataEt != emulEt) {
1145 if (dataPhi != emulPhi) {
1150 if (etGood && phiGood) {
1163 if (dataEt != emulEt)
1167 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1171 if (dataPhi != emulPhi)
1175 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataPhi <<
"\t" << emulPhi;
1187 && dataIt->hwPt() != 0) {
1189 dataEt = dataIt->hwPt();
1190 emulEt = emulIt->hwPt();
1192 dataPhi = dataIt->hwPhi();
1193 emulPhi = emulIt->hwPhi();
1195 if (dataEt != emulEt) {
1200 if (dataPhi != emulPhi) {
1205 if (etGood && phiGood) {
1218 if (dataEt != emulEt)
1222 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1225 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"METHF phi | ";
1226 if (dataPhi != emulPhi)
1230 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataPhi <<
"\t" << emulPhi;
1242 && dataIt->hwPt() != 0) {
1244 dataEt = dataIt->hwPt();
1245 emulEt = emulIt->hwPt();
1247 dataPhi = dataIt->hwPhi();
1248 emulPhi = emulIt->hwPhi();
1250 if (dataEt != emulEt) {
1255 if (!(etGood && dataEt == 0)) {
1256 if (dataPhi != emulPhi) {
1262 if (etGood && phiGood) {
1275 if (dataEt != emulEt)
1279 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1283 if (dataPhi != emulPhi)
1287 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataPhi <<
"\t" << emulPhi;
1299 && dataIt->hwPt() != 0) {
1301 dataEt = dataIt->hwPt();
1302 emulEt = emulIt->hwPt();
1304 dataPhi = dataIt->hwPhi();
1305 emulPhi = emulIt->hwPhi();
1307 if (dataEt != emulEt) {
1312 if (!(etGood && dataEt == 0)) {
1313 if (dataPhi != emulPhi) {
1319 if (etGood && phiGood) {
1332 if (dataEt != emulEt)
1336 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1339 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"MHTHF phi | ";
1340 if (dataPhi != emulPhi)
1344 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataPhi <<
"\t" << emulPhi;
1357 dataEt = dataIt->hwPt();
1358 emulEt = emulIt->hwPt();
1360 if (dataEt != emulEt) {
1372 if (dataEt != emulEt)
1376 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1389 dataEt = dataIt->hwPt();
1390 emulEt = emulIt->hwPt();
1392 if (dataEt != emulEt) {
1411 dataEt = dataIt->hwPt();
1412 emulEt = emulIt->hwPt();
1414 if (dataEt != emulEt) {
1426 if (dataEt != emulEt)
1430 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1443 dataEt = dataIt->hwPt();
1444 emulEt = emulIt->hwPt();
1448 if (dataEt != emulEt) {
1460 if (dataEt != emulEt)
1464 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1475 dataEt = dataIt->hwPt();
1476 emulEt = emulIt->hwPt();
1478 if (dataEt != emulEt) {
1489 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"TowCount | ";
1490 if (dataEt != emulEt)
1494 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1507 dataEt = dataIt->hwPt();
1508 emulEt = emulIt->hwPt();
1510 if (dataEt != emulEt) {
1521 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"AsymCount | ";
1522 if (dataEt != emulEt)
1527 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
1538 dataCentr = dataIt->hwPt();
1539 emulCentr = emulIt->hwPt();
1541 if (dataCentr != emulCentr) {
1545 for (
int i=0;
i<8;
i++)
1551 for (
int i=0;
i<8;
i++)
1562 edm::LogInfo(
"L1TdeStage2CaloLayer2") <<
"CentrCount | ";
1563 if (dataEt != emulEt)
1568 edm::LogInfo(
"L1TdeStage2CaloLayer2") << dataEt <<
"\t" << emulEt;
BXVector< EGamma > EGammaBxCollection
MonitorElement * ettHFEmul
edm::EDGetTokenT< l1t::JetBxCollection > calol2JetCollectionEmul
const_iterator end(int bx) const
MonitorElement * mbhfp1Data
edm::EDGetTokenT< l1t::TauBxCollection > calol2TauCollectionEmul
MonitorElement * jet2DEtaPhiData
MonitorElement * eg2DEtaPhiEmul
MonitorElement * metHFPhiEmul
MonitorElement * problemSummary
unsigned size(int bx) const
MonitorElement * centrCountData
MonitorElement * egEtData
MonitorElement * isoTauPhiEmul
MonitorElement * tauEtaEmul
L1TdeStage2CaloLayer2(const edm::ParameterSet &ps)
MonitorElement * httHFEmul
MonitorElement * towCountEmul
MonitorElement * metEtData
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * mhtEtEmul
MonitorElement * egPhiEmul
MonitorElement * metHFEtEmul
MonitorElement * mhtHFPhiData
MonitorElement * mhtEtData
MonitorElement * centrCountEmul
edm::EDGetTokenT< l1t::EtSumBxCollection > calol2EtSumCollectionData
MonitorElement * egSummary
MonitorElement * isoTauEtaData
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)
MonitorElement * agreementSummary
MonitorElement * jet2DEtaPhiEmul
bool compareSums(const edm::Handle< l1t::EtSumBxCollection > &dataCol, const edm::Handle< l1t::EtSumBxCollection > &emulCol)
MonitorElement * jetPhiEmul
bool isEmpty(int bx) const
MonitorElement * asymCountData
MonitorElement * mbhfm0Emul
BXVector< Tau > TauBxCollection
MonitorElement * egEtaData
MonitorElement * jetEtaData
MonitorElement * tauPhiEmul
MonitorElement * jetEtData
MonitorElement * isoEgEtData
edm::EDGetTokenT< l1t::EGammaBxCollection > calol2EGammaCollectionEmul
MonitorElement * isoEgPhiData
void bookHistograms(DQMStore::IBooker &, const edm::Run &, const edm::EventSetup &) override
MonitorElement * sumSummary
edm::EDGetTokenT< l1t::EGammaBxCollection > calol2EGammaCollectionData
BXVector< EtSum > EtSumBxCollection
MonitorElement * isoTauEtData
MonitorElement * isoEgEtaData
MonitorElement * tau2DEtaPhiEmul
MonitorElement * jetEtEmul
void setCurrentFolder(std::string const &fullpath)
MonitorElement * mbhfp0Emul
MonitorElement * metHFEtData
MonitorElement * tauSummary
MonitorElement * mhtHFPhiEmul
MonitorElement * isoTauEtEmul
MonitorElement * isoEgEtaEmul
bool compareEGs(const edm::Handle< l1t::EGammaBxCollection > &dataCol, const edm::Handle< l1t::EGammaBxCollection > &emulCol)
MonitorElement * book1D(Args &&...args)
MonitorElement * jetSummary
MonitorElement * egEtaEmul
MonitorElement * ettEmEmul
MonitorElement * mbhfm1Data
MonitorElement * isoTauEtaEmul
MonitorElement * jetPhiData
MonitorElement * mbhfm1Emul
MonitorElement * asymCountEmul
MonitorElement * mbhfp0Data
MonitorElement * mbhfp1Emul
BXVector< Jet > JetBxCollection
void analyze(const edm::Event &, const edm::EventSetup &) override
bool compareTaus(const edm::Handle< l1t::TauBxCollection > &dataCol, const edm::Handle< l1t::TauBxCollection > &emulCol)
MonitorElement * ettHFData
MonitorElement * book2D(Args &&...args)
MonitorElement * isoTau2DEtaPhiData
MonitorElement * mbhfm0Data
edm::EDGetTokenT< l1t::JetBxCollection > calol2JetCollectionData
MonitorElement * egPhiData
MonitorElement * metPhiEmul
MonitorElement * egEtEmul
MonitorElement * jetEtaEmul
MonitorElement * httHFData
bool compareJets(const edm::Handle< l1t::JetBxCollection > &dataCol, const edm::Handle< l1t::JetBxCollection > &emulCol)
MonitorElement * isoTauPhiData
MonitorElement * isoTau2DEtaPhiEmul
MonitorElement * tauEtEmul
MonitorElement * ettEmData
MonitorElement * isoEg2DEtaPhiData
MonitorElement * towCountData
MonitorElement * mhtPhiEmul
edm::EDGetTokenT< l1t::TauBxCollection > calol2TauCollectionData
MonitorElement * tau2DEtaPhiData
MonitorElement * metEtEmul
MonitorElement * isoEgEtEmul
MonitorElement * tauEtData
MonitorElement * mhtHFEtData
MonitorElement * tauEtaData
const unsigned int currBx
const_iterator begin(int bx) const
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
edm::EDGetTokenT< l1t::EtSumBxCollection > calol2EtSumCollectionEmul
MonitorElement * isoEgPhiEmul
MonitorElement * isoEg2DEtaPhiEmul
MonitorElement * tauPhiData
MonitorElement * mhtHFEtEmul
MonitorElement * mhtPhiData
MonitorElement * eg2DEtaPhiData
std::vector< T >::const_iterator const_iterator
MonitorElement * metHFPhiData
MonitorElement * metPhiData