Back to home page

Project CMSSW displayed by LXR
 
 

    

File indexing completed on 2024-05-10 02:21:37 

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()) {
0973           if (genParticle.mother() != nullptr and 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 
0981     math::XYZVector EleSigTrackMomentumAtVtx = ele.gsfTrack()->momentum();
0982     double EleSigTrackEtaAtVtx = ele.gsfTrack()->eta();
0983 
0984     double ele_sigTrkTime = -1;
0985     double ele_sigTrkTimeErr = -1;
0986     double ele_sigTrkMtdMva = -1;
0987 
0988     // if we found a track match, we add MTD timing information for it
0989     if (ele_TrkRef.isNonnull()) {
0990       // track pT/dz cuts
0991       bool Barrel_ele = ele.isEB();
0992       float min_pt_cut = Barrel_ele ? min_pt_cut_EB : min_pt_cut_EE;
0993       float max_dz_cut = Barrel_ele ? max_dz_cut_EB : max_dz_cut_EE;
0994 
0995       ele_sigTrkTime = t0Pid[ele_TrkRef];
0996       ele_sigTrkMtdMva = mtdQualMVA[ele_TrkRef];
0997       ele_sigTrkTimeErr = (ele_sigTrkMtdMva > min_track_mtd_mva_cut) ? Sigmat0Pid[ele_TrkRef] : -1;
0998 
0999       meEle_avg_error_SigTrk_check_->Fill(ele_sigTrkTimeErr);
1000 
1001       if (ele_Promt) {
1002         // For signal (promt)
1003         if (Barrel_ele) {
1004           // All selected electron information for efficiency plots later
1005           meEle_pt_tot_Sig_EB_->Fill(ele.pt());
1006           meEle_pt_sim_tot_Sig_EB_->Fill(ele_sim_pt);
1007           meEle_eta_tot_Sig_EB_->Fill(std::abs(ele.eta()));
1008           meEle_phi_tot_Sig_EB_->Fill(ele.phi());
1009         } else {
1010           // All selected electron information for efficiency plots later
1011           meEle_pt_tot_Sig_EE_->Fill(ele.pt());
1012           meEle_pt_sim_tot_Sig_EE_->Fill(ele_sim_pt);
1013           meEle_eta_tot_Sig_EE_->Fill(std::abs(ele.eta()));
1014           meEle_phi_tot_Sig_EE_->Fill(ele.phi());
1015         }
1016       } else {
1017         // For background (non-promt)
1018         if (Barrel_ele) {
1019           meEle_pt_tot_Bkg_EB_->Fill(ele.pt());
1020           meEle_pt_sim_tot_Bkg_EB_->Fill(ele_sim_pt);
1021           meEle_eta_tot_Bkg_EB_->Fill(std::abs(ele.eta()));
1022           meEle_phi_tot_Bkg_EB_->Fill(ele.phi());
1023         } else {
1024           meEle_pt_tot_Bkg_EE_->Fill(ele.pt());
1025           meEle_pt_sim_tot_Bkg_EE_->Fill(ele_sim_pt);
1026           meEle_eta_tot_Bkg_EE_->Fill(std::abs(ele.eta()));
1027           meEle_phi_tot_Bkg_EE_->Fill(ele.phi());
1028         }
1029       }
1030 
1031       int N_tracks_noMTD = 0;
1032       double pT_sum_noMTD = 0;
1033       double rel_pT_sum_noMTD = 0;
1034       std::vector<int> N_tracks_MTD{0, 0, 0, 0, 0, 0, 0};
1035       std::vector<double> pT_sum_MTD{0, 0, 0, 0, 0, 0, 0};
1036       std::vector<double> rel_pT_sum_MTD{0, 0, 0, 0, 0, 0, 0};
1037 
1038       std::vector<int> N_tracks_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1039       std::vector<double> pT_sum_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1040       std::vector<double> rel_pT_sum_sim_MTD{0, 0, 0, 0, 0, 0, 0};
1041       int N_tracks_gen = 0;
1042       double pT_sum_gen = 0;
1043       double rel_pT_sum_gen = 0;
1044 
1045       std::vector<int> N_tracks_MTD_significance{0, 0, 0};
1046       std::vector<double> pT_sum_MTD_significance{0, 0, 0};
1047       std::vector<double> rel_pT_sum_MTD_significance{0, 0, 0};
1048 
1049       std::vector<int> N_tracks_sim_MTD_significance{0, 0, 0};
1050       std::vector<double> pT_sum_sim_MTD_significance{0, 0, 0};
1051       std::vector<double> rel_pT_sum_sim_MTD_significance{0, 0, 0};
1052 
1053       int general_index = 0;
1054       for (const auto& trackGen : *GenRecTrackHandle) {
1055         const reco::TrackRef trackref_general(GenRecTrackHandle, general_index);
1056         general_index++;
1057 
1058         // Skip electron track
1059         if (trackref_general == ele_TrkRef)
1060           continue;
1061 
1062         if (trackGen.pt() < min_pt_cut) {
1063           continue;
1064         }
1065         if (std::abs(trackGen.vz() - ele.gsfTrack()->vz()) > max_dz_cut) {
1066           continue;
1067         }
1068 
1069         // cut for general track matching to PV
1070         if (track_match_PV_) {
1071           if (Vtx_chosen.trackWeight(trackref_general) < 0.5) {
1072             continue;
1073           }
1074         }
1075 
1076         double dR = reco::deltaR(trackGen.momentum(), EleSigTrackMomentumAtVtx);
1077         double deta = std::abs(trackGen.eta() - EleSigTrackEtaAtVtx);
1078 
1079         // restrict to tracks in the isolation cone
1080         if (dR < min_dR_cut || dR > max_dR_cut || deta < min_strip_cut)
1081           continue;
1082 
1083         // no MTD case
1084         ++N_tracks_noMTD;
1085         pT_sum_noMTD += trackGen.pt();
1086 
1087         // MTD case
1088         const reco::TrackBaseRef trkrefBase(trackref_general);
1089         auto TPmatched = r2s_->find(trkrefBase);
1090         double tsim_trk = -1.;
1091         double trk_ptSim = -1.;
1092         bool genMatched = false;
1093         if (TPmatched != r2s_->end()) {
1094           // reco track matched to a TP
1095           const auto& tp = (TPmatched->val)[0];
1096           tsim_trk = (tp.first)->parentVertex()->position().t() * 1e9;
1097           trk_ptSim = (tp.first)->pt();
1098           // check that the genParticle vector is not empty
1099           if (tp.first->status() != -99) {
1100             genMatched = true;
1101             meTrk_genMatch_check_->Fill(1);
1102           } else {
1103             meTrk_genMatch_check_->Fill(0);
1104           }
1105         }
1106 
1107         double TrkMTDTime = t0Pid[trackref_general];
1108         double TrkMTDMva = mtdQualMVA[trackref_general];
1109         double TrkMTDTimeErr = (TrkMTDMva > min_track_mtd_mva_cut) ? Sigmat0Pid[trackref_general] : -1;
1110 
1111         meEle_avg_error_PUTrk_check_->Fill(TrkMTDTimeErr);
1112 
1113         // MTD GEN case
1114         if (genMatched) {
1115           N_tracks_gen++;
1116           pT_sum_gen += trk_ptSim;
1117         }
1118 
1119         // dt with the track
1120         if (dt_sig_track_) {
1121           double dt_sigTrk = 0;
1122           double dt_sigTrk_signif = 0;
1123           double dt_sim_sigTrk = 0;
1124           double dt_sim_sigTrk_signif = 0;
1125 
1126           // MTD SIM CASE
1127           if (std::abs(tsim_trk) > 0 && std::abs(tsim_ele) > 0 && trk_ptSim > 0) {
1128             dt_sim_sigTrk = std::abs(tsim_trk - tsim_ele);
1129             dt_sim_sigTrk_signif = dt_sim_sigTrk / std::sqrt(avg_sim_PUtrack_t_err * avg_sim_PUtrack_t_err +
1130                                                              avg_sim_sigTrk_t_err * avg_sim_sigTrk_t_err);
1131 
1132             if (optionalPlots_) {
1133               // absolute timing cuts
1134               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1135                 if (dt_sim_sigTrk < max_dt_track_cut[i]) {
1136                   N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1137                   pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1138                 }
1139               }
1140             }
1141             // significance cuts
1142             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1143               if (dt_sim_sigTrk_signif < max_dt_significance_cut[i]) {
1144                 N_tracks_sim_MTD_significance[i]++;
1145                 pT_sum_sim_MTD_significance[i] += trk_ptSim;
1146               }
1147             }
1148 
1149           } else {
1150             // if there is no error for MTD information, we count the MTD isolation case same as noMTD
1151             if (optionalPlots_) {
1152               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1153                 N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1154                 pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1155               }
1156             }
1157             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1158               N_tracks_sim_MTD_significance[i]++;
1159               pT_sum_sim_MTD_significance[i] += trk_ptSim;
1160             }
1161           }
1162 
1163           // MTD reco case
1164           if (TrkMTDTimeErr > 0 && ele_sigTrkTimeErr > 0) {
1165             dt_sigTrk = std::abs(TrkMTDTime - ele_sigTrkTime);
1166             dt_sigTrk_signif =
1167                 dt_sigTrk / std::sqrt(TrkMTDTimeErr * TrkMTDTimeErr + ele_sigTrkTimeErr * ele_sigTrkTimeErr);
1168 
1169             meEle_no_dt_check_->Fill(1);
1170             if (optionalPlots_) {
1171               // absolute timing cuts
1172               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1173                 if (dt_sigTrk < max_dt_track_cut[i]) {
1174                   N_tracks_MTD[i] = N_tracks_MTD[i] + 1;
1175                   pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();
1176                 }
1177               }
1178             }
1179             // significance cuts
1180             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1181               if (dt_sigTrk_signif < max_dt_significance_cut[i]) {
1182                 N_tracks_MTD_significance[i]++;
1183                 pT_sum_MTD_significance[i] += trackGen.pt();
1184               }
1185             }
1186 
1187           } else {
1188             // if there is no error for MTD information, we count the MTD isolation case same as noMTD
1189             if (optionalPlots_) {
1190               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1191                 N_tracks_MTD[i] = N_tracks_MTD[i] + 1;          // N_tracks_noMTD
1192                 pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();  // pT sum
1193               }
1194             }
1195             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1196               N_tracks_MTD_significance[i]++;
1197               pT_sum_MTD_significance[i] += trackGen.pt();
1198             }
1199             meEle_no_dt_check_->Fill(0);
1200           }
1201 
1202           if (optionalPlots_) {
1203             for (long unsigned int i = 0; i < (pT_bins_dt_distrb.size() - 1); i++) {
1204               //stuff general pT
1205               if (ele.pt() > pT_bins_dt_distrb[i] && ele.pt() < pT_bins_dt_distrb[i + 1]) {
1206                 general_pT_list[i]->Fill(dt_sigTrk);
1207                 general_pT_Signif_list[i]->Fill(dt_sigTrk_signif);
1208               }
1209             }
1210 
1211             for (long unsigned int i = 0; i < (eta_bins_dt_distrib.size() - 1); i++) {
1212               //stuff general eta
1213               if (std::abs(ele.eta()) > eta_bins_dt_distrib[i] && std::abs(ele.eta()) < eta_bins_dt_distrib[i + 1]) {
1214                 general_eta_list[i]->Fill(dt_sigTrk);
1215                 general_eta_Signif_list[i]->Fill(dt_sigTrk_signif);
1216               }
1217             }
1218           }  // End of optional dt distributions plots
1219 
1220           // dt with the vertex
1221         } else {
1222           double dt_vtx = 0;  // dt regular track vs vtx
1223           double dt_vtx_signif = 0;
1224 
1225           double dt_sim_vtx = 0;  // dt regular track vs vtx
1226           double dt_sim_vtx_signif = 0;
1227 
1228           // MTD SIM case
1229           if (std::abs(tsim_trk) > 0 && Vtx_chosen.tError() > 0 && trk_ptSim > 0) {
1230             dt_sim_vtx = std::abs(tsim_trk - Vtx_chosen.t());
1231             dt_sim_vtx_signif = dt_sim_vtx / std::sqrt(avg_sim_PUtrack_t_err * avg_sim_PUtrack_t_err +
1232                                                        Vtx_chosen.tError() * Vtx_chosen.tError());
1233             if (optionalPlots_) {
1234               // absolute timing cuts
1235               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1236                 if (dt_sim_vtx < max_dt_vtx_cut[i]) {
1237                   N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;
1238                   pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;
1239                 }
1240               }
1241             }
1242             // significance timing cuts
1243             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1244               if (dt_sim_vtx_signif < max_dt_significance_cut[i]) {
1245                 N_tracks_sim_MTD_significance[i]++;
1246                 pT_sum_sim_MTD_significance[i] += trk_ptSim;
1247               }
1248             }
1249           } else {
1250             if (optionalPlots_) {
1251               for (long unsigned int i = 0; i < N_tracks_sim_MTD.size(); i++) {
1252                 N_tracks_sim_MTD[i] = N_tracks_sim_MTD[i] + 1;      // N_tracks_noMTD
1253                 pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] + trk_ptSim;  // pT_sum_noMTD
1254               }
1255             }
1256             for (long unsigned int i = 0; i < N_tracks_sim_MTD_significance.size(); i++) {
1257               N_tracks_sim_MTD_significance[i]++;
1258               pT_sum_sim_MTD_significance[i] += trk_ptSim;
1259             }
1260           }
1261 
1262           // MTD RECO case
1263           if (TrkMTDTimeErr > 0 && Vtx_chosen.tError() > 0) {
1264             dt_vtx = std::abs(TrkMTDTime - Vtx_chosen.t());
1265             dt_vtx_signif =
1266                 dt_vtx / std::sqrt(TrkMTDTimeErr * TrkMTDTimeErr + Vtx_chosen.tError() * Vtx_chosen.tError());
1267 
1268             meEle_no_dt_check_->Fill(1);
1269             meEle_avg_error_vtx_check_->Fill(Vtx_chosen.tError());
1270             if (optionalPlots_) {
1271               // absolute timing cuts
1272               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1273                 if (dt_vtx < max_dt_vtx_cut[i]) {
1274                   N_tracks_MTD[i] = N_tracks_MTD[i] + 1;
1275                   pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();
1276                 }
1277               }
1278             }
1279             // significance timing cuts
1280             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1281               if (dt_vtx_signif < max_dt_significance_cut[i]) {
1282                 N_tracks_MTD_significance[i]++;
1283                 pT_sum_MTD_significance[i] += trackGen.pt();
1284               }
1285             }
1286           } else {
1287             if (optionalPlots_) {
1288               for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1289                 N_tracks_MTD[i] = N_tracks_MTD[i] + 1;          // N_tracks_noMTD
1290                 pT_sum_MTD[i] = pT_sum_MTD[i] + trackGen.pt();  // pT_sum_noMTD
1291               }
1292             }
1293             for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1294               N_tracks_MTD_significance[i]++;
1295               pT_sum_MTD_significance[i] += trackGen.pt();
1296             }
1297             meEle_no_dt_check_->Fill(0);
1298           }
1299 
1300           // Optional dt distribution plots
1301           if (optionalPlots_) {
1302             for (long unsigned int i = 0; i < (pT_bins_dt_distrb.size() - 1); i++) {
1303               //stuff general pT
1304               if (ele.pt() > pT_bins_dt_distrb[i] && ele.pt() < pT_bins_dt_distrb[i + 1]) {
1305                 general_pT_list[i]->Fill(dt_vtx);
1306                 general_pT_Signif_list[i]->Fill(dt_vtx_signif);
1307               }
1308             }
1309 
1310             for (long unsigned int i = 0; i < (eta_bins_dt_distrib.size() - 1); i++) {
1311               //stuff general eta
1312               if (std::abs(ele.eta()) > eta_bins_dt_distrib[i] && std::abs(ele.eta()) < eta_bins_dt_distrib[i + 1]) {
1313                 general_eta_list[i]->Fill(dt_vtx);
1314                 general_eta_Signif_list[i]->Fill(dt_vtx_signif);
1315               }
1316             }
1317           }  // End of optional dt distributions plots
1318         }
1319       }
1320       rel_pT_sum_noMTD = pT_sum_noMTD / ele.gsfTrack()->pt();  // rel_ch_iso calculation
1321       if (optionalPlots_) {
1322         for (long unsigned int i = 0; i < N_tracks_MTD.size(); i++) {
1323           rel_pT_sum_MTD[i] = pT_sum_MTD[i] / ele.gsfTrack()->pt();
1324           rel_pT_sum_sim_MTD[i] = pT_sum_sim_MTD[i] / ele_sim_pt;
1325         }
1326         // now compute the isolation
1327         rel_pT_sum_noMTD = pT_sum_noMTD / ele.gsfTrack()->pt();
1328 
1329         rel_pT_sum_gen = pT_sum_gen / ele_sim_pt;
1330       }
1331 
1332       for (long unsigned int i = 0; i < N_tracks_MTD_significance.size(); i++) {
1333         rel_pT_sum_MTD_significance[i] = pT_sum_MTD_significance[i] / ele.gsfTrack()->pt();
1334         rel_pT_sum_sim_MTD_significance[i] = pT_sum_sim_MTD_significance[i] / ele_sim_pt;
1335       }
1336 
1337       if (ele_Promt) {  // promt part
1338         if (Barrel_ele) {
1339           meEleISO_Ntracks_Sig_EB_->Fill(N_tracks_noMTD);
1340           meEleISO_chIso_Sig_EB_->Fill(pT_sum_noMTD);
1341           meEleISO_rel_chIso_Sig_EB_->Fill(rel_pT_sum_noMTD);
1342           if (optionalPlots_) {
1343             for (long unsigned int j = 0; j < Ntracks_EB_list_Sig.size(); j++) {
1344               Ntracks_EB_list_Sig[j]->Fill(N_tracks_MTD[j]);
1345               ch_iso_EB_list_Sig[j]->Fill(pT_sum_MTD[j]);
1346               rel_ch_iso_EB_list_Sig[j]->Fill(rel_pT_sum_MTD[j]);
1347 
1348               Ntracks_sim_EB_list_Sig[j]->Fill(N_tracks_sim_MTD[j]);
1349               ch_iso_sim_EB_list_Sig[j]->Fill(pT_sum_sim_MTD[j]);
1350               rel_ch_iso_sim_EB_list_Sig[j]->Fill(rel_pT_sum_sim_MTD[j]);
1351             }
1352             meEleISO_Ntracks_gen_Sig_EB_->Fill(N_tracks_gen);
1353             meEleISO_chIso_gen_Sig_EB_->Fill(pT_sum_gen);
1354             meEleISO_rel_chIso_gen_Sig_EB_->Fill(rel_pT_sum_gen);
1355           }
1356 
1357           for (long unsigned int j = 0; j < Ntracks_EB_list_Significance_Sig.size(); j++) {
1358             Ntracks_EB_list_Significance_Sig[j]->Fill(N_tracks_MTD_significance[j]);
1359             ch_iso_EB_list_Significance_Sig[j]->Fill(pT_sum_MTD_significance[j]);
1360             rel_ch_iso_EB_list_Significance_Sig[j]->Fill(rel_pT_sum_MTD_significance[j]);
1361 
1362             if (optionalPlots_) {
1363               Ntracks_sim_EB_list_Significance_Sig[j]->Fill(N_tracks_sim_MTD_significance[j]);
1364               ch_iso_sim_EB_list_Significance_Sig[j]->Fill(pT_sum_sim_MTD_significance[j]);
1365               rel_ch_iso_sim_EB_list_Significance_Sig[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1366             }
1367           }
1368 
1369           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1370             meEle_pt_noMTD_Sig_EB_->Fill(ele.pt());
1371             meEle_eta_noMTD_Sig_EB_->Fill(std::abs(ele.eta()));
1372             meEle_phi_noMTD_Sig_EB_->Fill(ele.phi());
1373           }
1374           if (optionalPlots_) {
1375             for (long unsigned int k = 0; k < Ntracks_EB_list_Sig.size(); k++) {
1376               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1377                 Ele_pT_MTD_EB_list_Sig[k]->Fill(ele.pt());
1378                 Ele_eta_MTD_EB_list_Sig[k]->Fill(std::abs(ele.eta()));
1379                 Ele_phi_MTD_EB_list_Sig[k]->Fill(ele.phi());
1380 
1381                 Ele_pT_sim_MTD_EB_list_Sig[k]->Fill(ele_sim_pt);
1382               }
1383             }
1384             if (rel_pT_sum_gen < rel_iso_cut_) {
1385               meEle_pt_gen_Sig_EB_->Fill(ele_sim_pt);
1386               meEle_eta_gen_Sig_EB_->Fill(ele_sim_eta);
1387               meEle_phi_gen_Sig_EB_->Fill(ele_sim_phi);
1388             }
1389           }
1390 
1391           for (long unsigned int k = 0; k < Ntracks_EB_list_Significance_Sig.size(); k++) {
1392             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1393               Ele_pT_MTD_EB_list_Significance_Sig[k]->Fill(ele.pt());
1394               Ele_eta_MTD_EB_list_Significance_Sig[k]->Fill(std::abs(ele.eta()));
1395               Ele_phi_MTD_EB_list_Significance_Sig[k]->Fill(ele.phi());
1396             }
1397             if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1398               Ele_pT_sim_MTD_EB_list_Significance_Sig[k]->Fill(ele_sim_pt);
1399           }
1400 
1401         } else {  // for endcap
1402 
1403           meEleISO_Ntracks_Sig_EE_->Fill(N_tracks_noMTD);
1404           meEleISO_chIso_Sig_EE_->Fill(pT_sum_noMTD);
1405           meEleISO_rel_chIso_Sig_EE_->Fill(rel_pT_sum_noMTD);
1406           if (optionalPlots_) {
1407             for (long unsigned int j = 0; j < Ntracks_EE_list_Sig.size(); j++) {
1408               Ntracks_EE_list_Sig[j]->Fill(N_tracks_MTD[j]);
1409               ch_iso_EE_list_Sig[j]->Fill(pT_sum_MTD[j]);
1410               rel_ch_iso_EE_list_Sig[j]->Fill(rel_pT_sum_MTD[j]);
1411 
1412               Ntracks_sim_EE_list_Sig[j]->Fill(N_tracks_sim_MTD[j]);
1413               ch_iso_sim_EE_list_Sig[j]->Fill(pT_sum_sim_MTD[j]);
1414               rel_ch_iso_sim_EE_list_Sig[j]->Fill(rel_pT_sum_sim_MTD[j]);
1415             }
1416             meEleISO_Ntracks_gen_Sig_EE_->Fill(N_tracks_gen);
1417             meEleISO_chIso_gen_Sig_EE_->Fill(pT_sum_gen);
1418             meEleISO_rel_chIso_gen_Sig_EE_->Fill(rel_pT_sum_gen);
1419           }
1420 
1421           for (long unsigned int j = 0; j < Ntracks_EE_list_Significance_Sig.size(); j++) {
1422             Ntracks_EE_list_Significance_Sig[j]->Fill(N_tracks_MTD_significance[j]);
1423             ch_iso_EE_list_Significance_Sig[j]->Fill(pT_sum_MTD_significance[j]);
1424             rel_ch_iso_EE_list_Significance_Sig[j]->Fill(rel_pT_sum_MTD_significance[j]);
1425 
1426             if (optionalPlots_) {
1427               Ntracks_sim_EE_list_Significance_Sig[j]->Fill(N_tracks_sim_MTD_significance[j]);
1428               ch_iso_sim_EE_list_Significance_Sig[j]->Fill(pT_sum_sim_MTD_significance[j]);
1429               rel_ch_iso_sim_EE_list_Significance_Sig[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1430             }
1431           }
1432 
1433           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1434             meEle_pt_noMTD_Sig_EE_->Fill(ele.pt());
1435             meEle_eta_noMTD_Sig_EE_->Fill(std::abs(ele.eta()));
1436             meEle_phi_noMTD_Sig_EE_->Fill(ele.phi());
1437           }
1438           if (optionalPlots_) {
1439             for (long unsigned int k = 0; k < Ntracks_EE_list_Sig.size(); k++) {
1440               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1441                 Ele_pT_MTD_EE_list_Sig[k]->Fill(ele.pt());
1442                 Ele_eta_MTD_EE_list_Sig[k]->Fill(std::abs(ele.eta()));
1443                 Ele_phi_MTD_EE_list_Sig[k]->Fill(ele.phi());
1444 
1445                 Ele_pT_sim_MTD_EE_list_Sig[k]->Fill(ele_sim_pt);
1446               }
1447             }
1448             if (rel_pT_sum_gen < rel_iso_cut_) {
1449               meEle_pt_gen_Sig_EE_->Fill(ele_sim_pt);
1450               meEle_eta_gen_Sig_EE_->Fill(ele_sim_eta);
1451               meEle_phi_gen_Sig_EE_->Fill(ele_sim_phi);
1452             }
1453           }
1454           for (long unsigned int k = 0; k < Ntracks_EE_list_Significance_Sig.size(); k++) {
1455             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1456               Ele_pT_MTD_EE_list_Significance_Sig[k]->Fill(ele.pt());
1457               Ele_eta_MTD_EE_list_Significance_Sig[k]->Fill(std::abs(ele.eta()));
1458               Ele_phi_MTD_EE_list_Significance_Sig[k]->Fill(ele.phi());
1459 
1460               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1461                 Ele_pT_sim_MTD_EE_list_Significance_Sig[k]->Fill(ele_sim_pt);
1462             }
1463           }
1464         }
1465       } else {  // non-promt part
1466         if (Barrel_ele) {
1467           meEleISO_Ntracks_Bkg_EB_->Fill(N_tracks_noMTD);
1468           meEleISO_chIso_Bkg_EB_->Fill(pT_sum_noMTD);
1469           meEleISO_rel_chIso_Bkg_EB_->Fill(rel_pT_sum_noMTD);
1470           if (optionalPlots_) {
1471             for (long unsigned int j = 0; j < Ntracks_EB_list_Bkg.size(); j++) {
1472               Ntracks_EB_list_Bkg[j]->Fill(N_tracks_MTD[j]);
1473               ch_iso_EB_list_Bkg[j]->Fill(pT_sum_MTD[j]);
1474               rel_ch_iso_EB_list_Bkg[j]->Fill(rel_pT_sum_MTD[j]);
1475 
1476               Ntracks_sim_EB_list_Bkg[j]->Fill(N_tracks_sim_MTD[j]);
1477               ch_iso_sim_EB_list_Bkg[j]->Fill(pT_sum_sim_MTD[j]);
1478               rel_ch_iso_sim_EB_list_Bkg[j]->Fill(rel_pT_sum_sim_MTD[j]);
1479             }
1480             meEleISO_Ntracks_gen_Bkg_EB_->Fill(N_tracks_gen);
1481             meEleISO_chIso_gen_Bkg_EB_->Fill(pT_sum_gen);
1482             meEleISO_rel_chIso_gen_Bkg_EB_->Fill(rel_pT_sum_gen);
1483           }
1484 
1485           for (long unsigned int j = 0; j < Ntracks_EB_list_Significance_Bkg.size(); j++) {
1486             Ntracks_EB_list_Significance_Bkg[j]->Fill(N_tracks_MTD_significance[j]);
1487             ch_iso_EB_list_Significance_Bkg[j]->Fill(pT_sum_MTD_significance[j]);
1488             rel_ch_iso_EB_list_Significance_Bkg[j]->Fill(rel_pT_sum_MTD_significance[j]);
1489 
1490             if (optionalPlots_) {
1491               Ntracks_sim_EB_list_Significance_Bkg[j]->Fill(N_tracks_sim_MTD_significance[j]);
1492               ch_iso_sim_EB_list_Significance_Bkg[j]->Fill(pT_sum_sim_MTD_significance[j]);
1493               rel_ch_iso_sim_EB_list_Significance_Bkg[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1494             }
1495           }
1496 
1497           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1498             meEle_pt_noMTD_Bkg_EB_->Fill(ele.pt());
1499             meEle_eta_noMTD_Bkg_EB_->Fill(std::abs(ele.eta()));
1500             meEle_phi_noMTD_Bkg_EB_->Fill(ele.phi());
1501           }
1502           if (optionalPlots_) {
1503             for (long unsigned int k = 0; k < Ntracks_EB_list_Bkg.size(); k++) {
1504               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1505                 Ele_pT_MTD_EB_list_Bkg[k]->Fill(ele.pt());
1506                 Ele_eta_MTD_EB_list_Bkg[k]->Fill(std::abs(ele.eta()));
1507                 Ele_phi_MTD_EB_list_Bkg[k]->Fill(ele.phi());
1508 
1509                 Ele_pT_sim_MTD_EB_list_Bkg[k]->Fill(ele_sim_pt);
1510               }
1511             }
1512             if (rel_pT_sum_gen < rel_iso_cut_) {
1513               meEle_pt_gen_Bkg_EB_->Fill(ele_sim_pt);
1514               meEle_eta_gen_Bkg_EB_->Fill(ele_sim_eta);
1515               meEle_phi_gen_Bkg_EB_->Fill(ele_sim_phi);
1516             }
1517           }
1518           for (long unsigned int k = 0; k < Ntracks_EB_list_Significance_Bkg.size(); k++) {
1519             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1520               Ele_pT_MTD_EB_list_Significance_Bkg[k]->Fill(ele.pt());
1521               Ele_eta_MTD_EB_list_Significance_Bkg[k]->Fill(std::abs(ele.eta()));
1522               Ele_phi_MTD_EB_list_Significance_Bkg[k]->Fill(ele.phi());
1523 
1524               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1525                 Ele_pT_sim_MTD_EB_list_Significance_Bkg[k]->Fill(ele_sim_pt);
1526             }
1527           }
1528 
1529         } else {  // for endcap
1530           meEleISO_Ntracks_Bkg_EE_->Fill(N_tracks_noMTD);
1531           meEleISO_chIso_Bkg_EE_->Fill(pT_sum_noMTD);
1532           meEleISO_rel_chIso_Bkg_EE_->Fill(rel_pT_sum_noMTD);
1533           if (optionalPlots_) {
1534             for (long unsigned int j = 0; j < Ntracks_EE_list_Bkg.size(); j++) {
1535               Ntracks_EE_list_Bkg[j]->Fill(N_tracks_MTD[j]);
1536               ch_iso_EE_list_Bkg[j]->Fill(pT_sum_MTD[j]);
1537               rel_ch_iso_EE_list_Bkg[j]->Fill(rel_pT_sum_MTD[j]);
1538 
1539               Ntracks_sim_EE_list_Bkg[j]->Fill(N_tracks_sim_MTD[j]);
1540               ch_iso_sim_EE_list_Bkg[j]->Fill(pT_sum_sim_MTD[j]);
1541               rel_ch_iso_sim_EE_list_Bkg[j]->Fill(rel_pT_sum_sim_MTD[j]);
1542             }
1543             meEleISO_Ntracks_gen_Bkg_EE_->Fill(N_tracks_gen);
1544             meEleISO_chIso_gen_Bkg_EE_->Fill(pT_sum_gen);
1545             meEleISO_rel_chIso_gen_Bkg_EE_->Fill(rel_pT_sum_gen);
1546           }
1547 
1548           for (long unsigned int j = 0; j < Ntracks_EE_list_Significance_Bkg.size(); j++) {
1549             Ntracks_EE_list_Significance_Bkg[j]->Fill(N_tracks_MTD_significance[j]);
1550             ch_iso_EE_list_Significance_Bkg[j]->Fill(pT_sum_MTD_significance[j]);
1551             rel_ch_iso_EE_list_Significance_Bkg[j]->Fill(rel_pT_sum_MTD_significance[j]);
1552 
1553             if (optionalPlots_) {
1554               Ntracks_sim_EE_list_Significance_Bkg[j]->Fill(N_tracks_sim_MTD_significance[j]);
1555               ch_iso_sim_EE_list_Significance_Bkg[j]->Fill(pT_sum_sim_MTD_significance[j]);
1556               rel_ch_iso_sim_EE_list_Significance_Bkg[j]->Fill(rel_pT_sum_sim_MTD_significance[j]);
1557             }
1558           }
1559 
1560           if (rel_pT_sum_noMTD < rel_iso_cut_) {  // filling hists for iso efficiency calculations
1561             meEle_pt_noMTD_Bkg_EE_->Fill(ele.pt());
1562             meEle_eta_noMTD_Bkg_EE_->Fill(std::abs(ele.eta()));
1563             meEle_phi_noMTD_Bkg_EE_->Fill(ele.phi());
1564           }
1565           if (optionalPlots_) {
1566             for (long unsigned int k = 0; k < Ntracks_EE_list_Bkg.size(); k++) {
1567               if (rel_pT_sum_MTD[k] < rel_iso_cut_) {
1568                 Ele_pT_MTD_EE_list_Bkg[k]->Fill(ele.pt());
1569                 Ele_eta_MTD_EE_list_Bkg[k]->Fill(std::abs(ele.eta()));
1570                 Ele_phi_MTD_EE_list_Bkg[k]->Fill(ele.phi());
1571 
1572                 Ele_pT_sim_MTD_EE_list_Bkg[k]->Fill(ele_sim_pt);
1573               }
1574             }
1575             if (rel_pT_sum_gen < rel_iso_cut_) {
1576               meEle_pt_gen_Bkg_EE_->Fill(ele_sim_pt);
1577               meEle_eta_gen_Bkg_EE_->Fill(ele_sim_eta);
1578               meEle_phi_gen_Bkg_EE_->Fill(ele_sim_phi);
1579             }
1580           }
1581 
1582           for (long unsigned int k = 0; k < Ntracks_EE_list_Significance_Bkg.size(); k++) {
1583             if (rel_pT_sum_MTD_significance[k] < rel_iso_cut_) {
1584               Ele_pT_MTD_EE_list_Significance_Bkg[k]->Fill(ele.pt());
1585               Ele_eta_MTD_EE_list_Significance_Bkg[k]->Fill(std::abs(ele.eta()));
1586               Ele_phi_MTD_EE_list_Significance_Bkg[k]->Fill(ele.phi());
1587 
1588               if (optionalPlots_ and rel_pT_sum_sim_MTD_significance[k] < rel_iso_cut_)
1589                 Ele_pT_sim_MTD_EE_list_Significance_Bkg[k]->Fill(ele_sim_pt);
1590             }
1591           }
1592         }
1593       }
1594     }  // electron matched to a track
1595   }    // electron collection inside single event
1596 }
1597 
1598 // ------------ method for histogram booking ------------
1599 void MtdEleIsoValidation::bookHistograms(DQMStore::IBooker& ibook, edm::Run const& run, edm::EventSetup const& iSetup) {
1600   ibook.setCurrentFolder(folder_);
1601 
1602   // for regular Validation use a reduced binning, for detailed analysis and ROC curves use the larger one
1603   int nbin_1 = 40;
1604   int nbin_2 = 40;
1605   if (optionalPlots_) {
1606     nbin_1 = 1000;
1607     nbin_2 = 2000;
1608   }
1609 
1610   // histogram booking
1611 
1612   meEle_avg_error_SigTrk_check_ =
1613       ibook.book1D("SigTrk_avg_timing_err",
1614                    "Average signal electron track MTD timing uncertainty;Time Error (ns);Counts",
1615                    200,
1616                    0,
1617                    0.1);
1618   meEle_avg_error_PUTrk_check_ = ibook.book1D(
1619       "PUTrk_avg_timing_err", "Average PU track MTD timing uncertainty;Time Error (ns);Counts", 200, 0, 0.1);
1620   meEle_avg_error_vtx_check_ =
1621       ibook.book1D("Vtx_avg_timing_err", "Average vertex timing uncertainty;Time Error (ns);Counts", 200, 0, 0.1);
1622 
1623   meEle_no_dt_check_ =
1624       ibook.book1D("Track_dt_info_check",
1625                    "Tracks dt check - ratio between tracks with (value 1) and without (value 0) timing info",
1626                    2,
1627                    0,
1628                    2);
1629 
1630   meTrk_genMatch_check_ = ibook.book1D(
1631       "Track_genMatch_info_check", "Check on tracks matched with a GenParticle (matched=1, non matched=0)", 2, 0, 2);
1632 
1633   // signal
1634   meEleISO_Ntracks_Sig_EB_ = ibook.book1D("Ele_Iso_Ntracks_Sig_EB",
1635                                           "Number of tracks in isolation cone around electron track after basic cuts - "
1636                                           "Signal Barrel;Number of tracks;Counts",
1637                                           20,
1638                                           0,
1639                                           20);
1640 
1641   meEleISO_chIso_Sig_EB_ = ibook.book1D(
1642       "Ele_chIso_sum_Sig_EB",
1643       "Track pT sum in isolation cone around electron track after basic cuts - Signal Barrel;p_{T} (GeV);Counts",
1644       nbin_2,
1645       0,
1646       20);
1647 
1648   meEleISO_rel_chIso_Sig_EB_ = ibook.book1D(
1649       "Ele_rel_chIso_sum_Sig_EB",
1650       "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Barrel;Isolation;Counts",
1651       nbin_1,
1652       0,
1653       4);
1654   if (optionalPlots_) {
1655     meEleISO_Ntracks_MTD_1_Sig_EB_ =
1656         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Sig_EB",
1657                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1658                      20,
1659                      0,
1660                      20);
1661 
1662     meEleISO_chIso_MTD_1_Sig_EB_ = ibook.book1D(
1663         "Ele_chIso_sum_MTD_1_Sig_EB",
1664         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1665         nbin_2,
1666         0,
1667         20);
1668     meEleISO_rel_chIso_MTD_1_Sig_EB_ = ibook.book1D(
1669         "Ele_rel_chIso_sum_MTD_1_Sig_EB",
1670         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1671         nbin_1,
1672         0,
1673         4);
1674     // gen
1675     meEleISO_Ntracks_gen_Sig_EB_ = ibook.book1D("Ele_Iso_Ntracks_gen_Sig_EB",
1676                                                 "Number of tracks in isolation cone around electron track after basic "
1677                                                 "cuts using genInfo - Signal Barrel;Number of tracks;Counts",
1678                                                 20,
1679                                                 0,
1680                                                 20);
1681 
1682     meEleISO_chIso_gen_Sig_EB_ = ibook.book1D("Ele_chIso_sum_gen_Sig_EB",
1683                                               "Track pT sum in isolation cone around electron track after basic cuts "
1684                                               "using genInfo - Signal Barrel;p_{T} (GeV);Counts",
1685                                               nbin_2,
1686                                               0,
1687                                               20);
1688 
1689     meEleISO_rel_chIso_gen_Sig_EB_ = ibook.book1D("Ele_rel_chIso_sum_gen_Sig_EB",
1690                                                   "Track relative pT sum in isolation cone around electron track after "
1691                                                   "basic cuts using genInfo - Signal Barrel;Isolation;Counts",
1692                                                   nbin_1,
1693                                                   0,
1694                                                   4);
1695 
1696     meEleISO_Ntracks_MTD_2_Sig_EB_ =
1697         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Sig_EB",
1698                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1699                      20,
1700                      0,
1701                      20);
1702 
1703     meEleISO_chIso_MTD_2_Sig_EB_ = ibook.book1D(
1704         "Ele_chIso_sum_MTD_2_Sig_EB",
1705         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1706         nbin_2,
1707         0,
1708         20);
1709     meEleISO_rel_chIso_MTD_2_Sig_EB_ = ibook.book1D(
1710         "Ele_rel_chIso_sum_MTD_2_Sig_EB",
1711         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1712         nbin_1,
1713         0,
1714         4);
1715 
1716     meEleISO_Ntracks_MTD_3_Sig_EB_ =
1717         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Sig_EB",
1718                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1719                      20,
1720                      0,
1721                      20);
1722     meEleISO_chIso_MTD_3_Sig_EB_ = ibook.book1D(
1723         "Ele_chIso_sum_MTD_3_Sig_EB",
1724         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1725         nbin_2,
1726         0,
1727         20);
1728     meEleISO_rel_chIso_MTD_3_Sig_EB_ = ibook.book1D(
1729         "Ele_rel_chIso_sum_MTD_3_Sig_EB",
1730         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1731         nbin_1,
1732         0,
1733         4);
1734 
1735     meEleISO_Ntracks_MTD_4_Sig_EB_ =
1736         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Sig_EB",
1737                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1738                      20,
1739                      0,
1740                      20);
1741     meEleISO_chIso_MTD_4_Sig_EB_ = ibook.book1D(
1742         "Ele_chIso_sum_MTD_4_Sig_EB",
1743         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1744         nbin_2,
1745         0,
1746         20);
1747     meEleISO_rel_chIso_MTD_4_Sig_EB_ = ibook.book1D(
1748         "Ele_rel_chIso_sum_MTD_4_Sig_EB",
1749         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1750         nbin_1,
1751         0,
1752         4);
1753 
1754     meEleISO_Ntracks_MTD_5_Sig_EB_ =
1755         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Sig_EB",
1756                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1757                      20,
1758                      0,
1759                      20);
1760     meEleISO_chIso_MTD_5_Sig_EB_ = ibook.book1D(
1761         "Ele_chIso_sum_MTD_5_Sig_EB",
1762         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1763         nbin_2,
1764         0,
1765         20);
1766     meEleISO_rel_chIso_MTD_5_Sig_EB_ = ibook.book1D(
1767         "Ele_rel_chIso_sum_MTD_5_Sig_EB",
1768         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1769         nbin_1,
1770         0,
1771         4);
1772 
1773     meEleISO_Ntracks_MTD_6_Sig_EB_ =
1774         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Sig_EB",
1775                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1776                      20,
1777                      0,
1778                      20);
1779     meEleISO_chIso_MTD_6_Sig_EB_ = ibook.book1D(
1780         "Ele_chIso_sum_MTD_6_Sig_EB",
1781         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1782         nbin_2,
1783         0,
1784         20);
1785     meEleISO_rel_chIso_MTD_6_Sig_EB_ = ibook.book1D(
1786         "Ele_rel_chIso_sum_MTD_6_Sig_EB",
1787         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1788         nbin_1,
1789         0,
1790         4);
1791 
1792     meEleISO_Ntracks_MTD_7_Sig_EB_ =
1793         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Sig_EB",
1794                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1795                      20,
1796                      0,
1797                      20);
1798     meEleISO_chIso_MTD_7_Sig_EB_ = ibook.book1D(
1799         "Ele_chIso_sum_MTD_7_Sig_EB",
1800         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1801         nbin_2,
1802         0,
1803         20);
1804     meEleISO_rel_chIso_MTD_7_Sig_EB_ = ibook.book1D(
1805         "Ele_rel_chIso_sum_MTD_7_Sig_EB",
1806         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1807         nbin_1,
1808         0,
1809         4);
1810 
1811     meEleISO_Ntracks_MTD_sim_1_Sig_EB_ =
1812         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Sig_EB",
1813                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1814                      20,
1815                      0,
1816                      20);
1817 
1818     meEleISO_chIso_MTD_sim_1_Sig_EB_ = ibook.book1D(
1819         "Ele_chIso_sum_MTD_sim_1_Sig_EB",
1820         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1821         nbin_2,
1822         0,
1823         20);
1824     meEleISO_rel_chIso_MTD_sim_1_Sig_EB_ = ibook.book1D(
1825         "Ele_rel_chIso_sum_MTD_sim_1_Sig_EB",
1826         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1827         nbin_1,
1828         0,
1829         4);
1830 
1831     meEleISO_Ntracks_MTD_sim_2_Sig_EB_ =
1832         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Sig_EB",
1833                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1834                      20,
1835                      0,
1836                      20);
1837 
1838     meEleISO_chIso_MTD_sim_2_Sig_EB_ = ibook.book1D(
1839         "Ele_chIso_sum_MTD_sim_2_Sig_EB",
1840         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1841         nbin_2,
1842         0,
1843         20);
1844     meEleISO_rel_chIso_MTD_sim_2_Sig_EB_ = ibook.book1D(
1845         "Ele_rel_chIso_sum_MTD_sim_2_Sig_EB",
1846         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1847         nbin_1,
1848         0,
1849         4);
1850 
1851     meEleISO_Ntracks_MTD_sim_3_Sig_EB_ =
1852         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Sig_EB",
1853                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1854                      20,
1855                      0,
1856                      20);
1857     meEleISO_chIso_MTD_sim_3_Sig_EB_ = ibook.book1D(
1858         "Ele_chIso_sum_MTD_sim_3_Sig_EB",
1859         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1860         nbin_2,
1861         0,
1862         20);
1863     meEleISO_rel_chIso_MTD_sim_3_Sig_EB_ = ibook.book1D(
1864         "Ele_rel_chIso_sum_MTD_sim_3_Sig_EB",
1865         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1866         nbin_1,
1867         0,
1868         4);
1869 
1870     meEleISO_Ntracks_MTD_sim_4_Sig_EB_ =
1871         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Sig_EB",
1872                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1873                      20,
1874                      0,
1875                      20);
1876     meEleISO_chIso_MTD_sim_4_Sig_EB_ = ibook.book1D(
1877         "Ele_chIso_sum_MTD_sim_4_Sig_EB",
1878         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1879         nbin_2,
1880         0,
1881         20);
1882     meEleISO_rel_chIso_MTD_sim_4_Sig_EB_ = ibook.book1D(
1883         "Ele_rel_chIso_sum_MTD_sim_4_Sig_EB",
1884         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1885         nbin_1,
1886         0,
1887         4);
1888 
1889     meEleISO_Ntracks_MTD_sim_5_Sig_EB_ =
1890         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Sig_EB",
1891                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1892                      20,
1893                      0,
1894                      20);
1895     meEleISO_chIso_MTD_sim_5_Sig_EB_ = ibook.book1D(
1896         "Ele_chIso_sum_MTD_sim_5_Sig_EB",
1897         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1898         nbin_2,
1899         0,
1900         20);
1901     meEleISO_rel_chIso_MTD_sim_5_Sig_EB_ = ibook.book1D(
1902         "Ele_rel_chIso_sum_MTD_sim_5_Sig_EB",
1903         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1904         nbin_1,
1905         0,
1906         4);
1907 
1908     meEleISO_Ntracks_MTD_sim_6_Sig_EB_ =
1909         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Sig_EB",
1910                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1911                      20,
1912                      0,
1913                      20);
1914     meEleISO_chIso_MTD_sim_6_Sig_EB_ = ibook.book1D(
1915         "Ele_chIso_sum_MTD_sim_6_Sig_EB",
1916         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1917         nbin_2,
1918         0,
1919         20);
1920     meEleISO_rel_chIso_MTD_sim_6_Sig_EB_ = ibook.book1D(
1921         "Ele_rel_chIso_sum_MTD_sim_6_Sig_EB",
1922         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1923         nbin_1,
1924         0,
1925         4);
1926 
1927     meEleISO_Ntracks_MTD_sim_7_Sig_EB_ =
1928         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Sig_EB",
1929                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
1930                      20,
1931                      0,
1932                      20);
1933     meEleISO_chIso_MTD_sim_7_Sig_EB_ = ibook.book1D(
1934         "Ele_chIso_sum_MTD_sim_7_Sig_EB",
1935         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
1936         nbin_2,
1937         0,
1938         20);
1939     meEleISO_rel_chIso_MTD_sim_7_Sig_EB_ = ibook.book1D(
1940         "Ele_rel_chIso_sum_MTD_sim_7_Sig_EB",
1941         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
1942         nbin_1,
1943         0,
1944         4);
1945   }
1946   meEleISO_Ntracks_MTD_4sigma_Sig_EB_ =
1947       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Sig_EB",
1948                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma compatibiliy - "
1949                    "Signal Barrel;Number of tracks;Counts",
1950                    20,
1951                    0,
1952                    20);
1953   meEleISO_chIso_MTD_4sigma_Sig_EB_ =
1954       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Sig_EB",
1955                    "Track pT sum in isolation cone around electron track after basic "
1956                    "cuts with MTD - 4 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
1957                    nbin_2,
1958                    0,
1959                    20);
1960   meEleISO_rel_chIso_MTD_4sigma_Sig_EB_ =
1961       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Sig_EB",
1962                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
1963                    "compatibiliy - Signal Barrel;Isolation;Counts",
1964                    nbin_1,
1965                    0,
1966                    4);
1967 
1968   meEleISO_Ntracks_MTD_3sigma_Sig_EB_ =
1969       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Sig_EB",
1970                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma compatibiliy - "
1971                    "Signal Barrel;Number of tracks;Counts",
1972                    20,
1973                    0,
1974                    20);
1975   meEleISO_chIso_MTD_3sigma_Sig_EB_ =
1976       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Sig_EB",
1977                    "Track pT sum in isolation cone around electron track after basic "
1978                    "cuts with MTD - 3 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
1979                    nbin_2,
1980                    0,
1981                    20);
1982   meEleISO_rel_chIso_MTD_3sigma_Sig_EB_ = ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Sig_EB",
1983                                                        "Track relative pT sum in isolation cone around electron track "
1984                                                        "after basic cuts with MTD - 3 sigma;Isolation;Counts"
1985                                                        "compatibiliy - Signal Barrel",
1986                                                        nbin_1,
1987                                                        0,
1988                                                        4);
1989 
1990   meEleISO_Ntracks_MTD_2sigma_Sig_EB_ =
1991       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Sig_EB",
1992                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibiliy - "
1993                    "Signal Barrel;Number of tracks;Counts",
1994                    20,
1995                    0,
1996                    20);
1997   meEleISO_chIso_MTD_2sigma_Sig_EB_ =
1998       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Sig_EB",
1999                    "Track pT sum in isolation cone around electron track after basic "
2000                    "cuts with MTD - 2 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2001                    nbin_2,
2002                    0,
2003                    20);
2004   meEleISO_rel_chIso_MTD_2sigma_Sig_EB_ =
2005       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Sig_EB",
2006                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2007                    "compatibiliy - Signal Barrel;Isolation;Counts",
2008                    nbin_1,
2009                    0,
2010                    4);
2011 
2012   meEle_pt_tot_Sig_EB_ =
2013       ibook.book1D("Ele_pT_tot_Sig_EB", "Electron pT tot - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2014   meEle_pt_noMTD_Sig_EB_ =
2015       ibook.book1D("Ele_pT_noMTD_Sig_EB", "Electron pT noMTD - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2016 
2017   meEle_pt_sim_tot_Sig_EB_ =
2018       ibook.book1D("Ele_pT_sim_tot_Sig_EB", "Electron SIM pT tot - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2019 
2020   meEle_eta_tot_Sig_EB_ =
2021       ibook.book1D("Ele_eta_tot_Sig_EB", "Electron eta tot - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2022   meEle_eta_noMTD_Sig_EB_ =
2023       ibook.book1D("Ele_eta_noMTD_Sig_EB", "Electron eta noMTD - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2024 
2025   meEle_phi_tot_Sig_EB_ =
2026       ibook.book1D("Ele_phi_tot_Sig_EB", "Electron phi tot - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2027   meEle_phi_noMTD_Sig_EB_ =
2028       ibook.book1D("Ele_phi_noMTD_Sig_EB", "Electron phi noMTD - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2029 
2030   if (optionalPlots_) {
2031     meEleISO_Ntracks_MTD_sim_4sigma_Sig_EB_ =
2032         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Sig_EB",
2033                      "Number of tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2034                      "compatibiliy - Signal Barrel;Number of tracks;Counts",
2035                      20,
2036                      0,
2037                      20);
2038     meEleISO_chIso_MTD_sim_4sigma_Sig_EB_ =
2039         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Sig_EB",
2040                      "Track pT sum in isolation cone around electron track after "
2041                      "basic cuts with MTD - 4 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2042                      nbin_2,
2043                      0,
2044                      20);
2045     meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EB_ = ibook.book1D(
2046         "Ele_rel_chIso_sum_MTD_sim_4sigma_Sig_EB",
2047         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2048         "compatibiliy - Signal Barrel;Isolation;Counts",
2049         nbin_1,
2050         0,
2051         4);
2052 
2053     meEleISO_Ntracks_MTD_sim_3sigma_Sig_EB_ =
2054         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Sig_EB",
2055                      "Tracks in isolation cone around electron track after basic "
2056                      "cuts with MTD  - 3 sigma compatibiliy - Signal Barrel;Number of tracks;Counts",
2057                      20,
2058                      0,
2059                      20);
2060     meEleISO_chIso_MTD_sim_3sigma_Sig_EB_ =
2061         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Sig_EB",
2062                      "Track pT sum in isolation cone around electron track after "
2063                      "basic cuts with MTD - 3 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2064                      nbin_2,
2065                      0,
2066                      20);
2067     meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EB_ = ibook.book1D(
2068         "Ele_rel_chIso_sum_MTD_sim_3sigma_Sig_EB",
2069         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 3 sigma "
2070         "compatibiliy - Signal Barrel;Isolation;Counts",
2071         nbin_1,
2072         0,
2073         4);
2074 
2075     meEleISO_Ntracks_MTD_sim_2sigma_Sig_EB_ = ibook.book1D(
2076         "Ele_Iso_Ntracks_MTD_sim_2sigma_Sig_EB",
2077         "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibiliy - "
2078         "Signal Barrel;Number of tracks;Counts",
2079         20,
2080         0,
2081         20);
2082     meEleISO_chIso_MTD_sim_2sigma_Sig_EB_ =
2083         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Sig_EB",
2084                      "Track pT sum in isolation cone around electron track after "
2085                      "basic cuts with MTD - 2 sigma compatibiliy - Signal Barrel;p_{T} (GeV);Counts",
2086                      nbin_2,
2087                      0,
2088                      20);
2089     meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EB_ = ibook.book1D(
2090         "Ele_rel_chIso_sum_MTD_sim_2sigma_Sig_EB",
2091         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2092         "compatibiliy - Signal Barrel;Isolation;Counts",
2093         nbin_1,
2094         0,
2095         4);
2096 
2097     meEle_pt_gen_Sig_EB_ =
2098         ibook.book1D("Ele_pT_gen_Sig_EB", "Electron pT genInfo - Signal Barrel;p_{T} (GeV);Counts", 30, 10, 100);
2099     meEle_eta_gen_Sig_EB_ =
2100         ibook.book1D("Ele_eta_gen_Sig_EB", "Electron eta genInfo - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2101     meEle_phi_gen_Sig_EB_ =
2102         ibook.book1D("Ele_phi_gen_Sig_EB", "Electron phi genInfo - Signal Barrel;#phi;Counts", 64, -3.2, 3.2);
2103 
2104     meEle_pt_MTD_1_Sig_EB_ = ibook.book1D("Ele_pT_MTD_1_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2105     meEle_eta_MTD_1_Sig_EB_ = ibook.book1D("Ele_eta_MTD_1_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2106     meEle_phi_MTD_1_Sig_EB_ = ibook.book1D("Ele_phi_MTD_1_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2107 
2108     meEle_pt_MTD_2_Sig_EB_ = ibook.book1D("Ele_pT_MTD_2_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2109     meEle_eta_MTD_2_Sig_EB_ = ibook.book1D("Ele_eta_MTD_2_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2110     meEle_phi_MTD_2_Sig_EB_ = ibook.book1D("Ele_phi_MTD_2_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2111 
2112     meEle_pt_MTD_3_Sig_EB_ = ibook.book1D("Ele_pT_MTD_3_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2113     meEle_eta_MTD_3_Sig_EB_ = ibook.book1D("Ele_eta_MTD_3_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2114     meEle_phi_MTD_3_Sig_EB_ = ibook.book1D("Ele_phi_MTD_3_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2115 
2116     meEle_pt_MTD_4_Sig_EB_ = ibook.book1D("Ele_pT_MTD_4_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2117     meEle_eta_MTD_4_Sig_EB_ = ibook.book1D("Ele_eta_MTD_4_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2118     meEle_phi_MTD_4_Sig_EB_ = ibook.book1D("Ele_phi_MTD_4_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2119 
2120     meEle_pt_MTD_5_Sig_EB_ = ibook.book1D("Ele_pT_MTD_5_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2121     meEle_eta_MTD_5_Sig_EB_ = ibook.book1D("Ele_eta_MTD_5_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2122     meEle_phi_MTD_5_Sig_EB_ = ibook.book1D("Ele_phi_MTD_5_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2123 
2124     meEle_pt_MTD_6_Sig_EB_ = ibook.book1D("Ele_pT_MTD_6_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2125     meEle_eta_MTD_6_Sig_EB_ = ibook.book1D("Ele_eta_MTD_6_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2126     meEle_phi_MTD_6_Sig_EB_ = ibook.book1D("Ele_phi_MTD_6_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2127 
2128     meEle_pt_MTD_7_Sig_EB_ = ibook.book1D("Ele_pT_MTD_7_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2129     meEle_eta_MTD_7_Sig_EB_ = ibook.book1D("Ele_eta_MTD_7_Sig_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
2130     meEle_phi_MTD_7_Sig_EB_ = ibook.book1D("Ele_phi_MTD_7_Sig_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2131 
2132     meEle_pt_sim_MTD_1_Sig_EB_ =
2133         ibook.book1D("Ele_pT_sim_MTD_1_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2134     meEle_pt_sim_MTD_2_Sig_EB_ =
2135         ibook.book1D("Ele_pT_sim_MTD_2_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2136     meEle_pt_sim_MTD_3_Sig_EB_ =
2137         ibook.book1D("Ele_pT_sim_MTD_3_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2138     meEle_pt_sim_MTD_4_Sig_EB_ =
2139         ibook.book1D("Ele_pT_sim_MTD_4_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2140     meEle_pt_sim_MTD_5_Sig_EB_ =
2141         ibook.book1D("Ele_pT_sim_MTD_5_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2142     meEle_pt_sim_MTD_6_Sig_EB_ =
2143         ibook.book1D("Ele_pT_sim_MTD_6_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2144     meEle_pt_sim_MTD_7_Sig_EB_ =
2145         ibook.book1D("Ele_pT_sim_MTD_7_Sig_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2146   }
2147 
2148   meEle_pt_MTD_4sigma_Sig_EB_ =
2149       ibook.book1D("Ele_pT_MTD_4sigma_Sig_EB",
2150                    "Electron pT MTD - 4 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2151                    30,
2152                    10,
2153                    100);
2154   meEle_eta_MTD_4sigma_Sig_EB_ = ibook.book1D(
2155       "Ele_eta_MTD_4sigma_Sig_EB", "Electron eta MTD - 4 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2156   meEle_phi_MTD_4sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_4sigma_Sig_EB",
2157                                               "Electron phi MTD - 4 sigma compatibility - Signal Barrel;#phi;Counts",
2158                                               64,
2159                                               -3.2,
2160                                               3.2);
2161 
2162   meEle_pt_MTD_3sigma_Sig_EB_ =
2163       ibook.book1D("Ele_pT_MTD_3sigma_Sig_EB",
2164                    "Electron pT MTD - 3 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2165                    30,
2166                    10,
2167                    100);
2168   meEle_eta_MTD_3sigma_Sig_EB_ = ibook.book1D(
2169       "Ele_eta_MTD_3sigma_Sig_EB", "Electron eta MTD - 3 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2170   meEle_phi_MTD_3sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_3sigma_Sig_EB",
2171                                               "Electron phi MTD - 3 sigma compatibility - Signal Barrel;#phi;Counts",
2172                                               64,
2173                                               -3.2,
2174                                               3.2);
2175 
2176   meEle_pt_MTD_2sigma_Sig_EB_ =
2177       ibook.book1D("Ele_pT_MTD_2sigma_Sig_EB",
2178                    "Electron pT MTD - 2 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2179                    30,
2180                    10,
2181                    100);
2182   meEle_eta_MTD_2sigma_Sig_EB_ = ibook.book1D(
2183       "Ele_eta_MTD_2sigma_Sig_EB", "Electron eta MTD - 2 sigma compatibility - Signal Barrel;#eta;Counts", 32, 0., 1.6);
2184   meEle_phi_MTD_2sigma_Sig_EB_ = ibook.book1D("Ele_phi_MTD_2sigma_Sig_EB",
2185                                               "Electron phi MTD - 2 sigma compatibility - Signal Barrel;#phi;Counts",
2186                                               64,
2187                                               -3.2,
2188                                               3.2);
2189 
2190   meEleISO_Ntracks_Sig_EE_ = ibook.book1D("Ele_Iso_Ntracks_Sig_EE",
2191                                           "Number of tracks in isolation cone around electron track after basic cuts - "
2192                                           "Signal Endcap;Number of tracks;Counts",
2193                                           20,
2194                                           0,
2195                                           20);
2196   meEleISO_chIso_Sig_EE_ = ibook.book1D(
2197       "Ele_chIso_sum_Sig_EE",
2198       "Track pT sum in isolation cone around electron track after basic cuts - Signal Endcap;p_{T} (GeV);Counts",
2199       nbin_2,
2200       0,
2201       20);
2202   meEleISO_rel_chIso_Sig_EE_ = ibook.book1D(
2203       "Ele_rel_chIso_sum_Sig_EE",
2204       "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Endcap;Isolation;Counts",
2205       nbin_1,
2206       0,
2207       4);
2208 
2209   if (optionalPlots_) {
2210     meEle_pt_sim_MTD_4sigma_Sig_EB_ =
2211         ibook.book1D("Ele_pT_sim_MTD_4sigma_Sig_EB",
2212                      "Electron pT MTD SIM - 4 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2213                      30,
2214                      10,
2215                      100);
2216     meEle_pt_sim_MTD_3sigma_Sig_EB_ =
2217         ibook.book1D("Ele_pT_sim_MTD_3sigma_Sig_EB",
2218                      "Electron pT MTD SIM - 3 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2219                      30,
2220                      10,
2221                      100);
2222     meEle_pt_sim_MTD_2sigma_Sig_EB_ =
2223         ibook.book1D("Ele_pT_sim_MTD_2sigma_Sig_EB",
2224                      "Electron pT MTD SIM - 2 sigma compatibility - Signal Barrel;p_{T} (GeV);Counts",
2225                      30,
2226                      10,
2227                      100);
2228 
2229     meEleISO_Ntracks_MTD_1_Sig_EE_ =
2230         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Sig_EE",
2231                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2232                      20,
2233                      0,
2234                      20);
2235     meEleISO_chIso_MTD_1_Sig_EE_ = ibook.book1D(
2236         "Ele_chIso_sum_MTD_1_Sig_EE",
2237         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2238         nbin_2,
2239         0,
2240         20);
2241     meEleISO_rel_chIso_MTD_1_Sig_EE_ = ibook.book1D(
2242         "Ele_rel_chIso_sum_MTD_1_Sig_EE",
2243         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2244         nbin_1,
2245         0,
2246         4);
2247 
2248     meEleISO_Ntracks_MTD_2_Sig_EE_ =
2249         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Sig_EE",
2250                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2251                      20,
2252                      0,
2253                      20);
2254     meEleISO_chIso_MTD_2_Sig_EE_ = ibook.book1D(
2255         "Ele_chIso_sum_MTD_2_Sig_EE",
2256         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2257         nbin_2,
2258         0,
2259         20);
2260     meEleISO_rel_chIso_MTD_2_Sig_EE_ = ibook.book1D(
2261         "Ele_rel_chIso_sum_MTD_2_Sig_EE",
2262         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2263         nbin_1,
2264         0,
2265         4);
2266 
2267     meEleISO_Ntracks_gen_Sig_EE_ =
2268         ibook.book1D("Ele_Iso_Ntracks_gen_Sig_EE",
2269                      "Tracks in isolation cone around electron track after basic cuts using genInfo - Signal Endcap",
2270                      20,
2271                      0,
2272                      20);
2273     meEleISO_chIso_gen_Sig_EE_ =
2274         ibook.book1D("Ele_chIso_sum_gen_Sig_EE",
2275                      "Track pT sum in isolation cone around electron track after basic cuts - Signal Endcap",
2276                      nbin_2,
2277                      0,
2278                      20);
2279     meEleISO_rel_chIso_gen_Sig_EE_ =
2280         ibook.book1D("Ele_rel_chIso_sum_gen_Sig_EE",
2281                      "Track relative pT sum in isolation cone around electron track after basic cuts - Signal Endcap",
2282                      nbin_1,
2283                      0,
2284                      4);
2285 
2286     meEleISO_Ntracks_MTD_3_Sig_EE_ =
2287         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Sig_EE",
2288                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2289                      20,
2290                      0,
2291                      20);
2292     meEleISO_chIso_MTD_3_Sig_EE_ = ibook.book1D(
2293         "Ele_chIso_sum_MTD_3_Sig_EE",
2294         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2295         nbin_2,
2296         0,
2297         20);
2298     meEleISO_rel_chIso_MTD_3_Sig_EE_ = ibook.book1D(
2299         "Ele_rel_chIso_sum_MTD_3_Sig_EE",
2300         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2301         nbin_1,
2302         0,
2303         4);
2304 
2305     meEleISO_Ntracks_MTD_4_Sig_EE_ =
2306         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Sig_EE",
2307                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2308                      20,
2309                      0,
2310                      20);
2311     meEleISO_chIso_MTD_4_Sig_EE_ = ibook.book1D(
2312         "Ele_chIso_sum_MTD_4_Sig_EE",
2313         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2314         nbin_2,
2315         0,
2316         20);
2317     meEleISO_rel_chIso_MTD_4_Sig_EE_ = ibook.book1D(
2318         "Ele_rel_chIso_sum_MTD_4_Sig_EE",
2319         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2320         nbin_1,
2321         0,
2322         4);
2323 
2324     meEleISO_Ntracks_MTD_5_Sig_EE_ =
2325         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Sig_EE",
2326                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2327                      20,
2328                      0,
2329                      20);
2330     meEleISO_chIso_MTD_5_Sig_EE_ = ibook.book1D(
2331         "Ele_chIso_sum_MTD_5_Sig_EE",
2332         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2333         nbin_2,
2334         0,
2335         20);
2336     meEleISO_rel_chIso_MTD_5_Sig_EE_ = ibook.book1D(
2337         "Ele_rel_chIso_sum_MTD_5_Sig_EE",
2338         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2339         nbin_1,
2340         0,
2341         4);
2342 
2343     meEleISO_Ntracks_MTD_6_Sig_EE_ =
2344         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Sig_EE",
2345                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2346                      20,
2347                      0,
2348                      20);
2349     meEleISO_chIso_MTD_6_Sig_EE_ = ibook.book1D(
2350         "Ele_chIso_sum_MTD_6_Sig_EE",
2351         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2352         nbin_2,
2353         0,
2354         20);
2355     meEleISO_rel_chIso_MTD_6_Sig_EE_ = ibook.book1D(
2356         "Ele_rel_chIso_sum_MTD_6_Sig_EE",
2357         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2358         nbin_1,
2359         0,
2360         4);
2361 
2362     meEleISO_Ntracks_MTD_7_Sig_EE_ =
2363         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Sig_EE",
2364                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2365                      20,
2366                      0,
2367                      20);
2368     meEleISO_chIso_MTD_7_Sig_EE_ = ibook.book1D(
2369         "Ele_chIso_sum_MTD_7_Sig_EE",
2370         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2371         nbin_2,
2372         0,
2373         20);
2374     meEleISO_rel_chIso_MTD_7_Sig_EE_ = ibook.book1D(
2375         "Ele_rel_chIso_sum_MTD_7_Sig_EE",
2376         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2377         nbin_1,
2378         0,
2379         4);
2380 
2381     meEleISO_Ntracks_MTD_sim_1_Sig_EE_ =
2382         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Sig_EE",
2383                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2384                      20,
2385                      0,
2386                      20);
2387     meEleISO_chIso_MTD_sim_1_Sig_EE_ = ibook.book1D(
2388         "Ele_chIso_sum_MTD_sim_1_Sig_EE",
2389         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2390         nbin_2,
2391         0,
2392         20);
2393     meEleISO_rel_chIso_MTD_sim_1_Sig_EE_ = ibook.book1D(
2394         "Ele_rel_chIso_sum_MTD_sim_1_Sig_EE",
2395         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2396         nbin_1,
2397         0,
2398         4);
2399 
2400     meEleISO_Ntracks_MTD_sim_2_Sig_EE_ =
2401         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Sig_EE",
2402                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2403                      20,
2404                      0,
2405                      20);
2406     meEleISO_chIso_MTD_sim_2_Sig_EE_ = ibook.book1D(
2407         "Ele_chIso_sum_MTD_sim_2_Sig_EE",
2408         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2409         nbin_2,
2410         0,
2411         20);
2412     meEleISO_rel_chIso_MTD_sim_2_Sig_EE_ = ibook.book1D(
2413         "Ele_rel_chIso_sum_MTD_sim_2_Sig_EE",
2414         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2415         nbin_1,
2416         0,
2417         4);
2418 
2419     meEleISO_Ntracks_MTD_sim_3_Sig_EE_ =
2420         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Sig_EE",
2421                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2422                      20,
2423                      0,
2424                      20);
2425     meEleISO_chIso_MTD_sim_3_Sig_EE_ = ibook.book1D(
2426         "Ele_chIso_sum_MTD_sim_3_Sig_EE",
2427         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2428         nbin_2,
2429         0,
2430         20);
2431     meEleISO_rel_chIso_MTD_sim_3_Sig_EE_ = ibook.book1D(
2432         "Ele_rel_chIso_sum_MTD_sim_3_Sig_EE",
2433         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2434         nbin_1,
2435         0,
2436         4);
2437 
2438     meEleISO_Ntracks_MTD_sim_4_Sig_EE_ =
2439         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Sig_EE",
2440                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2441                      20,
2442                      0,
2443                      20);
2444     meEleISO_chIso_MTD_sim_4_Sig_EE_ = ibook.book1D(
2445         "Ele_chIso_sum_MTD_sim_4_Sig_EE",
2446         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2447         nbin_2,
2448         0,
2449         20);
2450     meEleISO_rel_chIso_MTD_sim_4_Sig_EE_ = ibook.book1D(
2451         "Ele_rel_chIso_sum_MTD_sim_4_Sig_EE",
2452         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2453         nbin_1,
2454         0,
2455         4);
2456 
2457     meEleISO_Ntracks_MTD_sim_5_Sig_EE_ =
2458         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Sig_EE",
2459                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2460                      20,
2461                      0,
2462                      20);
2463     meEleISO_chIso_MTD_sim_5_Sig_EE_ = ibook.book1D(
2464         "Ele_chIso_sum_MTD_sim_5_Sig_EE",
2465         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2466         nbin_2,
2467         0,
2468         20);
2469     meEleISO_rel_chIso_MTD_sim_5_Sig_EE_ = ibook.book1D(
2470         "Ele_rel_chIso_sum_MTD_sim_5_Sig_EE",
2471         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2472         nbin_1,
2473         0,
2474         4);
2475 
2476     meEleISO_Ntracks_MTD_sim_6_Sig_EE_ =
2477         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Sig_EE",
2478                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2479                      20,
2480                      0,
2481                      20);
2482     meEleISO_chIso_MTD_sim_6_Sig_EE_ = ibook.book1D(
2483         "Ele_chIso_sum_MTD_sim_6_Sig_EE",
2484         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2485         nbin_2,
2486         0,
2487         20);
2488     meEleISO_rel_chIso_MTD_sim_6_Sig_EE_ = ibook.book1D(
2489         "Ele_rel_chIso_sum_MTD_sim_6_Sig_EE",
2490         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2491         nbin_1,
2492         0,
2493         4);
2494 
2495     meEleISO_Ntracks_MTD_sim_7_Sig_EE_ =
2496         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Sig_EE",
2497                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2498                      20,
2499                      0,
2500                      20);
2501     meEleISO_chIso_MTD_sim_7_Sig_EE_ = ibook.book1D(
2502         "Ele_chIso_sum_MTD_sim_7_Sig_EE",
2503         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2504         nbin_2,
2505         0,
2506         20);
2507     meEleISO_rel_chIso_MTD_sim_7_Sig_EE_ = ibook.book1D(
2508         "Ele_rel_chIso_sum_MTD_sim_7_Sig_EE",
2509         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2510         nbin_1,
2511         0,
2512         4);
2513   }
2514   meEleISO_Ntracks_MTD_4sigma_Sig_EE_ =
2515       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Sig_EE",
2516                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma significance - "
2517                    "Signal Endcap;Number of tracks;Counts",
2518                    20,
2519                    0,
2520                    20);
2521   meEleISO_chIso_MTD_4sigma_Sig_EE_ =
2522       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Sig_EE",
2523                    "Track pT sum in isolation cone around electron track after basic "
2524                    "cuts with MTD - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2525                    nbin_2,
2526                    0,
2527                    20);
2528   meEleISO_rel_chIso_MTD_4sigma_Sig_EE_ =
2529       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Sig_EE",
2530                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 4 sigma "
2531                    "significance - Signal Endcap;Isolation;Counts",
2532                    nbin_1,
2533                    0,
2534                    4);
2535 
2536   meEleISO_Ntracks_MTD_3sigma_Sig_EE_ =
2537       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Sig_EE",
2538                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma significance - "
2539                    "Signal Endcap;Number of tracks;Counts",
2540                    20,
2541                    0,
2542                    20);
2543   meEleISO_chIso_MTD_3sigma_Sig_EE_ =
2544       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Sig_EE",
2545                    "Track pT sum in isolation cone around electron track after basic "
2546                    "cuts with MTD - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2547                    nbin_2,
2548                    0,
2549                    20);
2550   meEleISO_rel_chIso_MTD_3sigma_Sig_EE_ =
2551       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Sig_EE",
2552                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 3 sigma "
2553                    "significance - Signal Endcap;Isolation;Counts",
2554                    nbin_1,
2555                    0,
2556                    4);
2557 
2558   meEleISO_Ntracks_MTD_2sigma_Sig_EE_ =
2559       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Sig_EE",
2560                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma significance - "
2561                    "Signal Endcap;Number of tracks;Counts",
2562                    20,
2563                    0,
2564                    20);
2565   meEleISO_chIso_MTD_2sigma_Sig_EE_ =
2566       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Sig_EE",
2567                    "Track pT sum in isolation cone around electron track after basic "
2568                    "cuts with MTD - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2569                    nbin_2,
2570                    0,
2571                    20);
2572   meEleISO_rel_chIso_MTD_2sigma_Sig_EE_ =
2573       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Sig_EE",
2574                    "Track relative pT sum in isolation cone around electron track after basic cuts with MTD - 2 sigma "
2575                    "significance - Signal Endcap;Isolation;Counts",
2576                    nbin_1,
2577                    0,
2578                    4);
2579 
2580   meEle_pt_tot_Sig_EE_ =
2581       ibook.book1D("Ele_pT_tot_Sig_EE", "Electron pT tot - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2582   meEle_pt_noMTD_Sig_EE_ =
2583       ibook.book1D("Ele_pT_noMTD_Sig_EE", "Electron pT noMTD - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2584 
2585   meEle_pt_sim_tot_Sig_EE_ =
2586       ibook.book1D("Ele_pT_sim_tot_Sig_EE", "Electron pT tot - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2587 
2588   meEle_eta_tot_Sig_EE_ =
2589       ibook.book1D("Ele_eta_tot_Sig_EE", "Electron eta tot - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2590   meEle_eta_noMTD_Sig_EE_ =
2591       ibook.book1D("Ele_eta_noMTD_Sig_EE", "Electron eta noMTD - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2592 
2593   meEle_phi_tot_Sig_EE_ =
2594       ibook.book1D("Ele_phi_tot_Sig_EE", "Electron phi tot - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2595   meEle_phi_noMTD_Sig_EE_ =
2596       ibook.book1D("Ele_phi_noMTD_Sig_EE", "Electron phi noMTD - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2597 
2598   if (optionalPlots_) {
2599     meEleISO_Ntracks_MTD_sim_4sigma_Sig_EE_ =
2600         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Sig_EE",
2601                      "Number of tracks in isolation cone around electron track after basic cuts with MTD SIM - 4 sigma "
2602                      "significance - Signal Endcap;Number of tracks;Counts",
2603                      20,
2604                      0,
2605                      20);
2606     meEleISO_chIso_MTD_sim_4sigma_Sig_EE_ =
2607         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Sig_EE",
2608                      "Track pT sum in isolation cone around electron track after "
2609                      "basic cuts with MTD SIM - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2610                      nbin_2,
2611                      0,
2612                      20);
2613     meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EE_ =
2614         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Sig_EE",
2615                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
2616                      "sigma significance - Signal Endcap;Isolation;Counts",
2617                      nbin_1,
2618                      0,
2619                      4);
2620 
2621     meEleISO_Ntracks_MTD_sim_3sigma_Sig_EE_ =
2622         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Sig_EE",
2623                      "Tracks in isolation cone around electron track after basic "
2624                      "cuts with MTD SIM - 3 sigma significance - Signal Endcap;Number of tracks;Counts",
2625                      20,
2626                      0,
2627                      20);
2628     meEleISO_chIso_MTD_sim_3sigma_Sig_EE_ =
2629         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Sig_EE",
2630                      "Track pT sum in isolation cone around electron track after "
2631                      "basic cuts with MTD SIM - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2632                      nbin_2,
2633                      0,
2634                      20);
2635     meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EE_ =
2636         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Sig_EE",
2637                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
2638                      "sigma significance - Signal Endcap;Isolation;Counts",
2639                      nbin_1,
2640                      0,
2641                      4);
2642 
2643     meEleISO_Ntracks_MTD_sim_2sigma_Sig_EE_ =
2644         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Sig_EE",
2645                      "Tracks in isolation cone around electron track after basic "
2646                      "cuts with MTD SIM - 2 sigma significance - Signal Endcap;Number of tracks;Counts",
2647                      20,
2648                      0,
2649                      20);
2650     meEleISO_chIso_MTD_sim_2sigma_Sig_EE_ =
2651         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Sig_EE",
2652                      "Track pT sum in isolation cone around electron track after "
2653                      "basic cuts with MTD SIM - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2654                      nbin_2,
2655                      0,
2656                      20);
2657     meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EE_ =
2658         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Sig_EE",
2659                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
2660                      "sigma significance - Signal Endcap;Isolation;Counts",
2661                      nbin_1,
2662                      0,
2663                      4);
2664 
2665     meEle_pt_MTD_1_Sig_EE_ = ibook.book1D("Ele_pT_MTD_1_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2666     meEle_eta_MTD_1_Sig_EE_ = ibook.book1D("Ele_eta_MTD_1_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2667     meEle_phi_MTD_1_Sig_EE_ = ibook.book1D("Ele_phi_MTD_1_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2668     meEle_pt_gen_Sig_EE_ =
2669         ibook.book1D("Ele_pT_gen_Sig_EE", "Electron pT genInfo - Signal Endcap;p_{T} (GeV);Counts", 30, 10, 100);
2670     meEle_eta_gen_Sig_EE_ =
2671         ibook.book1D("Ele_eta_gen_Sig_EE", "Electron eta genInfo - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2672     meEle_phi_gen_Sig_EE_ =
2673         ibook.book1D("Ele_phi_gen_Sig_EE", "Electron phi genInfo - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2674 
2675     meEle_pt_MTD_2_Sig_EE_ = ibook.book1D("Ele_pT_MTD_2_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2676     meEle_eta_MTD_2_Sig_EE_ = ibook.book1D("Ele_eta_MTD_2_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2677     meEle_phi_MTD_2_Sig_EE_ = ibook.book1D("Ele_phi_MTD_2_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2678 
2679     meEle_pt_MTD_3_Sig_EE_ = ibook.book1D("Ele_pT_MTD_3_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2680     meEle_eta_MTD_3_Sig_EE_ = ibook.book1D("Ele_eta_MTD_3_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2681     meEle_phi_MTD_3_Sig_EE_ = ibook.book1D("Ele_phi_MTD_3_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2682 
2683     meEle_pt_MTD_4_Sig_EE_ = ibook.book1D("Ele_pT_MTD_4_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2684     meEle_eta_MTD_4_Sig_EE_ = ibook.book1D("Ele_eta_MTD_4_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2685     meEle_phi_MTD_4_Sig_EE_ = ibook.book1D("Ele_phi_MTD_4_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2686 
2687     meEle_pt_MTD_5_Sig_EE_ = ibook.book1D("Ele_pT_MTD_5_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2688     meEle_eta_MTD_5_Sig_EE_ = ibook.book1D("Ele_eta_MTD_5_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2689     meEle_phi_MTD_5_Sig_EE_ = ibook.book1D("Ele_phi_MTD_5_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2690 
2691     meEle_pt_MTD_6_Sig_EE_ = ibook.book1D("Ele_pT_MTD_6_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2692     meEle_eta_MTD_6_Sig_EE_ = ibook.book1D("Ele_eta_MTD_6_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2693     meEle_phi_MTD_6_Sig_EE_ = ibook.book1D("Ele_phi_MTD_6_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2694 
2695     meEle_pt_MTD_7_Sig_EE_ = ibook.book1D("Ele_pT_MTD_7_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2696     meEle_eta_MTD_7_Sig_EE_ = ibook.book1D("Ele_eta_MTD_7_Sig_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
2697     meEle_phi_MTD_7_Sig_EE_ = ibook.book1D("Ele_phi_MTD_7_Sig_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
2698 
2699     meEle_pt_sim_MTD_1_Sig_EE_ =
2700         ibook.book1D("Ele_pT_sim_MTD_1_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2701     meEle_pt_sim_MTD_2_Sig_EE_ =
2702         ibook.book1D("Ele_pT_sim_MTD_2_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2703     meEle_pt_sim_MTD_3_Sig_EE_ =
2704         ibook.book1D("Ele_pT_sim_MTD_3_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2705     meEle_pt_sim_MTD_4_Sig_EE_ =
2706         ibook.book1D("Ele_pT_sim_MTD_4_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2707     meEle_pt_sim_MTD_5_Sig_EE_ =
2708         ibook.book1D("Ele_pT_sim_MTD_5_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2709     meEle_pt_sim_MTD_6_Sig_EE_ =
2710         ibook.book1D("Ele_pT_sim_MTD_6_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2711     meEle_pt_sim_MTD_7_Sig_EE_ =
2712         ibook.book1D("Ele_pT_sim_MTD_7_Sig_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
2713 
2714     meEle_pt_sim_MTD_4sigma_Sig_EE_ =
2715         ibook.book1D("Ele_pT_sim_MTD_4sigma_Sig_EE",
2716                      "Electron pT MTD SIM - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2717                      30,
2718                      10,
2719                      100);
2720     meEle_pt_sim_MTD_3sigma_Sig_EE_ =
2721         ibook.book1D("Ele_pT_sim_MTD_3sigma_Sig_EE",
2722                      "Electron pT MTD SIM - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2723                      30,
2724                      10,
2725                      100);
2726     meEle_pt_sim_MTD_2sigma_Sig_EE_ =
2727         ibook.book1D("Ele_pT_sim_MTD_2sigma_Sig_EE",
2728                      "Electron pT MTD SIM - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2729                      30,
2730                      10,
2731                      100);
2732   }
2733 
2734   meEle_pt_MTD_4sigma_Sig_EE_ =
2735       ibook.book1D("Ele_pT_MTD_4sigma_Sig_EE",
2736                    "Electron pT MTD - 4 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2737                    30,
2738                    10,
2739                    100);
2740   meEle_eta_MTD_4sigma_Sig_EE_ = ibook.book1D(
2741       "Ele_eta_MTD_4sigma_Sig_EE", "Electron eta MTD - 4 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2742   meEle_phi_MTD_4sigma_Sig_EE_ = ibook.book1D(
2743       "Ele_phi_MTD_4sigma_Sig_EE", "Electron phi MTD - 4 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2744 
2745   meEle_pt_MTD_3sigma_Sig_EE_ =
2746       ibook.book1D("Ele_pT_MTD_3sigma_Sig_EE",
2747                    "Electron pT MTD - 3 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2748                    30,
2749                    10,
2750                    100);
2751   meEle_eta_MTD_3sigma_Sig_EE_ = ibook.book1D(
2752       "Ele_eta_MTD_3sigma_Sig_EE", "Electron eta MTD - 3 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2753   meEle_phi_MTD_3sigma_Sig_EE_ = ibook.book1D(
2754       "Ele_phi_MTD_3sigma_Sig_EE", "Electron phi MTD - 3 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2755 
2756   meEle_pt_MTD_2sigma_Sig_EE_ =
2757       ibook.book1D("Ele_pT_MTD_2sigma_Sig_EE",
2758                    "Electron pT MTD - 2 sigma significance - Signal Endcap;p_{T} (GeV);Counts",
2759                    30,
2760                    10,
2761                    100);
2762   meEle_eta_MTD_2sigma_Sig_EE_ = ibook.book1D(
2763       "Ele_eta_MTD_2sigma_Sig_EE", "Electron eta MTD - 2 sigma significance - Signal Endcap;#eta;Counts", 32, 1.6, 3.2);
2764   meEle_phi_MTD_2sigma_Sig_EE_ = ibook.book1D(
2765       "Ele_phi_MTD_2sigma_Sig_EE", "Electron phi MTD - 2 sigma significance - Signal Endcap;#phi;Counts", 64, -3.2, 3.2);
2766 
2767   // background
2768   meEleISO_Ntracks_Bkg_EB_ = ibook.book1D(
2769       "Ele_Iso_Ntracks_Bkg_EB",
2770       "Number of tracks in isolation cone around electron track after basic cuts - Bkg Barrel;Number of tracks;Counts",
2771       20,
2772       0,
2773       20);
2774   meEleISO_chIso_Bkg_EB_ = ibook.book1D(
2775       "Ele_chIso_sum_Bkg_EB",
2776       "Track pT sum in isolation cone around electron track after basic cuts - Bkg Barrel;p_{T} (GeV);Counts",
2777       nbin_2,
2778       0,
2779       20);
2780   meEleISO_rel_chIso_Bkg_EB_ = ibook.book1D(
2781       "Ele_rel_chIso_sum_Bkg_EB",
2782       "Track relative pT sum in isolation cone around electron track after basic cuts - Bkg Barrel;Isolation;Counts",
2783       nbin_1,
2784       0,
2785       4);
2786   if (optionalPlots_) {
2787     meEleISO_Ntracks_MTD_1_Bkg_EB_ =
2788         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Bkg_EB",
2789                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2790                      20,
2791                      0,
2792                      20);
2793     meEleISO_chIso_MTD_1_Bkg_EB_ = ibook.book1D(
2794         "Ele_chIso_sum_MTD_1_Bkg_EB",
2795         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2796         nbin_2,
2797         0,
2798         20);
2799     meEleISO_rel_chIso_MTD_1_Bkg_EB_ = ibook.book1D(
2800         "Ele_rel_chIso_sum_MTD_1_Bkg_EB",
2801         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2802         nbin_1,
2803         0,
2804         4);
2805 
2806     meEleISO_Ntracks_MTD_2_Bkg_EB_ =
2807         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Bkg_EB",
2808                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2809                      20,
2810                      0,
2811                      20);
2812     meEleISO_chIso_MTD_2_Bkg_EB_ = ibook.book1D(
2813         "Ele_chIso_sum_MTD_2_Bkg_EB",
2814         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2815         nbin_2,
2816         0,
2817         20);
2818     meEleISO_rel_chIso_MTD_2_Bkg_EB_ = ibook.book1D(
2819         "Ele_rel_chIso_sum_MTD_2_Bkg_EB",
2820         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2821         nbin_1,
2822         0,
2823         4);
2824     meEleISO_Ntracks_gen_Bkg_EB_ = ibook.book1D("Ele_Iso_Ntracks_gen_Bkg_EB",
2825                                                 "Tracks in isolation cone around electron track after basic cuts using "
2826                                                 "genInfo - Bkg Barrel;Number of tracks;Counts",
2827                                                 20,
2828                                                 0,
2829                                                 20);
2830     meEleISO_chIso_gen_Bkg_EB_ = ibook.book1D("Ele_chIso_sum_gen_Bkg_EB",
2831                                               "Track pT sum in isolation cone around electron track after basic cuts "
2832                                               "using genInfo - Bkg Barrel;p_{T} (GeV);Counts",
2833                                               nbin_2,
2834                                               0,
2835                                               20);
2836     meEleISO_rel_chIso_gen_Bkg_EB_ = ibook.book1D("Ele_rel_chIso_sum_gen_Bkg_EB",
2837                                                   "Track relative pT sum in isolation cone around electron track after "
2838                                                   "basic cuts using genInfo - Bkg Barrel;Isolation;Counts",
2839                                                   nbin_1,
2840                                                   0,
2841                                                   4);
2842     meEleISO_Ntracks_MTD_3_Bkg_EB_ =
2843         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Bkg_EB",
2844                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2845                      20,
2846                      0,
2847                      20);
2848     meEleISO_chIso_MTD_3_Bkg_EB_ = ibook.book1D(
2849         "Ele_chIso_sum_MTD_3_Bkg_EB",
2850         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2851         nbin_2,
2852         0,
2853         20);
2854     meEleISO_rel_chIso_MTD_3_Bkg_EB_ = ibook.book1D(
2855         "Ele_rel_chIso_sum_MTD_3_Bkg_EB",
2856         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2857         nbin_1,
2858         0,
2859         4);
2860 
2861     meEleISO_Ntracks_MTD_4_Bkg_EB_ =
2862         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Bkg_EB",
2863                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2864                      20,
2865                      0,
2866                      20);
2867     meEleISO_chIso_MTD_4_Bkg_EB_ = ibook.book1D(
2868         "Ele_chIso_sum_MTD_4_Bkg_EB",
2869         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2870         nbin_2,
2871         0,
2872         20);
2873     meEleISO_rel_chIso_MTD_4_Bkg_EB_ = ibook.book1D(
2874         "Ele_rel_chIso_sum_MTD_4_Bkg_EB",
2875         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2876         nbin_1,
2877         0,
2878         4);
2879 
2880     meEleISO_Ntracks_MTD_5_Bkg_EB_ =
2881         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Bkg_EB",
2882                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2883                      20,
2884                      0,
2885                      20);
2886     meEleISO_chIso_MTD_5_Bkg_EB_ = ibook.book1D(
2887         "Ele_chIso_sum_MTD_5_Bkg_EB",
2888         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2889         nbin_2,
2890         0,
2891         20);
2892     meEleISO_rel_chIso_MTD_5_Bkg_EB_ = ibook.book1D(
2893         "Ele_rel_chIso_sum_MTD_5_Bkg_EB",
2894         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2895         nbin_1,
2896         0,
2897         4);
2898 
2899     meEleISO_Ntracks_MTD_6_Bkg_EB_ =
2900         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Bkg_EB",
2901                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2902                      20,
2903                      0,
2904                      20);
2905     meEleISO_chIso_MTD_6_Bkg_EB_ = ibook.book1D(
2906         "Ele_chIso_sum_MTD_6_Bkg_EB",
2907         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2908         nbin_2,
2909         0,
2910         20);
2911     meEleISO_rel_chIso_MTD_6_Bkg_EB_ = ibook.book1D(
2912         "Ele_rel_chIso_sum_MTD_6_Bkg_EB",
2913         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2914         nbin_1,
2915         0,
2916         4);
2917 
2918     meEleISO_Ntracks_MTD_7_Bkg_EB_ =
2919         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Bkg_EB",
2920                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2921                      20,
2922                      0,
2923                      20);
2924     meEleISO_chIso_MTD_7_Bkg_EB_ = ibook.book1D(
2925         "Ele_chIso_sum_MTD_7_Bkg_EB",
2926         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2927         nbin_2,
2928         0,
2929         20);
2930     meEleISO_rel_chIso_MTD_7_Bkg_EB_ = ibook.book1D(
2931         "Ele_rel_chIso_sum_MTD_7_Bkg_EB",
2932         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2933         nbin_1,
2934         0,
2935         4);
2936 
2937     meEleISO_Ntracks_MTD_sim_1_Bkg_EB_ =
2938         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Bkg_EB",
2939                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2940                      20,
2941                      0,
2942                      20);
2943     meEleISO_chIso_MTD_sim_1_Bkg_EB_ = ibook.book1D(
2944         "Ele_chIso_sum_MTD_sim_1_Bkg_EB",
2945         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2946         nbin_2,
2947         0,
2948         20);
2949     meEleISO_rel_chIso_MTD_sim_1_Bkg_EB_ = ibook.book1D(
2950         "Ele_rel_chIso_sum_MTD_sim_1_Bkg_EB",
2951         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2952         nbin_1,
2953         0,
2954         4);
2955 
2956     meEleISO_Ntracks_MTD_sim_2_Bkg_EB_ =
2957         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Bkg_EB",
2958                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2959                      20,
2960                      0,
2961                      20);
2962     meEleISO_chIso_MTD_sim_2_Bkg_EB_ = ibook.book1D(
2963         "Ele_chIso_sum_MTD_sim_2_Bkg_EB",
2964         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2965         nbin_2,
2966         0,
2967         20);
2968     meEleISO_rel_chIso_MTD_sim_2_Bkg_EB_ = ibook.book1D(
2969         "Ele_rel_chIso_sum_MTD_sim_2_Bkg_EB",
2970         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2971         nbin_1,
2972         0,
2973         4);
2974 
2975     meEleISO_Ntracks_MTD_sim_3_Bkg_EB_ =
2976         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Bkg_EB",
2977                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2978                      20,
2979                      0,
2980                      20);
2981     meEleISO_chIso_MTD_sim_3_Bkg_EB_ = ibook.book1D(
2982         "Ele_chIso_sum_MTD_sim_3_Bkg_EB",
2983         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
2984         nbin_2,
2985         0,
2986         20);
2987     meEleISO_rel_chIso_MTD_sim_3_Bkg_EB_ = ibook.book1D(
2988         "Ele_rel_chIso_sum_MTD_sim_3_Bkg_EB",
2989         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
2990         nbin_1,
2991         0,
2992         4);
2993 
2994     meEleISO_Ntracks_MTD_sim_4_Bkg_EB_ =
2995         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Bkg_EB",
2996                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
2997                      20,
2998                      0,
2999                      20);
3000     meEleISO_chIso_MTD_sim_4_Bkg_EB_ = ibook.book1D(
3001         "Ele_chIso_sum_MTD_sim_4_Bkg_EB",
3002         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3003         nbin_2,
3004         0,
3005         20);
3006     meEleISO_rel_chIso_MTD_sim_4_Bkg_EB_ = ibook.book1D(
3007         "Ele_rel_chIso_sum_MTD_sim_4_Bkg_EB",
3008         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3009         nbin_1,
3010         0,
3011         4);
3012 
3013     meEleISO_Ntracks_MTD_sim_5_Bkg_EB_ =
3014         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Bkg_EB",
3015                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3016                      20,
3017                      0,
3018                      20);
3019     meEleISO_chIso_MTD_sim_5_Bkg_EB_ = ibook.book1D(
3020         "Ele_chIso_sum_MTD_sim_5_Bkg_EB",
3021         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3022         nbin_2,
3023         0,
3024         20);
3025     meEleISO_rel_chIso_MTD_sim_5_Bkg_EB_ = ibook.book1D(
3026         "Ele_rel_chIso_sum_MTD_sim_5_Bkg_EB",
3027         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3028         nbin_1,
3029         0,
3030         4);
3031 
3032     meEleISO_Ntracks_MTD_sim_6_Bkg_EB_ =
3033         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Bkg_EB",
3034                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3035                      20,
3036                      0,
3037                      20);
3038     meEleISO_chIso_MTD_sim_6_Bkg_EB_ = ibook.book1D(
3039         "Ele_chIso_sum_MTD_sim_6_Bkg_EB",
3040         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3041         nbin_2,
3042         0,
3043         20);
3044     meEleISO_rel_chIso_MTD_sim_6_Bkg_EB_ = ibook.book1D(
3045         "Ele_rel_chIso_sum_MTD_sim_6_Bkg_EB",
3046         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3047         nbin_1,
3048         0,
3049         4);
3050 
3051     meEleISO_Ntracks_MTD_sim_7_Bkg_EB_ =
3052         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Bkg_EB",
3053                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3054                      20,
3055                      0,
3056                      20);
3057     meEleISO_chIso_MTD_sim_7_Bkg_EB_ = ibook.book1D(
3058         "Ele_chIso_sum_MTD_sim_7_Bkg_EB",
3059         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3060         nbin_2,
3061         0,
3062         20);
3063     meEleISO_rel_chIso_MTD_sim_7_Bkg_EB_ = ibook.book1D(
3064         "Ele_rel_chIso_sum_MTD_sim_7_Bkg_EB",
3065         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3066         nbin_1,
3067         0,
3068         4);
3069   }
3070   meEleISO_Ntracks_MTD_4sigma_Bkg_EB_ =
3071       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Bkg_EB",
3072                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma significance - "
3073                    "Bkg Barrel;Number of tracks;Counts",
3074                    20,
3075                    0,
3076                    20);
3077   meEleISO_chIso_MTD_4sigma_Bkg_EB_ =
3078       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Bkg_EB",
3079                    "Track pT sum in isolation cone around electron track after basic "
3080                    "cuts with MTD - 4 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3081                    nbin_2,
3082                    0,
3083                    20);
3084   meEleISO_rel_chIso_MTD_4sigma_Bkg_EB_ =
3085       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Bkg_EB",
3086                    "Track relative pT sum in isolation cone around electron track "
3087                    "after basic cuts with MTD - 4 sigma significance - Bkg Barrel;Isolation;Counts",
3088                    nbin_1,
3089                    0,
3090                    4);
3091 
3092   meEleISO_Ntracks_MTD_3sigma_Bkg_EB_ =
3093       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Bkg_EB",
3094                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma significance - "
3095                    "Bkg Barrel;Number of tracks;Counts",
3096                    20,
3097                    0,
3098                    20);
3099   meEleISO_chIso_MTD_3sigma_Bkg_EB_ =
3100       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Bkg_EB",
3101                    "Track pT sum in isolation cone around electron track after basic "
3102                    "cuts with MTD - 3 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3103                    nbin_2,
3104                    0,
3105                    20);
3106   meEleISO_rel_chIso_MTD_3sigma_Bkg_EB_ =
3107       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Bkg_EB",
3108                    "Track relative pT sum in isolation cone around electron track "
3109                    "after basic cuts with MTD - 3 sigma significance - Bkg Barrel;Isolation;Counts",
3110                    nbin_1,
3111                    0,
3112                    4);
3113 
3114   meEleISO_Ntracks_MTD_2sigma_Bkg_EB_ =
3115       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Bkg_EB",
3116                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma significance - "
3117                    "Bkg Barrel;Number of tracks;Counts",
3118                    20,
3119                    0,
3120                    20);
3121   meEleISO_chIso_MTD_2sigma_Bkg_EB_ =
3122       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Bkg_EB",
3123                    "Track pT sum in isolation cone around electron track after basic "
3124                    "cuts with MTD - 2 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3125                    nbin_2,
3126                    0,
3127                    20);
3128   meEleISO_rel_chIso_MTD_2sigma_Bkg_EB_ =
3129       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Bkg_EB",
3130                    "Track relative pT sum in isolation cone around electron track "
3131                    "after basic cuts with MTD - 2 sigma significance - Bkg Barrel;Isolation;Counts",
3132                    nbin_1,
3133                    0,
3134                    4);
3135 
3136   meEle_pt_tot_Bkg_EB_ =
3137       ibook.book1D("Ele_pT_tot_Bkg_EB", "Electron pT tot - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3138   meEle_pt_noMTD_Bkg_EB_ =
3139       ibook.book1D("Ele_pT_noMTD_Bkg_EB", "Electron pT noMTD - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3140 
3141   meEle_pt_sim_tot_Bkg_EB_ =
3142       ibook.book1D("Ele_pT_sim_tot_Bkg_EB", "Electron pT tot - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3143 
3144   meEle_eta_tot_Bkg_EB_ = ibook.book1D("Ele_eta_tot_Bkg_EB", "Electron eta tot - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3145   meEle_eta_noMTD_Bkg_EB_ =
3146       ibook.book1D("Ele_eta_noMTD_Bkg_EB", "Electron eta noMTD - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3147 
3148   meEle_phi_tot_Bkg_EB_ =
3149       ibook.book1D("Ele_phi_tot_Bkg_EB", "Electron phi tot - Bkg Barrel;#phi;#Counts", 64, -3.2, 3.2);
3150   meEle_phi_noMTD_Bkg_EB_ =
3151       ibook.book1D("Ele_phi_noMTD_Bkg_EB", "Electron phi noMTD - Bkg Barrel;#phi;#Counts", 64, -3.2, 3.2);
3152 
3153   if (optionalPlots_) {
3154     meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EB_ =
3155         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Bkg_EB",
3156                      "Tracks in isolation cone around electron track after basic "
3157                      "cuts with MTD SIM - 4 sigma significance - Bkg Barrel;Number of tracks;Counts",
3158                      20,
3159                      0,
3160                      20);
3161     meEleISO_chIso_MTD_sim_4sigma_Bkg_EB_ =
3162         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Bkg_EB",
3163                      "Track pT sum in isolation cone around electron track after "
3164                      "basic cuts with MTD SIM - 4 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3165                      nbin_2,
3166                      0,
3167                      20);
3168     meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EB_ =
3169         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Bkg_EB",
3170                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
3171                      "sigma significance - Bkg Barrel;Isolation;Counts",
3172                      nbin_1,
3173                      0,
3174                      4);
3175 
3176     meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EB_ =
3177         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Bkg_EB",
3178                      "Tracks in isolation cone around electron track after basic "
3179                      "cuts with MTD SIM - 3 sigma significance - Bkg Barrel;Number of tracks;Counts",
3180                      20,
3181                      0,
3182                      20);
3183     meEleISO_chIso_MTD_sim_3sigma_Bkg_EB_ =
3184         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Bkg_EB",
3185                      "Track pT sum in isolation cone around electron track after "
3186                      "basic cuts with MTD SIM - 3 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3187                      nbin_2,
3188                      0,
3189                      20);
3190     meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EB_ =
3191         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Bkg_EB",
3192                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
3193                      "sigma significance - Bkg Barrel;Isolation;Counts",
3194                      nbin_1,
3195                      0,
3196                      4);
3197 
3198     meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EB_ =
3199         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Bkg_EB",
3200                      "Tracks in isolation cone around electron track after basic "
3201                      "cuts with MTD SIM - 3 sigma significance - Bkg Barrel;Number of tracks;Counts",
3202                      20,
3203                      0,
3204                      20);
3205     meEleISO_chIso_MTD_sim_2sigma_Bkg_EB_ =
3206         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Bkg_EB",
3207                      "Track pT sum in isolation cone around electron track after "
3208                      "basic cuts with MTD SIM - 2 sigma significance - Bkg Barrel;p_{T} (GeV);Counts",
3209                      nbin_2,
3210                      0,
3211                      20);
3212     meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EB_ =
3213         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Bkg_EB",
3214                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
3215                      "sigma significance - Bkg Barrel;Isolation;Counts",
3216                      nbin_1,
3217                      0,
3218                      4);
3219 
3220     meEle_pt_MTD_1_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_1_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3221     meEle_eta_MTD_1_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_1_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3222     meEle_phi_MTD_1_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_1_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3223     meEle_pt_gen_Bkg_EB_ =
3224         ibook.book1D("Ele_pT_gen_Bkg_EB", "Electron pT genInfo - Bkg Barrel;p_{T} (GeV);Counts", 30, 10, 100);
3225     meEle_eta_gen_Bkg_EB_ =
3226         ibook.book1D("Ele_eta_gen_Bkg_EB", "Electron eta genInfo - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3227     meEle_phi_gen_Bkg_EB_ =
3228         ibook.book1D("Ele_phi_gen_Bkg_EB", "Electron phi genInfo - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3229 
3230     meEle_pt_MTD_2_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_2_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3231     meEle_eta_MTD_2_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_2_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3232     meEle_phi_MTD_2_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_2_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3233 
3234     meEle_pt_MTD_3_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_3_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3235     meEle_eta_MTD_3_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_3_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3236     meEle_phi_MTD_3_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_3_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3237 
3238     meEle_pt_MTD_4_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_4_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3239     meEle_eta_MTD_4_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_4_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3240     meEle_phi_MTD_4_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_4_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3241 
3242     meEle_pt_MTD_5_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_5_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3243     meEle_eta_MTD_5_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_5_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3244     meEle_phi_MTD_5_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_5_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3245 
3246     meEle_pt_MTD_6_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_6_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3247     meEle_eta_MTD_6_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_6_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3248     meEle_phi_MTD_6_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_6_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3249 
3250     meEle_pt_MTD_7_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_7_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3251     meEle_eta_MTD_7_Bkg_EB_ = ibook.book1D("Ele_eta_MTD_7_Bkg_EB", "Electron eta MTD;#eta;Counts", 32, 0., 1.6);
3252     meEle_phi_MTD_7_Bkg_EB_ = ibook.book1D("Ele_phi_MTD_7_Bkg_EB", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3253 
3254     meEle_pt_sim_MTD_1_Bkg_EB_ =
3255         ibook.book1D("Ele_pT_sim_MTD_1_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3256     meEle_pt_sim_MTD_2_Bkg_EB_ =
3257         ibook.book1D("Ele_pT_sim_MTD_2_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3258     meEle_pt_sim_MTD_3_Bkg_EB_ =
3259         ibook.book1D("Ele_pT_sim_MTD_3_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3260     meEle_pt_sim_MTD_4_Bkg_EB_ =
3261         ibook.book1D("Ele_pT_sim_MTD_4_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3262     meEle_pt_sim_MTD_5_Bkg_EB_ =
3263         ibook.book1D("Ele_pT_sim_MTD_5_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3264     meEle_pt_sim_MTD_6_Bkg_EB_ =
3265         ibook.book1D("Ele_pT_sim_MTD_6_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3266     meEle_pt_sim_MTD_7_Bkg_EB_ =
3267         ibook.book1D("Ele_pT_sim_MTD_7_Bkg_EB", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3268   }
3269   meEle_pt_MTD_4sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_4sigma_Bkg_EB",
3270                                              "Electron pT MTD - 4 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3271                                              30,
3272                                              10,
3273                                              100);
3274   meEle_eta_MTD_4sigma_Bkg_EB_ = ibook.book1D(
3275       "Ele_eta_MTD_4sigma_Bkg_EB", "Electron eta MTD - 4 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3276   meEle_phi_MTD_4sigma_Bkg_EB_ = ibook.book1D(
3277       "Ele_phi_MTD_4sigma_Bkg_EB", "Electron phi MTD - 4 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3278 
3279   meEle_pt_MTD_3sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_3sigma_Bkg_EB",
3280                                              "Electron pT MTD - 3 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3281                                              30,
3282                                              10,
3283                                              100);
3284   meEle_eta_MTD_3sigma_Bkg_EB_ = ibook.book1D(
3285       "Ele_eta_MTD_3sigma_Bkg_EB", "Electron eta MTD - 3 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3286   meEle_phi_MTD_3sigma_Bkg_EB_ = ibook.book1D(
3287       "Ele_phi_MTD_3sigma_Bkg_EB", "Electron phi MTD - 3 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3288 
3289   meEle_pt_MTD_2sigma_Bkg_EB_ = ibook.book1D("Ele_pT_MTD_2sigma_Bkg_EB",
3290                                              "Electron pT MTD - 2 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3291                                              30,
3292                                              10,
3293                                              100);
3294   meEle_eta_MTD_2sigma_Bkg_EB_ = ibook.book1D(
3295       "Ele_eta_MTD_2sigma_Bkg_EB", "Electron eta MTD - 2 sigma compatibility - Bkg Barrel;#eta;Counts", 32, 0., 1.6);
3296   meEle_phi_MTD_2sigma_Bkg_EB_ = ibook.book1D(
3297       "Ele_phi_MTD_2sigma_Bkg_EB", "Electron phi MTD - 2 sigma compatibility - Bkg Barrel;#phi;Counts", 64, -3.2, 3.2);
3298 
3299   meEleISO_Ntracks_Bkg_EE_ = ibook.book1D(
3300       "Ele_Iso_Ntracks_Bkg_EE",
3301       "Number of tracks in isolation cone around electron track after basic cuts - Bkg Endcap;Number of tracks;Counts",
3302       20,
3303       0,
3304       20);
3305   meEleISO_chIso_Bkg_EE_ = ibook.book1D(
3306       "Ele_chIso_sum_Bkg_EE",
3307       "Track pT sum in isolation cone around electron track after basic cuts - Bkg Endcap;p_{T} (GeV);Counts",
3308       nbin_2,
3309       0,
3310       20);
3311   meEleISO_rel_chIso_Bkg_EE_ = ibook.book1D(
3312       "Ele_rel_chIso_sum_Bkg_EE",
3313       "Track relative pT sum in isolation cone around electron track after basic cuts - Bkg Endcap;Isolation;Counts",
3314       nbin_1,
3315       0,
3316       4);
3317   if (optionalPlots_) {
3318     meEle_pt_sim_MTD_4sigma_Bkg_EB_ =
3319         ibook.book1D("Ele_pT_sim_MTD_4sigma_Bkg_EB",
3320                      "Electron pT MTD SIM - 4 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3321                      30,
3322                      10,
3323                      100);
3324     meEle_pt_sim_MTD_3sigma_Bkg_EB_ =
3325         ibook.book1D("Ele_pT_sim_MTD_3sigma_Bkg_EB",
3326                      "Electron pT MTD SIM - 3 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3327                      30,
3328                      10,
3329                      100);
3330     meEle_pt_sim_MTD_2sigma_Bkg_EB_ =
3331         ibook.book1D("Ele_pT_sim_MTD_2sigma_Bkg_EB",
3332                      "Electron pT MTD SIM - 2 sigma compatibility - Bkg Barrel;p_{T} (GeV);Counts",
3333                      30,
3334                      10,
3335                      100);
3336 
3337     meEleISO_Ntracks_MTD_1_Bkg_EE_ =
3338         ibook.book1D("Ele_Iso_Ntracks_MTD_1_Bkg_EE",
3339                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3340                      20,
3341                      0,
3342                      20);
3343     meEleISO_chIso_MTD_1_Bkg_EE_ = ibook.book1D(
3344         "Ele_chIso_sum_MTD_1_Bkg_EE",
3345         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3346         nbin_2,
3347         0,
3348         20);
3349     meEleISO_rel_chIso_MTD_1_Bkg_EE_ = ibook.book1D(
3350         "Ele_rel_chIso_sum_MTD_1_Bkg_EE",
3351         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3352         nbin_1,
3353         0,
3354         4);
3355 
3356     meEleISO_Ntracks_MTD_2_Bkg_EE_ =
3357         ibook.book1D("Ele_Iso_Ntracks_MTD_2_Bkg_EE",
3358                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3359                      20,
3360                      0,
3361                      20);
3362     meEleISO_chIso_MTD_2_Bkg_EE_ = ibook.book1D(
3363         "Ele_chIso_sum_MTD_2_Bkg_EE",
3364         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3365         nbin_2,
3366         0,
3367         20);
3368     meEleISO_rel_chIso_MTD_2_Bkg_EE_ = ibook.book1D(
3369         "Ele_rel_chIso_sum_MTD_2_Bkg_EE",
3370         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3371         nbin_1,
3372         0,
3373         4);
3374     meEleISO_Ntracks_gen_Bkg_EE_ = ibook.book1D("Ele_Iso_Ntracks_gen_Bkg_EE",
3375                                                 "Tracks in isolation cone around electron track after basic cuts using "
3376                                                 "genInfo - Bkg Endcap;Number of tracks;Counts",
3377                                                 20,
3378                                                 0,
3379                                                 20);
3380     meEleISO_chIso_gen_Bkg_EE_ = ibook.book1D("Ele_chIso_sum_gen_Bkg_EE",
3381                                               "Track pT sum in isolation cone around electron track after basic cuts "
3382                                               "using genInfo - Bkg Endcap;p_{T} (GeV);Counts",
3383                                               nbin_2,
3384                                               0,
3385                                               20);
3386     meEleISO_rel_chIso_gen_Bkg_EE_ = ibook.book1D("Ele_rel_chIso_sum_gen_Bkg_EE",
3387                                                   "Track relative pT sum in isolation cone around electron track after "
3388                                                   "basic cuts using genInfo - Bkg Endcap;Isolation;Counts",
3389                                                   nbin_1,
3390                                                   0,
3391                                                   4);
3392 
3393     meEleISO_Ntracks_MTD_3_Bkg_EE_ =
3394         ibook.book1D("Ele_Iso_Ntracks_MTD_3_Bkg_EE",
3395                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3396                      20,
3397                      0,
3398                      20);
3399     meEleISO_chIso_MTD_3_Bkg_EE_ = ibook.book1D(
3400         "Ele_chIso_sum_MTD_3_Bkg_EE",
3401         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3402         nbin_2,
3403         0,
3404         20);
3405     meEleISO_rel_chIso_MTD_3_Bkg_EE_ = ibook.book1D(
3406         "Ele_rel_chIso_sum_MTD_3_Bkg_EE",
3407         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3408         nbin_1,
3409         0,
3410         4);
3411 
3412     meEleISO_Ntracks_MTD_4_Bkg_EE_ =
3413         ibook.book1D("Ele_Iso_Ntracks_MTD_4_Bkg_EE",
3414                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3415                      20,
3416                      0,
3417                      20);
3418     meEleISO_chIso_MTD_4_Bkg_EE_ = ibook.book1D(
3419         "Ele_chIso_sum_MTD_4_Bkg_EE",
3420         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3421         nbin_2,
3422         0,
3423         20);
3424     meEleISO_rel_chIso_MTD_4_Bkg_EE_ = ibook.book1D(
3425         "Ele_rel_chIso_sum_MTD_4_Bkg_EE",
3426         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3427         nbin_1,
3428         0,
3429         4);
3430 
3431     meEleISO_Ntracks_MTD_5_Bkg_EE_ =
3432         ibook.book1D("Ele_Iso_Ntracks_MTD_5_Bkg_EE",
3433                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3434                      20,
3435                      0,
3436                      20);
3437     meEleISO_chIso_MTD_5_Bkg_EE_ = ibook.book1D(
3438         "Ele_chIso_sum_MTD_5_Bkg_EE",
3439         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3440         nbin_2,
3441         0,
3442         20);
3443     meEleISO_rel_chIso_MTD_5_Bkg_EE_ = ibook.book1D(
3444         "Ele_rel_chIso_sum_MTD_5_Bkg_EE",
3445         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3446         nbin_1,
3447         0,
3448         4);
3449 
3450     meEleISO_Ntracks_MTD_6_Bkg_EE_ =
3451         ibook.book1D("Ele_Iso_Ntracks_MTD_6_Bkg_EE",
3452                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3453                      20,
3454                      0,
3455                      20);
3456     meEleISO_chIso_MTD_6_Bkg_EE_ = ibook.book1D(
3457         "Ele_chIso_sum_MTD_6_Bkg_EE",
3458         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3459         nbin_2,
3460         0,
3461         20);
3462     meEleISO_rel_chIso_MTD_6_Bkg_EE_ = ibook.book1D(
3463         "Ele_rel_chIso_sum_MTD_6_Bkg_EE",
3464         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3465         nbin_1,
3466         0,
3467         4);
3468 
3469     meEleISO_Ntracks_MTD_7_Bkg_EE_ =
3470         ibook.book1D("Ele_Iso_Ntracks_MTD_7_Bkg_EE",
3471                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3472                      20,
3473                      0,
3474                      20);
3475     meEleISO_chIso_MTD_7_Bkg_EE_ = ibook.book1D(
3476         "Ele_chIso_sum_MTD_7_Bkg_EE",
3477         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3478         nbin_2,
3479         0,
3480         20);
3481     meEleISO_rel_chIso_MTD_7_Bkg_EE_ = ibook.book1D(
3482         "Ele_rel_chIso_sum_MTD_7_Bkg_EE",
3483         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3484         nbin_1,
3485         0,
3486         4);
3487 
3488     meEleISO_Ntracks_MTD_sim_1_Bkg_EE_ =
3489         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_1_Bkg_EE",
3490                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3491                      20,
3492                      0,
3493                      20);
3494     meEleISO_chIso_MTD_sim_1_Bkg_EE_ = ibook.book1D(
3495         "Ele_chIso_sum_MTD_sim_1_Bkg_EE",
3496         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3497         nbin_2,
3498         0,
3499         20);
3500     meEleISO_rel_chIso_MTD_sim_1_Bkg_EE_ = ibook.book1D(
3501         "Ele_rel_chIso_sum_MTD_sim_1_Bkg_EE",
3502         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3503         nbin_1,
3504         0,
3505         4);
3506 
3507     meEleISO_Ntracks_MTD_sim_2_Bkg_EE_ =
3508         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2_Bkg_EE",
3509                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3510                      20,
3511                      0,
3512                      20);
3513     meEleISO_chIso_MTD_sim_2_Bkg_EE_ = ibook.book1D(
3514         "Ele_chIso_sum_MTD_sim_2_Bkg_EE",
3515         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3516         nbin_2,
3517         0,
3518         20);
3519     meEleISO_rel_chIso_MTD_sim_2_Bkg_EE_ = ibook.book1D(
3520         "Ele_rel_chIso_sum_MTD_sim_2_Bkg_EE",
3521         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3522         nbin_1,
3523         0,
3524         4);
3525 
3526     meEleISO_Ntracks_MTD_sim_3_Bkg_EE_ =
3527         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3_Bkg_EE",
3528                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3529                      20,
3530                      0,
3531                      20);
3532     meEleISO_chIso_MTD_sim_3_Bkg_EE_ = ibook.book1D(
3533         "Ele_chIso_sum_MTD_sim_3_Bkg_EE",
3534         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3535         nbin_2,
3536         0,
3537         20);
3538     meEleISO_rel_chIso_MTD_sim_3_Bkg_EE_ = ibook.book1D(
3539         "Ele_rel_chIso_sum_MTD_sim_3_Bkg_EE",
3540         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3541         nbin_1,
3542         0,
3543         4);
3544 
3545     meEleISO_Ntracks_MTD_sim_4_Bkg_EE_ =
3546         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4_Bkg_EE",
3547                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3548                      20,
3549                      0,
3550                      20);
3551     meEleISO_chIso_MTD_sim_4_Bkg_EE_ = ibook.book1D(
3552         "Ele_chIso_sum_MTD_sim_4_Bkg_EE",
3553         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3554         nbin_2,
3555         0,
3556         20);
3557     meEleISO_rel_chIso_MTD_sim_4_Bkg_EE_ = ibook.book1D(
3558         "Ele_rel_chIso_sum_MTD_sim_4_Bkg_EE",
3559         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3560         nbin_1,
3561         0,
3562         4);
3563 
3564     meEleISO_Ntracks_MTD_sim_5_Bkg_EE_ =
3565         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_5_Bkg_EE",
3566                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3567                      20,
3568                      0,
3569                      20);
3570     meEleISO_chIso_MTD_sim_5_Bkg_EE_ = ibook.book1D(
3571         "Ele_chIso_sum_MTD_sim_5_Bkg_EE",
3572         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3573         nbin_2,
3574         0,
3575         20);
3576     meEleISO_rel_chIso_MTD_sim_5_Bkg_EE_ = ibook.book1D(
3577         "Ele_rel_chIso_sum_MTD_sim_5_Bkg_EE",
3578         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3579         nbin_1,
3580         0,
3581         4);
3582 
3583     meEleISO_Ntracks_MTD_sim_6_Bkg_EE_ =
3584         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_6_Bkg_EE",
3585                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3586                      20,
3587                      0,
3588                      20);
3589     meEleISO_chIso_MTD_sim_6_Bkg_EE_ = ibook.book1D(
3590         "Ele_chIso_sum_MTD_sim_6_Bkg_EE",
3591         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3592         nbin_2,
3593         0,
3594         20);
3595     meEleISO_rel_chIso_MTD_sim_6_Bkg_EE_ = ibook.book1D(
3596         "Ele_rel_chIso_sum_MTD_sim_6_Bkg_EE",
3597         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3598         nbin_1,
3599         0,
3600         4);
3601 
3602     meEleISO_Ntracks_MTD_sim_7_Bkg_EE_ =
3603         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_7_Bkg_EE",
3604                      "Tracks in isolation cone around electron track after basic cuts with MTD;Number of tracks;Counts",
3605                      20,
3606                      0,
3607                      20);
3608     meEleISO_chIso_MTD_sim_7_Bkg_EE_ = ibook.book1D(
3609         "Ele_chIso_sum_MTD_sim_7_Bkg_EE",
3610         "Track pT sum in isolation cone around electron track after basic cuts with MTD;p_{T} (GeV);Counts",
3611         nbin_2,
3612         0,
3613         20);
3614     meEleISO_rel_chIso_MTD_sim_7_Bkg_EE_ = ibook.book1D(
3615         "Ele_rel_chIso_sum_MTD_sim_7_Bkg_EE",
3616         "Track relative pT sum in isolation cone around electron track after basic cuts with MTD;Isolation;Counts",
3617         nbin_1,
3618         0,
3619         4);
3620   }
3621   meEleISO_Ntracks_MTD_4sigma_Bkg_EE_ =
3622       ibook.book1D("Ele_Iso_Ntracks_MTD_4sigma_Bkg_EE",
3623                    "Tracks in isolation cone around electron track after basic cuts with MTD - 4 sigma compatibility - "
3624                    "Bkg Endcap;Number of tracks;Counts",
3625                    20,
3626                    0,
3627                    20);
3628   meEleISO_chIso_MTD_4sigma_Bkg_EE_ =
3629       ibook.book1D("Ele_chIso_sum_MTD_4sigma_Bkg_EE",
3630                    "Track pT sum in isolation cone around electron track after basic "
3631                    "cuts with MTD - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3632                    nbin_2,
3633                    0,
3634                    20);
3635   meEleISO_rel_chIso_MTD_4sigma_Bkg_EE_ =
3636       ibook.book1D("Ele_rel_chIso_sum_MTD_4sigma_Bkg_EE",
3637                    "Track relative pT sum in isolation cone around electron track "
3638                    "after basic cuts with MTD - 4 sigma compatibility - Bkg Endcap;Isolation;Counts",
3639                    nbin_1,
3640                    0,
3641                    4);
3642 
3643   meEleISO_Ntracks_MTD_3sigma_Bkg_EE_ =
3644       ibook.book1D("Ele_Iso_Ntracks_MTD_3sigma_Bkg_EE",
3645                    "Tracks in isolation cone around electron track after basic cuts with MTD - 3 sigma compatibility - "
3646                    "Bkg Endcap;Number of tracks;Counts",
3647                    20,
3648                    0,
3649                    20);
3650   meEleISO_chIso_MTD_3sigma_Bkg_EE_ =
3651       ibook.book1D("Ele_chIso_sum_MTD_3sigma_Bkg_EE",
3652                    "Track pT sum in isolation cone around electron track after basic "
3653                    "cuts with MTD - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3654                    nbin_2,
3655                    0,
3656                    20);
3657   meEleISO_rel_chIso_MTD_3sigma_Bkg_EE_ =
3658       ibook.book1D("Ele_rel_chIso_sum_MTD_3sigma_Bkg_EE",
3659                    "Track relative pT sum in isolation cone around electron track "
3660                    "after basic cuts with MTD - 3 sigma compatibility - Bkg Endcap;Isolation;Counts",
3661                    nbin_1,
3662                    0,
3663                    4);
3664 
3665   meEleISO_Ntracks_MTD_2sigma_Bkg_EE_ =
3666       ibook.book1D("Ele_Iso_Ntracks_MTD_2sigma_Bkg_EE",
3667                    "Tracks in isolation cone around electron track after basic cuts with MTD - 2 sigma compatibility - "
3668                    "Bkg Endcap;Number of tracks;Counts",
3669                    20,
3670                    0,
3671                    20);
3672   meEleISO_chIso_MTD_2sigma_Bkg_EE_ =
3673       ibook.book1D("Ele_chIso_sum_MTD_2sigma_Bkg_EE",
3674                    "Track pT sum in isolation cone around electron track after basic "
3675                    "cuts with MTD - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3676                    nbin_2,
3677                    0,
3678                    20);
3679   meEleISO_rel_chIso_MTD_2sigma_Bkg_EE_ =
3680       ibook.book1D("Ele_rel_chIso_sum_MTD_2sigma_Bkg_EE",
3681                    "Track relative pT sum in isolation cone around electron track "
3682                    "after basic cuts with MTD - 2 sigma compatibility - Bkg Endcap;Isolation;Counts",
3683                    nbin_1,
3684                    0,
3685                    4);
3686 
3687   meEle_pt_tot_Bkg_EE_ =
3688       ibook.book1D("Ele_pT_tot_Bkg_EE", "Electron pT tot - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3689   meEle_pt_noMTD_Bkg_EE_ =
3690       ibook.book1D("Ele_pT_noMTD_Bkg_EE", "Electron pT noMTD - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3691 
3692   meEle_pt_sim_tot_Bkg_EE_ =
3693       ibook.book1D("Ele_pT_sim_tot_Bkg_EE", "Electron pT tot - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3694 
3695   meEle_eta_tot_Bkg_EE_ = ibook.book1D("Ele_eta_tot_Bkg_EE", "Electron eta tot - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3696   meEle_eta_noMTD_Bkg_EE_ =
3697       ibook.book1D("Ele_eta_noMTD_Bkg_EE", "Electron eta noMTD - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3698 
3699   meEle_phi_tot_Bkg_EE_ =
3700       ibook.book1D("Ele_phi_tot_Bkg_EE", "Electron phi tot - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3701   meEle_phi_noMTD_Bkg_EE_ =
3702       ibook.book1D("Ele_phi_noMTD_Bkg_EE", "Electron phi noMTD - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3703   if (optionalPlots_) {
3704     meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EE_ =
3705         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_4sigma_Bkg_EE",
3706                      "Tracks in isolation cone around electron track after basic "
3707                      "cuts with MTD SIM - 4 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3708                      20,
3709                      0,
3710                      20);
3711     meEleISO_chIso_MTD_sim_4sigma_Bkg_EE_ =
3712         ibook.book1D("Ele_chIso_sum_MTD_sim_4sigma_Bkg_EE",
3713                      "Track pT sum in isolation cone around electron track after "
3714                      "basic cuts with MTD SIM - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3715                      nbin_2,
3716                      0,
3717                      20);
3718     meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EE_ =
3719         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_4sigma_Bkg_EE",
3720                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 4 "
3721                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3722                      nbin_1,
3723                      0,
3724                      4);
3725 
3726     meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EE_ =
3727         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_3sigma_Bkg_EE",
3728                      "Tracks in isolation cone around electron track after basic "
3729                      "cuts with MTD SIM - 3 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3730                      20,
3731                      0,
3732                      20);
3733     meEleISO_chIso_MTD_sim_3sigma_Bkg_EE_ =
3734         ibook.book1D("Ele_chIso_sum_MTD_sim_3sigma_Bkg_EE",
3735                      "Track pT sum in isolation cone around electron track after "
3736                      "basic cuts with MTD SIM - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3737                      nbin_2,
3738                      0,
3739                      20);
3740     meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EE_ =
3741         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_3sigma_Bkg_EE",
3742                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 3 "
3743                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3744                      nbin_1,
3745                      0,
3746                      4);
3747 
3748     meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EE_ =
3749         ibook.book1D("Ele_Iso_Ntracks_MTD_sim_2sigma_Bkg_EE",
3750                      "Tracks in isolation cone around electron track after basic "
3751                      "cuts with MTD SIM - 2 sigma compatibility - Bkg Endcap;Number of tracks;Counts",
3752                      20,
3753                      0,
3754                      20);
3755     meEleISO_chIso_MTD_sim_2sigma_Bkg_EE_ =
3756         ibook.book1D("Ele_chIso_sum_MTD_sim_2sigma_Bkg_EE",
3757                      "Track pT sum in isolation cone around electron track after "
3758                      "basic cuts with MTD SIM - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3759                      nbin_2,
3760                      0,
3761                      20);
3762     meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EE_ =
3763         ibook.book1D("Ele_rel_chIso_sum_MTD_sim_2sigma_Bkg_EE",
3764                      "Track relative pT sum in isolation cone around electron track after basic cuts with MTD SIM - 2 "
3765                      "sigma compatibility - Bkg Endcap;Isolation;Counts",
3766                      nbin_1,
3767                      0,
3768                      4);
3769 
3770     meEle_pt_MTD_1_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_1_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3771     meEle_eta_MTD_1_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_1_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3772     meEle_phi_MTD_1_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_1_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3773     meEle_pt_gen_Bkg_EE_ =
3774         ibook.book1D("Ele_pT_gen_Bkg_EE", "Electron pT genInfo - Bkg Endcap;p_{T} (GeV);Counts", 30, 10, 100);
3775     meEle_eta_gen_Bkg_EE_ =
3776         ibook.book1D("Ele_eta_gen_Bkg_EE", "Electron eta genInfo - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3777     meEle_phi_gen_Bkg_EE_ =
3778         ibook.book1D("Ele_phi_gen_Bkg_EE", "Electron phi genInfo - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3779 
3780     meEle_pt_MTD_2_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_2_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3781     meEle_eta_MTD_2_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_2_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3782     meEle_phi_MTD_2_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_2_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3783 
3784     meEle_pt_MTD_3_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_3_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3785     meEle_eta_MTD_3_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_3_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3786     meEle_phi_MTD_3_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_3_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3787 
3788     meEle_pt_MTD_4_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_4_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3789     meEle_eta_MTD_4_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_4_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3790     meEle_phi_MTD_4_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_4_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3791 
3792     meEle_pt_MTD_5_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_5_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3793     meEle_eta_MTD_5_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_5_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3794     meEle_phi_MTD_5_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_5_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3795 
3796     meEle_pt_MTD_6_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_6_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3797     meEle_eta_MTD_6_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_6_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3798     meEle_phi_MTD_6_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_6_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3799 
3800     meEle_pt_MTD_7_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_7_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3801     meEle_eta_MTD_7_Bkg_EE_ = ibook.book1D("Ele_eta_MTD_7_Bkg_EE", "Electron eta MTD;#eta;Counts", 32, 1.6, 3.2);
3802     meEle_phi_MTD_7_Bkg_EE_ = ibook.book1D("Ele_phi_MTD_7_Bkg_EE", "Electron phi MTD;#phi;Counts", 64, -3.2, 3.2);
3803 
3804     meEle_pt_sim_MTD_1_Bkg_EE_ =
3805         ibook.book1D("Ele_pT_sim_MTD_1_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3806     meEle_pt_sim_MTD_2_Bkg_EE_ =
3807         ibook.book1D("Ele_pT_sim_MTD_2_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3808     meEle_pt_sim_MTD_3_Bkg_EE_ =
3809         ibook.book1D("Ele_pT_sim_MTD_3_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3810     meEle_pt_sim_MTD_4_Bkg_EE_ =
3811         ibook.book1D("Ele_pT_sim_MTD_4_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3812     meEle_pt_sim_MTD_5_Bkg_EE_ =
3813         ibook.book1D("Ele_pT_sim_MTD_5_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3814     meEle_pt_sim_MTD_6_Bkg_EE_ =
3815         ibook.book1D("Ele_pT_sim_MTD_6_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3816     meEle_pt_sim_MTD_7_Bkg_EE_ =
3817         ibook.book1D("Ele_pT_sim_MTD_7_Bkg_EE", "Electron pT MTD;p_{T} (GeV);Counts", 30, 10, 100);
3818   }
3819 
3820   meEle_pt_MTD_4sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_4sigma_Bkg_EE",
3821                                              "Electron pT MTD - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3822                                              30,
3823                                              10,
3824                                              100);
3825   meEle_eta_MTD_4sigma_Bkg_EE_ = ibook.book1D(
3826       "Ele_eta_MTD_4sigma_Bkg_EE", "Electron eta MTD - 4 sigma compatibility - Bkg Endcapi;#eta;Counts", 32, 1.6, 3.2);
3827   meEle_phi_MTD_4sigma_Bkg_EE_ = ibook.book1D(
3828       "Ele_phi_MTD_4sigma_Bkg_EE", "Electron phi MTD - 4 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3829 
3830   meEle_pt_MTD_3sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_3sigma_Bkg_EE",
3831                                              "Electron pT MTD - 3 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3832                                              30,
3833                                              10,
3834                                              100);
3835   meEle_eta_MTD_3sigma_Bkg_EE_ = ibook.book1D(
3836       "Ele_eta_MTD_3sigma_Bkg_EE", "Electron eta MTD - 3 sigma compatibility - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3837   meEle_phi_MTD_3sigma_Bkg_EE_ = ibook.book1D(
3838       "Ele_phi_MTD_3sigma_Bkg_EE", "Electron phi MTD - 3 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3839 
3840   meEle_pt_MTD_2sigma_Bkg_EE_ = ibook.book1D("Ele_pT_MTD_2sigma_Bkg_EE",
3841                                              "Electron pT MTD - 2 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3842                                              30,
3843                                              10,
3844                                              100);
3845   meEle_eta_MTD_2sigma_Bkg_EE_ = ibook.book1D(
3846       "Ele_eta_MTD_2sigma_Bkg_EE", "Electron eta MTD - 2 sigma compatibility - Bkg Endcap;#eta;Counts", 32, 1.6, 3.2);
3847   meEle_phi_MTD_2sigma_Bkg_EE_ = ibook.book1D(
3848       "Ele_phi_MTD_2sigma_Bkg_EE", "Electron phi MTD - 2 sigma compatibility - Bkg Endcap;#phi;Counts", 64, -3.2, 3.2);
3849 
3850   if (optionalPlots_) {
3851     meEle_pt_sim_MTD_4sigma_Bkg_EE_ =
3852         ibook.book1D("Ele_pT_sim_MTD_4sigma_Bkg_EE",
3853                      "Electron pT MTD SIM - 4 sigma compatibility - Bkg Endcap;p_{T} (GeV);Counts",
3854                      30,
3855                      10,
3856                      100);
3857     meEle_pt_sim_MTD_3sigma_Bkg_EE_ =
3858         ibook.book1D("Ele_pT_sim_MTD_3sigma_Bkg_EE",
3859                      "Electron pT MTD SIM - 3 sigma compatibility - Bkg Endcap;#eta;Counts",
3860                      30,
3861                      10,
3862                      100);
3863     meEle_pt_sim_MTD_2sigma_Bkg_EE_ =
3864         ibook.book1D("Ele_pT_sim_MTD_2sigma_Bkg_EE",
3865                      "Electron pT MTD SIM - 2 sigma compatibility - Bkg Endcap;#phi;Counts",
3866                      30,
3867                      10,
3868                      100);
3869   }
3870 
3871   // defining vectors for more efficient hist filling
3872   // Promt part
3873   if (optionalPlots_) {
3874     Ntracks_EB_list_Sig = {meEleISO_Ntracks_MTD_1_Sig_EB_,
3875                            meEleISO_Ntracks_MTD_2_Sig_EB_,
3876                            meEleISO_Ntracks_MTD_3_Sig_EB_,
3877                            meEleISO_Ntracks_MTD_4_Sig_EB_,
3878                            meEleISO_Ntracks_MTD_5_Sig_EB_,
3879                            meEleISO_Ntracks_MTD_6_Sig_EB_,
3880                            meEleISO_Ntracks_MTD_7_Sig_EB_};
3881     ch_iso_EB_list_Sig = {meEleISO_chIso_MTD_1_Sig_EB_,
3882                           meEleISO_chIso_MTD_2_Sig_EB_,
3883                           meEleISO_chIso_MTD_3_Sig_EB_,
3884                           meEleISO_chIso_MTD_4_Sig_EB_,
3885                           meEleISO_chIso_MTD_5_Sig_EB_,
3886                           meEleISO_chIso_MTD_6_Sig_EB_,
3887                           meEleISO_chIso_MTD_7_Sig_EB_};
3888     rel_ch_iso_EB_list_Sig = {meEleISO_rel_chIso_MTD_1_Sig_EB_,
3889                               meEleISO_rel_chIso_MTD_2_Sig_EB_,
3890                               meEleISO_rel_chIso_MTD_3_Sig_EB_,
3891                               meEleISO_rel_chIso_MTD_4_Sig_EB_,
3892                               meEleISO_rel_chIso_MTD_5_Sig_EB_,
3893                               meEleISO_rel_chIso_MTD_6_Sig_EB_,
3894                               meEleISO_rel_chIso_MTD_7_Sig_EB_};
3895   }
3896   Ntracks_EB_list_Significance_Sig = {
3897       meEleISO_Ntracks_MTD_4sigma_Sig_EB_, meEleISO_Ntracks_MTD_3sigma_Sig_EB_, meEleISO_Ntracks_MTD_2sigma_Sig_EB_};
3898   ch_iso_EB_list_Significance_Sig = {
3899       meEleISO_chIso_MTD_4sigma_Sig_EB_, meEleISO_chIso_MTD_3sigma_Sig_EB_, meEleISO_chIso_MTD_2sigma_Sig_EB_};
3900   rel_ch_iso_EB_list_Significance_Sig = {meEleISO_rel_chIso_MTD_4sigma_Sig_EB_,
3901                                          meEleISO_rel_chIso_MTD_3sigma_Sig_EB_,
3902                                          meEleISO_rel_chIso_MTD_2sigma_Sig_EB_};
3903 
3904   if (optionalPlots_) {
3905     Ntracks_EE_list_Sig = {meEleISO_Ntracks_MTD_1_Sig_EE_,
3906                            meEleISO_Ntracks_MTD_2_Sig_EE_,
3907                            meEleISO_Ntracks_MTD_3_Sig_EE_,
3908                            meEleISO_Ntracks_MTD_4_Sig_EE_,
3909                            meEleISO_Ntracks_MTD_5_Sig_EE_,
3910                            meEleISO_Ntracks_MTD_6_Sig_EE_,
3911                            meEleISO_Ntracks_MTD_7_Sig_EE_};
3912     ch_iso_EE_list_Sig = {meEleISO_chIso_MTD_1_Sig_EE_,
3913                           meEleISO_chIso_MTD_2_Sig_EE_,
3914                           meEleISO_chIso_MTD_3_Sig_EE_,
3915                           meEleISO_chIso_MTD_4_Sig_EE_,
3916                           meEleISO_chIso_MTD_5_Sig_EE_,
3917                           meEleISO_chIso_MTD_6_Sig_EE_,
3918                           meEleISO_chIso_MTD_7_Sig_EE_};
3919     rel_ch_iso_EE_list_Sig = {meEleISO_rel_chIso_MTD_1_Sig_EE_,
3920                               meEleISO_rel_chIso_MTD_2_Sig_EE_,
3921                               meEleISO_rel_chIso_MTD_3_Sig_EE_,
3922                               meEleISO_rel_chIso_MTD_4_Sig_EE_,
3923                               meEleISO_rel_chIso_MTD_5_Sig_EE_,
3924                               meEleISO_rel_chIso_MTD_6_Sig_EE_,
3925                               meEleISO_rel_chIso_MTD_7_Sig_EE_};
3926   }
3927   Ntracks_EE_list_Significance_Sig = {
3928       meEleISO_Ntracks_MTD_4sigma_Sig_EE_, meEleISO_Ntracks_MTD_3sigma_Sig_EE_, meEleISO_Ntracks_MTD_2sigma_Sig_EE_};
3929   ch_iso_EE_list_Significance_Sig = {
3930       meEleISO_chIso_MTD_4sigma_Sig_EE_, meEleISO_chIso_MTD_3sigma_Sig_EE_, meEleISO_chIso_MTD_2sigma_Sig_EE_};
3931   rel_ch_iso_EE_list_Significance_Sig = {meEleISO_rel_chIso_MTD_4sigma_Sig_EE_,
3932                                          meEleISO_rel_chIso_MTD_3sigma_Sig_EE_,
3933                                          meEleISO_rel_chIso_MTD_2sigma_Sig_EE_};
3934 
3935   if (optionalPlots_) {
3936     Ele_pT_MTD_EB_list_Sig = {meEle_pt_MTD_1_Sig_EB_,
3937                               meEle_pt_MTD_2_Sig_EB_,
3938                               meEle_pt_MTD_3_Sig_EB_,
3939                               meEle_pt_MTD_4_Sig_EB_,
3940                               meEle_pt_MTD_5_Sig_EB_,
3941                               meEle_pt_MTD_6_Sig_EB_,
3942                               meEle_pt_MTD_7_Sig_EB_};
3943     Ele_eta_MTD_EB_list_Sig = {meEle_eta_MTD_1_Sig_EB_,
3944                                meEle_eta_MTD_2_Sig_EB_,
3945                                meEle_eta_MTD_3_Sig_EB_,
3946                                meEle_eta_MTD_4_Sig_EB_,
3947                                meEle_eta_MTD_5_Sig_EB_,
3948                                meEle_eta_MTD_6_Sig_EB_,
3949                                meEle_eta_MTD_7_Sig_EB_};
3950     Ele_phi_MTD_EB_list_Sig = {meEle_phi_MTD_1_Sig_EB_,
3951                                meEle_phi_MTD_2_Sig_EB_,
3952                                meEle_phi_MTD_3_Sig_EB_,
3953                                meEle_phi_MTD_4_Sig_EB_,
3954                                meEle_phi_MTD_5_Sig_EB_,
3955                                meEle_phi_MTD_6_Sig_EB_,
3956                                meEle_phi_MTD_7_Sig_EB_};
3957   }
3958 
3959   Ele_pT_MTD_EB_list_Significance_Sig = {
3960       meEle_pt_MTD_4sigma_Sig_EB_, meEle_pt_MTD_3sigma_Sig_EB_, meEle_pt_MTD_2sigma_Sig_EB_};
3961   Ele_eta_MTD_EB_list_Significance_Sig = {
3962       meEle_eta_MTD_4sigma_Sig_EB_, meEle_eta_MTD_3sigma_Sig_EB_, meEle_eta_MTD_2sigma_Sig_EB_};
3963   Ele_phi_MTD_EB_list_Significance_Sig = {
3964       meEle_phi_MTD_4sigma_Sig_EB_, meEle_phi_MTD_3sigma_Sig_EB_, meEle_phi_MTD_2sigma_Sig_EB_};
3965 
3966   if (optionalPlots_) {
3967     Ele_pT_MTD_EE_list_Sig = {meEle_pt_MTD_1_Sig_EE_,
3968                               meEle_pt_MTD_2_Sig_EE_,
3969                               meEle_pt_MTD_3_Sig_EE_,
3970                               meEle_pt_MTD_4_Sig_EE_,
3971                               meEle_pt_MTD_5_Sig_EE_,
3972                               meEle_pt_MTD_6_Sig_EE_,
3973                               meEle_pt_MTD_7_Sig_EE_};
3974     Ele_eta_MTD_EE_list_Sig = {meEle_eta_MTD_1_Sig_EE_,
3975                                meEle_eta_MTD_2_Sig_EE_,
3976                                meEle_eta_MTD_3_Sig_EE_,
3977                                meEle_eta_MTD_4_Sig_EE_,
3978                                meEle_eta_MTD_5_Sig_EE_,
3979                                meEle_eta_MTD_6_Sig_EE_,
3980                                meEle_eta_MTD_7_Sig_EE_};
3981     Ele_phi_MTD_EE_list_Sig = {meEle_phi_MTD_1_Sig_EE_,
3982                                meEle_phi_MTD_2_Sig_EE_,
3983                                meEle_phi_MTD_3_Sig_EE_,
3984                                meEle_phi_MTD_4_Sig_EE_,
3985                                meEle_phi_MTD_5_Sig_EE_,
3986                                meEle_phi_MTD_6_Sig_EE_,
3987                                meEle_phi_MTD_7_Sig_EE_};
3988   }
3989   Ele_pT_MTD_EE_list_Significance_Sig = {
3990       meEle_pt_MTD_4sigma_Sig_EE_, meEle_pt_MTD_3sigma_Sig_EE_, meEle_pt_MTD_2sigma_Sig_EE_};
3991   Ele_eta_MTD_EE_list_Significance_Sig = {
3992       meEle_eta_MTD_4sigma_Sig_EE_, meEle_eta_MTD_3sigma_Sig_EE_, meEle_eta_MTD_2sigma_Sig_EE_};
3993   Ele_phi_MTD_EE_list_Significance_Sig = {
3994       meEle_phi_MTD_4sigma_Sig_EE_, meEle_phi_MTD_3sigma_Sig_EE_, meEle_phi_MTD_2sigma_Sig_EE_};
3995 
3996   // For SIM CASE
3997   if (optionalPlots_) {
3998     Ntracks_sim_EB_list_Sig = {meEleISO_Ntracks_MTD_sim_1_Sig_EB_,
3999                                meEleISO_Ntracks_MTD_sim_2_Sig_EB_,
4000                                meEleISO_Ntracks_MTD_sim_3_Sig_EB_,
4001                                meEleISO_Ntracks_MTD_sim_4_Sig_EB_,
4002                                meEleISO_Ntracks_MTD_sim_5_Sig_EB_,
4003                                meEleISO_Ntracks_MTD_sim_6_Sig_EB_,
4004                                meEleISO_Ntracks_MTD_sim_7_Sig_EB_};
4005     ch_iso_sim_EB_list_Sig = {meEleISO_chIso_MTD_sim_1_Sig_EB_,
4006                               meEleISO_chIso_MTD_sim_2_Sig_EB_,
4007                               meEleISO_chIso_MTD_sim_3_Sig_EB_,
4008                               meEleISO_chIso_MTD_sim_4_Sig_EB_,
4009                               meEleISO_chIso_MTD_sim_5_Sig_EB_,
4010                               meEleISO_chIso_MTD_sim_6_Sig_EB_,
4011                               meEleISO_chIso_MTD_sim_7_Sig_EB_};
4012     rel_ch_iso_sim_EB_list_Sig = {meEleISO_rel_chIso_MTD_sim_1_Sig_EB_,
4013                                   meEleISO_rel_chIso_MTD_sim_2_Sig_EB_,
4014                                   meEleISO_rel_chIso_MTD_sim_3_Sig_EB_,
4015                                   meEleISO_rel_chIso_MTD_sim_4_Sig_EB_,
4016                                   meEleISO_rel_chIso_MTD_sim_5_Sig_EB_,
4017                                   meEleISO_rel_chIso_MTD_sim_6_Sig_EB_,
4018                                   meEleISO_rel_chIso_MTD_sim_7_Sig_EB_};
4019 
4020     Ntracks_sim_EB_list_Significance_Sig = {meEleISO_Ntracks_MTD_sim_4sigma_Sig_EB_,
4021                                             meEleISO_Ntracks_MTD_sim_3sigma_Sig_EB_,
4022                                             meEleISO_Ntracks_MTD_sim_2sigma_Sig_EB_};
4023     ch_iso_sim_EB_list_Significance_Sig = {meEleISO_chIso_MTD_sim_4sigma_Sig_EB_,
4024                                            meEleISO_chIso_MTD_sim_3sigma_Sig_EB_,
4025                                            meEleISO_chIso_MTD_sim_2sigma_Sig_EB_};
4026     rel_ch_iso_sim_EB_list_Significance_Sig = {meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EB_,
4027                                                meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EB_,
4028                                                meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EB_};
4029 
4030     Ntracks_sim_EE_list_Sig = {meEleISO_Ntracks_MTD_sim_1_Sig_EE_,
4031                                meEleISO_Ntracks_MTD_sim_2_Sig_EE_,
4032                                meEleISO_Ntracks_MTD_sim_3_Sig_EE_,
4033                                meEleISO_Ntracks_MTD_sim_4_Sig_EE_,
4034                                meEleISO_Ntracks_MTD_sim_5_Sig_EE_,
4035                                meEleISO_Ntracks_MTD_sim_6_Sig_EE_,
4036                                meEleISO_Ntracks_MTD_sim_7_Sig_EE_};
4037     ch_iso_sim_EE_list_Sig = {meEleISO_chIso_MTD_sim_1_Sig_EE_,
4038                               meEleISO_chIso_MTD_sim_2_Sig_EE_,
4039                               meEleISO_chIso_MTD_sim_3_Sig_EE_,
4040                               meEleISO_chIso_MTD_sim_4_Sig_EE_,
4041                               meEleISO_chIso_MTD_sim_5_Sig_EE_,
4042                               meEleISO_chIso_MTD_sim_6_Sig_EE_,
4043                               meEleISO_chIso_MTD_sim_7_Sig_EE_};
4044     rel_ch_iso_sim_EE_list_Sig = {meEleISO_rel_chIso_MTD_sim_1_Sig_EE_,
4045                                   meEleISO_rel_chIso_MTD_sim_2_Sig_EE_,
4046                                   meEleISO_rel_chIso_MTD_sim_3_Sig_EE_,
4047                                   meEleISO_rel_chIso_MTD_sim_4_Sig_EE_,
4048                                   meEleISO_rel_chIso_MTD_sim_5_Sig_EE_,
4049                                   meEleISO_rel_chIso_MTD_sim_6_Sig_EE_,
4050                                   meEleISO_rel_chIso_MTD_sim_7_Sig_EE_};
4051 
4052     Ntracks_sim_EE_list_Significance_Sig = {meEleISO_Ntracks_MTD_sim_4sigma_Sig_EE_,
4053                                             meEleISO_Ntracks_MTD_sim_3sigma_Sig_EE_,
4054                                             meEleISO_Ntracks_MTD_sim_2sigma_Sig_EE_};
4055     ch_iso_sim_EE_list_Significance_Sig = {meEleISO_chIso_MTD_sim_4sigma_Sig_EE_,
4056                                            meEleISO_chIso_MTD_sim_3sigma_Sig_EE_,
4057                                            meEleISO_chIso_MTD_sim_2sigma_Sig_EE_};
4058     rel_ch_iso_sim_EE_list_Significance_Sig = {meEleISO_rel_chIso_MTD_sim_4sigma_Sig_EE_,
4059                                                meEleISO_rel_chIso_MTD_sim_3sigma_Sig_EE_,
4060                                                meEleISO_rel_chIso_MTD_sim_2sigma_Sig_EE_};
4061 
4062     Ele_pT_sim_MTD_EB_list_Sig = {meEle_pt_sim_MTD_1_Sig_EB_,
4063                                   meEle_pt_sim_MTD_2_Sig_EB_,
4064                                   meEle_pt_sim_MTD_3_Sig_EB_,
4065                                   meEle_pt_sim_MTD_4_Sig_EB_,
4066                                   meEle_pt_sim_MTD_5_Sig_EB_,
4067                                   meEle_pt_sim_MTD_6_Sig_EB_,
4068                                   meEle_pt_sim_MTD_7_Sig_EB_};
4069 
4070     Ele_pT_sim_MTD_EB_list_Significance_Sig = {
4071         meEle_pt_sim_MTD_4sigma_Sig_EB_, meEle_pt_sim_MTD_3sigma_Sig_EB_, meEle_pt_sim_MTD_2sigma_Sig_EB_};
4072 
4073     Ele_pT_sim_MTD_EE_list_Sig = {meEle_pt_sim_MTD_1_Sig_EE_,
4074                                   meEle_pt_sim_MTD_2_Sig_EE_,
4075                                   meEle_pt_sim_MTD_3_Sig_EE_,
4076                                   meEle_pt_sim_MTD_4_Sig_EE_,
4077                                   meEle_pt_sim_MTD_5_Sig_EE_,
4078                                   meEle_pt_sim_MTD_6_Sig_EE_,
4079                                   meEle_pt_sim_MTD_7_Sig_EE_};
4080     Ele_pT_sim_MTD_EE_list_Significance_Sig = {
4081         meEle_pt_sim_MTD_4sigma_Sig_EE_, meEle_pt_sim_MTD_3sigma_Sig_EE_, meEle_pt_sim_MTD_2sigma_Sig_EE_};
4082   }
4083 
4084   // Non-promt part
4085   if (optionalPlots_) {
4086     Ntracks_EB_list_Bkg = {meEleISO_Ntracks_MTD_1_Bkg_EB_,
4087                            meEleISO_Ntracks_MTD_2_Bkg_EB_,
4088                            meEleISO_Ntracks_MTD_3_Bkg_EB_,
4089                            meEleISO_Ntracks_MTD_4_Bkg_EB_,
4090                            meEleISO_Ntracks_MTD_5_Bkg_EB_,
4091                            meEleISO_Ntracks_MTD_6_Bkg_EB_,
4092                            meEleISO_Ntracks_MTD_7_Bkg_EB_};
4093     ch_iso_EB_list_Bkg = {meEleISO_chIso_MTD_1_Bkg_EB_,
4094                           meEleISO_chIso_MTD_2_Bkg_EB_,
4095                           meEleISO_chIso_MTD_3_Bkg_EB_,
4096                           meEleISO_chIso_MTD_4_Bkg_EB_,
4097                           meEleISO_chIso_MTD_5_Bkg_EB_,
4098                           meEleISO_chIso_MTD_6_Bkg_EB_,
4099                           meEleISO_chIso_MTD_7_Bkg_EB_};
4100     rel_ch_iso_EB_list_Bkg = {meEleISO_rel_chIso_MTD_1_Bkg_EB_,
4101                               meEleISO_rel_chIso_MTD_2_Bkg_EB_,
4102                               meEleISO_rel_chIso_MTD_3_Bkg_EB_,
4103                               meEleISO_rel_chIso_MTD_4_Bkg_EB_,
4104                               meEleISO_rel_chIso_MTD_5_Bkg_EB_,
4105                               meEleISO_rel_chIso_MTD_6_Bkg_EB_,
4106                               meEleISO_rel_chIso_MTD_7_Bkg_EB_};
4107   }
4108   Ntracks_EB_list_Significance_Bkg = {
4109       meEleISO_Ntracks_MTD_4sigma_Bkg_EB_, meEleISO_Ntracks_MTD_3sigma_Bkg_EB_, meEleISO_Ntracks_MTD_2sigma_Bkg_EB_};
4110   ch_iso_EB_list_Significance_Bkg = {
4111       meEleISO_chIso_MTD_4sigma_Bkg_EB_, meEleISO_chIso_MTD_3sigma_Bkg_EB_, meEleISO_chIso_MTD_2sigma_Bkg_EB_};
4112   rel_ch_iso_EB_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_4sigma_Bkg_EB_,
4113                                          meEleISO_rel_chIso_MTD_3sigma_Bkg_EB_,
4114                                          meEleISO_rel_chIso_MTD_2sigma_Bkg_EB_};
4115 
4116   if (optionalPlots_) {
4117     Ntracks_EE_list_Bkg = {meEleISO_Ntracks_MTD_1_Bkg_EE_,
4118                            meEleISO_Ntracks_MTD_2_Bkg_EE_,
4119                            meEleISO_Ntracks_MTD_3_Bkg_EE_,
4120                            meEleISO_Ntracks_MTD_4_Bkg_EE_,
4121                            meEleISO_Ntracks_MTD_5_Bkg_EE_,
4122                            meEleISO_Ntracks_MTD_6_Bkg_EE_,
4123                            meEleISO_Ntracks_MTD_7_Bkg_EE_};
4124     ch_iso_EE_list_Bkg = {meEleISO_chIso_MTD_1_Bkg_EE_,
4125                           meEleISO_chIso_MTD_2_Bkg_EE_,
4126                           meEleISO_chIso_MTD_3_Bkg_EE_,
4127                           meEleISO_chIso_MTD_4_Bkg_EE_,
4128                           meEleISO_chIso_MTD_5_Bkg_EE_,
4129                           meEleISO_chIso_MTD_6_Bkg_EE_,
4130                           meEleISO_chIso_MTD_7_Bkg_EE_};
4131     rel_ch_iso_EE_list_Bkg = {meEleISO_rel_chIso_MTD_1_Bkg_EE_,
4132                               meEleISO_rel_chIso_MTD_2_Bkg_EE_,
4133                               meEleISO_rel_chIso_MTD_3_Bkg_EE_,
4134                               meEleISO_rel_chIso_MTD_4_Bkg_EE_,
4135                               meEleISO_rel_chIso_MTD_5_Bkg_EE_,
4136                               meEleISO_rel_chIso_MTD_6_Bkg_EE_,
4137                               meEleISO_rel_chIso_MTD_7_Bkg_EE_};
4138   }
4139   Ntracks_EE_list_Significance_Bkg = {
4140       meEleISO_Ntracks_MTD_4sigma_Bkg_EE_, meEleISO_Ntracks_MTD_3sigma_Bkg_EE_, meEleISO_Ntracks_MTD_2sigma_Bkg_EE_};
4141   ch_iso_EE_list_Significance_Bkg = {
4142       meEleISO_chIso_MTD_4sigma_Bkg_EE_, meEleISO_chIso_MTD_3sigma_Bkg_EE_, meEleISO_chIso_MTD_2sigma_Bkg_EE_};
4143   rel_ch_iso_EE_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_4sigma_Bkg_EE_,
4144                                          meEleISO_rel_chIso_MTD_3sigma_Bkg_EE_,
4145                                          meEleISO_rel_chIso_MTD_2sigma_Bkg_EE_};
4146   if (optionalPlots_) {
4147     Ele_pT_MTD_EB_list_Bkg = {meEle_pt_MTD_1_Bkg_EB_,
4148                               meEle_pt_MTD_2_Bkg_EB_,
4149                               meEle_pt_MTD_3_Bkg_EB_,
4150                               meEle_pt_MTD_4_Bkg_EB_,
4151                               meEle_pt_MTD_5_Bkg_EB_,
4152                               meEle_pt_MTD_6_Bkg_EB_,
4153                               meEle_pt_MTD_7_Bkg_EB_};
4154     Ele_eta_MTD_EB_list_Bkg = {meEle_eta_MTD_1_Bkg_EB_,
4155                                meEle_eta_MTD_2_Bkg_EB_,
4156                                meEle_eta_MTD_3_Bkg_EB_,
4157                                meEle_eta_MTD_4_Bkg_EB_,
4158                                meEle_eta_MTD_5_Bkg_EB_,
4159                                meEle_eta_MTD_6_Bkg_EB_,
4160                                meEle_eta_MTD_7_Bkg_EB_};
4161     Ele_phi_MTD_EB_list_Bkg = {meEle_phi_MTD_1_Bkg_EB_,
4162                                meEle_phi_MTD_2_Bkg_EB_,
4163                                meEle_phi_MTD_3_Bkg_EB_,
4164                                meEle_phi_MTD_4_Bkg_EB_,
4165                                meEle_phi_MTD_5_Bkg_EB_,
4166                                meEle_phi_MTD_6_Bkg_EB_,
4167                                meEle_phi_MTD_7_Bkg_EB_};
4168   }
4169   Ele_pT_MTD_EB_list_Significance_Bkg = {
4170       meEle_pt_MTD_4sigma_Bkg_EB_, meEle_pt_MTD_3sigma_Bkg_EB_, meEle_pt_MTD_2sigma_Bkg_EB_};
4171   Ele_eta_MTD_EB_list_Significance_Bkg = {
4172       meEle_eta_MTD_4sigma_Bkg_EB_, meEle_eta_MTD_3sigma_Bkg_EB_, meEle_eta_MTD_2sigma_Bkg_EB_};
4173   Ele_phi_MTD_EB_list_Significance_Bkg = {
4174       meEle_phi_MTD_4sigma_Bkg_EB_, meEle_phi_MTD_3sigma_Bkg_EB_, meEle_phi_MTD_2sigma_Bkg_EB_};
4175 
4176   if (optionalPlots_) {
4177     Ele_pT_MTD_EE_list_Bkg = {meEle_pt_MTD_1_Bkg_EE_,
4178                               meEle_pt_MTD_2_Bkg_EE_,
4179                               meEle_pt_MTD_3_Bkg_EE_,
4180                               meEle_pt_MTD_4_Bkg_EE_,
4181                               meEle_pt_MTD_5_Bkg_EE_,
4182                               meEle_pt_MTD_6_Bkg_EE_,
4183                               meEle_pt_MTD_7_Bkg_EE_};
4184     Ele_eta_MTD_EE_list_Bkg = {meEle_eta_MTD_1_Bkg_EE_,
4185                                meEle_eta_MTD_2_Bkg_EE_,
4186                                meEle_eta_MTD_3_Bkg_EE_,
4187                                meEle_eta_MTD_4_Bkg_EE_,
4188                                meEle_eta_MTD_5_Bkg_EE_,
4189                                meEle_eta_MTD_6_Bkg_EE_,
4190                                meEle_eta_MTD_7_Bkg_EE_};
4191     Ele_phi_MTD_EE_list_Bkg = {meEle_phi_MTD_1_Bkg_EE_,
4192                                meEle_phi_MTD_2_Bkg_EE_,
4193                                meEle_phi_MTD_3_Bkg_EE_,
4194                                meEle_phi_MTD_4_Bkg_EE_,
4195                                meEle_phi_MTD_5_Bkg_EE_,
4196                                meEle_phi_MTD_6_Bkg_EE_,
4197                                meEle_phi_MTD_7_Bkg_EE_};
4198   }
4199   Ele_pT_MTD_EE_list_Significance_Bkg = {
4200       meEle_pt_MTD_4sigma_Bkg_EE_, meEle_pt_MTD_3sigma_Bkg_EE_, meEle_pt_MTD_2sigma_Bkg_EE_};
4201   Ele_eta_MTD_EE_list_Significance_Bkg = {
4202       meEle_eta_MTD_4sigma_Bkg_EE_, meEle_eta_MTD_3sigma_Bkg_EE_, meEle_eta_MTD_2sigma_Bkg_EE_};
4203   Ele_phi_MTD_EE_list_Significance_Bkg = {
4204       meEle_phi_MTD_4sigma_Bkg_EE_, meEle_phi_MTD_3sigma_Bkg_EE_, meEle_phi_MTD_2sigma_Bkg_EE_};
4205 
4206   // SIM CASE
4207   if (optionalPlots_) {
4208     Ntracks_sim_EB_list_Bkg = {meEleISO_Ntracks_MTD_sim_1_Bkg_EB_,
4209                                meEleISO_Ntracks_MTD_sim_2_Bkg_EB_,
4210                                meEleISO_Ntracks_MTD_sim_3_Bkg_EB_,
4211                                meEleISO_Ntracks_MTD_sim_4_Bkg_EB_,
4212                                meEleISO_Ntracks_MTD_sim_5_Bkg_EB_,
4213                                meEleISO_Ntracks_MTD_sim_6_Bkg_EB_,
4214                                meEleISO_Ntracks_MTD_sim_7_Bkg_EB_};
4215     ch_iso_sim_EB_list_Bkg = {meEleISO_chIso_MTD_sim_1_Bkg_EB_,
4216                               meEleISO_chIso_MTD_sim_2_Bkg_EB_,
4217                               meEleISO_chIso_MTD_sim_3_Bkg_EB_,
4218                               meEleISO_chIso_MTD_sim_4_Bkg_EB_,
4219                               meEleISO_chIso_MTD_sim_5_Bkg_EB_,
4220                               meEleISO_chIso_MTD_sim_6_Bkg_EB_,
4221                               meEleISO_chIso_MTD_sim_7_Bkg_EB_};
4222     rel_ch_iso_sim_EB_list_Bkg = {meEleISO_rel_chIso_MTD_sim_1_Bkg_EB_,
4223                                   meEleISO_rel_chIso_MTD_sim_2_Bkg_EB_,
4224                                   meEleISO_rel_chIso_MTD_sim_3_Bkg_EB_,
4225                                   meEleISO_rel_chIso_MTD_sim_4_Bkg_EB_,
4226                                   meEleISO_rel_chIso_MTD_sim_5_Bkg_EB_,
4227                                   meEleISO_rel_chIso_MTD_sim_6_Bkg_EB_,
4228                                   meEleISO_rel_chIso_MTD_sim_7_Bkg_EB_};
4229     Ntracks_sim_EB_list_Significance_Bkg = {meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EB_,
4230                                             meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EB_,
4231                                             meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EB_};
4232     ch_iso_sim_EB_list_Significance_Bkg = {meEleISO_chIso_MTD_sim_4sigma_Bkg_EB_,
4233                                            meEleISO_chIso_MTD_sim_3sigma_Bkg_EB_,
4234                                            meEleISO_chIso_MTD_sim_2sigma_Bkg_EB_};
4235     rel_ch_iso_sim_EB_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EB_,
4236                                                meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EB_,
4237                                                meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EB_};
4238 
4239     Ntracks_sim_EE_list_Bkg = {meEleISO_Ntracks_MTD_sim_1_Bkg_EE_,
4240                                meEleISO_Ntracks_MTD_sim_2_Bkg_EE_,
4241                                meEleISO_Ntracks_MTD_sim_3_Bkg_EE_,
4242                                meEleISO_Ntracks_MTD_sim_4_Bkg_EE_,
4243                                meEleISO_Ntracks_MTD_sim_5_Bkg_EE_,
4244                                meEleISO_Ntracks_MTD_sim_6_Bkg_EE_,
4245                                meEleISO_Ntracks_MTD_sim_7_Bkg_EE_};
4246     ch_iso_sim_EE_list_Bkg = {meEleISO_chIso_MTD_sim_1_Bkg_EE_,
4247                               meEleISO_chIso_MTD_sim_2_Bkg_EE_,
4248                               meEleISO_chIso_MTD_sim_3_Bkg_EE_,
4249                               meEleISO_chIso_MTD_sim_4_Bkg_EE_,
4250                               meEleISO_chIso_MTD_sim_5_Bkg_EE_,
4251                               meEleISO_chIso_MTD_sim_6_Bkg_EE_,
4252                               meEleISO_chIso_MTD_sim_7_Bkg_EE_};
4253     rel_ch_iso_sim_EE_list_Bkg = {meEleISO_rel_chIso_MTD_sim_1_Bkg_EE_,
4254                                   meEleISO_rel_chIso_MTD_sim_2_Bkg_EE_,
4255                                   meEleISO_rel_chIso_MTD_sim_3_Bkg_EE_,
4256                                   meEleISO_rel_chIso_MTD_sim_4_Bkg_EE_,
4257                                   meEleISO_rel_chIso_MTD_sim_5_Bkg_EE_,
4258                                   meEleISO_rel_chIso_MTD_sim_6_Bkg_EE_,
4259                                   meEleISO_rel_chIso_MTD_sim_7_Bkg_EE_};
4260 
4261     Ntracks_sim_EE_list_Significance_Bkg = {meEleISO_Ntracks_MTD_sim_4sigma_Bkg_EE_,
4262                                             meEleISO_Ntracks_MTD_sim_3sigma_Bkg_EE_,
4263                                             meEleISO_Ntracks_MTD_sim_2sigma_Bkg_EE_};
4264     ch_iso_sim_EE_list_Significance_Bkg = {meEleISO_chIso_MTD_sim_4sigma_Bkg_EE_,
4265                                            meEleISO_chIso_MTD_sim_3sigma_Bkg_EE_,
4266                                            meEleISO_chIso_MTD_sim_2sigma_Bkg_EE_};
4267     rel_ch_iso_sim_EE_list_Significance_Bkg = {meEleISO_rel_chIso_MTD_sim_4sigma_Bkg_EE_,
4268                                                meEleISO_rel_chIso_MTD_sim_3sigma_Bkg_EE_,
4269                                                meEleISO_rel_chIso_MTD_sim_2sigma_Bkg_EE_};
4270 
4271     Ele_pT_sim_MTD_EB_list_Bkg = {meEle_pt_sim_MTD_1_Bkg_EB_,
4272                                   meEle_pt_sim_MTD_2_Bkg_EB_,
4273                                   meEle_pt_sim_MTD_3_Bkg_EB_,
4274                                   meEle_pt_sim_MTD_4_Bkg_EB_,
4275                                   meEle_pt_sim_MTD_5_Bkg_EB_,
4276                                   meEle_pt_sim_MTD_6_Bkg_EB_,
4277                                   meEle_pt_sim_MTD_7_Bkg_EB_};
4278 
4279     Ele_pT_sim_MTD_EB_list_Significance_Bkg = {
4280         meEle_pt_sim_MTD_4sigma_Bkg_EB_, meEle_pt_sim_MTD_3sigma_Bkg_EB_, meEle_pt_sim_MTD_2sigma_Bkg_EB_};
4281 
4282     Ele_pT_sim_MTD_EE_list_Bkg = {meEle_pt_sim_MTD_1_Bkg_EE_,
4283                                   meEle_pt_sim_MTD_2_Bkg_EE_,
4284                                   meEle_pt_sim_MTD_3_Bkg_EE_,
4285                                   meEle_pt_sim_MTD_4_Bkg_EE_,
4286                                   meEle_pt_sim_MTD_5_Bkg_EE_,
4287                                   meEle_pt_sim_MTD_6_Bkg_EE_,
4288                                   meEle_pt_sim_MTD_7_Bkg_EE_};
4289 
4290     Ele_pT_sim_MTD_EE_list_Significance_Bkg = {
4291         meEle_pt_sim_MTD_4sigma_Bkg_EE_, meEle_pt_sim_MTD_3sigma_Bkg_EE_, meEle_pt_sim_MTD_2sigma_Bkg_EE_};
4292   }
4293 }
4294 
4295 // ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
4296 
4297 void MtdEleIsoValidation::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
4298   edm::ParameterSetDescription desc;
4299 
4300   desc.add<std::string>("folder", "MTD/ElectronIso");
4301   desc.add<edm::InputTag>("inputTagG", edm::InputTag("generalTracks"));
4302   desc.add<edm::InputTag>("inputTag_vtx", edm::InputTag("offlinePrimaryVertices4D"));
4303   desc.add<edm::InputTag>("inputEle_EB", edm::InputTag("gedGsfElectrons"));
4304   desc.add<edm::InputTag>("inputEle_EE", edm::InputTag("ecalDrivenGsfElectronsHGC"));
4305   desc.add<edm::InputTag>("TPtoRecoTrackAssoc", edm::InputTag("trackingParticleRecoTrackAsssociation"));
4306   desc.add<edm::InputTag>("t0PID", edm::InputTag("tofPID:t0"));
4307   desc.add<edm::InputTag>("sigmat0PID", edm::InputTag("tofPID:sigmat0"));
4308   desc.add<edm::InputTag>("trackMVAQual", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
4309   desc.add<double>("trackMinimumPt", 1.0);  // [GeV]
4310   desc.add<double>("trackMinimumEta", 1.5);
4311   desc.add<double>("trackMaximumEta", 3.2);
4312   desc.add<double>("rel_iso_cut", 0.08);
4313   desc.add<bool>("optionTrackMatchToPV", false);
4314   desc.add<bool>("option_dtToTrack", true);  // default is dt with track, if false will do dt to vertex
4315   desc.add<bool>("option_plots", false);
4316   desc.add<double>("min_dR_cut", 0.01);
4317   desc.add<double>("max_dR_cut", 0.3);
4318   desc.add<double>("min_pt_cut_EB", 0.7);
4319   desc.add<double>("min_pt_cut_EE", 0.4);
4320   desc.add<double>("max_dz_cut_EB", 0.5);  // PARAM
4321   desc.add<double>("max_dz_cut_EE", 0.5);  // PARAM
4322   desc.add<double>("max_dz_vtx_cut", 0.5);
4323   desc.add<double>("max_dxy_vtx_cut", 0.2);
4324   desc.add<double>("min_strip_cut", 0.01);
4325   desc.add<double>("min_track_mtd_mva_cut", 0.5);
4326 
4327   descriptions.add("mtdEleIsoValid", desc);
4328 }
4329 
4330 DEFINE_FWK_MODULE(MtdEleIsoValidation);
4331 
4332 //*/
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.2.1 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!