16 std::cout<<
"I'm going to check the internal consistancy of EcalTPGScale transformation..."<<std::endl ;
38 std::vector<DetId>::const_iterator it ;
42 it = ebCells.begin() ;
46 for (
unsigned int ADC=0 ;
ADC<256 ;
ADC++) {
47 double gev = ecalScale.getTPGInGeV(
ADC, towidEB) ;
48 unsigned int tpgADC = ecalScale.getTPGInADC(gev, towidEB) ;
51 std::cout<<
" ERROR : with ADC = "<<
ADC<<
" getTPGInGeV = "<<gev<<
" getTPGInADC = "<<tpgADC<<std::endl ;
53 RCT_LUT_EB[
ADC] = ecalScale.getLinearizedTPG(
ADC, towidEB) ;
58 it = eeCells.begin() ;
62 for (
unsigned int ADC=0 ;
ADC<256 ;
ADC++) {
63 double gev = ecalScale.getTPGInGeV(
ADC, towidEE) ;
64 unsigned int tpgADC = ecalScale.getTPGInADC(gev, towidEE) ;
67 std::cout<<
" ERROR : with ADC = "<<
ADC<<
" getTPGInGeV = "<<gev<<
" getTPGInADC = "<<tpgADC<<std::endl ;
69 RCT_LUT_EE[
ADC] = ecalScale.getLinearizedTPG(
ADC, towidEE) ;
73 if (!error)
std::cout<<
" there is no error with EcalTPGScale internal consistancy "<<std::endl ;
void setEventSetup(const edm::EventSetup &evtSetup)
const CaloSubdetectorGeometry * theEndcapGeometry_
virtual const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const
Get a list of valid detector ids (for the given subdetector)
EcalTrigTowerDetId tower() const
get the HCAL/trigger iphi of this crystal
virtual void analyze(const edm::Event &evt, const edm::EventSetup &evtSetup)
testEcalTPGScale(edm::ParameterSet const &pSet)
const CaloSubdetectorGeometry * theBarrelGeometry_
edm::ESHandle< EcalTrigTowerConstituentsMap > eTTmap_