#include <AlCaElectronsTest.h>

Inheritance diagram for AlCaElectronsTest:

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

virtual void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)

virtual void	beginJob ()

virtual void	endJob ()

	~AlCaElectronsTest ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	fillAroundBarrel (const EcalRecHitCollection *recHits, int eta, int phi)

void	fillAroundEndcap (const EcalRecHitCollection *recHits, int ics, int ips)

EcalRecHit	getMaximum (const EcalRecHitCollection *recHits)

Private Attributes
edm::InputTag	m_barrelAlCa

TH2F *	m_barrelGlobalCrystalsEnergy
	ECAL Energy. More...

TH2F *	m_barrelGlobalCrystalsEnergyMap
	ECAL EnergyMap. More...

TH2F *	m_barrelGlobalCrystalsMap
	ECAL map. More...

TH2F *	m_barrelLocalCrystalsEnergy
	local Energy More...

TH2F *	m_barrelLocalCrystalsMap
	local map More...

edm::InputTag	m_endcapAlCa

TH2F *	m_endcapGlobalCrystalsEnergy
	ECAL Energy. More...

TH2F *	m_endcapGlobalCrystalsEnergyMap
	ECAL EnergyMap. More...

TH2F *	m_endcapGlobalCrystalsMap
	ECAL map. More...

TH2F *	m_endcapLocalCrystalsEnergy
	local Energy More...

TH2F *	m_endcapLocalCrystalsMap
	local map More...

std::string	m_outputFileName

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

typedef WorkerT< EDAnalyzer >	WorkerType

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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

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 32 of file AlCaElectronsTest.h.

Constructor & Destructor Documentation

AlCaElectronsTest::AlCaElectronsTest ( const edm::ParameterSet & iConfig )

explicit

Definition at line 80 of file AlCaElectronsTest.h.

                                                                     :
   m_barrelAlCa (iConfig.getParameter<edm::InputTag> ("alcaBarrelHitCollection")) ,
   m_endcapAlCa (iConfig.getParameter<edm::InputTag> ("alcaEndcapHitCollection")) ,
   m_outputFileName (iConfig.getUntrackedParameter<std::string>
                       ("HistOutFile",std::string ("AlCaElectronsTest.root"))) 
 {}

AlCaElectronsTest::~AlCaElectronsTest ( )

inline

Definition at line 35 of file AlCaElectronsTest.h.

35 {}

Member Function Documentation

void AlCaElectronsTest::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 135 of file AlCaElectronsTest.h.

References dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, CaloRecHit::energy(), fillAroundBarrel(), fillAroundEndcap(), edm::Event::getByLabel(), getMaximum(), edm::EventBase::id(), EcalRecHit::id(), edm::HandleBase::isValid(), m_barrelAlCa, m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsMap, m_endcapAlCa, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsMap, edm::Handle< T >::product(), and hitfit::return.

 {
   //FIXME replace with msg logger
   std::cout << "[AlCaElectronsTest] analysing event " 
             << iEvent.id () << std::endl ;
 
   //PG get the collections  
   // get Barrel RecHits
   edm::Handle<EBRecHitCollection> barrelRecHitsHandle ;
   iEvent.getByLabel (m_barrelAlCa, barrelRecHitsHandle) ;
   if (!barrelRecHitsHandle.isValid()) {
       std::cerr << "[AlCaElectronsTest] caught std::exception "
                 << " in rertieving " << m_barrelAlCa 
                 << std::endl ;
       return ;
   } else {
       const EBRecHitCollection* barrelHitsCollection = barrelRecHitsHandle.product () ;
       //PG fill the histo with the maximum
       EcalRecHit barrelMax = getMaximum (barrelHitsCollection) ;
       EBDetId barrelMaxId (barrelMax.id ()) ; 
       m_barrelGlobalCrystalsMap->Fill (
           barrelMaxId.ieta () ,
           barrelMaxId.iphi () 
         ) ;
       m_barrelGlobalCrystalsEnergy->Fill (
           barrelMaxId.ieta () ,
           barrelMaxId.iphi () ,
           barrelMax.energy ()
         ) ;
       fillAroundBarrel (
           barrelHitsCollection, 
           barrelMaxId.ieta (), 
           barrelMaxId.iphi ()
         ) ;
   }
   
   // get Endcap RecHits
   edm::Handle<EERecHitCollection> endcapRecHitsHandle ;
   iEvent.getByLabel (m_endcapAlCa,endcapRecHitsHandle) ;
   if (!endcapRecHitsHandle.isValid()) {
       std::cerr << "[AlCaElectronsTest] caught std::exception " 
                 << " in rertieving " << m_endcapAlCa 
                 << std::endl ;
       return ;
   } else {
       const EERecHitCollection* endcapHitsCollection = endcapRecHitsHandle.product () ;
       //PG fill the histo with the maximum
       EcalRecHit endcapMax = getMaximum (endcapHitsCollection) ;
       EEDetId endcapMaxId (endcapMax.id ()) ; 
       m_endcapGlobalCrystalsMap->Fill (
           endcapMaxId.ix () ,
           endcapMaxId.iy () 
         ) ;
       m_endcapGlobalCrystalsEnergy->Fill (
           endcapMaxId.ix () ,
           endcapMaxId.iy () ,
           endcapMax.energy ()
         ) ;
       fillAroundEndcap (
           endcapHitsCollection, 
           endcapMaxId.ix (), 
           endcapMaxId.iy ()
         ) ;
   }
 }

