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 ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 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)

Detailed Description

Description: SVSuite Super Cluster Validation

Implementation: \

Author

: Michael A. Balazs, Nov 2006

Definition at line 26 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_(0),
  hist_EE_RawSC_Size_(0),
  hist_EB_CorSC_Size_(0),
  hist_EE_CorSC_Size_(0),
  hist_EE_PreSC_Size_(0),
  hist_EB_RawSC_NumBC_(0),
  hist_EE_RawSC_NumBC_(0),
  hist_EB_CorSC_NumBC_(0),
  hist_EE_CorSC_NumBC_(0),
  hist_EE_PreSC_NumBC_(0),
  hist_EB_RawSC_ET_(0),
  hist_EE_RawSC_ET_(0),
  hist_EB_CorSC_ET_(0),
  hist_EE_CorSC_ET_(0),
  hist_EE_PreSC_ET_(0),
  hist_EB_RawSC_Eta_(0),
  hist_EE_RawSC_Eta_(0),
  hist_EB_CorSC_Eta_(0),
  hist_EE_CorSC_Eta_(0),
  hist_EE_PreSC_Eta_(0),
  hist_EB_RawSC_Phi_(0),
  hist_EE_RawSC_Phi_(0),
  hist_EB_CorSC_Phi_(0),
  hist_EE_CorSC_Phi_(0),
  hist_EE_PreSC_Phi_(0),
  hist_EB_RawSC_S1toS9_(0),
  hist_EE_RawSC_S1toS9_(0),
  hist_EB_CorSC_S1toS9_(0),
  hist_EE_CorSC_S1toS9_(0),
  hist_EE_PreSC_S1toS9_(0),
  hist_EB_RawSC_S25toE_(0),
  hist_EE_RawSC_S25toE_(0),
  hist_EB_CorSC_S25toE_(0),
  hist_EE_CorSC_S25toE_(0),
  hist_EE_PreSC_S25toE_(0),
  hist_EB_RawSC_EoverTruth_(0),
  hist_EE_RawSC_EoverTruth_(0),
  hist_EB_CorSC_EoverTruth_(0),
  hist_EE_CorSC_EoverTruth_(0),
  hist_EE_PreSC_EoverTruth_(0),
  hist_EB_RawSC_deltaR_(0),
  hist_EE_RawSC_deltaR_(0),
  hist_EB_CorSC_deltaR_(0),
  hist_EE_CorSC_deltaR_(0),
  hist_EE_PreSC_deltaR_(0),
  hist_EE_PreSC_preshowerE_(0),
  hist_EE_CorSC_preshowerE_(0),
  hist_EE_CorSC_phiWidth_(0),
  hist_EB_CorSC_phiWidth_(0),
  hist_EE_CorSC_etaWidth_(0),
  hist_EB_CorSC_etaWidth_(0),
  hist_EB_CorSC_ET_vs_Eta_(0),
  hist_EB_CorSC_ET_vs_Phi_(0),
  hist_EE_CorSC_ET_vs_Eta_(0),
  hist_EE_CorSC_ET_vs_Phi_(0),
  hist_EE_CorSC_ET_vs_R_(0)
{
}

EgammaSuperClusters::~EgammaSuperClusters

(

)

Definition at line 99 of file EgammaSuperClusters.cc.

{
}

Member Function Documentation

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

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 325 of file EgammaSuperClusters.cc.

References barrelCorSuperClusterCollectionToken_, barrelRawSuperClusterCollectionToken_, barrelRecHitCollectionToken_, closestMCParticle(), endcapCorSuperClusterCollectionToken_, endcapPreSuperClusterCollectionToken_, endcapRawSuperClusterCollectionToken_, endcapRecHitCollectionToken_, MonitorElement::Fill(), MCTruth::genEvent, edm::Event::getByToken(), 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, 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 104 of file EgammaSuperClusters.cc.

References HistSpec::bins, DQMStore::IBooker::book1D(), 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 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 *  genEvent, reco::SuperCluster const &  sc, double &  dRClosest, double &  energyClosest  ) const

private

Definition at line 578 of file EgammaSuperClusters.cc.

