#include <ESRecoSummary.h>

Inheritance diagram for ESRecoSummary:

Public Member Functions
	ESRecoSummary (const edm::ParameterSet &)

	~ESRecoSummary () override

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

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

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

	~DQMEDAnalyzer () override=default

Protected Attributes
edm::EDGetTokenT< reco::PreshowerClusterCollection >	esClusterCollectionX_

edm::EDGetTokenT< reco::PreshowerClusterCollection >	esClusterCollectionY_

edm::EDGetTokenT< ESRecHitCollection >	esRecHitCollection_

edm::EDGetTokenT< reco::SuperClusterCollection >	superClusterCollection_EE_

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

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

Private Attributes
MonitorElement *	h_esClusters_energy_plane1

MonitorElement *	h_esClusters_energy_plane2

MonitorElement *	h_esClusters_energy_ratio

MonitorElement *	h_recHits_ES_energyMax

MonitorElement *	h_recHits_ES_time

std::string	prefixME_

Detailed Description

Definition at line 15 of file ESRecoSummary.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 42 of file ESRecoSummary.cc.

References esClusterCollectionX_, esClusterCollectionY_, esRecHitCollection_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), prefixME_, AlCaHLTBitMon_QueryRunRegistry::string, and superClusterCollection_EE_.

                                                       {
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   // now do what ever initialization is needed
   esRecHitCollection_ = consumes<ESRecHitCollection>(ps.getParameter<edm::InputTag>("recHitCollection_ES"));
   esClusterCollectionX_ =
       consumes<reco::PreshowerClusterCollection>(ps.getParameter<edm::InputTag>("ClusterCollectionX_ES"));
   esClusterCollectionY_ =
       consumes<reco::PreshowerClusterCollection>(ps.getParameter<edm::InputTag>("ClusterCollectionY_ES"));
 
   superClusterCollection_EE_ =
       consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("superClusterCollection_EE"));
 }

ESRecoSummary::~ESRecoSummary ( )

inlineoverride

Definition at line 18 of file ESRecoSummary.h.

References analyze(), and bookHistograms().

18 {}

Member Function Documentation

