#include <AlCaElectronsTest.h>

Inheritance diagram for AlCaElectronsTest:

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

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

void	beginJob () override

	~AlCaElectronsTest ()=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () 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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 35 of file AlCaElectronsTest.h.

Constructor & Destructor Documentation

◆ AlCaElectronsTest()

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

explicit

Definition at line 75 of file AlCaElectronsTest.h.

References TFileService::kSharedResource.

     : 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"))) {
   usesResource(TFileService::kSharedResource);
 }

◆ ~AlCaElectronsTest()

AlCaElectronsTest::~AlCaElectronsTest ( )

default

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 111 of file AlCaElectronsTest.h.

References EcalRecHit::energy(), fillAroundBarrel(), fillAroundEndcap(), getMaximum(), EcalRecHit::id(), iEvent, edm::HandleBase::isValid(), SummaryClient_cfi::labels, edm::EDConsumerBase::labelsForToken(), m_barrelAlCa, m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsMap, m_endcapAlCa, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsMap, and edm::Handle< T >::product().

                                                                                    {
   //FIXME replace with msg logger
   edm::LogVerbatim("ElectronsTest") << "[AlCaElectronsTest] analysing event " << iEvent.id();
 
   //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);
     edm::LogError("ElectronsTest") << "[AlCaElectronsTest] caught std::exception in rertieving " << labels.module;
     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);
     edm::LogError("ElectronsTest") << "[AlCaElectronsTest] caught std::exception in rertieving " << labels.module;
     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());
   }
 }

◆ beginJob()

void AlCaElectronsTest::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 85 of file AlCaElectronsTest.h.

References compareTotals::fs, m_barrelGlobalCrystalsEnergy, m_barrelGlobalCrystalsEnergyMap, m_barrelGlobalCrystalsMap, m_barrelLocalCrystalsEnergy, m_barrelLocalCrystalsMap, m_endcapGlobalCrystalsEnergy, m_endcapGlobalCrystalsEnergyMap, m_endcapGlobalCrystalsMap, m_endcapLocalCrystalsEnergy, and m_endcapLocalCrystalsMap.

                                  {
   edm::Service<TFileService> fs;
   m_barrelGlobalCrystalsMap =
       fs->make<TH2F>("m_barrelGlobalCrystalsMap", "m_barrelGlobalCrystalsMap", 171, -85, 86, 360, 0, 360);
   m_barrelLocalCrystalsMap =
       fs->make<TH2F>("m_barrelLocalCrystalsMap", "m_barrelLocalCrystalsMap", 20, -10, 10, 20, -10, 10);
   m_endcapGlobalCrystalsMap =
       fs->make<TH2F>("m_endcapGlobalCrystalsMap", "m_endcapGlobalCrystalsMap", 100, 0, 100, 100, 0, 100);
   m_endcapLocalCrystalsMap =
       fs->make<TH2F>("m_endcapLocalCrystalsMap", "m_endcapLocalCrystalsMap", 20, -10, 10, 20, -10, 10);
   m_barrelGlobalCrystalsEnergy =
       fs->make<TH2F>("m_barrelGlobalCrystalsEnergy", "m_barrelGlobalCrystalsEnergy", 171, -85, 86, 360, 0, 360);
   m_barrelLocalCrystalsEnergy =
       fs->make<TH2F>("m_barrelLocalCrystalsEnergy", "m_barrelLocalCrystalsEnergy", 20, -10, 10, 20, -10, 10);
   m_endcapGlobalCrystalsEnergy =
       fs->make<TH2F>("m_endcapGlobalCrystalsEnergy", "m_endcapGlobalCrystalsEnergy", 100, 0, 100, 100, 0, 100);
   m_endcapLocalCrystalsEnergy =
       fs->make<TH2F>("m_endcapLocalCrystalsEnergy", "m_endcapLocalCrystalsEnergy", 20, -10, 10, 20, -10, 10);
   m_barrelGlobalCrystalsEnergyMap =
       fs->make<TH2F>("m_barrelGlobalCrystalsEnergyMap", "m_barrelGlobalCrystalsEnergyMap", 171, -85, 86, 360, 0, 360);
   m_endcapGlobalCrystalsEnergyMap =
       fs->make<TH2F>("m_endcapGlobalCrystalsEnergyMap", "m_endcapGlobalCrystalsEnergyMap", 100, 0, 100, 100, 0, 100);
 }

◆ fillAroundBarrel()

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

private

Definition at line 169 of file AlCaElectronsTest.h.

References PVValHelper::eta, EBDetId::ieta(), EBDetId::iphi(), m_barrelGlobalCrystalsEnergyMap, m_barrelLocalCrystalsEnergy, m_barrelLocalCrystalsMap, phi, and FastTrackerRecHitMaskProducer_cfi::recHits.

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

