HGCalGeomTools Class Reference

#include <HGCalGeomTools.h>

Public Member Functions
	HGCalGeomTools ()
std::pair< double, double >	shiftXY (int waferPosition, double waferSize) const
	~HGCalGeomTools ()

Static Public Member Functions
static void	radius (double zf, double zb, std::vector< double > const &zFront1, std::vector< double > const &rFront1, std::vector< double > const &slope1, std::vector< double > const &zFront2, std::vector< double > const &rFront2, std::vector< double > const &slope2, int flag, std::vector< double > &zz, std::vector< double > &rin, std::vector< double > &rout)
static double	radius (double z, std::vector< double > const &zFront, std::vector< double > const &rFront, std::vector< double > const &slope)
static double	radius (double z, int layer0, int layerf, std::vector< double > const &zFront, std::vector< double > const &rFront)
static double	slope (double z, std::vector< double > const &zFront, std::vector< double > const &slope)
static std::pair< int32_t, int32_t >	waferCorner (double xpos, double ypos, double r, double R, double rMin, double rMax, bool oldBug=false)
static std::pair< double, double >	zradius (double z1, double z2, std::vector< double > const &zFront, std::vector< double > const &rFront)

Static Public Attributes
static constexpr int	k_allCorners = 6
static constexpr int	k_fiveCorners = 5
static constexpr int	k_fourCorners = 4
static constexpr int	k_threeCorners = 3
static constexpr int	k_twoCorners = 2

Private Attributes
double	factor_

Static Private Attributes
static constexpr double	tol_ = 0.0001

Detailed Description

Definition at line 8 of file HGCalGeomTools.h.

Constructor & Destructor Documentation

HGCalGeomTools::HGCalGeomTools

(

)

Definition at line 10 of file HGCalGeomTools.cc.

: factor_(1.0 / std::sqrt(3.0)) {}

HGCalGeomTools::~HGCalGeomTools ( )

inline

Definition at line 11 of file HGCalGeomTools.h.

{}

Member Function Documentation

void HGCalGeomTools::radius ( double  zf, double  zb, std::vector< double > const &  zFront1, std::vector< double > const &  rFront1, std::vector< double > const &  slope1, std::vector< double > const &  zFront2, std::vector< double > const &  rFront2, std::vector< double > const &  slope2, int  flag, std::vector< double > &  zz, std::vector< double > &  rin, std::vector< double > &  rout  )

static

Definition at line 12 of file HGCalGeomTools.cc.

References funct::abs(), isotrackApplyRegressor::k, cuda_std::lower_bound(), SiStripPI::max, min(), findQualityFiles::rr, slope(), tol_, and findQualityFiles::v.