void ESRecoSummary::analyze	(	const edm::Event &	ev,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 74 of file ESRecoSummary.cc.

References edm::SortedCollection< T, SORT >::begin(), MessageLogger_cfi::cerr, DEFINE_FWK_MODULE, edm::SortedCollection< T, SORT >::end(), reco::CaloCluster::energy(), esClusterCollectionX_, esClusterCollectionY_, esRecHitCollection_, MonitorElement::Fill(), edm::Event::getByToken(), h_esClusters_energy_plane1, h_esClusters_energy_plane2, h_esClusters_energy_ratio, h_recHits_ES_energyMax, h_recHits_ES_time, edm::HandleBase::isValid(), edm::Handle< T >::product(), and superClusterCollection_EE_.

Referenced by ~ESRecoSummary().

                                                                      {
   // Preshower RecHits
   edm::Handle<ESRecHitCollection> recHitsES;
   ev.getByToken(esRecHitCollection_, recHitsES);
   const ESRecHitCollection *thePreShowerRecHits = recHitsES.product();
 
   if (!recHitsES.isValid()) {
     std::cerr << "ESRecoSummary::analyze --> recHitsES not found" << std::endl;
   }
 
   float maxRecHitEnergyES = -999.;
 
   for (ESRecHitCollection::const_iterator esItr = thePreShowerRecHits->begin(); esItr != thePreShowerRecHits->end();
        ++esItr) {
     h_recHits_ES_time->Fill(esItr->time());
     if (esItr->energy() > maxRecHitEnergyES)
       maxRecHitEnergyES = esItr->energy();
 
   }  // end loop over ES rec Hits
 
   h_recHits_ES_energyMax->Fill(maxRecHitEnergyES);
 
   // ES clusters in X plane
   edm::Handle<reco::PreshowerClusterCollection> esClustersX;
   ev.getByToken(esClusterCollectionX_, esClustersX);
   const reco::PreshowerClusterCollection *ESclustersX = esClustersX.product();
 
   // ES clusters in Y plane
   edm::Handle<reco::PreshowerClusterCollection> esClustersY;
   ev.getByToken(esClusterCollectionY_, esClustersY);
   const reco::PreshowerClusterCollection *ESclustersY = esClustersY.product();
 
   // ... endcap
   edm::Handle<reco::SuperClusterCollection> superClusters_EE_h;
   ev.getByToken(superClusterCollection_EE_, superClusters_EE_h);
   const reco::SuperClusterCollection *theEndcapSuperClusters = superClusters_EE_h.product();
   if (!superClusters_EE_h.isValid()) {
     std::cerr << "EcalRecHitSummary::analyze --> superClusters_EE_h not found" << std::endl;
   }
 
   // loop over all super clusters
   for (reco::SuperClusterCollection::const_iterator itSC = theEndcapSuperClusters->begin();
        itSC != theEndcapSuperClusters->end();
        ++itSC) {
     if (fabs(itSC->eta()) < 1.65 || fabs(itSC->eta()) > 2.6)
       continue;
 
     float ESenergyPlane1 = 0.;
     float ESenergyPlane2 = 0.;
 
     // Loop over all ECAL Basic clusters in the supercluster
     for (reco::CaloCluster_iterator ecalBasicCluster = itSC->clustersBegin(); ecalBasicCluster != itSC->clustersEnd();
          ecalBasicCluster++) {
       const reco::CaloClusterPtr ecalBasicClusterPtr = *(ecalBasicCluster);
 
       for (reco::PreshowerClusterCollection::const_iterator iESClus = ESclustersX->begin();
            iESClus != ESclustersX->end();
            ++iESClus) {
         const reco::CaloClusterPtr preshBasicCluster = iESClus->basicCluster();
         const reco::PreshowerCluster *esCluster = &*iESClus;
         if (preshBasicCluster == ecalBasicClusterPtr) {
           ESenergyPlane1 += esCluster->energy();
         }
       }  // end of x loop
 
       for (reco::PreshowerClusterCollection::const_iterator iESClus = ESclustersY->begin();
            iESClus != ESclustersY->end();
            ++iESClus) {
         const reco::CaloClusterPtr preshBasicCluster = iESClus->basicCluster();
         const reco::PreshowerCluster *esCluster = &*iESClus;
         if (preshBasicCluster == ecalBasicClusterPtr) {
           ESenergyPlane2 += esCluster->energy();
         }
       }  // end of y loop
     }    // end loop over all basic clusters in the supercluster
 
     // cout<<"DQM : "<<ESenergyPlane1<<" "<<ESenergyPlane2<<endl;
     h_esClusters_energy_plane1->Fill(ESenergyPlane1);
     h_esClusters_energy_plane2->Fill(ESenergyPlane2);
     if (ESenergyPlane1 > 0 && ESenergyPlane2 > 0)
       h_esClusters_energy_ratio->Fill(ESenergyPlane1 / ESenergyPlane2);
 
   }  // end loop over superclusters
 }

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

overrideprivate

Definition at line 56 of file ESRecoSummary.cc.

References DQMStore::IBooker::book1D(), h_esClusters_energy_plane1, h_esClusters_energy_plane2, h_esClusters_energy_ratio, h_recHits_ES_energyMax, h_recHits_ES_time, prefixME_, and DQMStore::IBooker::setCurrentFolder().

Referenced by ~ESRecoSummary().

                                                                                                   {
   // Monitor Elements (ex THXD)
   iBooker.setCurrentFolder(prefixME_ + "/ESRecoSummary");  // to organise the histos in folders
 
   // Preshower ----------------------------------------------
   h_recHits_ES_energyMax = iBooker.book1D("recHits_ES_energyMax", "recHits_ES_energyMax", 200, 0., 0.01);
   h_recHits_ES_time = iBooker.book1D("recHits_ES_time", "recHits_ES_time", 200, -100., 100.);
 
   h_esClusters_energy_plane1 = iBooker.book1D("esClusters_energy_plane1", "esClusters_energy_plane1", 200, 0., 0.01);
   h_esClusters_energy_plane2 = iBooker.book1D("esClusters_energy_plane2", "esClusters_energy_plane2", 200, 0., 0.01);
   h_esClusters_energy_ratio = iBooker.book1D("esClusters_energy_ratio", "esClusters_energy_ratio", 200, 0., 20.);
 }