EgammaSuperClusters.h

Go to the documentation of this file.

#ifndef Validation_EcalClusters_EgammaSuperClusters_h
#define Validation_EcalClusters_EgammaSuperClusters_h

//
//

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

#include "HistSpec.h"

#include <string>

class MonitorElement;

class EgammaSuperClusters : public DQMEDAnalyzer
{
 public:
  explicit EgammaSuperClusters( const edm::ParameterSet& );
  ~EgammaSuperClusters() override;
  
  void analyze( const edm::Event&, const edm::EventSetup& ) override;
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
  
 private:
  // mc truth
  edm::EDGetTokenT<edm::HepMCProduct> MCTruthCollectionToken_;
  
  // barrel clusters
  edm::EDGetTokenT<reco::SuperClusterCollection> barrelRawSuperClusterCollectionToken_;
  edm::EDGetTokenT<reco::SuperClusterCollection> barrelCorSuperClusterCollectionToken_;
  
  // endcap clusters
  edm::EDGetTokenT<reco::SuperClusterCollection> endcapRawSuperClusterCollectionToken_;
  edm::EDGetTokenT<reco::SuperClusterCollection> endcapPreSuperClusterCollectionToken_;
  edm::EDGetTokenT<reco::SuperClusterCollection> endcapCorSuperClusterCollectionToken_;
  
  // collections of hits
  edm::EDGetTokenT<EcalRecHitCollection> barrelRecHitCollectionToken_;
  edm::EDGetTokenT<EcalRecHitCollection> endcapRecHitCollectionToken_;

  HistSpec hsSize_;
  HistSpec hsNumBC_;
  HistSpec hsET_;
  HistSpec hsEta_;
  HistSpec hsPhi_;
  HistSpec hsS1toS9_;
  HistSpec hsS25toE_;
  HistSpec hsEoverTruth_;
  HistSpec hsdeltaR_;
  HistSpec hsphiWidth_;
  HistSpec hsetaWidth_;
  HistSpec hspreshowerE_;
  HistSpec hsR_;

  MonitorElement* hist_EB_RawSC_Size_;
  MonitorElement* hist_EE_RawSC_Size_;
  MonitorElement* hist_EB_CorSC_Size_;
  MonitorElement* hist_EE_CorSC_Size_;
  MonitorElement* hist_EE_PreSC_Size_;
  MonitorElement* hist_EB_RawSC_NumBC_;
  MonitorElement* hist_EE_RawSC_NumBC_;
  MonitorElement* hist_EB_CorSC_NumBC_;
  MonitorElement* hist_EE_CorSC_NumBC_;
  MonitorElement* hist_EE_PreSC_NumBC_;
  MonitorElement* hist_EB_RawSC_ET_;
  MonitorElement* hist_EE_RawSC_ET_;
  MonitorElement* hist_EB_CorSC_ET_;
  MonitorElement* hist_EE_CorSC_ET_;
  MonitorElement* hist_EE_PreSC_ET_;
  MonitorElement* hist_EB_RawSC_Eta_;
  MonitorElement* hist_EE_RawSC_Eta_;
  MonitorElement* hist_EB_CorSC_Eta_;
  MonitorElement* hist_EE_CorSC_Eta_;
  MonitorElement* hist_EE_PreSC_Eta_;
  MonitorElement* hist_EB_RawSC_Phi_;
  MonitorElement* hist_EE_RawSC_Phi_;
  MonitorElement* hist_EB_CorSC_Phi_;
  MonitorElement* hist_EE_CorSC_Phi_;
  MonitorElement* hist_EE_PreSC_Phi_;
  MonitorElement* hist_EB_RawSC_S1toS9_;
  MonitorElement* hist_EE_RawSC_S1toS9_;
  MonitorElement* hist_EB_CorSC_S1toS9_;
  MonitorElement* hist_EE_CorSC_S1toS9_;
  MonitorElement* hist_EE_PreSC_S1toS9_;
  MonitorElement* hist_EB_RawSC_S25toE_;
  MonitorElement* hist_EE_RawSC_S25toE_;
  MonitorElement* hist_EB_CorSC_S25toE_;
  MonitorElement* hist_EE_CorSC_S25toE_;
  MonitorElement* hist_EE_PreSC_S25toE_;
  MonitorElement* hist_EB_RawSC_EoverTruth_;
  MonitorElement* hist_EE_RawSC_EoverTruth_;
  MonitorElement* hist_EB_CorSC_EoverTruth_;
  MonitorElement* hist_EE_CorSC_EoverTruth_;
  MonitorElement* hist_EE_PreSC_EoverTruth_;
  MonitorElement* hist_EB_RawSC_deltaR_;
  MonitorElement* hist_EE_RawSC_deltaR_;
  MonitorElement* hist_EB_CorSC_deltaR_;
  MonitorElement* hist_EE_CorSC_deltaR_;
  MonitorElement* hist_EE_PreSC_deltaR_;
  MonitorElement* hist_EE_PreSC_preshowerE_;
  MonitorElement* hist_EE_CorSC_preshowerE_;
  MonitorElement* hist_EE_CorSC_phiWidth_;
  MonitorElement* hist_EB_CorSC_phiWidth_;
  MonitorElement* hist_EE_CorSC_etaWidth_;
  MonitorElement* hist_EB_CorSC_etaWidth_;
  MonitorElement* hist_EB_CorSC_ET_vs_Eta_;
  MonitorElement* hist_EB_CorSC_ET_vs_Phi_;
  MonitorElement* hist_EE_CorSC_ET_vs_Eta_;
  MonitorElement* hist_EE_CorSC_ET_vs_Phi_;
  MonitorElement* hist_EE_CorSC_ET_vs_R_;

  void closestMCParticle(HepMC::GenEvent const*, reco::SuperCluster const&,
                         double&, double&) const;

  float ecalEta(float, float, float) const;
};

#endif