Referenced by algorithm(), DDHGCalEEAlgo::constructLayers(), DDHGCalHEFileAlgo::constructLayers(), DDHGCalHEAlgo::constructLayers(), DDHGCalEEFileAlgo::constructLayers(), DDHGCalMixLayer::constructLayers(), DDHGCalSiliconModule::constructLayers(), HGCalEEAlgo::ConstructLayers(), HGCalEEFileAlgo::HGCalEEFileAlgo(), HGCalHEAlgo::HGCalHEAlgo(), HGCalHEFileAlgo::HGCalHEFileAlgo(), HGCalMixLayer::HGCalMixLayer(), HGCalSiliconModule::HGCalSiliconModule(), and HGCalDDDConstants::rangeR().

                                                     {
#ifdef EDM_ML_DEBUG
  std::ostringstream st1;
  st1 << "Z = " << zf << ":" << zb << "; zFront1 =";
  for (const auto& v : zFront1)
    st1 << " " << v;
  st1 << "; rFront1 =";
  for (const auto& v : rFront1)
    st1 << " " << v;
  st1 << "; slope1 =";
  for (const auto& v : slope1)
    st1 << " " << v;
  st1 << "; zFront2 =";
  for (const auto& v : zFront2)
    st1 << " " << v;
  st1 << "; rFront2 =";
  for (const auto& v : rFront2)
    st1 << " " << v;
  st1 << "; slope2 =";
  for (const auto& v : slope2)
    st1 << " " << v;
  edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX: " << st1.str();
  double dz2f(0), dz2b(0);
#endif
  double dz1f(0), dz1b(0);
  auto zf1 = std::lower_bound(zFront1.begin(), zFront1.end(), zf);
  if (zf1 != zFront1.begin())
    --zf1;
  if (((zf1 + 1) != zFront1.end()) && (std::abs(*(zf1 + 1) - zf) < tol_)) {
    ++zf1;
    dz1f = 2 * tol_;
  }
  auto zf2 = std::lower_bound(zFront2.begin(), zFront2.end(), zf);
  if (zf2 != zFront2.begin())
    --zf2;
  if (((zf2 + 1) != zFront2.end()) && (std::abs(*(zf2 + 1) - zf) < tol_)) {
    if (static_cast<int>(zf2 - zFront2.begin()) >= 2)
      ++zf2;
#ifdef EDM_ML_DEBUG
    dz2f = 2 * tol_;
#endif
  }
  auto zb1 = std::lower_bound(zFront1.begin(), zFront1.end(), zb);
  if (zb1 != zFront1.begin())
    --zb1;
  if ((zb1 != zFront1.begin()) && (std::abs(*zb1 - zb) < tol_)) {
    --zb1;
    dz1b = -2 * tol_;
  }
  if (((zb1 + 1) != zFront1.end()) && (std::abs(*(zb1 + 1) - zb) < tol_)) {
    dz1b = -2 * tol_;
  }
  auto zb2 = std::lower_bound(zFront2.begin(), zFront2.end(), zb);
  if (zb2 != zFront2.begin()) {
    --zb2;
  }
  if ((zb2 != zFront2.begin()) && (std::abs(*zb2 - zb) < tol_)) {
    if (static_cast<int>(zb2 - zFront2.begin()) > 2)
      --zb2;
#ifdef EDM_ML_DEBUG
    dz2b = -2 * tol_;
#endif
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX:zf " << zf << " : "
                                << static_cast<int>(zf1 - zFront1.begin()) << ":" << *zf1 << " : "
                                << static_cast<int>(zf2 - zFront2.begin()) << ":" << *zf2 << " zb " << zb << ":"
                                << static_cast<int>(zb1 - zFront1.begin()) << " : " << *zb1 << " : "
                                << static_cast<int>(zb2 - zFront2.begin()) << ":" << *zb2 << " Flag " << flag << ":"
                                << (zf1 == zb1) << ":" << (zf2 == zb2) << " dz " << dz1f << ":" << dz2f << ":" << dz1b
                                << ":" << dz2b;
#endif
  if ((zf1 == zb1) && (zf2 == zb2)) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX:Try 1: " << zf << ":" << zb;
#endif
    zz.emplace_back(zf);
    rin.emplace_back(radius(zf + dz1f, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zf, zFront2, rFront2, slope2));
    zz.emplace_back(zb);
    rin.emplace_back(radius(zb + dz1b, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zb, zFront2, rFront2, slope2));
  } else if (zf1 == zb1) {
#ifdef EDM_ML_DEBUG
    double z1 = std::max(*zf2, *zb2);
    edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX:Try 2:" << zf << ":" << *zb2 << " (" << z1 << ") : " << zb
                                  << " Test " << (slope(*zb2, zFront2, slope2) < tol_) << ":"
                                  << ((std::abs(*zb2 - zb) > tol_) && (std::abs(*zb2 - zf) > tol_));
#endif
    zz.emplace_back(zf);
    rin.emplace_back(radius(zf, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zf, zFront2, rFront2, slope2));
    if (slope(*zb2, zFront2, slope2) < tol_) {
      zz.emplace_back(*zb2);
      rin.emplace_back(radius(*zb2, zFront1, rFront1, slope1));
      rout.emplace_back(radius(*zb2 - tol_, zFront2, rFront2, slope2));
    }
    if ((std::abs(*zb2 - zb) > tol_) && (std::abs(*zb2 - zf) > tol_)) {
      zz.emplace_back(*zb2);
      rin.emplace_back(radius(*zb2, zFront1, rFront1, slope1));
      rout.emplace_back(radius(*zb2, zFront2, rFront2, slope2));
    }
    zz.emplace_back(zb);
    rin.emplace_back(radius(zb, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zb, zFront2, rFront2, slope2));
  } else if (zf2 == zb2) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX:Try 3: " << zf << ":" << *zb1 << ":" << zb << " Test "
                                  << (slope(*zb1, zFront1, slope1) < tol_) << ":"
                                  << ((std::abs(*zb1 - zb) > tol_) && (std::abs(*zb1 - zf) > tol_));
#endif
    zz.emplace_back(zf);
    rin.emplace_back(radius(zf, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zf, zFront2, rFront2, slope2));
    if (slope(*zb1, zFront1, slope1) < tol_) {
      zz.emplace_back(*zb1);
      rin.emplace_back(radius(*zb1 - tol_, zFront1, rFront1, slope1));
      rout.emplace_back(radius(*zb1, zFront2, rFront2, slope2));
    }
    if ((std::abs(*zb1 - zb) > tol_) && (std::abs(*zb1 - zf) > tol_)) {
      zz.emplace_back(*zb1);
      rin.emplace_back(radius(*zb1, zFront1, rFront1, slope1));
      rout.emplace_back(radius(*zb1, zFront2, rFront2, slope2));
    }
    zz.emplace_back(zb);
    rin.emplace_back(radius(zb, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zb, zFront2, rFront2, slope2));
  } else {
    double z1 = std::min(*zf2, *zb1);
    double z2 = std::max(*zf2, *zb1);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX:Try 4: " << zf << ":" << z1 << " : " << z2 << ":" << zb;
#endif
    zz.emplace_back(zf);
    rin.emplace_back(radius(zf, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zf, zFront2, rFront2, slope2));
    zz.emplace_back(z1);
    rin.emplace_back(radius(z1, zFront1, rFront1, slope1));
    rout.emplace_back(radius(z1, zFront2, rFront2, slope2));
    zz.emplace_back(z2);
    rin.emplace_back(radius(z2, zFront1, rFront1, slope1));
    rout.emplace_back(radius(z2, zFront2, rFront2, slope2));
    zz.emplace_back(zb);
    rin.emplace_back(radius(zb, zFront1, rFront1, slope1));
    rout.emplace_back(radius(zb, zFront2, rFront2, slope2));
  }
  double rmin = *(std::min_element(rout.begin(), rout.end()));
  if (flag > 1) {
    for (auto& rr : rout)
      rr = rmin;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radiusX has " << zz.size() << " sections: " << rmin;
  for (unsigned int k = 0; k < zz.size(); ++k)
    edm::LogVerbatim("HGCalGeom") << "[" << k << "] Z = " << zz[k] << " R = " << rin[k] << ":" << rout[k];
#endif
}

