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::EcalTPGScale

(

)

Definition at line 12 of file EcalTPGScale.cc.

{}

~EcalTPGScale()

EcalTPGScale::~EcalTPGScale

(

)

Definition at line 14 of file EcalTPGScale.cc.

{}

Member Function Documentation

getLinearizedTPG()

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

Definition at line 47 of file EcalTPGScale.cc.

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

References edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), get, EcalTPGGroups::getMap(), mps_fire::i, edm::ESHandle< T >::product(), setup_, and ecaldqm::towerId().

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

getTPGInADC()

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

Definition at line 75 of file EcalTPGScale.cc.

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

References DetId::Ecal, EcalBarrel, EcalEndcap, HCALHighEnergyHPDFilter_cfi::energy, edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), get, EcalTPGGroups::getMap(), EcalTPGPhysicsConst::getMap(), createfilelist::int, B2GTnPMonitor_cfi::item, edm::ESHandle< T >::product(), DetId::rawId(), setup_, and ecaldqm::towerId().

Referenced by testEcalTPGScale::analyze().

getTPGInGeV() [1/2]

double EcalTPGScale::getTPGInGeV

(

const EcalTriggerPrimitiveDigi &

tpDigi

)

Definition at line 18 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
	)

Definition at line 24 of file EcalTPGScale.cc.

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

References DetId::Ecal, EcalBarrel, EcalEndcap, edm::EventSetup::get(), get, getLinearizedTPG(), EcalTPGPhysicsConst::getMap(), B2GTnPMonitor_cfi::item, edm::ESHandle< T >::product(), DetId::rawId(), setup_, and ecaldqm::towerId().

setEventSetup()

void EcalTPGScale::setEventSetup

(

const edm::EventSetup &

evtSetup

)

Definition at line 16 of file EcalTPGScale.cc.

{ setup_ = &evtSetup; }

References setup_.

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

Member Data Documentation

setup_

const edm::EventSetup* EcalTPGScale::setup_

private

Definition at line 22 of file EcalTPGScale.h.

Referenced by getLinearizedTPG(), getTPGInADC(), getTPGInGeV(), setEventSetup(), and upgradeWorkflowComponents.UpgradeWorkflow::setup().