Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:32:21

0001 #include "Validation/EcalClusters/interface/EgammaSuperClusters.h"
0002 
0003 #include "FWCore/Framework/interface/Event.h"
0004 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0005 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0006 
0007 #include "DataFormats/Common/interface/Handle.h"
0008 #include "DataFormats/EgammaReco/interface/BasicCluster.h"
0009 #include "DataFormats/EgammaReco/interface/SuperCluster.h"
0010 #include "DataFormats/Math/interface/deltaPhi.h"
0011 
0012 #include "DQMServices/Core/interface/DQMStore.h"
0013 
0014 #include "FWCore/Framework/interface/MakerMacros.h"
0015 
0016 EgammaSuperClusters::EgammaSuperClusters(const edm::ParameterSet &ps)
0017     : MCTruthCollectionToken_(consumes<edm::HepMCProduct>(ps.getParameter<edm::InputTag>("MCTruthCollection"))),
0018       barrelRawSuperClusterCollectionToken_(
0019           consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelRawSuperClusterCollection"))),
0020       barrelCorSuperClusterCollectionToken_(
0021           consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelCorSuperClusterCollection"))),
0022       endcapRawSuperClusterCollectionToken_(
0023           consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapRawSuperClusterCollection"))),
0024       endcapPreSuperClusterCollectionToken_(
0025           consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapPreSuperClusterCollection"))),
0026       endcapCorSuperClusterCollectionToken_(
0027           consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapCorSuperClusterCollection"))),
0028       barrelRecHitCollectionToken_(
0029           consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelRecHitCollection"))),
0030       endcapRecHitCollectionToken_(
0031           consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapRecHitCollection"))),
0032       ecalClusterToolsESGetTokens_{consumesCollector()},
0033       hsSize_(ps, "Size"),
0034       hsNumBC_(ps, "NumBC"),
0035       hsET_(ps, "ET"),
0036       hsEta_(ps, "Eta"),
0037       hsPhi_(ps, "Phi"),
0038       hsS1toS9_(ps, "S1toS9"),
0039       hsS25toE_(ps, "S25toE"),
0040       hsEoverTruth_(ps, "EoverTruth"),
0041       hsdeltaR_(ps, "deltaR"),
0042       hsphiWidth_(ps, "phiWidth"),
0043       hsetaWidth_(ps, "etaWidth"),
0044       hspreshowerE_(ps, "preshowerE"),
0045       hsR_(ps, "R"),
0046       hist_EB_RawSC_Size_(nullptr),
0047       hist_EE_RawSC_Size_(nullptr),
0048       hist_EB_CorSC_Size_(nullptr),
0049       hist_EE_CorSC_Size_(nullptr),
0050       hist_EE_PreSC_Size_(nullptr),
0051       hist_EB_RawSC_NumBC_(nullptr),
0052       hist_EE_RawSC_NumBC_(nullptr),
0053       hist_EB_CorSC_NumBC_(nullptr),
0054       hist_EE_CorSC_NumBC_(nullptr),
0055       hist_EE_PreSC_NumBC_(nullptr),
0056       hist_EB_RawSC_ET_(nullptr),
0057       hist_EE_RawSC_ET_(nullptr),
0058       hist_EB_CorSC_ET_(nullptr),
0059       hist_EE_CorSC_ET_(nullptr),
0060       hist_EE_PreSC_ET_(nullptr),
0061       hist_EB_RawSC_Eta_(nullptr),
0062       hist_EE_RawSC_Eta_(nullptr),
0063       hist_EB_CorSC_Eta_(nullptr),
0064       hist_EE_CorSC_Eta_(nullptr),
0065       hist_EE_PreSC_Eta_(nullptr),
0066       hist_EB_RawSC_Phi_(nullptr),
0067       hist_EE_RawSC_Phi_(nullptr),
0068       hist_EB_CorSC_Phi_(nullptr),
0069       hist_EE_CorSC_Phi_(nullptr),
0070       hist_EE_PreSC_Phi_(nullptr),
0071       hist_EB_RawSC_S1toS9_(nullptr),
0072       hist_EE_RawSC_S1toS9_(nullptr),
0073       hist_EB_CorSC_S1toS9_(nullptr),
0074       hist_EE_CorSC_S1toS9_(nullptr),
0075       hist_EE_PreSC_S1toS9_(nullptr),
0076       hist_EB_RawSC_S25toE_(nullptr),
0077       hist_EE_RawSC_S25toE_(nullptr),
0078       hist_EB_CorSC_S25toE_(nullptr),
0079       hist_EE_CorSC_S25toE_(nullptr),
0080       hist_EE_PreSC_S25toE_(nullptr),
0081       hist_EB_RawSC_EoverTruth_(nullptr),
0082       hist_EE_RawSC_EoverTruth_(nullptr),
0083       hist_EB_CorSC_EoverTruth_(nullptr),
0084       hist_EE_CorSC_EoverTruth_(nullptr),
0085       hist_EE_PreSC_EoverTruth_(nullptr),
0086       hist_EB_RawSC_deltaR_(nullptr),
0087       hist_EE_RawSC_deltaR_(nullptr),
0088       hist_EB_CorSC_deltaR_(nullptr),
0089       hist_EE_CorSC_deltaR_(nullptr),
0090       hist_EE_PreSC_deltaR_(nullptr),
0091       hist_EE_PreSC_preshowerE_(nullptr),
0092       hist_EE_CorSC_preshowerE_(nullptr),
0093       hist_EE_CorSC_phiWidth_(nullptr),
0094       hist_EB_CorSC_phiWidth_(nullptr),
0095       hist_EE_CorSC_etaWidth_(nullptr),
0096       hist_EB_CorSC_etaWidth_(nullptr),
0097       hist_EB_CorSC_ET_vs_Eta_(nullptr),
0098       hist_EB_CorSC_ET_vs_Phi_(nullptr),
0099       hist_EE_CorSC_ET_vs_Eta_(nullptr),
0100       hist_EE_CorSC_ET_vs_Phi_(nullptr),
0101       hist_EE_CorSC_ET_vs_R_(nullptr) {}
0102 
0103 EgammaSuperClusters::~EgammaSuperClusters() {}
0104 
0105 void EgammaSuperClusters::bookHistograms(DQMStore::IBooker &_ibooker, edm::Run const &, edm::EventSetup const &) {
0106   _ibooker.setCurrentFolder("EcalClusterV/EcalSuperClusters/");
0107 
0108   // Number of SuperClusters
0109   //
0110   hist_EB_RawSC_Size_ =
0111       _ibooker.book1D("hist_EB_RawSC_Size_", "# Raw SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
0112   hist_EE_RawSC_Size_ =
0113       _ibooker.book1D("hist_EE_RawSC_Size_", "# Raw SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
0114   hist_EB_CorSC_Size_ = _ibooker.book1D(
0115       "hist_EB_CorSC_Size_", "# Corrected SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
0116   hist_EE_CorSC_Size_ = _ibooker.book1D(
0117       "hist_EE_CorSC_Size_", "# Corrected SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
0118   hist_EE_PreSC_Size_ = _ibooker.book1D(
0119       "hist_EE_PreSC_Size_", "# SuperClusters with Preshower in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
0120 
0121   // Number of BasicClusters in SuperCluster
0122   //
0123   hist_EB_RawSC_NumBC_ = _ibooker.book1D("hist_EB_RawSC_NumBC_",
0124                                          "# of Basic Clusters in Raw Super Clusters in Barrel",
0125                                          hsNumBC_.bins,
0126                                          hsNumBC_.min,
0127                                          hsNumBC_.max);
0128   hist_EE_RawSC_NumBC_ = _ibooker.book1D("hist_EE_RawSC_NumBC_",
0129                                          "# of Basic Clusters in Raw Super Clusters in Endcap",
0130                                          hsNumBC_.bins,
0131                                          hsNumBC_.min,
0132                                          hsNumBC_.max);
0133   hist_EB_CorSC_NumBC_ = _ibooker.book1D("hist_EB_CorSC_NumBC_",
0134                                          "# of Basic Clusters in Corrected SuperClusters in Barrel",
0135                                          hsNumBC_.bins,
0136                                          hsNumBC_.min,
0137                                          hsNumBC_.max);
0138   hist_EE_CorSC_NumBC_ = _ibooker.book1D("hist_EE_CorSC_NumBC_",
0139                                          "# of Basic Clusters in Corrected SuperClusters in Endcap",
0140                                          hsNumBC_.bins,
0141                                          hsNumBC_.min,
0142                                          hsNumBC_.max);
0143   hist_EE_PreSC_NumBC_ = _ibooker.book1D("hist_EE_PreSC_NumBC_",
0144                                          "# of Basic Clusters in SuperClusters with Preshower in Endcap",
0145                                          hsNumBC_.bins,
0146                                          hsNumBC_.min,
0147                                          hsNumBC_.max);
0148 
0149   // ET distribution of SuperClusters
0150   //
0151   hist_EB_RawSC_ET_ =
0152       _ibooker.book1D("hist_EB_RawSC_ET_", "ET of Raw SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
0153   hist_EE_RawSC_ET_ =
0154       _ibooker.book1D("hist_EE_RawSC_ET_", "ET of Raw SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
0155   hist_EB_CorSC_ET_ =
0156       _ibooker.book1D("hist_EB_CorSC_ET_", "ET of Corrected SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
0157   hist_EE_CorSC_ET_ =
0158       _ibooker.book1D("hist_EE_CorSC_ET_", "ET of Corrected SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
0159   hist_EE_PreSC_ET_ = _ibooker.book1D(
0160       "hist_EE_PreSC_ET_", "ET of SuperClusters with Preshower in Endcap", hsET_.bins, hsET_.min, hsET_.max);
0161 
0162   // Eta distribution of SuperClusters
0163   //
0164   hist_EB_RawSC_Eta_ =
0165       _ibooker.book1D("hist_EB_RawSC_Eta_", "Eta of Raw SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
0166   hist_EE_RawSC_Eta_ =
0167       _ibooker.book1D("hist_EE_RawSC_Eta_", "Eta of Raw SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
0168   hist_EB_CorSC_Eta_ = _ibooker.book1D(
0169       "hist_EB_CorSC_Eta_", "Eta of Corrected SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
0170   hist_EE_CorSC_Eta_ = _ibooker.book1D(
0171       "hist_EE_CorSC_Eta_", "Eta of Corrected SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
0172   hist_EE_PreSC_Eta_ = _ibooker.book1D(
0173       "hist_EE_PreSC_Eta_", "Eta of SuperClusters with Preshower in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
0174 
0175   // Phi distribution of SuperClusters
0176   //
0177   hist_EB_RawSC_Phi_ =
0178       _ibooker.book1D("hist_EB_RawSC_Phi_", "Phi of Raw SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
0179   hist_EE_RawSC_Phi_ =
0180       _ibooker.book1D("hist_EE_RawSC_Phi_", "Phi of Raw SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
0181   hist_EB_CorSC_Phi_ = _ibooker.book1D(
0182       "hist_EB_CorSC_Phi_", "Phi of Corrected SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
0183   hist_EE_CorSC_Phi_ = _ibooker.book1D(
0184       "hist_EE_CorSC_Phi_", "Phi of Corrected SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
0185   hist_EE_PreSC_Phi_ = _ibooker.book1D(
0186       "hist_EE_PreSC_Phi_", "Phi of SuperClusters with Preshower in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
0187 
0188   // S1/S9 distribution of SuperClusters
0189   //
0190   hist_EB_RawSC_S1toS9_ = _ibooker.book1D(
0191       "hist_EB_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Barrel", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
0192   hist_EE_RawSC_S1toS9_ = _ibooker.book1D(
0193       "hist_EE_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Endcap", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
0194   hist_EB_CorSC_S1toS9_ = _ibooker.book1D("hist_EB_CorSC_S1toS9_",
0195                                           "S1/S9 of Corrected SuperClusters in Barrel",
0196                                           hsS1toS9_.bins,
0197                                           hsS1toS9_.min,
0198                                           hsS1toS9_.max);
0199   hist_EE_CorSC_S1toS9_ = _ibooker.book1D("hist_EE_CorSC_S1toS9_",
0200                                           "S1/S9 of Corrected SuperClusters in Endcap",
0201                                           hsS1toS9_.bins,
0202                                           hsS1toS9_.min,
0203                                           hsS1toS9_.max);
0204   hist_EE_PreSC_S1toS9_ = _ibooker.book1D("hist_EE_PreSC_S1toS9_",
0205                                           "S1/S9 of SuperClusters with Preshower in Endcap",
0206                                           hsS1toS9_.bins,
0207                                           hsS1toS9_.min,
0208                                           hsS1toS9_.max);
0209 
0210   // S25/E distribution of SuperClusters
0211   //
0212   hist_EB_RawSC_S25toE_ = _ibooker.book1D(
0213       "hist_EB_RawSC_S25toE_", "S25/E of Raw Super Clusters in Barrel", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
0214   hist_EE_RawSC_S25toE_ = _ibooker.book1D(
0215       "hist_EE_RawSC_S25toE_", "S25/E of Raw Super Clusters in Endcap", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
0216   hist_EB_CorSC_S25toE_ = _ibooker.book1D("hist_EB_CorSC_S25toE_",
0217                                           "S25/E of Corrected SuperClusters in Barrel",
0218                                           hsS25toE_.bins,
0219                                           hsS25toE_.min,
0220                                           hsS25toE_.max);
0221   hist_EE_CorSC_S25toE_ = _ibooker.book1D("hist_EE_CorSC_S25toE_",
0222                                           "S25/E of Corrected SuperClusters in Endcap",
0223                                           hsS25toE_.bins,
0224                                           hsS25toE_.min,
0225                                           hsS25toE_.max);
0226   hist_EE_PreSC_S25toE_ = _ibooker.book1D("hist_EE_PreSC_S25toE_",
0227                                           "S25/E of SuperClusters with Preshower in Endcap",
0228                                           hsS25toE_.bins,
0229                                           hsS25toE_.min,
0230                                           hsS25toE_.max);
0231 
0232   // E/E(true) distribution of SuperClusters
0233   //
0234   hist_EB_RawSC_EoverTruth_ = _ibooker.book1D("hist_EB_RawSC_EoverTruth_",
0235                                               "E/True E of Raw SuperClusters in Barrel",
0236                                               hsEoverTruth_.bins,
0237                                               hsEoverTruth_.min,
0238                                               hsEoverTruth_.max);
0239   hist_EE_RawSC_EoverTruth_ = _ibooker.book1D("hist_EE_RawSC_EoverTruth_",
0240                                               "E/True E of Raw SuperClusters in Endcap",
0241                                               hsEoverTruth_.bins,
0242                                               hsEoverTruth_.min,
0243                                               hsEoverTruth_.max);
0244   hist_EB_CorSC_EoverTruth_ = _ibooker.book1D("hist_EB_CorSC_EoverTruth_",
0245                                               "E/True E of Corrected SuperClusters in Barrel",
0246                                               hsEoverTruth_.bins,
0247                                               hsEoverTruth_.min,
0248                                               hsEoverTruth_.max);
0249   hist_EE_CorSC_EoverTruth_ = _ibooker.book1D("hist_EE_CorSC_EoverTruth_",
0250                                               "E/True E of Corrected SuperClusters in Endcap",
0251                                               hsEoverTruth_.bins,
0252                                               hsEoverTruth_.min,
0253                                               hsEoverTruth_.max);
0254   hist_EE_PreSC_EoverTruth_ = _ibooker.book1D("hist_EE_PreSC_EoverTruth_",
0255                                               "E/True E of SuperClusters with Preshower in Endcap",
0256                                               hsEoverTruth_.bins,
0257                                               hsEoverTruth_.min,
0258                                               hsEoverTruth_.max);
0259 
0260   // dR distribution of SuperClusters from truth
0261   //
0262   hist_EB_RawSC_deltaR_ = _ibooker.book1D("hist_EB_RawSC_deltaR_",
0263                                           "dR to MC truth of Raw Super Clusters in Barrel",
0264                                           hsdeltaR_.bins,
0265                                           hsdeltaR_.min,
0266                                           hsdeltaR_.max);
0267   hist_EE_RawSC_deltaR_ = _ibooker.book1D("hist_EE_RawSC_deltaR_",
0268                                           "dR to MC truth of Raw Super Clusters in Endcap",
0269                                           hsdeltaR_.bins,
0270                                           hsdeltaR_.min,
0271                                           hsdeltaR_.max);
0272   hist_EB_CorSC_deltaR_ = _ibooker.book1D("hist_EB_CorSC_deltaR_",
0273                                           "dR to MC truth of Corrected SuperClusters in Barrel",
0274                                           hsdeltaR_.bins,
0275                                           hsdeltaR_.min,
0276                                           hsdeltaR_.max);
0277   hist_EE_CorSC_deltaR_ = _ibooker.book1D("hist_EE_CorSC_deltaR_",
0278                                           "dR to MC truth of Corrected SuperClusters in Endcap",
0279                                           hsdeltaR_.bins,
0280                                           hsdeltaR_.min,
0281                                           hsdeltaR_.max);
0282   hist_EE_PreSC_deltaR_ = _ibooker.book1D("hist_EE_PreSC_deltaR_",
0283                                           "dR to MC truth of SuperClusters with Preshower in Endcap",
0284                                           hsdeltaR_.bins,
0285                                           hsdeltaR_.min,
0286                                           hsdeltaR_.max);
0287 
0288   // phi width stored in corrected SuperClusters
0289   hist_EB_CorSC_phiWidth_ = _ibooker.book1D("hist_EB_CorSC_phiWidth_",
0290                                             "phiWidth of Corrected Super Clusters in Barrel",
0291                                             hsphiWidth_.bins,
0292                                             hsphiWidth_.min,
0293                                             hsphiWidth_.max);
0294   hist_EE_CorSC_phiWidth_ = _ibooker.book1D("hist_EE_CorSC_phiWidth_",
0295                                             "phiWidth of Corrected Super Clusters in Endcap",
0296                                             hsphiWidth_.bins,
0297                                             hsphiWidth_.min,
0298                                             hsphiWidth_.max);
0299 
0300   // eta width stored in corrected SuperClusters
0301   hist_EB_CorSC_etaWidth_ = _ibooker.book1D("hist_EB_CorSC_etaWidth_",
0302                                             "etaWidth of Corrected Super Clusters in Barrel",
0303                                             hsetaWidth_.bins,
0304                                             hsetaWidth_.min,
0305                                             hsetaWidth_.max);
0306   hist_EE_CorSC_etaWidth_ = _ibooker.book1D("hist_EE_CorSC_etaWidth_",
0307                                             "etaWidth of Corrected Super Clusters in Endcap",
0308                                             hsetaWidth_.bins,
0309                                             hsetaWidth_.min,
0310                                             hsetaWidth_.max);
0311 
0312   // preshower energy
0313   hist_EE_PreSC_preshowerE_ = _ibooker.book1D("hist_EE_PreSC_preshowerE_",
0314                                               "preshower energy in Super Clusters with Preshower in Endcap",
0315                                               hspreshowerE_.bins,
0316                                               hspreshowerE_.min,
0317                                               hspreshowerE_.max);
0318   hist_EE_CorSC_preshowerE_ = _ibooker.book1D("hist_EE_CorSC_preshowerE_",
0319                                               "preshower energy in Corrected Super Clusters with Preshower in Endcap",
0320                                               hspreshowerE_.bins,
0321                                               hspreshowerE_.min,
0322                                               hspreshowerE_.max);
0323 
0324   //
0325   hist_EB_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Eta_",
0326                                              "Corr Super Cluster ET versus Eta in Barrel",
0327                                              hsET_.bins,
0328                                              hsET_.min,
0329                                              hsET_.max,
0330                                              hsEta_.bins,
0331                                              hsEta_.min,
0332                                              hsEta_.max);
0333 
0334   hist_EB_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Phi_",
0335                                              "Corr Super Cluster ET versus Phi in Barrel",
0336                                              hsET_.bins,
0337                                              hsET_.min,
0338                                              hsET_.max,
0339                                              hsPhi_.bins,
0340                                              hsPhi_.min,
0341                                              hsPhi_.max);
0342 
0343   hist_EE_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Eta_",
0344                                              "Corr Super Cluster ET versus Eta in Endcap",
0345                                              hsET_.bins,
0346                                              hsET_.min,
0347                                              hsET_.max,
0348                                              hsEta_.bins,
0349                                              hsEta_.min,
0350                                              hsEta_.max);
0351 
0352   hist_EE_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Phi_",
0353                                              "Corr Super Cluster ET versus Phi in Endcap",
0354                                              hsET_.bins,
0355                                              hsET_.min,
0356                                              hsET_.max,
0357                                              hsPhi_.bins,
0358                                              hsPhi_.min,
0359                                              hsPhi_.max);
0360 
0361   hist_EE_CorSC_ET_vs_R_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_R_",
0362                                            "Corr Super Cluster ET versus Radius in Endcap",
0363                                            hsET_.bins,
0364                                            hsET_.min,
0365                                            hsET_.max,
0366                                            hsR_.bins,
0367                                            hsR_.min,
0368                                            hsR_.max);
0369 }
0370 
0371 void EgammaSuperClusters::analyze(const edm::Event &evt, const edm::EventSetup &es) {
0372   bool skipMC = false;
0373   bool skipBarrel = false;
0374   bool skipEndcap = false;
0375 
0376   //
0377   // Get MCTRUTH
0378   //
0379   edm::Handle<edm::HepMCProduct> pMCTruth;
0380   evt.getByToken(MCTruthCollectionToken_, pMCTruth);
0381   if (!pMCTruth.isValid()) {
0382     edm::LogError("EgammaSuperClusters") << "Error! can't get MC collection ";
0383     skipMC = true;
0384   }
0385   const HepMC::GenEvent *genEvent = pMCTruth->GetEvent();
0386 
0387   if (skipMC)
0388     return;
0389 
0390   //
0391   // Get the BARREL products
0392   //
0393   edm::Handle<reco::SuperClusterCollection> pBarrelRawSuperClusters;
0394   evt.getByToken(barrelRawSuperClusterCollectionToken_, pBarrelRawSuperClusters);
0395   if (!pBarrelRawSuperClusters.isValid()) {
0396     edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw SC";
0397     skipBarrel = true;
0398   }
0399 
0400   edm::Handle<reco::SuperClusterCollection> pBarrelCorSuperClusters;
0401   evt.getByToken(barrelCorSuperClusterCollectionToken_, pBarrelCorSuperClusters);
0402   if (!pBarrelCorSuperClusters.isValid()) {
0403     edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor SC";
0404     skipBarrel = true;
0405   }
0406 
0407   edm::Handle<EBRecHitCollection> pBarrelRecHitCollection;
0408   evt.getByToken(barrelRecHitCollectionToken_, pBarrelRecHitCollection);
0409   if (!pBarrelRecHitCollection.isValid()) {
0410     skipBarrel = true;
0411   }
0412   edm::Handle<EERecHitCollection> pEndcapRecHitCollection;
0413   evt.getByToken(endcapRecHitCollectionToken_, pEndcapRecHitCollection);
0414   if (!pEndcapRecHitCollection.isValid()) {
0415     skipEndcap = true;
0416   }
0417 
0418   if (skipBarrel || skipEndcap)
0419     return;
0420 
0421   EcalClusterLazyTools lazyTool(
0422       evt, ecalClusterToolsESGetTokens_.get(es), barrelRecHitCollectionToken_, endcapRecHitCollectionToken_);
0423 
0424   // Get the BARREL collections
0425   const reco::SuperClusterCollection *barrelRawSuperClusters = pBarrelRawSuperClusters.product();
0426   const reco::SuperClusterCollection *barrelCorSuperClusters = pBarrelCorSuperClusters.product();
0427 
0428   // Number of entries in collections
0429   hist_EB_RawSC_Size_->Fill(barrelRawSuperClusters->size());
0430   hist_EB_CorSC_Size_->Fill(barrelCorSuperClusters->size());
0431 
0432   // Do RAW BARREL SuperClusters
0433   for (reco::SuperClusterCollection::const_iterator aClus = barrelRawSuperClusters->begin();
0434        aClus != barrelRawSuperClusters->end();
0435        aClus++) {
0436     // kinematics
0437     hist_EB_RawSC_NumBC_->Fill(aClus->clustersSize());
0438     hist_EB_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
0439     hist_EB_RawSC_Eta_->Fill(aClus->position().eta());
0440     hist_EB_RawSC_Phi_->Fill(aClus->position().phi());
0441 
0442     // cluster shape
0443     const reco::CaloClusterPtr seed = aClus->seed();
0444     hist_EB_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
0445     hist_EB_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
0446 
0447     // truth
0448     double dRClosest = 999.9;
0449     double energyClosest = 0;
0450     closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
0451 
0452     if (dRClosest < 0.1) {
0453       hist_EB_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
0454       hist_EB_RawSC_deltaR_->Fill(dRClosest);
0455     }
0456   }
0457 
0458   // Do CORRECTED BARREL SuperClusters
0459   for (reco::SuperClusterCollection::const_iterator aClus = barrelCorSuperClusters->begin();
0460        aClus != barrelCorSuperClusters->end();
0461        aClus++) {
0462     // kinematics
0463     hist_EB_CorSC_NumBC_->Fill(aClus->clustersSize());
0464     hist_EB_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
0465     hist_EB_CorSC_Eta_->Fill(aClus->position().eta());
0466     hist_EB_CorSC_Phi_->Fill(aClus->position().phi());
0467 
0468     hist_EB_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
0469     hist_EB_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());
0470 
0471     // cluster shape
0472     const reco::CaloClusterPtr seed = aClus->seed();
0473     hist_EB_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
0474     hist_EB_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
0475 
0476     // correction variables
0477     hist_EB_CorSC_phiWidth_->Fill(aClus->phiWidth());
0478     hist_EB_CorSC_etaWidth_->Fill(aClus->etaWidth());
0479 
0480     // truth
0481     double dRClosest = 999.9;
0482     double energyClosest = 0;
0483     closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
0484 
0485     if (dRClosest < 0.1) {
0486       hist_EB_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
0487       hist_EB_CorSC_deltaR_->Fill(dRClosest);
0488     }
0489   }
0490 
0491   //
0492   // Get the ENDCAP products
0493   //
0494   edm::Handle<reco::SuperClusterCollection> pEndcapRawSuperClusters;
0495   evt.getByToken(endcapRawSuperClusterCollectionToken_, pEndcapRawSuperClusters);
0496   if (!pEndcapRawSuperClusters.isValid()) {
0497     edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw EE SC";
0498   }
0499 
0500   edm::Handle<reco::SuperClusterCollection> pEndcapPreSuperClusters;
0501   evt.getByToken(endcapPreSuperClusterCollectionToken_, pEndcapPreSuperClusters);
0502   if (!pEndcapPreSuperClusters.isValid()) {
0503     edm::LogError("EgammaSuperClusters") << "Error! can't get collection Pre EE SC";
0504   }
0505 
0506   edm::Handle<reco::SuperClusterCollection> pEndcapCorSuperClusters;
0507   evt.getByToken(endcapCorSuperClusterCollectionToken_, pEndcapCorSuperClusters);
0508   if (!pEndcapCorSuperClusters.isValid()) {
0509     edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor EE SC";
0510   }
0511 
0512   // Get the ENDCAP collections
0513   const reco::SuperClusterCollection *endcapRawSuperClusters = pEndcapRawSuperClusters.product();
0514   const reco::SuperClusterCollection *endcapPreSuperClusters = pEndcapPreSuperClusters.product();
0515   const reco::SuperClusterCollection *endcapCorSuperClusters = pEndcapCorSuperClusters.product();
0516 
0517   // Number of entries in collections
0518   hist_EE_RawSC_Size_->Fill(endcapRawSuperClusters->size());
0519   hist_EE_PreSC_Size_->Fill(endcapPreSuperClusters->size());
0520   hist_EE_CorSC_Size_->Fill(endcapCorSuperClusters->size());
0521 
0522   // Do RAW ENDCAP SuperClusters
0523   for (reco::SuperClusterCollection::const_iterator aClus = endcapRawSuperClusters->begin();
0524        aClus != endcapRawSuperClusters->end();
0525        aClus++) {
0526     hist_EE_RawSC_NumBC_->Fill(aClus->clustersSize());
0527     hist_EE_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
0528     hist_EE_RawSC_Eta_->Fill(aClus->position().eta());
0529     hist_EE_RawSC_Phi_->Fill(aClus->position().phi());
0530 
0531     const reco::CaloClusterPtr seed = aClus->seed();
0532     hist_EE_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
0533     hist_EE_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
0534 
0535     // truth
0536     double dRClosest = 999.9;
0537     double energyClosest = 0;
0538     closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
0539 
0540     if (dRClosest < 0.1) {
0541       hist_EE_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
0542       hist_EE_RawSC_deltaR_->Fill(dRClosest);
0543     }
0544   }
0545 
0546   // Do ENDCAP SuperClusters with PRESHOWER
0547   for (reco::SuperClusterCollection::const_iterator aClus = endcapPreSuperClusters->begin();
0548        aClus != endcapPreSuperClusters->end();
0549        aClus++) {
0550     hist_EE_PreSC_NumBC_->Fill(aClus->clustersSize());
0551     hist_EE_PreSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
0552     hist_EE_PreSC_Eta_->Fill(aClus->position().eta());
0553     hist_EE_PreSC_Phi_->Fill(aClus->position().phi());
0554     hist_EE_PreSC_preshowerE_->Fill(aClus->preshowerEnergy());
0555 
0556     const reco::CaloClusterPtr seed = aClus->seed();
0557     hist_EE_PreSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
0558     hist_EE_PreSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
0559 
0560     // truth
0561     double dRClosest = 999.9;
0562     double energyClosest = 0;
0563     closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
0564 
0565     if (dRClosest < 0.1) {
0566       hist_EE_PreSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
0567       hist_EE_PreSC_deltaR_->Fill(dRClosest);
0568     }
0569   }
0570 
0571   // Do CORRECTED ENDCAP SuperClusters
0572   for (reco::SuperClusterCollection::const_iterator aClus = endcapCorSuperClusters->begin();
0573        aClus != endcapCorSuperClusters->end();
0574        aClus++) {
0575     hist_EE_CorSC_NumBC_->Fill(aClus->clustersSize());
0576     hist_EE_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
0577     hist_EE_CorSC_Eta_->Fill(aClus->position().eta());
0578     hist_EE_CorSC_Phi_->Fill(aClus->position().phi());
0579     hist_EE_CorSC_preshowerE_->Fill(aClus->preshowerEnergy());
0580 
0581     hist_EE_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
0582     hist_EE_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());
0583     hist_EE_CorSC_ET_vs_R_->Fill(aClus->energy() / std::cosh(aClus->position().eta()),
0584                                  std::sqrt(std::pow(aClus->x(), 2) + std::pow(aClus->y(), 2)));
0585 
0586     // correction variables
0587     hist_EE_CorSC_phiWidth_->Fill(aClus->phiWidth());
0588     hist_EE_CorSC_etaWidth_->Fill(aClus->etaWidth());
0589 
0590     const reco::CaloClusterPtr seed = aClus->seed();
0591     hist_EE_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
0592     hist_EE_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
0593 
0594     // truth
0595     double dRClosest = 999.9;
0596     double energyClosest = 0;
0597     closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
0598 
0599     if (dRClosest < 0.1) {
0600       hist_EE_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
0601       hist_EE_CorSC_deltaR_->Fill(dRClosest);
0602     }
0603   }
0604 }
0605 
0606 //
0607 // Closest MC Particle
0608 //
0609 void EgammaSuperClusters::closestMCParticle(const HepMC::GenEvent *genEvent,
0610                                             const reco::SuperCluster &sc,
0611                                             double &dRClosest,
0612                                             double &energyClosest) const {
0613   // SuperCluster eta, phi
0614   double scEta = sc.eta();
0615   double scPhi = sc.phi();
0616 
0617   // initialize dRClosest to a large number
0618   dRClosest = 999.9;
0619 
0620   // loop over the MC truth particles to find the
0621   // closest to the superCluster in dR space
0622   for (HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
0623        currentParticle != genEvent->particles_end();
0624        currentParticle++) {
0625     if ((*currentParticle)->status() == 1) {
0626       // need GenParticle in ECAL co-ordinates
0627       HepMC::FourVector vtx = (*currentParticle)->production_vertex()->position();
0628       double phiTrue = (*currentParticle)->momentum().phi();
0629       double etaTrue = ecalEta((*currentParticle)->momentum().eta(), vtx.z() / 10., vtx.perp() / 10.);
0630 
0631       double dPhi = reco::deltaPhi(phiTrue, scPhi);
0632       double dEta = scEta - etaTrue;
0633       double deltaR = std::sqrt(dPhi * dPhi + dEta * dEta);
0634 
0635       if (deltaR < dRClosest) {
0636         dRClosest = deltaR;
0637         energyClosest = (*currentParticle)->momentum().e();
0638       }
0639 
0640     }  // end if stable particle
0641 
0642   }  // end loop on get particles
0643 }
0644 
0645 //
0646 // Compute Eta in the ECAL co-ordinate system
0647 //
0648 float EgammaSuperClusters::ecalEta(float EtaParticle, float Zvertex, float plane_Radius) const {
0649   const float R_ECAL = 136.5;
0650   const float Z_Endcap = 328.0;
0651   const float etaBarrelEndcap = 1.479;
0652 
0653   if (EtaParticle != 0.) {
0654     float Theta = 0.0;
0655     float ZEcal = (R_ECAL - plane_Radius) * sinh(EtaParticle) + Zvertex;
0656 
0657     if (ZEcal != 0.0)
0658       Theta = atan(R_ECAL / ZEcal);
0659     if (Theta < 0.0)
0660       Theta = Theta + Geom::pi();
0661 
0662     float ETA = -log(tan(0.5 * Theta));
0663 
0664     if (fabs(ETA) > etaBarrelEndcap) {
0665       float Zend = Z_Endcap;
0666       if (EtaParticle < 0.0)
0667         Zend = -Zend;
0668       float Zlen = Zend - Zvertex;
0669       float RR = Zlen / sinh(EtaParticle);
0670       Theta = atan((RR + plane_Radius) / Zend);
0671       if (Theta < 0.0)
0672         Theta = Theta + Geom::pi();
0673       ETA = -log(tan(0.5 * Theta));
0674     }
0675 
0676     return ETA;
0677   } else {
0678     edm::LogWarning("") << "[EgammaSuperClusters::ecalEta] Warning: Eta equals "
0679                            "to zero, not correcting";
0680     return EtaParticle;
0681   }
0682 }