Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-11-18 03:09:40

0001 #include <string>
0002 #include <numeric>
0003 #include <vector>
0004 #include <algorithm>
0005 #include <iostream>
0006 #include <memory>
0007 
0008 #include "FWCore/Framework/interface/Frameworkfwd.h"
0009 #include "FWCore/Framework/interface/Event.h"
0010 #include "FWCore/Framework/interface/MakerMacros.h"
0011 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0012 #include "DataFormats/Common/interface/ValueMap.h"
0013 
0014 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
0015 #include "FWCore/ParameterSet/interface/PluginDescription.h"
0016 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0017 #include "FWCore/Utilities/interface/ESGetToken.h"
0018 #include "FWCore/Framework/interface/ESHandle.h"
0019 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
0020 #include "FWCore/ServiceRegistry/interface/Service.h"
0021 
0022 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
0023 #include "DQMServices/Core/interface/DQMStore.h"
0024 
0025 #include "DataFormats/Common/interface/ValidHandle.h"
0026 #include "DataFormats/Math/interface/deltaR.h"
0027 #include "DataFormats/Math/interface/GeantUnits.h"
0028 #include "DataFormats/ForwardDetId/interface/ETLDetId.h"
0029 #include "DataFormats/ForwardDetId/interface/BTLDetId.h"
0030 
0031 #include "DataFormats/Common/interface/Ptr.h"
0032 #include "DataFormats/Common/interface/PtrVector.h"
0033 #include "DataFormats/Common/interface/RefProd.h"
0034 #include "DataFormats/Common/interface/Ref.h"
0035 #include "DataFormats/Common/interface/RefVector.h"
0036 
0037 #include "DataFormats/TrackReco/interface/Track.h"
0038 #include "DataFormats/VertexReco/interface/Vertex.h"
0039 #include "DataFormats/TrackReco/interface/TrackFwd.h"
0040 
0041 #include "Geometry/Records/interface/MTDDigiGeometryRecord.h"
0042 #include "Geometry/MTDGeometryBuilder/interface/MTDGeometry.h"
0043 
0044 #include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
0045 #include "SimDataFormats/Associations/interface/TrackToTrackingParticleAssociator.h"
0046 #include "SimDataFormats/TrackingAnalysis/interface/TrackingVertexContainer.h"
0047 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
0048 
0049 // Adding header files for electrons
0050 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
0051 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
0052 
0053 // eff vs PU test libraries
0054 #include "DataFormats/VertexReco/interface/VertexFwd.h"
0055 #include "RecoVertex/VertexPrimitives/interface/TransientVertex.h"
0056 #include "SimDataFormats/PileupSummaryInfo/interface/PileupSummaryInfo.h"
0057 #include "SimTracker/VertexAssociation/interface/calculateVertexSharedTracks.h"
0058 #include "RecoVertex/PrimaryVertexProducer/interface/TrackFilterForPVFinding.h"
0059 #include "SimDataFormats/Associations/interface/VertexToTrackingVertexAssociator.h"
0060 #include "DataFormats/Math/interface/LorentzVector.h"
0061 #include "DataFormats/Math/interface/Point3D.h"
0062 
0063 class MtdEleIsoValidation : public DQMEDAnalyzer {
0064 public:
0065   explicit MtdEleIsoValidation(const edm::ParameterSet&);
0066   ~MtdEleIsoValidation() override;
0067 
0068   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0069 
0070 private:
0071   void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
0072 
0073   void analyze(const edm::Event&, const edm::EventSetup&) override;
0074 
0075   // ------------ member data ------------
0076 
0077   const std::string folder_;
0078   const float trackMinPt_;
0079   const float trackMinEta_;
0080   const float trackMaxEta_;
0081   const double rel_iso_cut_;
0082 
0083   const bool track_match_PV_;
0084   const bool dt_sig_track_;
0085   const bool optionalPlots_;
0086 
0087   const float min_dR_cut;
0088   const float max_dR_cut;
0089   const float min_pt_cut_EB;
0090   const float min_pt_cut_EE;
0091   const float max_dz_cut_EB;
0092   const float max_dz_cut_EE;
0093   const float max_dz_vtx_cut;
0094   const float max_dxy_vtx_cut;
0095   const float min_strip_cut;
0096   const float min_track_mtd_mva_cut;
0097   const std::vector<double> max_dt_vtx_cut{0.30, 0.27, 0.24, 0.21, 0.18, 0.15, 0.12};
0098   const std::vector<double> max_dt_track_cut{0.30, 0.27, 0.24, 0.21, 0.18, 0.15, 0.12};
0099   const std::vector<double> max_dt_significance_cut{4.0, 3.0, 2.0};
0100   const std::vector<double> pT_bins_dt_distrb{10, 20, 30, 40, 50, 60, 70, 80, 90, 100};
0101   const std::vector<double> eta_bins_dt_distrib{0.0, 0.5, 1.0, 1.5, 2.0, 2.4, 2.7, 3};
0102   static constexpr double avg_sim_sigTrk_t_err = 0.03239;
0103   static constexpr double avg_sim_PUtrack_t_err = 0.03465;
0104 
0105   edm::EDGetTokenT<reco::TrackCollection> GenRecTrackToken_;
0106   edm::EDGetTokenT<std::vector<reco::Vertex>> RecVertexToken_;
0107 
0108   edm::EDGetTokenT<reco::GsfElectronCollection> GsfElectronToken_EB_;
0109   edm::EDGetTokenT<reco::GsfElectronCollection> GsfElectronToken_EE_;
0110   edm::EDGetTokenT<reco::GenParticleCollection> GenParticleToken_;
0111 
0112   edm::EDGetTokenT<edm::ValueMap<float>> t0PidToken_;
0113   edm::EDGetTokenT<edm::ValueMap<float>> Sigmat0PidToken_;
0114   edm::EDGetTokenT<edm::ValueMap<float>> trackMVAQualToken_;
0115 
0116   edm::EDGetTokenT<reco::RecoToSimCollection> recoToSimAssociationToken_;
0117 
0118   // Signal histograms
0119 
0120   MonitorElement* meEle_no_dt_check_;
0121   MonitorElement* meTrk_genMatch_check_;
0122 
0123   MonitorElement* meEle_avg_error_SigTrk_check_;
0124   MonitorElement* meEle_avg_error_PUTrk_check_;
0125   MonitorElement* meEle_avg_error_vtx_check_;
0126 
0127   // Adding histograms for barrel electrons
0128   MonitorElement* meEleISO_Ntracks_Sig_EB_;
0129   MonitorElement* meEleISO_chIso_Sig_EB_;
0130   MonitorElement* meEleISO_rel_chIso_Sig_EB_;
0131   MonitorElement* meEleISO_Ntracks_MTD_1_Sig_EB_;
0132   MonitorElement* meEleISO_chIso_MTD_1_Sig_EB_;
0133   MonitorElement* meEleISO_rel_chIso_MTD_1_Sig_EB_;
0134 
0135   MonitorElement* meEleISO_Ntracks_MTD_2_Sig_EB_;
0136   MonitorElement* meEleISO_chIso_MTD_2_Sig_EB_;
0137   MonitorElement* meEleISO_rel_chIso_MTD_2_Sig_EB_;
0138 
0139   MonitorElement* meEleISO_Ntracks_MTD_3_Sig_EB_;
0140   MonitorElement* meEleISO_chIso_MTD_3_Sig_EB_;
0141   MonitorElement* meEleISO_rel_chIso_MTD_3_Sig_EB_;
0142 
0143   MonitorElement* meEleISO_Ntracks_MTD_4_Sig_EB_;
0144   MonitorElement* meEleISO_chIso_MTD_4_Sig_EB_;
0145   MonitorElement* meEleISO_rel_chIso_MTD_4_Sig_EB_;
0146 
0147   MonitorElement* meEleISO_Ntracks_MTD_5_Sig_EB_;
0148   MonitorElement* meEleISO_chIso_MTD_5_Sig_EB_;
0149   MonitorElement* meEleISO_rel_chIso_MTD_5_Sig_EB_;
0150 
0151   MonitorElement* meEleISO_Ntracks_MTD_6_Sig_EB_;
0152   MonitorElement* meEleISO_chIso_MTD_6_Sig_EB_;
0153   MonitorElement* meEleISO_rel_chIso_MTD_6_Sig_EB_;
0154 
0155   MonitorElement* meEleISO_Ntracks_MTD_7_Sig_EB_;
0156   MonitorElement* meEleISO_chIso_MTD_7_Sig_EB_;
0157   MonitorElement* meEleISO_rel_chIso_MTD_7_Sig_EB_;
0158 
0159   MonitorElement* meEleISO_Ntracks_MTD_sim_1_Sig_EB_;
0160   MonitorElement* meEleISO_chIso_MTD_sim_1_Sig_EB_;
0161   MonitorElement* meEleISO_rel_chIso_MTD_sim_1_Sig_EB_;
0162 
0163   MonitorElement* meEleISO_Ntracks_MTD_sim_2_Sig_EB_;
0164   MonitorElement* meEleISO_chIso_MTD_sim_2_Sig_EB_;
0165   MonitorElement* meEleISO_rel_chIso_MTD_sim_2_Sig_EB_;
0166 
0167   MonitorElement* meEleISO_Ntracks_MTD_sim_3_Sig_EB_;
0168   MonitorElement* meEleISO_chIso_MTD_sim_3_Sig_EB_;
0169   MonitorElement* meEleISO_rel_chIso_MTD_sim_3_Sig_EB_;
0170 
0171   MonitorElement* meEleISO_Ntracks_MTD_sim_4_Sig_EB_;
0172   MonitorElement* meEleISO_chIso_MTD_sim_4_Sig_EB_;
0173   MonitorElement* meEleISO_rel_chIso_MTD_sim_4_Sig_EB_;
0174 
0175   MonitorElement* meEleISO_Ntracks_MTD_sim_5_Sig_EB_;
0176   MonitorElement* meEleISO_chIso_MTD_sim_5_Sig_EB_;
0177   MonitorElement* meEleISO_rel_chIso_MTD_sim_5_Sig_EB_;
0178 
0179   MonitorElement* meEleISO_Ntracks_MTD_sim_6_Sig_EB_;
0180   MonitorElement* meEleISO_chIso_MTD_sim_6_Sig_EB_;
0181   MonitorElement* meEleISO_rel_chIso_MTD_sim_6_Sig_EB_;
0182 
0183   MonitorElement* meEleISO_Ntracks_MTD_sim_7_Sig_EB_;
0184   MonitorElement* meEleISO_chIso_MTD_sim_7_Sig_EB_;
0185   MonitorElement* meEleISO_rel_chIso_MTD_sim_7_Sig_EB_;
0186 
0187   MonitorElement* meEleISO_Ntracks_gen_Sig_EB_;
0188   MonitorElement* meEleISO_chIso_gen_Sig_EB_;
0189   MonitorElement* meEleISO_rel_chIso_gen_Sig_EB_;
0190 
0191   MonitorElement* meEleISO_Ntracks_MTD_2sigma_Sig_EB_;
0192   MonitorElement* meEleISO_chIso_MTD_2sigma_Sig_EB_;
0193   MonitorElement* meEleISO_rel_chIso_MTD_2sigma_Sig_EB_;
0194 
0195   MonitorElement* meEleISO_Ntracks_MTD_3sigma_Sig_EB_;
0196   MonitorElement* meEleISO_chIso_MTD_3sigma_Sig_EB_;
0197   MonitorElement* meEleISO_rel_chIso_MTD_3sigma_Sig_EB_;
0198 
0199   MonitorElement* meEleISO_Ntracks_MTD_4sigma_Sig_EB_;
0200   MonitorElement* meEleISO_chIso_MTD_4sigma_Sig_EB_;
0201   MonitorElement* meEleISO_rel_chIso_MTD_4sigma_Sig_EB_;
0202 
0203   MonitorElement* meEleISO_Ntracks_MTD_sim_2sigma_Sig_EB_;
0204   MonitorElement* meEleISO_chIso_MTD_sim_2sigma_Sig_EB_;
0205   MonitorElement* meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EB_;
0206 
0207   MonitorElement* meEleISO_Ntracks_MTD_sim_3sigma_Sig_EB_;
0208   MonitorElement* meEleISO_chIso_MTD_sim_3sigma_Sig_EB_;
0209   MonitorElement* meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EB_;
0210 
0211   MonitorElement* meEleISO_Ntracks_MTD_sim_4sigma_Sig_EB_;
0212   MonitorElement* meEleISO_chIso_MTD_sim_4sigma_Sig_EB_;
0213   MonitorElement* meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EB_;
0214 
0215   MonitorElement* meEle_pt_tot_Sig_EB_;
0216   MonitorElement* meEle_pt_sim_tot_Sig_EB_;  // for GEN case is the same
0217   MonitorElement* meEle_eta_tot_Sig_EB_;
0218   MonitorElement* meEle_phi_tot_Sig_EB_;
0219   MonitorElement* meEle_pt_MTD_1_Sig_EB_;
0220   MonitorElement* meEle_pt_sim_MTD_1_Sig_EB_;
0221   MonitorElement* meEle_eta_MTD_1_Sig_EB_;
0222   MonitorElement* meEle_phi_MTD_1_Sig_EB_;
0223 
0224   MonitorElement* meEle_pt_MTD_2_Sig_EB_;
0225   MonitorElement* meEle_pt_sim_MTD_2_Sig_EB_;
0226   MonitorElement* meEle_eta_MTD_2_Sig_EB_;
0227   MonitorElement* meEle_phi_MTD_2_Sig_EB_;
0228 
0229   MonitorElement* meEle_pt_MTD_3_Sig_EB_;
0230   MonitorElement* meEle_pt_sim_MTD_3_Sig_EB_;
0231   MonitorElement* meEle_eta_MTD_3_Sig_EB_;
0232   MonitorElement* meEle_phi_MTD_3_Sig_EB_;
0233 
0234   MonitorElement* meEle_pt_MTD_4_Sig_EB_;
0235   MonitorElement* meEle_pt_sim_MTD_4_Sig_EB_;
0236   MonitorElement* meEle_eta_MTD_4_Sig_EB_;
0237   MonitorElement* meEle_phi_MTD_4_Sig_EB_;
0238 
0239   MonitorElement* meEle_pt_MTD_5_Sig_EB_;
0240   MonitorElement* meEle_pt_sim_MTD_5_Sig_EB_;
0241   MonitorElement* meEle_eta_MTD_5_Sig_EB_;
0242   MonitorElement* meEle_phi_MTD_5_Sig_EB_;
0243 
0244   MonitorElement* meEle_pt_MTD_6_Sig_EB_;
0245   MonitorElement* meEle_pt_sim_MTD_6_Sig_EB_;
0246   MonitorElement* meEle_eta_MTD_6_Sig_EB_;
0247   MonitorElement* meEle_phi_MTD_6_Sig_EB_;
0248 
0249   MonitorElement* meEle_pt_MTD_7_Sig_EB_;
0250   MonitorElement* meEle_pt_sim_MTD_7_Sig_EB_;
0251   MonitorElement* meEle_eta_MTD_7_Sig_EB_;
0252   MonitorElement* meEle_phi_MTD_7_Sig_EB_;
0253 
0254   MonitorElement* meEle_pt_noMTD_Sig_EB_;
0255   MonitorElement* meEle_eta_noMTD_Sig_EB_;
0256   MonitorElement* meEle_phi_noMTD_Sig_EB_;
0257 
0258   MonitorElement* meEle_pt_gen_Sig_EB_;
0259   MonitorElement* meEle_eta_gen_Sig_EB_;
0260   MonitorElement* meEle_phi_gen_Sig_EB_;
0261 
0262   MonitorElement* meEle_pt_MTD_2sigma_Sig_EB_;
0263   MonitorElement* meEle_pt_sim_MTD_2sigma_Sig_EB_;
0264   MonitorElement* meEle_eta_MTD_2sigma_Sig_EB_;
0265   MonitorElement* meEle_phi_MTD_2sigma_Sig_EB_;
0266 
0267   MonitorElement* meEle_pt_MTD_3sigma_Sig_EB_;
0268   MonitorElement* meEle_pt_sim_MTD_3sigma_Sig_EB_;
0269   MonitorElement* meEle_eta_MTD_3sigma_Sig_EB_;
0270   MonitorElement* meEle_phi_MTD_3sigma_Sig_EB_;
0271 
0272   MonitorElement* meEle_pt_MTD_4sigma_Sig_EB_;
0273   MonitorElement* meEle_pt_sim_MTD_4sigma_Sig_EB_;
0274   MonitorElement* meEle_eta_MTD_4sigma_Sig_EB_;
0275   MonitorElement* meEle_phi_MTD_4sigma_Sig_EB_;
0276 
0277   // Adding histograms for endcap electrons
0278   MonitorElement* meEleISO_Ntracks_Sig_EE_;
0279   MonitorElement* meEleISO_chIso_Sig_EE_;
0280   MonitorElement* meEleISO_rel_chIso_Sig_EE_;
0281   MonitorElement* meEleISO_Ntracks_MTD_1_Sig_EE_;
0282   MonitorElement* meEleISO_chIso_MTD_1_Sig_EE_;
0283   MonitorElement* meEleISO_rel_chIso_MTD_1_Sig_EE_;
0284 
0285   MonitorElement* meEleISO_Ntracks_MTD_2_Sig_EE_;
0286   MonitorElement* meEleISO_chIso_MTD_2_Sig_EE_;
0287   MonitorElement* meEleISO_rel_chIso_MTD_2_Sig_EE_;
0288 
0289   MonitorElement* meEleISO_Ntracks_MTD_3_Sig_EE_;
0290   MonitorElement* meEleISO_chIso_MTD_3_Sig_EE_;
0291   MonitorElement* meEleISO_rel_chIso_MTD_3_Sig_EE_;
0292 
0293   MonitorElement* meEleISO_Ntracks_MTD_4_Sig_EE_;
0294   MonitorElement* meEleISO_chIso_MTD_4_Sig_EE_;
0295   MonitorElement* meEleISO_rel_chIso_MTD_4_Sig_EE_;
0296 
0297   MonitorElement* meEleISO_Ntracks_MTD_5_Sig_EE_;
0298   MonitorElement* meEleISO_chIso_MTD_5_Sig_EE_;
0299   MonitorElement* meEleISO_rel_chIso_MTD_5_Sig_EE_;
0300 
0301   MonitorElement* meEleISO_Ntracks_MTD_6_Sig_EE_;
0302   MonitorElement* meEleISO_chIso_MTD_6_Sig_EE_;
0303   MonitorElement* meEleISO_rel_chIso_MTD_6_Sig_EE_;
0304 
0305   MonitorElement* meEleISO_Ntracks_MTD_7_Sig_EE_;
0306   MonitorElement* meEleISO_chIso_MTD_7_Sig_EE_;
0307   MonitorElement* meEleISO_rel_chIso_MTD_7_Sig_EE_;
0308 
0309   MonitorElement* meEleISO_Ntracks_MTD_sim_1_Sig_EE_;
0310   MonitorElement* meEleISO_chIso_MTD_sim_1_Sig_EE_;
0311   MonitorElement* meEleISO_rel_chIso_MTD_sim_1_Sig_EE_;
0312 
0313   MonitorElement* meEleISO_Ntracks_MTD_sim_2_Sig_EE_;
0314   MonitorElement* meEleISO_chIso_MTD_sim_2_Sig_EE_;
0315   MonitorElement* meEleISO_rel_chIso_MTD_sim_2_Sig_EE_;
0316 
0317   MonitorElement* meEleISO_Ntracks_MTD_sim_3_Sig_EE_;
0318   MonitorElement* meEleISO_chIso_MTD_sim_3_Sig_EE_;
0319   MonitorElement* meEleISO_rel_chIso_MTD_sim_3_Sig_EE_;
0320 
0321   MonitorElement* meEleISO_Ntracks_MTD_sim_4_Sig_EE_;
0322   MonitorElement* meEleISO_chIso_MTD_sim_4_Sig_EE_;
0323   MonitorElement* meEleISO_rel_chIso_MTD_sim_4_Sig_EE_;
0324 
0325   MonitorElement* meEleISO_Ntracks_MTD_sim_5_Sig_EE_;
0326   MonitorElement* meEleISO_chIso_MTD_sim_5_Sig_EE_;
0327   MonitorElement* meEleISO_rel_chIso_MTD_sim_5_Sig_EE_;
0328 
0329   MonitorElement* meEleISO_Ntracks_MTD_sim_6_Sig_EE_;
0330   MonitorElement* meEleISO_chIso_MTD_sim_6_Sig_EE_;
0331   MonitorElement* meEleISO_rel_chIso_MTD_sim_6_Sig_EE_;
0332 
0333   MonitorElement* meEleISO_Ntracks_MTD_sim_7_Sig_EE_;
0334   MonitorElement* meEleISO_chIso_MTD_sim_7_Sig_EE_;
0335   MonitorElement* meEleISO_rel_chIso_MTD_sim_7_Sig_EE_;
0336 
0337   MonitorElement* meEleISO_Ntracks_gen_Sig_EE_;
0338   MonitorElement* meEleISO_chIso_gen_Sig_EE_;
0339   MonitorElement* meEleISO_rel_chIso_gen_Sig_EE_;
0340 
0341   MonitorElement* meEleISO_Ntracks_MTD_2sigma_Sig_EE_;
0342   MonitorElement* meEleISO_chIso_MTD_2sigma_Sig_EE_;
0343   MonitorElement* meEleISO_rel_chIso_MTD_2sigma_Sig_EE_;
0344 
0345   MonitorElement* meEleISO_Ntracks_MTD_3sigma_Sig_EE_;
0346   MonitorElement* meEleISO_chIso_MTD_3sigma_Sig_EE_;
0347   MonitorElement* meEleISO_rel_chIso_MTD_3sigma_Sig_EE_;
0348 
0349   MonitorElement* meEleISO_Ntracks_MTD_4sigma_Sig_EE_;
0350   MonitorElement* meEleISO_chIso_MTD_4sigma_Sig_EE_;
0351   MonitorElement* meEleISO_rel_chIso_MTD_4sigma_Sig_EE_;
0352 
0353   MonitorElement* meEleISO_Ntracks_MTD_sim_2sigma_Sig_EE_;
0354   MonitorElement* meEleISO_chIso_MTD_sim_2sigma_Sig_EE_;
0355   MonitorElement* meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EE_;
0356 
0357   MonitorElement* meEleISO_Ntracks_MTD_sim_3sigma_Sig_EE_;
0358   MonitorElement* meEleISO_chIso_MTD_sim_3sigma_Sig_EE_;
0359   MonitorElement* meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EE_;
0360 
0361   MonitorElement* meEleISO_Ntracks_MTD_sim_4sigma_Sig_EE_;
0362   MonitorElement* meEleISO_chIso_MTD_sim_4sigma_Sig_EE_;
0363   MonitorElement* meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EE_;
0364 
0365   MonitorElement* meEle_pt_tot_Sig_EE_;
0366   MonitorElement* meEle_pt_sim_tot_Sig_EE_;
0367   MonitorElement* meEle_eta_tot_Sig_EE_;
0368   MonitorElement* meEle_phi_tot_Sig_EE_;
0369 
0370   MonitorElement* meEle_pt_MTD_1_Sig_EE_;
0371   MonitorElement* meEle_pt_sim_MTD_1_Sig_EE_;
0372   MonitorElement* meEle_eta_MTD_1_Sig_EE_;
0373   MonitorElement* meEle_phi_MTD_1_Sig_EE_;
0374 
0375   MonitorElement* meEle_pt_MTD_2_Sig_EE_;
0376   MonitorElement* meEle_pt_sim_MTD_2_Sig_EE_;
0377   MonitorElement* meEle_eta_MTD_2_Sig_EE_;
0378   MonitorElement* meEle_phi_MTD_2_Sig_EE_;
0379 
0380   MonitorElement* meEle_pt_MTD_3_Sig_EE_;
0381   MonitorElement* meEle_pt_sim_MTD_3_Sig_EE_;
0382   MonitorElement* meEle_eta_MTD_3_Sig_EE_;
0383   MonitorElement* meEle_phi_MTD_3_Sig_EE_;
0384 
0385   MonitorElement* meEle_pt_MTD_4_Sig_EE_;
0386   MonitorElement* meEle_pt_sim_MTD_4_Sig_EE_;
0387   MonitorElement* meEle_eta_MTD_4_Sig_EE_;
0388   MonitorElement* meEle_phi_MTD_4_Sig_EE_;
0389 
0390   MonitorElement* meEle_pt_MTD_5_Sig_EE_;
0391   MonitorElement* meEle_pt_sim_MTD_5_Sig_EE_;
0392   MonitorElement* meEle_eta_MTD_5_Sig_EE_;
0393   MonitorElement* meEle_phi_MTD_5_Sig_EE_;
0394 
0395   MonitorElement* meEle_pt_MTD_6_Sig_EE_;
0396   MonitorElement* meEle_pt_sim_MTD_6_Sig_EE_;
0397   MonitorElement* meEle_eta_MTD_6_Sig_EE_;
0398   MonitorElement* meEle_phi_MTD_6_Sig_EE_;
0399 
0400   MonitorElement* meEle_pt_MTD_7_Sig_EE_;
0401   MonitorElement* meEle_pt_sim_MTD_7_Sig_EE_;
0402   MonitorElement* meEle_eta_MTD_7_Sig_EE_;
0403   MonitorElement* meEle_phi_MTD_7_Sig_EE_;
0404 
0405   MonitorElement* meEle_pt_noMTD_Sig_EE_;
0406   MonitorElement* meEle_eta_noMTD_Sig_EE_;
0407   MonitorElement* meEle_phi_noMTD_Sig_EE_;
0408 
0409   MonitorElement* meEle_pt_gen_Sig_EE_;
0410   MonitorElement* meEle_eta_gen_Sig_EE_;
0411   MonitorElement* meEle_phi_gen_Sig_EE_;
0412 
0413   MonitorElement* meEle_pt_MTD_2sigma_Sig_EE_;
0414   MonitorElement* meEle_pt_sim_MTD_2sigma_Sig_EE_;
0415   MonitorElement* meEle_eta_MTD_2sigma_Sig_EE_;
0416   MonitorElement* meEle_phi_MTD_2sigma_Sig_EE_;
0417 
0418   MonitorElement* meEle_pt_MTD_3sigma_Sig_EE_;
0419   MonitorElement* meEle_pt_sim_MTD_3sigma_Sig_EE_;
0420   MonitorElement* meEle_eta_MTD_3sigma_Sig_EE_;
0421   MonitorElement* meEle_phi_MTD_3sigma_Sig_EE_;
0422 
0423   MonitorElement* meEle_pt_MTD_4sigma_Sig_EE_;
0424   MonitorElement* meEle_pt_sim_MTD_4sigma_Sig_EE_;
0425   MonitorElement* meEle_eta_MTD_4sigma_Sig_EE_;
0426   MonitorElement* meEle_phi_MTD_4sigma_Sig_EE_;
0427 
0428   // Signal histograms end
0429 
0430   // Background histograms
0431   // Adding histograms for barrel electrons
0432   MonitorElement* meEleISO_Ntracks_Bkg_EB_;
0433   MonitorElement* meEleISO_chIso_Bkg_EB_;
0434   MonitorElement* meEleISO_rel_chIso_Bkg_EB_;
0435   MonitorElement* meEleISO_Ntracks_MTD_1_Bkg_EB_;
0436   MonitorElement* meEleISO_chIso_MTD_1_Bkg_EB_;
0437   MonitorElement* meEleISO_rel_chIso_MTD_1_Bkg_EB_;
0438 
0439   MonitorElement* meEleISO_Ntracks_MTD_2_Bkg_EB_;
0440   MonitorElement* meEleISO_chIso_MTD_2_Bkg_EB_;
0441   MonitorElement* meEleISO_rel_chIso_MTD_2_Bkg_EB_;
0442 
0443   MonitorElement* meEleISO_Ntracks_MTD_3_Bkg_EB_;
0444   MonitorElement* meEleISO_chIso_MTD_3_Bkg_EB_;
0445   MonitorElement* meEleISO_rel_chIso_MTD_3_Bkg_EB_;
0446 
0447   MonitorElement* meEleISO_Ntracks_MTD_4_Bkg_EB_;
0448   MonitorElement* meEleISO_chIso_MTD_4_Bkg_EB_;
0449   MonitorElement* meEleISO_rel_chIso_MTD_4_Bkg_EB_;
0450 
0451   MonitorElement* meEleISO_Ntracks_MTD_5_Bkg_EB_;
0452   MonitorElement* meEleISO_chIso_MTD_5_Bkg_EB_;
0453   MonitorElement* meEleISO_rel_chIso_MTD_5_Bkg_EB_;
0454 
0455   MonitorElement* meEleISO_Ntracks_MTD_6_Bkg_EB_;
0456   MonitorElement* meEleISO_chIso_MTD_6_Bkg_EB_;
0457   MonitorElement* meEleISO_rel_chIso_MTD_6_Bkg_EB_;
0458 
0459   MonitorElement* meEleISO_Ntracks_MTD_7_Bkg_EB_;
0460   MonitorElement* meEleISO_chIso_MTD_7_Bkg_EB_;
0461   MonitorElement* meEleISO_rel_chIso_MTD_7_Bkg_EB_;
0462 
0463   MonitorElement* meEleISO_Ntracks_MTD_sim_1_Bkg_EB_;
0464   MonitorElement* meEleISO_chIso_MTD_sim_1_Bkg_EB_;
0465   MonitorElement* meEleISO_rel_chIso_MTD_sim_1_Bkg_EB_;
0466 
0467   MonitorElement* meEleISO_Ntracks_MTD_sim_2_Bkg_EB_;
0468   MonitorElement* meEleISO_chIso_MTD_sim_2_Bkg_EB_;
0469   MonitorElement* meEleISO_rel_chIso_MTD_sim_2_Bkg_EB_;
0470 
0471   MonitorElement* meEleISO_Ntracks_MTD_sim_3_Bkg_EB_;
0472   MonitorElement* meEleISO_chIso_MTD_sim_3_Bkg_EB_;
0473   MonitorElement* meEleISO_rel_chIso_MTD_sim_3_Bkg_EB_;
0474 
0475   MonitorElement* meEleISO_Ntracks_MTD_sim_4_Bkg_EB_;
0476   MonitorElement* meEleISO_chIso_MTD_sim_4_Bkg_EB_;
0477   MonitorElement* meEleISO_rel_chIso_MTD_sim_4_Bkg_EB_;
0478 
0479   MonitorElement* meEleISO_Ntracks_MTD_sim_5_Bkg_EB_;
0480   MonitorElement* meEleISO_chIso_MTD_sim_5_Bkg_EB_;
0481   MonitorElement* meEleISO_rel_chIso_MTD_sim_5_Bkg_EB_;
0482 
0483   MonitorElement* meEleISO_Ntracks_MTD_sim_6_Bkg_EB_;
0484   MonitorElement* meEleISO_chIso_MTD_sim_6_Bkg_EB_;
0485   MonitorElement* meEleISO_rel_chIso_MTD_sim_6_Bkg_EB_;
0486 
0487   MonitorElement* meEleISO_Ntracks_MTD_sim_7_Bkg_EB_;
0488   MonitorElement* meEleISO_chIso_MTD_sim_7_Bkg_EB_;
0489   MonitorElement* meEleISO_rel_chIso_MTD_sim_7_Bkg_EB_;
0490 
0491   MonitorElement* meEleISO_Ntracks_gen_Bkg_EB_;
0492   MonitorElement* meEleISO_chIso_gen_Bkg_EB_;
0493   MonitorElement* meEleISO_rel_chIso_gen_Bkg_EB_;
0494 
0495   MonitorElement* meEleISO_Ntracks_MTD_2sigma_Bkg_EB_;
0496   MonitorElement* meEleISO_chIso_MTD_2sigma_Bkg_EB_;
0497   MonitorElement* meEleISO_rel_chIso_MTD_2sigma_Bkg_EB_;
0498 
0499   MonitorElement* meEleISO_Ntracks_MTD_3sigma_Bkg_EB_;
0500   MonitorElement* meEleISO_chIso_MTD_3sigma_Bkg_EB_;
0501   MonitorElement* meEleISO_rel_chIso_MTD_3sigma_Bkg_EB_;
0502 
0503   MonitorElement* meEleISO_Ntracks_MTD_4sigma_Bkg_EB_;
0504   MonitorElement* meEleISO_chIso_MTD_4sigma_Bkg_EB_;
0505   MonitorElement* meEleISO_rel_chIso_MTD_4sigma_Bkg_EB_;
0506 
0507   MonitorElement* meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EB_;
0508   MonitorElement* meEleISO_chIso_MTD_sim_2sigma_Bkg_EB_;
0509   MonitorElement* meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EB_;
0510 
0511   MonitorElement* meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EB_;
0512   MonitorElement* meEleISO_chIso_MTD_sim_3sigma_Bkg_EB_;
0513   MonitorElement* meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EB_;
0514 
0515   MonitorElement* meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EB_;
0516   MonitorElement* meEleISO_chIso_MTD_sim_4sigma_Bkg_EB_;
0517   MonitorElement* meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EB_;
0518 
0519   MonitorElement* meEle_pt_tot_Bkg_EB_;
0520   MonitorElement* meEle_pt_sim_tot_Bkg_EB_;
0521   MonitorElement* meEle_eta_tot_Bkg_EB_;
0522   MonitorElement* meEle_phi_tot_Bkg_EB_;
0523   MonitorElement* meEle_pt_MTD_1_Bkg_EB_;
0524   MonitorElement* meEle_pt_sim_MTD_1_Bkg_EB_;
0525   MonitorElement* meEle_eta_MTD_1_Bkg_EB_;
0526   MonitorElement* meEle_phi_MTD_1_Bkg_EB_;
0527 
0528   MonitorElement* meEle_pt_MTD_2_Bkg_EB_;
0529   MonitorElement* meEle_pt_sim_MTD_2_Bkg_EB_;
0530   MonitorElement* meEle_eta_MTD_2_Bkg_EB_;
0531   MonitorElement* meEle_phi_MTD_2_Bkg_EB_;
0532 
0533   MonitorElement* meEle_pt_MTD_3_Bkg_EB_;
0534   MonitorElement* meEle_pt_sim_MTD_3_Bkg_EB_;
0535   MonitorElement* meEle_eta_MTD_3_Bkg_EB_;
0536   MonitorElement* meEle_phi_MTD_3_Bkg_EB_;
0537 
0538   MonitorElement* meEle_pt_MTD_4_Bkg_EB_;
0539   MonitorElement* meEle_pt_sim_MTD_4_Bkg_EB_;
0540   MonitorElement* meEle_eta_MTD_4_Bkg_EB_;
0541   MonitorElement* meEle_phi_MTD_4_Bkg_EB_;
0542 
0543   MonitorElement* meEle_pt_MTD_5_Bkg_EB_;
0544   MonitorElement* meEle_pt_sim_MTD_5_Bkg_EB_;
0545   MonitorElement* meEle_eta_MTD_5_Bkg_EB_;
0546   MonitorElement* meEle_phi_MTD_5_Bkg_EB_;
0547 
0548   MonitorElement* meEle_pt_MTD_6_Bkg_EB_;
0549   MonitorElement* meEle_pt_sim_MTD_6_Bkg_EB_;
0550   MonitorElement* meEle_eta_MTD_6_Bkg_EB_;
0551   MonitorElement* meEle_phi_MTD_6_Bkg_EB_;
0552 
0553   MonitorElement* meEle_pt_MTD_7_Bkg_EB_;
0554   MonitorElement* meEle_pt_sim_MTD_7_Bkg_EB_;
0555   MonitorElement* meEle_eta_MTD_7_Bkg_EB_;
0556   MonitorElement* meEle_phi_MTD_7_Bkg_EB_;
0557 
0558   MonitorElement* meEle_pt_noMTD_Bkg_EB_;
0559   MonitorElement* meEle_eta_noMTD_Bkg_EB_;
0560   MonitorElement* meEle_phi_noMTD_Bkg_EB_;
0561 
0562   MonitorElement* meEle_pt_gen_Bkg_EB_;
0563   MonitorElement* meEle_eta_gen_Bkg_EB_;
0564   MonitorElement* meEle_phi_gen_Bkg_EB_;
0565 
0566   MonitorElement* meEle_pt_MTD_2sigma_Bkg_EB_;
0567   MonitorElement* meEle_pt_sim_MTD_2sigma_Bkg_EB_;
0568   MonitorElement* meEle_eta_MTD_2sigma_Bkg_EB_;
0569   MonitorElement* meEle_phi_MTD_2sigma_Bkg_EB_;
0570 
0571   MonitorElement* meEle_pt_MTD_3sigma_Bkg_EB_;
0572   MonitorElement* meEle_pt_sim_MTD_3sigma_Bkg_EB_;
0573   MonitorElement* meEle_eta_MTD_3sigma_Bkg_EB_;
0574   MonitorElement* meEle_phi_MTD_3sigma_Bkg_EB_;
0575 
0576   MonitorElement* meEle_pt_MTD_4sigma_Bkg_EB_;
0577   MonitorElement* meEle_pt_sim_MTD_4sigma_Bkg_EB_;
0578   MonitorElement* meEle_eta_MTD_4sigma_Bkg_EB_;
0579   MonitorElement* meEle_phi_MTD_4sigma_Bkg_EB_;
0580 
0581   // Adding histograms for endcap electrons
0582   MonitorElement* meEleISO_Ntracks_Bkg_EE_;
0583   MonitorElement* meEleISO_chIso_Bkg_EE_;
0584   MonitorElement* meEleISO_rel_chIso_Bkg_EE_;
0585   MonitorElement* meEleISO_Ntracks_MTD_1_Bkg_EE_;
0586   MonitorElement* meEleISO_chIso_MTD_1_Bkg_EE_;
0587   MonitorElement* meEleISO_rel_chIso_MTD_1_Bkg_EE_;
0588 
0589   MonitorElement* meEleISO_Ntracks_MTD_2_Bkg_EE_;
0590   MonitorElement* meEleISO_chIso_MTD_2_Bkg_EE_;
0591   MonitorElement* meEleISO_rel_chIso_MTD_2_Bkg_EE_;
0592 
0593   MonitorElement* meEleISO_Ntracks_MTD_3_Bkg_EE_;
0594   MonitorElement* meEleISO_chIso_MTD_3_Bkg_EE_;
0595   MonitorElement* meEleISO_rel_chIso_MTD_3_Bkg_EE_;
0596 
0597   MonitorElement* meEleISO_Ntracks_MTD_4_Bkg_EE_;
0598   MonitorElement* meEleISO_chIso_MTD_4_Bkg_EE_;
0599   MonitorElement* meEleISO_rel_chIso_MTD_4_Bkg_EE_;
0600 
0601   MonitorElement* meEleISO_Ntracks_MTD_5_Bkg_EE_;
0602   MonitorElement* meEleISO_chIso_MTD_5_Bkg_EE_;
0603   MonitorElement* meEleISO_rel_chIso_MTD_5_Bkg_EE_;
0604 
0605   MonitorElement* meEleISO_Ntracks_MTD_6_Bkg_EE_;
0606   MonitorElement* meEleISO_chIso_MTD_6_Bkg_EE_;
0607   MonitorElement* meEleISO_rel_chIso_MTD_6_Bkg_EE_;
0608 
0609   MonitorElement* meEleISO_Ntracks_MTD_7_Bkg_EE_;
0610   MonitorElement* meEleISO_chIso_MTD_7_Bkg_EE_;
0611   MonitorElement* meEleISO_rel_chIso_MTD_7_Bkg_EE_;
0612 
0613   MonitorElement* meEleISO_Ntracks_MTD_sim_1_Bkg_EE_;
0614   MonitorElement* meEleISO_chIso_MTD_sim_1_Bkg_EE_;
0615   MonitorElement* meEleISO_rel_chIso_MTD_sim_1_Bkg_EE_;
0616 
0617   MonitorElement* meEleISO_Ntracks_MTD_sim_2_Bkg_EE_;
0618   MonitorElement* meEleISO_chIso_MTD_sim_2_Bkg_EE_;
0619   MonitorElement* meEleISO_rel_chIso_MTD_sim_2_Bkg_EE_;
0620 
0621   MonitorElement* meEleISO_Ntracks_MTD_sim_3_Bkg_EE_;
0622   MonitorElement* meEleISO_chIso_MTD_sim_3_Bkg_EE_;
0623   MonitorElement* meEleISO_rel_chIso_MTD_sim_3_Bkg_EE_;
0624 
0625   MonitorElement* meEleISO_Ntracks_MTD_sim_4_Bkg_EE_;
0626   MonitorElement* meEleISO_chIso_MTD_sim_4_Bkg_EE_;
0627   MonitorElement* meEleISO_rel_chIso_MTD_sim_4_Bkg_EE_;
0628 
0629   MonitorElement* meEleISO_Ntracks_MTD_sim_5_Bkg_EE_;
0630   MonitorElement* meEleISO_chIso_MTD_sim_5_Bkg_EE_;
0631   MonitorElement* meEleISO_rel_chIso_MTD_sim_5_Bkg_EE_;
0632 
0633   MonitorElement* meEleISO_Ntracks_MTD_sim_6_Bkg_EE_;
0634   MonitorElement* meEleISO_chIso_MTD_sim_6_Bkg_EE_;
0635   MonitorElement* meEleISO_rel_chIso_MTD_sim_6_Bkg_EE_;
0636 
0637   MonitorElement* meEleISO_Ntracks_MTD_sim_7_Bkg_EE_;
0638   MonitorElement* meEleISO_chIso_MTD_sim_7_Bkg_EE_;
0639   MonitorElement* meEleISO_rel_chIso_MTD_sim_7_Bkg_EE_;
0640 
0641   MonitorElement* meEleISO_Ntracks_gen_Bkg_EE_;
0642   MonitorElement* meEleISO_chIso_gen_Bkg_EE_;
0643   MonitorElement* meEleISO_rel_chIso_gen_Bkg_EE_;
0644 
0645   MonitorElement* meEleISO_Ntracks_MTD_2sigma_Bkg_EE_;
0646   MonitorElement* meEleISO_chIso_MTD_2sigma_Bkg_EE_;
0647   MonitorElement* meEleISO_rel_chIso_MTD_2sigma_Bkg_EE_;
0648 
0649   MonitorElement* meEleISO_Ntracks_MTD_3sigma_Bkg_EE_;
0650   MonitorElement* meEleISO_chIso_MTD_3sigma_Bkg_EE_;
0651   MonitorElement* meEleISO_rel_chIso_MTD_3sigma_Bkg_EE_;
0652 
0653   MonitorElement* meEleISO_Ntracks_MTD_4sigma_Bkg_EE_;
0654   MonitorElement* meEleISO_chIso_MTD_4sigma_Bkg_EE_;
0655   MonitorElement* meEleISO_rel_chIso_MTD_4sigma_Bkg_EE_;
0656 
0657   MonitorElement* meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EE_;
0658   MonitorElement* meEleISO_chIso_MTD_sim_2sigma_Bkg_EE_;
0659   MonitorElement* meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EE_;
0660 
0661   MonitorElement* meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EE_;
0662   MonitorElement* meEleISO_chIso_MTD_sim_3sigma_Bkg_EE_;
0663   MonitorElement* meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EE_;
0664 
0665   MonitorElement* meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EE_;
0666   MonitorElement* meEleISO_chIso_MTD_sim_4sigma_Bkg_EE_;
0667   MonitorElement* meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EE_;
0668 
0669   MonitorElement* meEle_pt_tot_Bkg_EE_;
0670   MonitorElement* meEle_pt_sim_tot_Bkg_EE_;
0671   MonitorElement* meEle_eta_tot_Bkg_EE_;
0672   MonitorElement* meEle_phi_tot_Bkg_EE_;
0673   MonitorElement* meEle_pt_MTD_1_Bkg_EE_;
0674   MonitorElement* meEle_pt_sim_MTD_1_Bkg_EE_;
0675   MonitorElement* meEle_eta_MTD_1_Bkg_EE_;
0676   MonitorElement* meEle_phi_MTD_1_Bkg_EE_;
0677 
0678   MonitorElement* meEle_pt_MTD_2_Bkg_EE_;
0679   MonitorElement* meEle_pt_sim_MTD_2_Bkg_EE_;
0680   MonitorElement* meEle_eta_MTD_2_Bkg_EE_;
0681   MonitorElement* meEle_phi_MTD_2_Bkg_EE_;
0682 
0683   MonitorElement* meEle_pt_MTD_3_Bkg_EE_;
0684   MonitorElement* meEle_pt_sim_MTD_3_Bkg_EE_;
0685   MonitorElement* meEle_eta_MTD_3_Bkg_EE_;
0686   MonitorElement* meEle_phi_MTD_3_Bkg_EE_;
0687 
0688   MonitorElement* meEle_pt_MTD_4_Bkg_EE_;
0689   MonitorElement* meEle_pt_sim_MTD_4_Bkg_EE_;
0690   MonitorElement* meEle_eta_MTD_4_Bkg_EE_;
0691   MonitorElement* meEle_phi_MTD_4_Bkg_EE_;
0692 
0693   MonitorElement* meEle_pt_MTD_5_Bkg_EE_;
0694   MonitorElement* meEle_pt_sim_MTD_5_Bkg_EE_;
0695   MonitorElement* meEle_eta_MTD_5_Bkg_EE_;
0696   MonitorElement* meEle_phi_MTD_5_Bkg_EE_;
0697 
0698   MonitorElement* meEle_pt_MTD_6_Bkg_EE_;
0699   MonitorElement* meEle_pt_sim_MTD_6_Bkg_EE_;
0700   MonitorElement* meEle_eta_MTD_6_Bkg_EE_;
0701   MonitorElement* meEle_phi_MTD_6_Bkg_EE_;
0702 
0703   MonitorElement* meEle_pt_MTD_7_Bkg_EE_;
0704   MonitorElement* meEle_pt_sim_MTD_7_Bkg_EE_;
0705   MonitorElement* meEle_eta_MTD_7_Bkg_EE_;
0706   MonitorElement* meEle_phi_MTD_7_Bkg_EE_;
0707 
0708   MonitorElement* meEle_pt_noMTD_Bkg_EE_;
0709   MonitorElement* meEle_eta_noMTD_Bkg_EE_;
0710   MonitorElement* meEle_phi_noMTD_Bkg_EE_;
0711 
0712   MonitorElement* meEle_pt_gen_Bkg_EE_;
0713   MonitorElement* meEle_eta_gen_Bkg_EE_;
0714   MonitorElement* meEle_phi_gen_Bkg_EE_;
0715 
0716   MonitorElement* meEle_pt_MTD_2sigma_Bkg_EE_;
0717   MonitorElement* meEle_pt_sim_MTD_2sigma_Bkg_EE_;
0718   MonitorElement* meEle_eta_MTD_2sigma_Bkg_EE_;
0719   MonitorElement* meEle_phi_MTD_2sigma_Bkg_EE_;
0720 
0721   MonitorElement* meEle_pt_MTD_3sigma_Bkg_EE_;
0722   MonitorElement* meEle_pt_sim_MTD_3sigma_Bkg_EE_;
0723   MonitorElement* meEle_eta_MTD_3sigma_Bkg_EE_;
0724   MonitorElement* meEle_phi_MTD_3sigma_Bkg_EE_;
0725 
0726   MonitorElement* meEle_pt_MTD_4sigma_Bkg_EE_;
0727   MonitorElement* meEle_pt_sim_MTD_4sigma_Bkg_EE_;
0728   MonitorElement* meEle_eta_MTD_4sigma_Bkg_EE_;
0729   MonitorElement* meEle_phi_MTD_4sigma_Bkg_EE_;
0730   // Background histograms end
0731 
0732   // promt part for histogram vectors
0733   std::vector<MonitorElement*> Ntracks_EB_list_Sig;
0734   std::vector<MonitorElement*> ch_iso_EB_list_Sig;
0735   std::vector<MonitorElement*> rel_ch_iso_EB_list_Sig;
0736 
0737   std::vector<MonitorElement*> Ntracks_EE_list_Sig;
0738   std::vector<MonitorElement*> ch_iso_EE_list_Sig;
0739   std::vector<MonitorElement*> rel_ch_iso_EE_list_Sig;
0740 
0741   std::vector<MonitorElement*> Ntracks_sim_EB_list_Sig;
0742   std::vector<MonitorElement*> ch_iso_sim_EB_list_Sig;
0743   std::vector<MonitorElement*> rel_ch_iso_sim_EB_list_Sig;
0744 
0745   std::vector<MonitorElement*> Ntracks_sim_EE_list_Sig;
0746   std::vector<MonitorElement*> ch_iso_sim_EE_list_Sig;
0747   std::vector<MonitorElement*> rel_ch_iso_sim_EE_list_Sig;
0748 
0749   std::vector<MonitorElement*> Ele_pT_MTD_EB_list_Sig;
0750   std::vector<MonitorElement*> Ele_pT_sim_MTD_EB_list_Sig;
0751   std::vector<MonitorElement*> Ele_eta_MTD_EB_list_Sig;
0752   std::vector<MonitorElement*> Ele_phi_MTD_EB_list_Sig;
0753 
0754   std::vector<MonitorElement*> Ele_pT_MTD_EE_list_Sig;
0755   std::vector<MonitorElement*> Ele_pT_sim_MTD_EE_list_Sig;
0756   std::vector<MonitorElement*> Ele_eta_MTD_EE_list_Sig;
0757   std::vector<MonitorElement*> Ele_phi_MTD_EE_list_Sig;
0758 
0759   std::vector<MonitorElement*> Ntracks_EB_list_Significance_Sig;
0760   std::vector<MonitorElement*> ch_iso_EB_list_Significance_Sig;
0761   std::vector<MonitorElement*> rel_ch_iso_EB_list_Significance_Sig;
0762 
0763   std::vector<MonitorElement*> Ntracks_EE_list_Significance_Sig;
0764   std::vector<MonitorElement*> ch_iso_EE_list_Significance_Sig;
0765   std::vector<MonitorElement*> rel_ch_iso_EE_list_Significance_Sig;
0766 
0767   std::vector<MonitorElement*> Ntracks_sim_EB_list_Significance_Sig;
0768   std::vector<MonitorElement*> ch_iso_sim_EB_list_Significance_Sig;
0769   std::vector<MonitorElement*> rel_ch_iso_sim_EB_list_Significance_Sig;
0770 
0771   std::vector<MonitorElement*> Ntracks_sim_EE_list_Significance_Sig;
0772   std::vector<MonitorElement*> ch_iso_sim_EE_list_Significance_Sig;
0773   std::vector<MonitorElement*> rel_ch_iso_sim_EE_list_Significance_Sig;
0774 
0775   std::vector<MonitorElement*> Ele_pT_MTD_EB_list_Significance_Sig;
0776   std::vector<MonitorElement*> Ele_pT_sim_MTD_EB_list_Significance_Sig;
0777   std::vector<MonitorElement*> Ele_eta_MTD_EB_list_Significance_Sig;
0778   std::vector<MonitorElement*> Ele_phi_MTD_EB_list_Significance_Sig;
0779 
0780   std::vector<MonitorElement*> Ele_pT_MTD_EE_list_Significance_Sig;
0781   std::vector<MonitorElement*> Ele_pT_sim_MTD_EE_list_Significance_Sig;
0782   std::vector<MonitorElement*> Ele_eta_MTD_EE_list_Significance_Sig;
0783   std::vector<MonitorElement*> Ele_phi_MTD_EE_list_Significance_Sig;
0784 
0785   // Non-promt part for histogram vectors
0786   std::vector<MonitorElement*> Ntracks_EB_list_Bkg;
0787   std::vector<MonitorElement*> ch_iso_EB_list_Bkg;
0788   std::vector<MonitorElement*> rel_ch_iso_EB_list_Bkg;
0789 
0790   std::vector<MonitorElement*> Ntracks_EE_list_Bkg;
0791   std::vector<MonitorElement*> ch_iso_EE_list_Bkg;
0792   std::vector<MonitorElement*> rel_ch_iso_EE_list_Bkg;
0793 
0794   std::vector<MonitorElement*> Ntracks_sim_EB_list_Bkg;
0795   std::vector<MonitorElement*> ch_iso_sim_EB_list_Bkg;
0796   std::vector<MonitorElement*> rel_ch_iso_sim_EB_list_Bkg;
0797 
0798   std::vector<MonitorElement*> Ntracks_sim_EE_list_Bkg;
0799   std::vector<MonitorElement*> ch_iso_sim_EE_list_Bkg;
0800   std::vector<MonitorElement*> rel_ch_iso_sim_EE_list_Bkg;
0801 
0802   std::vector<MonitorElement*> Ele_pT_MTD_EB_list_Bkg;
0803   std::vector<MonitorElement*> Ele_pT_sim_MTD_EB_list_Bkg;
0804   std::vector<MonitorElement*> Ele_eta_MTD_EB_list_Bkg;
0805   std::vector<MonitorElement*> Ele_phi_MTD_EB_list_Bkg;
0806 
0807   std::vector<MonitorElement*> Ele_pT_MTD_EE_list_Bkg;
0808   std::vector<MonitorElement*> Ele_pT_sim_MTD_EE_list_Bkg;
0809   std::vector<MonitorElement*> Ele_eta_MTD_EE_list_Bkg;
0810   std::vector<MonitorElement*> Ele_phi_MTD_EE_list_Bkg;
0811 
0812   std::vector<MonitorElement*> Ntracks_EB_list_Significance_Bkg;
0813   std::vector<MonitorElement*> ch_iso_EB_list_Significance_Bkg;
0814   std::vector<MonitorElement*> rel_ch_iso_EB_list_Significance_Bkg;
0815 
0816   std::vector<MonitorElement*> Ntracks_EE_list_Significance_Bkg;
0817   std::vector<MonitorElement*> ch_iso_EE_list_Significance_Bkg;
0818   std::vector<MonitorElement*> rel_ch_iso_EE_list_Significance_Bkg;
0819 
0820   std::vector<MonitorElement*> Ntracks_sim_EB_list_Significance_Bkg;
0821   std::vector<MonitorElement*> ch_iso_sim_EB_list_Significance_Bkg;
0822   std::vector<MonitorElement*> rel_ch_iso_sim_EB_list_Significance_Bkg;
0823 
0824   std::vector<MonitorElement*> Ntracks_sim_EE_list_Significance_Bkg;
0825   std::vector<MonitorElement*> ch_iso_sim_EE_list_Significance_Bkg;
0826   std::vector<MonitorElement*> rel_ch_iso_sim_EE_list_Significance_Bkg;
0827 
0828   std::vector<MonitorElement*> Ele_pT_MTD_EB_list_Significance_Bkg;
0829   std::vector<MonitorElement*> Ele_pT_sim_MTD_EB_list_Significance_Bkg;
0830   std::vector<MonitorElement*> Ele_eta_MTD_EB_list_Significance_Bkg;
0831   std::vector<MonitorElement*> Ele_phi_MTD_EB_list_Significance_Bkg;
0832 
0833   std::vector<MonitorElement*> Ele_pT_MTD_EE_list_Significance_Bkg;
0834   std::vector<MonitorElement*> Ele_pT_sim_MTD_EE_list_Significance_Bkg;
0835   std::vector<MonitorElement*> Ele_eta_MTD_EE_list_Significance_Bkg;
0836   std::vector<MonitorElement*> Ele_phi_MTD_EE_list_Significance_Bkg;
0837 
0838   // dt distribution part for histogram vectors
0839   std::vector<MonitorElement*> general_pT_list;
0840   std::vector<MonitorElement*> general_eta_list;
0841 
0842   std::vector<MonitorElement*> general_pT_Signif_list;
0843   std::vector<MonitorElement*> general_eta_Signif_list;
0844 };
0845 
0846 // ------------ constructor and destructor --------------
0847 
0848 MtdEleIsoValidation::MtdEleIsoValidation(const edm::ParameterSet& iConfig)
0849     : folder_(iConfig.getParameter<std::string>("folder")),
0850       trackMinPt_(iConfig.getParameter<double>("trackMinimumPt")),
0851       trackMinEta_(iConfig.getParameter<double>("trackMinimumEta")),
0852       trackMaxEta_(iConfig.getParameter<double>("trackMaximumEta")),
0853       rel_iso_cut_(iConfig.getParameter<double>("rel_iso_cut")),
0854       track_match_PV_(iConfig.getParameter<bool>("optionTrackMatchToPV")),
0855       dt_sig_track_(iConfig.getParameter<bool>("option_dtToTrack")),
0856       optionalPlots_(iConfig.getParameter<bool>("option_plots")),
0857       min_dR_cut(iConfig.getParameter<double>("min_dR_cut")),
0858       max_dR_cut(iConfig.getParameter<double>("max_dR_cut")),
0859       min_pt_cut_EB(iConfig.getParameter<double>("min_pt_cut_EB")),
0860       min_pt_cut_EE(iConfig.getParameter<double>("min_pt_cut_EE")),
0861       max_dz_cut_EB(iConfig.getParameter<double>("max_dz_cut_EB")),
0862       max_dz_cut_EE(iConfig.getParameter<double>("max_dz_cut_EE")),
0863       max_dz_vtx_cut(iConfig.getParameter<double>("max_dz_vtx_cut")),
0864       max_dxy_vtx_cut(iConfig.getParameter<double>("max_dxy_vtx_cut")),
0865       min_strip_cut(iConfig.getParameter<double>("min_strip_cut")),
0866       min_track_mtd_mva_cut(iConfig.getParameter<double>("min_track_mtd_mva_cut")) {
0867   GenRecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTagG"));
0868   RecVertexToken_ =
0869       consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("inputTag_vtx"));  // Vtx 4D collection
0870 
0871   GsfElectronToken_EB_ = consumes<reco::GsfElectronCollection>(
0872       iConfig.getParameter<edm::InputTag>("inputEle_EB"));  // Barrel electron collection input/token
0873   GsfElectronToken_EE_ = consumes<reco::GsfElectronCollection>(
0874       iConfig.getParameter<edm::InputTag>("inputEle_EE"));  // Endcap electron collection input/token
0875 
0876   t0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0PID"));
0877   Sigmat0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0PID"));
0878   trackMVAQualToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("trackMVAQual"));
0879 
0880   recoToSimAssociationToken_ =
0881       consumes<reco::RecoToSimCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
0882 }
0883 
0884 MtdEleIsoValidation::~MtdEleIsoValidation() {}
0885 
0886 // ------------ method called for each event  ------------
0887 void MtdEleIsoValidation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
0888   using namespace edm;
0889   using namespace geant_units::operators;
0890   using namespace std;
0891 
0892   auto GenRecTrackHandle = iEvent.getHandle(GenRecTrackToken_);
0893 
0894   auto VertexHandle = iEvent.getHandle(RecVertexToken_);
0895   std::vector<reco::Vertex> vertices = *VertexHandle;
0896 
0897   const auto& t0Pid = iEvent.get(t0PidToken_);
0898   const auto& Sigmat0Pid = iEvent.get(Sigmat0PidToken_);
0899   const auto& mtdQualMVA = iEvent.get(trackMVAQualToken_);
0900 
0901   auto eleHandle_EB = makeValid(iEvent.getHandle(GsfElectronToken_EB_));
0902   reco::GsfElectronCollection eleColl_EB = *(eleHandle_EB.product());
0903 
0904   auto eleHandle_EE = makeValid(iEvent.getHandle(GsfElectronToken_EE_));
0905   reco::GsfElectronCollection eleColl_EE = *(eleHandle_EE.product());
0906 
0907   auto recoToSimH = makeValid(iEvent.getHandle(recoToSimAssociationToken_));
0908   const reco::RecoToSimCollection* r2s_ = recoToSimH.product();
0909 
0910   // Creating combined electron collection
0911   std::vector<reco::GsfElectron> localEleCollection;
0912   localEleCollection.reserve(eleColl_EB.size() + eleColl_EE.size());
0913   for (const auto& ele_EB : eleColl_EB) {
0914     if (ele_EB.isEB()) {
0915       localEleCollection.emplace_back(ele_EB);
0916     }
0917   }
0918   for (const auto& ele_EE : eleColl_EE) {
0919     if (ele_EE.isEE()) {
0920       localEleCollection.emplace_back(ele_EE);
0921     }
0922   }
0923   localEleCollection.shrink_to_fit();
0924 
0925   reco::Vertex Vtx_chosen;
0926   // This part has to be included, because in ~1% of the events, the "good" vertex is the 1st one not the 0th one in the collection
0927   for (int iVtx = 0; iVtx < (int)vertices.size(); iVtx++) {
0928     const reco::Vertex& vertex = vertices.at(iVtx);
0929     if (!vertex.isFake() && vertex.ndof() >= 4) {
0930       Vtx_chosen = vertex;
0931       break;
0932     }
0933   }
0934 
0935   auto pdgCheck = [](int pdg) {
0936     pdg = std::abs(pdg);
0937     return (pdg == 23 or pdg == 24 or pdg == 15 or pdg == 11);  // some electrons are mothers to themselves?
0938   };
0939 
0940   for (const auto& ele : localEleCollection) {
0941     bool ele_Promt = false;
0942 
0943     float ele_track_source_dz = std::abs(ele.gsfTrack()->dz(Vtx_chosen.position()));
0944     float ele_track_source_dxy = std::abs(ele.gsfTrack()->dxy(Vtx_chosen.position()));
0945 
0946     const reco::TrackRef ele_TrkRef = ele.core()->ctfTrack();
0947     double tsim_ele = -1.;
0948     double ele_sim_pt = -1.;
0949     double ele_sim_phi = -1.;
0950     double ele_sim_eta = -1.;
0951 
0952     // selecting "good" RECO electrons
0953     // PARAM
0954     if (ele.pt() < 10 || std::abs(ele.eta()) > 2.4 || ele_track_source_dz > max_dz_vtx_cut ||
0955         ele_track_source_dxy > max_dxy_vtx_cut)
0956       continue;
0957 
0958     // association with tracking particle to have sim info
0959     const reco::TrackBaseRef trkrefb(ele_TrkRef);
0960     auto found = r2s_->find(trkrefb);
0961     if (found != r2s_->end()) {
0962       const auto& tp = (found->val)[0];
0963       tsim_ele = (tp.first)->parentVertex()->position().t() * 1e9;
0964       ele_sim_pt = (tp.first)->pt();
0965       ele_sim_phi = (tp.first)->phi();
0966       ele_sim_eta = (tp.first)->eta();
0967       // check that the genParticle vector is not empty
0968       if (tp.first->status() != -99) {
0969         const auto genParticle = *(tp.first->genParticles()[0]);
0970         // check if prompt (not from hadron, muon, or tau decay) and final state
0971         // or if is a direct decay product of a prompt tau and is final state
0972         if ((genParticle.isPromptFinalState() or genParticle.isDirectPromptTauDecayProductFinalState()) and
0973             pdgCheck(genParticle.mother()->pdgId())) {
0974           ele_Promt = true;
0975           // TODO get simtrackster from mtd, simtrack to tp and check that a recocluster was there
0976         }
0977       }
0978     }
0979 
0980     math::XYZVector EleSigTrackMomentumAtVtx = ele.gsfTrack()->momentum();
0981     double EleSigTrackEtaAtVtx = ele.gsfTrack()->eta();
0982 
0983     double ele_sigTrkTime = -1;
0984     double ele_sigTrkTimeErr = -1;
0985     double ele_sigTrkMtdMva = -1;
0986 
0987     // if we found a track match, we add MTD timing information for it
0988     if (ele_TrkRef.isNonnull()) {
0989       // track pT/dz cuts
0990       bool Barrel_ele = ele.isEB();
0991       float min_pt_cut = Barrel_ele ? min_pt_cut_EB : min_pt_cut_EE;
0992       float max_dz_cut = Barrel_ele ? max_dz_cut_EB : max_dz_cut_EE;
0993 
0994       ele_sigTrkTime = t0Pid[ele_TrkRef];
0995       ele_sigTrkMtdMva = mtdQualMVA[ele_TrkRef];
0996       ele_sigTrkTimeErr = (ele_sigTrkMtdMva > min_track_mtd_mva_cut) ? Sigmat0Pid[ele_TrkRef] : -1;
0997 
0998       meEle_avg_error_SigTrk_check_->Fill(ele_sigTrkTimeErr);
0999 
1000       if (ele_Promt) {
1001         // For signal (promt)
1002         if (Barrel_ele) {
1003           // All selected electron information for efficiency plots later
1004           meEle_pt_tot_Sig_EB_->Fill(ele.pt());
1005           meEle_pt_sim_tot_Sig_EB_->Fill(ele_sim_pt);
1006           meEle_eta_tot_Sig_EB_->Fill(std::abs(ele.eta()));
1007           meEle_phi_tot_Sig_EB_->Fill(ele.phi());
1008         } else {
1009           // All selected electron information for efficiency plots later
1010           meEle_pt_tot_Sig_EE_->Fill(ele.pt());
1011           meEle_pt_sim_tot_Sig_EE_->Fill(ele_sim_pt);
1012           meEle_eta_tot_Sig_EE_->Fill(std::abs(ele.eta()));
1013           meEle_phi_tot_Sig_EE_->Fill(ele.phi());
1014         }
1015       } else {
1016         // For background (non-promt)
1017         if (Barrel_ele) {
1018           meEle_pt_tot_Bkg_EB_->Fill(ele.pt());
1019           meEle_pt_sim_tot_Bkg_EB_->Fill(ele_sim_pt);
1020           meEle_eta_tot_Bkg_EB_->Fill(std::abs(ele.eta()));
1021           meEle_phi_tot_Bkg_EB_->Fill(ele.phi());
1022         } else {
1023           meEle_pt_tot_Bkg_EE_->Fill(ele.pt());
1024           meEle_pt_sim_tot_Bkg_EE_->Fill(ele_sim_pt);
1025           meEle_eta_tot_Bkg_EE_->Fill(std::abs(ele.eta()));
1026           meEle_phi_tot_Bkg_EE_->Fill(ele.phi());
1027         }
1028       }
1029 
1030       int N_tracks_noMTD = 0;
1031       double pT_sum_noMTD = 0;
1032       double rel_pT_sum_noMTD = 0;
1033       std::vector<int> N_tracks_MTD{0, 0, 0, 0, 0, 0, 0};
1034       std::vector<double> pT_sum_MTD{0, 0, 0, 0, 0, 0, 0};
1035       std::vector<double> rel_pT_sum_MTD{0, 0, 0, 0, 0, 0, 0};
1036 
1037       std::vector<int> N_tracks_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1038       std::vector<double> pT_sum_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1039       std::vector<double> rel_pT_sum_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1040       int N_tracks_gen = 0;
1041       double pT_sum_gen = 0;
1042       double rel_pT_sum_gen = 0;
1043 
1044       std::vector<int> N_tracks_MTD_significance{0, 0, 0};
1045       std::vector<double> pT_sum_MTD_significance{0, 0, 0};
1046       std::vector<double> rel_pT_sum_MTD_significance{0, 0, 0};
1047 
1048       std::vector<int> N_tracks_sim_MTD_significance{0, 0, 0};
1049       std::vector<double> pT_sum_sim_MTD_significance{0, 0, 0};
1050       std::vector<double> rel_pT_sum_sim_MTD_significance{0, 0, 0};
1051 
1052       int general_index = 0;
1053       for (const auto& trackGen : *GenRecTrackHandle) {
1054         const reco::TrackRef trackref_general(GenRecTrackHandle, general_index);
1055         general_index++;
1056 
1057         // Skip electron track
1058         if (trackref_general == ele_TrkRef)
1059           continue;
1060 
1061         if (trackGen.pt() < min_pt_cut) {
1062           continue;
1063         }
1064         if (std::abs(trackGen.vz() - ele.gsfTrack()->vz()) > max_dz_cut) {
1065           continue;
1066         }
1067 
1068         // cut for general track matching to PV
1069         if (track_match_PV_) {
1070           if (Vtx_chosen.trackWeight(trackref_general) < 0.5) {
1071             continue;
1072           }
1073         }
1074 
1075         double dR = reco::deltaR(trackGen.momentum(), EleSigTrackMomentumAtVtx);
1076         double deta = std::abs(trackGen.eta() - EleSigTrackEtaAtVtx);
1077 
1078         // restrict to tracks in the isolation cone
1079         if (dR < min_dR_cut || dR > max_dR_cut || deta < min_strip_cut)
1080           continue;
1081 
1082         // no MTD case
1083         ++N_tracks_noMTD;
1084         pT_sum_noMTD += trackGen.pt();
1085 
1086         // MTD case
1087         const reco::TrackBaseRef trkrefBase(trackref_general);
1088         auto TPmatched = r2s_->find(trkrefBase);
1089         double tsim_trk = -1.;
1090         double trk_ptSim = -1.;
1091         bool genMatched = false;
1092         if (TPmatched != r2s_->end()) {
1093           // reco track matched to a TP
1094           const auto& tp = (TPmatched->val)[0];
1095           tsim_trk = (tp.first)->parentVertex()->position().t() * 1e9;
1096           trk_ptSim = (tp.first)->pt();
1097           // check that the genParticle vector is not empty
1098           if (tp.first->status() != -99) {
1099             genMatched = true;
1100             meTrk_genMatch_check_->Fill(1);
1101           } else {
1102             meTrk_genMatch_check_->Fill(0);
1103           }
1104         }
1105 
1106         double TrkMTDTime = t0Pid[trackref_general];
1107         double TrkMTDMva = mtdQualMVA[trackref_general];
1108         double TrkMTDTimeErr = (TrkMTDMva > min_track_mtd_mva_cut) ? Sigmat0Pid[trackref_general] : -1;
1109 
1110         meEle_avg_error_PUTrk_check_->Fill(TrkMTDTimeErr);
1111 
1112         // MTD GEN case
1113         if (genMatched) {
1114           N_tracks_gen++;
1115           pT_sum_gen += trk_ptSim;
1116         }
1117 
1118         // dt with the track
1119         if (dt_sig_track_) {
1120           double dt_sigTrk = 0;
1121           double dt_sigTrk_signif = 0;
1122           double dt_sim_sigTrk = 0;
1123           double dt_sim_sigTrk_signif = 0;
1124 
1125           // MTD SIM CASE
1126           if (std::abs(tsim_trk) > 0 && std::abs(tsim_ele) > 0 && trk_ptSim > 0) {
1127             dt_sim_sigTrk = std::abs(tsim_trk - tsim_ele);
1128             dt_sim_sigTrk_signif = dt_sim_sigTrk / std::sqrt(avg_sim_PUtrack_t_err * avg_sim_PUtrack_t_err +
1129                                                              avg_sim_sigTrk_t_err * avg_sim_sigTrk_t_err);
1130 
1131             if (optionalPlots_) {
1132               // absolute timing cuts
1133               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1134                 if (dt_sim_sigTrk < max_dt_track_cut[i]) {
1135                   N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1136                   pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1137                 }
1138               }
1139             }
1140             // significance cuts
1141             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1142               if (dt_sim_sigTrk_signif < max_dt_significance_cut[i]) {
1143                 N_tracks_sim_MTD_significance[i]++;
1144                 pT_sum_sim_MTD_significance[i] += trk_ptSim;
1145               }
1146             }
1147 
1148           } else {
1149             // if there is no error for MTD information, we count the MTD isolation case same as noMTD
1150             if (optionalPlots_) {
1151               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1152                 N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1153                 pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1154               }
1155             }
1156             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1157               N_tracks_sim_MTD_significance[i]++;
1158               pT_sum_sim_MTD_significance[i] += trk_ptSim;
1159             }
1160           }
1161 
1162           // MTD reco case
1163           if (TrkMTDTimeErr > 0 && ele_sigTrkTimeErr > 0) {
1164             dt_sigTrk = std::abs(TrkMTDTime - ele_sigTrkTime);
1165             dt_sigTrk_signif =
1166                 dt_sigTrk / std::sqrt(TrkMTDTimeErr * TrkMTDTimeErr + ele_sigTrkTimeErr * ele_sigTrkTimeErr);
1167 
1168             meEle_no_dt_check_->Fill(1);
1169             if (optionalPlots_) {
1170               // absolute timing cuts
1171               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1172                 if (dt_sigTrk < max_dt_track_cut[i]) {
1173                   N_tracks_MTD[i] = N_tracks_MTD[i] + 1;
1174                   pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();
1175                 }
1176               }
1177             }
1178             // significance cuts
1179             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1180               if (dt_sigTrk_signif < max_dt_significance_cut[i]) {
1181                 N_tracks_MTD_significance[i]++;
1182                 pT_sum_MTD_significance[i] += trackGen.pt();
1183               }
1184             }
1185 
1186           } else {
1187             // if there is no error for MTD information, we count the MTD isolation case same as noMTD
1188             if (optionalPlots_) {
1189               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1190                 N_tracks_MTD[i] = N_tracks_MTD[i] + 1;          // N_tracks_noMTD
1191                 pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();  // pT sum
1192               }
1193             }
1194             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1195               N_tracks_MTD_significance[i]++;
1196               pT_sum_MTD_significance[i] += trackGen.pt();
1197             }
1198             meEle_no_dt_check_->Fill(0);
1199           }
1200 
1201           if (optionalPlots_) {
1202             for (long unsigned int i = 0; i < (pT_bins_dt_distrb.size() - 1); i++) {
1203               //stuff general pT
1204               if (ele.pt() > pT_bins_dt_distrb[i] && ele.pt() < pT_bins_dt_distrb[i + 1]) {
1205                 general_pT_list[i]->Fill(dt_sigTrk);
1206                 general_pT_Signif_list[i]->Fill(dt_sigTrk_signif);
1207               }
1208             }
1209 
1210             for (long unsigned int i = 0; i < (eta_bins_dt_distrib.size() - 1); i++) {
1211               //stuff general eta
1212               if (std::abs(ele.eta()) > eta_bins_dt_distrib[i] && std::abs(ele.eta()) < eta_bins_dt_distrib[i + 1]) {
1213                 general_eta_list[i]->Fill(dt_sigTrk);
1214                 general_eta_Signif_list[i]->Fill(dt_sigTrk_signif);
1215               }
1216             }
1217           }  // End of optional dt distributions plots
1218 
1219           // dt with the vertex
1220         } else {
1221           double dt_vtx = 0;  // dt regular track vs vtx
1222           double dt_vtx_signif = 0;
1223 
1224           double dt_sim_vtx = 0;  // dt regular track vs vtx
1225           double dt_sim_vtx_signif = 0;
1226 
1227           // MTD SIM case
1228           if (std::abs(tsim_trk) > 0 && Vtx_chosen.tError() > 0 && trk_ptSim > 0) {
1229             dt_sim_vtx = std::abs(tsim_trk - Vtx_chosen.t());
1230             dt_sim_vtx_signif = dt_sim_vtx / std::sqrt(avg_sim_PUtrack_t_err * avg_sim_PUtrack_t_err +
1231                                                        Vtx_chosen.tError() * Vtx_chosen.tError());
1232             if (optionalPlots_) {
1233               // absolute timing cuts
1234               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1235                 if (dt_sim_vtx < max_dt_vtx_cut[i]) {
1236                   N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1237                   pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1238                 }
1239               }
1240             }
1241             // significance timing cuts
1242             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1243               if (dt_sim_vtx_signif < max_dt_significance_cut[i]) {
1244                 N_tracks_sim_MTD_significance[i]++;
1245                 pT_sum_sim_MTD_significance[i] += trk_ptSim;
1246               }
1247             }
1248           } else {
1249             if (optionalPlots_) {
1250               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1251                 N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;      // N_tracks_noMTD
1252                 pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;  // pT_sum_noMTD
1253               }
1254             }
1255             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1256               N_tracks_sim_MTD_significance[i]++;
1257               pT_sum_sim_MTD_significance[i] += trk_ptSim;
1258             }
1259           }
1260 
1261           // MTD RECO case
1262           if (TrkMTDTimeErr > 0 && Vtx_chosen.tError() > 0) {
1263             dt_vtx = std::abs(TrkMTDTime - Vtx_chosen.t());
1264             dt_vtx_signif =
1265                 dt_vtx / std::sqrt(TrkMTDTimeErr * TrkMTDTimeErr + Vtx_chosen.tError() * Vtx_chosen.tError());
1266 
1267             meEle_no_dt_check_->Fill(1);
1268             meEle_avg_error_vtx_check_->Fill(Vtx_chosen.tError());
1269             if (optionalPlots_) {
1270               // absolute timing cuts
1271               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1272                 if (dt_vtx < max_dt_vtx_cut[i]) {
1273                   N_tracks_MTD[i] = N_tracks_MTD[i] + 1;
1274                   pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();
1275                 }
1276               }
1277             }
1278             // significance timing cuts
1279             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1280               if (dt_vtx_signif < max_dt_significance_cut[i]) {
1281                 N_tracks_MTD_significance[i]++;
1282                 pT_sum_MTD_significance[i] += trackGen.pt();
1283               }
1284             }
1285           } else {
1286             if (optionalPlots_) {
1287               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1288                 N_tracks_MTD[i] = N_tracks_MTD[i] + 1;          // N_tracks_noMTD
1289                 pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();  // pT_sum_noMTD
1290               }
1291             }
1292             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1293               N_tracks_MTD_significance[i]++;
1294               pT_sum_MTD_significance[i] += trackGen.pt();
1295             }
1296             meEle_no_dt_check_->Fill(0);
1297           }
1298 
1299           // Optional dt distribution plots
1300           if (optionalPlots_) {
1301             for (long unsigned int i = 0; i < (pT_bins_dt_distrb.size() - 1); i++) {
1302               //stuff general pT
1303               if (ele.pt() > pT_bins_dt_distrb[i] && ele.pt() < pT_bins_dt_distrb[i + 1]) {
1304                 general_pT_list[i]->Fill(dt_vtx);
1305                 general_pT_Signif_list[i]->Fill(dt_vtx_signif);
1306               }
1307             }
1308 
1309             for (long unsigned int i = 0; i < (eta_bins_dt_distrib.size() - 1); i++) {
1310               //stuff general eta
1311               if (std::abs(ele.eta()) > eta_bins_dt_distrib[i] && std::abs(ele.eta()) < eta_bins_dt_distrib[i + 1]) {
1312                 general_eta_list[i]->Fill(dt_vtx);
1313                 general_eta_Signif_list[i]->Fill(dt_vtx_signif);
1314               }
1315             }
1316           }  // End of optional dt distributions plots
1317         }
1318       }
1319       rel_pT_sum_noMTD = pT_sum_noMTD / ele.gsfTrack()->pt();  // rel_ch_iso calculation
1320       if (optionalPlots_) {
1321         for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1322           rel_pT_sum_MTD[i] = pT_sum_MTD[i] / ele.gsfTrack()->pt();
1323           rel_pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] / ele_sim_pt;
1324         }
1325         // now compute the isolation
1326         rel_pT_sum_noMTD = pT_sum_noMTD / ele.gsfTrack()->pt();
1327 
1328         rel_pT_sum_gen = pT_sum_gen / ele.gsfTrack()->pt();
1329       }
1330 
1331       for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1332         rel_pT_sum_MTD_significance[i] = pT_sum_MTD_significance[i] / ele.gsfTrack()->pt();
1333         rel_pT_sum_sim_MTD_significance[i] = pT_sum_sim_MTD_significance[i] / ele_sim_pt;
1334       }
1335 
1336       if (ele_Promt) {  // promt part
1337         if (Barrel_ele) {
1338           meEleISO_Ntracks_Sig_EB_->Fill(N_tracks_noMTD);
1339           meEleISO_chIso_Sig_EB_->Fill(pT_sum_noMTD);
1340           meEleISO_rel_chIso_Sig_EB_->Fill(rel_pT_sum_noMTD);
1341           if (optionalPlots_) {
1342             for (long unsigned int j = 0; j < Ntracks_EB_list_Sig.size(); j++) {
1343               Ntracks_EB_list_Sig[j]->Fill(N_tracks_MTD[j]);
1344               ch_iso_EB_list_Sig[j]->Fill(pT_sum_MTD[j]);
1345               rel_ch_iso_EB_list_Sig[j]->Fill(rel_pT_sum_MTD[j]);
1346 
1347               Ntracks_sim_EB_list_Sig[j]->Fill(N_tracks_sim_MTD[j]);
1348               ch_iso_sim_EB_list_Sig[j]->Fill(pT_sum_sim_MTD[j]);
1349               rel_ch_iso_sim_EB_list_Sig[j]->Fill(rel_pT_sum_sim_MTD[j]);
1350             }
1351             meEleISO_Ntracks_gen_Sig_EB_->Fill(N_tracks_gen);
1352             meEleISO_chIso_gen_Sig_EB_->Fill(pT_sum_gen);
1353             meEleISO_rel_chIso_gen_Sig_EB_->Fill(rel_pT_sum_gen);
1354           }
1355 
1356           for (long unsigned int j = 0; j < Ntracks_EB_list_Significance_Sig.size(); j++) {
1357             Ntracks_EB_list_Significance_Sig[j]->Fill(N_tracks_MTD_significance[j]);
1358             ch_iso_EB_list_Significance_Sig[j]->Fill(pT_sum_MTD_significance[j]);
1359             rel_ch_iso_EB_list_Significance_Sig[j]->Fill(rel_pT_sum_MTD_significance[j]);
1360 
1361             if (optionalPlots_) {
1362               Ntracks_sim_EB_list_Significance_Sig[j]->Fill(N_tracks_sim_MTD_significance[j]);
1363               ch_iso_sim_EB_list_Significance_Sig[j]->Fill(pT_sum_sim_MTD_significance[j]);
1364               rel_ch_iso_sim_EB_list_Significance_Sig[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1365             }
1366           }
1367 
1368           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1369             meEle_pt_noMTD_Sig_EB_->Fill(ele.pt());
1370             meEle_eta_noMTD_Sig_EB_->Fill(std::abs(ele.eta()));
1371             meEle_phi_noMTD_Sig_EB_->Fill(ele.phi());
1372           }
1373           if (optionalPlots_) {
1374             for (long unsigned int k = 0; k < Ntracks_EB_list_Sig.size(); k++) {
1375               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1376                 Ele_pT_MTD_EB_list_Sig[k]->Fill(ele.pt());
1377                 Ele_eta_MTD_EB_list_Sig[k]->Fill(std::abs(ele.eta()));
1378                 Ele_phi_MTD_EB_list_Sig[k]->Fill(ele.phi());
1379 
1380                 Ele_pT_sim_MTD_EB_list_Sig[k]->Fill(ele_sim_pt);
1381               }
1382             }
1383             if (rel_pT_sum_gen < rel_iso_cut_) {
1384               meEle_pt_gen_Sig_EB_->Fill(ele_sim_pt);
1385               meEle_eta_gen_Sig_EB_->Fill(ele_sim_eta);
1386               meEle_phi_gen_Sig_EB_->Fill(ele_sim_phi);
1387             }
1388           }
1389 
1390           for (long unsigned int k = 0; k < Ntracks_EB_list_Significance_Sig.size(); k++) {
1391             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1392               Ele_pT_MTD_EB_list_Significance_Sig[k]->Fill(ele.pt());
1393               Ele_eta_MTD_EB_list_Significance_Sig[k]->Fill(std::abs(ele.eta()));
1394               Ele_phi_MTD_EB_list_Significance_Sig[k]->Fill(ele.phi());
1395             }
1396             if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1397               Ele_pT_sim_MTD_EB_list_Significance_Sig[k]->Fill(ele_sim_pt);
1398           }
1399 
1400         } else {  // for endcap
1401 
1402           meEleISO_Ntracks_Sig_EE_->Fill(N_tracks_noMTD);
1403           meEleISO_chIso_Sig_EE_->Fill(pT_sum_noMTD);
1404           meEleISO_rel_chIso_Sig_EE_->Fill(rel_pT_sum_noMTD);
1405           if (optionalPlots_) {
1406             for (long unsigned int j = 0; j < Ntracks_EE_list_Sig.size(); j++) {
1407               Ntracks_EE_list_Sig[j]->Fill(N_tracks_MTD[j]);
1408               ch_iso_EE_list_Sig[j]->Fill(pT_sum_MTD[j]);
1409               rel_ch_iso_EE_list_Sig[j]->Fill(rel_pT_sum_MTD[j]);
1410 
1411               Ntracks_sim_EE_list_Sig[j]->Fill(N_tracks_sim_MTD[j]);
1412               ch_iso_sim_EE_list_Sig[j]->Fill(pT_sum_sim_MTD[j]);
1413               rel_ch_iso_sim_EE_list_Sig[j]->Fill(rel_pT_sum_sim_MTD[j]);
1414             }
1415             meEleISO_Ntracks_gen_Sig_EE_->Fill(N_tracks_gen);
1416             meEleISO_chIso_gen_Sig_EE_->Fill(pT_sum_gen);
1417             meEleISO_rel_chIso_gen_Sig_EE_->Fill(rel_pT_sum_gen);
1418           }
1419 
1420           for (long unsigned int j = 0; j < Ntracks_EE_list_Significance_Sig.size(); j++) {
1421             Ntracks_EE_list_Significance_Sig[j]->Fill(N_tracks_MTD_significance[j]);
1422             ch_iso_EE_list_Significance_Sig[j]->Fill(pT_sum_MTD_significance[j]);
1423             rel_ch_iso_EE_list_Significance_Sig[j]->Fill(rel_pT_sum_MTD_significance[j]);
1424 
1425             if (optionalPlots_) {
1426               Ntracks_sim_EE_list_Significance_Sig[j]->Fill(N_tracks_sim_MTD_significance[j]);
1427               ch_iso_sim_EE_list_Significance_Sig[j]->Fill(pT_sum_sim_MTD_significance[j]);
1428               rel_ch_iso_sim_EE_list_Significance_Sig[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1429             }
1430           }
1431 
1432           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1433             meEle_pt_noMTD_Sig_EE_->Fill(ele.pt());
1434             meEle_eta_noMTD_Sig_EE_->Fill(std::abs(ele.eta()));
1435             meEle_phi_noMTD_Sig_EE_->Fill(ele.phi());
1436           }
1437           if (optionalPlots_) {
1438             for (long unsigned int k = 0; k < Ntracks_EE_list_Sig.size(); k++) {
1439               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1440                 Ele_pT_MTD_EE_list_Sig[k]->Fill(ele.pt());
1441                 Ele_eta_MTD_EE_list_Sig[k]->Fill(std::abs(ele.eta()));
1442                 Ele_phi_MTD_EE_list_Sig[k]->Fill(ele.phi());
1443 
1444                 Ele_pT_sim_MTD_EE_list_Sig[k]->Fill(ele_sim_pt);
1445               }
1446             }
1447             if (rel_pT_sum_gen < rel_iso_cut_) {
1448               meEle_pt_gen_Sig_EE_->Fill(ele_sim_pt);
1449               meEle_eta_gen_Sig_EE_->Fill(ele_sim_eta);
1450               meEle_phi_gen_Sig_EE_->Fill(ele_sim_phi);
1451             }
1452           }
1453           for (long unsigned int k = 0; k < Ntracks_EE_list_Significance_Sig.size(); k++) {
1454             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1455               Ele_pT_MTD_EE_list_Significance_Sig[k]->Fill(ele.pt());
1456               Ele_eta_MTD_EE_list_Significance_Sig[k]->Fill(std::abs(ele.eta()));
1457               Ele_phi_MTD_EE_list_Significance_Sig[k]->Fill(ele.phi());
1458 
1459               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1460                 Ele_pT_sim_MTD_EE_list_Significance_Sig[k]->Fill(ele_sim_pt);
1461             }
1462           }
1463         }
1464       } else {  // non-promt part
1465         if (Barrel_ele) {
1466           meEleISO_Ntracks_Bkg_EB_->Fill(N_tracks_noMTD);
1467           meEleISO_chIso_Bkg_EB_->Fill(pT_sum_noMTD);
1468           meEleISO_rel_chIso_Bkg_EB_->Fill(rel_pT_sum_noMTD);
1469           if (optionalPlots_) {
1470             for (long unsigned int j = 0; j < Ntracks_EB_list_Bkg.size(); j++) {
1471               Ntracks_EB_list_Bkg[j]->Fill(N_tracks_MTD[j]);
1472               ch_iso_EB_list_Bkg[j]->Fill(pT_sum_MTD[j]);
1473               rel_ch_iso_EB_list_Bkg[j]->Fill(rel_pT_sum_MTD[j]);
1474 
1475               Ntracks_sim_EB_list_Bkg[j]->Fill(N_tracks_sim_MTD[j]);
1476               ch_iso_sim_EB_list_Bkg[j]->Fill(pT_sum_sim_MTD[j]);
1477               rel_ch_iso_sim_EB_list_Bkg[j]->Fill(rel_pT_sum_sim_MTD[j]);
1478             }
1479             meEleISO_Ntracks_gen_Bkg_EB_->Fill(N_tracks_gen);
1480             meEleISO_chIso_gen_Bkg_EB_->Fill(pT_sum_gen);
1481             meEleISO_rel_chIso_gen_Bkg_EB_->Fill(rel_pT_sum_gen);
1482           }
1483 
1484           for (long unsigned int j = 0; j < Ntracks_EB_list_Significance_Bkg.size(); j++) {
1485             Ntracks_EB_list_Significance_Bkg[j]->Fill(N_tracks_MTD_significance[j]);
1486             ch_iso_EB_list_Significance_Bkg[j]->Fill(pT_sum_MTD_significance[j]);
1487             rel_ch_iso_EB_list_Significance_Bkg[j]->Fill(rel_pT_sum_MTD_significance[j]);
1488 
1489             if (optionalPlots_) {
1490               Ntracks_sim_EB_list_Significance_Bkg[j]->Fill(N_tracks_sim_MTD_significance[j]);
1491               ch_iso_sim_EB_list_Significance_Bkg[j]->Fill(pT_sum_sim_MTD_significance[j]);
1492               rel_ch_iso_sim_EB_list_Significance_Bkg[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1493             }
1494           }
1495 
1496           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1497             meEle_pt_noMTD_Bkg_EB_->Fill(ele.pt());
1498             meEle_eta_noMTD_Bkg_EB_->Fill(std::abs(ele.eta()));
1499             meEle_phi_noMTD_Bkg_EB_->Fill(ele.phi());
1500           }
1501           if (optionalPlots_) {
1502             for (long unsigned int k = 0; k < Ntracks_EB_list_Bkg.size(); k++) {
1503               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1504                 Ele_pT_MTD_EB_list_Bkg[k]->Fill(ele.pt());
1505                 Ele_eta_MTD_EB_list_Bkg[k]->Fill(std::abs(ele.eta()));
1506                 Ele_phi_MTD_EB_list_Bkg[k]->Fill(ele.phi());
1507 
1508                 Ele_pT_sim_MTD_EB_list_Bkg[k]->Fill(ele_sim_pt);
1509               }
1510             }
1511             if (rel_pT_sum_gen < rel_iso_cut_) {
1512               meEle_pt_gen_Bkg_EB_->Fill(ele_sim_pt);
1513               meEle_eta_gen_Bkg_EB_->Fill(ele_sim_eta);
1514               meEle_phi_gen_Bkg_EB_->Fill(ele_sim_phi);
1515             }
1516           }
1517           for (long unsigned int k = 0; k < Ntracks_EB_list_Significance_Bkg.size(); k++) {
1518             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1519               Ele_pT_MTD_EB_list_Significance_Bkg[k]->Fill(ele.pt());
1520               Ele_eta_MTD_EB_list_Significance_Bkg[k]->Fill(std::abs(ele.eta()));
1521               Ele_phi_MTD_EB_list_Significance_Bkg[k]->Fill(ele.phi());
1522 
1523               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1524                 Ele_pT_sim_MTD_EB_list_Significance_Bkg[k]->Fill(ele_sim_pt);
1525             }
1526           }
1527 
1528         } else {  // for endcap
1529           meEleISO_Ntracks_Bkg_EE_->Fill(N_tracks_noMTD);
1530           meEleISO_chIso_Bkg_EE_->Fill(pT_sum_noMTD);
1531           meEleISO_rel_chIso_Bkg_EE_->Fill(rel_pT_sum_noMTD);
1532           if (optionalPlots_) {
1533             for (long unsigned int j = 0; j < Ntracks_EE_list_Bkg.size(); j++) {
1534               Ntracks_EE_list_Bkg[j]->Fill(N_tracks_MTD[j]);
1535               ch_iso_EE_list_Bkg[j]->Fill(pT_sum_MTD[j]);
1536               rel_ch_iso_EE_list_Bkg[j]->Fill(rel_pT_sum_MTD[j]);
1537 
1538               Ntracks_sim_EE_list_Bkg[j]->Fill(N_tracks_sim_MTD[j]);
1539               ch_iso_sim_EE_list_Bkg[j]->Fill(pT_sum_sim_MTD[j]);
1540               rel_ch_iso_sim_EE_list_Bkg[j]->Fill(rel_pT_sum_sim_MTD[j]);
1541             }
1542             meEleISO_Ntracks_gen_Bkg_EE_->Fill(N_tracks_gen);
1543             meEleISO_chIso_gen_Bkg_EE_->Fill(pT_sum_gen);
1544             meEleISO_rel_chIso_gen_Bkg_EE_->Fill(rel_pT_sum_gen);
1545           }
1546 
1547           for (long unsigned int j = 0; j < Ntracks_EE_list_Significance_Bkg.size(); j++) {
1548             Ntracks_EE_list_Significance_Bkg[j]->Fill(N_tracks_MTD_significance[j]);
1549             ch_iso_EE_list_Significance_Bkg[j]->Fill(pT_sum_MTD_significance[j]);
1550             rel_ch_iso_EE_list_Significance_Bkg[j]->Fill(rel_pT_sum_MTD_significance[j]);
1551 
1552             if (optionalPlots_) {
1553               Ntracks_sim_EE_list_Significance_Bkg[j]->Fill(N_tracks_sim_MTD_significance[j]);
1554               ch_iso_sim_EE_list_Significance_Bkg[j]->Fill(pT_sum_sim_MTD_significance[j]);
1555               rel_ch_iso_sim_EE_list_Significance_Bkg[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1556             }
1557           }
1558 
1559           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1560             meEle_pt_noMTD_Bkg_EE_->Fill(ele.pt());
1561             meEle_eta_noMTD_Bkg_EE_->Fill(std::abs(ele.eta()));
1562             meEle_phi_noMTD_Bkg_EE_->Fill(ele.phi());
1563           }
1564           if (optionalPlots_) {
1565             for (long unsigned int k = 0; k < Ntracks_EE_list_Bkg.size(); k++) {
1566               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1567                 Ele_pT_MTD_EE_list_Bkg[k]->Fill(ele.pt());
1568                 Ele_eta_MTD_EE_list_Bkg[k]->Fill(std::abs(ele.eta()));
1569                 Ele_phi_MTD_EE_list_Bkg[k]->Fill(ele.phi());
1570 
1571                 Ele_pT_sim_MTD_EE_list_Bkg[k]->Fill(ele_sim_pt);
1572               }
1573             }
1574             if (rel_pT_sum_gen < rel_iso_cut_) {
1575               meEle_pt_gen_Bkg_EE_->Fill(ele_sim_pt);
1576               meEle_eta_gen_Bkg_EE_->Fill(ele_sim_eta);
1577               meEle_phi_gen_Bkg_EE_->Fill(ele_sim_phi);
1578             }
1579           }
1580 
1581           for (long unsigned int k = 0; k < Ntracks_EE_list_Significance_Bkg.size(); k++) {
1582             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1583               Ele_pT_MTD_EE_list_Significance_Bkg[k]->Fill(ele.pt());
1584               Ele_eta_MTD_EE_list_Significance_Bkg[k]->Fill(std::abs(ele.eta()));
1585               Ele_phi_MTD_EE_list_Significance_Bkg[k]->Fill(ele.phi());
1586 
1587               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1588                 Ele_pT_sim_MTD_EE_list_Significance_Bkg[k]->Fill(ele_sim_pt);
1589             }
1590           }
1591         }
1592       }
1593     }  // electron matched to a track
1594   }    // electron collection inside single event
1595 }
1596 
1597 // ------------ method for histogram booking ------------
1598 void MtdEleIsoValidation::bookHistograms(DQMStore::IBooker& ibook, edm::Run const& run, edm::EventSetup const& iSetup) {
1599   ibook.setCurrentFolder(folder_);
1600 
1601   // for regular Validation use a reduced binning, for detailed analysis and ROC curves use the larger one
1602   int nbin_1 = 40;
1603   int nbin_2 = 40;
1604   if (optionalPlots_) {
1605     nbin_1 = 1000;
1606     nbin_2 = 2000;
1607   }
1608 
1609   // histogram booking
1610 
1611   meEle_avg_error_SigTrk_check_ =
1612       ibook.book1D("SigTrk_avg_timing_err",
1613                    "Average signal electron track MTD timing uncertainty;Time Error (ns);Counts",
1614                    200,
1615                    0,
1616                    0.1);
1617   meEle_avg_error_PUTrk_check_ = ibook.book1D(
1618       "PUTrk_avg_timing_err", "Average PU track MTD timing uncertainty;Time Error (ns);Counts", 200, 0, 0.1);
1619   meEle_avg_error_vtx_check_ =
1620       ibook.book1D("Vtx_avg_timing_err", "Average vertex timing uncertainty;Time Error (ns);Counts", 200, 0, 0.1);
1621 
1622   meEle_no_dt_check_ =
1623       ibook.book1D("Track_dt_info_check",
1624                    "Tracks dt check - ratio between tracks with (value 1) and without (value 0) timing info",
1625                    2,
1626                    0,
1627                    2);
1628 
1629   meTrk_genMatch_check_ = ibook.book1D(
1630       "Track_genMatch_info_check", "Check on tracks matched with a GenParticle (matched=1, non matched=0)", 2, 0, 2);
1631 
1632   // signal
1633   meEleISO_Ntracks_Sig_EB_ = ibook.book1D("Ele_Iso_Ntracks_Sig_EB",
1634                                           "Number of tracks in isolation cone around electron track after basic cuts - "
1635                                           "Signal Barrel;Number of tracks;Counts",
1636                                           20,
1637                                           0,
1638                                           20);
1639 
1640   meEleISO_chIso_Sig_EB_ = ibook.book1D(
1641       "Ele_chIso_sum_Sig_EB",
1642       "Track pT sum in isolation cone around electron track after basic cuts - Signal Barrel;p_{T} (GeV);Counts",
1643       nbin_2,
1644       0,
1645       20);
1646 
1647   meEleISO_rel_chIso_Sig_EB_ = ibook.book1D(
1648       "Ele_rel_chIso_sum_Sig_EB",
1649       "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Barrel;Isolation;Counts",
1650       nbin_1,
1651       0,
1652       4);
1653   if (optionalPlots_) {
1654     meEleISO_Ntracks_MTD_1_Sig_EB_ =
1655         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Sig_EB",
1656                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1657                      20,
1658                      0,
1659                      20);
1660 
1661     meEleISO_chIso_MTD_1_Sig_EB_ = ibook.book1D(
1662         "Ele_chIso_sum_MTD_1_Sig_EB",
1663         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1664         nbin_2,
1665         0,
1666         20);
1667     meEleISO_rel_chIso_MTD_1_Sig_EB_ = ibook.book1D(
1668         "Ele_rel_chIso_sum_MTD_1_Sig_EB",
1669         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1670         nbin_1,
1671         0,
1672         4);
1673     // gen
1674     meEleISO_Ntracks_gen_Sig_EB_ = ibook.book1D("Ele_Iso_Ntracks_gen_Sig_EB",
1675                                                 "Number of tracks in isolation cone around electron track after basic "
1676                                                 "cuts using genInfo - Signal Barrel;Number of tracks;Counts",
1677                                                 20,
1678                                                 0,
1679                                                 20);
1680 
1681     meEleISO_chIso_gen_Sig_EB_ = ibook.book1D("Ele_chIso_sum_gen_Sig_EB",
1682                                               "Track pT sum in isolation cone around electron track after basic cuts "
1683                                               "using genInfo - Signal Barrel;p_{T} (GeV);Counts",
1684                                               nbin_2,
1685                                               0,
1686                                               20);
1687 
1688     meEleISO_rel_chIso_gen_Sig_EB_ = ibook.book1D("Ele_rel_chIso_sum_gen_Sig_EB",
1689                                                   "Track relative pT sum in isolation cone around electron track after "
1690                                                   "basic cuts using genInfo - Signal Barrel;Isolation;Counts",
1691                                                   nbin_1,
1692                                                   0,
1693                                                   4);
1694 
1695     meEleISO_Ntracks_MTD_2_Sig_EB_ =
1696         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Sig_EB",
1697                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1698                      20,
1699                      0,
1700                      20);
1701 
1702     meEleISO_chIso_MTD_2_Sig_EB_ = ibook.book1D(
1703         "Ele_chIso_sum_MTD_2_Sig_EB",
1704         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1705         nbin_2,
1706         0,
1707         20);
1708     meEleISO_rel_chIso_MTD_2_Sig_EB_ = ibook.book1D(
1709         "Ele_rel_chIso_sum_MTD_2_Sig_EB",
1710         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1711         nbin_1,
1712         0,
1713         4);
1714 
1715     meEleISO_Ntracks_MTD_3_Sig_EB_ =
1716         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Sig_EB",
1717                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1718                      20,
1719                      0,
1720                      20);
1721     meEleISO_chIso_MTD_3_Sig_EB_ = ibook.book1D(
1722         "Ele_chIso_sum_MTD_3_Sig_EB",
1723         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1724         nbin_2,
1725         0,
1726         20);
1727     meEleISO_rel_chIso_MTD_3_Sig_EB_ = ibook.book1D(
1728         "Ele_rel_chIso_sum_MTD_3_Sig_EB",
1729         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1730         nbin_1,
1731         0,
1732         4);
1733 
1734     meEleISO_Ntracks_MTD_4_Sig_EB_ =
1735         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Sig_EB",
1736                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1737                      20,
1738                      0,
1739                      20);
1740     meEleISO_chIso_MTD_4_Sig_EB_ = ibook.book1D(
1741         "Ele_chIso_sum_MTD_4_Sig_EB",
1742         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1743         nbin_2,
1744         0,
1745         20);
1746     meEleISO_rel_chIso_MTD_4_Sig_EB_ = ibook.book1D(
1747         "Ele_rel_chIso_sum_MTD_4_Sig_EB",
1748         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1749         nbin_1,
1750         0,
1751         4);
1752 
1753     meEleISO_Ntracks_MTD_5_Sig_EB_ =
1754         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Sig_EB",
1755                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1756                      20,
1757                      0,
1758                      20);
1759     meEleISO_chIso_MTD_5_Sig_EB_ = ibook.book1D(
1760         "Ele_chIso_sum_MTD_5_Sig_EB",
1761         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1762         nbin_2,
1763         0,
1764         20);
1765     meEleISO_rel_chIso_MTD_5_Sig_EB_ = ibook.book1D(
1766         "Ele_rel_chIso_sum_MTD_5_Sig_EB",
1767         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1768         nbin_1,
1769         0,
1770         4);
1771 
1772     meEleISO_Ntracks_MTD_6_Sig_EB_ =
1773         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Sig_EB",
1774                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1775                      20,
1776                      0,
1777                      20);
1778     meEleISO_chIso_MTD_6_Sig_EB_ = ibook.book1D(
1779         "Ele_chIso_sum_MTD_6_Sig_EB",
1780         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1781         nbin_2,
1782         0,
1783         20);
1784     meEleISO_rel_chIso_MTD_6_Sig_EB_ = ibook.book1D(
1785         "Ele_rel_chIso_sum_MTD_6_Sig_EB",
1786         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1787         nbin_1,
1788         0,
1789         4);
1790 
1791     meEleISO_Ntracks_MTD_7_Sig_EB_ =
1792         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Sig_EB",
1793                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1794                      20,
1795                      0,
1796                      20);
1797     meEleISO_chIso_MTD_7_Sig_EB_ = ibook.book1D(
1798         "Ele_chIso_sum_MTD_7_Sig_EB",
1799         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1800         nbin_2,
1801         0,
1802         20);
1803     meEleISO_rel_chIso_MTD_7_Sig_EB_ = ibook.book1D(
1804         "Ele_rel_chIso_sum_MTD_7_Sig_EB",
1805         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1806         nbin_1,
1807         0,
1808         4);
1809 
1810     meEleISO_Ntracks_MTD_sim_1_Sig_EB_ =
1811         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Sig_EB",
1812                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1813                      20,
1814                      0,
1815                      20);
1816 
1817     meEleISO_chIso_MTD_sim_1_Sig_EB_ = ibook.book1D(
1818         "Ele_chIso_sum_MTD_sim_1_Sig_EB",
1819         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1820         nbin_2,
1821         0,
1822         20);
1823     meEleISO_rel_chIso_MTD_sim_1_Sig_EB_ = ibook.book1D(
1824         "Ele_rel_chIso_sum_MTD_sim_1_Sig_EB",
1825         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1826         nbin_1,
1827         0,
1828         4);
1829 
1830     meEleISO_Ntracks_MTD_sim_2_Sig_EB_ =
1831         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Sig_EB",
1832                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1833                      20,
1834                      0,
1835                      20);
1836 
1837     meEleISO_chIso_MTD_sim_2_Sig_EB_ = ibook.book1D(
1838         "Ele_chIso_sum_MTD_sim_2_Sig_EB",
1839         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1840         nbin_2,
1841         0,
1842         20);
1843     meEleISO_rel_chIso_MTD_sim_2_Sig_EB_ = ibook.book1D(
1844         "Ele_rel_chIso_sum_MTD_sim_2_Sig_EB",
1845         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1846         nbin_1,
1847         0,
1848         4);
1849 
1850     meEleISO_Ntracks_MTD_sim_3_Sig_EB_ =
1851         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Sig_EB",
1852                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1853                      20,
1854                      0,
1855                      20);
1856     meEleISO_chIso_MTD_sim_3_Sig_EB_ = ibook.book1D(
1857         "Ele_chIso_sum_MTD_sim_3_Sig_EB",
1858         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1859         nbin_2,
1860         0,
1861         20);
1862     meEleISO_rel_chIso_MTD_sim_3_Sig_EB_ = ibook.book1D(
1863         "Ele_rel_chIso_sum_MTD_sim_3_Sig_EB",
1864         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1865         nbin_1,
1866         0,
1867         4);
1868 
1869     meEleISO_Ntracks_MTD_sim_4_Sig_EB_ =
1870         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Sig_EB",
1871                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1872                      20,
1873                      0,
1874                      20);
1875     meEleISO_chIso_MTD_sim_4_Sig_EB_ = ibook.book1D(
1876         "Ele_chIso_sum_MTD_sim_4_Sig_EB",
1877         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1878         nbin_2,
1879         0,
1880         20);
1881     meEleISO_rel_chIso_MTD_sim_4_Sig_EB_ = ibook.book1D(
1882         "Ele_rel_chIso_sum_MTD_sim_4_Sig_EB",
1883         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1884         nbin_1,
1885         0,
1886         4);
1887 
1888     meEleISO_Ntracks_MTD_sim_5_Sig_EB_ =
1889         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Sig_EB",
1890                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1891                      20,
1892                      0,
1893                      20);
1894     meEleISO_chIso_MTD_sim_5_Sig_EB_ = ibook.book1D(
1895         "Ele_chIso_sum_MTD_sim_5_Sig_EB",
1896         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1897         nbin_2,
1898         0,
1899         20);
1900     meEleISO_rel_chIso_MTD_sim_5_Sig_EB_ = ibook.book1D(
1901         "Ele_rel_chIso_sum_MTD_sim_5_Sig_EB",
1902         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1903         nbin_1,
1904         0,
1905         4);
1906 
1907     meEleISO_Ntracks_MTD_sim_6_Sig_EB_ =
1908         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Sig_EB",
1909                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1910                      20,
1911                      0,
1912                      20);
1913     meEleISO_chIso_MTD_sim_6_Sig_EB_ = ibook.book1D(
1914         "Ele_chIso_sum_MTD_sim_6_Sig_EB",
1915         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1916         nbin_2,
1917         0,
1918         20);
1919     meEleISO_rel_chIso_MTD_sim_6_Sig_EB_ = ibook.book1D(
1920         "Ele_rel_chIso_sum_MTD_sim_6_Sig_EB",
1921         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1922         nbin_1,
1923         0,
1924         4);
1925 
1926     meEleISO_Ntracks_MTD_sim_7_Sig_EB_ =
1927         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Sig_EB",
1928                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1929                      20,
1930                      0,
1931                      20);
1932     meEleISO_chIso_MTD_sim_7_Sig_EB_ = ibook.book1D(
1933         "Ele_chIso_sum_MTD_sim_7_Sig_EB",
1934         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1935         nbin_2,
1936         0,
1937         20);
1938     meEleISO_rel_chIso_MTD_sim_7_Sig_EB_ = ibook.book1D(
1939         "Ele_rel_chIso_sum_MTD_sim_7_Sig_EB",
1940         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1941         nbin_1,
1942         0,
1943         4);
1944   }
1945   meEleISO_Ntracks_MTD_4sigma_Sig_EB_ =
1946       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Sig_EB",
1947                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma compatibiliy - "
1948                    "Signal Barrel;Number of tracks;Counts",
1949                    20,
1950                    0,
1951                    20);
1952   meEleISO_chIso_MTD_4sigma_Sig_EB_ =
1953       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Sig_EB",
1954                    "Track pT sum in isolation cone around electron track after basic "
1955                    "cuts with MTD - 4 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
1956                    nbin_2,
1957                    0,
1958                    20);
1959   meEleISO_rel_chIso_MTD_4sigma_Sig_EB_ =
1960       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Sig_EB",
1961                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
1962                    "compatibiliy - Signal Barrel;Isolation;Counts",
1963                    nbin_1,
1964                    0,
1965                    4);
1966 
1967   meEleISO_Ntracks_MTD_3sigma_Sig_EB_ =
1968       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Sig_EB",
1969                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma compatibiliy - "
1970                    "Signal Barrel;Number of tracks;Counts",
1971                    20,
1972                    0,
1973                    20);
1974   meEleISO_chIso_MTD_3sigma_Sig_EB_ =
1975       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Sig_EB",
1976                    "Track pT sum in isolation cone around electron track after basic "
1977                    "cuts with MTD - 3 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
1978                    nbin_2,
1979                    0,
1980                    20);
1981   meEleISO_rel_chIso_MTD_3sigma_Sig_EB_ = ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Sig_EB",
1982                                                        "Track relative pT sum in isolation cone around electron track "
1983                                                        "after basic cuts with MTD - 3 sigma;Isolation;Counts"
1984                                                        "compatibiliy - Signal Barrel",
1985                                                        nbin_1,
1986                                                        0,
1987                                                        4);
1988 
1989   meEleISO_Ntracks_MTD_2sigma_Sig_EB_ =
1990       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Sig_EB",
1991                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibiliy - "
1992                    "Signal Barrel;Number of tracks;Counts",
1993                    20,
1994                    0,
1995                    20);
1996   meEleISO_chIso_MTD_2sigma_Sig_EB_ =
1997       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Sig_EB",
1998                    "Track pT sum in isolation cone around electron track after basic "
1999                    "cuts with MTD - 2 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2000                    nbin_2,
2001                    0,
2002                    20);
2003   meEleISO_rel_chIso_MTD_2sigma_Sig_EB_ =
2004       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Sig_EB",
2005                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2006                    "compatibiliy - Signal Barrel;Isolation;Counts",
2007                    nbin_1,
2008                    0,
2009                    4);
2010 
2011   meEle_pt_tot_Sig_EB_ =
2012       ibook.book1D("Ele_pT_tot_Sig_EB", "Electron pT tot - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2013   meEle_pt_noMTD_Sig_EB_ =
2014       ibook.book1D("Ele_pT_noMTD_Sig_EB", "Electron pT noMTD - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2015 
2016   meEle_pt_sim_tot_Sig_EB_ =
2017       ibook.book1D("Ele_pT_sim_tot_Sig_EB", "Electron SIM pT tot - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2018 
2019   meEle_eta_tot_Sig_EB_ =
2020       ibook.book1D("Ele_eta_tot_Sig_EB", "Electron eta tot - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2021   meEle_eta_noMTD_Sig_EB_ =
2022       ibook.book1D("Ele_eta_noMTD_Sig_EB", "Electron eta noMTD - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2023 
2024   meEle_phi_tot_Sig_EB_ =
2025       ibook.book1D("Ele_phi_tot_Sig_EB", "Electron phi tot - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2026   meEle_phi_noMTD_Sig_EB_ =
2027       ibook.book1D("Ele_phi_noMTD_Sig_EB", "Electron phi noMTD - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2028 
2029   if (optionalPlots_) {
2030     meEleISO_Ntracks_MTD_sim_4sigma_Sig_EB_ =
2031         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Sig_EB",
2032                      "Number of tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2033                      "compatibiliy - Signal Barrel;Number of tracks;Counts",
2034                      20,
2035                      0,
2036                      20);
2037     meEleISO_chIso_MTD_sim_4sigma_Sig_EB_ =
2038         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Sig_EB",
2039                      "Track pT sum in isolation cone around electron track after "
2040                      "basic cuts with MTD - 4 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2041                      nbin_2,
2042                      0,
2043                      20);
2044     meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EB_ = ibook.book1D(
2045         "Ele_rel_chIso_sum_MTD_sim_4sigma_Sig_EB",
2046         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2047         "compatibiliy - Signal Barrel;Isolation;Counts",
2048         nbin_1,
2049         0,
2050         4);
2051 
2052     meEleISO_Ntracks_MTD_sim_3sigma_Sig_EB_ =
2053         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Sig_EB",
2054                      "Tracks in isolation cone around electron track after basic "
2055                      "cuts with MTD  - 3 sigma compatibiliy - Signal Barrel;Number of tracks;Counts",
2056                      20,
2057                      0,
2058                      20);
2059     meEleISO_chIso_MTD_sim_3sigma_Sig_EB_ =
2060         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Sig_EB",
2061                      "Track pT sum in isolation cone around electron track after "
2062                      "basic cuts with MTD - 3 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2063                      nbin_2,
2064                      0,
2065                      20);
2066     meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EB_ = ibook.book1D(
2067         "Ele_rel_chIso_sum_MTD_sim_3sigma_Sig_EB",
2068         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 3 sigma "
2069         "compatibiliy - Signal Barrel;Isolation;Counts",
2070         nbin_1,
2071         0,
2072         4);
2073 
2074     meEleISO_Ntracks_MTD_sim_2sigma_Sig_EB_ = ibook.book1D(
2075         "Ele_Iso_Ntracks_MTD_sim_2sigma_Sig_EB",
2076         "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibiliy - "
2077         "Signal Barrel;Number of tracks;Counts",
2078         20,
2079         0,
2080         20);
2081     meEleISO_chIso_MTD_sim_2sigma_Sig_EB_ =
2082         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Sig_EB",
2083                      "Track pT sum in isolation cone around electron track after "
2084                      "basic cuts with MTD - 2 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2085                      nbin_2,
2086                      0,
2087                      20);
2088     meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EB_ = ibook.book1D(
2089         "Ele_rel_chIso_sum_MTD_sim_2sigma_Sig_EB",
2090         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2091         "compatibiliy - Signal Barrel;Isolation;Counts",
2092         nbin_1,
2093         0,
2094         4);
2095 
2096     meEle_pt_gen_Sig_EB_ =
2097         ibook.book1D("Ele_pT_gen_Sig_EB", "Electron pT genInfo - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2098     meEle_eta_gen_Sig_EB_ =
2099         ibook.book1D("Ele_eta_gen_Sig_EB", "Electron eta genInfo - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2100     meEle_phi_gen_Sig_EB_ =
2101         ibook.book1D("Ele_phi_gen_Sig_EB", "Electron phi genInfo - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2102 
2103     meEle_pt_MTD_1_Sig_EB_ = ibook.book1D("Ele_pT_MTD_1_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2104     meEle_eta_MTD_1_Sig_EB_ = ibook.book1D("Ele_eta_MTD_1_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2105     meEle_phi_MTD_1_Sig_EB_ = ibook.book1D("Ele_phi_MTD_1_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2106 
2107     meEle_pt_MTD_2_Sig_EB_ = ibook.book1D("Ele_pT_MTD_2_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2108     meEle_eta_MTD_2_Sig_EB_ = ibook.book1D("Ele_eta_MTD_2_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2109     meEle_phi_MTD_2_Sig_EB_ = ibook.book1D("Ele_phi_MTD_2_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2110 
2111     meEle_pt_MTD_3_Sig_EB_ = ibook.book1D("Ele_pT_MTD_3_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2112     meEle_eta_MTD_3_Sig_EB_ = ibook.book1D("Ele_eta_MTD_3_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2113     meEle_phi_MTD_3_Sig_EB_ = ibook.book1D("Ele_phi_MTD_3_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2114 
2115     meEle_pt_MTD_4_Sig_EB_ = ibook.book1D("Ele_pT_MTD_4_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2116     meEle_eta_MTD_4_Sig_EB_ = ibook.book1D("Ele_eta_MTD_4_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2117     meEle_phi_MTD_4_Sig_EB_ = ibook.book1D("Ele_phi_MTD_4_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2118 
2119     meEle_pt_MTD_5_Sig_EB_ = ibook.book1D("Ele_pT_MTD_5_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2120     meEle_eta_MTD_5_Sig_EB_ = ibook.book1D("Ele_eta_MTD_5_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2121     meEle_phi_MTD_5_Sig_EB_ = ibook.book1D("Ele_phi_MTD_5_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2122 
2123     meEle_pt_MTD_6_Sig_EB_ = ibook.book1D("Ele_pT_MTD_6_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2124     meEle_eta_MTD_6_Sig_EB_ = ibook.book1D("Ele_eta_MTD_6_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2125     meEle_phi_MTD_6_Sig_EB_ = ibook.book1D("Ele_phi_MTD_6_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2126 
2127     meEle_pt_MTD_7_Sig_EB_ = ibook.book1D("Ele_pT_MTD_7_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2128     meEle_eta_MTD_7_Sig_EB_ = ibook.book1D("Ele_eta_MTD_7_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2129     meEle_phi_MTD_7_Sig_EB_ = ibook.book1D("Ele_phi_MTD_7_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2130 
2131     meEle_pt_sim_MTD_1_Sig_EB_ =
2132         ibook.book1D("Ele_pT_sim_MTD_1_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2133     meEle_pt_sim_MTD_2_Sig_EB_ =
2134         ibook.book1D("Ele_pT_sim_MTD_2_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2135     meEle_pt_sim_MTD_3_Sig_EB_ =
2136         ibook.book1D("Ele_pT_sim_MTD_3_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2137     meEle_pt_sim_MTD_4_Sig_EB_ =
2138         ibook.book1D("Ele_pT_sim_MTD_4_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2139     meEle_pt_sim_MTD_5_Sig_EB_ =
2140         ibook.book1D("Ele_pT_sim_MTD_5_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2141     meEle_pt_sim_MTD_6_Sig_EB_ =
2142         ibook.book1D("Ele_pT_sim_MTD_6_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2143     meEle_pt_sim_MTD_7_Sig_EB_ =
2144         ibook.book1D("Ele_pT_sim_MTD_7_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2145   }
2146 
2147   meEle_pt_MTD_4sigma_Sig_EB_ =
2148       ibook.book1D("Ele_pT_MTD_4sigma_Sig_EB",
2149                    "Electron pT MTD - 4 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2150                    30,
2151                    10,
2152                    100);
2153   meEle_eta_MTD_4sigma_Sig_EB_ = ibook.book1D(
2154       "Ele_eta_MTD_4sigma_Sig_EB", "Electron eta MTD - 4 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2155   meEle_phi_MTD_4sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_4sigma_Sig_EB",
2156                                               "Electron phi MTD - 4 sigma compatibility - Signal Barrel;#phi;Counts",
2157                                               64,
2158                                               -3.2,
2159                                               3.2);
2160 
2161   meEle_pt_MTD_3sigma_Sig_EB_ =
2162       ibook.book1D("Ele_pT_MTD_3sigma_Sig_EB",
2163                    "Electron pT MTD - 3 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2164                    30,
2165                    10,
2166                    100);
2167   meEle_eta_MTD_3sigma_Sig_EB_ = ibook.book1D(
2168       "Ele_eta_MTD_3sigma_Sig_EB", "Electron eta MTD - 3 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2169   meEle_phi_MTD_3sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_3sigma_Sig_EB",
2170                                               "Electron phi MTD - 3 sigma compatibility - Signal Barrel;#phi;Counts",
2171                                               64,
2172                                               -3.2,
2173                                               3.2);
2174 
2175   meEle_pt_MTD_2sigma_Sig_EB_ =
2176       ibook.book1D("Ele_pT_MTD_2sigma_Sig_EB",
2177                    "Electron pT MTD - 2 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2178                    30,
2179                    10,
2180                    100);
2181   meEle_eta_MTD_2sigma_Sig_EB_ = ibook.book1D(
2182       "Ele_eta_MTD_2sigma_Sig_EB", "Electron eta MTD - 2 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2183   meEle_phi_MTD_2sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_2sigma_Sig_EB",
2184                                               "Electron phi MTD - 2 sigma compatibility - Signal Barrel;#phi;Counts",
2185                                               64,
2186                                               -3.2,
2187                                               3.2);
2188 
2189   meEleISO_Ntracks_Sig_EE_ = ibook.book1D("Ele_Iso_Ntracks_Sig_EE",
2190                                           "Number of tracks in isolation cone around electron track after basic cuts - "
2191                                           "Signal Endcap;Number of tracks;Counts",
2192                                           20,
2193                                           0,
2194                                           20);
2195   meEleISO_chIso_Sig_EE_ = ibook.book1D(
2196       "Ele_chIso_sum_Sig_EE",
2197       "Track pT sum in isolation cone around electron track after basic cuts - Signal Endcap;p_{T} (GeV);Counts",
2198       nbin_2,
2199       0,
2200       20);
2201   meEleISO_rel_chIso_Sig_EE_ = ibook.book1D(
2202       "Ele_rel_chIso_sum_Sig_EE",
2203       "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Endcap;Isolation;Counts",
2204       nbin_1,
2205       0,
2206       4);
2207 
2208   if (optionalPlots_) {
2209     meEle_pt_sim_MTD_4sigma_Sig_EB_ =
2210         ibook.book1D("Ele_pT_sim_MTD_4sigma_Sig_EB",
2211                      "Electron pT MTD SIM - 4 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2212                      30,
2213                      10,
2214                      100);
2215     meEle_pt_sim_MTD_3sigma_Sig_EB_ =
2216         ibook.book1D("Ele_pT_sim_MTD_3sigma_Sig_EB",
2217                      "Electron pT MTD SIM - 3 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2218                      30,
2219                      10,
2220                      100);
2221     meEle_pt_sim_MTD_2sigma_Sig_EB_ =
2222         ibook.book1D("Ele_pT_sim_MTD_2sigma_Sig_EB",
2223                      "Electron pT MTD SIM - 2 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2224                      30,
2225                      10,
2226                      100);
2227 
2228     meEleISO_Ntracks_MTD_1_Sig_EE_ =
2229         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Sig_EE",
2230                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2231                      20,
2232                      0,
2233                      20);
2234     meEleISO_chIso_MTD_1_Sig_EE_ = ibook.book1D(
2235         "Ele_chIso_sum_MTD_1_Sig_EE",
2236         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2237         nbin_2,
2238         0,
2239         20);
2240     meEleISO_rel_chIso_MTD_1_Sig_EE_ = ibook.book1D(
2241         "Ele_rel_chIso_sum_MTD_1_Sig_EE",
2242         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2243         nbin_1,
2244         0,
2245         4);
2246 
2247     meEleISO_Ntracks_MTD_2_Sig_EE_ =
2248         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Sig_EE",
2249                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2250                      20,
2251                      0,
2252                      20);
2253     meEleISO_chIso_MTD_2_Sig_EE_ = ibook.book1D(
2254         "Ele_chIso_sum_MTD_2_Sig_EE",
2255         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2256         nbin_2,
2257         0,
2258         20);
2259     meEleISO_rel_chIso_MTD_2_Sig_EE_ = ibook.book1D(
2260         "Ele_rel_chIso_sum_MTD_2_Sig_EE",
2261         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2262         nbin_1,
2263         0,
2264         4);
2265 
2266     meEleISO_Ntracks_gen_Sig_EE_ =
2267         ibook.book1D("Ele_Iso_Ntracks_gen_Sig_EE",
2268                      "Tracks in isolation cone around electron track after basic cuts using genInfo - Signal Endcap",
2269                      20,
2270                      0,
2271                      20);
2272     meEleISO_chIso_gen_Sig_EE_ =
2273         ibook.book1D("Ele_chIso_sum_gen_Sig_EE",
2274                      "Track pT sum in isolation cone around electron track after basic cuts - Signal Endcap",
2275                      nbin_2,
2276                      0,
2277                      20);
2278     meEleISO_rel_chIso_gen_Sig_EE_ =
2279         ibook.book1D("Ele_rel_chIso_sum_gen_Sig_EE",
2280                      "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Endcap",
2281                      nbin_1,
2282                      0,
2283                      4);
2284 
2285     meEleISO_Ntracks_MTD_3_Sig_EE_ =
2286         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Sig_EE",
2287                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2288                      20,
2289                      0,
2290                      20);
2291     meEleISO_chIso_MTD_3_Sig_EE_ = ibook.book1D(
2292         "Ele_chIso_sum_MTD_3_Sig_EE",
2293         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2294         nbin_2,
2295         0,
2296         20);
2297     meEleISO_rel_chIso_MTD_3_Sig_EE_ = ibook.book1D(
2298         "Ele_rel_chIso_sum_MTD_3_Sig_EE",
2299         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2300         nbin_1,
2301         0,
2302         4);
2303 
2304     meEleISO_Ntracks_MTD_4_Sig_EE_ =
2305         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Sig_EE",
2306                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2307                      20,
2308                      0,
2309                      20);
2310     meEleISO_chIso_MTD_4_Sig_EE_ = ibook.book1D(
2311         "Ele_chIso_sum_MTD_4_Sig_EE",
2312         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2313         nbin_2,
2314         0,
2315         20);
2316     meEleISO_rel_chIso_MTD_4_Sig_EE_ = ibook.book1D(
2317         "Ele_rel_chIso_sum_MTD_4_Sig_EE",
2318         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2319         nbin_1,
2320         0,
2321         4);
2322 
2323     meEleISO_Ntracks_MTD_5_Sig_EE_ =
2324         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Sig_EE",
2325                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2326                      20,
2327                      0,
2328                      20);
2329     meEleISO_chIso_MTD_5_Sig_EE_ = ibook.book1D(
2330         "Ele_chIso_sum_MTD_5_Sig_EE",
2331         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2332         nbin_2,
2333         0,
2334         20);
2335     meEleISO_rel_chIso_MTD_5_Sig_EE_ = ibook.book1D(
2336         "Ele_rel_chIso_sum_MTD_5_Sig_EE",
2337         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2338         nbin_1,
2339         0,
2340         4);
2341 
2342     meEleISO_Ntracks_MTD_6_Sig_EE_ =
2343         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Sig_EE",
2344                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2345                      20,
2346                      0,
2347                      20);
2348     meEleISO_chIso_MTD_6_Sig_EE_ = ibook.book1D(
2349         "Ele_chIso_sum_MTD_6_Sig_EE",
2350         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2351         nbin_2,
2352         0,
2353         20);
2354     meEleISO_rel_chIso_MTD_6_Sig_EE_ = ibook.book1D(
2355         "Ele_rel_chIso_sum_MTD_6_Sig_EE",
2356         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2357         nbin_1,
2358         0,
2359         4);
2360 
2361     meEleISO_Ntracks_MTD_7_Sig_EE_ =
2362         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Sig_EE",
2363                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2364                      20,
2365                      0,
2366                      20);
2367     meEleISO_chIso_MTD_7_Sig_EE_ = ibook.book1D(
2368         "Ele_chIso_sum_MTD_7_Sig_EE",
2369         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2370         nbin_2,
2371         0,
2372         20);
2373     meEleISO_rel_chIso_MTD_7_Sig_EE_ = ibook.book1D(
2374         "Ele_rel_chIso_sum_MTD_7_Sig_EE",
2375         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2376         nbin_1,
2377         0,
2378         4);
2379 
2380     meEleISO_Ntracks_MTD_sim_1_Sig_EE_ =
2381         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Sig_EE",
2382                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2383                      20,
2384                      0,
2385                      20);
2386     meEleISO_chIso_MTD_sim_1_Sig_EE_ = ibook.book1D(
2387         "Ele_chIso_sum_MTD_sim_1_Sig_EE",
2388         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2389         nbin_2,
2390         0,
2391         20);
2392     meEleISO_rel_chIso_MTD_sim_1_Sig_EE_ = ibook.book1D(
2393         "Ele_rel_chIso_sum_MTD_sim_1_Sig_EE",
2394         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2395         nbin_1,
2396         0,
2397         4);
2398 
2399     meEleISO_Ntracks_MTD_sim_2_Sig_EE_ =
2400         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Sig_EE",
2401                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2402                      20,
2403                      0,
2404                      20);
2405     meEleISO_chIso_MTD_sim_2_Sig_EE_ = ibook.book1D(
2406         "Ele_chIso_sum_MTD_sim_2_Sig_EE",
2407         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2408         nbin_2,
2409         0,
2410         20);
2411     meEleISO_rel_chIso_MTD_sim_2_Sig_EE_ = ibook.book1D(
2412         "Ele_rel_chIso_sum_MTD_sim_2_Sig_EE",
2413         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2414         nbin_1,
2415         0,
2416         4);
2417 
2418     meEleISO_Ntracks_MTD_sim_3_Sig_EE_ =
2419         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Sig_EE",
2420                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2421                      20,
2422                      0,
2423                      20);
2424     meEleISO_chIso_MTD_sim_3_Sig_EE_ = ibook.book1D(
2425         "Ele_chIso_sum_MTD_sim_3_Sig_EE",
2426         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2427         nbin_2,
2428         0,
2429         20);
2430     meEleISO_rel_chIso_MTD_sim_3_Sig_EE_ = ibook.book1D(
2431         "Ele_rel_chIso_sum_MTD_sim_3_Sig_EE",
2432         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2433         nbin_1,
2434         0,
2435         4);
2436 
2437     meEleISO_Ntracks_MTD_sim_4_Sig_EE_ =
2438         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Sig_EE",
2439                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2440                      20,
2441                      0,
2442                      20);
2443     meEleISO_chIso_MTD_sim_4_Sig_EE_ = ibook.book1D(
2444         "Ele_chIso_sum_MTD_sim_4_Sig_EE",
2445         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2446         nbin_2,
2447         0,
2448         20);
2449     meEleISO_rel_chIso_MTD_sim_4_Sig_EE_ = ibook.book1D(
2450         "Ele_rel_chIso_sum_MTD_sim_4_Sig_EE",
2451         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2452         nbin_1,
2453         0,
2454         4);
2455 
2456     meEleISO_Ntracks_MTD_sim_5_Sig_EE_ =
2457         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Sig_EE",
2458                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2459                      20,
2460                      0,
2461                      20);
2462     meEleISO_chIso_MTD_sim_5_Sig_EE_ = ibook.book1D(
2463         "Ele_chIso_sum_MTD_sim_5_Sig_EE",
2464         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2465         nbin_2,
2466         0,
2467         20);
2468     meEleISO_rel_chIso_MTD_sim_5_Sig_EE_ = ibook.book1D(
2469         "Ele_rel_chIso_sum_MTD_sim_5_Sig_EE",
2470         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2471         nbin_1,
2472         0,
2473         4);
2474 
2475     meEleISO_Ntracks_MTD_sim_6_Sig_EE_ =
2476         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Sig_EE",
2477                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2478                      20,
2479                      0,
2480                      20);
2481     meEleISO_chIso_MTD_sim_6_Sig_EE_ = ibook.book1D(
2482         "Ele_chIso_sum_MTD_sim_6_Sig_EE",
2483         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2484         nbin_2,
2485         0,
2486         20);
2487     meEleISO_rel_chIso_MTD_sim_6_Sig_EE_ = ibook.book1D(
2488         "Ele_rel_chIso_sum_MTD_sim_6_Sig_EE",
2489         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2490         nbin_1,
2491         0,
2492         4);
2493 
2494     meEleISO_Ntracks_MTD_sim_7_Sig_EE_ =
2495         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Sig_EE",
2496                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2497                      20,
2498                      0,
2499                      20);
2500     meEleISO_chIso_MTD_sim_7_Sig_EE_ = ibook.book1D(
2501         "Ele_chIso_sum_MTD_sim_7_Sig_EE",
2502         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2503         nbin_2,
2504         0,
2505         20);
2506     meEleISO_rel_chIso_MTD_sim_7_Sig_EE_ = ibook.book1D(
2507         "Ele_rel_chIso_sum_MTD_sim_7_Sig_EE",
2508         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2509         nbin_1,
2510         0,
2511         4);
2512   }
2513   meEleISO_Ntracks_MTD_4sigma_Sig_EE_ =
2514       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Sig_EE",
2515                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma significance - "
2516                    "Signal Endcap;Number of tracks;Counts",
2517                    20,
2518                    0,
2519                    20);
2520   meEleISO_chIso_MTD_4sigma_Sig_EE_ =
2521       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Sig_EE",
2522                    "Track pT sum in isolation cone around electron track after basic "
2523                    "cuts with MTD - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2524                    nbin_2,
2525                    0,
2526                    20);
2527   meEleISO_rel_chIso_MTD_4sigma_Sig_EE_ =
2528       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Sig_EE",
2529                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2530                    "significance - Signal Endcap;Isolation;Counts",
2531                    nbin_1,
2532                    0,
2533                    4);
2534 
2535   meEleISO_Ntracks_MTD_3sigma_Sig_EE_ =
2536       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Sig_EE",
2537                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma significance - "
2538                    "Signal Endcap;Number of tracks;Counts",
2539                    20,
2540                    0,
2541                    20);
2542   meEleISO_chIso_MTD_3sigma_Sig_EE_ =
2543       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Sig_EE",
2544                    "Track pT sum in isolation cone around electron track after basic "
2545                    "cuts with MTD - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2546                    nbin_2,
2547                    0,
2548                    20);
2549   meEleISO_rel_chIso_MTD_3sigma_Sig_EE_ =
2550       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Sig_EE",
2551                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 3 sigma "
2552                    "significance - Signal Endcap;Isolation;Counts",
2553                    nbin_1,
2554                    0,
2555                    4);
2556 
2557   meEleISO_Ntracks_MTD_2sigma_Sig_EE_ =
2558       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Sig_EE",
2559                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma significance - "
2560                    "Signal Endcap;Number of tracks;Counts",
2561                    20,
2562                    0,
2563                    20);
2564   meEleISO_chIso_MTD_2sigma_Sig_EE_ =
2565       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Sig_EE",
2566                    "Track pT sum in isolation cone around electron track after basic "
2567                    "cuts with MTD - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2568                    nbin_2,
2569                    0,
2570                    20);
2571   meEleISO_rel_chIso_MTD_2sigma_Sig_EE_ =
2572       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Sig_EE",
2573                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2574                    "significance - Signal Endcap;Isolation;Counts",
2575                    nbin_1,
2576                    0,
2577                    4);
2578 
2579   meEle_pt_tot_Sig_EE_ =
2580       ibook.book1D("Ele_pT_tot_Sig_EE", "Electron pT tot - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2581   meEle_pt_noMTD_Sig_EE_ =
2582       ibook.book1D("Ele_pT_noMTD_Sig_EE", "Electron pT noMTD - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2583 
2584   meEle_pt_sim_tot_Sig_EE_ =
2585       ibook.book1D("Ele_pT_sim_tot_Sig_EE", "Electron pT tot - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2586 
2587   meEle_eta_tot_Sig_EE_ =
2588       ibook.book1D("Ele_eta_tot_Sig_EE", "Electron eta tot - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2589   meEle_eta_noMTD_Sig_EE_ =
2590       ibook.book1D("Ele_eta_noMTD_Sig_EE", "Electron eta noMTD - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2591 
2592   meEle_phi_tot_Sig_EE_ =
2593       ibook.book1D("Ele_phi_tot_Sig_EE", "Electron phi tot - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2594   meEle_phi_noMTD_Sig_EE_ =
2595       ibook.book1D("Ele_phi_noMTD_Sig_EE", "Electron phi noMTD - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2596 
2597   if (optionalPlots_) {
2598     meEleISO_Ntracks_MTD_sim_4sigma_Sig_EE_ =
2599         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Sig_EE",
2600                      "Number of tracks in isolation cone around electron track after basic cuts with MTD SIM - 4 sigma "
2601                      "significance - Signal Endcap;Number of tracks;Counts",
2602                      20,
2603                      0,
2604                      20);
2605     meEleISO_chIso_MTD_sim_4sigma_Sig_EE_ =
2606         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Sig_EE",
2607                      "Track pT sum in isolation cone around electron track after "
2608                      "basic cuts with MTD SIM - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2609                      nbin_2,
2610                      0,
2611                      20);
2612     meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EE_ =
2613         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Sig_EE",
2614                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
2615                      "sigma significance - Signal Endcap;Isolation;Counts",
2616                      nbin_1,
2617                      0,
2618                      4);
2619 
2620     meEleISO_Ntracks_MTD_sim_3sigma_Sig_EE_ =
2621         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Sig_EE",
2622                      "Tracks in isolation cone around electron track after basic "
2623                      "cuts with MTD SIM - 3 sigma significance - Signal Endcap;Number of tracks;Counts",
2624                      20,
2625                      0,
2626                      20);
2627     meEleISO_chIso_MTD_sim_3sigma_Sig_EE_ =
2628         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Sig_EE",
2629                      "Track pT sum in isolation cone around electron track after "
2630                      "basic cuts with MTD SIM - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2631                      nbin_2,
2632                      0,
2633                      20);
2634     meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EE_ =
2635         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Sig_EE",
2636                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
2637                      "sigma significance - Signal Endcap;Isolation;Counts",
2638                      nbin_1,
2639                      0,
2640                      4);
2641 
2642     meEleISO_Ntracks_MTD_sim_2sigma_Sig_EE_ =
2643         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Sig_EE",
2644                      "Tracks in isolation cone around electron track after basic "
2645                      "cuts with MTD SIM - 2 sigma significance - Signal Endcap;Number of tracks;Counts",
2646                      20,
2647                      0,
2648                      20);
2649     meEleISO_chIso_MTD_sim_2sigma_Sig_EE_ =
2650         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Sig_EE",
2651                      "Track pT sum in isolation cone around electron track after "
2652                      "basic cuts with MTD SIM - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2653                      nbin_2,
2654                      0,
2655                      20);
2656     meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EE_ =
2657         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Sig_EE",
2658                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
2659                      "sigma significance - Signal Endcap;Isolation;Counts",
2660                      nbin_1,
2661                      0,
2662                      4);
2663 
2664     meEle_pt_MTD_1_Sig_EE_ = ibook.book1D("Ele_pT_MTD_1_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2665     meEle_eta_MTD_1_Sig_EE_ = ibook.book1D("Ele_eta_MTD_1_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2666     meEle_phi_MTD_1_Sig_EE_ = ibook.book1D("Ele_phi_MTD_1_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2667     meEle_pt_gen_Sig_EE_ =
2668         ibook.book1D("Ele_pT_gen_Sig_EE", "Electron pT genInfo - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2669     meEle_eta_gen_Sig_EE_ =
2670         ibook.book1D("Ele_eta_gen_Sig_EE", "Electron eta genInfo - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2671     meEle_phi_gen_Sig_EE_ =
2672         ibook.book1D("Ele_phi_gen_Sig_EE", "Electron phi genInfo - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2673 
2674     meEle_pt_MTD_2_Sig_EE_ = ibook.book1D("Ele_pT_MTD_2_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2675     meEle_eta_MTD_2_Sig_EE_ = ibook.book1D("Ele_eta_MTD_2_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2676     meEle_phi_MTD_2_Sig_EE_ = ibook.book1D("Ele_phi_MTD_2_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2677 
2678     meEle_pt_MTD_3_Sig_EE_ = ibook.book1D("Ele_pT_MTD_3_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2679     meEle_eta_MTD_3_Sig_EE_ = ibook.book1D("Ele_eta_MTD_3_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2680     meEle_phi_MTD_3_Sig_EE_ = ibook.book1D("Ele_phi_MTD_3_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2681 
2682     meEle_pt_MTD_4_Sig_EE_ = ibook.book1D("Ele_pT_MTD_4_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2683     meEle_eta_MTD_4_Sig_EE_ = ibook.book1D("Ele_eta_MTD_4_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2684     meEle_phi_MTD_4_Sig_EE_ = ibook.book1D("Ele_phi_MTD_4_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2685 
2686     meEle_pt_MTD_5_Sig_EE_ = ibook.book1D("Ele_pT_MTD_5_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2687     meEle_eta_MTD_5_Sig_EE_ = ibook.book1D("Ele_eta_MTD_5_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2688     meEle_phi_MTD_5_Sig_EE_ = ibook.book1D("Ele_phi_MTD_5_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2689 
2690     meEle_pt_MTD_6_Sig_EE_ = ibook.book1D("Ele_pT_MTD_6_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2691     meEle_eta_MTD_6_Sig_EE_ = ibook.book1D("Ele_eta_MTD_6_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2692     meEle_phi_MTD_6_Sig_EE_ = ibook.book1D("Ele_phi_MTD_6_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2693 
2694     meEle_pt_MTD_7_Sig_EE_ = ibook.book1D("Ele_pT_MTD_7_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2695     meEle_eta_MTD_7_Sig_EE_ = ibook.book1D("Ele_eta_MTD_7_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2696     meEle_phi_MTD_7_Sig_EE_ = ibook.book1D("Ele_phi_MTD_7_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2697 
2698     meEle_pt_sim_MTD_1_Sig_EE_ =
2699         ibook.book1D("Ele_pT_sim_MTD_1_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2700     meEle_pt_sim_MTD_2_Sig_EE_ =
2701         ibook.book1D("Ele_pT_sim_MTD_2_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2702     meEle_pt_sim_MTD_3_Sig_EE_ =
2703         ibook.book1D("Ele_pT_sim_MTD_3_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2704     meEle_pt_sim_MTD_4_Sig_EE_ =
2705         ibook.book1D("Ele_pT_sim_MTD_4_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2706     meEle_pt_sim_MTD_5_Sig_EE_ =
2707         ibook.book1D("Ele_pT_sim_MTD_5_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2708     meEle_pt_sim_MTD_6_Sig_EE_ =
2709         ibook.book1D("Ele_pT_sim_MTD_6_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2710     meEle_pt_sim_MTD_7_Sig_EE_ =
2711         ibook.book1D("Ele_pT_sim_MTD_7_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2712 
2713     meEle_pt_sim_MTD_4sigma_Sig_EE_ =
2714         ibook.book1D("Ele_pT_sim_MTD_4sigma_Sig_EE",
2715                      "Electron pT MTD SIM - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2716                      30,
2717                      10,
2718                      100);
2719     meEle_pt_sim_MTD_3sigma_Sig_EE_ =
2720         ibook.book1D("Ele_pT_sim_MTD_3sigma_Sig_EE",
2721                      "Electron pT MTD SIM - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2722                      30,
2723                      10,
2724                      100);
2725     meEle_pt_sim_MTD_2sigma_Sig_EE_ =
2726         ibook.book1D("Ele_pT_sim_MTD_2sigma_Sig_EE",
2727                      "Electron pT MTD SIM - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2728                      30,
2729                      10,
2730                      100);
2731   }
2732 
2733   meEle_pt_MTD_4sigma_Sig_EE_ =
2734       ibook.book1D("Ele_pT_MTD_4sigma_Sig_EE",
2735                    "Electron pT MTD - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2736                    30,
2737                    10,
2738                    100);
2739   meEle_eta_MTD_4sigma_Sig_EE_ = ibook.book1D(
2740       "Ele_eta_MTD_4sigma_Sig_EE", "Electron eta MTD - 4 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2741   meEle_phi_MTD_4sigma_Sig_EE_ = ibook.book1D(
2742       "Ele_phi_MTD_4sigma_Sig_EE", "Electron phi MTD - 4 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2743 
2744   meEle_pt_MTD_3sigma_Sig_EE_ =
2745       ibook.book1D("Ele_pT_MTD_3sigma_Sig_EE",
2746                    "Electron pT MTD - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2747                    30,
2748                    10,
2749                    100);
2750   meEle_eta_MTD_3sigma_Sig_EE_ = ibook.book1D(
2751       "Ele_eta_MTD_3sigma_Sig_EE", "Electron eta MTD - 3 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2752   meEle_phi_MTD_3sigma_Sig_EE_ = ibook.book1D(
2753       "Ele_phi_MTD_3sigma_Sig_EE", "Electron phi MTD - 3 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2754 
2755   meEle_pt_MTD_2sigma_Sig_EE_ =
2756       ibook.book1D("Ele_pT_MTD_2sigma_Sig_EE",
2757                    "Electron pT MTD - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2758                    30,
2759                    10,
2760                    100);
2761   meEle_eta_MTD_2sigma_Sig_EE_ = ibook.book1D(
2762       "Ele_eta_MTD_2sigma_Sig_EE", "Electron eta MTD - 2 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2763   meEle_phi_MTD_2sigma_Sig_EE_ = ibook.book1D(
2764       "Ele_phi_MTD_2sigma_Sig_EE", "Electron phi MTD - 2 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2765 
2766   // background
2767   meEleISO_Ntracks_Bkg_EB_ = ibook.book1D(
2768       "Ele_Iso_Ntracks_Bkg_EB",
2769       "Number of tracks in isolation cone around electron track after basic cuts - Bkg Barrel;Number of tracks;Counts",
2770       20,
2771       0,
2772       20);
2773   meEleISO_chIso_Bkg_EB_ = ibook.book1D(
2774       "Ele_chIso_sum_Bkg_EB",
2775       "Track pT sum in isolation cone around electron track after basic cuts - Bkg Barrel;p_{T} (GeV);Counts",
2776       nbin_2,
2777       0,
2778       20);
2779   meEleISO_rel_chIso_Bkg_EB_ = ibook.book1D(
2780       "Ele_rel_chIso_sum_Bkg_EB",
2781       "Track relative pT sum in isolation cone around electron track after basic cuts - Bkg Barrel;Isolation;Counts",
2782       nbin_1,
2783       0,
2784       4);
2785   if (optionalPlots_) {
2786     meEleISO_Ntracks_MTD_1_Bkg_EB_ =
2787         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Bkg_EB",
2788                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2789                      20,
2790                      0,
2791                      20);
2792     meEleISO_chIso_MTD_1_Bkg_EB_ = ibook.book1D(
2793         "Ele_chIso_sum_MTD_1_Bkg_EB",
2794         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2795         nbin_2,
2796         0,
2797         20);
2798     meEleISO_rel_chIso_MTD_1_Bkg_EB_ = ibook.book1D(
2799         "Ele_rel_chIso_sum_MTD_1_Bkg_EB",
2800         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2801         nbin_1,
2802         0,
2803         4);
2804 
2805     meEleISO_Ntracks_MTD_2_Bkg_EB_ =
2806         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Bkg_EB",
2807                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2808                      20,
2809                      0,
2810                      20);
2811     meEleISO_chIso_MTD_2_Bkg_EB_ = ibook.book1D(
2812         "Ele_chIso_sum_MTD_2_Bkg_EB",
2813         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2814         nbin_2,
2815         0,
2816         20);
2817     meEleISO_rel_chIso_MTD_2_Bkg_EB_ = ibook.book1D(
2818         "Ele_rel_chIso_sum_MTD_2_Bkg_EB",
2819         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2820         nbin_1,
2821         0,
2822         4);
2823     meEleISO_Ntracks_gen_Bkg_EB_ = ibook.book1D("Ele_Iso_Ntracks_gen_Bkg_EB",
2824                                                 "Tracks in isolation cone around electron track after basic cuts using "
2825                                                 "genInfo - Bkg Barrel;Number of tracks;Counts",
2826                                                 20,
2827                                                 0,
2828                                                 20);
2829     meEleISO_chIso_gen_Bkg_EB_ = ibook.book1D("Ele_chIso_sum_gen_Bkg_EB",
2830                                               "Track pT sum in isolation cone around electron track after basic cuts "
2831                                               "using genInfo - Bkg Barrel;p_{T} (GeV);Counts",
2832                                               nbin_2,
2833                                               0,
2834                                               20);
2835     meEleISO_rel_chIso_gen_Bkg_EB_ = ibook.book1D("Ele_rel_chIso_sum_gen_Bkg_EB",
2836                                                   "Track relative pT sum in isolation cone around electron track after "
2837                                                   "basic cuts using genInfo - Bkg Barrel;Isolation;Counts",
2838                                                   nbin_1,
2839                                                   0,
2840                                                   4);
2841     meEleISO_Ntracks_MTD_3_Bkg_EB_ =
2842         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Bkg_EB",
2843                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2844                      20,
2845                      0,
2846                      20);
2847     meEleISO_chIso_MTD_3_Bkg_EB_ = ibook.book1D(
2848         "Ele_chIso_sum_MTD_3_Bkg_EB",
2849         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2850         nbin_2,
2851         0,
2852         20);
2853     meEleISO_rel_chIso_MTD_3_Bkg_EB_ = ibook.book1D(
2854         "Ele_rel_chIso_sum_MTD_3_Bkg_EB",
2855         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2856         nbin_1,
2857         0,
2858         4);
2859 
2860     meEleISO_Ntracks_MTD_4_Bkg_EB_ =
2861         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Bkg_EB",
2862                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2863                      20,
2864                      0,
2865                      20);
2866     meEleISO_chIso_MTD_4_Bkg_EB_ = ibook.book1D(
2867         "Ele_chIso_sum_MTD_4_Bkg_EB",
2868         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2869         nbin_2,
2870         0,
2871         20);
2872     meEleISO_rel_chIso_MTD_4_Bkg_EB_ = ibook.book1D(
2873         "Ele_rel_chIso_sum_MTD_4_Bkg_EB",
2874         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2875         nbin_1,
2876         0,
2877         4);
2878 
2879     meEleISO_Ntracks_MTD_5_Bkg_EB_ =
2880         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Bkg_EB",
2881                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2882                      20,
2883                      0,
2884                      20);
2885     meEleISO_chIso_MTD_5_Bkg_EB_ = ibook.book1D(
2886         "Ele_chIso_sum_MTD_5_Bkg_EB",
2887         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2888         nbin_2,
2889         0,
2890         20);
2891     meEleISO_rel_chIso_MTD_5_Bkg_EB_ = ibook.book1D(
2892         "Ele_rel_chIso_sum_MTD_5_Bkg_EB",
2893         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2894         nbin_1,
2895         0,
2896         4);
2897 
2898     meEleISO_Ntracks_MTD_6_Bkg_EB_ =
2899         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Bkg_EB",
2900                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2901                      20,
2902                      0,
2903                      20);
2904     meEleISO_chIso_MTD_6_Bkg_EB_ = ibook.book1D(
2905         "Ele_chIso_sum_MTD_6_Bkg_EB",
2906         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2907         nbin_2,
2908         0,
2909         20);
2910     meEleISO_rel_chIso_MTD_6_Bkg_EB_ = ibook.book1D(
2911         "Ele_rel_chIso_sum_MTD_6_Bkg_EB",
2912         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2913         nbin_1,
2914         0,
2915         4);
2916 
2917     meEleISO_Ntracks_MTD_7_Bkg_EB_ =
2918         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Bkg_EB",
2919                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2920                      20,
2921                      0,
2922                      20);
2923     meEleISO_chIso_MTD_7_Bkg_EB_ = ibook.book1D(
2924         "Ele_chIso_sum_MTD_7_Bkg_EB",
2925         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2926         nbin_2,
2927         0,
2928         20);
2929     meEleISO_rel_chIso_MTD_7_Bkg_EB_ = ibook.book1D(
2930         "Ele_rel_chIso_sum_MTD_7_Bkg_EB",
2931         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2932         nbin_1,
2933         0,
2934         4);
2935 
2936     meEleISO_Ntracks_MTD_sim_1_Bkg_EB_ =
2937         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Bkg_EB",
2938                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2939                      20,
2940                      0,
2941                      20);
2942     meEleISO_chIso_MTD_sim_1_Bkg_EB_ = ibook.book1D(
2943         "Ele_chIso_sum_MTD_sim_1_Bkg_EB",
2944         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2945         nbin_2,
2946         0,
2947         20);
2948     meEleISO_rel_chIso_MTD_sim_1_Bkg_EB_ = ibook.book1D(
2949         "Ele_rel_chIso_sum_MTD_sim_1_Bkg_EB",
2950         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2951         nbin_1,
2952         0,
2953         4);
2954 
2955     meEleISO_Ntracks_MTD_sim_2_Bkg_EB_ =
2956         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Bkg_EB",
2957                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2958                      20,
2959                      0,
2960                      20);
2961     meEleISO_chIso_MTD_sim_2_Bkg_EB_ = ibook.book1D(
2962         "Ele_chIso_sum_MTD_sim_2_Bkg_EB",
2963         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2964         nbin_2,
2965         0,
2966         20);
2967     meEleISO_rel_chIso_MTD_sim_2_Bkg_EB_ = ibook.book1D(
2968         "Ele_rel_chIso_sum_MTD_sim_2_Bkg_EB",
2969         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2970         nbin_1,
2971         0,
2972         4);
2973 
2974     meEleISO_Ntracks_MTD_sim_3_Bkg_EB_ =
2975         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Bkg_EB",
2976                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2977                      20,
2978                      0,
2979                      20);
2980     meEleISO_chIso_MTD_sim_3_Bkg_EB_ = ibook.book1D(
2981         "Ele_chIso_sum_MTD_sim_3_Bkg_EB",
2982         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2983         nbin_2,
2984         0,
2985         20);
2986     meEleISO_rel_chIso_MTD_sim_3_Bkg_EB_ = ibook.book1D(
2987         "Ele_rel_chIso_sum_MTD_sim_3_Bkg_EB",
2988         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2989         nbin_1,
2990         0,
2991         4);
2992 
2993     meEleISO_Ntracks_MTD_sim_4_Bkg_EB_ =
2994         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Bkg_EB",
2995                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2996                      20,
2997                      0,
2998                      20);
2999     meEleISO_chIso_MTD_sim_4_Bkg_EB_ = ibook.book1D(
3000         "Ele_chIso_sum_MTD_sim_4_Bkg_EB",
3001         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3002         nbin_2,
3003         0,
3004         20);
3005     meEleISO_rel_chIso_MTD_sim_4_Bkg_EB_ = ibook.book1D(
3006         "Ele_rel_chIso_sum_MTD_sim_4_Bkg_EB",
3007         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3008         nbin_1,
3009         0,
3010         4);
3011 
3012     meEleISO_Ntracks_MTD_sim_5_Bkg_EB_ =
3013         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Bkg_EB",
3014                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3015                      20,
3016                      0,
3017                      20);
3018     meEleISO_chIso_MTD_sim_5_Bkg_EB_ = ibook.book1D(
3019         "Ele_chIso_sum_MTD_sim_5_Bkg_EB",
3020         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3021         nbin_2,
3022         0,
3023         20);
3024     meEleISO_rel_chIso_MTD_sim_5_Bkg_EB_ = ibook.book1D(
3025         "Ele_rel_chIso_sum_MTD_sim_5_Bkg_EB",
3026         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3027         nbin_1,
3028         0,
3029         4);
3030 
3031     meEleISO_Ntracks_MTD_sim_6_Bkg_EB_ =
3032         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Bkg_EB",
3033                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3034                      20,
3035                      0,
3036                      20);
3037     meEleISO_chIso_MTD_sim_6_Bkg_EB_ = ibook.book1D(
3038         "Ele_chIso_sum_MTD_sim_6_Bkg_EB",
3039         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3040         nbin_2,
3041         0,
3042         20);
3043     meEleISO_rel_chIso_MTD_sim_6_Bkg_EB_ = ibook.book1D(
3044         "Ele_rel_chIso_sum_MTD_sim_6_Bkg_EB",
3045         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3046         nbin_1,
3047         0,
3048         4);
3049 
3050     meEleISO_Ntracks_MTD_sim_7_Bkg_EB_ =
3051         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Bkg_EB",
3052                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3053                      20,
3054                      0,
3055                      20);
3056     meEleISO_chIso_MTD_sim_7_Bkg_EB_ = ibook.book1D(
3057         "Ele_chIso_sum_MTD_sim_7_Bkg_EB",
3058         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3059         nbin_2,
3060         0,
3061         20);
3062     meEleISO_rel_chIso_MTD_sim_7_Bkg_EB_ = ibook.book1D(
3063         "Ele_rel_chIso_sum_MTD_sim_7_Bkg_EB",
3064         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3065         nbin_1,
3066         0,
3067         4);
3068   }
3069   meEleISO_Ntracks_MTD_4sigma_Bkg_EB_ =
3070       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Bkg_EB",
3071                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma significance - "
3072                    "Bkg Barrel;Number of tracks;Counts",
3073                    20,
3074                    0,
3075                    20);
3076   meEleISO_chIso_MTD_4sigma_Bkg_EB_ =
3077       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Bkg_EB",
3078                    "Track pT sum in isolation cone around electron track after basic "
3079                    "cuts with MTD - 4 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3080                    nbin_2,
3081                    0,
3082                    20);
3083   meEleISO_rel_chIso_MTD_4sigma_Bkg_EB_ =
3084       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Bkg_EB",
3085                    "Track relative pT sum in isolation cone around electron track "
3086                    "after basic cuts with MTD - 4 sigma significance - Bkg Barrel;Isolation;Counts",
3087                    nbin_1,
3088                    0,
3089                    4);
3090 
3091   meEleISO_Ntracks_MTD_3sigma_Bkg_EB_ =
3092       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Bkg_EB",
3093                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma significance - "
3094                    "Bkg Barrel;Number of tracks;Counts",
3095                    20,
3096                    0,
3097                    20);
3098   meEleISO_chIso_MTD_3sigma_Bkg_EB_ =
3099       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Bkg_EB",
3100                    "Track pT sum in isolation cone around electron track after basic "
3101                    "cuts with MTD - 3 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3102                    nbin_2,
3103                    0,
3104                    20);
3105   meEleISO_rel_chIso_MTD_3sigma_Bkg_EB_ =
3106       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Bkg_EB",
3107                    "Track relative pT sum in isolation cone around electron track "
3108                    "after basic cuts with MTD - 3 sigma significance - Bkg Barrel;Isolation;Counts",
3109                    nbin_1,
3110                    0,
3111                    4);
3112 
3113   meEleISO_Ntracks_MTD_2sigma_Bkg_EB_ =
3114       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Bkg_EB",
3115                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma significance - "
3116                    "Bkg Barrel;Number of tracks;Counts",
3117                    20,
3118                    0,
3119                    20);
3120   meEleISO_chIso_MTD_2sigma_Bkg_EB_ =
3121       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Bkg_EB",
3122                    "Track pT sum in isolation cone around electron track after basic "
3123                    "cuts with MTD - 2 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3124                    nbin_2,
3125                    0,
3126                    20);
3127   meEleISO_rel_chIso_MTD_2sigma_Bkg_EB_ =
3128       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Bkg_EB",
3129                    "Track relative pT sum in isolation cone around electron track "
3130                    "after basic cuts with MTD - 2 sigma significance - Bkg Barrel;Isolation;Counts",
3131                    nbin_1,
3132                    0,
3133                    4);
3134 
3135   meEle_pt_tot_Bkg_EB_ =
3136       ibook.book1D("Ele_pT_tot_Bkg_EB", "Electron pT tot - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3137   meEle_pt_noMTD_Bkg_EB_ =
3138       ibook.book1D("Ele_pT_noMTD_Bkg_EB", "Electron pT noMTD - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3139 
3140   meEle_pt_sim_tot_Bkg_EB_ =
3141       ibook.book1D("Ele_pT_sim_tot_Bkg_EB", "Electron pT tot - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3142 
3143   meEle_eta_tot_Bkg_EB_ = ibook.book1D("Ele_eta_tot_Bkg_EB", "Electron eta tot - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3144   meEle_eta_noMTD_Bkg_EB_ =
3145       ibook.book1D("Ele_eta_noMTD_Bkg_EB", "Electron eta noMTD - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3146 
3147   meEle_phi_tot_Bkg_EB_ =
3148       ibook.book1D("Ele_phi_tot_Bkg_EB", "Electron phi tot - Bkg Barrel;#phi;#Counts", 64, -3.2, 3.2);
3149   meEle_phi_noMTD_Bkg_EB_ =
3150       ibook.book1D("Ele_phi_noMTD_Bkg_EB", "Electron phi noMTD - Bkg Barrel;#phi;#Counts", 64, -3.2, 3.2);
3151 
3152   if (optionalPlots_) {
3153     meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EB_ =
3154         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Bkg_EB",
3155                      "Tracks in isolation cone around electron track after basic "
3156                      "cuts with MTD SIM - 4 sigma significance - Bkg Barrel;Number of tracks;Counts",
3157                      20,
3158                      0,
3159                      20);
3160     meEleISO_chIso_MTD_sim_4sigma_Bkg_EB_ =
3161         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Bkg_EB",
3162                      "Track pT sum in isolation cone around electron track after "
3163                      "basic cuts with MTD SIM - 4 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3164                      nbin_2,
3165                      0,
3166                      20);
3167     meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EB_ =
3168         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Bkg_EB",
3169                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
3170                      "sigma significance - Bkg Barrel;Isolation;Counts",
3171                      nbin_1,
3172                      0,
3173                      4);
3174 
3175     meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EB_ =
3176         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Bkg_EB",
3177                      "Tracks in isolation cone around electron track after basic "
3178                      "cuts with MTD SIM - 3 sigma significance - Bkg Barrel;Number of tracks;Counts",
3179                      20,
3180                      0,
3181                      20);
3182     meEleISO_chIso_MTD_sim_3sigma_Bkg_EB_ =
3183         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Bkg_EB",
3184                      "Track pT sum in isolation cone around electron track after "
3185                      "basic cuts with MTD SIM - 3 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3186                      nbin_2,
3187                      0,
3188                      20);
3189     meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EB_ =
3190         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Bkg_EB",
3191                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
3192                      "sigma significance - Bkg Barrel;Isolation;Counts",
3193                      nbin_1,
3194                      0,
3195                      4);
3196 
3197     meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EB_ =
3198         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Bkg_EB",
3199                      "Tracks in isolation cone around electron track after basic "
3200                      "cuts with MTD SIM - 3 sigma significance - Bkg Barrel;Number of tracks;Counts",
3201                      20,
3202                      0,
3203                      20);
3204     meEleISO_chIso_MTD_sim_2sigma_Bkg_EB_ =
3205         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Bkg_EB",
3206                      "Track pT sum in isolation cone around electron track after "
3207                      "basic cuts with MTD SIM - 2 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3208                      nbin_2,
3209                      0,
3210                      20);
3211     meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EB_ =
3212         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Bkg_EB",
3213                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
3214                      "sigma significance - Bkg Barrel;Isolation;Counts",
3215                      nbin_1,
3216                      0,
3217                      4);
3218 
3219     meEle_pt_MTD_1_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_1_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3220     meEle_eta_MTD_1_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_1_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3221     meEle_phi_MTD_1_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_1_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3222     meEle_pt_gen_Bkg_EB_ =
3223         ibook.book1D("Ele_pT_gen_Bkg_EB", "Electron pT genInfo - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3224     meEle_eta_gen_Bkg_EB_ =
3225         ibook.book1D("Ele_eta_gen_Bkg_EB", "Electron eta genInfo - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3226     meEle_phi_gen_Bkg_EB_ =
3227         ibook.book1D("Ele_phi_gen_Bkg_EB", "Electron phi genInfo - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3228 
3229     meEle_pt_MTD_2_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_2_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3230     meEle_eta_MTD_2_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_2_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3231     meEle_phi_MTD_2_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_2_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3232 
3233     meEle_pt_MTD_3_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_3_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3234     meEle_eta_MTD_3_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_3_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3235     meEle_phi_MTD_3_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_3_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3236 
3237     meEle_pt_MTD_4_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_4_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3238     meEle_eta_MTD_4_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_4_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3239     meEle_phi_MTD_4_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_4_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3240 
3241     meEle_pt_MTD_5_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_5_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3242     meEle_eta_MTD_5_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_5_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3243     meEle_phi_MTD_5_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_5_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3244 
3245     meEle_pt_MTD_6_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_6_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3246     meEle_eta_MTD_6_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_6_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3247     meEle_phi_MTD_6_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_6_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3248 
3249     meEle_pt_MTD_7_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_7_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3250     meEle_eta_MTD_7_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_7_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3251     meEle_phi_MTD_7_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_7_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3252 
3253     meEle_pt_sim_MTD_1_Bkg_EB_ =
3254         ibook.book1D("Ele_pT_sim_MTD_1_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3255     meEle_pt_sim_MTD_2_Bkg_EB_ =
3256         ibook.book1D("Ele_pT_sim_MTD_2_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3257     meEle_pt_sim_MTD_3_Bkg_EB_ =
3258         ibook.book1D("Ele_pT_sim_MTD_3_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3259     meEle_pt_sim_MTD_4_Bkg_EB_ =
3260         ibook.book1D("Ele_pT_sim_MTD_4_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3261     meEle_pt_sim_MTD_5_Bkg_EB_ =
3262         ibook.book1D("Ele_pT_sim_MTD_5_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3263     meEle_pt_sim_MTD_6_Bkg_EB_ =
3264         ibook.book1D("Ele_pT_sim_MTD_6_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3265     meEle_pt_sim_MTD_7_Bkg_EB_ =
3266         ibook.book1D("Ele_pT_sim_MTD_7_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3267   }
3268   meEle_pt_MTD_4sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_4sigma_Bkg_EB",
3269                                              "Electron pT MTD - 4 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3270                                              30,
3271                                              10,
3272                                              100);
3273   meEle_eta_MTD_4sigma_Bkg_EB_ = ibook.book1D(
3274       "Ele_eta_MTD_4sigma_Bkg_EB", "Electron eta MTD - 4 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3275   meEle_phi_MTD_4sigma_Bkg_EB_ = ibook.book1D(
3276       "Ele_phi_MTD_4sigma_Bkg_EB", "Electron phi MTD - 4 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3277 
3278   meEle_pt_MTD_3sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_3sigma_Bkg_EB",
3279                                              "Electron pT MTD - 3 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3280                                              30,
3281                                              10,
3282                                              100);
3283   meEle_eta_MTD_3sigma_Bkg_EB_ = ibook.book1D(
3284       "Ele_eta_MTD_3sigma_Bkg_EB", "Electron eta MTD - 3 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3285   meEle_phi_MTD_3sigma_Bkg_EB_ = ibook.book1D(
3286       "Ele_phi_MTD_3sigma_Bkg_EB", "Electron phi MTD - 3 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3287 
3288   meEle_pt_MTD_2sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_2sigma_Bkg_EB",
3289                                              "Electron pT MTD - 2 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3290                                              30,
3291                                              10,
3292                                              100);
3293   meEle_eta_MTD_2sigma_Bkg_EB_ = ibook.book1D(
3294       "Ele_eta_MTD_2sigma_Bkg_EB", "Electron eta MTD - 2 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3295   meEle_phi_MTD_2sigma_Bkg_EB_ = ibook.book1D(
3296       "Ele_phi_MTD_2sigma_Bkg_EB", "Electron phi MTD - 2 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3297 
3298   meEleISO_Ntracks_Bkg_EE_ = ibook.book1D(
3299       "Ele_Iso_Ntracks_Bkg_EE",
3300       "Number of tracks in isolation cone around electron track after basic cuts - Bkg Endcap;Number of tracks;Counts",
3301       20,
3302       0,
3303       20);
3304   meEleISO_chIso_Bkg_EE_ = ibook.book1D(
3305       "Ele_chIso_sum_Bkg_EE",
3306       "Track pT sum in isolation cone around electron track after basic cuts - Bkg Endcap;p_{T} (GeV);Counts",
3307       nbin_2,
3308       0,
3309       20);
3310   meEleISO_rel_chIso_Bkg_EE_ = ibook.book1D(
3311       "Ele_rel_chIso_sum_Bkg_EE",
3312       "Track relative pT sum in isolation cone around electron track after basic cuts - Bkg Endcap;Isolation;Counts",
3313       nbin_1,
3314       0,
3315       4);
3316   if (optionalPlots_) {
3317     meEle_pt_sim_MTD_4sigma_Bkg_EB_ =
3318         ibook.book1D("Ele_pT_sim_MTD_4sigma_Bkg_EB",
3319                      "Electron pT MTD SIM - 4 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3320                      30,
3321                      10,
3322                      100);
3323     meEle_pt_sim_MTD_3sigma_Bkg_EB_ =
3324         ibook.book1D("Ele_pT_sim_MTD_3sigma_Bkg_EB",
3325                      "Electron pT MTD SIM - 3 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3326                      30,
3327                      10,
3328                      100);
3329     meEle_pt_sim_MTD_2sigma_Bkg_EB_ =
3330         ibook.book1D("Ele_pT_sim_MTD_2sigma_Bkg_EB",
3331                      "Electron pT MTD SIM - 2 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3332                      30,
3333                      10,
3334                      100);
3335 
3336     meEleISO_Ntracks_MTD_1_Bkg_EE_ =
3337         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Bkg_EE",
3338                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3339                      20,
3340                      0,
3341                      20);
3342     meEleISO_chIso_MTD_1_Bkg_EE_ = ibook.book1D(
3343         "Ele_chIso_sum_MTD_1_Bkg_EE",
3344         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3345         nbin_2,
3346         0,
3347         20);
3348     meEleISO_rel_chIso_MTD_1_Bkg_EE_ = ibook.book1D(
3349         "Ele_rel_chIso_sum_MTD_1_Bkg_EE",
3350         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3351         nbin_1,
3352         0,
3353         4);
3354 
3355     meEleISO_Ntracks_MTD_2_Bkg_EE_ =
3356         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Bkg_EE",
3357                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3358                      20,
3359                      0,
3360                      20);
3361     meEleISO_chIso_MTD_2_Bkg_EE_ = ibook.book1D(
3362         "Ele_chIso_sum_MTD_2_Bkg_EE",
3363         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3364         nbin_2,
3365         0,
3366         20);
3367     meEleISO_rel_chIso_MTD_2_Bkg_EE_ = ibook.book1D(
3368         "Ele_rel_chIso_sum_MTD_2_Bkg_EE",
3369         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3370         nbin_1,
3371         0,
3372         4);
3373     meEleISO_Ntracks_gen_Bkg_EE_ = ibook.book1D("Ele_Iso_Ntracks_gen_Bkg_EE",
3374                                                 "Tracks in isolation cone around electron track after basic cuts using "
3375                                                 "genInfo - Bkg Endcap;Number of tracks;Counts",
3376                                                 20,
3377                                                 0,
3378                                                 20);
3379     meEleISO_chIso_gen_Bkg_EE_ = ibook.book1D("Ele_chIso_sum_gen_Bkg_EE",
3380                                               "Track pT sum in isolation cone around electron track after basic cuts "
3381                                               "using genInfo - Bkg Endcap;p_{T} (GeV);Counts",
3382                                               nbin_2,
3383                                               0,
3384                                               20);
3385     meEleISO_rel_chIso_gen_Bkg_EE_ = ibook.book1D("Ele_rel_chIso_sum_gen_Bkg_EE",
3386                                                   "Track relative pT sum in isolation cone around electron track after "
3387                                                   "basic cuts using genInfo - Bkg Endcap;Isolation;Counts",
3388                                                   nbin_1,
3389                                                   0,
3390                                                   4);
3391 
3392     meEleISO_Ntracks_MTD_3_Bkg_EE_ =
3393         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Bkg_EE",
3394                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3395                      20,
3396                      0,
3397                      20);
3398     meEleISO_chIso_MTD_3_Bkg_EE_ = ibook.book1D(
3399         "Ele_chIso_sum_MTD_3_Bkg_EE",
3400         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3401         nbin_2,
3402         0,
3403         20);
3404     meEleISO_rel_chIso_MTD_3_Bkg_EE_ = ibook.book1D(
3405         "Ele_rel_chIso_sum_MTD_3_Bkg_EE",
3406         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3407         nbin_1,
3408         0,
3409         4);
3410 
3411     meEleISO_Ntracks_MTD_4_Bkg_EE_ =
3412         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Bkg_EE",
3413                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3414                      20,
3415                      0,
3416                      20);
3417     meEleISO_chIso_MTD_4_Bkg_EE_ = ibook.book1D(
3418         "Ele_chIso_sum_MTD_4_Bkg_EE",
3419         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3420         nbin_2,
3421         0,
3422         20);
3423     meEleISO_rel_chIso_MTD_4_Bkg_EE_ = ibook.book1D(
3424         "Ele_rel_chIso_sum_MTD_4_Bkg_EE",
3425         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3426         nbin_1,
3427         0,
3428         4);
3429 
3430     meEleISO_Ntracks_MTD_5_Bkg_EE_ =
3431         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Bkg_EE",
3432                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3433                      20,
3434                      0,
3435                      20);
3436     meEleISO_chIso_MTD_5_Bkg_EE_ = ibook.book1D(
3437         "Ele_chIso_sum_MTD_5_Bkg_EE",
3438         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3439         nbin_2,
3440         0,
3441         20);
3442     meEleISO_rel_chIso_MTD_5_Bkg_EE_ = ibook.book1D(
3443         "Ele_rel_chIso_sum_MTD_5_Bkg_EE",
3444         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3445         nbin_1,
3446         0,
3447         4);
3448 
3449     meEleISO_Ntracks_MTD_6_Bkg_EE_ =
3450         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Bkg_EE",
3451                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3452                      20,
3453                      0,
3454                      20);
3455     meEleISO_chIso_MTD_6_Bkg_EE_ = ibook.book1D(
3456         "Ele_chIso_sum_MTD_6_Bkg_EE",
3457         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3458         nbin_2,
3459         0,
3460         20);
3461     meEleISO_rel_chIso_MTD_6_Bkg_EE_ = ibook.book1D(
3462         "Ele_rel_chIso_sum_MTD_6_Bkg_EE",
3463         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3464         nbin_1,
3465         0,
3466         4);
3467 
3468     meEleISO_Ntracks_MTD_7_Bkg_EE_ =
3469         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Bkg_EE",
3470                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3471                      20,
3472                      0,
3473                      20);
3474     meEleISO_chIso_MTD_7_Bkg_EE_ = ibook.book1D(
3475         "Ele_chIso_sum_MTD_7_Bkg_EE",
3476         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3477         nbin_2,
3478         0,
3479         20);
3480     meEleISO_rel_chIso_MTD_7_Bkg_EE_ = ibook.book1D(
3481         "Ele_rel_chIso_sum_MTD_7_Bkg_EE",
3482         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3483         nbin_1,
3484         0,
3485         4);
3486 
3487     meEleISO_Ntracks_MTD_sim_1_Bkg_EE_ =
3488         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Bkg_EE",
3489                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3490                      20,
3491                      0,
3492                      20);
3493     meEleISO_chIso_MTD_sim_1_Bkg_EE_ = ibook.book1D(
3494         "Ele_chIso_sum_MTD_sim_1_Bkg_EE",
3495         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3496         nbin_2,
3497         0,
3498         20);
3499     meEleISO_rel_chIso_MTD_sim_1_Bkg_EE_ = ibook.book1D(
3500         "Ele_rel_chIso_sum_MTD_sim_1_Bkg_EE",
3501         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3502         nbin_1,
3503         0,
3504         4);
3505 
3506     meEleISO_Ntracks_MTD_sim_2_Bkg_EE_ =
3507         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Bkg_EE",
3508                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3509                      20,
3510                      0,
3511                      20);
3512     meEleISO_chIso_MTD_sim_2_Bkg_EE_ = ibook.book1D(
3513         "Ele_chIso_sum_MTD_sim_2_Bkg_EE",
3514         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3515         nbin_2,
3516         0,
3517         20);
3518     meEleISO_rel_chIso_MTD_sim_2_Bkg_EE_ = ibook.book1D(
3519         "Ele_rel_chIso_sum_MTD_sim_2_Bkg_EE",
3520         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3521         nbin_1,
3522         0,
3523         4);
3524 
3525     meEleISO_Ntracks_MTD_sim_3_Bkg_EE_ =
3526         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Bkg_EE",
3527                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3528                      20,
3529                      0,
3530                      20);
3531     meEleISO_chIso_MTD_sim_3_Bkg_EE_ = ibook.book1D(
3532         "Ele_chIso_sum_MTD_sim_3_Bkg_EE",
3533         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3534         nbin_2,
3535         0,
3536         20);
3537     meEleISO_rel_chIso_MTD_sim_3_Bkg_EE_ = ibook.book1D(
3538         "Ele_rel_chIso_sum_MTD_sim_3_Bkg_EE",
3539         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3540         nbin_1,
3541         0,
3542         4);
3543 
3544     meEleISO_Ntracks_MTD_sim_4_Bkg_EE_ =
3545         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Bkg_EE",
3546                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3547                      20,
3548                      0,
3549                      20);
3550     meEleISO_chIso_MTD_sim_4_Bkg_EE_ = ibook.book1D(
3551         "Ele_chIso_sum_MTD_sim_4_Bkg_EE",
3552         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3553         nbin_2,
3554         0,
3555         20);
3556     meEleISO_rel_chIso_MTD_sim_4_Bkg_EE_ = ibook.book1D(
3557         "Ele_rel_chIso_sum_MTD_sim_4_Bkg_EE",
3558         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3559         nbin_1,
3560         0,
3561         4);
3562 
3563     meEleISO_Ntracks_MTD_sim_5_Bkg_EE_ =
3564         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Bkg_EE",
3565                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3566                      20,
3567                      0,
3568                      20);
3569     meEleISO_chIso_MTD_sim_5_Bkg_EE_ = ibook.book1D(
3570         "Ele_chIso_sum_MTD_sim_5_Bkg_EE",
3571         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3572         nbin_2,
3573         0,
3574         20);
3575     meEleISO_rel_chIso_MTD_sim_5_Bkg_EE_ = ibook.book1D(
3576         "Ele_rel_chIso_sum_MTD_sim_5_Bkg_EE",
3577         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3578         nbin_1,
3579         0,
3580         4);
3581 
3582     meEleISO_Ntracks_MTD_sim_6_Bkg_EE_ =
3583         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Bkg_EE",
3584                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3585                      20,
3586                      0,
3587                      20);
3588     meEleISO_chIso_MTD_sim_6_Bkg_EE_ = ibook.book1D(
3589         "Ele_chIso_sum_MTD_sim_6_Bkg_EE",
3590         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3591         nbin_2,
3592         0,
3593         20);
3594     meEleISO_rel_chIso_MTD_sim_6_Bkg_EE_ = ibook.book1D(
3595         "Ele_rel_chIso_sum_MTD_sim_6_Bkg_EE",
3596         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3597         nbin_1,
3598         0,
3599         4);
3600 
3601     meEleISO_Ntracks_MTD_sim_7_Bkg_EE_ =
3602         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Bkg_EE",
3603                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3604                      20,
3605                      0,
3606                      20);
3607     meEleISO_chIso_MTD_sim_7_Bkg_EE_ = ibook.book1D(
3608         "Ele_chIso_sum_MTD_sim_7_Bkg_EE",
3609         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3610         nbin_2,
3611         0,
3612         20);
3613     meEleISO_rel_chIso_MTD_sim_7_Bkg_EE_ = ibook.book1D(
3614         "Ele_rel_chIso_sum_MTD_sim_7_Bkg_EE",
3615         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3616         nbin_1,
3617         0,
3618         4);
3619   }
3620   meEleISO_Ntracks_MTD_4sigma_Bkg_EE_ =
3621       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Bkg_EE",
3622                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma compatibility - "
3623                    "Bkg Endcap;Number of tracks;Counts",
3624                    20,
3625                    0,
3626                    20);
3627   meEleISO_chIso_MTD_4sigma_Bkg_EE_ =
3628       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Bkg_EE",
3629                    "Track pT sum in isolation cone around electron track after basic "
3630                    "cuts with MTD - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3631                    nbin_2,
3632                    0,
3633                    20);
3634   meEleISO_rel_chIso_MTD_4sigma_Bkg_EE_ =
3635       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Bkg_EE",
3636                    "Track relative pT sum in isolation cone around electron track "
3637                    "after basic cuts with MTD - 4 sigma compatibility - Bkg Endcap;Isolation;Counts",
3638                    nbin_1,
3639                    0,
3640                    4);
3641 
3642   meEleISO_Ntracks_MTD_3sigma_Bkg_EE_ =
3643       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Bkg_EE",
3644                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma compatibility - "
3645                    "Bkg Endcap;Number of tracks;Counts",
3646                    20,
3647                    0,
3648                    20);
3649   meEleISO_chIso_MTD_3sigma_Bkg_EE_ =
3650       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Bkg_EE",
3651                    "Track pT sum in isolation cone around electron track after basic "
3652                    "cuts with MTD - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3653                    nbin_2,
3654                    0,
3655                    20);
3656   meEleISO_rel_chIso_MTD_3sigma_Bkg_EE_ =
3657       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Bkg_EE",
3658                    "Track relative pT sum in isolation cone around electron track "
3659                    "after basic cuts with MTD - 3 sigma compatibility - Bkg Endcap;Isolation;Counts",
3660                    nbin_1,
3661                    0,
3662                    4);
3663 
3664   meEleISO_Ntracks_MTD_2sigma_Bkg_EE_ =
3665       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Bkg_EE",
3666                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibility - "
3667                    "Bkg Endcap;Number of tracks;Counts",
3668                    20,
3669                    0,
3670                    20);
3671   meEleISO_chIso_MTD_2sigma_Bkg_EE_ =
3672       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Bkg_EE",
3673                    "Track pT sum in isolation cone around electron track after basic "
3674                    "cuts with MTD - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3675                    nbin_2,
3676                    0,
3677                    20);
3678   meEleISO_rel_chIso_MTD_2sigma_Bkg_EE_ =
3679       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Bkg_EE",
3680                    "Track relative pT sum in isolation cone around electron track "
3681                    "after basic cuts with MTD - 2 sigma compatibility - Bkg Endcap;Isolation;Counts",
3682                    nbin_1,
3683                    0,
3684                    4);
3685 
3686   meEle_pt_tot_Bkg_EE_ =
3687       ibook.book1D("Ele_pT_tot_Bkg_EE", "Electron pT tot - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3688   meEle_pt_noMTD_Bkg_EE_ =
3689       ibook.book1D("Ele_pT_noMTD_Bkg_EE", "Electron pT noMTD - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3690 
3691   meEle_pt_sim_tot_Bkg_EE_ =
3692       ibook.book1D("Ele_pT_sim_tot_Bkg_EE", "Electron pT tot - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3693 
3694   meEle_eta_tot_Bkg_EE_ = ibook.book1D("Ele_eta_tot_Bkg_EE", "Electron eta tot - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3695   meEle_eta_noMTD_Bkg_EE_ =
3696       ibook.book1D("Ele_eta_noMTD_Bkg_EE", "Electron eta noMTD - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3697 
3698   meEle_phi_tot_Bkg_EE_ =
3699       ibook.book1D("Ele_phi_tot_Bkg_EE", "Electron phi tot - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3700   meEle_phi_noMTD_Bkg_EE_ =
3701       ibook.book1D("Ele_phi_noMTD_Bkg_EE", "Electron phi noMTD - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3702   if (optionalPlots_) {
3703     meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EE_ =
3704         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Bkg_EE",
3705                      "Tracks in isolation cone around electron track after basic "
3706                      "cuts with MTD SIM - 4 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3707                      20,
3708                      0,
3709                      20);
3710     meEleISO_chIso_MTD_sim_4sigma_Bkg_EE_ =
3711         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Bkg_EE",
3712                      "Track pT sum in isolation cone around electron track after "
3713                      "basic cuts with MTD SIM - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3714                      nbin_2,
3715                      0,
3716                      20);
3717     meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EE_ =
3718         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Bkg_EE",
3719                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
3720                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3721                      nbin_1,
3722                      0,
3723                      4);
3724 
3725     meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EE_ =
3726         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Bkg_EE",
3727                      "Tracks in isolation cone around electron track after basic "
3728                      "cuts with MTD SIM - 3 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3729                      20,
3730                      0,
3731                      20);
3732     meEleISO_chIso_MTD_sim_3sigma_Bkg_EE_ =
3733         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Bkg_EE",
3734                      "Track pT sum in isolation cone around electron track after "
3735                      "basic cuts with MTD SIM - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3736                      nbin_2,
3737                      0,
3738                      20);
3739     meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EE_ =
3740         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Bkg_EE",
3741                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
3742                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3743                      nbin_1,
3744                      0,
3745                      4);
3746 
3747     meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EE_ =
3748         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Bkg_EE",
3749                      "Tracks in isolation cone around electron track after basic "
3750                      "cuts with MTD SIM - 2 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3751                      20,
3752                      0,
3753                      20);
3754     meEleISO_chIso_MTD_sim_2sigma_Bkg_EE_ =
3755         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Bkg_EE",
3756                      "Track pT sum in isolation cone around electron track after "
3757                      "basic cuts with MTD SIM - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3758                      nbin_2,
3759                      0,
3760                      20);
3761     meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EE_ =
3762         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Bkg_EE",
3763                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
3764                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3765                      nbin_1,
3766                      0,
3767                      4);
3768 
3769     meEle_pt_MTD_1_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_1_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3770     meEle_eta_MTD_1_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_1_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3771     meEle_phi_MTD_1_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_1_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3772     meEle_pt_gen_Bkg_EE_ =
3773         ibook.book1D("Ele_pT_gen_Bkg_EE", "Electron pT genInfo - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3774     meEle_eta_gen_Bkg_EE_ =
3775         ibook.book1D("Ele_eta_gen_Bkg_EE", "Electron eta genInfo - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3776     meEle_phi_gen_Bkg_EE_ =
3777         ibook.book1D("Ele_phi_gen_Bkg_EE", "Electron phi genInfo - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3778 
3779     meEle_pt_MTD_2_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_2_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3780     meEle_eta_MTD_2_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_2_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3781     meEle_phi_MTD_2_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_2_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3782 
3783     meEle_pt_MTD_3_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_3_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3784     meEle_eta_MTD_3_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_3_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3785     meEle_phi_MTD_3_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_3_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3786 
3787     meEle_pt_MTD_4_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_4_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3788     meEle_eta_MTD_4_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_4_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3789     meEle_phi_MTD_4_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_4_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3790 
3791     meEle_pt_MTD_5_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_5_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3792     meEle_eta_MTD_5_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_5_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3793     meEle_phi_MTD_5_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_5_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3794 
3795     meEle_pt_MTD_6_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_6_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3796     meEle_eta_MTD_6_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_6_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3797     meEle_phi_MTD_6_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_6_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3798 
3799     meEle_pt_MTD_7_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_7_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3800     meEle_eta_MTD_7_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_7_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3801     meEle_phi_MTD_7_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_7_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3802 
3803     meEle_pt_sim_MTD_1_Bkg_EE_ =
3804         ibook.book1D("Ele_pT_sim_MTD_1_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3805     meEle_pt_sim_MTD_2_Bkg_EE_ =
3806         ibook.book1D("Ele_pT_sim_MTD_2_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3807     meEle_pt_sim_MTD_3_Bkg_EE_ =
3808         ibook.book1D("Ele_pT_sim_MTD_3_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3809     meEle_pt_sim_MTD_4_Bkg_EE_ =
3810         ibook.book1D("Ele_pT_sim_MTD_4_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3811     meEle_pt_sim_MTD_5_Bkg_EE_ =
3812         ibook.book1D("Ele_pT_sim_MTD_5_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3813     meEle_pt_sim_MTD_6_Bkg_EE_ =
3814         ibook.book1D("Ele_pT_sim_MTD_6_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3815     meEle_pt_sim_MTD_7_Bkg_EE_ =
3816         ibook.book1D("Ele_pT_sim_MTD_7_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3817   }
3818 
3819   meEle_pt_MTD_4sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_4sigma_Bkg_EE",
3820                                              "Electron pT MTD - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3821                                              30,
3822                                              10,
3823                                              100);
3824   meEle_eta_MTD_4sigma_Bkg_EE_ = ibook.book1D(
3825       "Ele_eta_MTD_4sigma_Bkg_EE", "Electron eta MTD - 4 sigma compatibility - Bkg Endcapi;#eta;Counts", 32, 1.6, 3.2);
3826   meEle_phi_MTD_4sigma_Bkg_EE_ = ibook.book1D(
3827       "Ele_phi_MTD_4sigma_Bkg_EE", "Electron phi MTD - 4 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3828 
3829   meEle_pt_MTD_3sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_3sigma_Bkg_EE",
3830                                              "Electron pT MTD - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3831                                              30,
3832                                              10,
3833                                              100);
3834   meEle_eta_MTD_3sigma_Bkg_EE_ = ibook.book1D(
3835       "Ele_eta_MTD_3sigma_Bkg_EE", "Electron eta MTD - 3 sigma compatibility - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3836   meEle_phi_MTD_3sigma_Bkg_EE_ = ibook.book1D(
3837       "Ele_phi_MTD_3sigma_Bkg_EE", "Electron phi MTD - 3 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3838 
3839   meEle_pt_MTD_2sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_2sigma_Bkg_EE",
3840                                              "Electron pT MTD - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3841                                              30,
3842                                              10,
3843                                              100);
3844   meEle_eta_MTD_2sigma_Bkg_EE_ = ibook.book1D(
3845       "Ele_eta_MTD_2sigma_Bkg_EE", "Electron eta MTD - 2 sigma compatibility - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3846   meEle_phi_MTD_2sigma_Bkg_EE_ = ibook.book1D(
3847       "Ele_phi_MTD_2sigma_Bkg_EE", "Electron phi MTD - 2 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3848 
3849   if (optionalPlots_) {
3850     meEle_pt_sim_MTD_4sigma_Bkg_EE_ =
3851         ibook.book1D("Ele_pT_sim_MTD_4sigma_Bkg_EE",
3852                      "Electron pT MTD SIM - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3853                      30,
3854                      10,
3855                      100);
3856     meEle_pt_sim_MTD_3sigma_Bkg_EE_ =
3857         ibook.book1D("Ele_pT_sim_MTD_3sigma_Bkg_EE",
3858                      "Electron pT MTD SIM - 3 sigma compatibility - Bkg Endcap;#eta;Counts",
3859                      30,
3860                      10,
3861                      100);
3862     meEle_pt_sim_MTD_2sigma_Bkg_EE_ =
3863         ibook.book1D("Ele_pT_sim_MTD_2sigma_Bkg_EE",
3864                      "Electron pT MTD SIM - 2 sigma compatibility - Bkg Endcap;#phi;Counts",
3865                      30,
3866                      10,
3867                      100);
3868   }
3869 
3870   // defining vectors for more efficient hist filling
3871   // Promt part
3872   if (optionalPlots_) {
3873     Ntracks_EB_list_Sig = {meEleISO_Ntracks_MTD_1_Sig_EB_,
3874                            meEleISO_Ntracks_MTD_2_Sig_EB_,
3875                            meEleISO_Ntracks_MTD_3_Sig_EB_,
3876                            meEleISO_Ntracks_MTD_4_Sig_EB_,
3877                            meEleISO_Ntracks_MTD_5_Sig_EB_,
3878                            meEleISO_Ntracks_MTD_6_Sig_EB_,
3879                            meEleISO_Ntracks_MTD_7_Sig_EB_};
3880     ch_iso_EB_list_Sig = {meEleISO_chIso_MTD_1_Sig_EB_,
3881                           meEleISO_chIso_MTD_2_Sig_EB_,
3882                           meEleISO_chIso_MTD_3_Sig_EB_,
3883                           meEleISO_chIso_MTD_4_Sig_EB_,
3884                           meEleISO_chIso_MTD_5_Sig_EB_,
3885                           meEleISO_chIso_MTD_6_Sig_EB_,
3886                           meEleISO_chIso_MTD_7_Sig_EB_};
3887     rel_ch_iso_EB_list_Sig = {meEleISO_rel_chIso_MTD_1_Sig_EB_,
3888                               meEleISO_rel_chIso_MTD_2_Sig_EB_,
3889                               meEleISO_rel_chIso_MTD_3_Sig_EB_,
3890                               meEleISO_rel_chIso_MTD_4_Sig_EB_,
3891                               meEleISO_rel_chIso_MTD_5_Sig_EB_,
3892                               meEleISO_rel_chIso_MTD_6_Sig_EB_,
3893                               meEleISO_rel_chIso_MTD_7_Sig_EB_};
3894   }
3895   Ntracks_EB_list_Significance_Sig = {
3896       meEleISO_Ntracks_MTD_4sigma_Sig_EB_, meEleISO_Ntracks_MTD_3sigma_Sig_EB_, meEleISO_Ntracks_MTD_2sigma_Sig_EB_};
3897   ch_iso_EB_list_Significance_Sig = {
3898       meEleISO_chIso_MTD_4sigma_Sig_EB_, meEleISO_chIso_MTD_3sigma_Sig_EB_, meEleISO_chIso_MTD_2sigma_Sig_EB_};
3899   rel_ch_iso_EB_list_Significance_Sig = {meEleISO_rel_chIso_MTD_4sigma_Sig_EB_,
3900                                          meEleISO_rel_chIso_MTD_3sigma_Sig_EB_,
3901                                          meEleISO_rel_chIso_MTD_2sigma_Sig_EB_};
3902 
3903   if (optionalPlots_) {
3904     Ntracks_EE_list_Sig = {meEleISO_Ntracks_MTD_1_Sig_EE_,
3905                            meEleISO_Ntracks_MTD_2_Sig_EE_,
3906                            meEleISO_Ntracks_MTD_3_Sig_EE_,
3907                            meEleISO_Ntracks_MTD_4_Sig_EE_,
3908                            meEleISO_Ntracks_MTD_5_Sig_EE_,
3909                            meEleISO_Ntracks_MTD_6_Sig_EE_,
3910                            meEleISO_Ntracks_MTD_7_Sig_EE_};
3911     ch_iso_EE_list_Sig = {meEleISO_chIso_MTD_1_Sig_EE_,
3912                           meEleISO_chIso_MTD_2_Sig_EE_,
3913                           meEleISO_chIso_MTD_3_Sig_EE_,
3914                           meEleISO_chIso_MTD_4_Sig_EE_,
3915                           meEleISO_chIso_MTD_5_Sig_EE_,
3916                           meEleISO_chIso_MTD_6_Sig_EE_,
3917                           meEleISO_chIso_MTD_7_Sig_EE_};
3918     rel_ch_iso_EE_list_Sig = {meEleISO_rel_chIso_MTD_1_Sig_EE_,
3919                               meEleISO_rel_chIso_MTD_2_Sig_EE_,
3920                               meEleISO_rel_chIso_MTD_3_Sig_EE_,
3921                               meEleISO_rel_chIso_MTD_4_Sig_EE_,
3922                               meEleISO_rel_chIso_MTD_5_Sig_EE_,
3923                               meEleISO_rel_chIso_MTD_6_Sig_EE_,
3924                               meEleISO_rel_chIso_MTD_7_Sig_EE_};
3925   }
3926   Ntracks_EE_list_Significance_Sig = {
3927       meEleISO_Ntracks_MTD_4sigma_Sig_EE_, meEleISO_Ntracks_MTD_3sigma_Sig_EE_, meEleISO_Ntracks_MTD_2sigma_Sig_EE_};
3928   ch_iso_EE_list_Significance_Sig = {
3929       meEleISO_chIso_MTD_4sigma_Sig_EE_, meEleISO_chIso_MTD_3sigma_Sig_EE_, meEleISO_chIso_MTD_2sigma_Sig_EE_};
3930   rel_ch_iso_EE_list_Significance_Sig = {meEleISO_rel_chIso_MTD_4sigma_Sig_EE_,
3931                                          meEleISO_rel_chIso_MTD_3sigma_Sig_EE_,
3932                                          meEleISO_rel_chIso_MTD_2sigma_Sig_EE_};
3933 
3934   if (optionalPlots_) {
3935     Ele_pT_MTD_EB_list_Sig = {meEle_pt_MTD_1_Sig_EB_,
3936                               meEle_pt_MTD_2_Sig_EB_,
3937                               meEle_pt_MTD_3_Sig_EB_,
3938                               meEle_pt_MTD_4_Sig_EB_,
3939                               meEle_pt_MTD_5_Sig_EB_,
3940                               meEle_pt_MTD_6_Sig_EB_,
3941                               meEle_pt_MTD_7_Sig_EB_};
3942     Ele_eta_MTD_EB_list_Sig = {meEle_eta_MTD_1_Sig_EB_,
3943                                meEle_eta_MTD_2_Sig_EB_,
3944                                meEle_eta_MTD_3_Sig_EB_,
3945                                meEle_eta_MTD_4_Sig_EB_,
3946                                meEle_eta_MTD_5_Sig_EB_,
3947                                meEle_eta_MTD_6_Sig_EB_,
3948                                meEle_eta_MTD_7_Sig_EB_};
3949     Ele_phi_MTD_EB_list_Sig = {meEle_phi_MTD_1_Sig_EB_,
3950                                meEle_phi_MTD_2_Sig_EB_,
3951                                meEle_phi_MTD_3_Sig_EB_,
3952                                meEle_phi_MTD_4_Sig_EB_,
3953                                meEle_phi_MTD_5_Sig_EB_,
3954                                meEle_phi_MTD_6_Sig_EB_,
3955                                meEle_phi_MTD_7_Sig_EB_};
3956   }
3957 
3958   Ele_pT_MTD_EB_list_Significance_Sig = {
3959       meEle_pt_MTD_4sigma_Sig_EB_, meEle_pt_MTD_3sigma_Sig_EB_, meEle_pt_MTD_2sigma_Sig_EB_};
3960   Ele_eta_MTD_EB_list_Significance_Sig = {
3961       meEle_eta_MTD_4sigma_Sig_EB_, meEle_eta_MTD_3sigma_Sig_EB_, meEle_eta_MTD_2sigma_Sig_EB_};
3962   Ele_phi_MTD_EB_list_Significance_Sig = {
3963       meEle_phi_MTD_4sigma_Sig_EB_, meEle_phi_MTD_3sigma_Sig_EB_, meEle_phi_MTD_2sigma_Sig_EB_};
3964 
3965   if (optionalPlots_) {
3966     Ele_pT_MTD_EE_list_Sig = {meEle_pt_MTD_1_Sig_EE_,
3967                               meEle_pt_MTD_2_Sig_EE_,
3968                               meEle_pt_MTD_3_Sig_EE_,
3969                               meEle_pt_MTD_4_Sig_EE_,
3970                               meEle_pt_MTD_5_Sig_EE_,
3971                               meEle_pt_MTD_6_Sig_EE_,
3972                               meEle_pt_MTD_7_Sig_EE_};
3973     Ele_eta_MTD_EE_list_Sig = {meEle_eta_MTD_1_Sig_EE_,
3974                                meEle_eta_MTD_2_Sig_EE_,
3975                                meEle_eta_MTD_3_Sig_EE_,
3976                                meEle_eta_MTD_4_Sig_EE_,
3977                                meEle_eta_MTD_5_Sig_EE_,
3978                                meEle_eta_MTD_6_Sig_EE_,
3979                                meEle_eta_MTD_7_Sig_EE_};
3980     Ele_phi_MTD_EE_list_Sig = {meEle_phi_MTD_1_Sig_EE_,
3981                                meEle_phi_MTD_2_Sig_EE_,
3982                                meEle_phi_MTD_3_Sig_EE_,
3983                                meEle_phi_MTD_4_Sig_EE_,
3984                                meEle_phi_MTD_5_Sig_EE_,
3985                                meEle_phi_MTD_6_Sig_EE_,
3986                                meEle_phi_MTD_7_Sig_EE_};
3987   }
3988   Ele_pT_MTD_EE_list_Significance_Sig = {
3989       meEle_pt_MTD_4sigma_Sig_EE_, meEle_pt_MTD_3sigma_Sig_EE_, meEle_pt_MTD_2sigma_Sig_EE_};
3990   Ele_eta_MTD_EE_list_Significance_Sig = {
3991       meEle_eta_MTD_4sigma_Sig_EE_, meEle_eta_MTD_3sigma_Sig_EE_, meEle_eta_MTD_2sigma_Sig_EE_};
3992   Ele_phi_MTD_EE_list_Significance_Sig = {
3993       meEle_phi_MTD_4sigma_Sig_EE_, meEle_phi_MTD_3sigma_Sig_EE_, meEle_phi_MTD_2sigma_Sig_EE_};
3994 
3995   // For SIM CASE
3996   if (optionalPlots_) {
3997     Ntracks_sim_EB_list_Sig = {meEleISO_Ntracks_MTD_sim_1_Sig_EB_,
3998                                meEleISO_Ntracks_MTD_sim_2_Sig_EB_,
3999                                meEleISO_Ntracks_MTD_sim_3_Sig_EB_,
4000                                meEleISO_Ntracks_MTD_sim_4_Sig_EB_,
4001                                meEleISO_Ntracks_MTD_sim_5_Sig_EB_,
4002                                meEleISO_Ntracks_MTD_sim_6_Sig_EB_,
4003                                meEleISO_Ntracks_MTD_sim_7_Sig_EB_};
4004     ch_iso_sim_EB_list_Sig = {meEleISO_chIso_MTD_sim_1_Sig_EB_,
4005                               meEleISO_chIso_MTD_sim_2_Sig_EB_,
4006                               meEleISO_chIso_MTD_sim_3_Sig_EB_,
4007                               meEleISO_chIso_MTD_sim_4_Sig_EB_,
4008                               meEleISO_chIso_MTD_sim_5_Sig_EB_,
4009                               meEleISO_chIso_MTD_sim_6_Sig_EB_,
4010                               meEleISO_chIso_MTD_sim_7_Sig_EB_};
4011     rel_ch_iso_sim_EB_list_Sig = {meEleISO_rel_chIso_MTD_sim_1_Sig_EB_,
4012                                   meEleISO_rel_chIso_MTD_sim_2_Sig_EB_,
4013                                   meEleISO_rel_chIso_MTD_sim_3_Sig_EB_,
4014                                   meEleISO_rel_chIso_MTD_sim_4_Sig_EB_,
4015                                   meEleISO_rel_chIso_MTD_sim_5_Sig_EB_,
4016                                   meEleISO_rel_chIso_MTD_sim_6_Sig_EB_,
4017                                   meEleISO_rel_chIso_MTD_sim_7_Sig_EB_};
4018 
4019     Ntracks_sim_EB_list_Significance_Sig = {meEleISO_Ntracks_MTD_sim_4sigma_Sig_EB_,
4020                                             meEleISO_Ntracks_MTD_sim_3sigma_Sig_EB_,
4021                                             meEleISO_Ntracks_MTD_sim_2sigma_Sig_EB_};
4022     ch_iso_sim_EB_list_Significance_Sig = {meEleISO_chIso_MTD_sim_4sigma_Sig_EB_,
4023                                            meEleISO_chIso_MTD_sim_3sigma_Sig_EB_,
4024                                            meEleISO_chIso_MTD_sim_2sigma_Sig_EB_};
4025     rel_ch_iso_sim_EB_list_Significance_Sig = {meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EB_,
4026                                                meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EB_,
4027                                                meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EB_};
4028 
4029     Ntracks_sim_EE_list_Sig = {meEleISO_Ntracks_MTD_sim_1_Sig_EE_,
4030                                meEleISO_Ntracks_MTD_sim_2_Sig_EE_,
4031                                meEleISO_Ntracks_MTD_sim_3_Sig_EE_,
4032                                meEleISO_Ntracks_MTD_sim_4_Sig_EE_,
4033                                meEleISO_Ntracks_MTD_sim_5_Sig_EE_,
4034                                meEleISO_Ntracks_MTD_sim_6_Sig_EE_,
4035                                meEleISO_Ntracks_MTD_sim_7_Sig_EE_};
4036     ch_iso_sim_EE_list_Sig = {meEleISO_chIso_MTD_sim_1_Sig_EE_,
4037                               meEleISO_chIso_MTD_sim_2_Sig_EE_,
4038                               meEleISO_chIso_MTD_sim_3_Sig_EE_,
4039                               meEleISO_chIso_MTD_sim_4_Sig_EE_,
4040                               meEleISO_chIso_MTD_sim_5_Sig_EE_,
4041                               meEleISO_chIso_MTD_sim_6_Sig_EE_,
4042                               meEleISO_chIso_MTD_sim_7_Sig_EE_};
4043     rel_ch_iso_sim_EE_list_Sig = {meEleISO_rel_chIso_MTD_sim_1_Sig_EE_,
4044                                   meEleISO_rel_chIso_MTD_sim_2_Sig_EE_,
4045                                   meEleISO_rel_chIso_MTD_sim_3_Sig_EE_,
4046                                   meEleISO_rel_chIso_MTD_sim_4_Sig_EE_,
4047                                   meEleISO_rel_chIso_MTD_sim_5_Sig_EE_,
4048                                   meEleISO_rel_chIso_MTD_sim_6_Sig_EE_,
4049                                   meEleISO_rel_chIso_MTD_sim_7_Sig_EE_};
4050 
4051     Ntracks_sim_EE_list_Significance_Sig = {meEleISO_Ntracks_MTD_sim_4sigma_Sig_EE_,
4052                                             meEleISO_Ntracks_MTD_sim_3sigma_Sig_EE_,
4053                                             meEleISO_Ntracks_MTD_sim_2sigma_Sig_EE_};
4054     ch_iso_sim_EE_list_Significance_Sig = {meEleISO_chIso_MTD_sim_4sigma_Sig_EE_,
4055                                            meEleISO_chIso_MTD_sim_3sigma_Sig_EE_,
4056                                            meEleISO_chIso_MTD_sim_2sigma_Sig_EE_};
4057     rel_ch_iso_sim_EE_list_Significance_Sig = {meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EE_,
4058                                                meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EE_,
4059                                                meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EE_};
4060 
4061     Ele_pT_sim_MTD_EB_list_Sig = {meEle_pt_sim_MTD_1_Sig_EB_,
4062                                   meEle_pt_sim_MTD_2_Sig_EB_,
4063                                   meEle_pt_sim_MTD_3_Sig_EB_,
4064                                   meEle_pt_sim_MTD_4_Sig_EB_,
4065                                   meEle_pt_sim_MTD_5_Sig_EB_,
4066                                   meEle_pt_sim_MTD_6_Sig_EB_,
4067                                   meEle_pt_sim_MTD_7_Sig_EB_};
4068 
4069     Ele_pT_sim_MTD_EB_list_Significance_Sig = {
4070         meEle_pt_sim_MTD_4sigma_Sig_EB_, meEle_pt_sim_MTD_3sigma_Sig_EB_, meEle_pt_sim_MTD_2sigma_Sig_EB_};
4071 
4072     Ele_pT_sim_MTD_EE_list_Sig = {meEle_pt_sim_MTD_1_Sig_EE_,
4073                                   meEle_pt_sim_MTD_2_Sig_EE_,
4074                                   meEle_pt_sim_MTD_3_Sig_EE_,
4075                                   meEle_pt_sim_MTD_4_Sig_EE_,
4076                                   meEle_pt_sim_MTD_5_Sig_EE_,
4077                                   meEle_pt_sim_MTD_6_Sig_EE_,
4078                                   meEle_pt_sim_MTD_7_Sig_EE_};
4079     Ele_pT_sim_MTD_EE_list_Significance_Sig = {
4080         meEle_pt_sim_MTD_4sigma_Sig_EE_, meEle_pt_sim_MTD_3sigma_Sig_EE_, meEle_pt_sim_MTD_2sigma_Sig_EE_};
4081   }
4082 
4083   // Non-promt part
4084   if (optionalPlots_) {
4085     Ntracks_EB_list_Bkg = {meEleISO_Ntracks_MTD_1_Bkg_EB_,
4086                            meEleISO_Ntracks_MTD_2_Bkg_EB_,
4087                            meEleISO_Ntracks_MTD_3_Bkg_EB_,
4088                            meEleISO_Ntracks_MTD_4_Bkg_EB_,
4089                            meEleISO_Ntracks_MTD_5_Bkg_EB_,
4090                            meEleISO_Ntracks_MTD_6_Bkg_EB_,
4091                            meEleISO_Ntracks_MTD_7_Bkg_EB_};
4092     ch_iso_EB_list_Bkg = {meEleISO_chIso_MTD_1_Bkg_EB_,
4093                           meEleISO_chIso_MTD_2_Bkg_EB_,
4094                           meEleISO_chIso_MTD_3_Bkg_EB_,
4095                           meEleISO_chIso_MTD_4_Bkg_EB_,
4096                           meEleISO_chIso_MTD_5_Bkg_EB_,
4097                           meEleISO_chIso_MTD_6_Bkg_EB_,
4098                           meEleISO_chIso_MTD_7_Bkg_EB_};
4099     rel_ch_iso_EB_list_Bkg = {meEleISO_rel_chIso_MTD_1_Bkg_EB_,
4100                               meEleISO_rel_chIso_MTD_2_Bkg_EB_,
4101                               meEleISO_rel_chIso_MTD_3_Bkg_EB_,
4102                               meEleISO_rel_chIso_MTD_4_Bkg_EB_,
4103                               meEleISO_rel_chIso_MTD_5_Bkg_EB_,
4104                               meEleISO_rel_chIso_MTD_6_Bkg_EB_,
4105                               meEleISO_rel_chIso_MTD_7_Bkg_EB_};
4106   }
4107   Ntracks_EB_list_Significance_Bkg = {
4108       meEleISO_Ntracks_MTD_4sigma_Bkg_EB_, meEleISO_Ntracks_MTD_3sigma_Bkg_EB_, meEleISO_Ntracks_MTD_2sigma_Bkg_EB_};
4109   ch_iso_EB_list_Significance_Bkg = {
4110       meEleISO_chIso_MTD_4sigma_Bkg_EB_, meEleISO_chIso_MTD_3sigma_Bkg_EB_, meEleISO_chIso_MTD_2sigma_Bkg_EB_};
4111   rel_ch_iso_EB_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_4sigma_Bkg_EB_,
4112                                          meEleISO_rel_chIso_MTD_3sigma_Bkg_EB_,
4113                                          meEleISO_rel_chIso_MTD_2sigma_Bkg_EB_};
4114 
4115   if (optionalPlots_) {
4116     Ntracks_EE_list_Bkg = {meEleISO_Ntracks_MTD_1_Bkg_EE_,
4117                            meEleISO_Ntracks_MTD_2_Bkg_EE_,
4118                            meEleISO_Ntracks_MTD_3_Bkg_EE_,
4119                            meEleISO_Ntracks_MTD_4_Bkg_EE_,
4120                            meEleISO_Ntracks_MTD_5_Bkg_EE_,
4121                            meEleISO_Ntracks_MTD_6_Bkg_EE_,
4122                            meEleISO_Ntracks_MTD_7_Bkg_EE_};
4123     ch_iso_EE_list_Bkg = {meEleISO_chIso_MTD_1_Bkg_EE_,
4124                           meEleISO_chIso_MTD_2_Bkg_EE_,
4125                           meEleISO_chIso_MTD_3_Bkg_EE_,
4126                           meEleISO_chIso_MTD_4_Bkg_EE_,
4127                           meEleISO_chIso_MTD_5_Bkg_EE_,
4128                           meEleISO_chIso_MTD_6_Bkg_EE_,
4129                           meEleISO_chIso_MTD_7_Bkg_EE_};
4130     rel_ch_iso_EE_list_Bkg = {meEleISO_rel_chIso_MTD_1_Bkg_EE_,
4131                               meEleISO_rel_chIso_MTD_2_Bkg_EE_,
4132                               meEleISO_rel_chIso_MTD_3_Bkg_EE_,
4133                               meEleISO_rel_chIso_MTD_4_Bkg_EE_,
4134                               meEleISO_rel_chIso_MTD_5_Bkg_EE_,
4135                               meEleISO_rel_chIso_MTD_6_Bkg_EE_,
4136                               meEleISO_rel_chIso_MTD_7_Bkg_EE_};
4137   }
4138   Ntracks_EE_list_Significance_Bkg = {
4139       meEleISO_Ntracks_MTD_4sigma_Bkg_EE_, meEleISO_Ntracks_MTD_3sigma_Bkg_EE_, meEleISO_Ntracks_MTD_2sigma_Bkg_EE_};
4140   ch_iso_EE_list_Significance_Bkg = {
4141       meEleISO_chIso_MTD_4sigma_Bkg_EE_, meEleISO_chIso_MTD_3sigma_Bkg_EE_, meEleISO_chIso_MTD_2sigma_Bkg_EE_};
4142   rel_ch_iso_EE_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_4sigma_Bkg_EE_,
4143                                          meEleISO_rel_chIso_MTD_3sigma_Bkg_EE_,
4144                                          meEleISO_rel_chIso_MTD_2sigma_Bkg_EE_};
4145   if (optionalPlots_) {
4146     Ele_pT_MTD_EB_list_Bkg = {meEle_pt_MTD_1_Bkg_EB_,
4147                               meEle_pt_MTD_2_Bkg_EB_,
4148                               meEle_pt_MTD_3_Bkg_EB_,
4149                               meEle_pt_MTD_4_Bkg_EB_,
4150                               meEle_pt_MTD_5_Bkg_EB_,
4151                               meEle_pt_MTD_6_Bkg_EB_,
4152                               meEle_pt_MTD_7_Bkg_EB_};
4153     Ele_eta_MTD_EB_list_Bkg = {meEle_eta_MTD_1_Bkg_EB_,
4154                                meEle_eta_MTD_2_Bkg_EB_,
4155                                meEle_eta_MTD_3_Bkg_EB_,
4156                                meEle_eta_MTD_4_Bkg_EB_,
4157                                meEle_eta_MTD_5_Bkg_EB_,
4158                                meEle_eta_MTD_6_Bkg_EB_,
4159                                meEle_eta_MTD_7_Bkg_EB_};
4160     Ele_phi_MTD_EB_list_Bkg = {meEle_phi_MTD_1_Bkg_EB_,
4161                                meEle_phi_MTD_2_Bkg_EB_,
4162                                meEle_phi_MTD_3_Bkg_EB_,
4163                                meEle_phi_MTD_4_Bkg_EB_,
4164                                meEle_phi_MTD_5_Bkg_EB_,
4165                                meEle_phi_MTD_6_Bkg_EB_,
4166                                meEle_phi_MTD_7_Bkg_EB_};
4167   }
4168   Ele_pT_MTD_EB_list_Significance_Bkg = {
4169       meEle_pt_MTD_4sigma_Bkg_EB_, meEle_pt_MTD_3sigma_Bkg_EB_, meEle_pt_MTD_2sigma_Bkg_EB_};
4170   Ele_eta_MTD_EB_list_Significance_Bkg = {
4171       meEle_eta_MTD_4sigma_Bkg_EB_, meEle_eta_MTD_3sigma_Bkg_EB_, meEle_eta_MTD_2sigma_Bkg_EB_};
4172   Ele_phi_MTD_EB_list_Significance_Bkg = {
4173       meEle_phi_MTD_4sigma_Bkg_EB_, meEle_phi_MTD_3sigma_Bkg_EB_, meEle_phi_MTD_2sigma_Bkg_EB_};
4174 
4175   if (optionalPlots_) {
4176     Ele_pT_MTD_EE_list_Bkg = {meEle_pt_MTD_1_Bkg_EE_,
4177                               meEle_pt_MTD_2_Bkg_EE_,
4178                               meEle_pt_MTD_3_Bkg_EE_,
4179                               meEle_pt_MTD_4_Bkg_EE_,
4180                               meEle_pt_MTD_5_Bkg_EE_,
4181                               meEle_pt_MTD_6_Bkg_EE_,
4182                               meEle_pt_MTD_7_Bkg_EE_};
4183     Ele_eta_MTD_EE_list_Bkg = {meEle_eta_MTD_1_Bkg_EE_,
4184                                meEle_eta_MTD_2_Bkg_EE_,
4185                                meEle_eta_MTD_3_Bkg_EE_,
4186                                meEle_eta_MTD_4_Bkg_EE_,
4187                                meEle_eta_MTD_5_Bkg_EE_,
4188                                meEle_eta_MTD_6_Bkg_EE_,
4189                                meEle_eta_MTD_7_Bkg_EE_};
4190     Ele_phi_MTD_EE_list_Bkg = {meEle_phi_MTD_1_Bkg_EE_,
4191                                meEle_phi_MTD_2_Bkg_EE_,
4192                                meEle_phi_MTD_3_Bkg_EE_,
4193                                meEle_phi_MTD_4_Bkg_EE_,
4194                                meEle_phi_MTD_5_Bkg_EE_,
4195                                meEle_phi_MTD_6_Bkg_EE_,
4196                                meEle_phi_MTD_7_Bkg_EE_};
4197   }
4198   Ele_pT_MTD_EE_list_Significance_Bkg = {
4199       meEle_pt_MTD_4sigma_Bkg_EE_, meEle_pt_MTD_3sigma_Bkg_EE_, meEle_pt_MTD_2sigma_Bkg_EE_};
4200   Ele_eta_MTD_EE_list_Significance_Bkg = {
4201       meEle_eta_MTD_4sigma_Bkg_EE_, meEle_eta_MTD_3sigma_Bkg_EE_, meEle_eta_MTD_2sigma_Bkg_EE_};
4202   Ele_phi_MTD_EE_list_Significance_Bkg = {
4203       meEle_phi_MTD_4sigma_Bkg_EE_, meEle_phi_MTD_3sigma_Bkg_EE_, meEle_phi_MTD_2sigma_Bkg_EE_};
4204 
4205   // SIM CASE
4206   if (optionalPlots_) {
4207     Ntracks_sim_EB_list_Bkg = {meEleISO_Ntracks_MTD_sim_1_Bkg_EB_,
4208                                meEleISO_Ntracks_MTD_sim_2_Bkg_EB_,
4209                                meEleISO_Ntracks_MTD_sim_3_Bkg_EB_,
4210                                meEleISO_Ntracks_MTD_sim_4_Bkg_EB_,
4211                                meEleISO_Ntracks_MTD_sim_5_Bkg_EB_,
4212                                meEleISO_Ntracks_MTD_sim_6_Bkg_EB_,
4213                                meEleISO_Ntracks_MTD_sim_7_Bkg_EB_};
4214     ch_iso_sim_EB_list_Bkg = {meEleISO_chIso_MTD_sim_1_Bkg_EB_,
4215                               meEleISO_chIso_MTD_sim_2_Bkg_EB_,
4216                               meEleISO_chIso_MTD_sim_3_Bkg_EB_,
4217                               meEleISO_chIso_MTD_sim_4_Bkg_EB_,
4218                               meEleISO_chIso_MTD_sim_5_Bkg_EB_,
4219                               meEleISO_chIso_MTD_sim_6_Bkg_EB_,
4220                               meEleISO_chIso_MTD_sim_7_Bkg_EB_};
4221     rel_ch_iso_sim_EB_list_Bkg = {meEleISO_rel_chIso_MTD_sim_1_Bkg_EB_,
4222                                   meEleISO_rel_chIso_MTD_sim_2_Bkg_EB_,
4223                                   meEleISO_rel_chIso_MTD_sim_3_Bkg_EB_,
4224                                   meEleISO_rel_chIso_MTD_sim_4_Bkg_EB_,
4225                                   meEleISO_rel_chIso_MTD_sim_5_Bkg_EB_,
4226                                   meEleISO_rel_chIso_MTD_sim_6_Bkg_EB_,
4227                                   meEleISO_rel_chIso_MTD_sim_7_Bkg_EB_};
4228     Ntracks_sim_EB_list_Significance_Bkg = {meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EB_,
4229                                             meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EB_,
4230                                             meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EB_};
4231     ch_iso_sim_EB_list_Significance_Bkg = {meEleISO_chIso_MTD_sim_4sigma_Bkg_EB_,
4232                                            meEleISO_chIso_MTD_sim_3sigma_Bkg_EB_,
4233                                            meEleISO_chIso_MTD_sim_2sigma_Bkg_EB_};
4234     rel_ch_iso_sim_EB_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EB_,
4235                                                meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EB_,
4236                                                meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EB_};
4237 
4238     Ntracks_sim_EE_list_Bkg = {meEleISO_Ntracks_MTD_sim_1_Bkg_EE_,
4239                                meEleISO_Ntracks_MTD_sim_2_Bkg_EE_,
4240                                meEleISO_Ntracks_MTD_sim_3_Bkg_EE_,
4241                                meEleISO_Ntracks_MTD_sim_4_Bkg_EE_,
4242                                meEleISO_Ntracks_MTD_sim_5_Bkg_EE_,
4243                                meEleISO_Ntracks_MTD_sim_6_Bkg_EE_,
4244                                meEleISO_Ntracks_MTD_sim_7_Bkg_EE_};
4245     ch_iso_sim_EE_list_Bkg = {meEleISO_chIso_MTD_sim_1_Bkg_EE_,
4246                               meEleISO_chIso_MTD_sim_2_Bkg_EE_,
4247                               meEleISO_chIso_MTD_sim_3_Bkg_EE_,
4248                               meEleISO_chIso_MTD_sim_4_Bkg_EE_,
4249                               meEleISO_chIso_MTD_sim_5_Bkg_EE_,
4250                               meEleISO_chIso_MTD_sim_6_Bkg_EE_,
4251                               meEleISO_chIso_MTD_sim_7_Bkg_EE_};
4252     rel_ch_iso_sim_EE_list_Bkg = {meEleISO_rel_chIso_MTD_sim_1_Bkg_EE_,
4253                                   meEleISO_rel_chIso_MTD_sim_2_Bkg_EE_,
4254                                   meEleISO_rel_chIso_MTD_sim_3_Bkg_EE_,
4255                                   meEleISO_rel_chIso_MTD_sim_4_Bkg_EE_,
4256                                   meEleISO_rel_chIso_MTD_sim_5_Bkg_EE_,
4257                                   meEleISO_rel_chIso_MTD_sim_6_Bkg_EE_,
4258                                   meEleISO_rel_chIso_MTD_sim_7_Bkg_EE_};
4259 
4260     Ntracks_sim_EE_list_Significance_Bkg = {meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EE_,
4261                                             meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EE_,
4262                                             meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EE_};
4263     ch_iso_sim_EE_list_Significance_Bkg = {meEleISO_chIso_MTD_sim_4sigma_Bkg_EE_,
4264                                            meEleISO_chIso_MTD_sim_3sigma_Bkg_EE_,
4265                                            meEleISO_chIso_MTD_sim_2sigma_Bkg_EE_};
4266     rel_ch_iso_sim_EE_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EE_,
4267                                                meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EE_,
4268                                                meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EE_};
4269 
4270     Ele_pT_sim_MTD_EB_list_Bkg = {meEle_pt_sim_MTD_1_Bkg_EB_,
4271                                   meEle_pt_sim_MTD_2_Bkg_EB_,
4272                                   meEle_pt_sim_MTD_3_Bkg_EB_,
4273                                   meEle_pt_sim_MTD_4_Bkg_EB_,
4274                                   meEle_pt_sim_MTD_5_Bkg_EB_,
4275                                   meEle_pt_sim_MTD_6_Bkg_EB_,
4276                                   meEle_pt_sim_MTD_7_Bkg_EB_};
4277 
4278     Ele_pT_sim_MTD_EB_list_Significance_Bkg = {
4279         meEle_pt_sim_MTD_4sigma_Bkg_EB_, meEle_pt_sim_MTD_3sigma_Bkg_EB_, meEle_pt_sim_MTD_2sigma_Bkg_EB_};
4280 
4281     Ele_pT_sim_MTD_EE_list_Bkg = {meEle_pt_sim_MTD_1_Bkg_EE_,
4282                                   meEle_pt_sim_MTD_2_Bkg_EE_,
4283                                   meEle_pt_sim_MTD_3_Bkg_EE_,
4284                                   meEle_pt_sim_MTD_4_Bkg_EE_,
4285                                   meEle_pt_sim_MTD_5_Bkg_EE_,
4286                                   meEle_pt_sim_MTD_6_Bkg_EE_,
4287                                   meEle_pt_sim_MTD_7_Bkg_EE_};
4288 
4289     Ele_pT_sim_MTD_EE_list_Significance_Bkg = {
4290         meEle_pt_sim_MTD_4sigma_Bkg_EE_, meEle_pt_sim_MTD_3sigma_Bkg_EE_, meEle_pt_sim_MTD_2sigma_Bkg_EE_};
4291   }
4292 }
4293 
4294 // ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
4295 
4296 void MtdEleIsoValidation::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
4297   edm::ParameterSetDescription desc;
4298 
4299   desc.add<std::string>("folder", "MTD/ElectronIso");
4300   desc.add<edm::InputTag>("inputTagG", edm::InputTag("generalTracks"));
4301   desc.add<edm::InputTag>("inputTag_vtx", edm::InputTag("offlinePrimaryVertices4D"));
4302   desc.add<edm::InputTag>("inputEle_EB", edm::InputTag("gedGsfElectrons"));
4303   desc.add<edm::InputTag>("inputEle_EE", edm::InputTag("ecalDrivenGsfElectronsHGC"));
4304   desc.add<edm::InputTag>("TPtoRecoTrackAssoc", edm::InputTag("trackingParticleRecoTrackAsssociation"));
4305   desc.add<edm::InputTag>("t0PID", edm::InputTag("tofPID:t0"));
4306   desc.add<edm::InputTag>("sigmat0PID", edm::InputTag("tofPID:sigmat0"));
4307   desc.add<edm::InputTag>("trackMVAQual", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
4308   desc.add<double>("trackMinimumPt", 1.0);  // [GeV]
4309   desc.add<double>("trackMinimumEta", 1.5);
4310   desc.add<double>("trackMaximumEta", 3.2);
4311   desc.add<double>("rel_iso_cut", 0.08);
4312   desc.add<bool>("optionTrackMatchToPV", false);
4313   desc.add<bool>("option_dtToTrack", true);  // default is dt with track, if false will do dt to vertex
4314   desc.add<bool>("option_plots", false);
4315   desc.add<double>("min_dR_cut", 0.01);
4316   desc.add<double>("max_dR_cut", 0.3);
4317   desc.add<double>("min_pt_cut_EB", 0.7);
4318   desc.add<double>("min_pt_cut_EE", 0.4);
4319   desc.add<double>("max_dz_cut_EB", 0.5);  // PARAM
4320   desc.add<double>("max_dz_cut_EE", 0.5);  // PARAM
4321   desc.add<double>("max_dz_vtx_cut", 0.5);
4322   desc.add<double>("max_dxy_vtx_cut", 0.2);
4323   desc.add<double>("min_strip_cut", 0.01);
4324   desc.add<double>("min_track_mtd_mva_cut", 0.5);
4325 
4326   descriptions.add("mtdEleIsoValid", desc);
4327 }
4328 
4329 DEFINE_FWK_MODULE(MtdEleIsoValidation);
4330 
4331 //*/