double HGCalGeomTools::radius ( double  z, std::vector< double > const &  zFront, std::vector< double > const &  rFront, std::vector< double > const &  slope  )

static

Definition at line 181 of file HGCalGeomTools.cc.

References funct::abs(), cuda_std::lower_bound(), alignCSCRings::r, tol_, findQualityFiles::v, and z.

                                                              {
#ifdef EDM_ML_DEBUG
  std::ostringstream st1;
  st1 << "Z = " << z << "; zFront =";
  for (const auto& v : zFront)
    st1 << " " << v;
  st1 << "; rFront =";
  for (const auto& v : rFront)
    st1 << " " << v;
  st1 << "; slope =";
  for (const auto& v : slope)
    st1 << " " << v;
  edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::radius: " << st1.str();
#endif
  auto itrz = std::lower_bound(zFront.begin(), zFront.end(), z);
  if (itrz != zFront.begin())
    --itrz;
  unsigned int ik = static_cast<unsigned int>(itrz - zFront.begin());
  if ((ik + 1) < zFront.size() && std::abs(z - zFront[ik + 1]) < tol_)
    ++ik;
  double r = rFront[ik] + (z - zFront[ik]) * slope[ik];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalGeomTools: Z " << z << " k " << ik << " R " << r;
#endif
  return r;
}

double HGCalGeomTools::radius ( double  z, int  layer0, int  layerf, std::vector< double > const &  zFront, std::vector< double > const &  rFront  )

static

Definition at line 211 of file HGCalGeomTools.cc.

