#include <EcalTPGScale.h>

Classes
struct	Tokens

Public Member Functions
	EcalTPGScale (Tokens const &, const edm::EventSetup &evtSetup)

unsigned int	getLinearizedTPG (unsigned int ADC, const EcalTrigTowerDetId &towerId) const

unsigned int	getTPGInADC (double energy, const EcalTrigTowerDetId &towerId) const

double	getTPGInGeV (const EcalTriggerPrimitiveDigi &tpDigi) const

double	getTPGInGeV (unsigned int ADC, const EcalTrigTowerDetId &towerId) const

Private Attributes
EcalTPGLutIdMap const &	lut_

EcalTPGLutGroup const &	lutGrp_

EcalTPGPhysicsConst const &	phys_

Detailed Description

Definition at line 16 of file EcalTPGScale.h.

Constructor & Destructor Documentation

◆ EcalTPGScale()

EcalTPGScale::EcalTPGScale	(	Tokens const &	tokens,
		const edm::EventSetup &	evtSetup
	)

Definition at line 10 of file EcalTPGScale.cc.

     : phys_(evtSetup.getData(tokens.physConstToken_)),
       lutGrp_(evtSetup.getData(tokens.lutGrpToken_)),
       lut_(evtSetup.getData(tokens.lutMapToken_))
  
 {}

Member Function Documentation

◆ getLinearizedTPG()

unsigned int EcalTPGScale::getLinearizedTPG	(	unsigned int	ADC,
		const EcalTrigTowerDetId &	towerId
	)		const

Definition at line 44 of file EcalTPGScale.cc.

                                                                                                      {
   int tpg10bits = 0;
  
   const EcalTPGGroups::EcalTPGGroupsMap& lutGrpMap = lutGrp_.getMap();
   EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId.rawId());
   uint32_t lutGrp = 999;
   if (itgrp != lutGrpMap.end())
     lutGrp = itgrp->second;
  
   const EcalTPGLutIdMap::EcalTPGLutMap& lutMap = lut_.getMap();
   EcalTPGLutIdMap::EcalTPGLutMapItr itLut = lutMap.find(lutGrp);
   if (itLut != lutMap.end()) {
     const unsigned int* lut = (itLut->second).getLut();
     for (unsigned int i = 0; i < 1024; i++)
       if (ADC == (0xff & lut[i])) {
         tpg10bits = i;
         break;
       }
   }
  
   return tpg10bits;
 }

References EcalTPGLutIdMap::getMap(), EcalTPGGroups::getMap(), mps_fire::i, lut_, lutGrp_, and ecaldqm::towerId().

Referenced by testEcalTPGScale::analyze(), and getTPGInGeV().

◆ getTPGInADC()

unsigned int EcalTPGScale::getTPGInADC	(	double	energy,
		const EcalTrigTowerDetId &	towerId
	)		const

Definition at line 67 of file EcalTPGScale.cc.

                                                                                              {
   unsigned int tpgADC = 0;
  
   // 1. get lsb
  
   const EcalTPGPhysicsConstMap& physMap = phys_.getMap();
  
   uint32_t eb = DetId(DetId::Ecal, EcalBarrel).rawId();
   uint32_t ee = DetId(DetId::Ecal, EcalEndcap).rawId();
   EcalTPGPhysicsConstMapIterator it = physMap.end();
   if (towerId.subDet() == EcalBarrel)
     it = physMap.find(eb);
   else if (towerId.subDet() == EcalEndcap)
     it = physMap.find(ee);
   double lsb10bits = 0.;
   if (it != physMap.end()) {
     EcalTPGPhysicsConst::Item item = it->second;
     lsb10bits = item.EtSat / 1024.;
   }
  
   // 2. get compressed look-up table
  
   const EcalTPGGroups::EcalTPGGroupsMap& lutGrpMap = lutGrp_.getMap();
   EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId);
   uint32_t lutGrp = 0;
   if (itgrp != lutGrpMap.end())
     lutGrp = itgrp->second;
  
   const EcalTPGLutIdMap::EcalTPGLutMap& lutMap = lut_.getMap();
   EcalTPGLutIdMap::EcalTPGLutMapItr itLut = lutMap.find(lutGrp);
   if (itLut != lutMap.end()) {
     const unsigned int* lut = (itLut->second).getLut();
     if (lsb10bits > 0) {
       int tpgADC10b = int(energy / lsb10bits + 0.5);
       if (tpgADC10b >= 0 && tpgADC10b < 1024)
         tpgADC = (0xff & lut[tpgADC10b]);
       if (tpgADC10b >= 1024)
         tpgADC = 0xff;
     }
   }
  
   return tpgADC;
 }

References DetId::Ecal, EcalBarrel, EcalEndcap, HCALHighEnergyHPDFilter_cfi::energy, EcalTPGLutIdMap::getMap(), EcalTPGGroups::getMap(), EcalTPGPhysicsConst::getMap(), createfilelist::int, B2GTnPMonitor_cfi::item, lut_, lutGrp_, phys_, DetId::rawId(), and ecaldqm::towerId().

Referenced by testEcalTPGScale::analyze().

◆ getTPGInGeV() [1/2]

double EcalTPGScale::getTPGInGeV ( const EcalTriggerPrimitiveDigi & tpDigi ) const

Definition at line 17 of file EcalTPGScale.cc.

                                                                              {
   const EcalTrigTowerDetId& towerId = tpDigi.id();
   int ADC = tpDigi.compressedEt();
   return getTPGInGeV(ADC, towerId);
 }

References EcalTriggerPrimitiveDigi::compressedEt(), EcalTriggerPrimitiveDigi::id(), and ecaldqm::towerId().

Referenced by testEcalTPGScale::analyze(), L1CaloInputScalesGenerator::analyze(), L1CaloInputScaleTester::analyze(), EcalTrigPrimAnalyzer::analyze(), L1RCTSaveInput::analyze(), L1RCTLutWriter::analyze(), EcalRecHitWorkerRecover::estimateEnergy(), EcalRecHitWorkerRecover::run(), and EcalDeadCellTriggerPrimitiveFilter::setEvtTPstatus().

◆ getTPGInGeV() [2/2]

double EcalTPGScale::getTPGInGeV	(	unsigned int	ADC,
		const EcalTrigTowerDetId &	towerId
	)		const

Definition at line 23 of file EcalTPGScale.cc.

                                                                                           {
   // 1. get lsb
  
   const EcalTPGPhysicsConstMap& physMap = phys_.getMap();
   uint32_t eb = DetId(DetId::Ecal, EcalBarrel).rawId();
   uint32_t ee = DetId(DetId::Ecal, EcalEndcap).rawId();
   EcalTPGPhysicsConstMapIterator it = physMap.end();
   if (towerId.subDet() == EcalBarrel)
     it = physMap.find(eb);
   else if (towerId.subDet() == EcalEndcap)
     it = physMap.find(ee);
   double lsb10bits = 0.;
   if (it != physMap.end()) {
     EcalTPGPhysicsConst::Item item = it->second;
     lsb10bits = item.EtSat / 1024.;
   }
  
   // 2. linearized TPG
   return lsb10bits * getLinearizedTPG(ADC, towerId);
 }