HGCGuardRing Class Reference

#include <HGCGuardRing.h>

Public Member Functions
bool	exclude (G4ThreeVector &point, int zside, int frontBack, int layer, int waferU, int waferV)
	HGCGuardRing (const HGCalDDDConstants &hgc)

Static Public Member Functions
static bool	insidePolygon (double x, double y, const std::vector< std::pair< double, double > > &xyv)

Private Attributes
const double	guardRingOffset_
const HGCalDDDConstants &	hgcons_
const HGCalGeometryMode::GeometryMode	modeUV_
double	offset_
const double	sensorSizeOffset_
const double	waferSize_
double	xmax_
double	ymax_

Static Private Attributes
static constexpr double	sqrt3_ = 1.732050807568877

Detailed Description

Definition at line 9 of file HGCGuardRing.h.

Constructor & Destructor Documentation

HGCGuardRing()

HGCGuardRing::HGCGuardRing

(

const HGCalDDDConstants &

hgc

)

Definition at line 9 of file HGCGuardRing.cc.

References guardRingOffset_, modeUV_, offset_, sensorSizeOffset_, sqrt3_, waferSize_, xmax_, and ymax_.

    : hgcons_(hgc),
      modeUV_(hgcons_.geomMode()),
      waferSize_(hgcons_.waferSize(false)),
      sensorSizeOffset_(hgcons_.getParameter()->sensorSizeOffset_),
      guardRingOffset_(hgcons_.getParameter()->guardRingOffset_) {
  offset_ = sensorSizeOffset_ + 2.0 * guardRingOffset_;
  xmax_ = 0.5 * (waferSize_ - offset_);
  ymax_ = xmax_ / sqrt3_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "Creating HGCGuardRing with wafer size " << waferSize_ << ", Offsets "
                             << sensorSizeOffset_ << ":" << guardRingOffset_ << ":" << offset_ << ", and mode "
                             << modeUV_ << " xmax|ymax " << xmax_ << ":" << ymax_;
#endif
}

Member Function Documentation

exclude()

bool HGCGuardRing::exclude	(	G4ThreeVector &	point,
		int	zside,
		int	frontBack,
		int	layer,
		int	waferU,
		int	waferV
	)

Definition at line 25 of file HGCGuardRing.cc.

References funct::abs(), HGCalCell::cellPlacementIndex(), RPCNoise_example::check, PVValHelper::dx, PVValHelper::dy, dqmdumpme::first, HGCalWaferType::getOrient(), HGCalDDDConstants::getParameter(), HGCalWaferType::getPartial(), HGCalGeometryMode::Hexagon8Cassette, HGCalGeometryMode::Hexagon8Module, hgcons_, insidePolygon(), dqmdumpme::k, nano_mu_digi_cff::layer, modeUV_, offset_, hgcalPerformanceValidation::orient, point, edm::second(), sqrt3_, HGCalTypes::WaferFull, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, waferSize_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), HGCalWaferMask::waferXY(), xmax_, ymax_, and ecaldqm::zside().

                                                                                                            {
  bool check(false);
  if ((modeUV_ == HGCalGeometryMode::Hexagon8Module) || (modeUV_ == HGCalGeometryMode::Hexagon8Cassette)) {
    int index = HGCalWaferIndex::waferIndex(layer, waferU, waferV);
    int partial = HGCalWaferType::getPartial(index, hgcons_.getParameter()->waferInfoMap_);
    if (partial == HGCalTypes::WaferFull) {
      double dx = std::abs(point.x());
      double dy = std::abs(point.y());
      if (dx > xmax_) {
        check = true;
      } else if (dy > (2 * ymax_)) {
        check = true;
      } else {
        check = (dx > (sqrt3_ * (2 * ymax_ - dy)));
      }
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCSim") << "HGCGuardRing:: Point " << point << " zside " << zside << " layer " << layer
                                 << " wafer " << waferU << ":" << waferV << " partial type " << partial << ":"
                                 << HGCalTypes::WaferFull << " x " << dx << ":" << xmax_ << " y " << dy << ":" << ymax_
                                 << " check " << check;
#endif
    } else {
      int orient = HGCalWaferType::getOrient(index, hgcons_.getParameter()->waferInfoMap_);
      if (modeUV_ == HGCalGeometryMode::Hexagon8Module) {
        std::vector<std::pair<double, double> > wxy =
            HGCalWaferMask::waferXY(partial, orient, zside, waferSize_, offset_, 0.0, 0.0);
        check = !(insidePolygon(point.x(), point.y(), wxy));
#ifdef EDM_ML_DEBUG
        std::ostringstream st1;
        st1 << "HGCGuardRing:: Point " << point << " Partial/orient/zside/size/offset " << partial << ":" << orient
            << ":" << zside << ":" << waferSize_ << offset_ << " with " << wxy.size() << " points:";
        for (unsigned int k = 0; k < wxy.size(); ++k)
          st1 << " (" << wxy[k].first << ", " << wxy[k].second << ")";
        edm::LogVerbatim("HGCSim") << st1.str();
#endif
      } else {
        int placement = HGCalCell::cellPlacementIndex(zside, frontBack, orient);
        std::vector<std::pair<double, double> > wxy =
            HGCalWaferMask::waferXY(partial, placement, waferSize_, offset_, 0.0, 0.0);
        check = !(insidePolygon(point.x(), point.y(), wxy));
#ifdef EDM_ML_DEBUG
        std::ostringstream st1;
        st1 << "HGCGuardRing:: Point " << point << " Partial/frontback/orient/zside/placeemnt/size/offset " << partial
            << ":" << frontBack << ":" << orient << ":" << zside << ":" << placement << ":" << waferSize_ << offset_
            << " with " << wxy.size() << " points:";
        for (unsigned int k = 0; k < wxy.size(); ++k)
          st1 << " (" << wxy[k].first << ", " << wxy[k].second << ")";
        edm::LogVerbatim("HGCSim") << st1.str();
#endif
      }
    }
  }
  return check;
}

