#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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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_

EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_

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::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: SVSuite Super Cluster Validation

Implementation: \

Author: : Michael A. Balazs, Nov 2006

Definition at line 27 of file EgammaSuperClusters.h.

Constructor & Destructor Documentation

◆ EgammaSuperClusters()

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

explicit

Definition at line 16 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"))),
       ecalClusterToolsESGetTokens_{consumesCollector()},
       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::~EgammaSuperClusters ( )

override

Definition at line 103 of file EgammaSuperClusters.cc.

103 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 371 of file EgammaSuperClusters.cc.

References barrelCorSuperClusterCollectionToken_, barrelRawSuperClusterCollectionToken_, barrelRecHitCollectionToken_, closestMCParticle(), ecalClusterToolsESGetTokens_, endcapCorSuperClusterCollectionToken_, endcapPreSuperClusterCollectionToken_, endcapRawSuperClusterCollectionToken_, endcapRecHitCollectionToken_, dqm::impl::MonitorElement::Fill(), genWeightsTable_cfi::genEvent, EcalClusterLazyToolsBase::ESGetTokens::get(), 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(), fileCollector::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, ecalClusterToolsESGetTokens_.get(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);
     }
   }
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 105 of file EgammaSuperClusters.cc.

References HistSpec::bins, dqm::implementation::IBooker::book1D(), dqm::implementation::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::implementation::NavigatorBase::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);
 }

◆ closestMCParticle()

void EgammaSuperClusters::closestMCParticle	(	HepMC::GenEvent const *	genEvent,
		reco::SuperCluster const &	sc,
		double &	dRClosest,
		double &	energyClosest
	)		const

private

Definition at line 609 of file EgammaSuperClusters.cc.

References reco::deltaPhi(), PbPb_ZMuSkimMuonDPG_cff::deltaR, HLT_2022v12_cff::dEta, HLT_2022v12_cff::dPhi, ecalEta(), reco::CaloCluster::eta(), genWeightsTable_cfi::genEvent, reco::CaloCluster::phi(), mathSSE::sqrt(), and extraflags_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
 }

◆ ecalEta()

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(), Z_Endcap, and ZEcal.

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

◆ barrelCorSuperClusterCollectionToken_

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

private

Definition at line 41 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ barrelRawSuperClusterCollectionToken_

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

private

Definition at line 40 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ barrelRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::barrelRecHitCollectionToken_

private

Definition at line 49 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ ecalClusterToolsESGetTokens_

EcalClusterLazyTools::ESGetTokens EgammaSuperClusters::ecalClusterToolsESGetTokens_

private

Definition at line 52 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapCorSuperClusterCollectionToken_

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

private

Definition at line 46 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapPreSuperClusterCollectionToken_

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

private

Definition at line 45 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapRawSuperClusterCollectionToken_

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

private

Definition at line 44 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::endcapRecHitCollectionToken_

private

Definition at line 50 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ hist_EB_CorSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_deltaR_

private

Definition at line 110 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_EoverTruth_

private

Definition at line 105 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_

private

Definition at line 80 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_vs_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Eta_

private

Definition at line 119 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_vs_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Phi_

private

Definition at line 120 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Eta_

private

Definition at line 85 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_etaWidth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_etaWidth_

private

Definition at line 118 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_NumBC_

private

Definition at line 75 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Phi_

private

Definition at line 90 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_phiWidth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_phiWidth_

private

Definition at line 116 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S1toS9_

private

Definition at line 95 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S25toE_

private

Definition at line 100 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Size_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Size_

private

Definition at line 70 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_deltaR_

private

Definition at line 108 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_EoverTruth_

private

Definition at line 103 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_ET_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_ET_

private

Definition at line 78 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Eta_

private

Definition at line 83 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_NumBC_

private

Definition at line 73 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Phi_

private

Definition at line 88 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S1toS9_

private

Definition at line 93 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S25toE_

private

Definition at line 98 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Size_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Size_

private

Definition at line 68 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_deltaR_

private

Definition at line 111 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_EoverTruth_

private

Definition at line 106 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_

private

Definition at line 81 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Eta_

private

Definition at line 121 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Phi_

private

Definition at line 122 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_R_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_R_

private

Definition at line 123 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Eta_

private

Definition at line 86 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_etaWidth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_etaWidth_

private

Definition at line 117 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_NumBC_

private

Definition at line 76 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Phi_

private

Definition at line 91 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_phiWidth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_phiWidth_

private

Definition at line 115 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_preshowerE_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_preshowerE_

private

Definition at line 114 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S1toS9_

private

Definition at line 96 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S25toE_

private

Definition at line 101 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Size_

private

Definition at line 71 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_deltaR_

private

Definition at line 112 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_EoverTruth_

private

Definition at line 107 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_ET_

private

Definition at line 82 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Eta_

private

Definition at line 87 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_NumBC_

private

Definition at line 77 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Phi_

private

Definition at line 92 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_preshowerE_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_preshowerE_

private

Definition at line 113 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S1toS9_

private

Definition at line 97 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S25toE_

private

Definition at line 102 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Size_

private

Definition at line 72 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_deltaR_

private

Definition at line 109 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_EoverTruth_

private

Definition at line 104 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_ET_

private

Definition at line 79 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Eta_

private

Definition at line 84 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_NumBC_

private

Definition at line 74 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Phi_

private

Definition at line 89 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S1toS9_

private

Definition at line 94 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S25toE_

private

Definition at line 99 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Size_

private

Definition at line 69 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hsdeltaR_

HistSpec EgammaSuperClusters::hsdeltaR_

private

Definition at line 62 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsEoverTruth_

HistSpec EgammaSuperClusters::hsEoverTruth_

private

Definition at line 61 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsET_

HistSpec EgammaSuperClusters::hsET_

private

Definition at line 56 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsEta_

HistSpec EgammaSuperClusters::hsEta_

private

Definition at line 57 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsetaWidth_

HistSpec EgammaSuperClusters::hsetaWidth_

private

Definition at line 64 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsNumBC_

HistSpec EgammaSuperClusters::hsNumBC_

private

Definition at line 55 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsPhi_

HistSpec EgammaSuperClusters::hsPhi_

private

Definition at line 58 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsphiWidth_

HistSpec EgammaSuperClusters::hsphiWidth_

private

Definition at line 63 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hspreshowerE_

HistSpec EgammaSuperClusters::hspreshowerE_

private

Definition at line 65 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsR_

HistSpec EgammaSuperClusters::hsR_

private

Definition at line 66 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsS1toS9_

HistSpec EgammaSuperClusters::hsS1toS9_

private

Definition at line 59 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsS25toE_

HistSpec EgammaSuperClusters::hsS25toE_

private

Definition at line 60 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsSize_

HistSpec EgammaSuperClusters::hsSize_

private

Definition at line 54 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ MCTruthCollectionToken_

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

private

Definition at line 37 of file EgammaSuperClusters.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EgammaSuperClusters()

◆ ~EgammaSuperClusters()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ closestMCParticle()

◆ ecalEta()

Member Data Documentation

◆ barrelCorSuperClusterCollectionToken_

◆ barrelRawSuperClusterCollectionToken_

◆ barrelRecHitCollectionToken_

◆ ecalClusterToolsESGetTokens_

◆ endcapCorSuperClusterCollectionToken_

◆ endcapPreSuperClusterCollectionToken_

◆ endcapRawSuperClusterCollectionToken_

◆ endcapRecHitCollectionToken_

◆ 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_