EgammaSuperClusters Class Reference

#include <EgammaSuperClusters.h>

Inheritance diagram for EgammaSuperClusters:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	EgammaSuperClusters (const edm::ParameterSet &)
	~EgammaSuperClusters () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	closestMCParticle (HepMC::GenEvent const *, reco::SuperCluster const &, double &, double &) const
float	ecalEta (float, float, float) const

Private Attributes
edm::EDGetTokenT< reco::SuperClusterCollection >	barrelCorSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection >	barrelRawSuperClusterCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection >	barrelRecHitCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection >	endcapCorSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection >	endcapPreSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection >	endcapRawSuperClusterCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection >	endcapRecHitCollectionToken_
MonitorElement *	hist_EB_CorSC_deltaR_
MonitorElement *	hist_EB_CorSC_EoverTruth_
MonitorElement *	hist_EB_CorSC_ET_
MonitorElement *	hist_EB_CorSC_ET_vs_Eta_
MonitorElement *	hist_EB_CorSC_ET_vs_Phi_
MonitorElement *	hist_EB_CorSC_Eta_
MonitorElement *	hist_EB_CorSC_etaWidth_
MonitorElement *	hist_EB_CorSC_NumBC_
MonitorElement *	hist_EB_CorSC_Phi_
MonitorElement *	hist_EB_CorSC_phiWidth_
MonitorElement *	hist_EB_CorSC_S1toS9_
MonitorElement *	hist_EB_CorSC_S25toE_
MonitorElement *	hist_EB_CorSC_Size_
MonitorElement *	hist_EB_RawSC_deltaR_
MonitorElement *	hist_EB_RawSC_EoverTruth_
MonitorElement *	hist_EB_RawSC_ET_
MonitorElement *	hist_EB_RawSC_Eta_
MonitorElement *	hist_EB_RawSC_NumBC_
MonitorElement *	hist_EB_RawSC_Phi_
MonitorElement *	hist_EB_RawSC_S1toS9_
MonitorElement *	hist_EB_RawSC_S25toE_
MonitorElement *	hist_EB_RawSC_Size_
MonitorElement *	hist_EE_CorSC_deltaR_
MonitorElement *	hist_EE_CorSC_EoverTruth_
MonitorElement *	hist_EE_CorSC_ET_
MonitorElement *	hist_EE_CorSC_ET_vs_Eta_
MonitorElement *	hist_EE_CorSC_ET_vs_Phi_
MonitorElement *	hist_EE_CorSC_ET_vs_R_
MonitorElement *	hist_EE_CorSC_Eta_
MonitorElement *	hist_EE_CorSC_etaWidth_
MonitorElement *	hist_EE_CorSC_NumBC_
MonitorElement *	hist_EE_CorSC_Phi_
MonitorElement *	hist_EE_CorSC_phiWidth_
MonitorElement *	hist_EE_CorSC_preshowerE_
MonitorElement *	hist_EE_CorSC_S1toS9_
MonitorElement *	hist_EE_CorSC_S25toE_
MonitorElement *	hist_EE_CorSC_Size_
MonitorElement *	hist_EE_PreSC_deltaR_
MonitorElement *	hist_EE_PreSC_EoverTruth_
MonitorElement *	hist_EE_PreSC_ET_
MonitorElement *	hist_EE_PreSC_Eta_
MonitorElement *	hist_EE_PreSC_NumBC_
MonitorElement *	hist_EE_PreSC_Phi_
MonitorElement *	hist_EE_PreSC_preshowerE_
MonitorElement *	hist_EE_PreSC_S1toS9_
MonitorElement *	hist_EE_PreSC_S25toE_
MonitorElement *	hist_EE_PreSC_Size_
MonitorElement *	hist_EE_RawSC_deltaR_
MonitorElement *	hist_EE_RawSC_EoverTruth_
MonitorElement *	hist_EE_RawSC_ET_
MonitorElement *	hist_EE_RawSC_Eta_
MonitorElement *	hist_EE_RawSC_NumBC_
MonitorElement *	hist_EE_RawSC_Phi_
MonitorElement *	hist_EE_RawSC_S1toS9_
MonitorElement *	hist_EE_RawSC_S25toE_
MonitorElement *	hist_EE_RawSC_Size_
HistSpec	hsdeltaR_
HistSpec	hsEoverTruth_
HistSpec	hsET_
HistSpec	hsEta_
HistSpec	hsetaWidth_
HistSpec	hsNumBC_
HistSpec	hsPhi_
HistSpec	hsphiWidth_
HistSpec	hspreshowerE_
HistSpec	hsR_
HistSpec	hsS1toS9_
HistSpec	hsS25toE_
HistSpec	hsSize_
edm::EDGetTokenT< edm::HepMCProduct >	MCTruthCollectionToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Description: SVSuite Super Cluster Validation

Implementation: \

Author

: Michael A. Balazs, Nov 2006

Definition at line 25 of file EgammaSuperClusters.h.

Constructor & Destructor Documentation

EgammaSuperClusters::EgammaSuperClusters ( const edm::ParameterSet &  ps )

explicit