References reco::deltaPhi(), deltaR(), dPhi(), ecalEta(), reco::CaloCluster::eta(), reco::CaloCluster::phi(), and mathSSE::sqrt().

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 622 of file EgammaSuperClusters.cc.

References ETA, etaBarrelEndcap, cmsBatch::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 41 of file EgammaSuperClusters.h.

Referenced by analyze().

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

private

Definition at line 40 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::barrelRecHitCollectionToken_

private

Definition at line 49 of file EgammaSuperClusters.h.

Referenced by analyze().

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

private

Definition at line 46 of file EgammaSuperClusters.h.

Referenced by analyze().

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

private

Definition at line 45 of file EgammaSuperClusters.h.

Referenced by analyze().

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

private

Definition at line 44 of file EgammaSuperClusters.h.

Referenced by analyze().

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::endcapRecHitCollectionToken_

private

Definition at line 50 of file EgammaSuperClusters.h.

Referenced by analyze().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_deltaR_

private

Definition at line 108 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_EoverTruth_

private

Definition at line 103 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_

private

Definition at line 78 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Eta_

private

Definition at line 117 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Phi_

private

Definition at line 118 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Eta_

private

Definition at line 83 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_etaWidth_

private

Definition at line 116 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_NumBC_

private

Definition at line 73 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Phi_

private

Definition at line 88 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_phiWidth_

private

Definition at line 114 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S1toS9_

private

Definition at line 93 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S25toE_

private

Definition at line 98 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Size_

private

Definition at line 68 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_deltaR_

private

Definition at line 106 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_EoverTruth_

private

Definition at line 101 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_ET_

private

Definition at line 76 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Eta_

private

Definition at line 81 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_NumBC_

private

Definition at line 71 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Phi_

private

Definition at line 86 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S1toS9_

private

Definition at line 91 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S25toE_

private

Definition at line 96 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Size_

private

Definition at line 66 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_deltaR_

private

Definition at line 109 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_EoverTruth_

private

Definition at line 104 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_

private

Definition at line 79 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Eta_

private

Definition at line 119 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Phi_

private

Definition at line 120 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_R_

private

Definition at line 121 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Eta_

private

Definition at line 84 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_etaWidth_

private

Definition at line 115 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_NumBC_

private

Definition at line 74 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Phi_

private

Definition at line 89 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_phiWidth_

private

Definition at line 113 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_preshowerE_

private

Definition at line 112 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S1toS9_

private

Definition at line 94 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S25toE_

private

Definition at line 99 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Size_

private

Definition at line 69 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_deltaR_

private

Definition at line 110 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_EoverTruth_

private

Definition at line 105 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_ET_

private

Definition at line 80 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Eta_

private

Definition at line 85 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_NumBC_

private

Definition at line 75 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Phi_

private

Definition at line 90 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_preshowerE_

private

Definition at line 111 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S1toS9_

private

Definition at line 95 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S25toE_

private

Definition at line 100 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Size_

private

Definition at line 70 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_deltaR_

private

Definition at line 107 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_EoverTruth_

private

Definition at line 102 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_ET_

private

Definition at line 77 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Eta_

private

Definition at line 82 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_NumBC_

private

Definition at line 72 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Phi_

private

Definition at line 87 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S1toS9_

private

Definition at line 92 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S25toE_

private

Definition at line 97 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Size_

private

Definition at line 67 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

HistSpec EgammaSuperClusters::hsdeltaR_

private

Definition at line 60 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsEoverTruth_

private

Definition at line 59 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsET_

private

Definition at line 54 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsEta_

private

Definition at line 55 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsetaWidth_

private

Definition at line 62 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsNumBC_

private

Definition at line 53 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsPhi_

private

Definition at line 56 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsphiWidth_

private

Definition at line 61 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hspreshowerE_

private

Definition at line 63 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsR_

private

Definition at line 64 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsS1toS9_

private

Definition at line 57 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsS25toE_

private

Definition at line 58 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

HistSpec EgammaSuperClusters::hsSize_

private

Definition at line 52 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

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

private

Definition at line 37 of file EgammaSuperClusters.h.

Referenced by analyze().