#include <EgammaBasicClusters.h>

Inheritance diagram for EgammaBasicClusters:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	EgammaBasicClusters (const edm::ParameterSet &)

	~EgammaBasicClusters () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Attributes
edm::EDGetTokenT< reco::BasicClusterCollection >	barrelBasicClusterCollection_

edm::EDGetTokenT< reco::BasicClusterCollection >	endcapBasicClusterCollection_

MonitorElement *	hist_EB_BC_ET_

MonitorElement *	hist_EB_BC_ET_vs_Eta_

MonitorElement *	hist_EB_BC_ET_vs_Phi_

MonitorElement *	hist_EB_BC_Eta_

MonitorElement *	hist_EB_BC_NumRecHits_

MonitorElement *	hist_EB_BC_Phi_

MonitorElement *	hist_EB_BC_Size_

MonitorElement *	hist_EE_BC_ET_

MonitorElement *	hist_EE_BC_ET_vs_Eta_

MonitorElement *	hist_EE_BC_ET_vs_Phi_

MonitorElement *	hist_EE_BC_ET_vs_R_

MonitorElement *	hist_EE_BC_Eta_

MonitorElement *	hist_EE_BC_NumRecHits_

MonitorElement *	hist_EE_BC_Phi_

MonitorElement *	hist_EE_BC_Size_

HistSpec	hsET_

HistSpec	hsEta_

HistSpec	hsNumRecHits_

HistSpec	hsPhi_

HistSpec	hsR_

HistSpec	hsSize_

Detailed Description

Description: SVSuite Basic Cluster Validation

Implementation: \

Author: : Michael A. Balazs, Nov 2006

Definition at line 24 of file EgammaBasicClusters.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file EgammaBasicClusters.cc.

                                                                     :
   barrelBasicClusterCollection_(consumes<reco::BasicClusterCollection>(ps.getParameter<edm::InputTag>("barrelBasicClusterCollection"))),
   endcapBasicClusterCollection_(consumes<reco::BasicClusterCollection>(ps.getParameter<edm::InputTag>("endcapBasicClusterCollection"))),
   hsSize_(ps, "Size"),
   hsNumRecHits_(ps, "NumRecHits"),
   hsET_(ps, "ET"),
   hsEta_(ps, "Eta"),
   hsPhi_(ps, "Phi"),
   hsR_(ps, "R"),
   hist_EB_BC_Size_(nullptr),
   hist_EE_BC_Size_(nullptr),
   hist_EB_BC_NumRecHits_(nullptr),
   hist_EE_BC_NumRecHits_(nullptr),
   hist_EB_BC_ET_(nullptr),
   hist_EE_BC_ET_(nullptr),
   hist_EB_BC_Eta_(nullptr),
   hist_EE_BC_Eta_(nullptr),
   hist_EB_BC_Phi_(nullptr),
   hist_EE_BC_Phi_(nullptr),
   hist_EB_BC_ET_vs_Eta_(nullptr),
   hist_EB_BC_ET_vs_Phi_(nullptr),
   hist_EE_BC_ET_vs_Eta_(nullptr),
   hist_EE_BC_ET_vs_Phi_(nullptr),
   hist_EE_BC_ET_vs_R_(nullptr)
 {
 }

EgammaBasicClusters::~EgammaBasicClusters ( )

override

Definition at line 40 of file EgammaBasicClusters.cc.

41 {

42 }

Member Function Documentation

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

override

Definition at line 112 of file EgammaBasicClusters.cc.

