#include <EcalTrigPrimAnalyzer.h>

Inheritance diagram for EcalTrigPrimAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	EcalTrigPrimAnalyzer (const edm::ParameterSet &)

	~EcalTrigPrimAnalyzer () override=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 Attributes
edm::ESGetToken< CaloSubdetectorGeometry, EcalBarrelGeometryRecord >	barrelGeomToken_

const edm::EDGetTokenT< EcalRecHitCollection >	ebToken_

TH1I *	ecal_et_ [2]

TH1I *	ecal_fgvb_ [2]

std::vector< std::string >	ecal_parts_

TH1I *	ecal_tt_ [2]

const edm::EDGetTokenT< EcalRecHitCollection >	eeToken_

edm::ESGetToken< CaloSubdetectorGeometry, EcalEndcapGeometryRecord >	endcapGeomToken_

float	eRec_

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	eTTmapToken_

int	fg_

TH1I *	histBar

TH1I *	histEndc

TH1F *	hTPoverRechit_

TH2F *	hTPvsRechit_

int	ieta_

int	iphi_

const edm::InputTag	label_

const bool	recHits_

const edm::InputTag	rechits_labelEB_

const edm::InputTag	rechits_labelEE_

EcalTPGScale::Tokens	tokens_

int	tpgADC_

float	tpgGeV_

const edm::EDGetTokenT< EcalTrigPrimDigiCollection >	tpToken_

TTree *	tree_

int	ttf_

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

Description: test of the output of EcalTrigPrimProducer

Description: rereads the result of the EcalTrigPrimProducer

Definition at line 50 of file EcalTrigPrimAnalyzer.h.

Constructor & Destructor Documentation

◆ EcalTrigPrimAnalyzer()

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

explicit

Definition at line 33 of file EcalTrigPrimAnalyzer.cc.

References barrelGeomToken_, ecal_et_, ecal_fgvb_, ecal_parts_, ecal_tt_, endcapGeomToken_, eRec_, eTTmapToken_, fg_, compareTotals::fs, hTPoverRechit_, hTPvsRechit_, mps_fire::i, ieta_, iphi_, TFileService::kSharedResource, recHits_, runGCPTkAlMap::title, tpgADC_, tpgGeV_, tree_, and ttf_.

     : recHits_(iConfig.getParameter<bool>("AnalyzeRecHits")),
       label_(iConfig.getParameter<edm::InputTag>("inputTP")),
       rechits_labelEB_(iConfig.getParameter<edm::InputTag>("inputRecHitsEB")),
       rechits_labelEE_(iConfig.getParameter<edm::InputTag>("inputRecHitsEE")),
       tpToken_(consumes<EcalTrigPrimDigiCollection>(label_)),
       ebToken_(consumes<EcalRecHitCollection>(rechits_labelEB_)),
       eeToken_(consumes<EcalRecHitCollection>(rechits_labelEE_)),
       tokens_(consumesCollector()) {
   usesResource(TFileService::kSharedResource);
 
   ecal_parts_.push_back("Barrel");
   ecal_parts_.push_back("Endcap");
 
   edm::Service<TFileService> fs;
   tree_ = fs->make<TTree>("TPGtree", "TPGtree");
   tree_->Branch("iphi", &iphi_, "iphi/I");
   tree_->Branch("ieta", &ieta_, "ieta/I");
   tree_->Branch("eRec", &eRec_, "eRec/F");
   tree_->Branch("tpgADC", &tpgADC_, "tpgADC/I");
   tree_->Branch("tpgGeV", &tpgGeV_, "tpgGeV/F");
   tree_->Branch("ttf", &ttf_, "ttf/I");
   tree_->Branch("fg", &fg_, "fg/I");
   for (unsigned int i = 0; i < 2; ++i) {
     ecal_et_[i] = fs->make<TH1I>(ecal_parts_[i].c_str(), "Et", 255, 0, 255);
     char title[30];
     sprintf(title, "%s_ttf", ecal_parts_[i].c_str());
     ecal_tt_[i] = fs->make<TH1I>(title, "TTF", 10, 0, 10);
     sprintf(title, "%s_fgvb", ecal_parts_[i].c_str());
     ecal_fgvb_[i] = fs->make<TH1I>(title, "FGVB", 10, 0, 10);
   }
 
   if (recHits_) {
     hTPvsRechit_ = fs->make<TH2F>("TP_vs_RecHit", "TP vs  rechit", 256, -1, 255, 255, 0, 255);
     hTPoverRechit_ = fs->make<TH1F>("TP_over_RecHit", "TP over rechit", 500, 0, 4);
     endcapGeomToken_ = esConsumes<CaloSubdetectorGeometry, EcalEndcapGeometryRecord>(edm::ESInputTag("", "EcalEndcap"));
     barrelGeomToken_ = esConsumes<CaloSubdetectorGeometry, EcalBarrelGeometryRecord>(edm::ESInputTag("", "EcalBarrel"));
     eTTmapToken_ = esConsumes<EcalTrigTowerConstituentsMap, IdealGeometryRecord>();
   }
 }