void AlCaElectronsTest::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 92 of file AlCaElectronsTest.h.

References m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsEnergyMap, m_barrelGlobalCrystalsMap, m_barrelLocalCrystalsEnergy, m_barrelLocalCrystalsMap, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsEnergyMap, m_endcapGlobalCrystalsMap, m_endcapLocalCrystalsEnergy, m_endcapLocalCrystalsMap, and hitfit::return.

 {
   m_barrelGlobalCrystalsMap = new TH2F ("m_barrelGlobalCrystalsMap","m_barrelGlobalCrystalsMap",171,-85,86,360,0,360) ;
   m_barrelLocalCrystalsMap = new TH2F ("m_barrelLocalCrystalsMap","m_barrelLocalCrystalsMap",20,-10,10,20,-10,10) ;
   m_endcapGlobalCrystalsMap = new TH2F ("m_endcapGlobalCrystalsMap","m_endcapGlobalCrystalsMap",100,0,100,100,0,100) ;
   m_endcapLocalCrystalsMap = new TH2F ("m_endcapLocalCrystalsMap","m_endcapLocalCrystalsMap",20,-10,10,20,-10,10) ;
   m_barrelGlobalCrystalsEnergy = new TH2F ("m_barrelGlobalCrystalsEnergy","m_barrelGlobalCrystalsEnergy",171,-85,86,360,0,360) ;
   m_barrelLocalCrystalsEnergy = new TH2F ("m_barrelLocalCrystalsEnergy","m_barrelLocalCrystalsEnergy",20,-10,10,20,-10,10) ;
   m_endcapGlobalCrystalsEnergy = new TH2F ("m_endcapGlobalCrystalsEnergy","m_endcapGlobalCrystalsEnergy",100,0,100,100,0,100) ;
   m_endcapLocalCrystalsEnergy = new TH2F ("m_endcapLocalCrystalsEnergy","m_endcapLocalCrystalsEnergy",20,-10,10,20,-10,10) ;
   m_barrelGlobalCrystalsEnergyMap = new TH2F ("m_barrelGlobalCrystalsEnergyMap","m_barrelGlobalCrystalsEnergyMap",171,-85,86,360,0,360) ;
   m_endcapGlobalCrystalsEnergyMap = new TH2F ("m_endcapGlobalCrystalsEnergyMap","m_endcapGlobalCrystalsEnergyMap",100,0,100,100,0,100) ;
    return ;
 }

void AlCaElectronsTest::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 112 of file AlCaElectronsTest.h.