References barrelBasicClusterCollection_, endcapBasicClusterCollection_, MonitorElement::Fill(), edm::Event::getByToken(), hist_EB_BC_ET_, hist_EB_BC_ET_vs_Eta_, hist_EB_BC_ET_vs_Phi_, hist_EB_BC_Eta_, hist_EB_BC_NumRecHits_, hist_EB_BC_Phi_, hist_EB_BC_Size_, hist_EE_BC_ET_, hist_EE_BC_ET_vs_Eta_, hist_EE_BC_ET_vs_Phi_, hist_EE_BC_ET_vs_R_, hist_EE_BC_Eta_, hist_EE_BC_NumRecHits_, hist_EE_BC_Phi_, hist_EE_BC_Size_, edm::HandleBase::isValid(), checklumidiff::l, funct::pow(), edm::Handle< T >::product(), and mathSSE::sqrt().

 {
   edm::Handle<reco::BasicClusterCollection> pBarrelBasicClusters;
   evt.getByToken(barrelBasicClusterCollection_, pBarrelBasicClusters);
   if (!pBarrelBasicClusters.isValid()) {
 
     Labels l;
     labelsForToken(barrelBasicClusterCollection_,l);
     edm::LogError("EgammaBasicClusters") << "Error! can't get collection with label " 
                                          << l.module;
   }
 
   const reco::BasicClusterCollection* barrelBasicClusters = pBarrelBasicClusters.product();
   hist_EB_BC_Size_->Fill(barrelBasicClusters->size());
 
   for(reco::BasicClusterCollection::const_iterator aClus = barrelBasicClusters->begin(); 
       aClus != barrelBasicClusters->end(); aClus++)
     {
       hist_EB_BC_NumRecHits_->Fill(aClus->size());
       hist_EB_BC_ET_->Fill(aClus->energy()/std::cosh(aClus->position().eta()));
       hist_EB_BC_Eta_->Fill(aClus->position().eta());
       hist_EB_BC_Phi_->Fill(aClus->position().phi());
         
       hist_EB_BC_ET_vs_Eta_->Fill( aClus->energy()/std::cosh(aClus->position().eta()), aClus->eta() );
       hist_EB_BC_ET_vs_Phi_->Fill( aClus->energy()/std::cosh(aClus->position().eta()), aClus->phi() );
     }
 
   edm::Handle<reco::BasicClusterCollection> pEndcapBasicClusters;
 
   evt.getByToken(endcapBasicClusterCollection_, pEndcapBasicClusters);
   if (!pEndcapBasicClusters.isValid()) {
 
     Labels l;
     labelsForToken(endcapBasicClusterCollection_,l);
     edm::LogError("EgammaBasicClusters") << "Error! can't get collection with label " 
                                          << l.module;
   }
 
   const reco::BasicClusterCollection* endcapBasicClusters = pEndcapBasicClusters.product();
   hist_EE_BC_Size_->Fill(endcapBasicClusters->size());
 
   for(reco::BasicClusterCollection::const_iterator aClus = endcapBasicClusters->begin(); 
       aClus != endcapBasicClusters->end(); aClus++)
     {
       hist_EE_BC_NumRecHits_->Fill(aClus->size());
       hist_EE_BC_ET_->Fill(aClus->energy()/std::cosh(aClus->position().eta()));
       hist_EE_BC_Eta_->Fill(aClus->position().eta());
       hist_EE_BC_Phi_->Fill(aClus->position().phi());
 
       hist_EE_BC_ET_vs_Eta_->Fill( aClus->energy()/std::cosh(aClus->position().eta()), aClus->eta() );
       hist_EE_BC_ET_vs_Phi_->Fill( aClus->energy()/std::cosh(aClus->position().eta()), aClus->phi() );
       hist_EE_BC_ET_vs_R_->Fill( aClus->energy()/std::cosh(aClus->position().eta()), 
                                  std::sqrt( std::pow(aClus->x(),2) + std::pow(aClus->y(),2) ) );
 
     }
 }

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

override

Definition at line 45 of file EgammaBasicClusters.cc.

