#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. More...

static short	getModuleIndex (DetId aDetId)

static short	getRingIndex (DetId aDetId)
	Retrieve the phi-ring index corresponding to a DetId. More...

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

Detailed Description

Definition at line 20 of file EcalRingCalibrationTools.h.

Constructor & Destructor Documentation

EcalRingCalibrationTools::EcalRingCalibrationTools ( )

inline

Definition at line 23 of file EcalRingCalibrationTools.h.

23 {};

EcalRingCalibrationTools::~EcalRingCalibrationTools ( )

inline

Definition at line 24 of file EcalRingCalibrationTools.h.

24 {};

Member Function Documentation

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, EBDetId::validDetId(), and EEDetId::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 ( )

staticprivate

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