Definition at line 18 of file EgammaSuperClusters.cc.

    : MCTruthCollectionToken_(consumes<edm::HepMCProduct>(ps.getParameter<edm::InputTag>("MCTruthCollection"))),
      barrelRawSuperClusterCollectionToken_(
          consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelRawSuperClusterCollection"))),
      barrelCorSuperClusterCollectionToken_(
          consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelCorSuperClusterCollection"))),
      endcapRawSuperClusterCollectionToken_(
          consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapRawSuperClusterCollection"))),
      endcapPreSuperClusterCollectionToken_(
          consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapPreSuperClusterCollection"))),
      endcapCorSuperClusterCollectionToken_(
          consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapCorSuperClusterCollection"))),
      barrelRecHitCollectionToken_(
          consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelRecHitCollection"))),
      endcapRecHitCollectionToken_(
          consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapRecHitCollection"))),
      hsSize_(ps, "Size"),
      hsNumBC_(ps, "NumBC"),
      hsET_(ps, "ET"),
      hsEta_(ps, "Eta"),
      hsPhi_(ps, "Phi"),
      hsS1toS9_(ps, "S1toS9"),
      hsS25toE_(ps, "S25toE"),
      hsEoverTruth_(ps, "EoverTruth"),
      hsdeltaR_(ps, "deltaR"),
      hsphiWidth_(ps, "phiWidth"),
      hsetaWidth_(ps, "etaWidth"),
      hspreshowerE_(ps, "preshowerE"),
      hsR_(ps, "R"),
      hist_EB_RawSC_Size_(nullptr),
      hist_EE_RawSC_Size_(nullptr),
      hist_EB_CorSC_Size_(nullptr),
      hist_EE_CorSC_Size_(nullptr),
      hist_EE_PreSC_Size_(nullptr),
      hist_EB_RawSC_NumBC_(nullptr),
      hist_EE_RawSC_NumBC_(nullptr),
      hist_EB_CorSC_NumBC_(nullptr),
      hist_EE_CorSC_NumBC_(nullptr),
      hist_EE_PreSC_NumBC_(nullptr),
      hist_EB_RawSC_ET_(nullptr),
      hist_EE_RawSC_ET_(nullptr),
      hist_EB_CorSC_ET_(nullptr),
      hist_EE_CorSC_ET_(nullptr),
      hist_EE_PreSC_ET_(nullptr),
      hist_EB_RawSC_Eta_(nullptr),
      hist_EE_RawSC_Eta_(nullptr),
      hist_EB_CorSC_Eta_(nullptr),
      hist_EE_CorSC_Eta_(nullptr),
      hist_EE_PreSC_Eta_(nullptr),
      hist_EB_RawSC_Phi_(nullptr),
      hist_EE_RawSC_Phi_(nullptr),
      hist_EB_CorSC_Phi_(nullptr),
      hist_EE_CorSC_Phi_(nullptr),
      hist_EE_PreSC_Phi_(nullptr),
      hist_EB_RawSC_S1toS9_(nullptr),
      hist_EE_RawSC_S1toS9_(nullptr),
      hist_EB_CorSC_S1toS9_(nullptr),
      hist_EE_CorSC_S1toS9_(nullptr),
      hist_EE_PreSC_S1toS9_(nullptr),
      hist_EB_RawSC_S25toE_(nullptr),
      hist_EE_RawSC_S25toE_(nullptr),
      hist_EB_CorSC_S25toE_(nullptr),
      hist_EE_CorSC_S25toE_(nullptr),
      hist_EE_PreSC_S25toE_(nullptr),
      hist_EB_RawSC_EoverTruth_(nullptr),
      hist_EE_RawSC_EoverTruth_(nullptr),
      hist_EB_CorSC_EoverTruth_(nullptr),
      hist_EE_CorSC_EoverTruth_(nullptr),
      hist_EE_PreSC_EoverTruth_(nullptr),
      hist_EB_RawSC_deltaR_(nullptr),
      hist_EE_RawSC_deltaR_(nullptr),
      hist_EB_CorSC_deltaR_(nullptr),
      hist_EE_CorSC_deltaR_(nullptr),
      hist_EE_PreSC_deltaR_(nullptr),
      hist_EE_PreSC_preshowerE_(nullptr),
      hist_EE_CorSC_preshowerE_(nullptr),
      hist_EE_CorSC_phiWidth_(nullptr),
      hist_EB_CorSC_phiWidth_(nullptr),
      hist_EE_CorSC_etaWidth_(nullptr),
      hist_EB_CorSC_etaWidth_(nullptr),
      hist_EB_CorSC_ET_vs_Eta_(nullptr),
      hist_EB_CorSC_ET_vs_Phi_(nullptr),
      hist_EE_CorSC_ET_vs_Eta_(nullptr),
      hist_EE_CorSC_ET_vs_Phi_(nullptr),
      hist_EE_CorSC_ET_vs_R_(nullptr) {}

EgammaSuperClusters::~EgammaSuperClusters ( )

override

Definition at line 104 of file EgammaSuperClusters.cc.

{}

Member Function Documentation

void EgammaSuperClusters::analyze ( const edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 372 of file EgammaSuperClusters.cc.

References barrelCorSuperClusterCollectionToken_, barrelRawSuperClusterCollectionToken_, barrelRecHitCollectionToken_, closestMCParticle(), endcapCorSuperClusterCollectionToken_, endcapPreSuperClusterCollectionToken_, endcapRawSuperClusterCollectionToken_, endcapRecHitCollectionToken_, dqm::impl::MonitorElement::Fill(), nano_cff::genEvent, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hist_EB_CorSC_deltaR_, hist_EB_CorSC_EoverTruth_, hist_EB_CorSC_ET_, hist_EB_CorSC_ET_vs_Eta_, hist_EB_CorSC_ET_vs_Phi_, hist_EB_CorSC_Eta_, hist_EB_CorSC_etaWidth_, hist_EB_CorSC_NumBC_, hist_EB_CorSC_Phi_, hist_EB_CorSC_phiWidth_, hist_EB_CorSC_S1toS9_, hist_EB_CorSC_S25toE_, hist_EB_CorSC_Size_, hist_EB_RawSC_deltaR_, hist_EB_RawSC_EoverTruth_, hist_EB_RawSC_ET_, hist_EB_RawSC_Eta_, hist_EB_RawSC_NumBC_, hist_EB_RawSC_Phi_, hist_EB_RawSC_S1toS9_, hist_EB_RawSC_S25toE_, hist_EB_RawSC_Size_, hist_EE_CorSC_deltaR_, hist_EE_CorSC_EoverTruth_, hist_EE_CorSC_ET_, hist_EE_CorSC_ET_vs_Eta_, hist_EE_CorSC_ET_vs_Phi_, hist_EE_CorSC_ET_vs_R_, hist_EE_CorSC_Eta_, hist_EE_CorSC_etaWidth_, hist_EE_CorSC_NumBC_, hist_EE_CorSC_Phi_, hist_EE_CorSC_phiWidth_, hist_EE_CorSC_preshowerE_, hist_EE_CorSC_S1toS9_, hist_EE_CorSC_S25toE_, hist_EE_CorSC_Size_, hist_EE_PreSC_deltaR_, hist_EE_PreSC_EoverTruth_, hist_EE_PreSC_ET_, hist_EE_PreSC_Eta_, hist_EE_PreSC_NumBC_, hist_EE_PreSC_Phi_, hist_EE_PreSC_preshowerE_, hist_EE_PreSC_S1toS9_, hist_EE_PreSC_S25toE_, hist_EE_PreSC_Size_, hist_EE_RawSC_deltaR_, hist_EE_RawSC_EoverTruth_, hist_EE_RawSC_ET_, hist_EE_RawSC_Eta_, hist_EE_RawSC_NumBC_, hist_EE_RawSC_Phi_, hist_EE_RawSC_S1toS9_, hist_EE_RawSC_S25toE_, hist_EE_RawSC_Size_, edm::HandleBase::isValid(), MCTruthCollectionToken_, funct::pow(), edm::Handle< T >::product(), SurveyInfoScenario_cff::seed, and mathSSE::sqrt().

                                                                              {
  bool skipMC = false;
  bool skipBarrel = false;
  bool skipEndcap = false;

  //
  // Get MCTRUTH
  //
  edm::Handle<edm::HepMCProduct> pMCTruth;
  evt.getByToken(MCTruthCollectionToken_, pMCTruth);
  if (!pMCTruth.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get MC collection ";
    skipMC = true;
  }
  const HepMC::GenEvent *genEvent = pMCTruth->GetEvent();

  if (skipMC)
    return;

  //
  // Get the BARREL products
  //
  edm::Handle<reco::SuperClusterCollection> pBarrelRawSuperClusters;
  evt.getByToken(barrelRawSuperClusterCollectionToken_, pBarrelRawSuperClusters);
  if (!pBarrelRawSuperClusters.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw SC";
    skipBarrel = true;
  }

  edm::Handle<reco::SuperClusterCollection> pBarrelCorSuperClusters;
  evt.getByToken(barrelCorSuperClusterCollectionToken_, pBarrelCorSuperClusters);
  if (!pBarrelCorSuperClusters.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor SC";
    skipBarrel = true;
  }

  edm::Handle<EBRecHitCollection> pBarrelRecHitCollection;
  evt.getByToken(barrelRecHitCollectionToken_, pBarrelRecHitCollection);
  if (!pBarrelRecHitCollection.isValid()) {
    skipBarrel = true;
  }
  edm::Handle<EERecHitCollection> pEndcapRecHitCollection;
  evt.getByToken(endcapRecHitCollectionToken_, pEndcapRecHitCollection);
  if (!pEndcapRecHitCollection.isValid()) {
    skipEndcap = true;
  }

  if (skipBarrel || skipEndcap)
    return;

  EcalClusterLazyTools lazyTool(evt, es, barrelRecHitCollectionToken_, endcapRecHitCollectionToken_);

  // Get the BARREL collections
  const reco::SuperClusterCollection *barrelRawSuperClusters = pBarrelRawSuperClusters.product();
  const reco::SuperClusterCollection *barrelCorSuperClusters = pBarrelCorSuperClusters.product();

  // Number of entries in collections
  hist_EB_RawSC_Size_->Fill(barrelRawSuperClusters->size());
  hist_EB_CorSC_Size_->Fill(barrelCorSuperClusters->size());

  // Do RAW BARREL SuperClusters
  for (reco::SuperClusterCollection::const_iterator aClus = barrelRawSuperClusters->begin();
       aClus != barrelRawSuperClusters->end();
       aClus++) {
    // kinematics
    hist_EB_RawSC_NumBC_->Fill(aClus->clustersSize());
    hist_EB_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
    hist_EB_RawSC_Eta_->Fill(aClus->position().eta());
    hist_EB_RawSC_Phi_->Fill(aClus->position().phi());

    // cluster shape
    const reco::CaloClusterPtr seed = aClus->seed();
    hist_EB_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
    hist_EB_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());

    // truth
    double dRClosest = 999.9;
    double energyClosest = 0;
    closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);

    if (dRClosest < 0.1) {
      hist_EB_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
      hist_EB_RawSC_deltaR_->Fill(dRClosest);
    }
  }

  // Do CORRECTED BARREL SuperClusters
  for (reco::SuperClusterCollection::const_iterator aClus = barrelCorSuperClusters->begin();
       aClus != barrelCorSuperClusters->end();
       aClus++) {
    // kinematics
    hist_EB_CorSC_NumBC_->Fill(aClus->clustersSize());
    hist_EB_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
    hist_EB_CorSC_Eta_->Fill(aClus->position().eta());
    hist_EB_CorSC_Phi_->Fill(aClus->position().phi());

    hist_EB_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
    hist_EB_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());

    // cluster shape
    const reco::CaloClusterPtr seed = aClus->seed();
    hist_EB_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
    hist_EB_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());

    // correction variables
    hist_EB_CorSC_phiWidth_->Fill(aClus->phiWidth());
    hist_EB_CorSC_etaWidth_->Fill(aClus->etaWidth());

    // truth
    double dRClosest = 999.9;
    double energyClosest = 0;
    closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);

    if (dRClosest < 0.1) {
      hist_EB_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
      hist_EB_CorSC_deltaR_->Fill(dRClosest);
    }
  }

  //
  // Get the ENDCAP products
  //
  edm::Handle<reco::SuperClusterCollection> pEndcapRawSuperClusters;
  evt.getByToken(endcapRawSuperClusterCollectionToken_, pEndcapRawSuperClusters);
  if (!pEndcapRawSuperClusters.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw EE SC";
  }

  edm::Handle<reco::SuperClusterCollection> pEndcapPreSuperClusters;
  evt.getByToken(endcapPreSuperClusterCollectionToken_, pEndcapPreSuperClusters);
  if (!pEndcapPreSuperClusters.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get collection Pre EE SC";
  }

  edm::Handle<reco::SuperClusterCollection> pEndcapCorSuperClusters;
  evt.getByToken(endcapCorSuperClusterCollectionToken_, pEndcapCorSuperClusters);
  if (!pEndcapCorSuperClusters.isValid()) {
    edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor EE SC";
  }

  // Get the ENDCAP collections
  const reco::SuperClusterCollection *endcapRawSuperClusters = pEndcapRawSuperClusters.product();
  const reco::SuperClusterCollection *endcapPreSuperClusters = pEndcapPreSuperClusters.product();
  const reco::SuperClusterCollection *endcapCorSuperClusters = pEndcapCorSuperClusters.product();

  // Number of entries in collections
  hist_EE_RawSC_Size_->Fill(endcapRawSuperClusters->size());
  hist_EE_PreSC_Size_->Fill(endcapPreSuperClusters->size());
  hist_EE_CorSC_Size_->Fill(endcapCorSuperClusters->size());

  // Do RAW ENDCAP SuperClusters
  for (reco::SuperClusterCollection::const_iterator aClus = endcapRawSuperClusters->begin();
       aClus != endcapRawSuperClusters->end();
       aClus++) {
    hist_EE_RawSC_NumBC_->Fill(aClus->clustersSize());
    hist_EE_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
    hist_EE_RawSC_Eta_->Fill(aClus->position().eta());
    hist_EE_RawSC_Phi_->Fill(aClus->position().phi());

    const reco::CaloClusterPtr seed = aClus->seed();
    hist_EE_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
    hist_EE_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());

    // truth
    double dRClosest = 999.9;
    double energyClosest = 0;
    closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);

    if (dRClosest < 0.1) {
      hist_EE_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
      hist_EE_RawSC_deltaR_->Fill(dRClosest);
    }
  }

  // Do ENDCAP SuperClusters with PRESHOWER
  for (reco::SuperClusterCollection::const_iterator aClus = endcapPreSuperClusters->begin();
       aClus != endcapPreSuperClusters->end();
       aClus++) {
    hist_EE_PreSC_NumBC_->Fill(aClus->clustersSize());
    hist_EE_PreSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
    hist_EE_PreSC_Eta_->Fill(aClus->position().eta());
    hist_EE_PreSC_Phi_->Fill(aClus->position().phi());
    hist_EE_PreSC_preshowerE_->Fill(aClus->preshowerEnergy());

    const reco::CaloClusterPtr seed = aClus->seed();
    hist_EE_PreSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
    hist_EE_PreSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());

    // truth
    double dRClosest = 999.9;
    double energyClosest = 0;
    closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);

    if (dRClosest < 0.1) {
      hist_EE_PreSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
      hist_EE_PreSC_deltaR_->Fill(dRClosest);
    }
  }

  // Do CORRECTED ENDCAP SuperClusters
  for (reco::SuperClusterCollection::const_iterator aClus = endcapCorSuperClusters->begin();
       aClus != endcapCorSuperClusters->end();
       aClus++) {
    hist_EE_CorSC_NumBC_->Fill(aClus->clustersSize());
    hist_EE_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
    hist_EE_CorSC_Eta_->Fill(aClus->position().eta());
    hist_EE_CorSC_Phi_->Fill(aClus->position().phi());
    hist_EE_CorSC_preshowerE_->Fill(aClus->preshowerEnergy());

    hist_EE_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
    hist_EE_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());
    hist_EE_CorSC_ET_vs_R_->Fill(aClus->energy() / std::cosh(aClus->position().eta()),
                                 std::sqrt(std::pow(aClus->x(), 2) + std::pow(aClus->y(), 2)));

    // correction variables
    hist_EE_CorSC_phiWidth_->Fill(aClus->phiWidth());
    hist_EE_CorSC_etaWidth_->Fill(aClus->etaWidth());

    const reco::CaloClusterPtr seed = aClus->seed();
    hist_EE_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
    hist_EE_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());

    // truth
    double dRClosest = 999.9;
    double energyClosest = 0;
    closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);

    if (dRClosest < 0.1) {
      hist_EE_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
      hist_EE_CorSC_deltaR_->Fill(dRClosest);
    }
  }
}