References HistSpec::bins, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), hist_EB_BC_ET_, hist_EB_BC_ET_vs_Eta_, hist_EB_BC_ET_vs_Phi_, hist_EB_BC_Eta_, hist_EB_BC_NumRecHits_, hist_EB_BC_Phi_, hist_EB_BC_Size_, hist_EE_BC_ET_, hist_EE_BC_ET_vs_Eta_, hist_EE_BC_ET_vs_Phi_, hist_EE_BC_ET_vs_R_, hist_EE_BC_Eta_, hist_EE_BC_NumRecHits_, hist_EE_BC_Phi_, hist_EE_BC_Size_, hsET_, hsEta_, hsNumRecHits_, hsPhi_, hsR_, hsSize_, HistSpec::max, HistSpec::min, and DQMStore::IBooker::setCurrentFolder().

 {
   _ibooker.setCurrentFolder("EcalClusterV/EcalBasicClusters/");
 
   hist_EB_BC_Size_ 
     = _ibooker.book1D("hist_EB_BC_Size_","# Basic Clusters in Barrel",
                       hsSize_.bins, hsSize_.min, hsSize_.max);
   hist_EE_BC_Size_ 
     = _ibooker.book1D("hist_EE_BC_Size_","# Basic Clusters in Endcap",
                       hsSize_.bins, hsSize_.min, hsSize_.max);
 
   hist_EB_BC_NumRecHits_ 
     = _ibooker.book1D("hist_EB_BC_NumRecHits_","# of RecHits in Basic Clusters in Barrel",
                       hsNumRecHits_.bins, hsNumRecHits_.min, hsNumRecHits_.max);
   hist_EE_BC_NumRecHits_ 
     = _ibooker.book1D("hist_EE_BC_NumRecHits_","# of RecHits in Basic Clusters in Endcap",
                       hsNumRecHits_.bins, hsNumRecHits_.min, hsNumRecHits_.max);
 
   hist_EB_BC_ET_ 
     = _ibooker.book1D("hist_EB_BC_ET_","ET of Basic Clusters in Barrel",
                       hsET_.bins, hsET_.min, hsET_.max);
   hist_EE_BC_ET_ 
     = _ibooker.book1D("hist_EE_BC_ET_","ET of Basic Clusters in Endcap",
                       hsET_.bins, hsET_.min, hsET_.max);
 
   hist_EB_BC_Eta_ 
     = _ibooker.book1D("hist_EB_BC_Eta_","Eta of Basic Clusters in Barrel",
                       hsEta_.bins, hsEta_.min, hsEta_.max);
   hist_EE_BC_Eta_ 
     = _ibooker.book1D("hist_EE_BC_Eta_","Eta of Basic Clusters in Endcap",
                       hsEta_.bins, hsEta_.min, hsEta_.max);
 
   hist_EB_BC_Phi_
     = _ibooker.book1D("hist_EB_BC_Phi_","Phi of Basic Clusters in Barrel",
                       hsPhi_.bins, hsPhi_.min, hsPhi_.max);
   hist_EE_BC_Phi_ 
     = _ibooker.book1D("hist_EE_BC_Phi_","Phi of Basic Clusters in Endcap",
                       hsPhi_.bins, hsPhi_.min, hsPhi_.max);
 
     
   hist_EB_BC_ET_vs_Eta_
     = _ibooker.book2D( "hist_EB_BC_ET_vs_Eta_", "Basic Cluster ET versus Eta in Barrel", 
                        hsET_.bins, hsET_.min, hsET_.max,
                        hsEta_.bins, hsEta_.min, hsEta_.max );
 
   hist_EB_BC_ET_vs_Phi_
     = _ibooker.book2D( "hist_EB_BC_ET_vs_Phi_", "Basic Cluster ET versus Phi in Barrel", 
                        hsET_.bins, hsET_.min, hsET_.max,
                        hsPhi_.bins, hsPhi_.min, hsPhi_.max );
 
   hist_EE_BC_ET_vs_Eta_
     = _ibooker.book2D( "hist_EE_BC_ET_vs_Eta_", "Basic Cluster ET versus Eta in Endcap", 
                        hsET_.bins, hsET_.min, hsET_.max,
                        hsEta_.bins, hsEta_.min, hsEta_.max );
 
   hist_EE_BC_ET_vs_Phi_
     = _ibooker.book2D( "hist_EE_BC_ET_vs_Phi_", "Basic Cluster ET versus Phi in Endcap", 
                        hsET_.bins, hsET_.min, hsET_.max,
                        hsPhi_.bins, hsPhi_.min, hsPhi_.max );
 
   hist_EE_BC_ET_vs_R_
     = _ibooker.book2D( "hist_EE_BC_ET_vs_R_", "Basic Cluster ET versus Radius in Endcap", 
                        hsET_.bins, hsET_.min, hsET_.max,
                        hsR_.bins, hsR_.min, hsR_.max );
 }