259 std::stringstream
text;
266 for (
int ibx = gmuons->getFirstBX(); ibx <= gmuons->getLastBX(); ++ibx) {
267 for (
auto itr = gmuons->begin(ibx); itr != gmuons->end(ibx); ++itr) {
269 het_.at(
gMU)->Fill(itr->hwPt());
270 heta_.at(
gMU)->Fill(itr->hwEtaAtVtx());
271 hphi_.at(
gMU)->Fill(itr->hwPhiAtVtx());
272 hetaphi_.at(
gMU)->Fill(itr->hwEtaAtVtx(), itr->hwPhiAtVtx(), itr->hwPt());
275 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEtaAtVtx()
276 <<
" iphi=" << itr->hwPhiAtVtx() << std::endl;
286 for (
int ibx = dmxegs->getFirstBX(); ibx <= dmxegs->getLastBX(); ++ibx) {
287 for (
auto itr = dmxegs->begin(ibx); itr != dmxegs->end(ibx); ++itr) {
292 hetaphi_.at(
DmxEG)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
295 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
306 for (
int ibx = dmxtaus->getFirstBX(); ibx <= dmxtaus->getLastBX(); ++ibx) {
307 for (
auto itr = dmxtaus->begin(ibx); itr != dmxtaus->end(ibx); ++itr) {
312 hetaphi_.at(
DmxTau)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
315 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
326 for (
int ibx = dmxjets->getFirstBX(); ibx <= dmxjets->getLastBX(); ++ibx) {
327 for (
auto itr = dmxjets->begin(ibx); itr != dmxjets->end(ibx); ++itr) {
332 hetaphi_.at(
DmxJet)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
335 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
346 for (
int ibx = dmxsums->getFirstBX(); ibx <= dmxsums->getLastBX(); ++ibx) {
347 for (
auto itr = dmxsums->begin(ibx); itr != dmxsums->end(ibx); ++itr) {
352 hetaphi_.at(
DmxSum)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
355 <<
" type=" << itr->getType() <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta()
356 <<
" iphi=" << itr->hwPhi() << std::endl;
366 for (
int ibx =
muons->getFirstBX(); ibx <=
muons->getLastBX(); ++ibx) {
367 for (
auto itr =
muons->begin(ibx); itr !=
muons->end(ibx); ++itr) {
369 het_.at(
MU)->Fill(itr->hwPt());
370 heta_.at(
MU)->Fill(itr->hwEtaAtVtx());
371 hphi_.at(
MU)->Fill(itr->hwPhiAtVtx());
372 hetaphi_.at(
MU)->Fill(itr->hwEtaAtVtx(), itr->hwPhiAtVtx(), itr->hwPt());
375 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEtaAtVtx()
376 <<
" iphi=" << itr->hwPhiAtVtx() << std::endl;
386 for (
int ibx = egs->getFirstBX(); ibx <= egs->getLastBX(); ++ibx) {
387 for (
auto itr = egs->begin(ibx); itr != egs->end(ibx); ++itr) {
389 het_.at(
EG)->Fill(itr->hwPt());
390 heta_.at(
EG)->Fill(itr->hwEta());
391 hphi_.at(
EG)->Fill(itr->hwPhi());
392 hetaphi_.at(
EG)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
395 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
406 for (
int ibx =
taus->getFirstBX(); ibx <=
taus->getLastBX(); ++ibx) {
407 for (
auto itr =
taus->begin(ibx); itr !=
taus->end(ibx); ++itr) {
409 het_.at(
Tau)->Fill(itr->hwPt());
412 hetaphi_.at(
Tau)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
415 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
426 for (
int ibx =
jets->getFirstBX(); ibx <=
jets->getLastBX(); ++ibx) {
427 for (
auto itr =
jets->begin(ibx); itr !=
jets->end(ibx); ++itr) {
429 het_.at(
Jet)->Fill(itr->hwPt());
432 hetaphi_.at(
Jet)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
435 <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta() <<
" iphi=" << itr->hwPhi()
446 for (
int ibx = sums->getFirstBX(); ibx <= sums->getLastBX(); ++ibx) {
447 for (
auto itr = sums->begin(ibx); itr != sums->end(ibx); ++itr) {
449 het_.at(
Sum)->Fill(itr->hwPt());
452 hetaphi_.at(
Sum)->Fill(itr->hwEta(), itr->hwPhi(), itr->hwPt());
454 <<
" type=" << itr->getType() <<
" BX=" << ibx <<
" ipt=" << itr->hwPt() <<
" ieta=" << itr->hwEta()
455 <<
" iphi=" << itr->hwPhi() << std::endl;
465 for (
int ibx = algs->getFirstBX(); ibx <= algs->getLastBX(); ++ibx) {
466 for (
auto itr = algs->begin(ibx); itr != algs->end(ibx); ++itr) {
467 for (
int algBit = 0; algBit < 128;
469 if (itr->getAlgoDecisionFinal(algBit)) {
471 text <<
"HW Fired Alg Bit : " << algBit << std::endl;
483 for (
int ibx = emulGtAlgs->getFirstBX(); ibx <= emulGtAlgs->getLastBX(); ++ibx) {
484 for (
auto itr = emulGtAlgs->begin(ibx); itr != emulGtAlgs->end(ibx); ++itr) {
485 for (
int algBit = 0; algBit < 128;
487 if (itr->getAlgoDecisionFinal(algBit)) {
489 text <<
"Emul w/ GTInput Fired Alg Bit : " << algBit << std::endl;
501 for (
int ibx = emulDxAlgs->getFirstBX(); ibx <= emulDxAlgs->getLastBX(); ++ibx) {
502 for (
auto itr = emulDxAlgs->begin(ibx); itr != emulDxAlgs->end(ibx); ++itr) {
503 for (
int algBit = 0; algBit < 128;
505 if (itr->getAlgoDecisionFinal(algBit)) {
507 text <<
"Emul w/ Demux output Fired Alg Bit : " << algBit << std::endl;
522 for (
int ibx = mus->getFirstBX(); ibx <= mus->getLastBX(); ++ibx) {
524 for (
unsigned int igtMU = 0; igtMU < mus->size(ibx); igtMU++) {
525 double gtMUEt = mus->at(ibx, igtMU).hwPt();
526 double gmtMUEt = 0.0;
527 if (gmtmus->size(ibx) > igtMU)
528 gmtMUEt = gmtmus->at(ibx, igtMU).hwPt();
531 double gtMUEta = mus->at(ibx, igtMU).hwEtaAtVtx();
532 double gmtMUEta = 0.0;
533 if (gmtmus->size(ibx) > igtMU)
534 gmtMUEta = gmtmus->at(ibx, igtMU).hwEtaAtVtx();
537 double gtMUPhi = mus->at(ibx, igtMU).hwPhiAtVtx();
538 double gmtMUPhi = 0.0;
539 if (gmtmus->size(ibx) > igtMU)
540 gmtMUPhi = gmtmus->at(ibx, igtMU).hwPhiAtVtx();
544 for (
unsigned int igmtMU = mus->size(ibx); igmtMU < gmtmus->size(ibx); igmtMU++) {
546 double gmtMUEt = gmtmus->at(ibx, igmtMU).hwPt();
549 double gtMUEta = 0.0;
550 double gmtMUEta = gmtmus->at(ibx, igmtMU).hwEtaAtVtx();
553 double gtMUPhi = 0.0;
554 double gmtMUPhi = gmtmus->at(ibx, igmtMU).hwPhiAtVtx();
568 for (
int ibx = egs->getFirstBX(); ibx <= egs->getLastBX(); ++ibx) {
570 for (
unsigned int igtEG = 0; igtEG < egs->size(ibx); igtEG++) {
571 double gtEGEt = egs->at(ibx, igtEG).hwPt();
572 double dmxEGEt = 0.0;
573 if (dmxegs->size(ibx) > igtEG)
574 dmxEGEt = dmxegs->at(ibx, igtEG).hwPt();
577 double gtEGEta = egs->at(ibx, igtEG).hwEta();
578 double dmxEGEta = 0.0;
579 if (dmxegs->size(ibx) > igtEG)
580 dmxEGEta = dmxegs->at(ibx, igtEG).hwEta();
583 double gtEGPhi = egs->at(ibx, igtEG).hwPhi();
584 double dmxEGPhi = 0.0;
585 if (dmxegs->size(ibx) > igtEG)
586 dmxEGPhi = dmxegs->at(ibx, igtEG).hwPhi();
590 for (
unsigned int idmEG = egs->size(ibx); idmEG < dmxegs->size(ibx); idmEG++) {
592 double dmxEGEt = dmxegs->at(ibx, idmEG).hwPt();
595 double gtEGEta = 0.0;
596 double dmxEGEta = dmxegs->at(ibx, idmEG).hwEta();
599 double gtEGPhi = 0.0;
600 double dmxEGPhi = dmxegs->at(ibx, idmEG).hwPhi();
614 for (
int ibx =
taus->getFirstBX(); ibx <=
taus->getLastBX(); ++ibx) {
616 for (
unsigned int igtTau = 0; igtTau <
taus->size(ibx); igtTau++) {
617 double gtTauEt =
taus->at(ibx, igtTau).hwPt();
618 double dmxTauEt = 0.0;
619 if (dmxtaus->size(ibx) > igtTau)
620 dmxTauEt = dmxtaus->at(ibx, igtTau).hwPt();
623 double gtTauEta =
taus->at(ibx, igtTau).hwEta();
624 double dmxTauEta = 0.0;
625 if (dmxtaus->size(ibx) > igtTau)
626 dmxTauEta = dmxtaus->at(ibx, igtTau).hwEta();
629 double gtTauPhi =
taus->at(ibx, igtTau).hwPhi();
630 double dmxTauPhi = 0.0;
631 if (dmxtaus->size(ibx) > igtTau)
632 dmxTauPhi = dmxtaus->at(ibx, igtTau).hwPhi();
636 for (
unsigned int idmTau =
taus->size(ibx); idmTau < dmxtaus->size(ibx); idmTau++) {
637 double gtTauEt = 0.0;
638 double dmxTauEt = dmxtaus->at(ibx, idmTau).hwPt();
641 double gtTauEta = 0.0;
642 double dmxTauEta = dmxtaus->at(ibx, idmTau).hwEta();
645 double gtTauPhi = 0.0;
646 double dmxTauPhi = dmxtaus->at(ibx, idmTau).hwPhi();
660 for (
int ibx =
jets->getFirstBX(); ibx <=
jets->getLastBX(); ++ibx) {
662 for (
unsigned int igtJ = 0; igtJ <
jets->size(ibx); igtJ++) {
663 double gtJetEt =
jets->at(ibx, igtJ).hwPt();
664 double dmxJetEt = 0.0;
665 if (dmxjets->size(ibx) > igtJ)
666 dmxJetEt = dmxjets->at(ibx, igtJ).hwPt();
669 double gtJetEta =
jets->at(ibx, igtJ).hwEta();
670 double dmxJetEta = 0.0;
671 if (dmxjets->size(ibx) > igtJ)
672 dmxJetEta = dmxjets->at(ibx, igtJ).hwEta();
675 double gtJetPhi =
jets->at(ibx, igtJ).hwPhi();
676 double dmxJetPhi = 0.0;
677 if (dmxjets->size(ibx) > igtJ)
678 dmxJetPhi = dmxjets->at(ibx, igtJ).hwPhi();
682 for (
unsigned int idmJ =
jets->size(ibx); idmJ < dmxjets->size(ibx); idmJ++) {
683 double gtJetEt = 0.0;
684 double dmxJetEt = dmxjets->at(ibx, idmJ).hwPt();
687 double gtJetEta = 0.0;
688 double dmxJetEta = dmxjets->at(ibx, idmJ).hwEta();
691 double gtJetPhi = 0.0;
692 double dmxJetPhi = dmxjets->at(ibx, idmJ).hwPhi();
706 for (
int ibx = sums->getFirstBX(); ibx <= sums->getLastBX(); ++ibx) {
708 for (
unsigned int igtS = 0; igtS < sums->size(ibx); igtS++) {
709 double gtSumEt = sums->at(ibx, igtS).hwPt();
710 double dmxSumEt = 0.0;
711 if (dmxSums->size(ibx) > igtS)
712 dmxSumEt = dmxSums->at(ibx, igtS).hwPt();
714 double gtSumPhi = sums->at(ibx, igtS).hwPhi();
715 double dmxSumPhi = 0.0;
716 if (dmxSums->size(ibx) > igtS)
717 dmxSumPhi = dmxSums->at(ibx, igtS).hwPhi();
719 if (sums->at(ibx, igtS).getType() == dmxSums->at(ibx, igtS).getType()) {
720 switch (sums->at(ibx, igtS).getType()) {
761 text <<
"WARNING: EtSum Types do not line up between DeMux and uGT " << std::endl;
765 for (
unsigned int idmS = sums->size(ibx); idmS < dmxSums->size(ibx); idmS++) {
766 double gtSumEt = -1.0;
767 double dmxSumEt = dmxSums->at(ibx, idmS).hwPt();
769 double gtSumPhi = -1.0;
770 double dmxSumPhi = dmxSums->at(ibx, idmS).hwPhi();
772 switch (dmxSums->at(ibx, idmS).getType()) {
815 auto itr = hwalgs->begin(ibx);
816 auto itr_emul = emulAlgs->begin(ibx);
818 for (
int algBit = 0; algBit < 128;
820 bool hw = itr->getAlgoDecisionFinal(algBit);
821 bool em = itr_emul->getAlgoDecisionFinal(algBit);
824 }
else if (hw & !em) {
826 text <<
"WARNING: HW Fnd Alg Bit " << algBit <<
" but emulation did not " << std::endl;
827 }
else if (!hw & em) {
829 text <<
"WARNING: Emul. Fnd Alg Bit " << algBit <<
" but hardware did not " << std::endl;
847 auto itr = hwalgs->begin(ibx);
848 auto itr_emul = emulAlgs->begin(ibx);
850 for (
int algBit = 0; algBit < 128;
852 bool hw = itr->getAlgoDecisionFinal(algBit);
853 bool em = itr_emul->getAlgoDecisionFinal(algBit);
856 }
else if (hw & !em) {
858 }
else if (!hw & em) {
TH2F * hAlgoBitsEmulGtVsHw_
TH2F * hDmxVsGTSumEt_ETMHF_
edm::EDGetToken m_emulDxAlgToken
std::map< ObjectType, TH1F * > het_
edm::EDGetToken m_dmxEGToken
TH2F * hDmxVsGTSumEt_HTM_
edm::EDGetToken m_egToken
TH2F * hDmxVsGTSumPhi_ETM_
std::map< ObjectType, TH1F * > hphi_
muons
the two sets of parameters below are mutually exclusive, depending if RECO or ALCARECO is used the us...
edm::EDGetToken m_emulGtAlgToken
edm::EDGetToken m_muToken
TH2F * hDmxVsGTSumEt_HFP1_
TH2F * hDmxVsGTSumEt_ETT_
edm::EDGetToken m_gmuToken
TH2F * hDmxVsGTSumPhi_ETMHF_
edm::EDGetToken m_dmxTauToken
edm::EDGetToken m_gtAlgToken
TH2F * hDmxVsGTSumEt_ETM_
edm::EDGetToken m_sumToken
std::map< ObjectType, TH1F * > hbx_
Log< level::Info, false > LogInfo
std::map< ObjectType, TH2F * > hetaphi_
TH2F * hDmxVsGTSumEt_HTT_
TH2F * hAlgoBitsEmulDxVsHw_
TH2F * hDmxVsGTSumEt_HFM0_
TH2F * hDmxVsGTSumEt_HFM1_
edm::EDGetToken m_dmxSumToken
edm::EDGetToken m_tauToken
TH2F * hDmxVsGTSumEt_ETTem_
TH2F * hDmxVsGTSumPhi_HTM_
edm::EDGetToken m_jetToken
edm::EDGetToken m_dmxJetToken
TH2F * hDmxVsGTSumEt_HFP0_
std::map< ObjectType, TH1F * > heta_