References m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsEnergyMap, m_barrelGlobalCrystalsMap, m_barrelLocalCrystalsEnergy, m_barrelLocalCrystalsMap, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsEnergyMap, m_endcapGlobalCrystalsMap, m_endcapLocalCrystalsEnergy, m_endcapLocalCrystalsMap, m_outputFileName, convertSQLitetoXML_cfg::output, and hitfit::return.

 {      
    TFile output (m_outputFileName.c_str (),"recreate") ;
    m_barrelGlobalCrystalsMap->Write () ;
    m_barrelLocalCrystalsMap->Write () ;
    m_endcapGlobalCrystalsMap->Write () ;
    m_endcapLocalCrystalsMap->Write () ;   
    m_barrelGlobalCrystalsEnergy->Write () ;
    m_barrelLocalCrystalsEnergy->Write () ;
    m_endcapGlobalCrystalsEnergy->Write () ;
    m_endcapLocalCrystalsEnergy->Write () ;   
    m_barrelGlobalCrystalsEnergyMap->Write () ;
    m_endcapGlobalCrystalsEnergyMap->Write () ;
    output.Close () ;
    //PG save root things
    return ;
 }

void AlCaElectronsTest::fillAroundBarrel	(	const EcalRecHitCollection *	recHits,
		int	eta,
		int	phi
	)

private

Definition at line 229 of file AlCaElectronsTest.h.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), eta(), EBDetId::ieta(), EBDetId::iphi(), m_barrelGlobalCrystalsEnergyMap, m_barrelLocalCrystalsEnergy, m_barrelLocalCrystalsMap, phi, and hitfit::return.

Referenced by analyze().

 {
   for (EcalRecHitCollection::const_iterator elem = recHits->begin () ;
        elem != recHits->end () ;
        ++elem)
     {
       EBDetId elementId = elem->id () ; 
       m_barrelLocalCrystalsMap->Fill (
         elementId.ieta () - eta ,
         elementId.iphi () - phi 
       ) ;
       m_barrelLocalCrystalsEnergy->Fill (
         elementId.ieta () - eta ,
         elementId.iphi () - phi ,
         elem->energy ()
       ) ;
      m_barrelGlobalCrystalsEnergyMap->Fill (
         elementId.ieta () ,
         elementId.iphi () ,
         elem->energy ()
       ) ;
 
     }   
   return ;
 }

void AlCaElectronsTest::fillAroundEndcap	(	const EcalRecHitCollection *	recHits,
		int	ics,
		int	ips
	)

private

Definition at line 260 of file AlCaElectronsTest.h.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), h2_mapping_cfi::ics, EEDetId::ix(), EEDetId::iy(), m_endcapGlobalCrystalsEnergyMap, m_endcapLocalCrystalsEnergy, m_endcapLocalCrystalsMap, and hitfit::return.

Referenced by analyze().

 {
   for (EcalRecHitCollection::const_iterator elem = recHits->begin () ;
        elem != recHits->end () ;
        ++elem)
     {
       EEDetId elementId = elem->id () ; 
       m_endcapLocalCrystalsMap->Fill (
         elementId.ix () - ics ,
         elementId.iy () - ips 
       ) ;
       m_endcapLocalCrystalsEnergy->Fill (
         elementId.ix () - ics ,
         elementId.iy () - ips ,
         elem->energy ()
       ) ;
       m_endcapGlobalCrystalsEnergyMap->Fill (
         elementId.ix () ,
         elementId.iy () ,
         elem->energy ()
       ) ;
     }   
   return ;
 }

EcalRecHit AlCaElectronsTest::getMaximum ( const EcalRecHitCollection * recHits )

private

Definition at line 207 of file AlCaElectronsTest.h.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, and max().

Referenced by analyze().

 {
   double energy = 0. ;
   EcalRecHit max ;
   for (EcalRecHitCollection::const_iterator elem = recHits->begin () ;
        elem != recHits->end () ;
        ++elem)
     {
       if (elem->energy () > energy)
         {
           energy = elem->energy () ;
           max = *elem ;
         }        
     }   
   return max ;
 }