References isotrackApplyRegressor::k, alignCSCRings::r, and tol_.

                                                                                                        {
  double r = rFront[0];
#ifdef EDM_ML_DEBUG
  unsigned int ik(0);
#endif
  for (unsigned int k = 0; k < rFront.size(); ++k) {
    int k1 = layerf - layer0 + (int)(k);
    if (k1 < (int)(zFront.size())) {
      r = rFront[k];
#ifdef EDM_ML_DEBUG
      ik = k;
#endif
      if (z < zFront[k1] + tol_)
        break;
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalGeomTools: Z " << z << ":" << ik << " R " << r;
#endif
  return r;
}

std::pair< double, double > HGCalGeomTools::shiftXY	(	int	waferPosition,
		double	waferSize
	)		const

Definition at line 234 of file HGCalGeomTools.cc.

References HGCalTypes::CornerCenterXm, HGCalTypes::CornerCenterXp, HGCalTypes::CornerCenterYm, HGCalTypes::CornerCenterYp, PVValHelper::dx, PVValHelper::dy, and factor_.

Referenced by HGCalGeomParameters::loadWaferHexagon8(), and HGCalDDDConstants::locateCell().

                                                                                         {
  double dx(0), dy(0);
  switch (waferPosition) {
    case (HGCalTypes::CornerCenterYp): {
      dy = factor_ * waferSize;
      break;
    }
    case (HGCalTypes::CornerCenterYm): {
      dy = -factor_ * waferSize;
      break;
    }
    case (HGCalTypes::CornerCenterXp): {
      dx = factor_ * waferSize;
      break;
    }
    case (HGCalTypes::CornerCenterXm): {
      dx = -factor_ * waferSize;
      break;
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Shift for " << waferPosition << " is (" << dx << ":" << dy << ")";
#endif
  return std::make_pair(dx, dy);
}

double HGCalGeomTools::slope ( double  z, std::vector< double > const &  zFront, std::vector< double > const &  slope  )

static

Definition at line 260 of file HGCalGeomTools.cc.

References cuda_std::lower_bound(), and z.

Referenced by radius().

                                                                                                      {
  auto itrz = std::lower_bound(zFront.begin(), zFront.end(), z);
  if (itrz != zFront.begin())
    --itrz;
  unsigned int ik = static_cast<unsigned int>(itrz - zFront.begin());
  //  if (ik < zFront.size() && std::abs(z-zFront[ik+1]) < tol_) ++ik;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalGeomTools::slope:z " << z << " k " << ik << " Slope " << slope[ik];
#endif
  return slope[ik];
}

std::pair< int32_t, int32_t > HGCalGeomTools::waferCorner ( double  xpos, double  ypos, double  r, double  R, double  rMin, double  rMax, bool  oldBug = false  )

static

Definition at line 287 of file HGCalGeomTools.cc.

References spr::find(), isotrackApplyRegressor::k, HGCalParameters::k_CornerSize, convertSQLiteXML::ok, dttmaxenums::R, alignCSCRings::r, mathSSE::sqrt(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by algorithm(), HGCalGeomParameters::loadWaferHexagon(), HGCalGeomParameters::loadWaferHexagon8(), DDHGCalEEAlgo::positionSensitive(), DDHGCalModule::positionSensitive(), DDHGCalModuleAlgo::positionSensitive(), DDHGCalEEFileAlgo::positionSensitive(), DDHGCalHEFileAlgo::positionSensitive(), DDHGCalHEAlgo::positionSensitive(), HGCalEEFileAlgo::positionSensitive(), HGCalEEAlgo::PositionSensitive(), HGCalHEFileAlgo::positionSensitive(), HGCalHEAlgo::positionSensitive(), and HGCalDDDConstants::waferInLayerTest().

                                                                                         {
  double xc[HGCalParameters::k_CornerSize], yc[HGCalParameters::k_CornerSize];
  xc[0] = xpos;
  yc[0] = ypos + R;
  xc[1] = xpos - r;
  yc[1] = ypos + 0.5 * R;
  if (oldBug) {
    xc[2] = xpos + r;
    yc[2] = ypos - 0.5 * R;
  } else {
    xc[2] = xpos - r;
    yc[2] = ypos - 0.5 * R;
  }
  xc[3] = xpos;
  yc[3] = ypos - R;
  xc[4] = xpos + r;
  yc[4] = ypos - 0.5 * R;
  xc[5] = xpos + r;
  yc[5] = ypos + 0.5 * R;
  int32_t nCorner(0), firstCorner(-1), firstMiss(-1);
#ifdef EDM_ML_DEBUG
  std::vector<uint32_t> corners;
#endif
  for (uint32_t k = 0; k < HGCalParameters::k_CornerSize; ++k) {
    double rpos = sqrt(xc[k] * xc[k] + yc[k] * yc[k]);
    if ((rpos <= rMax) && (rpos >= rMin)) {
#ifdef EDM_ML_DEBUG
      corners.emplace_back(k);
#endif
      if (firstCorner < 0)
        firstCorner = k;
      ++nCorner;
    } else {
      if (firstMiss < 0)
        firstMiss = k;
    }
  }
  if ((nCorner > 1) && (firstCorner == 0) && (firstMiss < nCorner)) {
    firstCorner = firstMiss + HGCalParameters::k_CornerSize - nCorner;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "waferCorner:: R " << rMin << ":" << rMax << nCorner << " corners; first corner "
                                << firstCorner;
  for (uint32_t k = 0; k < HGCalParameters::k_CornerSize; ++k) {
    double rpos = std::sqrt(xc[k] * xc[k] + yc[k] * yc[k]);
    std::string ok = (std::find(corners.begin(), corners.end(), k) != corners.end()) ? " In" : " Out";
    edm::LogVerbatim("HGCalGeom") << "Corner[" << k << "] x " << xc[k] << " y " << yc[k] << " R " << rpos << ok;
  }
#endif
  return std::make_pair(nCorner, firstCorner);
}

std::pair< double, double > HGCalGeomTools::zradius ( double  z1, double  z2, std::vector< double > const &  zFront, std::vector< double > const &  rFront  )

static

Definition at line 272 of file HGCalGeomTools.cc.

References isotrackApplyRegressor::k, alignCSCRings::r, tol_, and z.

                                                                               {
  double z(-1), r(-1);
  for (unsigned int k = 0; k < rF.size(); ++k) {
    if ((z1 > zF[k] - tol_) && (z2 < zF[k] + tol_)) {
      z = zF[k];
      r = rF[k];
      break;
    }
  }
  return std::make_pair(z, r);
}

Member Data Documentation

double HGCalGeomTools::factor_

private

Definition at line 48 of file HGCalGeomTools.h.

Referenced by shiftXY().

constexpr int HGCalGeomTools::k_allCorners = 6

static

Definition at line 13 of file HGCalGeomTools.h.

Referenced by HGCalWaferMask::getTypeMode(), HGCalDDDConstants::waferType(), and HGCalDDDConstants::waferTypeRotation().

constexpr int HGCalGeomTools::k_fiveCorners = 5

static

Definition at line 14 of file HGCalGeomTools.h.

Referenced by HGCalWaferMask::getTypeMode(), HGCalDDDConstants::waferType(), and HGCalDDDConstants::waferTypeRotation().

constexpr int HGCalGeomTools::k_fourCorners = 4

static

Definition at line 15 of file HGCalGeomTools.h.

Referenced by HGCalWaferMask::getTypeMode(), HGCalWaferMask::maskCell(), HGCalDDDConstants::waferType(), and HGCalDDDConstants::waferTypeRotation().

constexpr int HGCalGeomTools::k_threeCorners = 3

static

Definition at line 16 of file HGCalGeomTools.h.

Referenced by HGCalWaferMask::getTypeMode(), HGCalWaferMask::maskCell(), HGCalDDDConstants::waferType(), and HGCalDDDConstants::waferTypeRotation().

constexpr int HGCalGeomTools::k_twoCorners = 2

static

Definition at line 17 of file HGCalGeomTools.h.

Referenced by HGCalWaferMask::getTypeMode().

constexpr double HGCalGeomTools::tol_ = 0.0001

staticprivate

Definition at line 47 of file HGCalGeomTools.h.

Referenced by radius(), and zradius().