#include <EBRecoSummary.h>

Inheritance diagram for EBRecoSummary:

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

	~EBRecoSummary ()

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_EB_

double	ethrEB_

edm::EDGetTokenT < EcalRecHitCollection >	recHitCollection_EB_

edm::EDGetTokenT < EcalRecHitCollection >	redRecHitCollection_EB_

double	scEtThrEB_

edm::EDGetTokenT < reco::SuperClusterCollection >	superClusterCollection_EB_

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

virtual void	beginJob ()

void	convxtalid (int &, int &)

int	diff_neta_s (int, int)

int	diff_nphi_s (int, int)

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 EBRecoSummary.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file EBRecoSummary.cc.

References basicClusterCollection_EB_, gather_cfg::cout, dqmStore_, ethrEB_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), h_basicClusters_recHits_EB_recoFlag, h_recHits_EB_Chi2, h_recHits_EB_E1oE4, h_recHits_EB_energyMax, h_recHits_EB_recoFlag, h_recHits_EB_time, h_redRecHits_EB_recoFlag, h_superClusters_EB_E1oE4, h_superClusters_EB_nBC, h_superClusters_EB_phi, h_superClusters_eta, cppFunctionSkipper::operator, prefixME_, recHitCollection_EB_, redRecHitCollection_EB_, scEtThrEB_, AlCaHLTBitMon_QueryRunRegistry::string, and superClusterCollection_EB_.

 {
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   //now do what ever initialization is needed
   recHitCollection_EB_       = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("recHitCollection_EB"));
   redRecHitCollection_EB_    = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("redRecHitCollection_EB"));
   basicClusterCollection_EB_ = consumes<edm::View<reco::CaloCluster> >(ps.getParameter<edm::InputTag>("basicClusterCollection_EB"));
   superClusterCollection_EB_ = consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("superClusterCollection_EB"));
 
   std::cout << "EBRecoSummary " << ps.getParameter<edm::InputTag>("basicClusterCollection_EB") << std::endl;
 
   ethrEB_                    = ps.getParameter<double>("ethrEB");
 
   scEtThrEB_                 = ps.getParameter<double>("scEtThrEB");
 
   // DQM Store -------------------
   dqmStore_= edm::Service<DQMStore>().operator->();
 
   // Monitor Elements (ex THXD)
   dqmStore_->setCurrentFolder(prefixME_ + "/EBRecoSummary"); // to organise the histos in folders
      
   // ReducedRecHits ----------------------------------------------
   // ... barrel 
   h_redRecHits_EB_recoFlag = dqmStore_->book1D("redRecHits_EB_recoFlag","redRecHits_EB_recoFlag",16,-0.5,15.5);  
 
   // RecHits ---------------------------------------------- 
   // ... barrel
   h_recHits_EB_energyMax     = dqmStore_->book1D("recHits_EB_energyMax","recHitsEB_energyMax",110,-10,100);
   h_recHits_EB_Chi2          = dqmStore_->book1D("recHits_EB_Chi2","recHits_EB_Chi2",200,0,100);
   h_recHits_EB_time          = dqmStore_->book1D("recHits_EB_time","recHits_EB_time",200,-50,50);
   h_recHits_EB_E1oE4         = dqmStore_->book1D("recHits_EB_E1oE4","recHitsEB_E1oE4",150, 0, 1.5);
   h_recHits_EB_recoFlag      = dqmStore_->book1D("recHits_EB_recoFlag","recHits_EB_recoFlag",16,-0.5,15.5);  
 
   // Basic Clusters ----------------------------------------------    
   // ... associated barrel rec hits
   h_basicClusters_recHits_EB_recoFlag = dqmStore_->book1D("basicClusters_recHits_EB_recoFlag","basicClusters_recHits_EB_recoFlag",16,-0.5,15.5);  
 
   // Super Clusters ----------------------------------------------
   // ... barrel
   h_superClusters_EB_nBC     = dqmStore_->book1D("superClusters_EB_nBC","superClusters_EB_nBC",100,0.,100.);
   h_superClusters_EB_E1oE4   = dqmStore_->book1D("superClusters_EB_E1oE4","superClusters_EB_E1oE4",150,0,1.5);
 
   h_superClusters_eta        = dqmStore_->book1D("superClusters_eta","superClusters_eta",150,-3.,3.);
   h_superClusters_EB_phi     = dqmStore_->book1D("superClusters_EB_phi","superClusters_EB_phi",360,-3.1415927,3.1415927);
   
 }