◆ ~EcalTrigPrimAnalyzer()

EcalTrigPrimAnalyzer::~EcalTrigPrimAnalyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 79 of file EcalTrigPrimAnalyzer.cc.

References barrelGeomToken_, submitPVResolutionJobs::count, ztail::d, ebToken_, ecal_et_, ecal_fgvb_, ecal_tt_, eeToken_, endcapGeomToken_, eRec_, eTTmapToken_, fg_, edm::EventSetup::getData(), EcalTPGScale::getTPGInGeV(), hTPoverRechit_, hTPvsRechit_, mps_fire::i, EcalTrigTowerDetId::ieta(), ieta_, iEvent, EcalTrigTowerDetId::iphi(), iphi_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, recHits_, funct::sin(), edm::SortedCollection< T, SORT >::size(), theta(), tokens_, EBDetId::tower(), cmsswSequenceInfo::tp, tpgADC_, tpgGeV_, tpToken_, tree_, and ttf_.

                                                                                       {
   // Get input
   const auto &tp = iEvent.get(tpToken_);
   for (unsigned int i = 0; i < tp.size(); i++) {
     const EcalTriggerPrimitiveDigi &d = tp[i];
     int subdet = d.id().subDet() - 1;
     if (subdet == 0) {
       ecal_et_[subdet]->Fill(d.compressedEt());
     } else {
       if (d.id().ietaAbs() == 27 || d.id().ietaAbs() == 28) {
         if (i % 2)
           ecal_et_[subdet]->Fill(d.compressedEt() * 2.);
       } else
         ecal_et_[subdet]->Fill(d.compressedEt());
     }
     ecal_tt_[subdet]->Fill(d.ttFlag());
     ecal_fgvb_[subdet]->Fill(d.fineGrain());
   }
   if (!recHits_)
     return;
 
   // comparison with RecHits
   const EcalRecHitCollection &rechit_EB_col = iEvent.get(ebToken_);
   const EcalRecHitCollection &rechit_EE_col = iEvent.get(eeToken_);
 
   const auto &theEndcapGeometry = iSetup.getData(endcapGeomToken_);
   const auto &theBarrelGeometry = iSetup.getData(barrelGeomToken_);
   const auto &eTTmap = iSetup.getData(eTTmapToken_);
 
   std::map<EcalTrigTowerDetId, float> mapTow_Et;
 
   for (unsigned int i = 0; i < rechit_EB_col.size(); i++) {
     const EBDetId &myid1 = rechit_EB_col[i].id();
     EcalTrigTowerDetId towid1 = myid1.tower();
     float theta = theBarrelGeometry.getGeometry(myid1)->getPosition().theta();
     float Etsum = rechit_EB_col[i].energy() * sin(theta);
     bool test_alreadyin = false;
     std::map<EcalTrigTowerDetId, float>::iterator ittest = mapTow_Et.find(towid1);
     if (ittest != mapTow_Et.end())
       test_alreadyin = true;
     if (test_alreadyin)
       continue;
     unsigned int j = i + 1;
     bool loopend = false;
     unsigned int count = 0;
     while (j < rechit_EB_col.size() && !loopend) {
       count++;
       const EBDetId &myid2 = rechit_EB_col[j].id();
       EcalTrigTowerDetId towid2 = myid2.tower();
       if (towid1 == towid2) {
         float theta = theBarrelGeometry.getGeometry(myid2)->getPosition().theta();
         Etsum += rechit_EB_col[j].energy() * sin(theta);
       }
       j++;
       if (count > 1800)
         loopend = true;
     }
     mapTow_Et.insert(std::pair<EcalTrigTowerDetId, float>(towid1, Etsum));
   }
 
   for (unsigned int i = 0; i < rechit_EE_col.size(); i++) {
     const EEDetId &myid1 = rechit_EE_col[i].id();
     EcalTrigTowerDetId towid1 = eTTmap.towerOf(myid1);
     float theta = theEndcapGeometry.getGeometry(myid1)->getPosition().theta();
     float Etsum = rechit_EE_col[i].energy() * sin(theta);
     bool test_alreadyin = false;
     std::map<EcalTrigTowerDetId, float>::iterator ittest = mapTow_Et.find(towid1);
     if (ittest != mapTow_Et.end())
       test_alreadyin = true;
     if (test_alreadyin)
       continue;
     unsigned int j = i + 1;
     bool loopend = false;
     unsigned int count = 0;
     while (j < rechit_EE_col.size() && !loopend) {
       const EEDetId &myid2 = rechit_EE_col[j].id();
       EcalTrigTowerDetId towid2 = eTTmap.towerOf(myid2);
       if (towid1 == towid2) {
         float theta = theEndcapGeometry.getGeometry(myid2)->getPosition().theta();
         Etsum += rechit_EE_col[j].energy() * sin(theta);
       }
       //  else loopend=true;
       j++;
       if (count > 500)
         loopend = true;
     }
     //    alreadyin_EE.push_back(towid1);
     mapTow_Et.insert(std::pair<EcalTrigTowerDetId, float>(towid1, Etsum));
   }
 
   EcalTPGScale ecalScale(tokens_, iSetup);
   for (unsigned int i = 0; i < tp.size(); i++) {
     const EcalTriggerPrimitiveDigi &d = tp[i];
     const EcalTrigTowerDetId TPtowid = d.id();
     std::map<EcalTrigTowerDetId, float>::iterator it = mapTow_Et.find(TPtowid);
     float Et = ecalScale.getTPGInGeV(d.compressedEt(), TPtowid);
     if (d.id().ietaAbs() == 27 || d.id().ietaAbs() == 28)
       Et *= 2;
     iphi_ = TPtowid.iphi();
     ieta_ = TPtowid.ieta();
     tpgADC_ = d.compressedEt();
     tpgGeV_ = Et;
     ttf_ = d.ttFlag();
     fg_ = d.fineGrain();
     if (it != mapTow_Et.end()) {
       hTPvsRechit_->Fill(it->second, Et);
       hTPoverRechit_->Fill(Et / it->second);
       eRec_ = it->second;
     }
     tree_->Fill();
   }
 }