◆ fillAroundEndcap()

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

private

Definition at line 181 of file AlCaElectronsTest.h.

References Ecal07UnpackerData_cfi::ics, EEDetId::ix(), EEDetId::iy(), m_endcapGlobalCrystalsEnergyMap, m_endcapLocalCrystalsEnergy, m_endcapLocalCrystalsMap, and FastTrackerRecHitMaskProducer_cfi::recHits.

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

◆ getMaximum()

EcalRecHit AlCaElectronsTest::getMaximum ( const EcalRecHitCollection * recHits )

private

Definition at line 155 of file AlCaElectronsTest.h.

References hcalRecHitTable_cff::energy, SiStripPI::max, and FastTrackerRecHitMaskProducer_cfi::recHits.

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

Member Data Documentation

◆ m_barrelAlCa

edm::EDGetTokenT<EBRecHitCollection> AlCaElectronsTest::m_barrelAlCa

private

Definition at line 47 of file AlCaElectronsTest.h.

Referenced by analyze().

◆ m_barrelGlobalCrystalsEnergy

TH2F* AlCaElectronsTest::m_barrelGlobalCrystalsEnergy

private

ECAL Energy.

Definition at line 60 of file AlCaElectronsTest.h.

Referenced by analyze(), and beginJob().

◆ m_barrelGlobalCrystalsEnergyMap

TH2F* AlCaElectronsTest::m_barrelGlobalCrystalsEnergyMap

private

ECAL EnergyMap.

Definition at line 68 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundBarrel().

◆ m_barrelGlobalCrystalsMap

TH2F* AlCaElectronsTest::m_barrelGlobalCrystalsMap

private

ECAL map.

Definition at line 52 of file AlCaElectronsTest.h.

Referenced by analyze(), and beginJob().

◆ m_barrelLocalCrystalsEnergy

TH2F* AlCaElectronsTest::m_barrelLocalCrystalsEnergy

private

local Energy

Definition at line 62 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundBarrel().

◆ m_barrelLocalCrystalsMap

TH2F* AlCaElectronsTest::m_barrelLocalCrystalsMap

private

local map

Definition at line 54 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundBarrel().

◆ m_endcapAlCa

edm::EDGetTokenT<EERecHitCollection> AlCaElectronsTest::m_endcapAlCa

private

Definition at line 48 of file AlCaElectronsTest.h.

Referenced by analyze().

◆ m_endcapGlobalCrystalsEnergy

TH2F* AlCaElectronsTest::m_endcapGlobalCrystalsEnergy

private

ECAL Energy.

Definition at line 64 of file AlCaElectronsTest.h.

Referenced by analyze(), and beginJob().

◆ m_endcapGlobalCrystalsEnergyMap

TH2F* AlCaElectronsTest::m_endcapGlobalCrystalsEnergyMap

private

ECAL EnergyMap.

Definition at line 70 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundEndcap().

◆ m_endcapGlobalCrystalsMap

TH2F* AlCaElectronsTest::m_endcapGlobalCrystalsMap

private

ECAL map.

Definition at line 56 of file AlCaElectronsTest.h.

Referenced by analyze(), and beginJob().

◆ m_endcapLocalCrystalsEnergy

TH2F* AlCaElectronsTest::m_endcapLocalCrystalsEnergy

private

local Energy

Definition at line 66 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundEndcap().

◆ m_endcapLocalCrystalsMap

TH2F* AlCaElectronsTest::m_endcapLocalCrystalsMap

private

local map

Definition at line 58 of file AlCaElectronsTest.h.

Referenced by beginJob(), and fillAroundEndcap().

◆ m_outputFileName

std::string AlCaElectronsTest::m_outputFileName

private

Definition at line 49 of file AlCaElectronsTest.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AlCaElectronsTest()

◆ ~AlCaElectronsTest()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ fillAroundBarrel()

◆ fillAroundEndcap()

◆ getMaximum()

Member Data Documentation

◆ m_barrelAlCa

◆ m_barrelGlobalCrystalsEnergy

◆ m_barrelGlobalCrystalsEnergyMap

◆ m_barrelGlobalCrystalsMap

◆ m_barrelLocalCrystalsEnergy

◆ m_barrelLocalCrystalsMap

◆ m_endcapAlCa

◆ m_endcapGlobalCrystalsEnergy

◆ m_endcapGlobalCrystalsEnergyMap

◆ m_endcapGlobalCrystalsMap

◆ m_endcapLocalCrystalsEnergy

◆ m_endcapLocalCrystalsMap

◆ m_outputFileName