EBRecoSummary::~EBRecoSummary ( )

Definition at line 115 of file EBRecoSummary.cc.

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

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 127 of file EBRecoSummary.cc.

References funct::abs(), basicClusterCollection_EB_, edm::View< T >::begin(), edm::SortedCollection< T, SORT >::begin(), jptDQMConfig_cff::eMax, edm::View< T >::end(), edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, eta(), ethrEB_, create_public_lumi_plots::exp, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::SortedCollection< T, SORT >::find(), edm::EventSetup::get(), edm::Event::getByToken(), h_basicClusters_recHits_EB_recoFlag, h_recHits_EB_Chi2, h_recHits_EB_E1oE4, h_recHits_EB_energyMax, h_recHits_EB_recoFlag, h_recHits_EB_time, h_redRecHits_EB_recoFlag, h_superClusters_EB_E1oE4, h_superClusters_EB_nBC, h_superClusters_EB_phi, h_superClusters_eta, EBDetId::ieta(), edm::HandleBase::isValid(), phi, edm::Handle< T >::product(), edm::ESHandle< class >::product(), recHitCollection_EB_, redRecHitCollection_EB_, scEtThrEB_, funct::sin(), superClusterCollection_EB_, EcalTools::swissCross(), and cond::rpcobgas::time.

 {
   
   //Get the magnetic field
   edm::ESHandle<MagneticField> theMagField;
   iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
 
   // calo topology
   edm::ESHandle<CaloTopology> pTopology;
   iSetup.get<CaloTopologyRecord>().get(pTopology);
   const CaloTopology *topology = pTopology.product();
 
   // --- REDUCED REC HITS ------------------------------------------------------------------------------------- 
   edm::Handle<EcalRecHitCollection> redRecHitsEB;
   ev.getByToken( redRecHitCollection_EB_, redRecHitsEB );
   const EcalRecHitCollection* theBarrelEcalredRecHits = redRecHitsEB.product () ;
   if ( ! redRecHitsEB.isValid() ) {
     edm::LogWarning("EBRecoSummary") << "redRecHitsEB not found"; 
   }
   
   for ( EcalRecHitCollection::const_iterator itr = theBarrelEcalredRecHits->begin () ;
         itr != theBarrelEcalredRecHits->end () ;++itr)
   {
       
     h_redRecHits_EB_recoFlag->Fill( itr -> recoFlag() );
   
   }
 
   // --- REC HITS ------------------------------------------------------------------------------------- 
   
   // ... barrel
   edm::Handle<EcalRecHitCollection> recHitsEB;
   ev.getByToken( recHitCollection_EB_, recHitsEB );
   const EcalRecHitCollection* theBarrelEcalRecHits = recHitsEB.product () ;
   if ( ! recHitsEB.isValid() ) {
     edm::LogWarning("EBRecoSummary") << "recHitsEB not found"; 
   }
 
   float maxRecHitEnergyEB = -999.;
   
   EBDetId ebid_MrecHitEB;
 
   for ( EcalRecHitCollection::const_iterator itr = theBarrelEcalRecHits->begin () ;
     itr != theBarrelEcalRecHits->end () ;++itr)
     {
 
       EBDetId ebid( itr -> id() );
       
       h_recHits_EB_recoFlag      -> Fill( itr -> recoFlag() );
 
       // max E rec hit
       if (itr -> energy() > maxRecHitEnergyEB ){
     maxRecHitEnergyEB = itr -> energy() ;
       }       
 
       if ( itr -> energy() > ethrEB_ ){
     h_recHits_EB_Chi2          -> Fill( itr -> chi2() );
         h_recHits_EB_time          -> Fill( itr -> time() );
       }
 
       float R4 = EcalTools::swissCross( ebid, *theBarrelEcalRecHits, 0. );
       
       if ( itr -> energy() > 3. && abs(ebid.ieta())!=85 )  h_recHits_EB_E1oE4-> Fill( R4 );
       
     }
   
   h_recHits_EB_energyMax         -> Fill( maxRecHitEnergyEB );
   
   //--- BASIC CLUSTERS --------------------------------------------------------------
 
   // ... barrel
   edm::Handle<edm::View<reco::CaloCluster> > basicClusters_EB_h;
   if(ev.getByToken( basicClusterCollection_EB_, basicClusters_EB_h )){
 
     const edm::View<reco::CaloCluster>* theBarrelBasicClusters = basicClusters_EB_h.product () ;
 
     for (edm::View<reco::CaloCluster>::const_iterator itBC = theBarrelBasicClusters->begin(); 
          itBC != theBarrelBasicClusters->end(); ++itBC ) {
          
       //Get the associated RecHits
       const std::vector<std::pair<DetId,float> > & hits= itBC->hitsAndFractions();
       for (std::vector<std::pair<DetId,float> > ::const_iterator rh = hits.begin(); rh!=hits.end(); ++rh){
       
         EBRecHitCollection::const_iterator itrechit = theBarrelEcalRecHits->find((*rh).first);
         if (itrechit==theBarrelEcalRecHits->end()) continue;
         h_basicClusters_recHits_EB_recoFlag -> Fill ( itrechit -> recoFlag() );
     
       }
   
     }
   }
   else{
     //    edm::LogWarning("EBRecoSummary") << "basicClusters_EB_h not found"; 
   }
  
   // Super Clusters
   // ... barrel
   edm::Handle<reco::SuperClusterCollection> superClusters_EB_h;
   ev.getByToken( superClusterCollection_EB_, superClusters_EB_h );
   const reco::SuperClusterCollection* theBarrelSuperClusters = superClusters_EB_h.product () ;
   if ( ! superClusters_EB_h.isValid() ) {
     edm::LogWarning("EBRecoSummary") << "superClusters_EB_h not found"; 
   }
 
   for (reco::SuperClusterCollection::const_iterator itSC = theBarrelSuperClusters->begin(); 
        itSC != theBarrelSuperClusters->end(); ++itSC ) {
     
     double scEt = itSC -> energy() * sin(2.*atan( exp(- itSC->position().eta() )));
     
     if (scEt < scEtThrEB_ ) continue;
 
     h_superClusters_EB_nBC    -> Fill( itSC -> clustersSize());
     h_superClusters_eta       -> Fill( itSC -> eta() );
     h_superClusters_EB_phi    -> Fill( itSC -> phi() );
  
     float E1 = EcalClusterTools::eMax   ( *itSC, theBarrelEcalRecHits);
     float E4 = EcalClusterTools::eTop   ( *itSC, theBarrelEcalRecHits, topology )+
                EcalClusterTools::eRight ( *itSC, theBarrelEcalRecHits, topology )+
                EcalClusterTools::eBottom( *itSC, theBarrelEcalRecHits, topology )+
                EcalClusterTools::eLeft  ( *itSC, theBarrelEcalRecHits, topology );
 
     if ( E1 > 3. ) h_superClusters_EB_E1oE4  -> Fill( 1.- E4/E1);
     
   }
 
 }