insidePolygon()

bool HGCGuardRing::insidePolygon ( double  x, double  y, const std::vector< std::pair< double, double > > &  xyv  )

static

Definition at line 80 of file HGCGuardRing.cc.

References counter, dqmdumpme::first, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, SiStripPI::max, SiStripPI::min, edm::second(), x, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, y, testProducerWithPsetDescEmpty_cfi::y1, and testProducerWithPsetDescEmpty_cfi::y2.

Referenced by HGCalTestGuardRing::analyze(), and exclude().

                                                                                                   {
  int counter(0);
  double x1(xyv[0].first), y1(xyv[0].second);
  for (unsigned i1 = 1; i1 <= xyv.size(); i1++) {
    unsigned i2 = (i1 % xyv.size());
    double x2(xyv[i2].first), y2(xyv[i2].second);
    if (y > std::min(y1, y2)) {
      if (y <= std::max(y1, y2)) {
        if (x <= std::max(x1, x2)) {
          if (y1 != y2) {
            double xinter = (y - y1) * (x2 - x1) / (y2 - y1) + x1;
            if ((x1 == x2) || (x <= xinter))
              ++counter;
          }
        }
      }
    }
    x1 = x2;
    y1 = y2;
  }

  if (counter % 2 == 0)
    return false;
  else
    return true;
}

Member Data Documentation

guardRingOffset_

const double HGCGuardRing::guardRingOffset_

private

Definition at line 19 of file HGCGuardRing.h.

Referenced by HGCGuardRing().

hgcons_

const HGCalDDDConstants& HGCGuardRing::hgcons_

private

Definition at line 17 of file HGCGuardRing.h.

Referenced by exclude().

modeUV_

const HGCalGeometryMode::GeometryMode HGCGuardRing::modeUV_

private

Definition at line 18 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().

offset_

double HGCGuardRing::offset_

private

Definition at line 20 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().

sensorSizeOffset_

const double HGCGuardRing::sensorSizeOffset_

private

Definition at line 19 of file HGCGuardRing.h.

Referenced by HGCGuardRing().

sqrt3_

constexpr double HGCGuardRing::sqrt3_ = 1.732050807568877

staticprivate

Definition at line 16 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().

waferSize_

const double HGCGuardRing::waferSize_

private

Definition at line 19 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().

xmax_

double HGCGuardRing::xmax_

private

Definition at line 20 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().

ymax_

double HGCGuardRing::ymax_

private

Definition at line 20 of file HGCGuardRing.h.

Referenced by exclude(), and HGCGuardRing().