|
ATLAS Offline Software
|
Go to the documentation of this file.
16 #include "GaudiKernel/MsgStream.h"
17 #include "GaudiKernel/StatusCode.h"
28 #include "CaloDetDescr/CaloDetDescrElement.h"
31 #include "Identifier/Identifier.h"
49 const std::string &
name,
52 m_ttTool(
"LVL1::L1TriggerTowerTool/L1TriggerTowerTool"),
53 m_errorTool(
"LVL1::TrigT1CaloMonErrorTool/TrigT1CaloMonErrorTool"),
54 m_histTool(
"LVL1::TrigT1CaloLWHistogramTool/TrigT1CaloLWHistogramTool"),
55 m_caloTool(
"LVL1::L1CaloMonitoringCaloTool/L1CaloMonitoringCaloTool"),
57 m_h_CaloCell_tile_phi(0),
58 m_h_CaloCell_tile_eta(0),
59 m_h_CaloCell_tile_HitMap(0),
60 m_h_CaloCell_tile_et(0),
61 m_h_CaloCell_lar_phi(0),
62 m_h_CaloCell_lar_eta(0),
63 m_h_CaloCell_lar_HitMap(0),
64 m_h_CaloCell_lar_et(0),
65 m_h_CaloTT_HitMap_emLUT_Thresh0(0),
66 m_h_CaloTT_HitMap_emLUT_Thresh1(0),
67 m_h_CaloTT_HitMap_emLUT_Thresh2(0),
68 m_h_CaloTT_HitMap_hadLUT_Thresh0(0),
69 m_h_CaloTT_HitMap_hadLUT_Thresh1(0),
70 m_h_CaloTT_HitMap_hadLUT_Thresh2(0),
72 m_h_CaloTT_emLUT_eta(0),
73 m_h_CaloTT_emLUT_phi(0),
75 m_h_CaloTT_hadLUT_eta(0),
76 m_h_CaloTT_hadLUT_phi(0),
77 m_h_emEnergy_Match_CaloTT_TT_LArEMB(0),
78 m_h_emEnergy_Match_CaloTT_TT_LArOverlap(0),
79 m_h_emEnergy_Match_CaloTT_TT_LArEMEC(0),
80 m_h_emEnergy_Match_CaloTT_TT_LArFCAL1(0),
81 m_h_hadEnergy_Match_CaloTT_TT_TileLB(0),
82 m_h_hadEnergy_Match_CaloTT_TT_TileEB(0),
83 m_h_hadEnergy_Match_CaloTT_TT_LArHEC(0),
84 m_h_hadEnergy_Match_CaloTT_TT_LArFCAL23(0),
85 m_h_emEnergy_Match_CaloTT_TT_LArEMB_CP(0),
86 m_h_emEnergy_Match_CaloTT_TT_LArOverlap_CP(0),
87 m_h_emEnergy_Match_CaloTT_TT_LArEMEC_CP(0),
88 m_h_emEnergy_Match_CaloTT_TT_LArFCAL1_CP(0),
89 m_h_hadEnergy_Match_CaloTT_TT_TileLB_CP(0),
90 m_h_hadEnergy_Match_CaloTT_TT_TileEB_CP(0),
91 m_h_hadEnergy_Match_CaloTT_TT_LArHEC_CP(0),
92 m_h_hadEnergy_Match_CaloTT_TT_LArFCAL23_CP(0),
93 m_h_em_profile_Match(0),
94 m_h_had_profile_Match(0),
95 m_h_em_profile_Match_CP(0),
96 m_h_had_profile_Match_CP(0),
105 m_h_em_1D_Match_LArEMB(0),
106 m_h_em_1D_Match_LArOverlap(0),
107 m_h_em_1D_Match_LArEMEC(0),
108 m_h_em_1D_Match_LArFCAL1(0),
109 m_h_had_1D_Match_TileLB(0),
110 m_h_had_1D_Match_TileEB(0),
111 m_h_had_1D_Match_LArHEC(0),
112 m_h_had_1D_Match_LArFCAL23(0),
113 m_h_em_1D_Match_LArEMB_CP(0),
114 m_h_em_1D_Match_LArOverlap_CP(0),
115 m_h_em_1D_Match_LArEMEC_CP(0),
116 m_h_em_1D_Match_LArFCAL1_CP(0),
117 m_h_had_1D_Match_TileLB_CP(0),
118 m_h_had_1D_Match_TileEB_CP(0),
119 m_h_had_1D_Match_LArHEC_CP(0),
120 m_h_had_1D_Match_LArFCAL23_CP(0),
121 m_h_emEnergy_Calo_L1CaloSat_LArEMB(0),
122 m_h_emEnergy_Calo_L1CaloSat_LArOverlap(0),
123 m_h_emEnergy_Calo_L1CaloSat_LArEMEC(0),
124 m_h_emEnergy_Calo_L1CaloSat_LArFCAL1(0),
125 m_h_hadEnergy_Calo_L1CaloSat_TileLB(0),
126 m_h_hadEnergy_Calo_L1CaloSat_TileEB(0),
127 m_h_hadEnergy_Calo_L1CaloSat_LArHEC(0),
128 m_h_hadEnergy_Calo_L1CaloSat_LArFCAL23(0),
129 m_h_emEnergy_Calo_L1CaloSat_LArEMB_CP(0),
130 m_h_emEnergy_Calo_L1CaloSat_LArOverlap_CP(0),
131 m_h_emEnergy_Calo_L1CaloSat_LArEMEC_CP(0),
132 m_h_emEnergy_Calo_L1CaloSat_LArFCAL1_CP(0),
133 m_h_hadEnergy_Calo_L1CaloSat_TileLB_CP(0),
134 m_h_hadEnergy_Calo_L1CaloSat_TileEB_CP(0),
135 m_h_hadEnergy_Calo_L1CaloSat_LArHEC_CP(0),
136 m_h_hadEnergy_Calo_L1CaloSat_LArFCAL23_CP(0),
137 m_h_em_profile_etaRegion(0),
138 m_h_had_profile_etaRegion(0),
139 m_h_em_profile_phiRegion(0),
140 m_h_had_profile_phiRegion(0),
141 m_h_em_profile_etaRegion_CP(0),
142 m_h_had_profile_etaRegion_CP(0),
143 m_h_em_profile_phiRegion_CP(0),
144 m_h_had_profile_phiRegion_CP(0),
145 m_h_average_emDE_map(0),
146 m_h_average_hadDE_map(0),
147 m_h_average_emDE_map_CP(0),
148 m_h_average_hadDE_map_CP(0),
149 m_h_emTTME_etaphi(0),
150 m_h_hadTTME_etaphi(0),
151 m_h_emTTME_etaphi_CP(0),
152 m_h_hadTTME_etaphi_CP(0),
153 m_h_emcaloME_etaphi(0),
154 m_h_hadcaloME_etaphi(0)
185 if (
sc.isFailure())
return sc;
188 if (
sc.isFailure()) {
189 msg(MSG::ERROR) <<
"Unable to locate tool L1TriggerTowerTool" <<
endmsg;
194 if(
sc.isFailure() ) {
195 msg(MSG::ERROR) <<
"Unable to locate Tool TrigT1CaloMonErrorTool"
201 if (
sc.isFailure()) {
202 msg(MSG::ERROR) <<
"Unable to locate Tool TrigT1CaloLWHistogramTool"
208 if (
sc.isFailure()) {
209 msg(MSG::ERROR) <<
"Unable to locate Tool L1CaloMonitoringCaloTool"
214 return StatusCode::SUCCESS;
222 return StatusCode::SUCCESS;
258 "Distribution of Tile CaloCell Energy;energy",
261 "#eta - Distribution of Tile CaloCells;eta", 98, -4.9, 4.9);
263 "#phi - Distribution of Tile CaloCells;phi", 64, 0., 2*
M_PI);
265 "#eta - #phi Map of Tile CaloCells;eta;phi", 98, -4.9, 4.9, 64, 0., 2*
M_PI);
268 "Distribution of LAr CaloCell Energy;energy",
271 "#eta - Distribution of LAr CaloCells;eta", 98, -4.9, 4.9);
273 "#phi - Distribution of LAr CaloCells;phi", 64, 0., 2*
M_PI);
275 "#eta - #phi Map of LAr CaloCells;eta;phi", 98, -4.9, 4.9, 64, 0., 2*
M_PI);
285 "CaloTTEmLUT"+
buffer.str(),
286 "#eta - #phi Map of Calo EM E_{T} > "+
buffer.str());
291 "CaloTTEmLUT"+
buffer.str(),
292 "#eta - #phi Map of Calo EM E_{T} > "+
buffer.str());
297 "CaloTTEmLUT"+
buffer.str(),
298 "#eta - #phi Map of Calo EM E_{T} > "+
buffer.str());
303 "CaloTTHadLUT"+
buffer.str(),
304 "#eta - #phi Map of Calo Had E_{T} > "+
buffer.str());
309 "CaloTTHadLUT"+
buffer.str(),
310 "#eta - #phi Map of Calo Had E_{T} > "+
buffer.str());
315 "CaloTTHadLUT"+
buffer.str(),
316 "#eta - #phi Map of Calo Had E_{T} > "+
buffer.str());
325 "Calo EM #eta Distribution");
327 "Calo EM #phi Distribution;phi", 64, 0., 2*
M_PI);
332 "Calo HAD #eta Distribution");
334 "Calo HAD #phi Distribution;phi",64,0,2*
M_PI);
341 "emEnergy_CompTT2CaloTT_LArEMB",
342 "LArEMB Match between Calo and L1Calo em E_{T} (JEP);L1Calo;Calo",
346 "emEnergy_CompTT2CaloTT_LArOverlap",
347 "LArOverlap Match between Calo and L1Calo em E_{T} (JEP);L1Calo;Calo",
351 "emEnergy_CompTT2CaloTT_LArEMEC",
352 "LArEMEC Match between Calo and L1Calo em E_{T} (JEP);L1Calo;Calo",
356 "emEnergy_CompTT2CaloTT_LArFCAL1",
357 "LArFCAL1 Match between Calo and L1Calo em E_{T} (JEP);L1Calo;Calo",
362 "hadEnergy_CompTT2CaloTT_TileLB",
363 "TileLB Match between Calo and L1Calo had E_{T} (JEP);L1Calo;Calo",
367 "hadEnergy_CompTT2CaloTT_TileEB",
368 "TileEB Match between Calo and L1Calo had E_{T} (JEP);L1Calo;Calo",
372 "hadEnergy_CompTT2CaloTT_LArHEC",
373 "LArHEC Match between Calo and L1Calo had E_{T} (JEP);L1Calo;Calo",
377 "hadEnergy_CompTT2CaloTT_LArFCAL23",
378 "LArFCAL23 Match between Calo and L1Calo had E_{T} (JEP);L1Calo;Calo",
383 "emEnergy_CompTT2CaloTT_LArEMB_CP",
384 "LArEMB Match between Calo and L1Calo em E_{T} (CP);L1Calo;Calo",
388 "emEnergy_CompTT2CaloTT_LArOverlap_CP",
389 "LArOverlap Match between Calo and L1Calo em E_{T} (CP);L1Calo;Calo",
393 "emEnergy_CompTT2CaloTT_LArEMEC_CP",
394 "LArEMEC Match between Calo and L1Calo em E_{T} (CP);L1Calo;Calo",
398 "emEnergy_CompTT2CaloTT_LArFCAL1_CP",
399 "LArFCAL1 Match between Calo and L1Calo em E_{T} (CP);L1Calo;Calo",
404 "hadEnergy_CompTT2CaloTT_TileLB_CP",
405 "TileLB Match between Calo and L1Calo had E_{T} (CP);L1Calo;Calo",
409 "hadEnergy_CompTT2CaloTT_TileEB_CP",
410 "TileEB Match between Calo and L1Calo had E_{T} (CP);L1Calo;Calo",
414 "hadEnergy_CompTT2CaloTT_LArHEC_CP",
415 "LArHEC Match between Calo and L1Calo had E_{T} (CP);L1Calo;Calo",
419 "hadEnergy_CompTT2CaloTT_LArFCAL23_CP",
420 "LArFCAL23 Match between Calo and L1Calo had E_{T} (CP);L1Calo;Calo",
425 "emEnergy_eta_DiffTT2CaloTT",
426 "Relative difference in eta: (L1Calo (JEP) - Calo) / Calo em E_{T};eta;(L1Calo-Calo)/Calo Et", 98, -4.9, 4.9);
428 "hadEnergy_eta_DiffTT2CaloTT",
429 "Relative difference in eta: (L1Calo (JEP) - Calo) / Calo had E_{T};eta;(L1Calo-Calo)/Calo Et", 98, -4.9, 4.9);
432 "emEnergy_eta_DiffTT2CaloTT_CP",
433 "Relative difference in eta: (L1Calo (CP) - Calo) / Calo em E_{T};eta;(L1Calo-Calo)/Calo Et", 98, -4.9, 4.9);
435 "hadEnergy_eta_DiffTT2CaloTT_CP",
436 "Relative difference in eta: (L1Calo (CP) - Calo) / Calo had E_{T};eta;(L1Calo-Calo)/Calo Et", 98, -4.9, 4.9);
465 "emEnergy_DiffCaloL1CaloE_LArEMB",
466 "LArEMB relative difference between L1Calo (JEP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
468 "emEnergy_DiffCaloL1CaloE_LArOverlap",
469 "LArOverlap relative difference between L1Calo (JEP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
471 "emEnergy_DiffCaloL1CaloE_LArEMEC",
472 "LArEMEC relative difference between L1Calo (JEP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
474 "emEnergy_DiffCaloL1CaloE_LArFCAL1",
475 "LArFCAL1 relative difference between L1Calo (JEP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
478 "hadEnergy_DiffCaloL1CaloE_TileLB",
479 "TileLB relative difference between L1Calo (JEP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
481 "hadEnergy_DiffCaloL1CaloE_TileEB",
482 "TileEB relative difference between L1Calo (JEP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
484 "hadEnergy_DiffCaloL1CaloE_LArHEC",
485 "LArHEC relative difference between L1Calo (JEP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
487 "hadEnergy_DiffCaloL1CaloE_LArFCAL23",
488 "LArFCAL23 relative difference between L1Calo (JEP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
491 "emEnergy_DiffCaloL1CaloE_LArEMB_CP",
492 "LArEMB relative difference between L1Calo (CP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
494 "emEnergy_DiffCaloL1CaloE_LArOverlap_CP",
495 "LArOverlap relative difference between L1Calo (CP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
497 "emEnergy_DiffCaloL1CaloE_LArEMEC_CP",
498 "LArEMEC relative difference between L1Calo (CP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
500 "emEnergy_DiffCaloL1CaloE_LArFCAL1_CP",
501 "LArFCAL1 relative difference between L1Calo (CP) and Calo em E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
504 "hadEnergy_DiffCaloL1CaloE_TileLB_CP",
505 "TileLB relative difference between L1Calo (CP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
507 "hadEnergy_DiffCaloL1CaloE_TileEB_CP",
508 "TileEB relative difference between L1Calo (CP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
510 "hadEnergy_DiffCaloL1CaloE_LArHEC_CP",
511 "LArHEC relative difference between L1Calo (CP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
513 "hadEnergy_DiffCaloL1CaloE_LArFCAL23_CP",
514 "LArFCAL23 relative difference between L1Calo (CP) and Calo had E_{T};(L1Calo-Calo)/Calo Et", 24, -0.6, 0.6);
517 "Relative difference between L1Calo (JEP) and Calo em E_{T};L1Calo Et;(L1Calo-Calo)/Calo Et",
520 "Relative difference between L1Calo (JEP) and Calo had E_{T};L1Calo Et;(L1Calo-Calo)/Calo Et",
523 "emEnergy_phi_DiffTT2CaloTT",
524 "Relative difference in phi: L1Calo (JEP) and Calo em E_{T};phi;(L1Calo-Calo)/Calo Et", 64, 0., 2.*
M_PI);
526 "hadEnergy_phi_DiffTT2CaloTT",
527 "Relative difference in phi: L1Calo (JEP) and Calo had E_{T};phi;(L1Calo-Calo)/Calo Et", 64, 0., 2.*
M_PI);
530 "Relative difference between L1Calo (CP) and Calo em E_{T};L1Calo Et;(L1Calo-Calo)/Calo Et",
533 "Relative difference between L1Calo (CP) and Calo had E_{T};L1Calo Et;(L1Calo-Calo)/Calo Et",
536 "emEnergy_phi_DiffTT2CaloTT_CP",
537 "Relative difference in phi: L1Calo (CP) and Calo em E_{T};phi;(L1Calo-Calo)/Calo Et", 64, 0., 2.*
M_PI);
539 "hadEnergy_phi_DiffTT2CaloTT_CP",
540 "Relative difference in phi: L1Calo (CP) and Calo had E_{T};phi;(L1Calo-Calo)/Calo Et", 64, 0., 2.*
M_PI);
544 "#eta-#phi map of average (L1Calo (JEP) - Calo) / Calo em E_{T} relative difference");
546 "hadDeltaEnergy_Map",
547 "#eta-#phi map of average (L1Calo (JEP) - Calo) / Calo had E_{T} relative difference");
550 "emDeltaEnergy_Map_CP",
551 "#eta-#phi map of average (L1Calo (CP) - Calo) / Calo em E_{T} relative difference");
553 "hadDeltaEnergy_Map_CP",
554 "#eta-#phi map of average (L1Calo (CP) - Calo) / Calo had E_{T} relative difference");
557 "Distribution of LArEMB Calo em Et for L1Calo JEP Et = 255;energy",
560 "Distribution of LArOverlap Calo em Et for L1Calo JEP Et = 255;energy",
563 "Distribution of LArEMEC Calo em Et for L1Calo JEP Et = 255;energy",
566 "Distribution of LArFCAL1 Calo em Et for L1Calo JEP Et = 255;energy",
569 "Distribution of TileLB Calo had Et for L1Calo JEP Et = 255;energy",
572 "Distribution of TileEB Calo had Et for L1Calo JEP Et = 255;energy",
575 "Distribution of LArHEC Calo had Et for L1Calo JEP Et = 255;energy",
578 "Distribution of LArFCAL23 Calo had Et for L1Calo JEP Et = 255;energy",
582 "Distribution of LArEMB Calo em Et for L1Calo CP Et = 255;energy",
585 "Distribution of LArOverlap Calo em Et for L1Calo CP Et = 255;energy",
588 "Distribution of LArEMEC Calo em Et for L1Calo CP Et = 255;energy",
591 "Distribution of LArFCAL1 Calo em Et for L1Calo CP Et = 255;energy",
594 "Distribution of TileLB Calo had Et for L1Calo CP Et = 255;energy",
597 "Distribution of TileEB Calo had Et for L1Calo CP Et = 255;energy",
600 "Distribution of LArHEC Calo had Et for L1Calo CP Et = 255;energy",
603 "Distribution of LArFCAL23 Calo had Et for L1Calo CP Et = 255;energy",
607 "em_MostEnergetic_TTMap",
"#eta-#phi map of em most energetic TT (JEP)");
609 "had_MostEnergetic_TTMap",
"#eta-#phi map of had most energetic TT (JEP)");
612 "em_MostEnergetic_TTMap_CP",
"#eta-#phi map of em most energetic TT (CP)");
614 "had_MostEnergetic_TTMap_CP",
"#eta-#phi map of had most energetic TT (CP)");
617 "em_MostEnergetic_caloMap",
"#eta-#phi map of em most energetic caloTT");
619 "had_MostEnergetic_caloMap",
620 "#eta-#phi map of had most energetic caloTT");
626 if(!
sc.isSuccess()) {
627 msg(MSG::WARNING) <<
"Conditions not retrieved " <<
endmsg;
632 return StatusCode::SUCCESS;
645 return StatusCode::SUCCESS;
652 return StatusCode::SUCCESS;
661 if(!
sc.isSuccess() || !caloCellContainer) {
662 msg(MSG::WARNING) <<
"No CaloCellContainer found at AllCalo" <<
endmsg;
663 return StatusCode::SUCCESS;
666 if(!
sc.isSuccess()) {
667 msg(MSG::WARNING) <<
"Could not load CaloCells" <<
endmsg;
668 return StatusCode::SUCCESS;
681 for (; CaloCellIterator != CaloCellIteratorEnd; ++CaloCellIterator) {
684 double energy = (*CaloCellIterator)->energy()*1
e-3;
712 if(
sc==StatusCode::FAILURE || !triggerTowerTES) {
713 msg(MSG::INFO) <<
"No xAODTriggerTower found in TES at "
715 return StatusCode::SUCCESS;
719 triggerTowerTES->
begin();
721 triggerTowerTES->
end();
723 double em_maximum = 0;
724 double had_maximum = 0;
725 double em_maximum_CP = 0;
726 double had_maximum_CP = 0;
727 double caloem_maximum = 0;
728 double calohad_maximum = 0;
730 double eta_max_had = -10;
731 double eta_max_em = -10;
732 double phi_max_had = -10;
733 double phi_max_em = -10;
734 double eta_max_had_CP = -10;
735 double eta_max_em_CP = -10;
736 double phi_max_had_CP = -10;
737 double phi_max_em_CP = -10;
738 double eta_max_calohad = -10;
739 double eta_max_caloem = -10;
740 double phi_max_calohad = -10;
741 double phi_max_caloem = -10;
747 double had_caloE = 0;
748 double em_L1CaloE = 0;
749 double had_L1CaloE = 0;
750 double em_L1CaloE_CP = 0;
751 double had_L1CaloE_CP = 0;
754 CaloCellEnergyByLayerAcc(
"CaloCellEnergyByLayer");
756 CaloCellETByLayerAcc(
"CaloCellETByLayer");
758 for (; ttIterator != ttIteratorEnd; ++ttIterator) {
759 double eta = (*ttIterator)->eta();
761 double phi = (*ttIterator)->phi();
762 const int layer = (*ttIterator)->layer();
772 if (!(
tt->jepET() ||
tt->cpET()))
continue;
779 em_L1CaloE =
int(
tt->jepET());
780 em_L1CaloE_CP =
int(0.5*(
tt->cpET()));
785 for (
float e : CaloCellEnergyByLayerAcc(*
tt)) {
788 em_caloE = em_caloE / cosh(
eta);
792 for (
float et : CaloCellETByLayerAcc(*
tt)) {
798 if (em_caloE < 0) em_caloE = 0;
799 em_caloE =
int(em_caloE+0.5);
800 if (em_caloE > 255) em_caloE = 255;
804 double emRelDiffE = 0.;
807 if (em_L1CaloE > 5 && em_caloE > 5) {
809 bool em_disabled =
m_ttTool->disabledChannel(em_coolId);
812 emDiffE = em_L1CaloE-em_caloE;
815 emRelDiffE = emDiffE / em_caloE;
832 if (em_L1CaloE >= 255) {
838 }
else if (
absEta < 1.5) {
839 if (em_L1CaloE >= 255) {
845 }
else if (
absEta < 3.2) {
846 if (em_L1CaloE >= 255) {
852 }
else if (
absEta < 4.9) {
853 if (em_L1CaloE >= 255) {
865 double emDiffE_CP = 0.;
866 double emRelDiffE_CP = 0.;
869 if (em_L1CaloE_CP > 5 && em_caloE > 5) {
871 bool em_disabled =
m_ttTool->disabledChannel(em_coolId);
874 emDiffE_CP = em_L1CaloE_CP-em_caloE;
877 emRelDiffE_CP = emDiffE_CP / em_caloE;
895 if (em_L1CaloE_CP >= 127) {
901 }
else if (
absEta < 1.5) {
902 if (em_L1CaloE_CP >= 127) {
908 }
else if (
absEta < 3.2) {
909 if (em_L1CaloE_CP >= 127) {
915 }
else if (
absEta < 4.9) {
916 if (em_L1CaloE_CP >= 127) {
926 if (em_maximum < em_L1CaloE) {
927 em_maximum = em_L1CaloE;
932 if (em_maximum_CP < em_L1CaloE_CP) {
933 em_maximum_CP = em_L1CaloE_CP;
938 if(caloem_maximum < em_caloE) {
939 caloem_maximum = em_caloE;
940 eta_max_caloem =
eta;
941 phi_max_caloem =
phi;
966 had_L1CaloE =
int(
tt->jepET());
967 had_L1CaloE_CP =
int(0.5*(
tt->cpET()));
971 for (
float e : CaloCellEnergyByLayerAcc(*
tt)) {
974 had_caloE = had_caloE / cosh(
eta);
978 for (
float et : CaloCellETByLayerAcc(*
tt)) {
984 if (had_caloE < 0) had_caloE = 0;
985 had_caloE =
int(had_caloE+0.5);
986 if (had_caloE > 255) had_caloE = 255;
990 if (had_L1CaloE > 5 && had_caloE > 5) {
992 bool had_disabled =
m_ttTool->disabledChannel(had_coolId);
995 double hadDiffE = had_L1CaloE-had_caloE;
996 double hadRelDiffE = 0.;
998 if (had_caloE != 0) {
999 hadRelDiffE = hadDiffE / had_caloE;
1016 if (had_L1CaloE >= 255) {
1022 }
else if (
absEta < 1.5) {
1023 if (had_L1CaloE >= 255) {
1029 }
else if (
absEta < 3.2) {
1030 if (had_L1CaloE >= 255) {
1036 }
else if (
absEta < 4.9) {
1037 if (had_L1CaloE >= 255) {
1049 if (had_L1CaloE_CP > 5 && had_caloE > 5) {
1051 bool had_disabled =
m_ttTool->disabledChannel(had_coolId);
1052 if (!had_disabled) {
1054 double hadDiffE_CP = had_L1CaloE_CP-had_caloE;
1055 double hadRelDiffE_CP = 0;
1057 if (had_caloE != 0) {
1058 hadRelDiffE_CP = hadDiffE_CP / had_caloE;
1074 if (had_L1CaloE_CP >= 127) {
1080 }
else if (
absEta < 1.5) {
1081 if (had_L1CaloE_CP >= 127) {
1087 }
else if (
absEta < 3.2) {
1088 if (had_L1CaloE_CP >= 127) {
1094 }
else if (
absEta < 4.9) {
1095 if (had_L1CaloE_CP >= 127) {
1105 if (had_maximum < had_L1CaloE) {
1106 had_maximum = had_L1CaloE;
1111 if (had_maximum_CP < had_L1CaloE_CP) {
1112 had_maximum_CP = had_L1CaloE_CP;
1113 eta_max_had_CP =
eta;
1114 phi_max_had_CP =
phi;
1117 if (calohad_maximum < had_caloE) {
1118 calohad_maximum = had_caloE;
1119 eta_max_calohad =
eta;
1120 phi_max_calohad =
phi;
1123 if (had_caloE > 0) {
1143 if (eta_max_caloem != -10) {
1147 if (eta_max_calohad != -10) {
1151 if (eta_max_em != -10) {
1154 if (eta_max_had != -10) {
1157 if (eta_max_em_CP != -10) {
1160 if (eta_max_had_CP != -10) {
1164 return StatusCode::SUCCESS;
1174 return StatusCode::SUCCESS;
TH1F_LW * m_h_had_1D_Match_LArFCAL23
LArFCAL23 Relative difference between L1Calo (JEP) and Calo had Et.
virtual void Fill(const double &x, const double &y) override
std::string m_caloCellContainerName
CaloCellContainer StoreGate key.
TH1F_LW * m_h_had_1D_Match_TileEB
TileEB Relative difference between L1Calo (JEP) and Calo had Et.
TH1F_LW * m_h_had_1D_Match_TileLB
TileLB Relative difference between L1Calo (JEP) and Calo had Et.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_TileEB
TileEB Match between Calo and L1Calo (JEP) had Et.
ToolHandle< LVL1::IL1TriggerTowerTool > m_ttTool
Tool for identifiers and disabled channels.
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_TileEB_CP
TileEB Et for Calo TT with L1Calo TT CP ET = 255.
Extra patterns decribing particle interation process.
TProfile_LW * m_h_em_profile_etaRegion_CP
Relative difference in eta between L1Calo (CP) and Calo em Et.
Const iterator class for DataVector/DataList.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_LArFCAL23_CP
LArFCAL23 Match between Calo and L1Calo (CP) had Et.
Scalar phi() const
phi method
virtual void Fill(const double &x) override
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArEMB
LArEMB Match between Calo and L1Calo (JEP) em Et
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArFCAL1
LArFCAL1 Match between Calo and L1Calo (JEP) em Et.
TH2F_LW * m_h_CaloTT_HitMap_hadLUT_Thresh1
eta - phi Map of Calo Had Et > Thresh1
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArEMB_CP
LArEMB Match between Calo and L1Calo (CP) em Et.
TH1F_LW * m_h_had_1D_Match_LArHEC
LArHEC Relative difference between L1Calo (JEP) and Calo had Et.
TH1F_LW * m_h_em_1D_Match_LArFCAL1
LArFCAL1 Relative difference between LCalo (JEP) and Calo em Et.
bool m_histBooked
Histograms booked flag.
TProfile2D_LW * m_h_average_emDE_map
eta-phi map of average L1Calo (JEP) - Calo em Et relative difference
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArEMEC_CP
LArEMEC Match between Calo and L1Calo (CP) em Et.
TH1F_LW * m_h_CaloTT_emLUT_eta
Calo EM eta Distribution.
TH2F_LW * m_h_hadcaloME_etaphi
eta-phi map of had most energetic caloTT
TH1F_LW * m_h_had_1D_Match_LArFCAL23_CP
LArFCAL23 Relative difference between L1Calo (CP) and Calo had Et.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_TileLB_CP
TileLB Match between Calo and L1Calo (CP) had Et.
TProfile_LW * m_h_em_profile_etaRegion
Relative difference in eta between L1Calo (JEP) and Calo em Et.
Scalar eta() const
pseudorapidity method
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
TProfile_LW * m_h_had_profile_etaRegion_CP
Relative difference in eta between L1Calo (CP) and Calo had Et.
double m_CaloThreshold
Minimum CaloCell energy in GeV.
TH2F_LW * m_h_emcaloME_etaphi
eta-phi map of em most energetic caloTT
ToolHandle< ITrigT1CaloMonErrorTool > m_errorTool
Event veto tool.
TH1F_LW * m_h_had_1D_Match_TileEB_CP
TileEB Relative difference between L1Calo (CP) and Calo had Et.
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_LArFCAL23
LArFCAL23 Et for Calo TT with L1Calo TT JEP ET = 255.
This class groups all DetDescr information related to a CaloCell. Provides a generic interface for al...
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_LArFCAL23_CP
LArFCAL23 Et for Calo TT with L1Calo TT CP ET = 255.
TProfile_LW * m_h_em_profile_Match
Relative difference between L1Calo (JEP) - Calo em Et.
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_LArHEC_CP
LArHEC Et for Calo TT with L1Calo TT CP ET = 255.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_LArHEC
LArHEC Match between Calo and L1Calo (JEP) had Et.
TProfile_LW * m_h_em_profile_phiRegion_CP
Relative difference in phi between L1Calo (CP) and Calo em Et.
bool msgLvl(const MSG::Level lvl) const
TH1F_LW * m_h_CaloCell_tile_phi
phi - Distribution of Tile CaloCells
TH1F_LW * m_h_CaloCell_lar_phi
phi - Distribution of LAr CaloCells
ToolHandle< TrigT1CaloLWHistogramTool > m_histTool
Histogram helper.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_LArFCAL23
LArFCAL23 Match between Calo and L1Calo (JEP) had Et.
Helper class to provide constant type-safe access to aux data.
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArFCAL1_CP
LArFCAL1 Match between Calo and L1Calo (CP) em Et.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArFCAL1_CP
LArFCAL1 Et for Calo TT with L1Calo TT CP ET = 255.
TH1F_LW * m_h_em_1D_Match_LArEMB
LArEMB Relative difference between L1Calo (JEP) and Calo em Et.
eFexTowerBuilder creates xAOD::eFexTowerContainer from supercells (LATOME) and triggerTowers (TREX) i...
TH1F_LW * m_h_em_1D_Match_LArOverlap_CP
LArOverlap Relative difference between L1Calo (CP) and Calo em Et.
TH1F_LW * m_h_CaloTT_emLUT
Calo EM Et Distribution.
int m_CaloTT_HitMap_Thresh0
Hitmaps per threshold first value.
TProfile2D_LW * m_h_average_hadDE_map_CP
eta-phi map of average L1Calo (CP) - Calo had Et relative difference
float eta_raw() const
cell eta_raw
bool usedInTT(const CaloDetDescrElement *caloDDE)
Return true if Calo cell is used in TT.
TH1F_LW * m_h_em_1D_Match_LArOverlap
LArOverlap Relative difference between L1Calo (JEP) and Calo em Et.
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
std::string m_PathInRootFile
Root directory.
virtual StatusCode finalize()
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_TileLB
TileLB Et for Calo TT with L1Calo TT JEP ET = 255.
TH1F_LW * m_h_CaloCell_lar_eta
eta - Distribution of LAr CaloCells
TH1F_LW * m_h_CaloTT_hadLUT_phi
Calo HAD phi Distribution.
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArOverlap_CP
LArOverlap Match between Calo and L1Calo (CP) em Et.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_TileLB
TileLB Match between Calo and L1Calo (JEP) had Et.
::StatusCode StatusCode
StatusCode definition for legacy code.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArOverlap
LAr Overlap Et for Calo TT with L1Calo TT JEP ET = 255.
std::string m_xAODTriggerTowerContainerName
TriggerTower collection StoreGate key.
virtual ~CalorimeterL1CaloMon()
TH2F_LW * m_h_CaloTT_HitMap_emLUT_Thresh0
eta - phi Map of Calo EM Et > Thresh0
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_TileEB
TileEB Et for Calo TT with L1Calo TT JEP ET = 255.
Description of TriggerTower_v2.
TH1F_LW * m_h_CaloCell_lar_et
Distribution of LAr CaloCell Energy.
TH1F_LW * m_h_had_1D_Match_TileLB_CP
TileLB Relative difference between L1Calo (CP) and Calo had Et.
TProfile_LW * m_h_had_profile_etaRegion
Relative difference in eta between L1Calo (JEP) and Calo had Et.
TProfile_LW * m_h_em_profile_phiRegion
Relative difference in phi between L1Calo (JEP) and Calo em Et.
virtual StatusCode initialize()
CalorimeterL1CaloMon(const std::string &type, const std::string &name, const IInterface *parent)
TH1F_LW * m_h_CaloTT_emLUT_phi
Calo EM phi Distribution.
TProfile_LW * m_h_had_profile_Match
Relative difference between L1Calo (JEP) - Calo had Et.
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArEMEC
LArEMEC Match between Calo and L1Calo (JEP) em Et.
TH2F_LW * m_h_emTTME_etaphi
eta-phi map of em most energetic TT
TH1F_LW * m_h_CaloCell_tile_et
Distribution of Tile CaloCell Energy.
TH2F_LW * m_h_CaloTT_HitMap_hadLUT_Thresh0
eta - phi Map of Calo Had Et > Thresh0
void Fill(const double &x, const double &y)
int m_CaloTT_HitMap_Thresh2
Hitmaps per threshold third value.
bool is_tile() const
cell belongs to Tile
TProfile_LW * m_h_had_profile_phiRegion_CP
Relative difference in phi between L1Calo (CP) and Calo had Et.
TProfile2D_LW * m_h_average_emDE_map_CP
eta-phi map of average L1Calo (CP) - Calo em Et relative difference
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArEMEC
LAr EMEC Et for Calo TT with L1Calo TT JEP ET = 255.
TH2F_LW * m_h_CaloCell_lar_HitMap
eta - phi Map of LAr CaloCells
virtual StatusCode bookHistogramsRecurrent()
An inheriting class should either override this function, bookHists() or bookHistograms().
TH2F_LW * m_h_hadTTME_etaphi_CP
eta-phi map of had most energetic TT (CP)
TProfile_LW * m_h_had_profile_Match_CP
Relative difference between L1Calo (CP) - Calo had Et.
TH1F_LW * m_h_CaloTT_hadLUT
Calo HAD Et Distribution.
TProfile_LW * m_h_em_profile_Match_CP
Relative difference between L1Calo (CP) - Calo em Et.
TProfile2D_LW * m_h_average_hadDE_map
eta-phi map of average L1Calo (JEP) - Calo had Et relative difference
int m_CaloTT_HitMap_Thresh1
Hitmaps per threshold second value.
ToolHandle< LVL1::IL1CaloMonitoringCaloTool > m_caloTool
Tool for CaloTT energies.
TH2F_LW * m_h_CaloTT_HitMap_emLUT_Thresh2
eta - phi Map of Calo EM Et > Thresh2
Container class for CaloCell.
static const std::string xAODTriggerTowerLocation
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_TileLB_CP
TileLB Et for Calo TT with L1Calo TT CP ET = 255.
TH2F_LW * m_h_emTTME_etaphi_CP
eta-phi map of em most energetic TT (CP)
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
TProfile_LW * m_h_had_profile_phiRegion
Relative difference in phi between L1Calo (JEP) and Calo had Et.
TH1F_LW * m_h_CaloCell_tile_eta
eta - Distribution of Tile CaloCells
CaloCell_ID::CaloSample getSampling() const
cell sampling
TH1F_LW * m_h_em_1D_Match_LArFCAL1_CP
LArFCAL1 Relative difference between L1Calo (CP) and Calo em Et.
TH2F_LW * m_h_emEnergy_Match_CaloTT_TT_LArOverlap
LArOverlap Match between Calo and L1Calo (JEP) em Et.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArFCAL1
LArFCAL1 Et for Calo TT with L1Calo TT JEP ET = 255.
bool absEta(const xAOD::TauJet &tau, double &out)
TH1F_LW * m_h_had_1D_Match_LArHEC_CP
LArHEC Relative difference between L1Calo (CP) and Calo had Et.
TH1F_LW * m_h_em_1D_Match_LArEMEC_CP
LArEMEC Relative difference between L1Calo (CP) and Calo em Et.
TH1F_LW * m_h_em_1D_Match_LArEMEC
LArEMEC Relative difference between L1Calo (JEP)and Calo em Et.
TH2F_LW * m_h_CaloCell_tile_HitMap
eta - phi Map of Tile CaloCells
TH2F_LW * m_h_hadTTME_etaphi
eta-phi map of had most energetic TT
bool isAvailable(const ELT &e) const
Test to see if this variable exists in the store.
TH1F_LW * m_h_CaloTT_hadLUT_eta
Calo HAD eta Distribution.
TH2F_LW * m_h_CaloTT_HitMap_emLUT_Thresh1
eta - phi Map of Calo EM Et > Thresh1
Helper class to provide constant type-safe access to aux data.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArOverlap_CP
LAr Overlap Et for Calo TT with L1Calo TT CP ET = 255.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArEMB_CP
LAr EMB Et for Calo TT with L1Calo TT CP ET = 255.
TH1F_LW * m_h_hadEnergy_Calo_L1CaloSat_LArHEC
LArHEC Et for Calo TT with L1Calo TT JEP ET = 255.
TH1F_LW * m_h_em_1D_Match_LArEMB_CP
LArEMB Relative difference between L1Calo (CP) and Calo em Et.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_LArHEC_CP
LArHEC Match between Calo and L1Calo (CP) had Et.
int m_MaxEnergyRange
Maximim energy in plots.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArEMEC_CP
LAr EMEC Et for Calo TT with L1Calo TT CP ET = 255.
TH2F_LW * m_h_hadEnergy_Match_CaloTT_TT_TileEB_CP
TileEB Match between Calo and L1Calo (CP) had Et.
virtual StatusCode procHistograms()
An inheriting class should either override this function or finalHists().
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
TH1F_LW * m_h_emEnergy_Calo_L1CaloSat_LArEMB
LAr EMB Et for Calo TT with L1Calo TT JEP ET = 255.
float phi_raw() const
cell phi_raw
virtual StatusCode fillHistograms()
An inheriting class should either override this function or fillHists().
TH2F_LW * m_h_CaloTT_HitMap_hadLUT_Thresh2
eta - phi Map of Calo Had Et > Thresh2