#include <EERecoSummary.h>

Inheritance diagram for EERecoSummary:

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

	~EERecoSummary ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Attributes
edm::EDGetTokenT< edm::View < reco::CaloCluster > >	basicClusterCollection_EE_

double	ethrEE_

edm::EDGetTokenT < EcalRecHitCollection >	recHitCollection_EE_

edm::EDGetTokenT < EcalRecHitCollection >	redRecHitCollection_EE_

double	scEtThrEE_

edm::EDGetTokenT < reco::SuperClusterCollection >	superClusterCollection_EE_

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

virtual void	beginJob ()

virtual void	endJob ()

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

Definition at line 40 of file EERecoSummary.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file EERecoSummary.cc.

References basicClusterCollection_EE_, dqmStore_, ethrEE_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), h_basicClusters_recHits_EE_recoFlag, h_recHits_EE_recoFlag, h_recHits_EEM_Chi2, h_recHits_EEM_energyMax, h_recHits_EEM_time, h_recHits_EEP_Chi2, h_recHits_EEP_energyMax, h_recHits_EEP_time, h_redRecHits_EE_recoFlag, h_superClusters_EE_phi, h_superClusters_EEM_nBC, h_superClusters_EEP_nBC, h_superClusters_eta, cppFunctionSkipper::operator, prefixME_, recHitCollection_EE_, redRecHitCollection_EE_, scEtThrEE_, AlCaHLTBitMon_QueryRunRegistry::string, and superClusterCollection_EE_.

 {
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   //now do what ever initialization is needed
   recHitCollection_EE_       = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("recHitCollection_EE"));
   redRecHitCollection_EE_    = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("redRecHitCollection_EE"));
   basicClusterCollection_EE_ = consumes<edm::View<reco::CaloCluster> >(ps.getParameter<edm::InputTag>("basicClusterCollection_EE"));
   superClusterCollection_EE_ = consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("superClusterCollection_EE"));
 
   ethrEE_                    = ps.getParameter<double>("ethrEE");
 
   scEtThrEE_                 = ps.getParameter<double>("scEtThrEE");
 
   // DQM Store -------------------
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   // Monitor Elements (ex THXD)
   dqmStore_->setCurrentFolder(prefixME_ + "/EERecoSummary"); // to organise the histos in folders
      
   // ReducedRecHits ----------------------------------------------
   // ... endcap 
   h_redRecHits_EE_recoFlag = dqmStore_->book1D("redRecHits_EE_recoFlag","redRecHits_EE_recoFlag",16,-0.5,15.5);  
 
   // RecHits ---------------------------------------------- 
   // ... endcap
   h_recHits_EE_recoFlag = dqmStore_->book1D("recHits_EE_recoFlag","recHits_EE_recoFlag",16,-0.5,15.5);  
 
   // ... endcap +
   h_recHits_EEP_energyMax     = dqmStore_->book1D("recHits_EEP_energyMax","recHitsEEP_energyMax",110,-10,100);
   h_recHits_EEP_Chi2          = dqmStore_->book1D("recHits_EEP_Chi2","recHits_EEP_Chi2",200,0,100);
   h_recHits_EEP_time          = dqmStore_->book1D("recHits_EEP_time","recHits_EEP_time",200,-50,50);
 
   // ... endcap -
   h_recHits_EEM_energyMax     = dqmStore_->book1D("recHits_EEM_energyMax","recHits_EEM_energyMax",110,-10,100);
   h_recHits_EEM_Chi2          = dqmStore_->book1D("recHits_EEM_Chi2","recHits_EEM_Chi2",200,0,100);
   h_recHits_EEM_time          = dqmStore_->book1D("recHits_EEM_time","recHits_EEM_time",200,-50,50);
 
   // Basic Clusters ----------------------------------------------    
   // ... associated endcap rec hits
   h_basicClusters_recHits_EE_recoFlag = dqmStore_->book1D("basicClusters_recHits_EE_recoFlag","basicClusters_recHits_EE_recoFlag",16,-0.5,15.5);  
 
   // Super Clusters ----------------------------------------------
   // ... endcap
   h_superClusters_EEP_nBC    = dqmStore_->book1D("superClusters_EEP_nBC","superClusters_EEP_nBC",100,0.,100.);
   h_superClusters_EEM_nBC    = dqmStore_->book1D("superClusters_EEM_nBC","superClusters_EEM_nBC",100,0.,100.);
 
   h_superClusters_eta        = dqmStore_->book1D("superClusters_eta","superClusters_eta",150,-3.,3.);
   h_superClusters_EE_phi     = dqmStore_->book1D("superClusters_EE_phi","superClusters_EE_phi",360,-3.1415927,3.1415927);
   
 }

EERecoSummary::~EERecoSummary ( )

