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 constexpr short	N_MODULES_BARREL = 144
static constexpr short	N_RING_BARREL = 170
static constexpr short	N_RING_ENDCAP = 78
static constexpr 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::EcalRingCalibrationTools ( )

inline

Definition at line 24 of file EcalRingCalibrationTools.h.

:
  EcalRingCalibrationTools(){};

~EcalRingCalibrationTools()

EcalRingCalibrationTools::~EcalRingCalibrationTools ( )

inline

Definition at line 25 of file EcalRingCalibrationTools.h.

{};

Member Function Documentation

getDetIdsInECAL()

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

static

Definition at line 102 of file EcalRingCalibrationTools.cc.

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

References LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, 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().

getDetIdsInModule()

std::vector< DetId > EcalRingCalibrationTools::getDetIdsInModule ( short int  )

static

Definition at line 125 of file EcalRingCalibrationTools.cc.

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

References gather_cfg::cout, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and N_MODULES_BARREL.

Referenced by ZeeCalibration::duringLoop(), and ZeeCalibration::endOfLoop().

getDetIdsInRing()

std::vector< DetId > EcalRingCalibrationTools::getDetIdsInRing ( short  aRingIndex )

static

Retrieve the DetIds in a phi-ring.

Definition at line 66 of file EcalRingCalibrationTools.cc.

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

References endcapRingIndex_, LEDCalibrationChannels::iphi, isInitializedFromGeometry_, EEDetId::IX_MAX, EEDetId::IY_MAX, dqmdumpme::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().

getModuleIndex()

short EcalRingCalibrationTools::getModuleIndex ( DetId  aDetId )

static

Definition at line 49 of file EcalRingCalibrationTools.cc.

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

References DetId::Ecal, EcalBarrel, EcalEndcap, EBDetId::im(), and EBDetId::ism().

Referenced by calib::CalibElectron::getCalibModulesWeights().

getRingIndex()

short EcalRingCalibrationTools::getRingIndex ( DetId  aDetId )

static

Retrieve the phi-ring index corresponding to a DetId.

Definition at line 24 of file EcalRingCalibrationTools.cc.

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

References DetId::Ecal, EcalBarrel, EcalEndcap, runTauDisplay::eid, endcapRingIndex_, EBDetId::ieta(), LEDCalibrationChannels::ieta, isInitializedFromGeometry_, N_RING_BARREL, and N_RING_ENDCAP.

Referenced by calib::CalibElectron::getCalibModulesWeights(), and PFRecHitQTestECALMultiThreshold::pass().

initializeFromGeometry()

void EcalRingCalibrationTools::initializeFromGeometry ( CaloGeometry const *  geometry )

staticprivate

Definition at line 185 of file EcalRingCalibrationTools.cc.

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

References DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, Ecal07UnpackerData_cfi::ics, isInitializedFromGeometry_, EEDetId::ix(), EEDetId::IX_MAX, EEDetId::iy(), EEDetId::IY_MAX, N_RING_ENDCAP, relativeConstraints::ring, and EEDetId::zside().

Referenced by setCaloGeometry().

setCaloGeometry()

void EcalRingCalibrationTools::setCaloGeometry ( const CaloGeometry *  geometry )

static

Definition at line 181 of file EcalRingCalibrationTools.cc.

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

References initializeFromGeometry(), and once_.

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

Member Data Documentation

endcapRingIndex_

short EcalRingCalibrationTools::endcapRingIndex_

staticprivate

Definition at line 51 of file EcalRingCalibrationTools.h.

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

isInitializedFromGeometry_

std::atomic< bool > EcalRingCalibrationTools::isInitializedFromGeometry_

staticprivate

Definition at line 48 of file EcalRingCalibrationTools.h.

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

N_MODULES_BARREL

constexpr short EcalRingCalibrationTools::N_MODULES_BARREL = 144

staticconstexpr

Definition at line 41 of file EcalRingCalibrationTools.h.

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

N_RING_BARREL

constexpr short EcalRingCalibrationTools::N_RING_BARREL = 170

staticconstexpr

Definition at line 38 of file EcalRingCalibrationTools.h.

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

N_RING_ENDCAP

constexpr short EcalRingCalibrationTools::N_RING_ENDCAP = 78

staticconstexpr

Definition at line 39 of file EcalRingCalibrationTools.h.

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

N_RING_TOTAL

constexpr short EcalRingCalibrationTools::N_RING_TOTAL = 248

staticconstexpr

Definition at line 37 of file EcalRingCalibrationTools.h.

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

once_

std::once_flag EcalRingCalibrationTools::once_

staticprivate

Definition at line 53 of file EcalRingCalibrationTools.h.

Referenced by setCaloGeometry().