Member Data Documentation

◆ barrelGeomToken_

edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> EcalTrigPrimAnalyzer::barrelGeomToken_

private

Definition at line 79 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ ebToken_

const edm::EDGetTokenT<EcalRecHitCollection> EcalTrigPrimAnalyzer::ebToken_

private

Definition at line 76 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze().

◆ ecal_et_

TH1I* EcalTrigPrimAnalyzer::ecal_et_[2]

private

Definition at line 60 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ ecal_fgvb_

TH1I* EcalTrigPrimAnalyzer::ecal_fgvb_[2]

private

Definition at line 62 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ ecal_parts_

std::vector<std::string> EcalTrigPrimAnalyzer::ecal_parts_

private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

Referenced by EcalTrigPrimAnalyzer().

◆ ecal_tt_

TH1I* EcalTrigPrimAnalyzer::ecal_tt_[2]

private

Definition at line 61 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ eeToken_

const edm::EDGetTokenT<EcalRecHitCollection> EcalTrigPrimAnalyzer::eeToken_

private

Definition at line 77 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze().

◆ endcapGeomToken_

edm::ESGetToken<CaloSubdetectorGeometry, EcalEndcapGeometryRecord> EcalTrigPrimAnalyzer::endcapGeomToken_

private

Definition at line 78 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ eRec_