Definition at line 119 of file EERecoSummary.cc.

 {
         // do anything here that needs to be done at desctruction time
         // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 131 of file EERecoSummary.cc.

References basicClusterCollection_EE_, edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, eta(), ethrEE_, create_public_lumi_plots::exp, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::SortedCollection< T, SORT >::find(), edm::EventSetup::get(), edm::Event::getByToken(), h_basicClusters_recHits_EE_recoFlag, h_recHits_EE_recoFlag, h_recHits_EEM_Chi2, h_recHits_EEM_energyMax, h_recHits_EEM_time, h_recHits_EEP_Chi2, h_recHits_EEP_energyMax, h_recHits_EEP_time, h_redRecHits_EE_recoFlag, h_superClusters_EE_phi, h_superClusters_EEM_nBC, h_superClusters_EEP_nBC, h_superClusters_eta, edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), phi, edm::Handle< T >::product(), recHitCollection_EE_, redRecHitCollection_EE_, scEtThrEE_, funct::sin(), superClusterCollection_EE_, cond::rpcobgas::time, detailsBasic3DVector::z, and EEDetId::zside().

 {
   
   //Get the magnetic field
   edm::ESHandle<MagneticField> theMagField;
   iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
 
   // --- REDUCED REC HITS ------------------------------------------------------------------------------------- 
   // ... endcap
   edm::Handle<EcalRecHitCollection> redRecHitsEE;
   ev.getByToken( redRecHitCollection_EE_, redRecHitsEE );
   const EcalRecHitCollection* theEndcapEcalredRecHits = redRecHitsEE.product () ;
   if ( ! redRecHitsEE.isValid() ) {
     edm::LogWarning("EERecoSummary") << "redRecHitsEE not found"; 
   }
   
   for ( EcalRecHitCollection::const_iterator itr = theEndcapEcalredRecHits->begin () ;
         itr != theEndcapEcalredRecHits->end () ; ++itr)
   {
       
     EEDetId eeid( itr -> id() );
 
     h_redRecHits_EE_recoFlag->Fill( itr -> recoFlag() );
 
   }
 
   // --- REC HITS ------------------------------------------------------------------------------------- 
   
   // ... endcap
   edm::Handle<EcalRecHitCollection> recHitsEE;
   ev.getByToken( recHitCollection_EE_, recHitsEE );
   const EcalRecHitCollection* theEndcapEcalRecHits = recHitsEE.product () ;
   if ( ! recHitsEE.isValid() ) {
     edm::LogWarning("EERecoSummary") << "recHitsEE not found"; 
   }
 
   float maxRecHitEnergyEEP = -999.;
   float maxRecHitEnergyEEM = -999.;
 
   EEDetId eeid_MrecHitEEM;
   EEDetId eeid_MrecHitEEP;
 
   for ( EcalRecHitCollection::const_iterator itr = theEndcapEcalRecHits->begin () ;
     itr != theEndcapEcalRecHits->end () ; ++itr)
     {
       
       EEDetId eeid( itr -> id() );
 
       // EE+
       if ( eeid.zside() > 0 ){
 
         h_recHits_EE_recoFlag       -> Fill( itr -> recoFlag() );
 
     // max E rec hit
     if (itr -> energy() > maxRecHitEnergyEEP && 
         !(eeid.ix()>=41 && eeid.ix()<=60 && eeid.iy()>=41 && eeid.iy()<=60) ) {
       maxRecHitEnergyEEP = itr -> energy() ;
     }
     
     // only channels above noise
     if (  itr -> energy() > ethrEE_ ){
       h_recHits_EEP_Chi2          -> Fill( itr -> chi2() );
           h_recHits_EEP_time          -> Fill( itr -> time() );
     }
       }
       
       // EE-
       if ( eeid.zside() < 0 ){
     
         h_recHits_EE_recoFlag       -> Fill( itr -> recoFlag() );
     
     // max E rec hit
     if (itr -> energy() > maxRecHitEnergyEEM && 
         !(eeid.ix()>=41 && eeid.ix()<=60 && eeid.iy()>=41 && eeid.iy()<=60) ) {
       maxRecHitEnergyEEM = itr -> energy() ;
     }
     
     // only channels above noise
     if (  itr -> energy() > ethrEE_ ) {
       h_recHits_EEM_Chi2          -> Fill( itr -> chi2() );
           h_recHits_EEM_time          -> Fill( itr -> time() );
     }
 
       }
     } // end loop over EE rec hits
 
   // energy
   h_recHits_EEP_energyMax -> Fill( maxRecHitEnergyEEP );
   h_recHits_EEM_energyMax -> Fill( maxRecHitEnergyEEM );
 
   //--- BASIC CLUSTERS --------------------------------------------------------------
 
   // ... endcap
   edm::Handle<edm::View<reco::CaloCluster> > basicClusters_EE_h;
   if(ev.getByToken( basicClusterCollection_EE_, basicClusters_EE_h )){
 
     for (unsigned int icl = 0; icl < basicClusters_EE_h->size(); ++icl) {
     
       //Get the associated RecHits
       const std::vector<std::pair<DetId,float> > & hits= (*basicClusters_EE_h)[icl].hitsAndFractions();
       for (std::vector<std::pair<DetId,float> > ::const_iterator rh = hits.begin(); rh!=hits.end(); ++rh){
       
         EBRecHitCollection::const_iterator itrechit = theEndcapEcalRecHits->find((*rh).first);
         if (itrechit==theEndcapEcalRecHits->end()) continue;
         h_basicClusters_recHits_EE_recoFlag -> Fill ( itrechit -> recoFlag() );
       }
   
     }
   }
   else{
     //    edm::LogWarning("EERecoSummary") << "basicClusters_EE_h not found"; 
   }
 
   // Super Clusters
   // ... 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() ) {
     edm::LogWarning("EERecoSummary") << "superClusters_EE_h not found"; 
   }
 
   for (reco::SuperClusterCollection::const_iterator itSC = theEndcapSuperClusters->begin(); 
        itSC != theEndcapSuperClusters->end(); ++itSC ) {
 
     double scEt = itSC -> energy() * sin(2.*atan( exp(- itSC->position().eta() )));
 
     if (scEt < scEtThrEE_ ) continue;
 
     h_superClusters_eta       -> Fill( itSC -> eta() );
     h_superClusters_EE_phi    -> Fill( itSC -> phi() );
     
     if  ( itSC -> z() > 0 ){
       h_superClusters_EEP_nBC    -> Fill( itSC -> clustersSize() );      
     }
 
     if  ( itSC -> z() < 0 ){
       h_superClusters_EEM_nBC    -> Fill( itSC -> clustersSize() );      
     }
   }
 
   //--------------------------------------------------------
  
 }

