EcalRingCalibrationTools Class Reference

#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 (CaloGeometry const *geometry)

Static Private Attributes
static short	endcapRingIndex_ [EEDetId::IX_MAX][EEDetId::IY_MAX]
static std::atomic< bool >	isInitializedFromGeometry_
static std::once_flag	once_

Detailed Description

Definition at line 22 of file EcalRingCalibrationTools.h.

Constructor & Destructor Documentation

EcalRingCalibrationTools::EcalRingCalibrationTools ( )

inline

Definition at line 25 of file EcalRingCalibrationTools.h.

{};

EcalRingCalibrationTools::~EcalRingCalibrationTools ( )

inline

Definition at line 26 of file EcalRingCalibrationTools.h.

{};

Member Function Documentation

std::vector< DetId > EcalRingCalibrationTools::getDetIdsInECAL ( )

static

Definition at line 112 of file EcalRingCalibrationTools.cc.

References isInitializedFromGeometry_, EEDetId::IX_MAX, EEDetId::IY_MAX, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, EBDetId::validDetId(), EEDetId::validDetId(), and ecaldqm::zside().

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 (not isInitializedFromGeometry_)
    throw std::logic_error("EcalRingCalibrationTools::initializeFromGeometry Ecal Endcap geometry is not initialized");
  
  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 137 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_, isInitializedFromGeometry_, EEDetId::IX_MAX, EEDetId::IY_MAX, relval_2017::k, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, N_RING_BARREL, N_RING_ENDCAP, N_RING_TOTAL, EBDetId::validDetId(), and ecaldqm::zside().

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 (not isInitializedFromGeometry_)
        throw std::logic_error("EcalRingCalibrationTools::initializeFromGeometry Ecal Endcap geometry is not initialized");
      
      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 19 of file EcalRingCalibrationTools.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, EBDetId::ieta(), 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 (not isInitializedFromGeometry_)
        throw std::logic_error("EcalRingCalibrationTools::initializeFromGeometry Ecal Endcap geometry is not initialized");

      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 ( CaloGeometry const *  geometry )

staticprivate

Definition at line 206 of file EcalRingCalibrationTools.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, PV3DBase< T, PVType, FrameType >::eta(), CaloCellGeometry::getPosition(), 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 setCaloGeometry().

{
  if (not geometry)
    throw std::invalid_argument("EcalRingCalibrationTools::initializeFromGeometry called with a nullptr argument");

  float 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)
    {
      cellPosEta[ix][iy] = -1.;
      endcapRingIndex_[ix][iy]=-9;
    }
  
  CaloSubdetectorGeometry const* endcapGeometry = geometry->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
  if (not endcapGeometry)
    throw std::logic_error("EcalRingCalibrationTools::initializeFromGeometry Ecal Endcap geometry not found");

  std::vector<DetId> const& endcapCells = geometry->getValidDetIds(DetId::Ecal, EcalEndcap);

  for (std::vector<DetId>::const_iterator endcapIt = endcapCells.begin();
       endcapIt!=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 ;
      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]=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 (cellPosEta[ix][iy]>etaBoundary[ring] && cellPosEta[ix][iy]<etaBoundary[ring+1])
    {
      endcapRingIndex_[ix][iy]=ring;
      //std::cout<<"endcapRing_["<<ix+1<<"]["<<iy+1<<"] = "<<ring<<";"<<std::endl;  
    }
  }

  std::vector<DetId> const& barrelCells = geometry->getValidDetIds(DetId::Ecal, EcalBarrel);

  for (std::vector<DetId>::const_iterator barrelIt = barrelCells.begin();
       barrelIt!=barrelCells.end();
       ++barrelIt)
    {
      EBDetId eb(*barrelIt);
    }

  //EB

  isInitializedFromGeometry_ = true;

}

void EcalRingCalibrationTools::setCaloGeometry ( const CaloGeometry *  geometry )

static

Definition at line 201 of file EcalRingCalibrationTools.cc.

References initializeFromGeometry(), and once_.

Referenced by ZeeCalibration::duringLoop(), and HLTEcalPhiSymFilter::filter().

{ 
  std::call_once(once_, EcalRingCalibrationTools::initializeFromGeometry, geometry);
}

Member Data Documentation

short EcalRingCalibrationTools::endcapRingIndex_

staticprivate

Definition at line 52 of file EcalRingCalibrationTools.h.

Referenced by getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().

std::atomic< bool > EcalRingCalibrationTools::isInitializedFromGeometry_

staticprivate

Definition at line 49 of file EcalRingCalibrationTools.h.

Referenced by getDetIdsInECAL(), getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().

const short EcalRingCalibrationTools::N_MODULES_BARREL = 144

static

Definition at line 42 of file EcalRingCalibrationTools.h.

Referenced by calib::CalibElectron::getCalibModulesWeights(), getDetIdsInModule(), and ZIterativeAlgorithmWithFit::ZIterativeAlgorithmWithFit().

const short EcalRingCalibrationTools::N_RING_BARREL = 170

static

Definition at line 39 of file EcalRingCalibrationTools.h.

Referenced by HLTEcalPhiSymFilter::filter(), getDetIdsInRing(), and getRingIndex().

const short EcalRingCalibrationTools::N_RING_ENDCAP = 78

static

Definition at line 40 of file EcalRingCalibrationTools.h.

Referenced by HLTEcalPhiSymFilter::filter(), getDetIdsInRing(), getRingIndex(), and initializeFromGeometry().

const short EcalRingCalibrationTools::N_RING_TOTAL = 248

static

Definition at line 38 of file EcalRingCalibrationTools.h.

Referenced by calib::CalibElectron::getCalibModulesWeights(), getDetIdsInRing(), and ZIterativeAlgorithmWithFit::ZIterativeAlgorithmWithFit().

std::once_flag EcalRingCalibrationTools::once_

staticprivate

Definition at line 54 of file EcalRingCalibrationTools.h.

Referenced by setCaloGeometry().