void EgammaSuperClusters::bookHistograms ( DQMStore::IBooker &  _ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 106 of file EgammaSuperClusters.cc.

References HistSpec::bins, dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), hist_EB_CorSC_deltaR_, hist_EB_CorSC_EoverTruth_, hist_EB_CorSC_ET_, hist_EB_CorSC_ET_vs_Eta_, hist_EB_CorSC_ET_vs_Phi_, hist_EB_CorSC_Eta_, hist_EB_CorSC_etaWidth_, hist_EB_CorSC_NumBC_, hist_EB_CorSC_Phi_, hist_EB_CorSC_phiWidth_, hist_EB_CorSC_S1toS9_, hist_EB_CorSC_S25toE_, hist_EB_CorSC_Size_, hist_EB_RawSC_deltaR_, hist_EB_RawSC_EoverTruth_, hist_EB_RawSC_ET_, hist_EB_RawSC_Eta_, hist_EB_RawSC_NumBC_, hist_EB_RawSC_Phi_, hist_EB_RawSC_S1toS9_, hist_EB_RawSC_S25toE_, hist_EB_RawSC_Size_, hist_EE_CorSC_deltaR_, hist_EE_CorSC_EoverTruth_, hist_EE_CorSC_ET_, hist_EE_CorSC_ET_vs_Eta_, hist_EE_CorSC_ET_vs_Phi_, hist_EE_CorSC_ET_vs_R_, hist_EE_CorSC_Eta_, hist_EE_CorSC_etaWidth_, hist_EE_CorSC_NumBC_, hist_EE_CorSC_Phi_, hist_EE_CorSC_phiWidth_, hist_EE_CorSC_preshowerE_, hist_EE_CorSC_S1toS9_, hist_EE_CorSC_S25toE_, hist_EE_CorSC_Size_, hist_EE_PreSC_deltaR_, hist_EE_PreSC_EoverTruth_, hist_EE_PreSC_ET_, hist_EE_PreSC_Eta_, hist_EE_PreSC_NumBC_, hist_EE_PreSC_Phi_, hist_EE_PreSC_preshowerE_, hist_EE_PreSC_S1toS9_, hist_EE_PreSC_S25toE_, hist_EE_PreSC_Size_, hist_EE_RawSC_deltaR_, hist_EE_RawSC_EoverTruth_, hist_EE_RawSC_ET_, hist_EE_RawSC_Eta_, hist_EE_RawSC_NumBC_, hist_EE_RawSC_Phi_, hist_EE_RawSC_S1toS9_, hist_EE_RawSC_S25toE_, hist_EE_RawSC_Size_, hsdeltaR_, hsEoverTruth_, hsET_, hsEta_, hsetaWidth_, hsNumBC_, hsPhi_, hsphiWidth_, hspreshowerE_, hsR_, hsS1toS9_, hsS25toE_, hsSize_, HistSpec::max, HistSpec::min, and dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder().

                                                                                                         {
  _ibooker.setCurrentFolder("EcalClusterV/EcalSuperClusters/");

  // Number of SuperClusters
  //
  hist_EB_RawSC_Size_ =
      _ibooker.book1D("hist_EB_RawSC_Size_", "# Raw SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
  hist_EE_RawSC_Size_ =
      _ibooker.book1D("hist_EE_RawSC_Size_", "# Raw SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
  hist_EB_CorSC_Size_ = _ibooker.book1D(
      "hist_EB_CorSC_Size_", "# Corrected SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
  hist_EE_CorSC_Size_ = _ibooker.book1D(
      "hist_EE_CorSC_Size_", "# Corrected SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
  hist_EE_PreSC_Size_ = _ibooker.book1D(
      "hist_EE_PreSC_Size_", "# SuperClusters with Preshower in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);

  // Number of BasicClusters in SuperCluster
  //
  hist_EB_RawSC_NumBC_ = _ibooker.book1D("hist_EB_RawSC_NumBC_",
                                         "# of Basic Clusters in Raw Super Clusters in Barrel",
                                         hsNumBC_.bins,
                                         hsNumBC_.min,
                                         hsNumBC_.max);
  hist_EE_RawSC_NumBC_ = _ibooker.book1D("hist_EE_RawSC_NumBC_",
                                         "# of Basic Clusters in Raw Super Clusters in Endcap",
                                         hsNumBC_.bins,
                                         hsNumBC_.min,
                                         hsNumBC_.max);
  hist_EB_CorSC_NumBC_ = _ibooker.book1D("hist_EB_CorSC_NumBC_",
                                         "# of Basic Clusters in Corrected SuperClusters in Barrel",
                                         hsNumBC_.bins,
                                         hsNumBC_.min,
                                         hsNumBC_.max);
  hist_EE_CorSC_NumBC_ = _ibooker.book1D("hist_EE_CorSC_NumBC_",
                                         "# of Basic Clusters in Corrected SuperClusters in Endcap",
                                         hsNumBC_.bins,
                                         hsNumBC_.min,
                                         hsNumBC_.max);
  hist_EE_PreSC_NumBC_ = _ibooker.book1D("hist_EE_PreSC_NumBC_",
                                         "# of Basic Clusters in SuperClusters with Preshower in Endcap",
                                         hsNumBC_.bins,
                                         hsNumBC_.min,
                                         hsNumBC_.max);

  // ET distribution of SuperClusters
  //
  hist_EB_RawSC_ET_ =
      _ibooker.book1D("hist_EB_RawSC_ET_", "ET of Raw SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
  hist_EE_RawSC_ET_ =
      _ibooker.book1D("hist_EE_RawSC_ET_", "ET of Raw SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
  hist_EB_CorSC_ET_ =
      _ibooker.book1D("hist_EB_CorSC_ET_", "ET of Corrected SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
  hist_EE_CorSC_ET_ =
      _ibooker.book1D("hist_EE_CorSC_ET_", "ET of Corrected SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
  hist_EE_PreSC_ET_ = _ibooker.book1D(
      "hist_EE_PreSC_ET_", "ET of SuperClusters with Preshower in Endcap", hsET_.bins, hsET_.min, hsET_.max);

  // Eta distribution of SuperClusters
  //
  hist_EB_RawSC_Eta_ =
      _ibooker.book1D("hist_EB_RawSC_Eta_", "Eta of Raw SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
  hist_EE_RawSC_Eta_ =
      _ibooker.book1D("hist_EE_RawSC_Eta_", "Eta of Raw SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
  hist_EB_CorSC_Eta_ = _ibooker.book1D(
      "hist_EB_CorSC_Eta_", "Eta of Corrected SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
  hist_EE_CorSC_Eta_ = _ibooker.book1D(
      "hist_EE_CorSC_Eta_", "Eta of Corrected SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
  hist_EE_PreSC_Eta_ = _ibooker.book1D(
      "hist_EE_PreSC_Eta_", "Eta of SuperClusters with Preshower in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);

  // Phi distribution of SuperClusters
  //
  hist_EB_RawSC_Phi_ =
      _ibooker.book1D("hist_EB_RawSC_Phi_", "Phi of Raw SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
  hist_EE_RawSC_Phi_ =
      _ibooker.book1D("hist_EE_RawSC_Phi_", "Phi of Raw SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
  hist_EB_CorSC_Phi_ = _ibooker.book1D(
      "hist_EB_CorSC_Phi_", "Phi of Corrected SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
  hist_EE_CorSC_Phi_ = _ibooker.book1D(
      "hist_EE_CorSC_Phi_", "Phi of Corrected SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
  hist_EE_PreSC_Phi_ = _ibooker.book1D(
      "hist_EE_PreSC_Phi_", "Phi of SuperClusters with Preshower in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);

  // S1/S9 distribution of SuperClusters
  //
  hist_EB_RawSC_S1toS9_ = _ibooker.book1D(
      "hist_EB_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Barrel", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
  hist_EE_RawSC_S1toS9_ = _ibooker.book1D(
      "hist_EE_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Endcap", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
  hist_EB_CorSC_S1toS9_ = _ibooker.book1D("hist_EB_CorSC_S1toS9_",
                                          "S1/S9 of Corrected SuperClusters in Barrel",
                                          hsS1toS9_.bins,
                                          hsS1toS9_.min,
                                          hsS1toS9_.max);
  hist_EE_CorSC_S1toS9_ = _ibooker.book1D("hist_EE_CorSC_S1toS9_",
                                          "S1/S9 of Corrected SuperClusters in Endcap",
                                          hsS1toS9_.bins,
                                          hsS1toS9_.min,
                                          hsS1toS9_.max);
  hist_EE_PreSC_S1toS9_ = _ibooker.book1D("hist_EE_PreSC_S1toS9_",
                                          "S1/S9 of SuperClusters with Preshower in Endcap",
                                          hsS1toS9_.bins,
                                          hsS1toS9_.min,
                                          hsS1toS9_.max);

  // S25/E distribution of SuperClusters
  //
  hist_EB_RawSC_S25toE_ = _ibooker.book1D(
      "hist_EB_RawSC_S25toE_", "S25/E of Raw Super Clusters in Barrel", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
  hist_EE_RawSC_S25toE_ = _ibooker.book1D(
      "hist_EE_RawSC_S25toE_", "S25/E of Raw Super Clusters in Endcap", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
  hist_EB_CorSC_S25toE_ = _ibooker.book1D("hist_EB_CorSC_S25toE_",
                                          "S25/E of Corrected SuperClusters in Barrel",
                                          hsS25toE_.bins,
                                          hsS25toE_.min,
                                          hsS25toE_.max);
  hist_EE_CorSC_S25toE_ = _ibooker.book1D("hist_EE_CorSC_S25toE_",
                                          "S25/E of Corrected SuperClusters in Endcap",
                                          hsS25toE_.bins,
                                          hsS25toE_.min,
                                          hsS25toE_.max);
  hist_EE_PreSC_S25toE_ = _ibooker.book1D("hist_EE_PreSC_S25toE_",
                                          "S25/E of SuperClusters with Preshower in Endcap",
                                          hsS25toE_.bins,
                                          hsS25toE_.min,
                                          hsS25toE_.max);

  // E/E(true) distribution of SuperClusters
  //
  hist_EB_RawSC_EoverTruth_ = _ibooker.book1D("hist_EB_RawSC_EoverTruth_",
                                              "E/True E of Raw SuperClusters in Barrel",
                                              hsEoverTruth_.bins,
                                              hsEoverTruth_.min,
                                              hsEoverTruth_.max);
  hist_EE_RawSC_EoverTruth_ = _ibooker.book1D("hist_EE_RawSC_EoverTruth_",
                                              "E/True E of Raw SuperClusters in Endcap",
                                              hsEoverTruth_.bins,
                                              hsEoverTruth_.min,
                                              hsEoverTruth_.max);
  hist_EB_CorSC_EoverTruth_ = _ibooker.book1D("hist_EB_CorSC_EoverTruth_",
                                              "E/True E of Corrected SuperClusters in Barrel",
                                              hsEoverTruth_.bins,
                                              hsEoverTruth_.min,
                                              hsEoverTruth_.max);
  hist_EE_CorSC_EoverTruth_ = _ibooker.book1D("hist_EE_CorSC_EoverTruth_",
                                              "E/True E of Corrected SuperClusters in Endcap",
                                              hsEoverTruth_.bins,
                                              hsEoverTruth_.min,
                                              hsEoverTruth_.max);
  hist_EE_PreSC_EoverTruth_ = _ibooker.book1D("hist_EE_PreSC_EoverTruth_",
                                              "E/True E of SuperClusters with Preshower in Endcap",
                                              hsEoverTruth_.bins,
                                              hsEoverTruth_.min,
                                              hsEoverTruth_.max);

  // dR distribution of SuperClusters from truth
  //
  hist_EB_RawSC_deltaR_ = _ibooker.book1D("hist_EB_RawSC_deltaR_",
                                          "dR to MC truth of Raw Super Clusters in Barrel",
                                          hsdeltaR_.bins,
                                          hsdeltaR_.min,
                                          hsdeltaR_.max);
  hist_EE_RawSC_deltaR_ = _ibooker.book1D("hist_EE_RawSC_deltaR_",
                                          "dR to MC truth of Raw Super Clusters in Endcap",
                                          hsdeltaR_.bins,
                                          hsdeltaR_.min,
                                          hsdeltaR_.max);
  hist_EB_CorSC_deltaR_ = _ibooker.book1D("hist_EB_CorSC_deltaR_",
                                          "dR to MC truth of Corrected SuperClusters in Barrel",
                                          hsdeltaR_.bins,
                                          hsdeltaR_.min,
                                          hsdeltaR_.max);
  hist_EE_CorSC_deltaR_ = _ibooker.book1D("hist_EE_CorSC_deltaR_",
                                          "dR to MC truth of Corrected SuperClusters in Endcap",
                                          hsdeltaR_.bins,
                                          hsdeltaR_.min,
                                          hsdeltaR_.max);
  hist_EE_PreSC_deltaR_ = _ibooker.book1D("hist_EE_PreSC_deltaR_",
                                          "dR to MC truth of SuperClusters with Preshower in Endcap",
                                          hsdeltaR_.bins,
                                          hsdeltaR_.min,
                                          hsdeltaR_.max);

  // phi width stored in corrected SuperClusters
  hist_EB_CorSC_phiWidth_ = _ibooker.book1D("hist_EB_CorSC_phiWidth_",
                                            "phiWidth of Corrected Super Clusters in Barrel",
                                            hsphiWidth_.bins,
                                            hsphiWidth_.min,
                                            hsphiWidth_.max);
  hist_EE_CorSC_phiWidth_ = _ibooker.book1D("hist_EE_CorSC_phiWidth_",
                                            "phiWidth of Corrected Super Clusters in Endcap",
                                            hsphiWidth_.bins,
                                            hsphiWidth_.min,
                                            hsphiWidth_.max);

  // eta width stored in corrected SuperClusters
  hist_EB_CorSC_etaWidth_ = _ibooker.book1D("hist_EB_CorSC_etaWidth_",
                                            "etaWidth of Corrected Super Clusters in Barrel",
                                            hsetaWidth_.bins,
                                            hsetaWidth_.min,
                                            hsetaWidth_.max);
  hist_EE_CorSC_etaWidth_ = _ibooker.book1D("hist_EE_CorSC_etaWidth_",
                                            "etaWidth of Corrected Super Clusters in Endcap",
                                            hsetaWidth_.bins,
                                            hsetaWidth_.min,
                                            hsetaWidth_.max);

  // preshower energy
  hist_EE_PreSC_preshowerE_ = _ibooker.book1D("hist_EE_PreSC_preshowerE_",
                                              "preshower energy in Super Clusters with Preshower in Endcap",
                                              hspreshowerE_.bins,
                                              hspreshowerE_.min,
                                              hspreshowerE_.max);
  hist_EE_CorSC_preshowerE_ = _ibooker.book1D("hist_EE_CorSC_preshowerE_",
                                              "preshower energy in Corrected Super Clusters with Preshower in Endcap",
                                              hspreshowerE_.bins,
                                              hspreshowerE_.min,
                                              hspreshowerE_.max);

  //
  hist_EB_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Eta_",
                                             "Corr Super Cluster ET versus Eta in Barrel",
                                             hsET_.bins,
                                             hsET_.min,
                                             hsET_.max,
                                             hsEta_.bins,
                                             hsEta_.min,
                                             hsEta_.max);

  hist_EB_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Phi_",
                                             "Corr Super Cluster ET versus Phi in Barrel",
                                             hsET_.bins,
                                             hsET_.min,
                                             hsET_.max,
                                             hsPhi_.bins,
                                             hsPhi_.min,
                                             hsPhi_.max);

  hist_EE_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Eta_",
                                             "Corr Super Cluster ET versus Eta in Endcap",
                                             hsET_.bins,
                                             hsET_.min,
                                             hsET_.max,
                                             hsEta_.bins,
                                             hsEta_.min,
                                             hsEta_.max);

  hist_EE_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Phi_",
                                             "Corr Super Cluster ET versus Phi in Endcap",
                                             hsET_.bins,
                                             hsET_.min,
                                             hsET_.max,
                                             hsPhi_.bins,
                                             hsPhi_.min,
                                             hsPhi_.max);

  hist_EE_CorSC_ET_vs_R_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_R_",
                                           "Corr Super Cluster ET versus Radius in Endcap",
                                           hsET_.bins,
                                           hsET_.min,
                                           hsET_.max,
                                           hsR_.bins,
                                           hsR_.min,
                                           hsR_.max);
}

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

private

Definition at line 609 of file EgammaSuperClusters.cc.

References reco::deltaPhi(), PbPb_ZMuSkimMuonDPG_cff::deltaR, HLT_2018_cff::dEta, HLT_2018_cff::dPhi, ecalEta(), reco::CaloCluster::eta(), reco::CaloCluster::phi(), mathSSE::sqrt(), and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by analyze().

                                                                         {
  // SuperCluster eta, phi
  double scEta = sc.eta();
  double scPhi = sc.phi();

  // initialize dRClosest to a large number
  dRClosest = 999.9;

  // loop over the MC truth particles to find the
  // closest to the superCluster in dR space
  for (HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
       currentParticle != genEvent->particles_end();
       currentParticle++) {
    if ((*currentParticle)->status() == 1) {
      // need GenParticle in ECAL co-ordinates
      HepMC::FourVector vtx = (*currentParticle)->production_vertex()->position();
      double phiTrue = (*currentParticle)->momentum().phi();
      double etaTrue = ecalEta((*currentParticle)->momentum().eta(), vtx.z() / 10., vtx.perp() / 10.);

      double dPhi = reco::deltaPhi(phiTrue, scPhi);
      double dEta = scEta - etaTrue;
      double deltaR = std::sqrt(dPhi * dPhi + dEta * dEta);

      if (deltaR < dRClosest) {
        dRClosest = deltaR;
        energyClosest = (*currentParticle)->momentum().e();
      }

    }  // end if stable particle

  }  // end loop on get particles
}

float EgammaSuperClusters::ecalEta ( float  EtaParticle, float  Zvertex, float  plane_Radius  ) const

private

Definition at line 648 of file EgammaSuperClusters.cc.

References ETA, etaBarrelEndcap, dqm-mbProfile::log, Geom::pi(), R_ECAL, funct::tan(), and Z_Endcap.

Referenced by closestMCParticle().

                                                                                             {
  const float R_ECAL = 136.5;
  const float Z_Endcap = 328.0;
  const float etaBarrelEndcap = 1.479;

  if (EtaParticle != 0.) {
    float Theta = 0.0;
    float ZEcal = (R_ECAL - plane_Radius) * sinh(EtaParticle) + Zvertex;

    if (ZEcal != 0.0)
      Theta = atan(R_ECAL / ZEcal);
    if (Theta < 0.0)
      Theta = Theta + Geom::pi();

    float ETA = -log(tan(0.5 * Theta));

    if (fabs(ETA) > etaBarrelEndcap) {
      float Zend = Z_Endcap;
      if (EtaParticle < 0.0)
        Zend = -Zend;
      float Zlen = Zend - Zvertex;
      float RR = Zlen / sinh(EtaParticle);
      Theta = atan((RR + plane_Radius) / Zend);
      if (Theta < 0.0)
        Theta = Theta + Geom::pi();
      ETA = -log(tan(0.5 * Theta));
    }

    return ETA;
  } else {
    edm::LogWarning("") << "[EgammaSuperClusters::ecalEta] Warning: Eta equals "
                           "to zero, not correcting";
    return EtaParticle;
  }
}

Member Data Documentation

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::barrelCorSuperClusterCollectionToken_

private

Definition at line 39 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::barrelRawSuperClusterCollectionToken_

private

Definition at line 38 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::barrelRecHitCollectionToken_

private

Definition at line 47 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapCorSuperClusterCollectionToken_

private

Definition at line 44 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapPreSuperClusterCollectionToken_

private

Definition at line 43 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapRawSuperClusterCollectionToken_

private

Definition at line 42 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::endcapRecHitCollectionToken_

private

Definition at line 48 of file EgammaSuperClusters.h.

Referenced by analyze().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_deltaR_

private

Definition at line 106 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_EoverTruth_

private

Definition at line 101 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_

private

Definition at line 76 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Eta_

private

Definition at line 115 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Phi_

private

Definition at line 116 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Eta_

private

Definition at line 81 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_etaWidth_

private

Definition at line 114 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_NumBC_

private

Definition at line 71 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Phi_

private

Definition at line 86 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_phiWidth_

private

Definition at line 112 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S1toS9_

private

Definition at line 91 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S25toE_

private

Definition at line 96 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Size_

private

Definition at line 66 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_deltaR_

private

Definition at line 104 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_EoverTruth_

private

Definition at line 99 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_ET_

private

Definition at line 74 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Eta_

private

Definition at line 79 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_NumBC_

private

Definition at line 69 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Phi_

private

Definition at line 84 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S1toS9_

private

Definition at line 89 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S25toE_

private

Definition at line 94 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Size_

private

Definition at line 64 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_deltaR_

private

Definition at line 107 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_EoverTruth_

private

Definition at line 102 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_

private

Definition at line 77 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Eta_

private

Definition at line 117 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Phi_

private

Definition at line 118 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_R_

private

Definition at line 119 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Eta_

private

Definition at line 82 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_etaWidth_

private

Definition at line 113 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_NumBC_

private

Definition at line 72 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Phi_

private

Definition at line 87 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_phiWidth_

private

Definition at line 111 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_preshowerE_

private

Definition at line 110 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S1toS9_

private

Definition at line 92 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S25toE_

private

Definition at line 97 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Size_

private

Definition at line 67 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_deltaR_

private

Definition at line 108 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_EoverTruth_

private

Definition at line 103 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_ET_

private

Definition at line 78 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Eta_

private

Definition at line 83 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_NumBC_

private

Definition at line 73 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Phi_

private

Definition at line 88 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_preshowerE_

private

Definition at line 109 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S1toS9_

private

Definition at line 93 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S25toE_

private

Definition at line 98 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Size_

private

Definition at line 68 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_deltaR_

private

Definition at line 105 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_EoverTruth_

private

Definition at line 100 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_ET_

private

Definition at line 75 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Eta_

private

Definition at line 80 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_NumBC_

private

Definition at line 70 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Phi_

private

Definition at line 85 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S1toS9_

private

Definition at line 90 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S25toE_

private

Definition at line 95 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Size_

private

Definition at line 65 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

HistSpec EgammaSuperClusters::hsdeltaR_

private

Definition at line 58 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsEoverTruth_

private

Definition at line 57 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsET_

private

Definition at line 52 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsEta_

private

Definition at line 53 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsetaWidth_

private

Definition at line 60 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsNumBC_

private

Definition at line 51 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsPhi_

private

Definition at line 54 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsphiWidth_

private

Definition at line 59 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hspreshowerE_

private

Definition at line 61 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsR_

private

Definition at line 62 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsS1toS9_

private

Definition at line 55 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsS25toE_

private

Definition at line 56 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsSize_

private

Definition at line 50 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

edm::EDGetTokenT<edm::HepMCProduct> EgammaSuperClusters::MCTruthCollectionToken_

private

Definition at line 35 of file EgammaSuperClusters.h.

Referenced by analyze().