void EERecoSummary::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 279 of file EERecoSummary.cc.

280 {

281 }

void EERecoSummary::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 285 of file EERecoSummary.cc.

286 {}

Member Data Documentation

edm::EDGetTokenT<edm::View<reco::CaloCluster> > EERecoSummary::basicClusterCollection_EE_

protected

Definition at line 92 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

DQMStore* EERecoSummary::dqmStore_

private

Definition at line 54 of file EERecoSummary.h.

Referenced by EERecoSummary().

double EERecoSummary::ethrEE_

protected

Definition at line 95 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_basicClusters_recHits_EE_recoFlag

private

Definition at line 77 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EE_recoFlag

private

Definition at line 66 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEM_Chi2

private

Definition at line 73 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEM_energyMax

private

Definition at line 72 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEM_time

private

Definition at line 74 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEP_Chi2

private

Definition at line 69 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEP_energyMax

private

Definition at line 68 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_recHits_EEP_time

private

Definition at line 70 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_redRecHits_EE_recoFlag

private

Definition at line 62 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_superClusters_EE_phi

private

Definition at line 85 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_superClusters_EEM_nBC

private

Definition at line 82 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_superClusters_EEP_nBC

private

Definition at line 81 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

MonitorElement* EERecoSummary::h_superClusters_eta

private

Definition at line 84 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

std::string EERecoSummary::prefixME_

private

Definition at line 56 of file EERecoSummary.h.

Referenced by EERecoSummary().

edm::EDGetTokenT<EcalRecHitCollection> EERecoSummary::recHitCollection_EE_

protected

Definition at line 90 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

edm::EDGetTokenT<EcalRecHitCollection> EERecoSummary::redRecHitCollection_EE_

protected

Definition at line 91 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

double EERecoSummary::scEtThrEE_

protected

Definition at line 97 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

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

protected

Definition at line 93 of file EERecoSummary.h.

Referenced by analyze(), and EERecoSummary().

Private Attributes
DQMStore *	dqmStore_

MonitorElement *	h_basicClusters_recHits_EE_recoFlag

MonitorElement *	h_recHits_EE_recoFlag

MonitorElement *	h_recHits_EEM_Chi2

MonitorElement *	h_recHits_EEM_energyMax

MonitorElement *	h_recHits_EEM_time

MonitorElement *	h_recHits_EEP_Chi2

MonitorElement *	h_recHits_EEP_energyMax

MonitorElement *	h_recHits_EEP_time

MonitorElement *	h_redRecHits_EE_recoFlag

MonitorElement *	h_superClusters_EE_phi

MonitorElement *	h_superClusters_EEM_nBC

MonitorElement *	h_superClusters_EEP_nBC

MonitorElement *	h_superClusters_eta

std::string	prefixME_

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation