Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-10-30 00:11:33

0001 /*
0002  * \file GEMTnPEfficiencyTask.cc
0003  * \author Qianying
0004  *
0005  * \interited from the TnP framework of  
0006  * \author L. Lunerti - INFN Bologna
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   /// Constructor
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   /// Destructor
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   /// Analyze
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   // Bad VFAT Mask + Empty OH Mask
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   // Bad VFAT Mask + Empty OH Mask + EvtNF
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   // Bad VFAT Mask + Empty OH Mask + OHerrorMask
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       //oh_warning = oh_warning | (oh_status->warnings()); // If doing oh warning masking
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))) {  // will not work for GE21
1316           oh_warning = oh_warning | (1 << 2);
1317         }
1318       }
1319     }  // range
1320   }  // collection
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     //if (oh_status_collem_tion.isValid()) {
1355     oh_status = *oh_status_collection;
1356   } else {
1357     LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask") << "failed to get GEMOHStatusCollection" << std::endl;
1358     return;
1359   }
1360 
1361   // event.getByToken(kGEMVFATStatusCollectionToken_, vfat_status_collection);
1362   // if (vfat_status_collection.isValid()) {
1363   //   vfat_status = *vfat_status_collection;
1364   // } else {
1365   //   LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask") << "failed to get GEMVFATStatusCollection" << std::endl;
1366   //   return;
1367   // }
1368   // }
1369   edm::Handle<reco::MuonCollection> muons;
1370   event.getByToken(m_muToken, muons);
1371 
1372   //GE11 variables
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   //GE21 variables
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;  // ME0 to 0b0001, GE11 to 0b0010, GE21 to 0b0100
1397 
1398   //ME0 variables
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   //Fill probe dx + subdetector coordinates
1416   for (const auto i : probe_indices) {
1417     //GE11 variables
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     //GE21 variables
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     //std::vector<float> probe_GEM_dx_seg;
1440     uint8_t GEM_stationMatching = 0;
1441     //ME0 variables
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;  // fill detailed plots only for probes matching GEM
1455 
1456     for (const auto& chambMatch : (*muons).at(i).matches()) {
1457       // look in GEMs
1458       bool hit_matched = false;  // true if chambermatch has at least one hit (GE11, GE21) or segment (ME0)
1459       if (chambMatch.detector() == MuonSubdetId::GEM) {
1460         if (chambMatch.edgeX < m_borderCut && chambMatch.edgeY < m_borderCut) {
1461           gem_matched = true;  //fill detailed plots if at least one GEM probe match
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     }  //loop over chamber matches
1579 
1580     //Fill detailed plots
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       //for(int ii=0; i<probe_GEM_dx.size(); ii++)
1587       //{
1588       //    m_histos.find("GEMhit_dx")->second->Fill(probe_GEM_dx[ii]);
1589       //    m_histos.find("GEMseg_dx")->second->Fill(probe_GEM_dx_seg[ii]);
1590       //}
1591     }
1592 
1593     //Fill GEM variables
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     //Fill GE21 variables
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     //Fill ME0 variables
1620     probe_coll_ME0_region.push_back(probe_ME0_region);
1621     probe_coll_ME0_roll.push_back(probe_ME0_roll);  // same as ieta
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   }  //loop over probe collection
1633 
1634   //Loop over probes
1635   for (unsigned i = 0; i < probe_indices.size(); ++i) {
1636     //uint8_t GEM_matchPatt = probe_coll_GEM_staMatch.at(i);  // ME0 to 0b0001, GE11 to 0b0010, GE21 to 0b0100
1637 
1638     //Loop over ME0 matches
1639     unsigned nME0_matches = probe_coll_ME0_region.at(i).size();
1640     for (unsigned j = 0; j < nME0_matches; ++j) {
1641       //ME0 variables
1642       int ME0_region = probe_coll_ME0_region.at(i).at(j);
1643       //int ME0_roll   = probe_coll_ME0_roll.at(i).at(j);
1644       //int ME0_sta = probe_coll_ME0_sta.at(i).at(j);
1645       //int ME0_lay    = probe_coll_ME0_lay.at(i).at(j);
1646       int ME0_chamber = probe_coll_ME0_chamber.at(i).at(j);
1647       //float ME0_pt   = probe_coll_ME0_pt.at(i).at(j);
1648       float ME0_dx = probe_coll_ME0_dx.at(i).at(j);
1649       //float ME0_eta   = probe_coll_ME0_eta.at(i).at(j);
1650       //float ME0_phi   = probe_coll_ME0_phi.at(i).at(j);
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     //Loop over GE11 matches
1669     unsigned nGE11_matches = probe_coll_GE11_region.at(i).size();
1670     for (unsigned j = 0; j < nGE11_matches; ++j) {
1671       //GEM variables
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       //Fill GEM plots
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     //Loop over GE21 matches
2041     unsigned nGE21_matches = probe_coll_GE21_region.at(i).size();
2042     for (unsigned j = 0; j < nGE21_matches; ++j) {
2043       //GEM variables
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       //Fill GEM plots
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);