float EcalTrigPrimAnalyzer::eRec_

private

Definition at line 69 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ eTTmapToken_

edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> EcalTrigPrimAnalyzer::eTTmapToken_

private

Definition at line 80 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ fg_

int EcalTrigPrimAnalyzer::fg_

private

Definition at line 68 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ histBar

TH1I * EcalTrigPrimAnalyzer::histBar

private

Definition at line 63 of file EcalTrigPrimAnalyzer.h.

◆ histEndc

TH1I* EcalTrigPrimAnalyzer::histEndc

private

Definition at line 63 of file EcalTrigPrimAnalyzer.h.

◆ hTPoverRechit_

TH1F* EcalTrigPrimAnalyzer::hTPoverRechit_

private

Definition at line 65 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ hTPvsRechit_

TH2F* EcalTrigPrimAnalyzer::hTPvsRechit_

private

Definition at line 64 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ ieta_

int EcalTrigPrimAnalyzer::ieta_

private

Definition at line 68 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ iphi_

int EcalTrigPrimAnalyzer::iphi_

private

Definition at line 68 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ label_

const edm::InputTag EcalTrigPrimAnalyzer::label_

private

Definition at line 72 of file EcalTrigPrimAnalyzer.h.

Referenced by Modules.SwitchProducer::__addParameter(), SequenceTypes.DummyModule::__repr__(), Modules.SwitchProducer::__setattr__(), Mixins._Labelable::_findDependencies(), and Mixins._Labelable::setLabel().

◆ recHits_

const bool EcalTrigPrimAnalyzer::recHits_

private

Definition at line 71 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ rechits_labelEB_

const edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEB_

private

Definition at line 73 of file EcalTrigPrimAnalyzer.h.

◆ rechits_labelEE_

const edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEE_

private

Definition at line 74 of file EcalTrigPrimAnalyzer.h.

◆ tokens_

EcalTPGScale::Tokens EcalTrigPrimAnalyzer::tokens_

private

Definition at line 81 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze().

◆ tpgADC_

int EcalTrigPrimAnalyzer::tpgADC_

private

Definition at line 68 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ tpgGeV_

float EcalTrigPrimAnalyzer::tpgGeV_

private

Definition at line 69 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ tpToken_

const edm::EDGetTokenT<EcalTrigPrimDigiCollection> EcalTrigPrimAnalyzer::tpToken_

private

Definition at line 75 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze().

◆ tree_

TTree* EcalTrigPrimAnalyzer::tree_

private

Definition at line 66 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

◆ ttf_

int EcalTrigPrimAnalyzer::ttf_

private

Definition at line 68 of file EcalTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalTrigPrimAnalyzer().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalTrigPrimAnalyzer()

◆ ~EcalTrigPrimAnalyzer()

Member Function Documentation

◆ analyze()

Member Data Documentation

◆ barrelGeomToken_

◆ ebToken_

◆ ecal_et_

◆ ecal_fgvb_

◆ ecal_parts_

◆ ecal_tt_

◆ eeToken_

◆ endcapGeomToken_

◆ eRec_

◆ eTTmapToken_

◆ fg_

◆ histBar

◆ histEndc

◆ hTPoverRechit_

◆ hTPvsRechit_

◆ ieta_

◆ iphi_

◆ label_

◆ recHits_

◆ rechits_labelEB_

◆ rechits_labelEE_

◆ tokens_

◆ tpgADC_

◆ tpgGeV_

◆ tpToken_

◆ tree_

◆ ttf_