void EBRecoSummary::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 257 of file EBRecoSummary.cc.

258 {

259 }

void EBRecoSummary::convxtalid	(	int &	,
		int &
	)

private

Definition at line 268 of file EBRecoSummary.cc.

 {
   // Barrel only
   // Output nphi 0...359; neta 0...84; nside=+1 (for eta>0), or 0 (for eta<0).
   // neta will be [-85,-1] , or [0,84], the minus sign indicates the z<0 side.
   
   if(neta > 0) neta -= 1;
   if(nphi > 359) nphi=nphi-360;
   
 } //end of convxtalid

int EBRecoSummary::diff_neta_s	(	int	,
		int
	)

private

Definition at line 279 of file EBRecoSummary.cc.

                                                       {
   Int_t mdiff;
   mdiff=(neta1-neta2);
   return mdiff;
 }

int EBRecoSummary::diff_nphi_s	(	int	,
		int
	)

private

Definition at line 286 of file EBRecoSummary.cc.

References funct::abs().

                                                       {
   Int_t mdiff;
   if(abs(nphi1-nphi2) < (360-abs(nphi1-nphi2))) {
     mdiff=nphi1-nphi2;
   }
   else {
     mdiff=360-abs(nphi1-nphi2);
     if(nphi1>nphi2) mdiff=-mdiff;
   }
   return mdiff;
 }

