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 23 of file EcalRingCalibrationTools.h.

Constructor & Destructor Documentation

EcalRingCalibrationTools()

EcalRingCalibrationTools::EcalRingCalibrationTools ( )

inline

Definition at line 25 of file EcalRingCalibrationTools.h.

{}

~EcalRingCalibrationTools()

EcalRingCalibrationTools::~EcalRingCalibrationTools ( )

inline

Definition at line 26 of file EcalRingCalibrationTools.h.

{}

Member Function Documentation

getDetIdsInECAL()

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

static

Definition at line 102 of file EcalRingCalibrationTools.cc.

References hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, isInitializedFromGeometry_, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::IX_MAX, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), 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;
}

getDetIdsInModule()

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

static

Definition at line 125 of file EcalRingCalibrationTools.cc.

References gather_cfg::cout, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, N_MODULES_BARREL, and ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::sm().

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

getDetIdsInRing()

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

static

Retrieve the DetIds in a phi-ring.

Definition at line 66 of file EcalRingCalibrationTools.cc.

References endcapRingIndex_, hcalRecHitTable_cff::iphi, isInitializedFromGeometry_, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::IX_MAX, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), EEDetId::IY_MAX, isotrackApplyRegressor::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;
}

getModuleIndex()

short EcalRingCalibrationTools::getModuleIndex ( DetId  aDetId )

static

Definition at line 49 of file EcalRingCalibrationTools.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, EBDetId::im(), EBDetId::ism(), and callgraph::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;
}

getRingIndex()

short EcalRingCalibrationTools::getRingIndex ( DetId  aDetId )

static

Retrieve the phi-ring index corresponding to a DetId.

Definition at line 24 of file EcalRingCalibrationTools.cc.

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

Referenced by calib::CalibElectron::getCalibModulesWeights(), EcalPhiSymRecHitProducerBase::initializeConfigCache(), EcalPhiSymRecHitProducerBase::initializePhiSymCache(), PFRecHitQTestECALMultiThreshold::pass(), and EcalPhiSymRecHitProducerBase::processEvent().

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

initializeFromGeometry()

void EcalRingCalibrationTools::initializeFromGeometry ( CaloGeometry const *  geometry )

staticprivate

Definition at line 185 of file EcalRingCalibrationTools.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, endcapRingIndex_, unpacker2007h4_cff::ics, isInitializedFromGeometry_, EEDetId::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::IX_MAX, EEDetId::iy(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::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
    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;
}

setCaloGeometry()

void EcalRingCalibrationTools::setCaloGeometry ( const CaloGeometry *  geometry )

static

Definition at line 181 of file EcalRingCalibrationTools.cc.

References initializeFromGeometry(), and once_.

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

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

Member Data Documentation

endcapRingIndex_

short EcalRingCalibrationTools::endcapRingIndex_

staticprivate

Definition at line 52 of file EcalRingCalibrationTools.h.

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

isInitializedFromGeometry_

std::atomic< bool > EcalRingCalibrationTools::isInitializedFromGeometry_

staticprivate

Definition at line 49 of file EcalRingCalibrationTools.h.

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

N_MODULES_BARREL

constexpr short EcalRingCalibrationTools::N_MODULES_BARREL = 144

static

Definition at line 42 of file EcalRingCalibrationTools.h.

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

N_RING_BARREL

constexpr short EcalRingCalibrationTools::N_RING_BARREL = 170

static

Definition at line 39 of file EcalRingCalibrationTools.h.

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

N_RING_ENDCAP

constexpr short EcalRingCalibrationTools::N_RING_ENDCAP = 78

static

Definition at line 40 of file EcalRingCalibrationTools.h.

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

N_RING_TOTAL

constexpr short EcalRingCalibrationTools::N_RING_TOTAL = 248

static

Definition at line 38 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 54 of file EcalRingCalibrationTools.h.

Referenced by setCaloGeometry().