|
|
|
#include <EcalTPGScale.h>
Public Member Functions | |
| EcalTPGScale () | |
| virtual double | et_bin_high (const EcalTrigTowerDetId &id, const EcalTriggerPrimitiveSample &s) const =0 |
| virtual double | et_bin_low (const EcalTrigTowerDetId &id, const EcalTriggerPrimitiveSample &s) const =0 |
| smallest ET value which would resolve into this sample | |
| virtual double | et_RCT (const EcalTrigTowerDetId &id, const EcalTriggerPrimitiveSample &s) const =0 |
| nominal ET value of this sample/code as to be used in RCT LUT creation | |
| 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) |
| virtual void | releaseSetup () const |
| Release any objects obtained from the EventSetup. | |
| void | setEventSetup (const edm::EventSetup &evtSetup) |
| virtual void | setup (const edm::EventSetup &es) const |
| Get any needed information from the event setup. | |
| ~EcalTPGScale () | |
Private Attributes | |
| const edm::EventSetup * | setup_ |
To be stored in the CaloTPGRecord
Definition at line 8 of file EcalTPGScale.h.
| EcalTPGScale::EcalTPGScale | ( | ) |
Definition at line 13 of file EcalTPGScale.cc.
{ }
| EcalTPGScale::~EcalTPGScale | ( | ) |
Definition at line 16 of file EcalTPGScale.cc.
{ }
| virtual double EcalTPGScale::et_bin_high | ( | const EcalTrigTowerDetId & | id, |
| const EcalTriggerPrimitiveSample & | s | ||
| ) | const [pure virtual] |
smallest ET value which would not resolve into sample (too big)
| virtual double EcalTPGScale::et_bin_low | ( | const EcalTrigTowerDetId & | id, |
| const EcalTriggerPrimitiveSample & | s | ||
| ) | const [pure virtual] |
smallest ET value which would resolve into this sample
| virtual double EcalTPGScale::et_RCT | ( | const EcalTrigTowerDetId & | id, |
| const EcalTriggerPrimitiveSample & | s | ||
| ) | const [pure virtual] |
nominal ET value of this sample/code as to be used in RCT LUT creation
| unsigned int EcalTPGScale::getLinearizedTPG | ( | unsigned int | ADC, |
| const EcalTrigTowerDetId & | towerId | ||
| ) |
Definition at line 54 of file EcalTPGScale.cc.
References edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), i, lumiPlot::lut, edm::ESHandle< T >::product(), DetId::rawId(), and setup_.
Referenced by 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::EventSetupRecord::get(), lumiPlot::lut, edm::ESHandle< T >::product(), DetId::rawId(), setup_, and EcalTrigTowerDetId::subDet().
{
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 ecaldqm::towerId().
Referenced by L1RCTLutWriter::analyze(), L1RCTSaveInput::analyze(), EcalSelectiveReadoutValidation::analyzeTP(), 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(), 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) ;
}
| virtual void EcalTPGScale::releaseSetup | ( | ) | const [inline, virtual] |
| void EcalTPGScale::setEventSetup | ( | const edm::EventSetup & | evtSetup | ) |
Definition at line 19 of file EcalTPGScale.cc.
References setup_.
Referenced by L1RCTLutWriter::analyze(), L1RCTSaveInput::analyze(), testEcalTPGScale::analyze(), EcalTrigPrimAnalyzer::analyze(), EcalSelectiveReadoutValidation::analyzeTP(), EcalDeadCellTriggerPrimitiveFilter::envSet(), EcalDeadCellDeltaRFilter::envSet(), and EcalRecHitWorkerRecover::set().
{
setup_ = &evtSetup ;
}
| virtual void EcalTPGScale::setup | ( | const edm::EventSetup & | es | ) | const [inline, virtual] |
const edm::EventSetup* EcalTPGScale::setup_ [private] |
Definition at line 23 of file EcalTPGScale.h.
Referenced by getLinearizedTPG(), getTPGInADC(), getTPGInGeV(), and setEventSetup().
1.7.3