File indexing completed on 2024-10-30 00:11:33
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010 #include "FWCore/Framework/interface/MakerMacros.h"
0011
0012 #include "DataFormats/MuonReco/interface/MuonSegmentMatch.h"
0013 #include "DataFormats/MuonReco/interface/MuonGEMHitMatch.h"
0014
0015 #include "DQMOffline/MuonDPG/interface/BaseTnPEfficiencyTask.h"
0016
0017 #include "DataFormats/GEMDigi/interface/GEMVFATStatusCollection.h"
0018 #include "DataFormats/GEMDigi/interface/GEMOHStatusCollection.h"
0019 #include "DataFormats/GEMDigi/interface/GEMAMCStatusCollection.h"
0020 #include "Geometry/GEMGeometry/interface/GEMGeometry.h"
0021 #include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
0022 #include "FWCore/Framework/interface/ConsumesCollector.h"
0023 #include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
0024
0025 class GEMTnPEfficiencyTask : public BaseTnPEfficiencyTask {
0026 public:
0027
0028 GEMTnPEfficiencyTask(const edm::ParameterSet& config);
0029 uint16_t maskChamberWithError(const GEMDetId& chamber_id,
0030 const GEMOHStatusCollection*,
0031 const GEMVFATStatusCollection*);
0032 bool checkBounds(const GeomDet* geomDet, const GlobalPoint& global_position, const float bordercut);
0033
0034 ~GEMTnPEfficiencyTask() override;
0035
0036 const bool m_maskChamberWithError_;
0037 const edm::EDGetTokenT<GEMOHStatusCollection> m_GEMOHStatusCollectionToken_;
0038 const edm::EDGetTokenT<GEMVFATStatusCollection> m_GEMVFATStatusCollectionToken_;
0039 const edm::EDGetTokenT<GEMAMCStatusCollection> m_GEMAMCStatusCollectionToken_;
0040 std::unique_ptr<MuonServiceProxy> muon_service_;
0041
0042 protected:
0043 std::string topFolder() const override;
0044
0045 void bookHistograms(DQMStore::IBooker& iBooker, edm::Run const& run, edm::EventSetup const& context) override;
0046
0047
0048 void analyze(const edm::Event& event, const edm::EventSetup& context) override;
0049 };
0050
0051 GEMTnPEfficiencyTask::GEMTnPEfficiencyTask(const edm::ParameterSet& config)
0052 : BaseTnPEfficiencyTask(config),
0053 m_maskChamberWithError_((config.getUntrackedParameter<bool>("maskChamberWithError"))),
0054 m_GEMOHStatusCollectionToken_(
0055 consumes<GEMOHStatusCollection>(config.getUntrackedParameter<edm::InputTag>("ohStatusTag"))),
0056 m_GEMVFATStatusCollectionToken_(
0057 consumes<GEMVFATStatusCollection>(config.getUntrackedParameter<edm::InputTag>("vfatStatusTag"))),
0058 m_GEMAMCStatusCollectionToken_(
0059 consumes<GEMAMCStatusCollection>(config.getUntrackedParameter<edm::InputTag>("amcStatusTag"))) {
0060 LogTrace("DQMOffline|MuonDPG|GEMTnPEfficiencyTask") << "[GEMTnPEfficiencyTask]: Constructor" << std::endl;
0061 muon_service_ = std::make_unique<MuonServiceProxy>(config.getParameter<edm::ParameterSet>("ServiceParameters"),
0062 consumesCollector());
0063 }
0064
0065 GEMTnPEfficiencyTask::~GEMTnPEfficiencyTask() {
0066 LogTrace("DQMOffline|MuonDPG|GEMTnPEfficiencyTask")
0067 << "[GEMTnPEfficiencyTask]: analyzed " << m_nEvents << " events" << std::endl;
0068 }
0069
0070 void GEMTnPEfficiencyTask::bookHistograms(DQMStore::IBooker& iBooker,
0071 edm::Run const& run,
0072 edm::EventSetup const& context) {
0073 BaseTnPEfficiencyTask::bookHistograms(iBooker, run, context);
0074
0075 LogTrace("DQMOffline|MuonDPG|GEMTnPEfficiencyTask") << "[GEMTnPEfficiencyTask]: bookHistograms" << std::endl;
0076
0077 auto baseDir = topFolder() + "Task/";
0078 iBooker.setCurrentFolder(baseDir);
0079 MonitorElement* me_GE11_pass_Ch_region =
0080 iBooker.book2D("GE11_nPassingProbe_Ch_region", "GE11_nPassingProbe_Ch_region", 2, -1.5, 1.5, 36, 1, 37);
0081 MonitorElement* me_GE11_fail_Ch_region =
0082 iBooker.book2D("GE11_nFailingProbe_Ch_region", "GE11_nFailingProbe_Ch_region", 2, -1.5, 1.5, 36, 1, 37);
0083 MonitorElement* me_GE21_pass_Ch_region =
0084 iBooker.book2D("GE21_nPassingProbe_Ch_region", "GE21_nPassingProbe_Ch_region", 2, -1.5, 1.5, 36, 1, 37);
0085 MonitorElement* me_GE21_fail_Ch_region =
0086 iBooker.book2D("GE21_nFailingProbe_Ch_region", "GE21_nFailingProbe_Ch_region", 2, -1.5, 1.5, 36, 1, 37);
0087 MonitorElement* me_GEM_pass_Ch_region_GE1 =
0088 iBooker.book2D("GEM_nPassingProbe_Ch_region_GE1", "GEM_nPassingProbe_Ch_region_GE1", 4, 0, 4, 36, 1, 37);
0089 MonitorElement* me_GEM_fail_Ch_region_GE1 =
0090 iBooker.book2D("GEM_nFailingProbe_Ch_region_GE1", "GEM_nFailingProbe_Ch_region_GE1", 4, 0, 4, 36, 1, 37);
0091 MonitorElement* me_GEM_pass_Ch_region_GE1_NoL =
0092 iBooker.book2D("GEM_nPassingProbe_Ch_region_GE1_NoL", "GEM_nPassingProbe_Ch_region_GE1_NoL", 2, 0, 2, 36, 1, 37);
0093 MonitorElement* me_GEM_fail_Ch_region_GE1_NoL =
0094 iBooker.book2D("GEM_nFailingProbe_Ch_region_GE1_NoL", "GEM_nFailingProbe_Ch_region_GE1_NoL", 2, 0, 2, 36, 1, 37);
0095 MonitorElement* me_GE11_pass_Ch_ieta =
0096 iBooker.book2D("GE11_nPassingProbe_Ch_ieta", "GE11_nPassingProbe_Ch_ieta", 8, 1, 9, 36, 1, 37);
0097 MonitorElement* me_GE11_fail_Ch_ieta =
0098 iBooker.book2D("GE11_nFailingProbe_Ch_ieta", "GE11_nFailingProbe_Ch_ieta", 8, 1, 9, 36, 1, 37);
0099 MonitorElement* me_GE11_pass_Ch_phi = iBooker.book2D(
0100 "GE11_nPassingProbe_Ch_phi", "GE11_nPassingProbe_Ch_phi", 20, -TMath::Pi(), TMath::Pi(), 36, 1, 37);
0101 MonitorElement* me_GE11_fail_Ch_phi = iBooker.book2D(
0102 "GE11_nFailingProbe_Ch_phi", "GE11_nFailingProbe_Ch_phi", 20, -TMath::Pi(), TMath::Pi(), 36, 1, 37);
0103 MonitorElement* me_GE11_pass_allCh_1D =
0104 iBooker.book1D("GE11_nPassingProbe_allCh_1D", "GE11_nPassingProbe_allCh_1D", 2, -1.5, 1.5);
0105 MonitorElement* me_GE11_fail_allCh_1D =
0106 iBooker.book1D("GE11_nFailingProbe_allCh_1D", "GE11_nFailingProbe_allCh_1D", 2, -1.5, 1.5);
0107 MonitorElement* me_GE11_pass_chamber_1D =
0108 iBooker.book1D("GE11_nPassingProbe_chamber_1D", "GE11_nPassingProbe_chamber_1D", 36, 1, 37);
0109 MonitorElement* me_GE11_fail_chamber_1D =
0110 iBooker.book1D("GE11_nFailingProbe_chamber_1D", "GE11_nFailingProbe_chamber_1D", 36, 1, 37);
0111 MonitorElement* me_GE21_pass_Ch_ieta =
0112 iBooker.book2D("GE21_nPassingProbe_Ch_ieta", "GE21_nPassingProbe_Ch_ieta", 16, 1, 17, 18, 1, 19);
0113 MonitorElement* me_GE21_fail_Ch_ieta =
0114 iBooker.book2D("GE21_nFailingProbe_Ch_ieta", "GE21_nFailingProbe_Ch_ieta", 16, 1, 17, 18, 1, 19);
0115 MonitorElement* me_GE21_pass_Ch_phi = iBooker.book2D(
0116 "GE21_nPassingProbe_Ch_phi", "GE21_nPassingProbe_Ch_phi", 20, -TMath::Pi(), TMath::Pi(), 18, 1, 19);
0117 MonitorElement* me_GE21_fail_Ch_phi = iBooker.book2D(
0118 "GE21_nFailingProbe_Ch_phi", "GE21_nFailingProbe_Ch_phi", 20, -TMath::Pi(), TMath::Pi(), 18, 1, 19);
0119 MonitorElement* me_GE21_pass_allCh_1D =
0120 iBooker.book1D("GE21_nPassingProbe_allCh_1D", "GE21_nPassingProbe_allCh_1D", 2, -1.5, 1.5);
0121 MonitorElement* me_GE21_fail_allCh_1D =
0122 iBooker.book1D("GE21_nFailingProbe_allCh_1D", "GE21_nFailingProbe_allCh_1D", 2, -1.5, 1.5);
0123 MonitorElement* me_GE21_pass_chamber_1D =
0124 iBooker.book1D("GE21_nPassingProbe_chamber_1D", "GE21_nPassingProbe_chamber_1D", 18, 1, 19);
0125 MonitorElement* me_GE21_fail_chamber_1D =
0126 iBooker.book1D("GE21_nFailingProbe_chamber_1D", "GE21_nFailingProbe_chamber_1D", 18, 1, 19);
0127 MonitorElement* me_GEM_pass_chamber_p1_1D =
0128 iBooker.book1D("GEM_nPassingProbe_chamber_p1_1D", "GEM_nPassingProbe_chamber_p1_1D", 36, 1, 37);
0129 MonitorElement* me_GEM_fail_chamber_p1_1D =
0130 iBooker.book1D("GEM_nFailingProbe_chamber_p1_1D", "GEM_nFailingProbe_chamber_p1_1D", 36, 1, 37);
0131 MonitorElement* me_GEM_pass_chamber_p2_1D =
0132 iBooker.book1D("GEM_nPassingProbe_chamber_p2_1D", "GEM_nPassingProbe_chamber_p2_1D", 36, 1, 37);
0133 MonitorElement* me_GEM_fail_chamber_p2_1D =
0134 iBooker.book1D("GEM_nFailingProbe_chamber_p2_1D", "GEM_nFailingProbe_chamber_p2_1D", 36, 1, 37);
0135 MonitorElement* me_GEM_pass_chamber_n1_1D =
0136 iBooker.book1D("GEM_nPassingProbe_chamber_n1_1D", "GEM_nPassingProbe_chamber_n1_1D", 36, 1, 37);
0137 MonitorElement* me_GEM_fail_chamber_n1_1D =
0138 iBooker.book1D("GEM_nFailingProbe_chamber_n1_1D", "GEM_nFailingProbe_chamber_n1_1D", 36, 1, 37);
0139 MonitorElement* me_GEM_pass_chamber_n2_1D =
0140 iBooker.book1D("GEM_nPassingProbe_chamber_n2_1D", "GEM_nPassingProbe_chamber_n2_1D", 36, 1, 37);
0141 MonitorElement* me_GEM_fail_chamber_n2_1D =
0142 iBooker.book1D("GEM_nFailingProbe_chamber_n2_1D", "GEM_nFailingProbe_chamber_n2_1D", 36, 1, 37);
0143
0144 MonitorElement* me_GEM_pass_pt_1D = iBooker.book1D("GEM_nPassingProbe_pt_1D", "GEM_nPassingProbe_pt_1D", 20, 0, 100);
0145 MonitorElement* me_GEM_fail_pt_1D = iBooker.book1D("GEM_nFailingProbe_pt_1D", "GEM_nFailingProbe_pt_1D", 20, 0, 100);
0146 MonitorElement* me_GEM_pass_eta_1D =
0147 iBooker.book1D("GEM_nPassingProbe_eta_1D", "GEM_nPassingProbe_eta_1D", 24, 0, 2.4);
0148 MonitorElement* me_GEM_fail_eta_1D =
0149 iBooker.book1D("GEM_nFailingProbe_eta_1D", "GEM_nFailingProbe_eta_1D", 24, 0, 2.4);
0150 MonitorElement* me_GEM_pass_phi_1D =
0151 iBooker.book1D("GEM_nPassingProbe_phi_1D", "GEM_nPassingProbe_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0152 MonitorElement* me_GEM_fail_phi_1D =
0153 iBooker.book1D("GEM_nFailingProbe_phi_1D", "GEM_nFailingProbe_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0154
0155 MonitorElement* me_GEM_pass_pt_p1_1D =
0156 iBooker.book1D("GEM_nPassingProbe_pt_p1_1D", "GEM_nPassingProbe_pt_p1_1D", 20, 0, 100);
0157 MonitorElement* me_GEM_fail_pt_p1_1D =
0158 iBooker.book1D("GEM_nFailingProbe_pt_p1_1D", "GEM_nFailingProbe_pt_p1_1D", 20, 0, 100);
0159 MonitorElement* me_GEM_pass_eta_p1_1D =
0160 iBooker.book1D("GEM_nPassingProbe_eta_p1_1D", "GEM_nPassingProbe_eta_p1_1D", 24, 0, 2.4);
0161 MonitorElement* me_GEM_fail_eta_p1_1D =
0162 iBooker.book1D("GEM_nFailingProbe_eta_p1_1D", "GEM_nFailingProbe_eta_p1_1D", 24, 0, 2.4);
0163 MonitorElement* me_GEM_pass_phi_p1_1D =
0164 iBooker.book1D("GEM_nPassingProbe_phi_p1_1D", "GEM_nPassingProbe_phi_p1_1D", 20, -TMath::Pi(), TMath::Pi());
0165 MonitorElement* me_GEM_fail_phi_p1_1D =
0166 iBooker.book1D("GEM_nFailingProbe_phi_p1_1D", "GEM_nFailingProbe_phi_p1_1D", 20, -TMath::Pi(), TMath::Pi());
0167 MonitorElement* me_GEM_pass_pt_p2_1D =
0168 iBooker.book1D("GEM_nPassingProbe_pt_p2_1D", "GEM_nPassingProbe_pt_p2_1D", 20, 0, 100);
0169 MonitorElement* me_GEM_fail_pt_p2_1D =
0170 iBooker.book1D("GEM_nFailingProbe_pt_p2_1D", "GEM_nFailingProbe_pt_p2_1D", 20, 0, 100);
0171 MonitorElement* me_GEM_pass_eta_p2_1D =
0172 iBooker.book1D("GEM_nPassingProbe_eta_p2_1D", "GEM_nPassingProbe_eta_p2_1D", 24, 0, 2.4);
0173 MonitorElement* me_GEM_fail_eta_p2_1D =
0174 iBooker.book1D("GEM_nFailingProbe_eta_p2_1D", "GEM_nFailingProbe_eta_p2_1D", 24, 0, 2.4);
0175 MonitorElement* me_GEM_pass_phi_p2_1D =
0176 iBooker.book1D("GEM_nPassingProbe_phi_p2_1D", "GEM_nPassingProbe_phi_p2_1D", 20, -TMath::Pi(), TMath::Pi());
0177 MonitorElement* me_GEM_fail_phi_p2_1D =
0178 iBooker.book1D("GEM_nFailingProbe_phi_p2_1D", "GEM_nFailingProbe_phi_p2_1D", 20, -TMath::Pi(), TMath::Pi());
0179 MonitorElement* me_GEM_pass_pt_n1_1D =
0180 iBooker.book1D("GEM_nPassingProbe_pt_n1_1D", "GEM_nPassingProbe_pt_n1_1D", 20, 0, 100);
0181 MonitorElement* me_GEM_fail_pt_n1_1D =
0182 iBooker.book1D("GEM_nFailingProbe_pt_n1_1D", "GEM_nFailingProbe_pt_n1_1D", 20, 0, 100);
0183 MonitorElement* me_GEM_pass_eta_n1_1D =
0184 iBooker.book1D("GEM_nPassingProbe_eta_n1_1D", "GEM_nPassingProbe_eta_n1_1D", 24, 0, 2.4);
0185 MonitorElement* me_GEM_fail_eta_n1_1D =
0186 iBooker.book1D("GEM_nFailingProbe_eta_n1_1D", "GEM_nFailingProbe_eta_n1_1D", 24, 0, 2.4);
0187 MonitorElement* me_GEM_pass_phi_n1_1D =
0188 iBooker.book1D("GEM_nPassingProbe_phi_n1_1D", "GEM_nPassingProbe_phi_n1_1D", 20, -TMath::Pi(), TMath::Pi());
0189 MonitorElement* me_GEM_fail_phi_n1_1D =
0190 iBooker.book1D("GEM_nFailingProbe_phi_n1_1D", "GEM_nFailingProbe_phi_n1_1D", 20, -TMath::Pi(), TMath::Pi());
0191 MonitorElement* me_GEM_pass_pt_n2_1D =
0192 iBooker.book1D("GEM_nPassingProbe_pt_n2_1D", "GEM_nPassingProbe_pt_n2_1D", 20, 0, 100);
0193 MonitorElement* me_GEM_fail_pt_n2_1D =
0194 iBooker.book1D("GEM_nFailingProbe_pt_n2_1D", "GEM_nFailingProbe_pt_n2_1D", 20, 0, 100);
0195 MonitorElement* me_GEM_pass_eta_n2_1D =
0196 iBooker.book1D("GEM_nPassingProbe_eta_n2_1D", "GEM_nPassingProbe_eta_n2_1D", 24, 0, 2.4);
0197 MonitorElement* me_GEM_fail_eta_n2_1D =
0198 iBooker.book1D("GEM_nFailingProbe_eta_n2_1D", "GEM_nFailingProbe_eta_n2_1D", 24, 0, 2.4);
0199 MonitorElement* me_GEM_pass_phi_n2_1D =
0200 iBooker.book1D("GEM_nPassingProbe_phi_n2_1D", "GEM_nPassingProbe_phi_n2_1D", 20, -TMath::Pi(), TMath::Pi());
0201 MonitorElement* me_GEM_fail_phi_n2_1D =
0202 iBooker.book1D("GEM_nFailingProbe_phi_n2_1D", "GEM_nFailingProbe_phi_n2_1D", 20, -TMath::Pi(), TMath::Pi());
0203 MonitorElement* me_GE11_pass_Ch_ieta_p1 =
0204 iBooker.book2D("GE11_nPassingProbe_Ch_ieta_p1", "GE11_nPassingProbe_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0205 MonitorElement* me_GE11_fail_Ch_ieta_p1 =
0206 iBooker.book2D("GE11_nFailingProbe_Ch_ieta_p1", "GE11_nFailingProbe_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0207 MonitorElement* me_GE11_pass_Ch_ieta_p2 =
0208 iBooker.book2D("GE11_nPassingProbe_Ch_ieta_p2", "GE11_nPassingProbe_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0209 MonitorElement* me_GE11_fail_Ch_ieta_p2 =
0210 iBooker.book2D("GE11_nFailingProbe_Ch_ieta_p2", "GE11_nFailingProbe_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0211 MonitorElement* me_GE11_pass_Ch_ieta_n1 =
0212 iBooker.book2D("GE11_nPassingProbe_Ch_ieta_n1", "GE11_nPassingProbe_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0213 MonitorElement* me_GE11_fail_Ch_ieta_n1 =
0214 iBooker.book2D("GE11_nFailingProbe_Ch_ieta_n1", "GE11_nFailingProbe_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0215 MonitorElement* me_GE11_pass_Ch_ieta_n2 =
0216 iBooker.book2D("GE11_nPassingProbe_Ch_ieta_n2", "GE11_nPassingProbe_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0217 MonitorElement* me_GE11_fail_Ch_ieta_n2 =
0218 iBooker.book2D("GE11_nFailingProbe_Ch_ieta_n2", "GE11_nFailingProbe_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0219
0220
0221 MonitorElement* me_ME0_pass_chamber_1D =
0222 iBooker.book1D("ME0_nPassingProbe_chamber_1D", "ME0_nPassingProbe_chamber_1D", 18, 1, 19);
0223 MonitorElement* me_ME0_fail_chamber_1D =
0224 iBooker.book1D("ME0_nFailingProbe_chamber_1D", "ME0_nFailingProbe_chamber_1D", 18, 1, 19);
0225 MonitorElement* me_GEM_pass_Ch_region_layer_phase2 = iBooker.book2D(
0226 "GEM_nPassingProbe_Ch_region_layer_phase2", "GEM_nPassingProbe_Ch_region_layer_phase2", 10, 0, 10, 36, 1, 37);
0227 MonitorElement* me_GEM_fail_Ch_region_layer_phase2 = iBooker.book2D(
0228 "GEM_nFailingProbe_Ch_region_layer_phase2", "GEM_nFailingProbe_Ch_region_layer_phase2", 10, 0, 10, 36, 1, 37);
0229
0230 MonitorElement* OHmissing_me_GE11_pass_allCh_1D =
0231 iBooker.book1D("GE11_nPassingProbe_OHmissing_allCh_1D", "GE11_nPassingProbe_OHmissing_allCh_1D", 2, -1.5, 1.5);
0232 MonitorElement* OHmissing_me_GE11_fail_allCh_1D =
0233 iBooker.book1D("GE11_nFailingProbe_OHmissing_allCh_1D", "GE11_nFailingProbe_OHmissing_allCh_1D", 2, -1.5, 1.5);
0234 MonitorElement* OHmissing_me_GEM_pass_chamber_p1_1D = iBooker.book1D(
0235 "GEM_nPassingProbe_OHmissing_chamber_p1_1D", "GEM_nPassingProbe_OHmissing_chamber_p1_1D", 36, 1, 37);
0236 MonitorElement* OHmissing_me_GEM_fail_chamber_p1_1D = iBooker.book1D(
0237 "GEM_nFailingProbe_OHmissing_chamber_p1_1D", "GEM_nFailingProbe_OHmissing_chamber_p1_1D", 36, 1, 37);
0238 MonitorElement* OHmissing_me_GEM_pass_chamber_p2_1D = iBooker.book1D(
0239 "GEM_nPassingProbe_OHmissing_chamber_p2_1D", "GEM_nPassingProbe_OHmissing_chamber_p2_1D", 36, 1, 37);
0240 MonitorElement* OHmissing_me_GEM_fail_chamber_p2_1D = iBooker.book1D(
0241 "GEM_nFailingProbe_OHmissing_chamber_p2_1D", "GEM_nFailingProbe_OHmissing_chamber_p2_1D", 36, 1, 37);
0242 MonitorElement* OHmissing_me_GEM_pass_chamber_n1_1D = iBooker.book1D(
0243 "GEM_nPassingProbe_OHmissing_chamber_n1_1D", "GEM_nPassingProbe_OHmissing_chamber_n1_1D", 36, 1, 37);
0244 MonitorElement* OHmissing_me_GEM_fail_chamber_n1_1D = iBooker.book1D(
0245 "GEM_nFailingProbe_OHmissing_chamber_n1_1D", "GEM_nFailingProbe_OHmissing_chamber_n1_1D", 36, 1, 37);
0246 MonitorElement* OHmissing_me_GEM_pass_chamber_n2_1D = iBooker.book1D(
0247 "GEM_nPassingProbe_OHmissing_chamber_n2_1D", "GEM_nPassingProbe_OHmissing_chamber_n2_1D", 36, 1, 37);
0248 MonitorElement* OHmissing_me_GEM_fail_chamber_n2_1D = iBooker.book1D(
0249 "GEM_nFailingProbe_OHmissing_chamber_n2_1D", "GEM_nFailingProbe_OHmissing_chamber_n2_1D", 36, 1, 37);
0250 MonitorElement* OHmissing_me_GEM_pass_pt_1D =
0251 iBooker.book1D("GEM_nPassingProbe_OHmissing_pt_1D", "GEM_nPassingProbe_OHmissing_pt_1D", 20, 0, 100);
0252 MonitorElement* OHmissing_me_GEM_fail_pt_1D =
0253 iBooker.book1D("GEM_nFailingProbe_OHmissing_pt_1D", "GEM_nFailingProbe_OHmissing_pt_1D", 20, 0, 100);
0254 MonitorElement* OHmissing_me_GEM_pass_eta_1D =
0255 iBooker.book1D("GEM_nPassingProbe_OHmissing_eta_1D", "GEM_nPassingProbe_OHmissing_eta_1D", 24, 0, 2.4);
0256 MonitorElement* OHmissing_me_GEM_fail_eta_1D =
0257 iBooker.book1D("GEM_nFailingProbe_OHmissing_eta_1D", "GEM_nFailingProbe_OHmissing_eta_1D", 24, 0, 2.4);
0258 MonitorElement* OHmissing_me_GEM_pass_phi_1D = iBooker.book1D(
0259 "GEM_nPassingProbe_OHmissing_phi_1D", "GEM_nPassingProbe_OHmissing_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0260 MonitorElement* OHmissing_me_GEM_fail_phi_1D = iBooker.book1D(
0261 "GEM_nFailingProbe_OHmissing_phi_1D", "GEM_nFailingProbe_OHmissing_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0262 MonitorElement* OHmissing_me_GE11_pass_Ch_ieta_p1 = iBooker.book2D(
0263 "GE11_nPassingProbe_OHmissing_Ch_ieta_p1", "GE11_nPassingProbe_OHmissing_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0264 MonitorElement* OHmissing_me_GE11_fail_Ch_ieta_p1 = iBooker.book2D(
0265 "GE11_nFailingProbe_OHmissing_Ch_ieta_p1", "GE11_nFailingProbe_OHmissing_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0266 MonitorElement* OHmissing_me_GE11_pass_Ch_ieta_p2 = iBooker.book2D(
0267 "GE11_nPassingProbe_OHmissing_Ch_ieta_p2", "GE11_nPassingProbe_OHmissing_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0268 MonitorElement* OHmissing_me_GE11_fail_Ch_ieta_p2 = iBooker.book2D(
0269 "GE11_nFailingProbe_OHmissing_Ch_ieta_p2", "GE11_nFailingProbe_OHmissing_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0270 MonitorElement* OHmissing_me_GE11_pass_Ch_ieta_n1 = iBooker.book2D(
0271 "GE11_nPassingProbe_OHmissing_Ch_ieta_n1", "GE11_nPassingProbe_OHmissing_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0272 MonitorElement* OHmissing_me_GE11_fail_Ch_ieta_n1 = iBooker.book2D(
0273 "GE11_nFailingProbe_OHmissing_Ch_ieta_n1", "GE11_nFailingProbe_OHmissing_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0274 MonitorElement* OHmissing_me_GE11_pass_Ch_ieta_n2 = iBooker.book2D(
0275 "GE11_nPassingProbe_OHmissing_Ch_ieta_n2", "GE11_nPassingProbe_OHmissing_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0276 MonitorElement* OHmissing_me_GE11_fail_Ch_ieta_n2 = iBooker.book2D(
0277 "GE11_nFailingProbe_OHmissing_Ch_ieta_n2", "GE11_nFailingProbe_OHmissing_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0278
0279 OHmissing_me_GE11_pass_allCh_1D->setBinLabel(1, "GE-11", 1);
0280 OHmissing_me_GE11_pass_allCh_1D->setBinLabel(2, "GE+11", 1);
0281 OHmissing_me_GE11_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
0282
0283 OHmissing_me_GE11_fail_allCh_1D->setBinLabel(1, "GE-11", 1);
0284 OHmissing_me_GE11_fail_allCh_1D->setBinLabel(2, "GE+11", 1);
0285 OHmissing_me_GE11_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
0286
0287 OHmissing_me_GEM_pass_chamber_p1_1D->setAxisTitle("Chamber", 1);
0288 OHmissing_me_GEM_pass_chamber_p1_1D->setAxisTitle("Number of passing probes", 2);
0289 OHmissing_me_GEM_fail_chamber_p1_1D->setAxisTitle("Chamber", 1);
0290 OHmissing_me_GEM_fail_chamber_p1_1D->setAxisTitle("Number of failing probes", 2);
0291
0292 OHmissing_me_GEM_pass_chamber_p2_1D->setAxisTitle("Chamber", 1);
0293 OHmissing_me_GEM_pass_chamber_p2_1D->setAxisTitle("Number of passing probes", 2);
0294 OHmissing_me_GEM_fail_chamber_p2_1D->setAxisTitle("Chamber", 1);
0295 OHmissing_me_GEM_fail_chamber_p2_1D->setAxisTitle("Number of failing probes", 2);
0296
0297 OHmissing_me_GEM_pass_chamber_n1_1D->setAxisTitle("Chamber", 1);
0298 OHmissing_me_GEM_pass_chamber_n1_1D->setAxisTitle("Number of passing probes", 2);
0299 OHmissing_me_GEM_fail_chamber_n1_1D->setAxisTitle("Chamber", 1);
0300 OHmissing_me_GEM_fail_chamber_n1_1D->setAxisTitle("Number of failing probes", 2);
0301
0302 OHmissing_me_GEM_pass_chamber_n2_1D->setAxisTitle("Chamber", 1);
0303 OHmissing_me_GEM_pass_chamber_n2_1D->setAxisTitle("Number of passing probes", 2);
0304 OHmissing_me_GEM_fail_chamber_n2_1D->setAxisTitle("Chamber", 1);
0305 OHmissing_me_GEM_fail_chamber_n2_1D->setAxisTitle("Number of failing probes", 2);
0306
0307 OHmissing_me_GEM_pass_pt_1D->setAxisTitle("P_{T}", 1);
0308 OHmissing_me_GEM_pass_pt_1D->setAxisTitle("Number of passing probes", 2);
0309 OHmissing_me_GEM_fail_pt_1D->setAxisTitle("P_{T}", 1);
0310 OHmissing_me_GEM_fail_pt_1D->setAxisTitle("Number of failing probes", 2);
0311
0312 OHmissing_me_GEM_pass_eta_1D->setAxisTitle("#eta", 1);
0313 OHmissing_me_GEM_pass_eta_1D->setAxisTitle("Number of passing probes", 2);
0314 OHmissing_me_GEM_fail_eta_1D->setAxisTitle("#eta", 1);
0315 OHmissing_me_GEM_fail_eta_1D->setAxisTitle("Number of failing probes", 2);
0316
0317 OHmissing_me_GEM_pass_phi_1D->setAxisTitle("#phi", 1);
0318 OHmissing_me_GEM_pass_phi_1D->setAxisTitle("Number of passing probes", 2);
0319 OHmissing_me_GEM_fail_phi_1D->setAxisTitle("#phi", 1);
0320 OHmissing_me_GEM_fail_phi_1D->setAxisTitle("Number of failing probes", 2);
0321
0322 for (int i = 1; i < 37; ++i) {
0323 OHmissing_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0324 }
0325 for (int i = 1; i < 9; ++i) {
0326 OHmissing_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0327 }
0328 OHmissing_me_GE11_fail_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0329 OHmissing_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0330 OHmissing_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Number of failing probes", 3);
0331
0332 for (int i = 1; i < 37; ++i) {
0333 OHmissing_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0334 }
0335 for (int i = 1; i < 9; ++i) {
0336 OHmissing_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0337 }
0338 OHmissing_me_GE11_pass_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0339 OHmissing_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0340 OHmissing_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Number of passing probes", 3);
0341
0342 for (int i = 1; i < 37; ++i) {
0343 OHmissing_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0344 }
0345 for (int i = 1; i < 9; ++i) {
0346 OHmissing_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0347 }
0348 OHmissing_me_GE11_fail_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0349 OHmissing_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0350 OHmissing_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Number of failing probes", 3);
0351
0352 for (int i = 1; i < 37; ++i) {
0353 OHmissing_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0354 }
0355 for (int i = 1; i < 9; ++i) {
0356 OHmissing_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0357 }
0358 OHmissing_me_GE11_pass_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0359 OHmissing_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0360 OHmissing_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Number of passing probes", 3);
0361
0362 for (int i = 1; i < 37; ++i) {
0363 OHmissing_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0364 }
0365 for (int i = 1; i < 9; ++i) {
0366 OHmissing_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0367 }
0368 OHmissing_me_GE11_fail_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0369 OHmissing_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0370 OHmissing_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Number of failing probes", 3);
0371
0372 for (int i = 1; i < 37; ++i) {
0373 OHmissing_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0374 }
0375 for (int i = 1; i < 9; ++i) {
0376 OHmissing_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0377 }
0378 OHmissing_me_GE11_pass_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0379 OHmissing_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0380 OHmissing_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Number of passing probes", 3);
0381
0382 for (int i = 1; i < 37; ++i) {
0383 OHmissing_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0384 }
0385 for (int i = 1; i < 9; ++i) {
0386 OHmissing_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0387 }
0388 OHmissing_me_GE11_fail_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0389 OHmissing_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0390 OHmissing_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Number of failing probes", 3);
0391
0392 for (int i = 1; i < 37; ++i) {
0393 OHmissing_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0394 }
0395 for (int i = 1; i < 9; ++i) {
0396 OHmissing_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0397 }
0398 OHmissing_me_GE11_pass_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0399 OHmissing_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0400 OHmissing_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Number of passing probes", 3);
0401
0402
0403
0404 MonitorElement* OHerror_me_GE11_pass_allCh_1D =
0405 iBooker.book1D("GE11_nPassingProbe_OHerror_allCh_1D", "GE11_nPassingProbe_OHerror_allCh_1D", 2, -1.5, 1.5);
0406 MonitorElement* OHerror_me_GE11_fail_allCh_1D =
0407 iBooker.book1D("GE11_nFailingProbe_OHerror_allCh_1D", "GE11_nFailingProbe_OHerror_allCh_1D", 2, -1.5, 1.5);
0408 MonitorElement* OHerror_me_GEM_pass_chamber_p1_1D =
0409 iBooker.book1D("GEM_nPassingProbe_OHerror_chamber_p1_1D", "GEM_nPassingProbe_OHerror_chamber_p1_1D", 36, 1, 37);
0410 MonitorElement* OHerror_me_GEM_fail_chamber_p1_1D =
0411 iBooker.book1D("GEM_nFailingProbe_OHerror_chamber_p1_1D", "GEM_nFailingProbe_OHerror_chamber_p1_1D", 36, 1, 37);
0412 MonitorElement* OHerror_me_GEM_pass_chamber_p2_1D =
0413 iBooker.book1D("GEM_nPassingProbe_OHerror_chamber_p2_1D", "GEM_nPassingProbe_OHerror_chamber_p2_1D", 36, 1, 37);
0414 MonitorElement* OHerror_me_GEM_fail_chamber_p2_1D =
0415 iBooker.book1D("GEM_nFailingProbe_OHerror_chamber_p2_1D", "GEM_nFailingProbe_OHerror_chamber_p2_1D", 36, 1, 37);
0416 MonitorElement* OHerror_me_GEM_pass_chamber_n1_1D =
0417 iBooker.book1D("GEM_nPassingProbe_OHerror_chamber_n1_1D", "GEM_nPassingProbe_OHerror_chamber_n1_1D", 36, 1, 37);
0418 MonitorElement* OHerror_me_GEM_fail_chamber_n1_1D =
0419 iBooker.book1D("GEM_nFailingProbe_OHerror_chamber_n1_1D", "GEM_nFailingProbe_OHerror_chamber_n1_1D", 36, 1, 37);
0420 MonitorElement* OHerror_me_GEM_pass_chamber_n2_1D =
0421 iBooker.book1D("GEM_nPassingProbe_OHerror_chamber_n2_1D", "GEM_nPassingProbe_OHerror_chamber_n2_1D", 36, 1, 37);
0422 MonitorElement* OHerror_me_GEM_fail_chamber_n2_1D =
0423 iBooker.book1D("GEM_nFailingProbe_OHerror_chamber_n2_1D", "GEM_nFailingProbe_OHerror_chamber_n2_1D", 36, 1, 37);
0424 MonitorElement* OHerror_me_GEM_pass_pt_1D =
0425 iBooker.book1D("GEM_nPassingProbe_OHerror_pt_1D", "GEM_nPassingProbe_OHerror_pt_1D", 20, 0, 100);
0426 MonitorElement* OHerror_me_GEM_fail_pt_1D =
0427 iBooker.book1D("GEM_nFailingProbe_OHerror_pt_1D", "GEM_nFailingProbe_OHerror_pt_1D", 20, 0, 100);
0428 MonitorElement* OHerror_me_GEM_pass_eta_1D =
0429 iBooker.book1D("GEM_nPassingProbe_OHerror_eta_1D", "GEM_nPassingProbe_OHerror_eta_1D", 24, 0, 2.4);
0430 MonitorElement* OHerror_me_GEM_fail_eta_1D =
0431 iBooker.book1D("GEM_nFailingProbe_OHerror_eta_1D", "GEM_nFailingProbe_OHerror_eta_1D", 24, 0, 2.4);
0432 MonitorElement* OHerror_me_GEM_pass_phi_1D = iBooker.book1D(
0433 "GEM_nPassingProbe_OHerror_phi_1D", "GEM_nPassingProbe_OHerror_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0434 MonitorElement* OHerror_me_GEM_fail_phi_1D = iBooker.book1D(
0435 "GEM_nFailingProbe_OHerror_phi_1D", "GEM_nFailingProbe_OHerror_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0436 MonitorElement* OHerror_me_GE11_pass_Ch_ieta_p1 = iBooker.book2D(
0437 "GE11_nPassingProbe_OHerror_Ch_ieta_p1", "GE11_nPassingProbe_OHerror_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0438 MonitorElement* OHerror_me_GE11_fail_Ch_ieta_p1 = iBooker.book2D(
0439 "GE11_nFailingProbe_OHerror_Ch_ieta_p1", "GE11_nFailingProbe_OHerror_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0440 MonitorElement* OHerror_me_GE11_pass_Ch_ieta_p2 = iBooker.book2D(
0441 "GE11_nPassingProbe_OHerror_Ch_ieta_p2", "GE11_nPassingProbe_OHerror_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0442 MonitorElement* OHerror_me_GE11_fail_Ch_ieta_p2 = iBooker.book2D(
0443 "GE11_nFailingProbe_OHerror_Ch_ieta_p2", "GE11_nFailingProbe_OHerror_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0444 MonitorElement* OHerror_me_GE11_pass_Ch_ieta_n1 = iBooker.book2D(
0445 "GE11_nPassingProbe_OHerror_Ch_ieta_n1", "GE11_nPassingProbe_OHerror_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0446 MonitorElement* OHerror_me_GE11_fail_Ch_ieta_n1 = iBooker.book2D(
0447 "GE11_nFailingProbe_OHerror_Ch_ieta_n1", "GE11_nFailingProbe_OHerror_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0448 MonitorElement* OHerror_me_GE11_pass_Ch_ieta_n2 = iBooker.book2D(
0449 "GE11_nPassingProbe_OHerror_Ch_ieta_n2", "GE11_nPassingProbe_OHerror_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0450 MonitorElement* OHerror_me_GE11_fail_Ch_ieta_n2 = iBooker.book2D(
0451 "GE11_nFailingProbe_OHerror_Ch_ieta_n2", "GE11_nFailingProbe_OHerror_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0452
0453 OHerror_me_GE11_pass_allCh_1D->setBinLabel(1, "GE-11", 1);
0454 OHerror_me_GE11_pass_allCh_1D->setBinLabel(2, "GE+11", 1);
0455 OHerror_me_GE11_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
0456
0457 OHerror_me_GE11_fail_allCh_1D->setBinLabel(1, "GE-11", 1);
0458 OHerror_me_GE11_fail_allCh_1D->setBinLabel(2, "GE+11", 1);
0459 OHerror_me_GE11_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
0460
0461 OHerror_me_GEM_pass_chamber_p1_1D->setAxisTitle("Chamber", 1);
0462 OHerror_me_GEM_pass_chamber_p1_1D->setAxisTitle("Number of passing probes", 2);
0463 OHerror_me_GEM_fail_chamber_p1_1D->setAxisTitle("Chamber", 1);
0464 OHerror_me_GEM_fail_chamber_p1_1D->setAxisTitle("Number of failing probes", 2);
0465
0466 OHerror_me_GEM_pass_chamber_p2_1D->setAxisTitle("Chamber", 1);
0467 OHerror_me_GEM_pass_chamber_p2_1D->setAxisTitle("Number of passing probes", 2);
0468 OHerror_me_GEM_fail_chamber_p2_1D->setAxisTitle("Chamber", 1);
0469 OHerror_me_GEM_fail_chamber_p2_1D->setAxisTitle("Number of failing probes", 2);
0470
0471 OHerror_me_GEM_pass_chamber_n1_1D->setAxisTitle("Chamber", 1);
0472 OHerror_me_GEM_pass_chamber_n1_1D->setAxisTitle("Number of passing probes", 2);
0473 OHerror_me_GEM_fail_chamber_n1_1D->setAxisTitle("Chamber", 1);
0474 OHerror_me_GEM_fail_chamber_n1_1D->setAxisTitle("Number of failing probes", 2);
0475
0476 OHerror_me_GEM_pass_chamber_n2_1D->setAxisTitle("Chamber", 1);
0477 OHerror_me_GEM_pass_chamber_n2_1D->setAxisTitle("Number of passing probes", 2);
0478 OHerror_me_GEM_fail_chamber_n2_1D->setAxisTitle("Chamber", 1);
0479 OHerror_me_GEM_fail_chamber_n2_1D->setAxisTitle("Number of failing probes", 2);
0480
0481 OHerror_me_GEM_pass_pt_1D->setAxisTitle("P_{T}", 1);
0482 OHerror_me_GEM_pass_pt_1D->setAxisTitle("Number of passing probes", 2);
0483 OHerror_me_GEM_fail_pt_1D->setAxisTitle("P_{T}", 1);
0484 OHerror_me_GEM_fail_pt_1D->setAxisTitle("Number of failing probes", 2);
0485
0486 OHerror_me_GEM_pass_eta_1D->setAxisTitle("#eta", 1);
0487 OHerror_me_GEM_pass_eta_1D->setAxisTitle("Number of passing probes", 2);
0488 OHerror_me_GEM_fail_eta_1D->setAxisTitle("#eta", 1);
0489 OHerror_me_GEM_fail_eta_1D->setAxisTitle("Number of failing probes", 2);
0490
0491 OHerror_me_GEM_pass_phi_1D->setAxisTitle("#phi", 1);
0492 OHerror_me_GEM_pass_phi_1D->setAxisTitle("Number of passing probes", 2);
0493 OHerror_me_GEM_fail_phi_1D->setAxisTitle("#phi", 1);
0494 OHerror_me_GEM_fail_phi_1D->setAxisTitle("Number of failing probes", 2);
0495 for (int i = 1; i < 37; ++i) {
0496 OHerror_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0497 }
0498 for (int i = 1; i < 9; ++i) {
0499 OHerror_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0500 }
0501 OHerror_me_GE11_fail_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0502 OHerror_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0503 OHerror_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Number of failing probes", 3);
0504
0505 for (int i = 1; i < 37; ++i) {
0506 OHerror_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0507 }
0508 for (int i = 1; i < 9; ++i) {
0509 OHerror_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0510 }
0511 OHerror_me_GE11_pass_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0512 OHerror_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0513 OHerror_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Number of passing probes", 3);
0514
0515 for (int i = 1; i < 37; ++i) {
0516 OHerror_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0517 }
0518 for (int i = 1; i < 9; ++i) {
0519 OHerror_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0520 }
0521 OHerror_me_GE11_fail_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0522 OHerror_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0523 OHerror_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Number of failing probes", 3);
0524
0525 for (int i = 1; i < 37; ++i) {
0526 OHerror_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0527 }
0528 for (int i = 1; i < 9; ++i) {
0529 OHerror_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0530 }
0531 OHerror_me_GE11_pass_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0532 OHerror_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0533 OHerror_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Number of passing probes", 3);
0534
0535 for (int i = 1; i < 37; ++i) {
0536 OHerror_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0537 }
0538 for (int i = 1; i < 9; ++i) {
0539 OHerror_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0540 }
0541 OHerror_me_GE11_fail_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0542 OHerror_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0543 OHerror_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Number of failing probes", 3);
0544
0545 for (int i = 1; i < 37; ++i) {
0546 OHerror_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0547 }
0548 for (int i = 1; i < 9; ++i) {
0549 OHerror_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0550 }
0551 OHerror_me_GE11_pass_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0552 OHerror_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0553 OHerror_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Number of passing probes", 3);
0554
0555 for (int i = 1; i < 37; ++i) {
0556 OHerror_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0557 }
0558 for (int i = 1; i < 9; ++i) {
0559 OHerror_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0560 }
0561 OHerror_me_GE11_fail_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0562 OHerror_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0563 OHerror_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Number of failing probes", 3);
0564
0565 for (int i = 1; i < 37; ++i) {
0566 OHerror_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0567 }
0568 for (int i = 1; i < 9; ++i) {
0569 OHerror_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0570 }
0571 OHerror_me_GE11_pass_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0572 OHerror_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0573 OHerror_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Number of passing probes", 3);
0574
0575
0576
0577 MonitorElement* VFATMask_me_GE11_pass_allCh_1D =
0578 iBooker.book1D("GE11_nPassingProbe_VFATMask_allCh_1D", "GE11_nPassingProbe_VFATMask_allCh_1D", 2, -1.5, 1.5);
0579 MonitorElement* VFATMask_me_GE11_fail_allCh_1D =
0580 iBooker.book1D("GE11_nFailingProbe_VFATMask_allCh_1D", "GE11_nFailingProbe_VFATMask_allCh_1D", 2, -1.5, 1.5);
0581 MonitorElement* VFATMask_me_GEM_pass_chamber_p1_1D =
0582 iBooker.book1D("GEM_nPassingProbe_VFATMask_chamber_p1_1D", "GEM_nPassingProbe_VFATMask_chamber_p1_1D", 36, 1, 37);
0583 MonitorElement* VFATMask_me_GEM_fail_chamber_p1_1D =
0584 iBooker.book1D("GEM_nFailingProbe_VFATMask_chamber_p1_1D", "GEM_nFailingProbe_VFATMask_chamber_p1_1D", 36, 1, 37);
0585 MonitorElement* VFATMask_me_GEM_pass_chamber_p2_1D =
0586 iBooker.book1D("GEM_nPassingProbe_VFATMask_chamber_p2_1D", "GEM_nPassingProbe_VFATMask_chamber_p2_1D", 36, 1, 37);
0587 MonitorElement* VFATMask_me_GEM_fail_chamber_p2_1D =
0588 iBooker.book1D("GEM_nFailingProbe_VFATMask_chamber_p2_1D", "GEM_nFailingProbe_VFATMask_chamber_p2_1D", 36, 1, 37);
0589 MonitorElement* VFATMask_me_GEM_pass_chamber_n1_1D =
0590 iBooker.book1D("GEM_nPassingProbe_VFATMask_chamber_n1_1D", "GEM_nPassingProbe_VFATMask_chamber_n1_1D", 36, 1, 37);
0591 MonitorElement* VFATMask_me_GEM_fail_chamber_n1_1D =
0592 iBooker.book1D("GEM_nFailingProbe_VFATMask_chamber_n1_1D", "GEM_nFailingProbe_VFATMask_chamber_n1_1D", 36, 1, 37);
0593 MonitorElement* VFATMask_me_GEM_pass_chamber_n2_1D =
0594 iBooker.book1D("GEM_nPassingProbe_VFATMask_chamber_n2_1D", "GEM_nPassingProbe_VFATMask_chamber_n2_1D", 36, 1, 37);
0595 MonitorElement* VFATMask_me_GEM_fail_chamber_n2_1D =
0596 iBooker.book1D("GEM_nFailingProbe_VFATMask_chamber_n2_1D", "GEM_nFailingProbe_VFATMask_chamber_n2_1D", 36, 1, 37);
0597 MonitorElement* VFATMask_me_GEM_pass_pt_1D =
0598 iBooker.book1D("GEM_nPassingProbe_VFATMask_pt_1D", "GEM_nPassingProbe_VFATMask_pt_1D", 20, 0, 100);
0599 MonitorElement* VFATMask_me_GEM_fail_pt_1D =
0600 iBooker.book1D("GEM_nFailingProbe_VFATMask_pt_1D", "GEM_nFailingProbe_VFATMask_pt_1D", 20, 0, 100);
0601 MonitorElement* VFATMask_me_GEM_pass_eta_1D =
0602 iBooker.book1D("GEM_nPassingProbe_VFATMask_eta_1D", "GEM_nPassingProbe_VFATMask_eta_1D", 24, 0, 2.4);
0603 MonitorElement* VFATMask_me_GEM_fail_eta_1D =
0604 iBooker.book1D("GEM_nFailingProbe_VFATMask_eta_1D", "GEM_nFailingProbe_VFATMask_eta_1D", 24, 0, 2.4);
0605 MonitorElement* VFATMask_me_GEM_pass_phi_1D = iBooker.book1D(
0606 "GEM_nPassingProbe_VFATMask_phi_1D", "GEM_nPassingProbe_VFATMask_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0607 MonitorElement* VFATMask_me_GEM_fail_phi_1D = iBooker.book1D(
0608 "GEM_nFailingProbe_VFATMask_phi_1D", "GEM_nFailingProbe_VFATMask_phi_1D", 20, -TMath::Pi(), TMath::Pi());
0609 MonitorElement* VFATMask_me_GE11_pass_Ch_ieta_p1 = iBooker.book2D(
0610 "GE11_nPassingProbe_VFATMask_Ch_ieta_p1", "GE11_nPassingProbe_VFATMask_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0611 MonitorElement* VFATMask_me_GE11_fail_Ch_ieta_p1 = iBooker.book2D(
0612 "GE11_nFailingProbe_VFATMask_Ch_ieta_p1", "GE11_nFailingProbe_VFATMask_Ch_ieta_p1", 8, 1, 9, 36, 1, 37);
0613 MonitorElement* VFATMask_me_GE11_pass_Ch_ieta_p2 = iBooker.book2D(
0614 "GE11_nPassingProbe_VFATMask_Ch_ieta_p2", "GE11_nPassingProbe_VFATMask_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0615 MonitorElement* VFATMask_me_GE11_fail_Ch_ieta_p2 = iBooker.book2D(
0616 "GE11_nFailingProbe_VFATMask_Ch_ieta_p2", "GE11_nFailingProbe_VFATMask_Ch_ieta_p2", 8, 1, 9, 36, 1, 37);
0617 MonitorElement* VFATMask_me_GE11_pass_Ch_ieta_n1 = iBooker.book2D(
0618 "GE11_nPassingProbe_VFATMask_Ch_ieta_n1", "GE11_nPassingProbe_VFATMask_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0619 MonitorElement* VFATMask_me_GE11_fail_Ch_ieta_n1 = iBooker.book2D(
0620 "GE11_nFailingProbe_VFATMask_Ch_ieta_n1", "GE11_nFailingProbe_VFATMask_Ch_ieta_n1", 8, 1, 9, 36, 1, 37);
0621 MonitorElement* VFATMask_me_GE11_pass_Ch_ieta_n2 = iBooker.book2D(
0622 "GE11_nPassingProbe_VFATMask_Ch_ieta_n2", "GE11_nPassingProbe_VFATMask_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0623 MonitorElement* VFATMask_me_GE11_fail_Ch_ieta_n2 = iBooker.book2D(
0624 "GE11_nFailingProbe_VFATMask_Ch_ieta_n2", "GE11_nFailingProbe_VFATMask_Ch_ieta_n2", 8, 1, 9, 36, 1, 37);
0625
0626 VFATMask_me_GE11_pass_allCh_1D->setBinLabel(1, "GE-11", 1);
0627 VFATMask_me_GE11_pass_allCh_1D->setBinLabel(2, "GE+11", 1);
0628 VFATMask_me_GE11_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
0629
0630 VFATMask_me_GE11_fail_allCh_1D->setBinLabel(1, "GE-11", 1);
0631 VFATMask_me_GE11_fail_allCh_1D->setBinLabel(2, "GE+11", 1);
0632 VFATMask_me_GE11_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
0633
0634 VFATMask_me_GEM_pass_chamber_p1_1D->setAxisTitle("Chamber", 1);
0635 VFATMask_me_GEM_pass_chamber_p1_1D->setAxisTitle("Number of passing probes", 2);
0636 VFATMask_me_GEM_fail_chamber_p1_1D->setAxisTitle("Chamber", 1);
0637 VFATMask_me_GEM_fail_chamber_p1_1D->setAxisTitle("Number of failing probes", 2);
0638
0639 VFATMask_me_GEM_pass_chamber_p2_1D->setAxisTitle("Chamber", 1);
0640 VFATMask_me_GEM_pass_chamber_p2_1D->setAxisTitle("Number of passing probes", 2);
0641 VFATMask_me_GEM_fail_chamber_p2_1D->setAxisTitle("Chamber", 1);
0642 VFATMask_me_GEM_fail_chamber_p2_1D->setAxisTitle("Number of failing probes", 2);
0643
0644 VFATMask_me_GEM_pass_chamber_n1_1D->setAxisTitle("Chamber", 1);
0645 VFATMask_me_GEM_pass_chamber_n1_1D->setAxisTitle("Number of passing probes", 2);
0646 VFATMask_me_GEM_fail_chamber_n1_1D->setAxisTitle("Chamber", 1);
0647 VFATMask_me_GEM_fail_chamber_n1_1D->setAxisTitle("Number of failing probes", 2);
0648
0649 VFATMask_me_GEM_pass_chamber_n2_1D->setAxisTitle("Chamber", 1);
0650 VFATMask_me_GEM_pass_chamber_n2_1D->setAxisTitle("Number of passing probes", 2);
0651 VFATMask_me_GEM_fail_chamber_n2_1D->setAxisTitle("Chamber", 1);
0652 VFATMask_me_GEM_fail_chamber_n2_1D->setAxisTitle("Number of failing probes", 2);
0653
0654 VFATMask_me_GEM_pass_pt_1D->setAxisTitle("P_{T}", 1);
0655 VFATMask_me_GEM_pass_pt_1D->setAxisTitle("Number of passing probes", 2);
0656 VFATMask_me_GEM_fail_pt_1D->setAxisTitle("P_{T}", 1);
0657 VFATMask_me_GEM_fail_pt_1D->setAxisTitle("Number of failing probes", 2);
0658
0659 VFATMask_me_GEM_pass_eta_1D->setAxisTitle("#eta", 1);
0660 VFATMask_me_GEM_pass_eta_1D->setAxisTitle("Number of passing probes", 2);
0661 VFATMask_me_GEM_fail_eta_1D->setAxisTitle("#eta", 1);
0662 VFATMask_me_GEM_fail_eta_1D->setAxisTitle("Number of failing probes", 2);
0663
0664 VFATMask_me_GEM_pass_phi_1D->setAxisTitle("#phi", 1);
0665 VFATMask_me_GEM_pass_phi_1D->setAxisTitle("Number of passing probes", 2);
0666 VFATMask_me_GEM_fail_phi_1D->setAxisTitle("#phi", 1);
0667 VFATMask_me_GEM_fail_phi_1D->setAxisTitle("Number of failing probes", 2);
0668
0669 for (int i = 1; i < 37; ++i) {
0670 VFATMask_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0671 }
0672 for (int i = 1; i < 9; ++i) {
0673 VFATMask_me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0674 }
0675 VFATMask_me_GE11_fail_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0676 VFATMask_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0677 VFATMask_me_GE11_fail_Ch_ieta_p1->setAxisTitle("Number of failing probes", 3);
0678
0679 for (int i = 1; i < 37; ++i) {
0680 VFATMask_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
0681 }
0682 for (int i = 1; i < 9; ++i) {
0683 VFATMask_me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
0684 }
0685 VFATMask_me_GE11_pass_Ch_ieta_p1->setAxisTitle("#ieta", 1);
0686 VFATMask_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Chamber", 2);
0687 VFATMask_me_GE11_pass_Ch_ieta_p1->setAxisTitle("Number of passing probes", 3);
0688
0689 for (int i = 1; i < 37; ++i) {
0690 VFATMask_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0691 }
0692 for (int i = 1; i < 9; ++i) {
0693 VFATMask_me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0694 }
0695 VFATMask_me_GE11_fail_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0696 VFATMask_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0697 VFATMask_me_GE11_fail_Ch_ieta_p2->setAxisTitle("Number of failing probes", 3);
0698
0699 for (int i = 1; i < 37; ++i) {
0700 VFATMask_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
0701 }
0702 for (int i = 1; i < 9; ++i) {
0703 VFATMask_me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
0704 }
0705 VFATMask_me_GE11_pass_Ch_ieta_p2->setAxisTitle("#ieta", 1);
0706 VFATMask_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Chamber", 2);
0707 VFATMask_me_GE11_pass_Ch_ieta_p2->setAxisTitle("Number of passing probes", 3);
0708
0709 for (int i = 1; i < 37; ++i) {
0710 VFATMask_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0711 }
0712 for (int i = 1; i < 9; ++i) {
0713 VFATMask_me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0714 }
0715 VFATMask_me_GE11_fail_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0716 VFATMask_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0717 VFATMask_me_GE11_fail_Ch_ieta_n1->setAxisTitle("Number of failing probes", 3);
0718
0719 for (int i = 1; i < 37; ++i) {
0720 VFATMask_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
0721 }
0722 for (int i = 1; i < 9; ++i) {
0723 VFATMask_me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
0724 }
0725 VFATMask_me_GE11_pass_Ch_ieta_n1->setAxisTitle("#ieta", 1);
0726 VFATMask_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Chamber", 2);
0727 VFATMask_me_GE11_pass_Ch_ieta_n1->setAxisTitle("Number of passing probes", 3);
0728
0729 for (int i = 1; i < 37; ++i) {
0730 VFATMask_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0731 }
0732 for (int i = 1; i < 9; ++i) {
0733 VFATMask_me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0734 }
0735 VFATMask_me_GE11_fail_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0736 VFATMask_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0737 VFATMask_me_GE11_fail_Ch_ieta_n2->setAxisTitle("Number of failing probes", 3);
0738
0739 for (int i = 1; i < 37; ++i) {
0740 VFATMask_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
0741 }
0742 for (int i = 1; i < 9; ++i) {
0743 VFATMask_me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
0744 }
0745 VFATMask_me_GE11_pass_Ch_ieta_n2->setAxisTitle("#ieta", 1);
0746 VFATMask_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Chamber", 2);
0747 VFATMask_me_GE11_pass_Ch_ieta_n2->setAxisTitle("Number of passing probes", 3);
0748
0749 me_GE11_pass_allCh_1D->setBinLabel(1, "GE-11", 1);
0750 me_GE11_pass_allCh_1D->setBinLabel(2, "GE+11", 1);
0751 me_GE11_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
0752
0753 me_GE11_fail_allCh_1D->setBinLabel(1, "GE-11", 1);
0754 me_GE11_fail_allCh_1D->setBinLabel(2, "GE+11", 1);
0755 me_GE11_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
0756
0757 me_GE11_pass_chamber_1D->setAxisTitle("Chamber", 1);
0758 me_GE11_pass_chamber_1D->setAxisTitle("Number of passing probes", 2);
0759 me_GE11_fail_chamber_1D->setAxisTitle("Chamber", 1);
0760 me_GE11_fail_chamber_1D->setAxisTitle("Number of failing probes", 2);
0761
0762 me_GE21_pass_allCh_1D->setBinLabel(1, "GE-21", 1);
0763 me_GE21_pass_allCh_1D->setBinLabel(2, "GE+21", 1);
0764 me_GE21_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
0765
0766 me_GE21_fail_allCh_1D->setBinLabel(1, "GE-21", 1);
0767 me_GE21_fail_allCh_1D->setBinLabel(2, "GE+21", 1);
0768 me_GE21_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
0769
0770 me_GE21_pass_chamber_1D->setAxisTitle("Chamber", 1);
0771 me_GE21_pass_chamber_1D->setAxisTitle("Number of passing probes", 2);
0772 me_GE21_fail_chamber_1D->setAxisTitle("Chamber", 1);
0773 me_GE21_fail_chamber_1D->setAxisTitle("Number of failing probes", 2);
0774
0775 me_GEM_pass_chamber_p1_1D->setAxisTitle("Chamber", 1);
0776 me_GEM_pass_chamber_p1_1D->setAxisTitle("Number of passing probes", 2);
0777 me_GEM_fail_chamber_p1_1D->setAxisTitle("Chamber", 1);
0778 me_GEM_fail_chamber_p1_1D->setAxisTitle("Number of failing probes", 2);
0779
0780 me_GEM_pass_chamber_p2_1D->setAxisTitle("Chamber", 1);
0781 me_GEM_pass_chamber_p2_1D->setAxisTitle("Number of passing probes", 2);
0782 me_GEM_fail_chamber_p2_1D->setAxisTitle("Chamber", 1);
0783 me_GEM_fail_chamber_p2_1D->setAxisTitle("Number of failing probes", 2);
0784
0785 me_GEM_pass_chamber_n1_1D->setAxisTitle("Chamber", 1);
0786 me_GEM_pass_chamber_n1_1D->setAxisTitle("Number of passing probes", 2);
0787 me_GEM_fail_chamber_n1_1D->setAxisTitle("Chamber", 1);
0788 me_GEM_fail_chamber_n1_1D->setAxisTitle("Number of failing probes", 2);
0789
0790 me_GEM_pass_chamber_n2_1D->setAxisTitle("Chamber", 1);
0791 me_GEM_pass_chamber_n2_1D->setAxisTitle("Number of passing probes", 2);
0792 me_GEM_fail_chamber_n2_1D->setAxisTitle("Chamber", 1);
0793 me_GEM_fail_chamber_n2_1D->setAxisTitle("Number of failing probes", 2);
0794
0795 me_GEM_pass_pt_1D->setAxisTitle("P_{T}", 1);
0796 me_GEM_pass_pt_1D->setAxisTitle("Number of passing probes", 2);
0797 me_GEM_fail_pt_1D->setAxisTitle("P_{T}", 1);
0798 me_GEM_fail_pt_1D->setAxisTitle("Number of failing probes", 2);
0799
0800 me_GEM_pass_eta_1D->setAxisTitle("#eta", 1);
0801 me_GEM_pass_eta_1D->setAxisTitle("Number of passing probes", 2);
0802 me_GEM_fail_eta_1D->setAxisTitle("#eta", 1);
0803 me_GEM_fail_eta_1D->setAxisTitle("Number of failing probes", 2);
0804
0805 me_GEM_pass_phi_1D->setAxisTitle("#phi", 1);
0806 me_GEM_pass_phi_1D->setAxisTitle("Number of passing probes", 2);
0807 me_GEM_fail_phi_1D->setAxisTitle("#phi", 1);
0808 me_GEM_fail_phi_1D->setAxisTitle("Number of failing probes", 2);
0809
0810 me_GEM_pass_pt_p1_1D->setAxisTitle("P_{T}", 1);
0811 me_GEM_pass_pt_p1_1D->setAxisTitle("Number of passing probes", 2);
0812 me_GEM_fail_pt_p1_1D->setAxisTitle("P_{T}", 1);
0813 me_GEM_fail_pt_p1_1D->setAxisTitle("Number of failing probes", 2);
0814
0815 me_GEM_pass_eta_p1_1D->setAxisTitle("#eta", 1);
0816 me_GEM_pass_eta_p1_1D->setAxisTitle("Number of passing probes", 2);
0817 me_GEM_fail_eta_p1_1D->setAxisTitle("#eta", 1);
0818 me_GEM_fail_eta_p1_1D->setAxisTitle("Number of failing probes", 2);
0819
0820 me_GEM_pass_phi_p1_1D->setAxisTitle("#phi", 1);
0821 me_GEM_pass_phi_p1_1D->setAxisTitle("Number of passing probes", 2);
0822 me_GEM_fail_phi_p1_1D->setAxisTitle("#phi", 1);
0823 me_GEM_fail_phi_p1_1D->setAxisTitle("Number of failing probes", 2);
0824
0825 me_GEM_pass_pt_p2_1D->setAxisTitle("P_{T}", 1);
0826 me_GEM_pass_pt_p2_1D->setAxisTitle("Number of passing probes", 2);
0827 me_GEM_fail_pt_p2_1D->setAxisTitle("P_{T}", 1);
0828 me_GEM_fail_pt_p2_1D->setAxisTitle("Number of failing probes", 2);
0829
0830 me_GEM_pass_eta_p2_1D->setAxisTitle("#eta", 1);
0831 me_GEM_pass_eta_p2_1D->setAxisTitle("Number of passing probes", 2);
0832 me_GEM_fail_eta_p2_1D->setAxisTitle("#eta", 1);
0833 me_GEM_fail_eta_p2_1D->setAxisTitle("Number of failing probes", 2);
0834
0835 me_GEM_pass_phi_p2_1D->setAxisTitle("#phi", 1);
0836 me_GEM_pass_phi_p2_1D->setAxisTitle("Number of passing probes", 2);
0837 me_GEM_fail_phi_p2_1D->setAxisTitle("#phi", 1);
0838 me_GEM_fail_phi_p2_1D->setAxisTitle("Number of failing probes", 2);
0839
0840 me_GEM_pass_pt_n1_1D->setAxisTitle("P_{T}", 1);
0841 me_GEM_pass_pt_n1_1D->setAxisTitle("Number of passing probes", 2);
0842 me_GEM_fail_pt_n1_1D->setAxisTitle("P_{T}", 1);
0843 me_GEM_fail_pt_n1_1D->setAxisTitle("Number of failing probes", 2);
0844
0845 me_GEM_pass_eta_n1_1D->setAxisTitle("#eta", 1);
0846 me_GEM_pass_eta_n1_1D->setAxisTitle("Number of passing probes", 2);
0847 me_GEM_fail_eta_n1_1D->setAxisTitle("#eta", 1);
0848 me_GEM_fail_eta_n1_1D->setAxisTitle("Number of failing probes", 2);
0849
0850 me_GEM_pass_phi_n1_1D->setAxisTitle("#phi", 1);
0851 me_GEM_pass_phi_n1_1D->setAxisTitle("Number of passing probes", 2);
0852 me_GEM_fail_phi_n1_1D->setAxisTitle("#phi", 1);
0853 me_GEM_fail_phi_n1_1D->setAxisTitle("Number of failing probes", 2);
0854
0855 me_GEM_pass_pt_n2_1D->setAxisTitle("P_{T}", 1);
0856 me_GEM_pass_pt_n2_1D->setAxisTitle("Number of passing probes", 2);
0857 me_GEM_fail_pt_n2_1D->setAxisTitle("P_{T}", 1);
0858 me_GEM_fail_pt_n2_1D->setAxisTitle("Number of failing probes", 2);
0859
0860 me_GEM_pass_eta_n2_1D->setAxisTitle("#eta", 1);
0861 me_GEM_pass_eta_n2_1D->setAxisTitle("Number of passing probes", 2);
0862 me_GEM_fail_eta_n2_1D->setAxisTitle("#eta", 1);
0863 me_GEM_fail_eta_n2_1D->setAxisTitle("Number of failing probes", 2);
0864
0865 me_GEM_pass_phi_n2_1D->setAxisTitle("#phi", 1);
0866 me_GEM_pass_phi_n2_1D->setAxisTitle("Number of passing probes", 2);
0867 me_GEM_fail_phi_n2_1D->setAxisTitle("#phi", 1);
0868 me_GEM_fail_phi_n2_1D->setAxisTitle("Number of failing probes", 2);
0869
0870 me_GE11_fail_Ch_region->setBinLabel(1, "GE-11", 1);
0871 me_GE11_fail_Ch_region->setBinLabel(2, "GE+11", 1);
0872 for (int i = 1; i < 37; ++i) {
0873 me_GE11_fail_Ch_region->setBinLabel(i, std::to_string(i), 2);
0874 }
0875 me_GE11_fail_Ch_region->setAxisTitle("Chamber", 2);
0876 me_GE11_fail_Ch_region->setAxisTitle("Number of failing probes", 3);
0877
0878 me_GE11_pass_Ch_region->setBinLabel(1, "GE-11", 1);
0879 me_GE11_pass_Ch_region->setBinLabel(2, "GE+11", 1);
0880 for (int i = 1; i < 37; ++i) {
0881 me_GE11_pass_Ch_region->setBinLabel(i, std::to_string(i), 2);
0882 }
0883 me_GE11_pass_Ch_region->setAxisTitle("Chamber", 2);
0884 me_GE11_pass_Ch_region->setAxisTitle("Number of passing probes", 3);
0885
0886 me_GE21_fail_Ch_region->setBinLabel(1, "GE-21", 1);
0887 me_GE21_fail_Ch_region->setBinLabel(2, "GE+21", 1);
0888 for (int i = 1; i < 19; ++i) {
0889 me_GE21_fail_Ch_region->setBinLabel(i, std::to_string(i), 2);
0890 }
0891 me_GE21_fail_Ch_region->setAxisTitle("Chamber", 2);
0892 me_GE21_fail_Ch_region->setAxisTitle("Number of failing probes", 3);
0893
0894 me_GE21_pass_Ch_region->setBinLabel(1, "GE-21", 1);
0895 me_GE21_pass_Ch_region->setBinLabel(2, "GE+21", 1);
0896 for (int i = 1; i < 19; ++i) {
0897 me_GE21_pass_Ch_region->setBinLabel(i, std::to_string(i), 2);
0898 }
0899 me_GE21_pass_Ch_region->setAxisTitle("Chamber", 2);
0900 me_GE21_pass_Ch_region->setAxisTitle("Number of passing probes", 3);
0901
0902 me_GEM_fail_Ch_region_GE1->setBinLabel(1, "GE-1/1_L2", 1);
0903 me_GEM_fail_Ch_region_GE1->setBinLabel(2, "GE-1/1_L1", 1);
0904 me_GEM_fail_Ch_region_GE1->setBinLabel(3, "GE+1/1_L1", 1);
0905 me_GEM_fail_Ch_region_GE1->setBinLabel(4, "GE+1/1_L2", 1);
0906 for (int i = 1; i < 37; ++i) {
0907 me_GEM_fail_Ch_region_GE1->setBinLabel(i, std::to_string(i), 2);
0908 }
0909 me_GEM_fail_Ch_region_GE1->setAxisTitle("Chamber", 2);
0910 me_GEM_fail_Ch_region_GE1->setAxisTitle("Number of passing probes", 3);
0911
0912 me_GEM_pass_Ch_region_GE1->setBinLabel(1, "GE-1/1_L2", 1);
0913 me_GEM_pass_Ch_region_GE1->setBinLabel(2, "GE-1/1_L1", 1);
0914 me_GEM_pass_Ch_region_GE1->setBinLabel(3, "GE+1/1_L1", 1);
0915 me_GEM_pass_Ch_region_GE1->setBinLabel(4, "GE+1/1_L2", 1);
0916 for (int i = 1; i < 37; ++i) {
0917 me_GEM_pass_Ch_region_GE1->setBinLabel(i, std::to_string(i), 2);
0918 }
0919 me_GEM_pass_Ch_region_GE1->setAxisTitle("Chamber", 2);
0920 me_GEM_pass_Ch_region_GE1->setAxisTitle("Number of passing probes", 3);
0921
0922 me_GEM_fail_Ch_region_GE1_NoL->setBinLabel(1, "GE-1", 1);
0923 me_GEM_fail_Ch_region_GE1_NoL->setBinLabel(2, "GE+1", 1);
0924 for (int i = 1; i < 37; ++i) {
0925 me_GEM_fail_Ch_region_GE1_NoL->setBinLabel(i, std::to_string(i), 2);
0926 }
0927 me_GEM_fail_Ch_region_GE1_NoL->setAxisTitle("Chamber", 2);
0928 me_GEM_fail_Ch_region_GE1_NoL->setAxisTitle("Number of passing probes", 3);
0929
0930 me_GEM_pass_Ch_region_GE1_NoL->setBinLabel(1, "GE-1", 1);
0931 me_GEM_pass_Ch_region_GE1_NoL->setBinLabel(2, "GE+1", 1);
0932 for (int i = 1; i < 37; ++i) {
0933 me_GEM_pass_Ch_region_GE1_NoL->setBinLabel(i, std::to_string(i), 2);
0934 }
0935 me_GEM_pass_Ch_region_GE1_NoL->setAxisTitle("Chamber", 2);
0936 me_GEM_pass_Ch_region_GE1_NoL->setAxisTitle("Number of passing probes", 3);
0937
0938 for (int i = 1; i < 37; ++i) {
0939 me_GE11_fail_Ch_ieta->setBinLabel(i, std::to_string(i), 2);
0940 }
0941 for (int i = 1; i < 9; ++i) {
0942 me_GE11_fail_Ch_ieta->setBinLabel(i, std::to_string(i), 1);
0943 }
0944 me_GE11_fail_Ch_ieta->setAxisTitle("#ieta", 1);
0945 me_GE11_fail_Ch_ieta->setAxisTitle("Chamber", 2);
0946 me_GE11_fail_Ch_ieta->setAxisTitle("Number of failing probes", 3);
0947
0948 for (int i = 1; i < 37; ++i) {
0949 me_GE11_pass_Ch_ieta->setBinLabel(i, std::to_string(i), 2);
0950 }
0951 for (int i = 1; i < 9; ++i) {
0952 me_GE11_pass_Ch_ieta->setBinLabel(i, std::to_string(i), 1);
0953 }
0954 me_GE11_pass_Ch_ieta->setAxisTitle("#ieta", 1);
0955 me_GE11_pass_Ch_ieta->setAxisTitle("Chamber", 2);
0956 me_GE11_pass_Ch_ieta->setAxisTitle("Number of passing probes", 3);
0957
0958 for (int i = 1; i < 37; ++i) {
0959 me_GE11_fail_Ch_phi->setBinLabel(i, std::to_string(i), 2);
0960 }
0961 me_GE11_fail_Ch_phi->setAxisTitle("#phi", 1);
0962 me_GE11_fail_Ch_phi->setAxisTitle("Chamber", 2);
0963 me_GE11_fail_Ch_phi->setAxisTitle("Number of failing probes", 3);
0964
0965 for (int i = 1; i < 37; ++i) {
0966 me_GE11_pass_Ch_phi->setBinLabel(i, std::to_string(i), 2);
0967 }
0968 me_GE11_pass_Ch_phi->setAxisTitle("#phi", 1);
0969 me_GE11_pass_Ch_phi->setAxisTitle("Chamber", 2);
0970 me_GE11_pass_Ch_phi->setAxisTitle("Number of passing probes", 3);
0971
0972 for (int i = 1; i < 19; ++i) {
0973 me_GE21_fail_Ch_ieta->setBinLabel(i, std::to_string(i), 2);
0974 }
0975 for (int i = 1; i < 17; ++i) {
0976 me_GE21_fail_Ch_ieta->setBinLabel(i, std::to_string(i), 1);
0977 }
0978 me_GE21_fail_Ch_ieta->setAxisTitle("#ieta", 1);
0979 me_GE21_fail_Ch_ieta->setAxisTitle("Chamber", 2);
0980 me_GE21_fail_Ch_ieta->setAxisTitle("Number of failing probes", 3);
0981
0982 for (int i = 1; i < 19; ++i) {
0983 me_GE21_pass_Ch_ieta->setBinLabel(i, std::to_string(i), 2);
0984 }
0985 for (int i = 1; i < 17; ++i) {
0986 me_GE21_pass_Ch_ieta->setBinLabel(i, std::to_string(i), 1);
0987 }
0988 me_GE21_pass_Ch_ieta->setAxisTitle("#ieta", 1);
0989 me_GE21_pass_Ch_ieta->setAxisTitle("Chamber", 2);
0990 me_GE21_pass_Ch_ieta->setAxisTitle("Number of passing probes", 3);
0991
0992 for (int i = 1; i < 19; ++i) {
0993 me_GE21_fail_Ch_phi->setBinLabel(i, std::to_string(i), 2);
0994 }
0995 me_GE21_fail_Ch_phi->setAxisTitle("#phi", 1);
0996 me_GE21_fail_Ch_phi->setAxisTitle("Chamber", 2);
0997 me_GE21_fail_Ch_phi->setAxisTitle("Number of failing probes", 3);
0998
0999 for (int i = 1; i < 19; ++i) {
1000 me_GE21_pass_Ch_phi->setBinLabel(i, std::to_string(i), 2);
1001 }
1002 me_GE21_pass_Ch_phi->setAxisTitle("#phi", 1);
1003 me_GE21_pass_Ch_phi->setAxisTitle("Chamber", 2);
1004 me_GE21_pass_Ch_phi->setAxisTitle("Number of passing probes", 3);
1005
1006 for (int i = 1; i < 19; ++i) {
1007 me_ME0_pass_chamber_1D->setBinLabel(i, std::to_string(i), 1);
1008 }
1009 me_ME0_pass_chamber_1D->setAxisTitle("Chamber", 1);
1010 me_ME0_pass_chamber_1D->setAxisTitle("Number of passing probes", 2);
1011 for (int i = 1; i < 19; ++i) {
1012 me_ME0_fail_chamber_1D->setBinLabel(i, std::to_string(i), 1);
1013 }
1014 me_ME0_fail_chamber_1D->setAxisTitle("Chamber", 1);
1015 me_ME0_fail_chamber_1D->setAxisTitle("Number of failing probes", 2);
1016
1017 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(1, "GE-2/1_L2", 1);
1018 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(2, "GE-2/1_L1", 1);
1019 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(3, "GE-1/1_L2", 1);
1020 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(4, "GE-1/1_L1", 1);
1021 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(5, "GE0-", 1);
1022 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(6, "GE0+", 1);
1023 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(7, "GE+1/1_L1", 1);
1024 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(8, "GE+1/1_L2", 1);
1025 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(9, "GE+2/1_L1", 1);
1026 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(10, "GE+2/1_L2", 1);
1027 for (int i = 1; i < 37; ++i) {
1028 me_GEM_fail_Ch_region_layer_phase2->setBinLabel(i, std::to_string(i), 2);
1029 }
1030 me_GEM_fail_Ch_region_layer_phase2->setAxisTitle("Chamber", 2);
1031 me_GEM_fail_Ch_region_layer_phase2->setAxisTitle("Number of passing probes", 3);
1032
1033 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(1, "GE-2/1_L2", 1);
1034 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(2, "GE-2/1_L1", 1);
1035 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(3, "GE-1/1_L2", 1);
1036 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(4, "GE-1/1_L1", 1);
1037 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(5, "GE0-", 1);
1038 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(6, "GE0+", 1);
1039 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(7, "GE+1/1_L1", 1);
1040 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(8, "GE+1/1_L2", 1);
1041 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(9, "GE+2/1_L1", 1);
1042 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(10, "GE+2/1_L2", 1);
1043
1044 for (int i = 1; i < 37; ++i) {
1045 me_GEM_pass_Ch_region_layer_phase2->setBinLabel(i, std::to_string(i), 2);
1046 }
1047 me_GEM_pass_Ch_region_layer_phase2->setAxisTitle("Chamber", 2);
1048 me_GEM_pass_Ch_region_layer_phase2->setAxisTitle("Number of passing probes", 3);
1049
1050 for (int i = 1; i < 37; ++i) {
1051 me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
1052 }
1053 for (int i = 1; i < 9; ++i) {
1054 me_GE11_fail_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
1055 }
1056 me_GE11_fail_Ch_ieta_p1->setAxisTitle("#ieta", 1);
1057 me_GE11_fail_Ch_ieta_p1->setAxisTitle("Chamber", 2);
1058 me_GE11_fail_Ch_ieta_p1->setAxisTitle("Number of failing probes", 3);
1059
1060 for (int i = 1; i < 37; ++i) {
1061 me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 2);
1062 }
1063 for (int i = 1; i < 9; ++i) {
1064 me_GE11_pass_Ch_ieta_p1->setBinLabel(i, std::to_string(i), 1);
1065 }
1066 me_GE11_pass_Ch_ieta_p1->setAxisTitle("#ieta", 1);
1067 me_GE11_pass_Ch_ieta_p1->setAxisTitle("Chamber", 2);
1068 me_GE11_pass_Ch_ieta_p1->setAxisTitle("Number of passing probes", 3);
1069
1070 for (int i = 1; i < 37; ++i) {
1071 me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
1072 }
1073 for (int i = 1; i < 9; ++i) {
1074 me_GE11_fail_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
1075 }
1076 me_GE11_fail_Ch_ieta_p2->setAxisTitle("#ieta", 1);
1077 me_GE11_fail_Ch_ieta_p2->setAxisTitle("Chamber", 2);
1078 me_GE11_fail_Ch_ieta_p2->setAxisTitle("Number of failing probes", 3);
1079
1080 for (int i = 1; i < 37; ++i) {
1081 me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 2);
1082 }
1083 for (int i = 1; i < 9; ++i) {
1084 me_GE11_pass_Ch_ieta_p2->setBinLabel(i, std::to_string(i), 1);
1085 }
1086 me_GE11_pass_Ch_ieta_p2->setAxisTitle("#ieta", 1);
1087 me_GE11_pass_Ch_ieta_p2->setAxisTitle("Chamber", 2);
1088 me_GE11_pass_Ch_ieta_p2->setAxisTitle("Number of passing probes", 3);
1089
1090 for (int i = 1; i < 37; ++i) {
1091 me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
1092 }
1093 for (int i = 1; i < 9; ++i) {
1094 me_GE11_fail_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
1095 }
1096 me_GE11_fail_Ch_ieta_n1->setAxisTitle("#ieta", 1);
1097 me_GE11_fail_Ch_ieta_n1->setAxisTitle("Chamber", 2);
1098 me_GE11_fail_Ch_ieta_n1->setAxisTitle("Number of failing probes", 3);
1099
1100 for (int i = 1; i < 37; ++i) {
1101 me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 2);
1102 }
1103 for (int i = 1; i < 9; ++i) {
1104 me_GE11_pass_Ch_ieta_n1->setBinLabel(i, std::to_string(i), 1);
1105 }
1106 me_GE11_pass_Ch_ieta_n1->setAxisTitle("#ieta", 1);
1107 me_GE11_pass_Ch_ieta_n1->setAxisTitle("Chamber", 2);
1108 me_GE11_pass_Ch_ieta_n1->setAxisTitle("Number of passing probes", 3);
1109
1110 for (int i = 1; i < 37; ++i) {
1111 me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
1112 }
1113 for (int i = 1; i < 9; ++i) {
1114 me_GE11_fail_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
1115 }
1116 me_GE11_fail_Ch_ieta_n2->setAxisTitle("#ieta", 1);
1117 me_GE11_fail_Ch_ieta_n2->setAxisTitle("Chamber", 2);
1118 me_GE11_fail_Ch_ieta_n2->setAxisTitle("Number of failing probes", 3);
1119
1120 for (int i = 1; i < 37; ++i) {
1121 me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 2);
1122 }
1123 for (int i = 1; i < 9; ++i) {
1124 me_GE11_pass_Ch_ieta_n2->setBinLabel(i, std::to_string(i), 1);
1125 }
1126 me_GE11_pass_Ch_ieta_n2->setAxisTitle("#ieta", 1);
1127 me_GE11_pass_Ch_ieta_n2->setAxisTitle("Chamber", 2);
1128 me_GE11_pass_Ch_ieta_n2->setAxisTitle("Number of passing probes", 3);
1129
1130 m_histos["GE11_nPassingProbe_Ch_region"] = me_GE11_pass_Ch_region;
1131 m_histos["GE11_nFailingProbe_Ch_region"] = me_GE11_fail_Ch_region;
1132 m_histos["GE21_nPassingProbe_Ch_region"] = me_GE21_pass_Ch_region;
1133 m_histos["GE21_nFailingProbe_Ch_region"] = me_GE21_fail_Ch_region;
1134 m_histos["GEM_nPassingProbe_Ch_region_GE1"] = me_GEM_pass_Ch_region_GE1;
1135 m_histos["GEM_nFailingProbe_Ch_region_GE1"] = me_GEM_fail_Ch_region_GE1;
1136 m_histos["GEM_nPassingProbe_Ch_region_GE1_NoL"] = me_GEM_pass_Ch_region_GE1_NoL;
1137 m_histos["GEM_nFailingProbe_Ch_region_GE1_NoL"] = me_GEM_fail_Ch_region_GE1_NoL;
1138 m_histos["GE11_nPassingProbe_Ch_ieta"] = me_GE11_pass_Ch_ieta;
1139 m_histos["GE11_nFailingProbe_Ch_ieta"] = me_GE11_fail_Ch_ieta;
1140 m_histos["GE11_nPassingProbe_Ch_phi"] = me_GE11_pass_Ch_phi;
1141 m_histos["GE11_nFailingProbe_Ch_phi"] = me_GE11_fail_Ch_phi;
1142 m_histos["GE21_nPassingProbe_Ch_ieta"] = me_GE21_pass_Ch_ieta;
1143 m_histos["GE21_nFailingProbe_Ch_ieta"] = me_GE21_fail_Ch_ieta;
1144 m_histos["GE21_nPassingProbe_Ch_phi"] = me_GE21_pass_Ch_phi;
1145 m_histos["GE21_nFailingProbe_Ch_phi"] = me_GE21_fail_Ch_phi;
1146 m_histos["GE11_nPassingProbe_allCh_1D"] = me_GE11_pass_allCh_1D;
1147 m_histos["GE11_nFailingProbe_allCh_1D"] = me_GE11_fail_allCh_1D;
1148 m_histos["GE21_nPassingProbe_allCh_1D"] = me_GE21_pass_allCh_1D;
1149 m_histos["GE21_nFailingProbe_allCh_1D"] = me_GE21_fail_allCh_1D;
1150 m_histos["GE11_nPassingProbe_chamber_1D"] = me_GE11_pass_chamber_1D;
1151 m_histos["GE11_nFailingProbe_chamber_1D"] = me_GE11_fail_chamber_1D;
1152 m_histos["GE21_nPassingProbe_chamber_1D"] = me_GE21_pass_chamber_1D;
1153 m_histos["GE21_nFailingProbe_chamber_1D"] = me_GE21_fail_chamber_1D;
1154 m_histos["GEM_nPassingProbe_chamber_p1_1D"] = me_GEM_pass_chamber_p1_1D;
1155 m_histos["GEM_nFailingProbe_chamber_p1_1D"] = me_GEM_fail_chamber_p1_1D;
1156 m_histos["GEM_nPassingProbe_chamber_p2_1D"] = me_GEM_pass_chamber_p2_1D;
1157 m_histos["GEM_nFailingProbe_chamber_p2_1D"] = me_GEM_fail_chamber_p2_1D;
1158 m_histos["GEM_nPassingProbe_chamber_n1_1D"] = me_GEM_pass_chamber_n1_1D;
1159 m_histos["GEM_nFailingProbe_chamber_n1_1D"] = me_GEM_fail_chamber_n1_1D;
1160 m_histos["GEM_nPassingProbe_chamber_n2_1D"] = me_GEM_pass_chamber_n2_1D;
1161 m_histos["GEM_nFailingProbe_chamber_n2_1D"] = me_GEM_fail_chamber_n2_1D;
1162 m_histos["GEM_nPassingProbe_pt_1D"] = me_GEM_pass_pt_1D;
1163 m_histos["GEM_nFailingProbe_pt_1D"] = me_GEM_fail_pt_1D;
1164 m_histos["GEM_nPassingProbe_eta_1D"] = me_GEM_pass_eta_1D;
1165 m_histos["GEM_nFailingProbe_eta_1D"] = me_GEM_fail_eta_1D;
1166 m_histos["GEM_nPassingProbe_phi_1D"] = me_GEM_pass_phi_1D;
1167 m_histos["GEM_nFailingProbe_phi_1D"] = me_GEM_fail_phi_1D;
1168 m_histos["GEM_nPassingProbe_pt_p1_1D"] = me_GEM_pass_pt_p1_1D;
1169 m_histos["GEM_nFailingProbe_pt_p1_1D"] = me_GEM_fail_pt_p1_1D;
1170 m_histos["GEM_nPassingProbe_eta_p1_1D"] = me_GEM_pass_eta_p1_1D;
1171 m_histos["GEM_nFailingProbe_eta_p1_1D"] = me_GEM_fail_eta_p1_1D;
1172 m_histos["GEM_nPassingProbe_phi_p1_1D"] = me_GEM_pass_phi_p1_1D;
1173 m_histos["GEM_nFailingProbe_phi_p1_1D"] = me_GEM_fail_phi_p1_1D;
1174 m_histos["GEM_nPassingProbe_pt_p2_1D"] = me_GEM_pass_pt_p2_1D;
1175 m_histos["GEM_nFailingProbe_pt_p2_1D"] = me_GEM_fail_pt_p2_1D;
1176 m_histos["GEM_nPassingProbe_eta_p2_1D"] = me_GEM_pass_eta_p2_1D;
1177 m_histos["GEM_nFailingProbe_eta_p2_1D"] = me_GEM_fail_eta_p2_1D;
1178 m_histos["GEM_nPassingProbe_phi_p2_1D"] = me_GEM_pass_phi_p2_1D;
1179 m_histos["GEM_nFailingProbe_phi_p2_1D"] = me_GEM_fail_phi_p2_1D;
1180 m_histos["GEM_nPassingProbe_pt_n1_1D"] = me_GEM_pass_pt_n1_1D;
1181 m_histos["GEM_nFailingProbe_pt_n1_1D"] = me_GEM_fail_pt_n1_1D;
1182 m_histos["GEM_nPassingProbe_eta_n1_1D"] = me_GEM_pass_eta_n1_1D;
1183 m_histos["GEM_nFailingProbe_eta_n1_1D"] = me_GEM_fail_eta_n1_1D;
1184 m_histos["GEM_nPassingProbe_phi_n1_1D"] = me_GEM_pass_phi_n1_1D;
1185 m_histos["GEM_nFailingProbe_phi_n1_1D"] = me_GEM_fail_phi_n1_1D;
1186 m_histos["GEM_nPassingProbe_pt_n2_1D"] = me_GEM_pass_pt_n2_1D;
1187 m_histos["GEM_nFailingProbe_pt_n2_1D"] = me_GEM_fail_pt_n2_1D;
1188 m_histos["GEM_nPassingProbe_eta_n2_1D"] = me_GEM_pass_eta_n2_1D;
1189 m_histos["GEM_nFailingProbe_eta_n2_1D"] = me_GEM_fail_eta_n2_1D;
1190 m_histos["GEM_nPassingProbe_phi_n2_1D"] = me_GEM_pass_phi_n2_1D;
1191 m_histos["GEM_nFailingProbe_phi_n2_1D"] = me_GEM_fail_phi_n2_1D;
1192 m_histos["ME0_nPassingProbe_chamber_1D"] = me_ME0_pass_chamber_1D;
1193 m_histos["ME0_nFailingProbe_chamber_1D"] = me_ME0_fail_chamber_1D;
1194 m_histos["GEM_nPassingProbe_Ch_region_layer_phase2"] = me_GEM_pass_Ch_region_layer_phase2;
1195 m_histos["GEM_nFailingProbe_Ch_region_layer_phase2"] = me_GEM_fail_Ch_region_layer_phase2;
1196 m_histos["GE11_nPassingProbe_Ch_ieta_p1"] = me_GE11_pass_Ch_ieta_p1;
1197 m_histos["GE11_nFailingProbe_Ch_ieta_p1"] = me_GE11_fail_Ch_ieta_p1;
1198 m_histos["GE11_nPassingProbe_Ch_ieta_p2"] = me_GE11_pass_Ch_ieta_p2;
1199 m_histos["GE11_nFailingProbe_Ch_ieta_p2"] = me_GE11_fail_Ch_ieta_p2;
1200 m_histos["GE11_nPassingProbe_Ch_ieta_n1"] = me_GE11_pass_Ch_ieta_n1;
1201 m_histos["GE11_nFailingProbe_Ch_ieta_n1"] = me_GE11_fail_Ch_ieta_n1;
1202 m_histos["GE11_nPassingProbe_Ch_ieta_n2"] = me_GE11_pass_Ch_ieta_n2;
1203 m_histos["GE11_nFailingProbe_Ch_ieta_n2"] = me_GE11_fail_Ch_ieta_n2;
1204
1205 m_histos["GE11_nPassingProbe_OHmissing_allCh_1D"] = OHmissing_me_GE11_pass_allCh_1D;
1206 m_histos["GE11_nFailingProbe_OHmissing_allCh_1D"] = OHmissing_me_GE11_fail_allCh_1D;
1207 m_histos["GEM_nPassingProbe_OHmissing_chamber_p1_1D"] = OHmissing_me_GEM_pass_chamber_p1_1D;
1208 m_histos["GEM_nFailingProbe_OHmissing_chamber_p1_1D"] = OHmissing_me_GEM_fail_chamber_p1_1D;
1209 m_histos["GEM_nPassingProbe_OHmissing_chamber_p2_1D"] = OHmissing_me_GEM_pass_chamber_p2_1D;
1210 m_histos["GEM_nFailingProbe_OHmissing_chamber_p2_1D"] = OHmissing_me_GEM_fail_chamber_p2_1D;
1211 m_histos["GEM_nPassingProbe_OHmissing_chamber_n1_1D"] = OHmissing_me_GEM_pass_chamber_n1_1D;
1212 m_histos["GEM_nFailingProbe_OHmissing_chamber_n1_1D"] = OHmissing_me_GEM_fail_chamber_n1_1D;
1213 m_histos["GEM_nPassingProbe_OHmissing_chamber_n2_1D"] = OHmissing_me_GEM_pass_chamber_n2_1D;
1214 m_histos["GEM_nFailingProbe_OHmissing_chamber_n2_1D"] = OHmissing_me_GEM_fail_chamber_n2_1D;
1215 m_histos["GEM_nPassingProbe_OHmissing_pt_1D"] = OHmissing_me_GEM_pass_pt_1D;
1216 m_histos["GEM_nFailingProbe_OHmissing_pt_1D"] = OHmissing_me_GEM_fail_pt_1D;
1217 m_histos["GEM_nPassingProbe_OHmissing_eta_1D"] = OHmissing_me_GEM_pass_eta_1D;
1218 m_histos["GEM_nFailingProbe_OHmissing_eta_1D"] = OHmissing_me_GEM_fail_eta_1D;
1219 m_histos["GEM_nPassingProbe_OHmissing_phi_1D"] = OHmissing_me_GEM_pass_phi_1D;
1220 m_histos["GEM_nFailingProbe_OHmissing_phi_1D"] = OHmissing_me_GEM_fail_phi_1D;
1221 m_histos["GE11_nPassingProbe_OHmissing_Ch_ieta_p1"] = OHmissing_me_GE11_pass_Ch_ieta_p1;
1222 m_histos["GE11_nFailingProbe_OHmissing_Ch_ieta_p1"] = OHmissing_me_GE11_fail_Ch_ieta_p1;
1223 m_histos["GE11_nPassingProbe_OHmissing_Ch_ieta_p2"] = OHmissing_me_GE11_pass_Ch_ieta_p2;
1224 m_histos["GE11_nFailingProbe_OHmissing_Ch_ieta_p2"] = OHmissing_me_GE11_fail_Ch_ieta_p2;
1225 m_histos["GE11_nPassingProbe_OHmissing_Ch_ieta_n1"] = OHmissing_me_GE11_pass_Ch_ieta_n1;
1226 m_histos["GE11_nFailingProbe_OHmissing_Ch_ieta_n1"] = OHmissing_me_GE11_fail_Ch_ieta_n1;
1227 m_histos["GE11_nPassingProbe_OHmissing_Ch_ieta_n2"] = OHmissing_me_GE11_pass_Ch_ieta_n2;
1228 m_histos["GE11_nFailingProbe_OHmissing_Ch_ieta_n2"] = OHmissing_me_GE11_fail_Ch_ieta_n2;
1229
1230 m_histos["GE11_nPassingProbe_OHerror_allCh_1D"] = OHerror_me_GE11_pass_allCh_1D;
1231 m_histos["GE11_nFailingProbe_OHerror_allCh_1D"] = OHerror_me_GE11_fail_allCh_1D;
1232 m_histos["GEM_nPassingProbe_OHerror_chamber_p1_1D"] = OHerror_me_GEM_pass_chamber_p1_1D;
1233 m_histos["GEM_nFailingProbe_OHerror_chamber_p1_1D"] = OHerror_me_GEM_fail_chamber_p1_1D;
1234 m_histos["GEM_nPassingProbe_OHerror_chamber_p2_1D"] = OHerror_me_GEM_pass_chamber_p2_1D;
1235 m_histos["GEM_nFailingProbe_OHerror_chamber_p2_1D"] = OHerror_me_GEM_fail_chamber_p2_1D;
1236 m_histos["GEM_nPassingProbe_OHerror_chamber_n1_1D"] = OHerror_me_GEM_pass_chamber_n1_1D;
1237 m_histos["GEM_nFailingProbe_OHerror_chamber_n1_1D"] = OHerror_me_GEM_fail_chamber_n1_1D;
1238 m_histos["GEM_nPassingProbe_OHerror_chamber_n2_1D"] = OHerror_me_GEM_pass_chamber_n2_1D;
1239 m_histos["GEM_nFailingProbe_OHerror_chamber_n2_1D"] = OHerror_me_GEM_fail_chamber_n2_1D;
1240 m_histos["GEM_nPassingProbe_OHerror_pt_1D"] = OHerror_me_GEM_pass_pt_1D;
1241 m_histos["GEM_nFailingProbe_OHerror_pt_1D"] = OHerror_me_GEM_fail_pt_1D;
1242 m_histos["GEM_nPassingProbe_OHerror_eta_1D"] = OHerror_me_GEM_pass_eta_1D;
1243 m_histos["GEM_nFailingProbe_OHerror_eta_1D"] = OHerror_me_GEM_fail_eta_1D;
1244 m_histos["GEM_nPassingProbe_OHerror_phi_1D"] = OHerror_me_GEM_pass_phi_1D;
1245 m_histos["GEM_nFailingProbe_OHerror_phi_1D"] = OHerror_me_GEM_fail_phi_1D;
1246 m_histos["GE11_nPassingProbe_OHerror_Ch_ieta_p1"] = OHerror_me_GE11_pass_Ch_ieta_p1;
1247 m_histos["GE11_nFailingProbe_OHerror_Ch_ieta_p1"] = OHerror_me_GE11_fail_Ch_ieta_p1;
1248 m_histos["GE11_nPassingProbe_OHerror_Ch_ieta_p2"] = OHerror_me_GE11_pass_Ch_ieta_p2;
1249 m_histos["GE11_nFailingProbe_OHerror_Ch_ieta_p2"] = OHerror_me_GE11_fail_Ch_ieta_p2;
1250 m_histos["GE11_nPassingProbe_OHerror_Ch_ieta_n1"] = OHerror_me_GE11_pass_Ch_ieta_n1;
1251 m_histos["GE11_nFailingProbe_OHerror_Ch_ieta_n1"] = OHerror_me_GE11_fail_Ch_ieta_n1;
1252 m_histos["GE11_nPassingProbe_OHerror_Ch_ieta_n2"] = OHerror_me_GE11_pass_Ch_ieta_n2;
1253 m_histos["GE11_nFailingProbe_OHerror_Ch_ieta_n2"] = OHerror_me_GE11_fail_Ch_ieta_n2;
1254
1255 m_histos["GE11_nPassingProbe_VFATMask_allCh_1D"] = VFATMask_me_GE11_pass_allCh_1D;
1256 m_histos["GE11_nFailingProbe_VFATMask_allCh_1D"] = VFATMask_me_GE11_fail_allCh_1D;
1257 m_histos["GEM_nPassingProbe_VFATMask_chamber_p1_1D"] = VFATMask_me_GEM_pass_chamber_p1_1D;
1258 m_histos["GEM_nFailingProbe_VFATMask_chamber_p1_1D"] = VFATMask_me_GEM_fail_chamber_p1_1D;
1259 m_histos["GEM_nPassingProbe_VFATMask_chamber_p2_1D"] = VFATMask_me_GEM_pass_chamber_p2_1D;
1260 m_histos["GEM_nFailingProbe_VFATMask_chamber_p2_1D"] = VFATMask_me_GEM_fail_chamber_p2_1D;
1261 m_histos["GEM_nPassingProbe_VFATMask_chamber_n1_1D"] = VFATMask_me_GEM_pass_chamber_n1_1D;
1262 m_histos["GEM_nFailingProbe_VFATMask_chamber_n1_1D"] = VFATMask_me_GEM_fail_chamber_n1_1D;
1263 m_histos["GEM_nPassingProbe_VFATMask_chamber_n2_1D"] = VFATMask_me_GEM_pass_chamber_n2_1D;
1264 m_histos["GEM_nFailingProbe_VFATMask_chamber_n2_1D"] = VFATMask_me_GEM_fail_chamber_n2_1D;
1265 m_histos["GEM_nPassingProbe_VFATMask_pt_1D"] = VFATMask_me_GEM_pass_pt_1D;
1266 m_histos["GEM_nFailingProbe_VFATMask_pt_1D"] = VFATMask_me_GEM_fail_pt_1D;
1267 m_histos["GEM_nPassingProbe_VFATMask_eta_1D"] = VFATMask_me_GEM_pass_eta_1D;
1268 m_histos["GEM_nFailingProbe_VFATMask_eta_1D"] = VFATMask_me_GEM_fail_eta_1D;
1269 m_histos["GEM_nPassingProbe_VFATMask_phi_1D"] = VFATMask_me_GEM_pass_phi_1D;
1270 m_histos["GEM_nFailingProbe_VFATMask_phi_1D"] = VFATMask_me_GEM_fail_phi_1D;
1271 m_histos["GE11_nPassingProbe_VFATMask_Ch_ieta_p1"] = VFATMask_me_GE11_pass_Ch_ieta_p1;
1272 m_histos["GE11_nFailingProbe_VFATMask_Ch_ieta_p1"] = VFATMask_me_GE11_fail_Ch_ieta_p1;
1273 m_histos["GE11_nPassingProbe_VFATMask_Ch_ieta_p2"] = VFATMask_me_GE11_pass_Ch_ieta_p2;
1274 m_histos["GE11_nFailingProbe_VFATMask_Ch_ieta_p2"] = VFATMask_me_GE11_fail_Ch_ieta_p2;
1275 m_histos["GE11_nPassingProbe_VFATMask_Ch_ieta_n1"] = VFATMask_me_GE11_pass_Ch_ieta_n1;
1276 m_histos["GE11_nFailingProbe_VFATMask_Ch_ieta_n1"] = VFATMask_me_GE11_fail_Ch_ieta_n1;
1277 m_histos["GE11_nPassingProbe_VFATMask_Ch_ieta_n2"] = VFATMask_me_GE11_pass_Ch_ieta_n2;
1278 m_histos["GE11_nFailingProbe_VFATMask_Ch_ieta_n2"] = VFATMask_me_GE11_fail_Ch_ieta_n2;
1279
1280 std::string baseDir_ = topFolder() + "/detailed/";
1281 iBooker.setCurrentFolder(baseDir_);
1282 m_histos["GEMseg_dx_ME0"] = iBooker.book1D("GEMseg_dx_ME0", "GEMseg_dx;probe dx [cm];Events", 100, 0., 20.);
1283 m_histos["GEMhit_dx_GE1"] = iBooker.book1D("GEMhit_dx_GE1", "GEMhit_dx;probe dx [cm];Events", 100, 0., 10.);
1284 m_histos["GEMhit_dx_GE2"] = iBooker.book1D("GEMhit_dx_GE2", "GEMhit_dx;probe dx [cm];Events", 100, 0., 10.);
1285
1286 m_histos["GEMseg_x_ME0"] = iBooker.book1D("GEMhit_x_ME0", "GEMhit_x;probe x [cm];Events", 100, -10., 10.);
1287 m_histos["GEMhit_x_GE1"] = iBooker.book1D("GEMhit_x_GE1", "GEMhit_x;probe x [cm];Events", 100, -10., 10.);
1288 m_histos["GEMhit_x_GE2"] = iBooker.book1D("GEMhit_x_GE2", "GEMhit_x;probe x [cm];Events", 100, -10., 10.);
1289 m_histos["Cham_x_ME0"] = iBooker.book1D("Cham_x_ME0", "Cham_x;probe x [cm];Events", 100, -10., 10.);
1290 m_histos["Cham_x_GE1"] = iBooker.book1D("Cham_x_GE1", "Cham_x;probe x [cm];Events", 100, -10., 10.);
1291 m_histos["Cham_x_GE2"] = iBooker.book1D("Cham_x_GE2", "Cham_x;probe x [cm];Events", 100, -10., 10.);
1292 m_histos["xyErr_GE1"] = iBooker.book2D("xyErr_GE1", "xyErr_GE1", 50, 0., 5., 50, 0., 5.);
1293 }
1294
1295 uint16_t GEMTnPEfficiencyTask::maskChamberWithError(const GEMDetId& chamber_id,
1296 const GEMOHStatusCollection* oh_status_collection,
1297 const GEMVFATStatusCollection* vfat_status_collection) {
1298 uint16_t oh_warning = 0;
1299 bool oh_exists = false;
1300 for (auto iter = oh_status_collection->begin(); iter != oh_status_collection->end(); iter++) {
1301 const auto [oh_id, range] = (*iter);
1302 if (chamber_id.chamberId() != oh_id) {
1303 continue;
1304 }
1305 for (auto oh_status = range.first; oh_status != range.second; oh_status++) {
1306 oh_exists = true;
1307 if (oh_status->isBad()) {
1308 oh_warning = oh_warning | (1 << 1);
1309 }
1310
1311 uint32_t vfatmask = oh_status->vfatMask();
1312 if (vfatmask != 16777215) {
1313 int ieta = chamber_id.ieta();
1314 if (!((vfatmask >> (8 - ieta) & 1) && (vfatmask >> (16 - ieta) & 1) &&
1315 (vfatmask >> (24 - ieta) & 1))) {
1316 oh_warning = oh_warning | (1 << 2);
1317 }
1318 }
1319 }
1320 }
1321 if (!oh_exists) {
1322 oh_warning = oh_warning | 1;
1323 }
1324 return oh_warning;
1325 }
1326
1327 bool GEMTnPEfficiencyTask::checkBounds(const GeomDet* geomDet,
1328 const GlobalPoint& global_position,
1329 const float bordercut) {
1330 const TrapezoidalPlaneBounds* bounds = dynamic_cast<const TrapezoidalPlaneBounds*>(&geomDet->surface().bounds());
1331 LocalPoint localPoint = geomDet->surface().toLocal(global_position);
1332 float wideWidth = bounds->width();
1333 float narrowWidth = 2.f * bounds->widthAtHalfLength() - wideWidth;
1334 float length = bounds->length();
1335 float tangent = (wideWidth - narrowWidth) / (2.f * length);
1336 float halfWidthAtY = tangent * localPoint.y() + 0.25f * (narrowWidth + wideWidth);
1337 float distanceY = std::abs(localPoint.y()) - 0.5f * length;
1338 float distanceX = std::abs(localPoint.x()) - halfWidthAtY;
1339 if (distanceX < bordercut && distanceY < bordercut) {
1340 return true;
1341 }
1342 return false;
1343 }
1344
1345 void GEMTnPEfficiencyTask::analyze(const edm::Event& event, const edm::EventSetup& context) {
1346 BaseTnPEfficiencyTask::analyze(event, context);
1347 GEMOHStatusCollection oh_status;
1348 GEMVFATStatusCollection vfat_status;
1349 edm::Handle<GEMOHStatusCollection> oh_status_collection;
1350 edm::Handle<GEMVFATStatusCollection> vfat_status_collection;
1351 muon_service_->update(context);
1352 if (m_maskChamberWithError_) {
1353 event.getByToken(m_GEMOHStatusCollectionToken_, oh_status_collection);
1354
1355 oh_status = *oh_status_collection;
1356 } else {
1357 LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask") << "failed to get GEMOHStatusCollection" << std::endl;
1358 return;
1359 }
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369 edm::Handle<reco::MuonCollection> muons;
1370 event.getByToken(m_muToken, muons);
1371
1372
1373 std::vector<std::vector<int>> probe_coll_GE11_region;
1374 std::vector<std::vector<int>> probe_coll_GE11_lay;
1375 std::vector<std::vector<int>> probe_coll_GE11_chamber;
1376 std::vector<std::vector<float>> probe_coll_GE11_pt;
1377 std::vector<std::vector<float>> probe_coll_GE11_eta;
1378 std::vector<std::vector<float>> probe_coll_GE11_ieta;
1379 std::vector<std::vector<float>> probe_coll_GE11_phi;
1380 std::vector<std::vector<int>> probe_coll_GE11_sta;
1381 std::vector<std::vector<float>> probe_coll_GE11_dx;
1382 std::vector<std::vector<uint16_t>> probe_coll_GE11_warnings;
1383
1384
1385 std::vector<std::vector<int>> probe_coll_GE21_region;
1386 std::vector<std::vector<int>> probe_coll_GE21_lay;
1387 std::vector<std::vector<int>> probe_coll_GE21_chamber;
1388 std::vector<std::vector<float>> probe_coll_GE21_pt;
1389 std::vector<std::vector<float>> probe_coll_GE21_eta;
1390 std::vector<std::vector<float>> probe_coll_GE21_ieta;
1391 std::vector<std::vector<float>> probe_coll_GE21_phi;
1392 std::vector<std::vector<int>> probe_coll_GE21_sta;
1393 std::vector<std::vector<float>> probe_coll_GE21_dx;
1394 std::vector<std::vector<uint16_t>> probe_coll_GE21_warnings;
1395
1396 std::vector<uint8_t> probe_coll_GEM_staMatch;
1397
1398
1399 std::vector<std::vector<int>> probe_coll_ME0_region;
1400 std::vector<std::vector<int>> probe_coll_ME0_roll;
1401 std::vector<std::vector<int>> probe_coll_ME0_lay;
1402 std::vector<std::vector<int>> probe_coll_ME0_chamber;
1403 std::vector<std::vector<float>> probe_coll_ME0_pt;
1404 std::vector<std::vector<float>> probe_coll_ME0_eta;
1405 std::vector<std::vector<float>> probe_coll_ME0_ieta;
1406 std::vector<std::vector<float>> probe_coll_ME0_phi;
1407 std::vector<std::vector<int>> probe_coll_ME0_sta;
1408 std::vector<std::vector<float>> probe_coll_ME0_dx;
1409 std::vector<std::vector<uint16_t>> probe_coll_ME0_warnings;
1410
1411 std::vector<unsigned> probe_indices;
1412 if (!m_probeIndices.empty())
1413 probe_indices = m_probeIndices.back();
1414
1415
1416 for (const auto i : probe_indices) {
1417
1418 std::vector<int> probe_GE11_region;
1419 std::vector<int> probe_GE11_sta;
1420 std::vector<int> probe_GE11_lay;
1421 std::vector<int> probe_GE11_chamber;
1422 std::vector<float> probe_GE11_pt;
1423 std::vector<float> probe_GE11_eta;
1424 std::vector<float> probe_GE11_ieta;
1425 std::vector<float> probe_GE11_phi;
1426 std::vector<float> probe_GE11_dx;
1427 std::vector<uint16_t> probe_GE11_warnings;
1428
1429 std::vector<int> probe_GE21_region;
1430 std::vector<int> probe_GE21_sta;
1431 std::vector<int> probe_GE21_lay;
1432 std::vector<int> probe_GE21_chamber;
1433 std::vector<float> probe_GE21_pt;
1434 std::vector<float> probe_GE21_eta;
1435 std::vector<float> probe_GE21_ieta;
1436 std::vector<float> probe_GE21_phi;
1437 std::vector<float> probe_GE21_dx;
1438 std::vector<uint16_t> probe_GE21_warnings;
1439
1440 uint8_t GEM_stationMatching = 0;
1441
1442 std::vector<int> probe_ME0_region;
1443 std::vector<int> probe_ME0_roll;
1444 std::vector<int> probe_ME0_sta;
1445 std::vector<int> probe_ME0_lay;
1446 std::vector<int> probe_ME0_chamber;
1447 std::vector<float> probe_ME0_pt;
1448 std::vector<float> probe_ME0_eta;
1449 std::vector<float> probe_ME0_ieta;
1450 std::vector<float> probe_ME0_phi;
1451 std::vector<float> probe_ME0_dx;
1452 std::vector<uint16_t> probe_ME0_warnings;
1453
1454 bool gem_matched = false;
1455
1456 for (const auto& chambMatch : (*muons).at(i).matches()) {
1457
1458 bool hit_matched = false;
1459 if (chambMatch.detector() == MuonSubdetId::GEM) {
1460 if (chambMatch.edgeX < m_borderCut && chambMatch.edgeY < m_borderCut) {
1461 gem_matched = true;
1462
1463 GEMDetId chId(chambMatch.id.rawId());
1464 const uint16_t warnings = maskChamberWithError(chId, &oh_status, &vfat_status);
1465 const int roll = chId.roll();
1466 const int region = chId.region();
1467 const int station = chId.station();
1468 const int layer = chId.layer();
1469 const int chamber = chId.chamber();
1470 const float pt = (*muons).at(i).pt();
1471 const float eta = (*muons).at(i).eta();
1472 const float phi = (*muons).at(i).phi();
1473 int ieta = 0;
1474 GEM_stationMatching = GEM_stationMatching | (1 << station);
1475
1476 const GeomDet* geomDet = muon_service_->trackingGeometry()->idToDet(chId);
1477 LocalPoint pos(chambMatch.x, chambMatch.y);
1478
1479 const GlobalPoint& global_position = geomDet->toGlobal(pos);
1480
1481 if (const GEMChamber* gemChamber = dynamic_cast<const GEMChamber*>(geomDet)) {
1482 for (const GEMEtaPartition* eta_partition : gemChamber->etaPartitions())
1483 if (checkBounds(eta_partition, global_position, m_borderCut)) {
1484 ieta = eta_partition->id().ieta();
1485 break;
1486 }
1487 }
1488
1489 if (station == 1 || station == 2) {
1490 reco::MuonGEMHitMatch closest_matchedHit;
1491 double smallestDx = 99999.;
1492 double matched_GEMHit_x = 99999.;
1493
1494 for (auto& gemHit : chambMatch.gemHitMatches) {
1495 float dx = std::abs(chambMatch.x - gemHit.x);
1496 if (dx < smallestDx) {
1497 smallestDx = dx;
1498 closest_matchedHit = gemHit;
1499 matched_GEMHit_x = gemHit.x;
1500 hit_matched = true;
1501 }
1502 }
1503
1504 if (station == 1) {
1505 probe_GE11_region.push_back(region);
1506 probe_GE11_sta.push_back(station);
1507 probe_GE11_lay.push_back(layer);
1508 probe_GE11_chamber.push_back(chamber);
1509 probe_GE11_ieta.push_back(ieta);
1510 probe_GE11_pt.push_back(pt);
1511 probe_GE11_eta.push_back(eta);
1512 probe_GE11_phi.push_back(phi);
1513 probe_GE11_dx.push_back(smallestDx);
1514 probe_GE11_warnings.push_back(warnings);
1515 }
1516
1517 if (station == 2) {
1518 probe_GE21_region.push_back(region);
1519 probe_GE21_sta.push_back(station);
1520 probe_GE21_lay.push_back(layer);
1521 probe_GE21_chamber.push_back(chamber);
1522 probe_GE21_ieta.push_back(ieta);
1523 probe_GE21_pt.push_back(pt);
1524 probe_GE21_eta.push_back(eta);
1525 probe_GE21_phi.push_back(phi);
1526 probe_GE21_dx.push_back(smallestDx);
1527 probe_GE21_warnings.push_back(warnings);
1528 }
1529 if (m_detailedAnalysis && hit_matched) {
1530 if (station == 1) {
1531 m_histos.find("GEMhit_dx_GE1")->second->Fill(smallestDx);
1532 m_histos.find("GEMhit_x_GE1")->second->Fill(matched_GEMHit_x);
1533 m_histos.find("Cham_x_GE1")->second->Fill(chambMatch.x);
1534 m_histos.find("xyErr_GE1")->second->Fill(chambMatch.xErr, chambMatch.yErr);
1535 }
1536 if (station == 2) {
1537 m_histos.find("GEMhit_dx_GE2")->second->Fill(smallestDx);
1538 m_histos.find("GEMhit_x_GE2")->second->Fill(matched_GEMHit_x);
1539 m_histos.find("Cham_x_GE2")->second->Fill(chambMatch.x);
1540 }
1541 }
1542 }
1543
1544 if (station == 0) {
1545 reco::MuonSegmentMatch closest_matchedSegment;
1546 double smallestDx_seg = 99999.;
1547
1548 for (auto& seg : chambMatch.gemMatches) {
1549 float dx_seg = std::abs(chambMatch.x - seg.x);
1550 if (dx_seg < smallestDx_seg) {
1551 smallestDx_seg = dx_seg;
1552 closest_matchedSegment = seg;
1553 hit_matched = true;
1554 }
1555 }
1556
1557 probe_ME0_region.push_back(region);
1558 probe_ME0_roll.push_back(roll);
1559 probe_ME0_sta.push_back(station);
1560 probe_ME0_lay.push_back(layer);
1561 probe_ME0_chamber.push_back(chamber);
1562 probe_ME0_ieta.push_back(ieta);
1563 probe_ME0_pt.push_back(pt);
1564 probe_ME0_eta.push_back(eta);
1565 probe_ME0_phi.push_back(phi);
1566 probe_ME0_dx.push_back(smallestDx_seg);
1567 probe_ME0_warnings.push_back(warnings);
1568
1569 if (m_detailedAnalysis && hit_matched) {
1570 m_histos.find("GEMseg_dx_ME0")->second->Fill(smallestDx_seg);
1571 m_histos.find("GEMseg_x_ME0")->second->Fill(closest_matchedSegment.x);
1572 m_histos.find("Cham_x_ME0")->second->Fill(chambMatch.x);
1573 }
1574 }
1575 }
1576 } else
1577 continue;
1578 }
1579
1580
1581 if (m_detailedAnalysis && gem_matched) {
1582 m_histos.find("probeEta")->second->Fill((*muons).at(i).eta());
1583 m_histos.find("probePhi")->second->Fill((*muons).at(i).phi());
1584 m_histos.find("probeNumberOfMatchedStations")->second->Fill((*muons).at(i).numberOfMatchedStations());
1585 m_histos.find("probePt")->second->Fill((*muons).at(i).pt());
1586
1587
1588
1589
1590
1591 }
1592
1593
1594 probe_coll_GE11_region.push_back(probe_GE11_region);
1595 probe_coll_GE11_sta.push_back(probe_GE11_sta);
1596 probe_coll_GE11_lay.push_back(probe_GE11_lay);
1597 probe_coll_GE11_chamber.push_back(probe_GE11_chamber);
1598 probe_coll_GE11_ieta.push_back(probe_GE11_ieta);
1599 probe_coll_GE11_pt.push_back(probe_GE11_pt);
1600 probe_coll_GE11_eta.push_back(probe_GE11_eta);
1601 probe_coll_GE11_phi.push_back(probe_GE11_phi);
1602 probe_coll_GE11_dx.push_back(probe_GE11_dx);
1603 probe_coll_GE11_warnings.push_back(probe_GE11_warnings);
1604
1605 probe_coll_GEM_staMatch.push_back(GEM_stationMatching);
1606
1607
1608 probe_coll_GE21_region.push_back(probe_GE21_region);
1609 probe_coll_GE21_sta.push_back(probe_GE21_sta);
1610 probe_coll_GE21_lay.push_back(probe_GE21_lay);
1611 probe_coll_GE21_chamber.push_back(probe_GE21_chamber);
1612 probe_coll_GE21_ieta.push_back(probe_GE21_ieta);
1613 probe_coll_GE21_pt.push_back(probe_GE21_pt);
1614 probe_coll_GE21_eta.push_back(probe_GE21_eta);
1615 probe_coll_GE21_phi.push_back(probe_GE21_phi);
1616 probe_coll_GE21_dx.push_back(probe_GE21_dx);
1617 probe_coll_GE21_warnings.push_back(probe_GE21_warnings);
1618
1619
1620 probe_coll_ME0_region.push_back(probe_ME0_region);
1621 probe_coll_ME0_roll.push_back(probe_ME0_roll);
1622 probe_coll_ME0_sta.push_back(probe_ME0_sta);
1623 probe_coll_ME0_lay.push_back(probe_ME0_lay);
1624 probe_coll_ME0_chamber.push_back(probe_ME0_chamber);
1625 probe_coll_ME0_ieta.push_back(probe_ME0_ieta);
1626 probe_coll_ME0_pt.push_back(probe_ME0_pt);
1627 probe_coll_ME0_eta.push_back(probe_ME0_eta);
1628 probe_coll_ME0_phi.push_back(probe_ME0_phi);
1629 probe_coll_ME0_dx.push_back(probe_ME0_dx);
1630 probe_coll_ME0_warnings.push_back(probe_ME0_warnings);
1631
1632 }
1633
1634
1635 for (unsigned i = 0; i < probe_indices.size(); ++i) {
1636
1637
1638
1639 unsigned nME0_matches = probe_coll_ME0_region.at(i).size();
1640 for (unsigned j = 0; j < nME0_matches; ++j) {
1641
1642 int ME0_region = probe_coll_ME0_region.at(i).at(j);
1643
1644
1645
1646 int ME0_chamber = probe_coll_ME0_chamber.at(i).at(j);
1647
1648 float ME0_dx = probe_coll_ME0_dx.at(i).at(j);
1649
1650
1651
1652 if (ME0_dx < m_dxCut) {
1653 m_histos.find("ME0_nPassingProbe_chamber_1D")->second->Fill(ME0_chamber);
1654 if (ME0_region < 0)
1655 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(4, ME0_chamber);
1656 else if (ME0_region > 0)
1657 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(5, ME0_chamber);
1658 } else {
1659 m_histos.find("ME0_nFailingProbe_chamber_1D")->second->Fill(ME0_chamber);
1660 if (ME0_region < 0)
1661 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(4, ME0_chamber);
1662 else if (ME0_region > 0)
1663 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(5, ME0_chamber);
1664 }
1665 }
1666
1667
1668
1669 unsigned nGE11_matches = probe_coll_GE11_region.at(i).size();
1670 for (unsigned j = 0; j < nGE11_matches; ++j) {
1671
1672 int GEM_region = probe_coll_GE11_region.at(i).at(j);
1673 int GEM_sta = probe_coll_GE11_sta.at(i).at(j);
1674 int GEM_lay = probe_coll_GE11_lay.at(i).at(j);
1675 int GEM_chamber = probe_coll_GE11_chamber.at(i).at(j);
1676 int GEM_ieta = probe_coll_GE11_ieta.at(i).at(j);
1677 float GEM_pt = probe_coll_GE11_pt.at(i).at(j);
1678 float GEM_dx = probe_coll_GE11_dx.at(i).at(j);
1679 float GEM_eta = probe_coll_GE11_eta.at(i).at(j);
1680 float GEM_phi = probe_coll_GE11_phi.at(i).at(j);
1681
1682 uint16_t GEM_warning = probe_coll_GE11_warnings.at(i).at(j);
1683
1684 if (GEM_dx < m_dxCut) {
1685 m_histos.find("GE11_nPassingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
1686 m_histos.find("GE11_nPassingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
1687 m_histos.find("GE11_nPassingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
1688 m_histos.find("GE11_nPassingProbe_allCh_1D")->second->Fill(GEM_region);
1689 if (~GEM_warning & 1) {
1690 m_histos.find("GE11_nPassingProbe_OHmissing_allCh_1D")->second->Fill(GEM_region);
1691 if (~GEM_warning >> 1 & 1) {
1692 m_histos.find("GE11_nPassingProbe_OHerror_allCh_1D")->second->Fill(GEM_region);
1693 if (~GEM_warning >> 2 & 1) {
1694 m_histos.find("GE11_nPassingProbe_VFATMask_allCh_1D")->second->Fill(GEM_region);
1695 }
1696 }
1697 }
1698
1699 m_histos.find("GE11_nPassingProbe_chamber_1D")->second->Fill(GEM_chamber);
1700 if (GEM_region < 0) {
1701 if (GEM_lay == 2)
1702 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(2, GEM_chamber);
1703 else if (GEM_lay == 1)
1704 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(3, GEM_chamber);
1705 }
1706 if (GEM_region > 0) {
1707 if (GEM_lay == 1)
1708 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(6, GEM_chamber);
1709 else if (GEM_lay == 2)
1710 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(7, GEM_chamber);
1711 }
1712 if (GEM_region == -1) {
1713 m_histos.find("GEM_nPassingProbe_Ch_region_GE1_NoL")->second->Fill(0, GEM_chamber);
1714 } else if (GEM_region == 1) {
1715 m_histos.find("GEM_nPassingProbe_Ch_region_GE1_NoL")->second->Fill(1, GEM_chamber);
1716 }
1717
1718 if (GEM_region == 1 && GEM_lay == 1) {
1719 m_histos.find("GEM_nPassingProbe_chamber_p1_1D")->second->Fill(GEM_chamber);
1720 if (~GEM_warning & 1) {
1721 m_histos.find("GEM_nPassingProbe_OHmissing_chamber_p1_1D")->second->Fill(GEM_chamber);
1722 if (~GEM_warning >> 1 & 1) {
1723 m_histos.find("GEM_nPassingProbe_OHerror_chamber_p1_1D")->second->Fill(GEM_chamber);
1724 if (~GEM_warning >> 2 & 1) {
1725 m_histos.find("GEM_nPassingProbe_VFATMask_chamber_p1_1D")->second->Fill(GEM_chamber);
1726 }
1727 }
1728 }
1729 if (GEM_sta == 1) {
1730 m_histos.find("GE11_nPassingProbe_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1731 if (~GEM_warning & 1) {
1732 m_histos.find("GE11_nPassingProbe_OHmissing_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1733 if (~GEM_warning >> 1 & 1) {
1734 m_histos.find("GE11_nPassingProbe_OHerror_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1735 if (~GEM_warning >> 2 & 1) {
1736 m_histos.find("GE11_nPassingProbe_VFATMask_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1737 }
1738 }
1739 }
1740 }
1741 m_histos.find("GEM_nPassingProbe_Ch_region_GE1")->second->Fill(2, GEM_chamber);
1742 m_histos.find("GEM_nPassingProbe_pt_p1_1D")->second->Fill(GEM_pt);
1743 m_histos.find("GEM_nPassingProbe_eta_p1_1D")->second->Fill(abs(GEM_eta));
1744 m_histos.find("GEM_nPassingProbe_phi_p1_1D")->second->Fill(GEM_phi);
1745 } else if (GEM_region == 1 && GEM_lay == 2) {
1746 m_histos.find("GEM_nPassingProbe_chamber_p2_1D")->second->Fill(GEM_chamber);
1747 if (~GEM_warning & 1) {
1748 m_histos.find("GEM_nPassingProbe_OHmissing_chamber_p2_1D")->second->Fill(GEM_chamber);
1749 if (~GEM_warning >> 1 & 1) {
1750 m_histos.find("GEM_nPassingProbe_OHerror_chamber_p2_1D")->second->Fill(GEM_chamber);
1751 if (~GEM_warning >> 2 & 1) {
1752 m_histos.find("GEM_nPassingProbe_VFATMask_chamber_p2_1D")->second->Fill(GEM_chamber);
1753 }
1754 }
1755 }
1756 if (GEM_sta == 1) {
1757 m_histos.find("GE11_nPassingProbe_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1758 if (~GEM_warning & 1) {
1759 m_histos.find("GE11_nPassingProbe_OHmissing_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1760 if (~GEM_warning >> 1 & 1) {
1761 m_histos.find("GE11_nPassingProbe_OHerror_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1762 if (~GEM_warning >> 2 & 1) {
1763 m_histos.find("GE11_nPassingProbe_VFATMask_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1764 }
1765 }
1766 }
1767 }
1768 m_histos.find("GEM_nPassingProbe_Ch_region_GE1")->second->Fill(3, GEM_chamber);
1769 m_histos.find("GEM_nPassingProbe_pt_p2_1D")->second->Fill(GEM_pt);
1770 m_histos.find("GEM_nPassingProbe_eta_p2_1D")->second->Fill(abs(GEM_eta));
1771 m_histos.find("GEM_nPassingProbe_phi_p2_1D")->second->Fill(GEM_phi);
1772 } else if (GEM_region == -1 && GEM_lay == 1) {
1773 m_histos.find("GEM_nPassingProbe_chamber_n1_1D")->second->Fill(GEM_chamber);
1774 if (~GEM_warning & 1) {
1775 m_histos.find("GEM_nPassingProbe_OHmissing_chamber_n1_1D")->second->Fill(GEM_chamber);
1776 if (~GEM_warning >> 1 & 1) {
1777 m_histos.find("GEM_nPassingProbe_OHerror_chamber_n1_1D")->second->Fill(GEM_chamber);
1778 if (~GEM_warning >> 2 & 1) {
1779 m_histos.find("GEM_nPassingProbe_VFATMask_chamber_n1_1D")->second->Fill(GEM_chamber);
1780 }
1781 }
1782 }
1783 if (GEM_sta == 1) {
1784 m_histos.find("GE11_nPassingProbe_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1785 if (~GEM_warning & 1) {
1786 m_histos.find("GE11_nPassingProbe_OHmissing_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1787 if (~GEM_warning >> 1 & 1) {
1788 m_histos.find("GE11_nPassingProbe_OHerror_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1789 if (~GEM_warning >> 2 & 1) {
1790 m_histos.find("GE11_nPassingProbe_VFATMask_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1791 }
1792 }
1793 }
1794 }
1795 m_histos.find("GEM_nPassingProbe_Ch_region_GE1")->second->Fill(1, GEM_chamber);
1796 m_histos.find("GEM_nPassingProbe_pt_n1_1D")->second->Fill(GEM_pt);
1797 m_histos.find("GEM_nPassingProbe_eta_n1_1D")->second->Fill(abs(GEM_eta));
1798 m_histos.find("GEM_nPassingProbe_phi_n1_1D")->second->Fill(GEM_phi);
1799 } else if (GEM_region == -1 && GEM_lay == 2) {
1800 m_histos.find("GEM_nPassingProbe_chamber_n2_1D")->second->Fill(GEM_chamber);
1801 if (~GEM_warning & 1) {
1802 m_histos.find("GEM_nPassingProbe_OHmissing_chamber_n2_1D")->second->Fill(GEM_chamber);
1803 if (~GEM_warning >> 1 & 1) {
1804 m_histos.find("GEM_nPassingProbe_OHerror_chamber_n2_1D")->second->Fill(GEM_chamber);
1805 if (~GEM_warning >> 2 & 1) {
1806 m_histos.find("GEM_nPassingProbe_VFATMask_chamber_n2_1D")->second->Fill(GEM_chamber);
1807 }
1808 }
1809 }
1810 if (GEM_sta == 1) {
1811 m_histos.find("GE11_nPassingProbe_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1812 if (~GEM_warning & 1) {
1813 m_histos.find("GE11_nPassingProbe_OHmissing_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1814 if (~GEM_warning >> 1 & 1) {
1815 m_histos.find("GE11_nPassingProbe_OHerror_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1816 if (~GEM_warning >> 2 & 1) {
1817 m_histos.find("GE11_nPassingProbe_VFATMask_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1818 }
1819 }
1820 }
1821 }
1822 m_histos.find("GEM_nPassingProbe_Ch_region_GE1")->second->Fill(0, GEM_chamber);
1823 m_histos.find("GEM_nPassingProbe_pt_n2_1D")->second->Fill(GEM_pt);
1824 m_histos.find("GEM_nPassingProbe_eta_n2_1D")->second->Fill(abs(GEM_eta));
1825 m_histos.find("GEM_nPassingProbe_phi_n2_1D")->second->Fill(GEM_phi);
1826 }
1827 m_histos.find("GEM_nPassingProbe_pt_1D")->second->Fill(GEM_pt);
1828 if (~GEM_warning & 1) {
1829 m_histos.find("GEM_nPassingProbe_OHmissing_pt_1D")->second->Fill(GEM_pt);
1830 if (~GEM_warning >> 1 & 1) {
1831 m_histos.find("GEM_nPassingProbe_OHerror_pt_1D")->second->Fill(GEM_pt);
1832 if (~GEM_warning >> 2 & 1) {
1833 m_histos.find("GEM_nPassingProbe_VFATMask_pt_1D")->second->Fill(GEM_pt);
1834 }
1835 }
1836 }
1837 m_histos.find("GEM_nPassingProbe_eta_1D")->second->Fill(abs(GEM_eta));
1838 if (~GEM_warning & 1) {
1839 m_histos.find("GEM_nPassingProbe_OHmissing_eta_1D")->second->Fill(GEM_eta);
1840 if (~GEM_warning >> 1 & 1) {
1841 m_histos.find("GEM_nPassingProbe_OHerror_eta_1D")->second->Fill(GEM_eta);
1842 if (~GEM_warning >> 2 & 1) {
1843 m_histos.find("GEM_nPassingProbe_VFATMask_eta_1D")->second->Fill(GEM_eta);
1844 }
1845 }
1846 }
1847 m_histos.find("GEM_nPassingProbe_phi_1D")->second->Fill(GEM_phi);
1848 if (~GEM_warning & 1) {
1849 m_histos.find("GEM_nPassingProbe_OHmissing_phi_1D")->second->Fill(GEM_phi);
1850 if (~GEM_warning >> 1 & 1) {
1851 m_histos.find("GEM_nPassingProbe_OHerror_phi_1D")->second->Fill(GEM_phi);
1852 if (~GEM_warning >> 2 & 1) {
1853 m_histos.find("GEM_nPassingProbe_VFATMask_phi_1D")->second->Fill(GEM_phi);
1854 }
1855 }
1856 }
1857 } else {
1858 m_histos.find("GE11_nFailingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
1859 m_histos.find("GE11_nFailingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
1860 m_histos.find("GE11_nFailingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
1861 m_histos.find("GE11_nFailingProbe_allCh_1D")->second->Fill(GEM_region);
1862 if (~GEM_warning & 1) {
1863 m_histos.find("GE11_nFailingProbe_OHmissing_allCh_1D")->second->Fill(GEM_region);
1864 if (~GEM_warning >> 1 & 1) {
1865 m_histos.find("GE11_nFailingProbe_OHerror_allCh_1D")->second->Fill(GEM_region);
1866 if (~GEM_warning >> 2 & 1) {
1867 m_histos.find("GE11_nFailingProbe_VFATMask_allCh_1D")->second->Fill(GEM_region);
1868 }
1869 }
1870 }
1871 m_histos.find("GE11_nFailingProbe_chamber_1D")->second->Fill(GEM_chamber);
1872 if (GEM_region < 0) {
1873 if (GEM_sta == 2 and GEM_lay == 2)
1874 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
1875 else if (GEM_sta == 2 and GEM_lay == 1)
1876 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
1877 else if (GEM_sta == 1 and GEM_lay == 2)
1878 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(2, GEM_chamber);
1879 else if (GEM_sta == 1 and GEM_lay == 1)
1880 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(3, GEM_chamber);
1881 }
1882 if (GEM_region > 0) {
1883 if (GEM_sta == 1 and GEM_lay == 1)
1884 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(6, GEM_chamber);
1885 else if (GEM_sta == 1 and GEM_lay == 2)
1886 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(7, GEM_chamber);
1887 else if (GEM_sta == 2 and GEM_lay == 1)
1888 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
1889 else if (GEM_sta == 2 and GEM_lay == 2)
1890 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
1891 }
1892 if (GEM_region == -1) {
1893 m_histos.find("GEM_nFailingProbe_Ch_region_GE1_NoL")->second->Fill(0, GEM_chamber);
1894 } else if (GEM_region == 1) {
1895 m_histos.find("GEM_nFailingProbe_Ch_region_GE1_NoL")->second->Fill(1, GEM_chamber);
1896 }
1897
1898 if (GEM_region == 1 && GEM_lay == 1) {
1899 m_histos.find("GEM_nFailingProbe_chamber_p1_1D")->second->Fill(GEM_chamber);
1900 if (~GEM_warning & 1) {
1901 m_histos.find("GEM_nFailingProbe_OHmissing_chamber_p1_1D")->second->Fill(GEM_chamber);
1902 if (~GEM_warning >> 1 & 1) {
1903 m_histos.find("GEM_nFailingProbe_OHerror_chamber_p1_1D")->second->Fill(GEM_chamber);
1904 if (~GEM_warning >> 2 & 1) {
1905 m_histos.find("GEM_nFailingProbe_VFATMask_chamber_p1_1D")->second->Fill(GEM_chamber);
1906 }
1907 }
1908 }
1909 m_histos.find("GEM_nFailingProbe_Ch_region_GE1")->second->Fill(2, GEM_chamber);
1910 m_histos.find("GEM_nFailingProbe_pt_p1_1D")->second->Fill(GEM_pt);
1911 m_histos.find("GEM_nFailingProbe_eta_p1_1D")->second->Fill(abs(GEM_eta));
1912 m_histos.find("GEM_nFailingProbe_phi_p1_1D")->second->Fill(GEM_phi);
1913 if (GEM_sta == 1) {
1914 m_histos.find("GE11_nFailingProbe_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1915 if (~GEM_warning & 1) {
1916 m_histos.find("GE11_nFailingProbe_OHmissing_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1917 if (~GEM_warning >> 1 & 1) {
1918 m_histos.find("GE11_nFailingProbe_OHerror_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1919 if (~GEM_warning >> 2 & 1) {
1920 m_histos.find("GE11_nFailingProbe_VFATMask_Ch_ieta_p1")->second->Fill(GEM_ieta, GEM_chamber);
1921 }
1922 }
1923 }
1924 }
1925 } else if (GEM_region == 1 && GEM_lay == 2) {
1926 m_histos.find("GEM_nFailingProbe_chamber_p2_1D")->second->Fill(GEM_chamber);
1927 if (~GEM_warning & 1) {
1928 m_histos.find("GEM_nFailingProbe_OHmissing_chamber_p2_1D")->second->Fill(GEM_chamber);
1929 if (~GEM_warning >> 1 & 1) {
1930 m_histos.find("GEM_nFailingProbe_OHerror_chamber_p2_1D")->second->Fill(GEM_chamber);
1931 if (~GEM_warning >> 2 & 1) {
1932 m_histos.find("GEM_nFailingProbe_VFATMask_chamber_p2_1D")->second->Fill(GEM_chamber);
1933 }
1934 }
1935 }
1936 m_histos.find("GEM_nFailingProbe_Ch_region_GE1")->second->Fill(3, GEM_chamber);
1937 m_histos.find("GEM_nFailingProbe_pt_p2_1D")->second->Fill(GEM_pt);
1938 m_histos.find("GEM_nFailingProbe_eta_p2_1D")->second->Fill(abs(GEM_eta));
1939 m_histos.find("GEM_nFailingProbe_phi_p2_1D")->second->Fill(GEM_phi);
1940 if (GEM_sta == 1) {
1941 m_histos.find("GE11_nFailingProbe_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1942 if (~GEM_warning & 1) {
1943 m_histos.find("GE11_nFailingProbe_OHmissing_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1944 if (~GEM_warning >> 1 & 1) {
1945 m_histos.find("GE11_nFailingProbe_OHerror_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1946 if (~GEM_warning >> 2 & 1) {
1947 m_histos.find("GE11_nFailingProbe_VFATMask_Ch_ieta_p2")->second->Fill(GEM_ieta, GEM_chamber);
1948 }
1949 }
1950 }
1951 }
1952 } else if (GEM_region == -1 && GEM_lay == 1) {
1953 m_histos.find("GEM_nFailingProbe_chamber_n1_1D")->second->Fill(GEM_chamber);
1954 if (~GEM_warning & 1) {
1955 m_histos.find("GEM_nFailingProbe_OHmissing_chamber_n1_1D")->second->Fill(GEM_chamber);
1956 if (~GEM_warning >> 1 & 1) {
1957 m_histos.find("GEM_nFailingProbe_OHerror_chamber_n1_1D")->second->Fill(GEM_chamber);
1958 if (~GEM_warning >> 2 & 1) {
1959 m_histos.find("GEM_nFailingProbe_VFATMask_chamber_n1_1D")->second->Fill(GEM_chamber);
1960 }
1961 }
1962 }
1963 m_histos.find("GEM_nFailingProbe_Ch_region_GE1")->second->Fill(1, GEM_chamber);
1964 m_histos.find("GEM_nFailingProbe_pt_n1_1D")->second->Fill(GEM_pt);
1965 m_histos.find("GEM_nFailingProbe_eta_n1_1D")->second->Fill(abs(GEM_eta));
1966 m_histos.find("GEM_nFailingProbe_phi_n1_1D")->second->Fill(GEM_phi);
1967 if (GEM_sta == 1) {
1968 m_histos.find("GE11_nFailingProbe_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1969 if (~GEM_warning & 1) {
1970 m_histos.find("GE11_nFailingProbe_OHmissing_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1971 if (~GEM_warning >> 1 & 1) {
1972 m_histos.find("GE11_nFailingProbe_OHerror_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1973 if (~GEM_warning >> 2 & 1) {
1974 m_histos.find("GE11_nFailingProbe_VFATMask_Ch_ieta_n1")->second->Fill(GEM_ieta, GEM_chamber);
1975 }
1976 }
1977 }
1978 }
1979 } else if (GEM_region == -1 && GEM_lay == 2) {
1980 m_histos.find("GEM_nFailingProbe_chamber_n2_1D")->second->Fill(GEM_chamber);
1981 if (~GEM_warning & 1) {
1982 m_histos.find("GEM_nFailingProbe_OHmissing_chamber_n2_1D")->second->Fill(GEM_chamber);
1983 if (~GEM_warning >> 1 & 1) {
1984 m_histos.find("GEM_nFailingProbe_OHerror_chamber_n2_1D")->second->Fill(GEM_chamber);
1985 if (~GEM_warning >> 2 & 1) {
1986 m_histos.find("GEM_nFailingProbe_VFATMask_chamber_n2_1D")->second->Fill(GEM_chamber);
1987 }
1988 }
1989 }
1990 m_histos.find("GEM_nFailingProbe_Ch_region_GE1")->second->Fill(0, GEM_chamber);
1991 m_histos.find("GEM_nFailingProbe_pt_n2_1D")->second->Fill(GEM_pt);
1992 m_histos.find("GEM_nFailingProbe_eta_n2_1D")->second->Fill(abs(GEM_eta));
1993 m_histos.find("GEM_nFailingProbe_phi_n2_1D")->second->Fill(GEM_phi);
1994 if (GEM_sta == 1) {
1995 m_histos.find("GE11_nFailingProbe_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1996 if (~GEM_warning & 1) {
1997 m_histos.find("GE11_nFailingProbe_OHmissing_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
1998 if (~GEM_warning >> 1 & 1) {
1999 m_histos.find("GE11_nFailingProbe_OHerror_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
2000 if (~GEM_warning >> 2 & 1) {
2001 m_histos.find("GE11_nFailingProbe_VFATMask_Ch_ieta_n2")->second->Fill(GEM_ieta, GEM_chamber);
2002 }
2003 }
2004 }
2005 }
2006 }
2007 m_histos.find("GEM_nFailingProbe_pt_1D")->second->Fill(GEM_pt);
2008 if (~GEM_warning & 1) {
2009 m_histos.find("GEM_nFailingProbe_OHmissing_pt_1D")->second->Fill(GEM_pt);
2010 if (~GEM_warning >> 1 & 1) {
2011 m_histos.find("GEM_nFailingProbe_OHerror_pt_1D")->second->Fill(GEM_pt);
2012 if (~GEM_warning >> 2 & 1) {
2013 m_histos.find("GEM_nFailingProbe_VFATMask_pt_1D")->second->Fill(GEM_pt);
2014 }
2015 }
2016 }
2017 m_histos.find("GEM_nFailingProbe_eta_1D")->second->Fill(abs(GEM_eta));
2018 if (~GEM_warning & 1) {
2019 m_histos.find("GEM_nFailingProbe_OHmissing_eta_1D")->second->Fill(GEM_eta);
2020 if (~GEM_warning >> 1 & 1) {
2021 m_histos.find("GEM_nFailingProbe_OHerror_eta_1D")->second->Fill(GEM_eta);
2022 if (~GEM_warning >> 2 & 1) {
2023 m_histos.find("GEM_nFailingProbe_VFATMask_eta_1D")->second->Fill(GEM_eta);
2024 }
2025 }
2026 }
2027 m_histos.find("GEM_nFailingProbe_phi_1D")->second->Fill(GEM_phi);
2028 if (~GEM_warning & 1) {
2029 m_histos.find("GEM_nFailingProbe_OHmissing_phi_1D")->second->Fill(GEM_phi);
2030 if (~GEM_warning >> 1 & 1) {
2031 m_histos.find("GEM_nFailingProbe_OHerror_phi_1D")->second->Fill(GEM_phi);
2032 if (~GEM_warning >> 2 & 1) {
2033 m_histos.find("GEM_nFailingProbe_VFATMask_phi_1D")->second->Fill(GEM_phi);
2034 }
2035 }
2036 }
2037 }
2038 }
2039
2040
2041 unsigned nGE21_matches = probe_coll_GE21_region.at(i).size();
2042 for (unsigned j = 0; j < nGE21_matches; ++j) {
2043
2044 int GEM_region = probe_coll_GE21_region.at(i).at(j);
2045 int GEM_lay = probe_coll_GE21_lay.at(i).at(j);
2046 int GEM_chamber = probe_coll_GE21_chamber.at(i).at(j);
2047 float GEM_ieta = probe_coll_GE21_ieta.at(i).at(j);
2048 float GEM_dx = probe_coll_GE21_dx.at(i).at(j);
2049 float GEM_phi = probe_coll_GE21_phi.at(i).at(j);
2050
2051
2052 if (GEM_dx < m_dxCut) {
2053 m_histos.find("GE21_nPassingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
2054 m_histos.find("GE21_nPassingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
2055 m_histos.find("GE21_nPassingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
2056 m_histos.find("GE21_nPassingProbe_allCh_1D")->second->Fill(GEM_region);
2057 m_histos.find("GE21_nPassingProbe_chamber_1D")->second->Fill(GEM_chamber);
2058 if (GEM_region < 0) {
2059 if (GEM_lay == 2)
2060 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
2061 else if (GEM_lay == 1)
2062 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
2063 }
2064 if (GEM_region > 0) {
2065 if (GEM_lay == 1)
2066 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
2067 else if (GEM_lay == 2)
2068 m_histos.find("GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
2069 }
2070 } else {
2071 m_histos.find("GE21_nFailingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
2072 m_histos.find("GE21_nFailingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
2073 m_histos.find("GE21_nFailingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
2074 m_histos.find("GE21_nFailingProbe_allCh_1D")->second->Fill(GEM_region);
2075 m_histos.find("GE21_nFailingProbe_chamber_1D")->second->Fill(GEM_chamber);
2076 if (GEM_region < 0) {
2077 if (GEM_lay == 2)
2078 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
2079 else if (GEM_lay == 1)
2080 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
2081 }
2082 if (GEM_region > 0) {
2083 if (GEM_lay == 1)
2084 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
2085 else if (GEM_lay == 2)
2086 m_histos.find("GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
2087 }
2088 }
2089 }
2090 }
2091 }
2092
2093 std::string GEMTnPEfficiencyTask::topFolder() const { return "GEM/Segment_TnP/"; };
2094
2095 DEFINE_FWK_MODULE(GEMTnPEfficiencyTask);