void EBRecoSummary::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 263 of file EBRecoSummary.cc.

264 {}

Member Data Documentation

edm::EDGetTokenT<edm::View<reco::CaloCluster> > EBRecoSummary::basicClusterCollection_EB_

protected

Definition at line 93 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

DQMStore* EBRecoSummary::dqmStore_

private

Definition at line 59 of file EBRecoSummary.h.

Referenced by EBRecoSummary().

double EBRecoSummary::ethrEB_

protected

Definition at line 96 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_basicClusters_recHits_EB_recoFlag

private

Definition at line 78 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_recHits_EB_Chi2

private

Definition at line 72 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_recHits_EB_E1oE4

private

Definition at line 74 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_recHits_EB_energyMax

private

Definition at line 71 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_recHits_EB_recoFlag

private

Definition at line 75 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_recHits_EB_time

private

Definition at line 73 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_redRecHits_EB_recoFlag

private

Definition at line 67 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_superClusters_EB_E1oE4

private

Definition at line 83 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_superClusters_EB_nBC

private

Definition at line 82 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_superClusters_EB_phi

private

Definition at line 86 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

MonitorElement* EBRecoSummary::h_superClusters_eta

private

Definition at line 85 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

std::string EBRecoSummary::prefixME_

private

Definition at line 61 of file EBRecoSummary.h.

Referenced by EBRecoSummary().

edm::EDGetTokenT<EcalRecHitCollection> EBRecoSummary::recHitCollection_EB_

protected

Definition at line 91 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

edm::EDGetTokenT<EcalRecHitCollection> EBRecoSummary::redRecHitCollection_EB_

protected

Definition at line 92 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

double EBRecoSummary::scEtThrEB_

protected

Definition at line 98 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

edm::EDGetTokenT<reco::SuperClusterCollection> EBRecoSummary::superClusterCollection_EB_

protected

Definition at line 94 of file EBRecoSummary.h.

Referenced by analyze(), and EBRecoSummary().

Private Attributes
DQMStore *	dqmStore_

MonitorElement *	h_basicClusters_recHits_EB_recoFlag

MonitorElement *	h_recHits_EB_Chi2

MonitorElement *	h_recHits_EB_E1oE4

MonitorElement *	h_recHits_EB_energyMax

MonitorElement *	h_recHits_EB_recoFlag

MonitorElement *	h_recHits_EB_time

MonitorElement *	h_redRecHits_EB_recoFlag

MonitorElement *	h_superClusters_EB_E1oE4

MonitorElement *	h_superClusters_EB_nBC

MonitorElement *	h_superClusters_EB_phi

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