EcalTPGScale Class Reference

#include <EcalTPGScale.h>

Public Member Functions
	EcalTPGScale ()
unsigned int	getLinearizedTPG (unsigned int ADC, const EcalTrigTowerDetId &towerId)
unsigned int	getTPGInADC (double energy, const EcalTrigTowerDetId &towerId)
double	getTPGInGeV (const EcalTriggerPrimitiveDigi &tpDigi)
double	getTPGInGeV (unsigned int ADC, const EcalTrigTowerDetId &towerId)
void	setEventSetup (const edm::EventSetup &evtSetup)
	~EcalTPGScale ()

Private Attributes
const edm::EventSetup *	setup_

Detailed Description

Definition at line 8 of file EcalTPGScale.h.

Constructor & Destructor Documentation

EcalTPGScale::EcalTPGScale

(

)

Definition at line 13 of file EcalTPGScale.cc.

{ }

EcalTPGScale::~EcalTPGScale

(

)

Definition at line 16 of file EcalTPGScale.cc.

{ }

Member Function Documentation

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

Definition at line 54 of file EcalTPGScale.cc.

References edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), EcalTPGGroups::getMap(), mps_fire::i, lumiPlot::lut, edm::ESHandle< T >::product(), DetId::rawId(), and setup_.

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

{
  int tpg10bits = 0 ;
 
  // Get compressed look-up table
  edm::ESHandle<EcalTPGLutGroup> lutGrpHandle;
  setup_->get<EcalTPGLutGroupRcd>().get( lutGrpHandle );
  const EcalTPGGroups::EcalTPGGroupsMap & lutGrpMap = lutGrpHandle.product()->getMap() ;  
  EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId.rawId()) ;
  uint32_t lutGrp = 999 ;
  if (itgrp != lutGrpMap.end()) lutGrp = itgrp->second ;

  edm::ESHandle<EcalTPGLutIdMap> lutMapHandle;
  setup_->get<EcalTPGLutIdMapRcd>().get( lutMapHandle );
  const EcalTPGLutIdMap::EcalTPGLutMap & lutMap = lutMapHandle.product()->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 ;
}

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

Definition at line 82 of file EcalTPGScale.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, EcalTPGPhysicsConst::Item::EtSat, edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), EcalTPGGroups::getMap(), EcalTPGPhysicsConst::getMap(), createfilelist::int, lumiPlot::lut, edm::ESHandle< T >::product(), DetId::rawId(), setup_, and EcalTrigTowerDetId::subDet().

Referenced by testEcalTPGScale::analyze().

{
  unsigned int tpgADC = 0 ;

  // 1. get lsb
  edm::ESHandle<EcalTPGPhysicsConst> physHandle;
  setup_->get<EcalTPGPhysicsConstRcd>().get( physHandle );
  const EcalTPGPhysicsConstMap & physMap = physHandle.product()->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
  edm::ESHandle<EcalTPGLutGroup> lutGrpHandle;
  setup_->get<EcalTPGLutGroupRcd>().get( lutGrpHandle );
  const EcalTPGGroups::EcalTPGGroupsMap & lutGrpMap = lutGrpHandle.product()->getMap() ;  
  EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId) ;
  uint32_t lutGrp = 0 ;
  if (itgrp != lutGrpMap.end()) lutGrp = itgrp->second ;
  
  edm::ESHandle<EcalTPGLutIdMap> lutMapHandle;
  setup_->get<EcalTPGLutIdMapRcd>().get( lutMapHandle );
  const EcalTPGLutIdMap::EcalTPGLutMap & lutMap = lutMapHandle.product()->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 ;
}

double EcalTPGScale::getTPGInGeV

(

const EcalTriggerPrimitiveDigi &

tpDigi

)

Definition at line 24 of file EcalTPGScale.cc.

References EcalTriggerPrimitiveDigi::compressedEt(), EcalTriggerPrimitiveDigi::id(), and DataFormats::L1THGCal::towerId.

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

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

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

Definition at line 31 of file EcalTPGScale.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, EcalTPGPhysicsConst::Item::EtSat, edm::EventSetup::get(), getLinearizedTPG(), EcalTPGPhysicsConst::getMap(), edm::ESHandle< T >::product(), DetId::rawId(), setup_, and EcalTrigTowerDetId::subDet().

{ 
  // 1. get lsb
  edm::ESHandle<EcalTPGPhysicsConst> physHandle;
  setup_->get<EcalTPGPhysicsConstRcd>().get( physHandle );
  const EcalTPGPhysicsConstMap & physMap = physHandle.product()->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) ;

}

void EcalTPGScale::setEventSetup

(

const edm::EventSetup &

evtSetup

)

Definition at line 19 of file EcalTPGScale.cc.

References setup_.

Referenced by testEcalTPGScale::analyze(), EcalEBTrigPrimAnalyzer::analyze(), L1CaloInputScalesGenerator::analyze(), L1CaloInputScaleTester::analyze(), EcalTrigPrimAnalyzer::analyze(), L1RCTSaveInput::analyze(), L1RCTLutWriter::analyze(), EcalDeadCellTriggerPrimitiveFilter::filter(), and EcalRecHitWorkerRecover::set().

{
  setup_ = &evtSetup ;
}

Member Data Documentation

const edm::EventSetup* EcalTPGScale::setup_

private

Definition at line 23 of file EcalTPGScale.h.

Referenced by getLinearizedTPG(), getTPGInADC(), getTPGInGeV(), and setEventSetup().