ESRecoSummary Class Reference

#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 53 of file ESRecoSummary.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  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().

{}

Member Function Documentation

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

overrideprivate

Definition at line 86 of file ESRecoSummary.cc.

References edm::SortedCollection< T, SORT >::begin(), MessageLogger_cfi::cerr, DEFINE_FWK_MODULE, edm::SortedCollection< T, SORT >::end(), reco::CaloCluster::energy(), HcalObjRepresent::Fill(), edm::Event::getByToken(), edm::HandleBase::isValid(), and edm::Handle< T >::product().

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 66 of file ESRecoSummary.cc.

References DQMStore::IBooker::book1D(), 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.);
}

Member Data Documentation

edm::EDGetTokenT<reco::PreshowerClusterCollection> ESRecoSummary::esClusterCollectionX_

protected

Definition at line 40 of file ESRecoSummary.h.

edm::EDGetTokenT<reco::PreshowerClusterCollection> ESRecoSummary::esClusterCollectionY_

protected

Definition at line 41 of file ESRecoSummary.h.

edm::EDGetTokenT<ESRecHitCollection> ESRecoSummary::esRecHitCollection_

protected

Definition at line 39 of file ESRecoSummary.h.

MonitorElement* ESRecoSummary::h_esClusters_energy_plane1

private

Definition at line 30 of file ESRecoSummary.h.

MonitorElement* ESRecoSummary::h_esClusters_energy_plane2

private

Definition at line 31 of file ESRecoSummary.h.

MonitorElement* ESRecoSummary::h_esClusters_energy_ratio

private

Definition at line 32 of file ESRecoSummary.h.

MonitorElement* ESRecoSummary::h_recHits_ES_energyMax

private

Definition at line 27 of file ESRecoSummary.h.

MonitorElement* ESRecoSummary::h_recHits_ES_time

private

Definition at line 28 of file ESRecoSummary.h.

std::string ESRecoSummary::prefixME_

private

Definition at line 24 of file ESRecoSummary.h.

edm::EDGetTokenT<reco::SuperClusterCollection> ESRecoSummary::superClusterCollection_EE_

protected

Definition at line 38 of file ESRecoSummary.h.