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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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::EDGetTokenT< EBRecHitCollection >	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::EDGetTokenT< EERecHitCollection >	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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Constructor & Destructor Documentation

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

explicit

Definition at line 81 of file AlCaElectronsTest.h.

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

AlCaElectronsTest::~AlCaElectronsTest ( )

inline

Definition at line 36 of file AlCaElectronsTest.h.

References analyze(), beginJob(), endJob(), PVValHelper::eta, fillAroundBarrel(), fillAroundEndcap(), getMaximum(), iEvent, and phi.

36 {}

Member Function Documentation

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

virtual

Definition at line 136 of file AlCaElectronsTest.h.

References MessageLogger_cfi::cerr, gather_cfg::cout, EcalRecHit::energy(), fillAroundBarrel(), fillAroundEndcap(), edm::Event::getByToken(), getMaximum(), edm::EventBase::id(), EcalRecHit::id(), edm::HandleBase::isValid(), tablePrinter::labels, edm::EDConsumerBase::labelsForToken(), m_barrelAlCa, m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsMap, m_endcapAlCa, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsMap, edm::ProductLabels::module, edm::Handle< T >::product(), and mathSSE::return().

Referenced by ~AlCaElectronsTest().

 {
   //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.getByToken (m_barrelAlCa, barrelRecHitsHandle) ;
   if (!barrelRecHitsHandle.isValid()) {
       edm::EDConsumerBase::Labels labels;
       labelsForToken(m_barrelAlCa, labels);
       std::cerr << "[AlCaElectronsTest] caught std::exception "
                 << " in rertieving " << labels.module
                 << 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.getByToken (m_endcapAlCa,endcapRecHitsHandle) ;
   if (!endcapRecHitsHandle.isValid()) {
       edm::EDConsumerBase::Labels labels;
       labelsForToken(m_endcapAlCa, labels);
       std::cerr << "[AlCaElectronsTest] caught std::exception " 
                 << " in rertieving " << labels.module
                 << 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 93 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 mathSSE::return().

Referenced by ~AlCaElectronsTest().

 {
   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 113 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 mathSSE::return().

Referenced by o2olib.O2ORunMgr::executeJob(), and ~AlCaElectronsTest().

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

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

Referenced by analyze(), and ~AlCaElectronsTest().

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

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

Referenced by analyze(), and ~AlCaElectronsTest().

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

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

Referenced by analyze(), and ~AlCaElectronsTest().

 {
   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 ;
 }