|
|
|
#include <EcalRingCalibrationTools.h>
Public Member Functions | |
| EcalRingCalibrationTools () | |
| ~EcalRingCalibrationTools () | |
Static Public Member Functions | |
| static std::vector< DetId > | getDetIdsInECAL () |
| static std::vector< DetId > | getDetIdsInModule (short int) |
| static std::vector< DetId > | getDetIdsInRing (short aRingIndex) |
| Retrieve the DetIds in a phi-ring. | |
| static short | getModuleIndex (DetId aDetId) |
| static short | getRingIndex (DetId aDetId) |
| Retrieve the phi-ring index corresponding to a DetId. | |
| static void | setCaloGeometry (const CaloGeometry *geometry) |
Static Public Attributes | |
| static const short | N_MODULES_BARREL = 144 |
| static const short | N_RING_BARREL = 170 |
| static const short | N_RING_ENDCAP = 78 |
| static const short | N_RING_TOTAL = 248 |
Static Private Member Functions | |
| static void | initializeFromGeometry () |
Static Private Attributes | |
| static const CaloGeometry * | caloGeometry_ = 0 |
| static short | endcapRingIndex_ [EEDetId::IX_MAX][EEDetId::IY_MAX] |
| static bool | isInitializedFromGeometry_ = false |
Definition at line 20 of file EcalRingCalibrationTools.h.
| EcalRingCalibrationTools::EcalRingCalibrationTools | ( | ) | [inline] |
Definition at line 23 of file EcalRingCalibrationTools.h.
{};
| EcalRingCalibrationTools::~EcalRingCalibrationTools | ( | ) | [inline] |
Definition at line 24 of file EcalRingCalibrationTools.h.
{};
| std::vector< DetId > EcalRingCalibrationTools::getDetIdsInECAL | ( | ) | [static] |
Definition at line 113 of file EcalRingCalibrationTools.cc.
References initializeFromGeometry(), isInitializedFromGeometry_, EEDetId::IX_MAX, EEDetId::IY_MAX, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EEDetId::validDetId(), and EBDetId::validDetId().
Referenced by ZeeCalibration::duringLoop(), and ZeeCalibration::endOfLoop().
{
std::vector<DetId> ringIds;
for(int ieta= - EBDetId::MAX_IETA; ieta<=EBDetId::MAX_IETA; ++ieta)
for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId::MAX_IPHI; ++iphi)
if (EBDetId::validDetId(ieta,iphi))
ringIds.push_back(EBDetId(ieta,iphi));
//needed only for the EE, it can be replaced at some point maybe with something smarter
if (!isInitializedFromGeometry_)
initializeFromGeometry();
for (int ix=0;ix<EEDetId::IX_MAX;++ix)
for (int iy=0;iy<EEDetId::IY_MAX;++iy)
for(int zside = -1; zside<2; zside += 2)
if ( EEDetId::validDetId(ix+1,iy+1,zside) )
ringIds.push_back( EEDetId(ix+1,iy+1,zside) );
// std::cout<<" [EcalRingCalibrationTools::getDetIdsInECAL()] DetId.size() is "<<ringIds.size()<<std::endl;
return ringIds;
}
| std::vector< DetId > EcalRingCalibrationTools::getDetIdsInModule | ( | short | int | ) | [static] |
Definition at line 139 of file EcalRingCalibrationTools.cc.
References gather_cfg::cout, and N_MODULES_BARREL.
Referenced by ZeeCalibration::duringLoop(), and ZeeCalibration::endOfLoop().
{
std::vector<DetId> ringIds;
if (moduleIndex < 0)
return ringIds;
short moduleBound[5] = {1, 26, 46, 66, 86};
if ( moduleIndex < EcalRingCalibrationTools::N_MODULES_BARREL)
{
short sm, moduleInSm, zsm;
short minModuleiphi, maxModuleiphi, minModuleieta=360,maxModuleieta=0;
// if(moduleIndex%4 != 0 )
sm = moduleIndex / 4 + 1;
// else
//sm = moduleIndex/4;//i.e. module 8 belongs to sm=3, not sm=3
//if(moduleIndex%4 != 0 )
moduleInSm = moduleIndex%4;
//else
//moduleInSm = 4;//moduleInSm is [1,2,3,4]
if(moduleIndex > 71)
zsm = -1;
else
zsm = 1;
minModuleiphi = ( (sm - 1) %18 + 1 ) *20 - 19;
maxModuleiphi = ( (sm - 1) %18 + 1 ) * 20;
if(zsm == 1)
{
minModuleieta = moduleBound[ moduleInSm ];
maxModuleieta = moduleBound[ moduleInSm + 1 ] - 1;
}
else if(zsm == -1){
minModuleieta = - moduleBound[ moduleInSm + 1 ] + 1;
maxModuleieta = - moduleBound[ moduleInSm ];
}
std::cout<<"Called moduleIndex "<<moduleIndex<<std::endl;
std::cout<<"minModuleieta "<<minModuleieta<<" maxModuleieta "<<maxModuleieta<<" minModuleiphi "<<minModuleiphi<<" maxModuleiphi "<<maxModuleiphi<<std::endl;
for(int ieta = minModuleieta; ieta <= maxModuleieta; ++ieta){
for(int iphi = minModuleiphi; iphi<= maxModuleiphi; ++iphi){
ringIds.push_back(EBDetId(ieta,iphi));
// std::cout<<"Putting Xtal with ieta: "<<ieta<<" iphi "<<iphi<<" of SM "<<sm<<" into Module "<<moduleIndex<<std::endl;
}//close loop on phi
}//close loop on eta
}//close if ( moduleInstallationNumber < 144)
return ringIds;
}
| std::vector< DetId > EcalRingCalibrationTools::getDetIdsInRing | ( | short | aRingIndex | ) | [static] |
Retrieve the DetIds in a phi-ring.
Definition at line 71 of file EcalRingCalibrationTools.cc.
References endcapRingIndex_, initializeFromGeometry(), isInitializedFromGeometry_, EEDetId::IX_MAX, EEDetId::IY_MAX, gen::k, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, N_RING_BARREL, N_RING_ENDCAP, N_RING_TOTAL, and EBDetId::validDetId().
Referenced by ZeeCalibration::duringLoop(), and ZeeCalibration::endOfLoop().
{
std::vector<DetId> ringIds;
if (etaIndex < 0)
return ringIds;
if (etaIndex < N_RING_BARREL)
{
int k =0;
if (etaIndex<85)
k=-85 + etaIndex;
else
k= etaIndex - 84;
for(int iphi=EBDetId::MIN_IPHI; iphi<=EBDetId::MAX_IPHI; ++iphi)
if (EBDetId::validDetId(k,iphi))
ringIds.push_back(EBDetId(k,iphi));
}
else if (etaIndex < N_RING_TOTAL)
{
//needed only for the EE, it can be replaced at some point maybe with something smarter
if (!isInitializedFromGeometry_)
initializeFromGeometry();
int zside= (etaIndex < N_RING_BARREL + (N_RING_ENDCAP/2) ) ? -1 : 1;
short eeEtaIndex = (etaIndex - N_RING_BARREL)%(N_RING_ENDCAP/2);
for (int ix=0;ix<EEDetId::IX_MAX;++ix)
for (int iy=0;iy<EEDetId::IY_MAX;++iy)
if (endcapRingIndex_[ix][iy] == eeEtaIndex)
ringIds.push_back(EEDetId(ix+1,iy+1,zside));
}
return ringIds;
}
| short EcalRingCalibrationTools::getModuleIndex | ( | DetId | aDetId | ) | [static] |
Definition at line 46 of file EcalRingCalibrationTools.cc.
References DetId::Ecal, EcalBarrel, EcalEndcap, EBDetId::im(), EBDetId::ism(), and python::rootplot::argparse::module.
Referenced by calib::CalibElectron::getCalibModulesWeights().
{
if (id.det() != DetId::Ecal)
return -1;
if (id.subdetId() == EcalBarrel)
{
short module = 4*( EBDetId(id).ism() -1 ) + EBDetId(id).im() -1;
// std::cout<<"SM construction # : "<<EBDetId(id).ism() <<" SM installation # : "<< installationSMNumber[ EBDetId(id).ism() -1 ]<< "Xtal is Module :"<<module<< std::endl;
return module;
}
if (id.subdetId() == EcalEndcap)
{
return -1;
}
return -1;
}
| short EcalRingCalibrationTools::getRingIndex | ( | DetId | aDetId | ) | [static] |
Retrieve the phi-ring index corresponding to a DetId.
Definition at line 20 of file EcalRingCalibrationTools.cc.
References DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, EBDetId::ieta(), initializeFromGeometry(), isInitializedFromGeometry_, EEDetId::ix(), EEDetId::iy(), N_RING_BARREL, N_RING_ENDCAP, and EEDetId::zside().
Referenced by calib::CalibElectron::getCalibModulesWeights().
{
if (id.det() != DetId::Ecal)
return -1;
if (id.subdetId() == EcalBarrel)
{
if(EBDetId(id).ieta()<0)
return EBDetId(id).ieta() + 85;
else
return EBDetId(id).ieta() + 84;
}
if (id.subdetId() == EcalEndcap)
{
//needed only for the EE, it can be replaced at some point with something smarter
if (!isInitializedFromGeometry_)
initializeFromGeometry();
EEDetId eid(id);
short endcapRingIndex = endcapRingIndex_[eid.ix()-1][eid.iy()-1] + N_RING_BARREL;
if (eid.zside() == 1) endcapRingIndex += N_RING_ENDCAP/2;
return endcapRingIndex;
}
return -1;
}
| void EcalRingCalibrationTools::initializeFromGeometry | ( | ) | [static, private] |
Definition at line 203 of file EcalRingCalibrationTools.cc.
References caloGeometry_, DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, PV3DBase< T, PVType, FrameType >::eta(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), CaloGeometry::getValidDetIds(), h2_mapping_cfi::ics, isInitializedFromGeometry_, EEDetId::ix(), EEDetId::IX_MAX, EEDetId::iy(), EEDetId::IY_MAX, N_RING_ENDCAP, relativeConstraints::ring, and EEDetId::zside().
Referenced by getDetIdsInECAL(), getDetIdsInRing(), and getRingIndex().
{
if (!caloGeometry_)
{
edm::LogError("EcalRingCalibrationTools") << "BIG ERROR::Initializing without geometry handle" ;
return;
}
float m_cellPosEta[EEDetId::IX_MAX][EEDetId::IY_MAX];
for (int ix=0; ix<EEDetId::IX_MAX; ++ix)
for (int iy=0; iy<EEDetId::IY_MAX; ++iy)
{
m_cellPosEta[ix][iy] = -1.;
endcapRingIndex_[ix][iy]=-9;
}
const CaloSubdetectorGeometry *endcapGeometry = caloGeometry_->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
if (!endcapGeometry)
{
edm::LogError("EcalRingCalibrationTools") << "BIG ERROR::Ecal Endcap geometry not found" ;
return;
}
const std::vector<DetId>& m_endcapCells= caloGeometry_->getValidDetIds(DetId::Ecal, EcalEndcap);
for (std::vector<DetId>::const_iterator endcapIt = m_endcapCells.begin();
endcapIt!=m_endcapCells.end();
++endcapIt)
{
EEDetId ee(*endcapIt);
if (ee.zside() == -1) continue; //Just using +side to fill absEta x,y map
const CaloCellGeometry *cellGeometry = endcapGeometry->getGeometry(*endcapIt) ;
int ics=ee.ix() - 1 ;
int ips=ee.iy() - 1 ;
m_cellPosEta[ics][ips] = fabs(cellGeometry->getPosition().eta());
//std::cout<<"EE Xtal, |eta| is "<<fabs(cellGeometry->getPosition().eta())<<std::endl;
}
float eta_ring[N_RING_ENDCAP/2];
for (int ring=0; ring<N_RING_ENDCAP/2; ++ring)
eta_ring[ring]=m_cellPosEta[ring][50];
double etaBoundary[N_RING_ENDCAP/2 + 1];
etaBoundary[0]=1.47;
etaBoundary[N_RING_ENDCAP/2]=4.0;
for (int ring=1; ring<N_RING_ENDCAP/2; ++ring)
etaBoundary[ring]=(eta_ring[ring]+eta_ring[ring-1])/2.;
for (int ring=0; ring<N_RING_ENDCAP/2; ring++){
// std::cout<<"***********************EE ring: "<<ring<<" eta "<<(etaBoundary[ring] + etaBoundary[ring+1])/2.<<std::endl;
for (int ix=0; ix<EEDetId::IX_MAX; ix++)
for (int iy=0; iy<EEDetId::IY_MAX; iy++)
if (m_cellPosEta[ix][iy]>etaBoundary[ring] && m_cellPosEta[ix][iy]<etaBoundary[ring+1])
{
endcapRingIndex_[ix][iy]=ring;
//std::cout<<"endcapRing_["<<ix+1<<"]["<<iy+1<<"] = "<<ring<<";"<<std::endl;
}
}
const std::vector<DetId>& m_barrelCells= caloGeometry_->getValidDetIds(DetId::Ecal, EcalBarrel);
for (std::vector<DetId>::const_iterator barrelIt = m_barrelCells.begin();
barrelIt!=m_barrelCells.end();
++barrelIt)
{
EBDetId eb(*barrelIt);
}
//EB
isInitializedFromGeometry_ = true;
}
| static void EcalRingCalibrationTools::setCaloGeometry | ( | const CaloGeometry * | geometry | ) | [inline, static] |
Definition at line 42 of file EcalRingCalibrationTools.h.
References caloGeometry_, and geometry.
Referenced by ZeeCalibration::duringLoop().
{ caloGeometry_ = geometry; };
const CaloGeometry * EcalRingCalibrationTools::caloGeometry_ = 0 [static, private] |
Definition at line 50 of file EcalRingCalibrationTools.h.
Referenced by initializeFromGeometry(), and setCaloGeometry().
short EcalRingCalibrationTools::endcapRingIndex_ [static, private] |
Definition at line 49 of file EcalRingCalibrationTools.h.
Referenced by getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().
bool EcalRingCalibrationTools::isInitializedFromGeometry_ = false [static, private] |
Definition at line 48 of file EcalRingCalibrationTools.h.
Referenced by getDetIdsInECAL(), getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().
const short EcalRingCalibrationTools::N_MODULES_BARREL = 144 [static] |
Definition at line 40 of file EcalRingCalibrationTools.h.
Referenced by calib::CalibElectron::getCalibModulesWeights(), getDetIdsInModule(), and ZIterativeAlgorithmWithFit::ZIterativeAlgorithmWithFit().
const short EcalRingCalibrationTools::N_RING_BARREL = 170 [static] |
Definition at line 37 of file EcalRingCalibrationTools.h.
Referenced by getDetIdsInRing(), and getRingIndex().
const short EcalRingCalibrationTools::N_RING_ENDCAP = 78 [static] |
Definition at line 38 of file EcalRingCalibrationTools.h.
Referenced by getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().
const short EcalRingCalibrationTools::N_RING_TOTAL = 248 [static] |
Definition at line 36 of file EcalRingCalibrationTools.h.
Referenced by calib::CalibElectron::getCalibModulesWeights(), getDetIdsInRing(), and ZIterativeAlgorithmWithFit::ZIterativeAlgorithmWithFit().
1.7.3