|
| 1 | +#include "FWCore/MessageLogger/interface/MessageLogger.h" |
| 2 | +#include "Geometry/HGCalCommonData/interface/HGCalWaferMask.h" |
| 3 | +#include "Geometry/HGCalCommonData/interface/HGCalWaferType.h" |
| 4 | +#include "Geometry/HGCalCommonData/interface/HGCGuardRing.h" |
| 5 | +#include <iostream> |
| 6 | + |
| 7 | +//#define EDM_ML_DEBUG |
| 8 | + |
| 9 | +HGCGuardRing::HGCGuardRing(const HGCalDDDConstants& hgc) |
| 10 | + : hgcons_(hgc), |
| 11 | + modeUV_(hgcons_.geomMode()), |
| 12 | + v17OrLess_(hgcons_.v17OrLess()), |
| 13 | + waferSize_(hgcons_.waferSize(false)), |
| 14 | + sensorSizeOffset_(hgcons_.sensorSizeOffset(false)), |
| 15 | + guardRingOffset_(hgcons_.guardRingOffset(false)) { |
| 16 | + offset_ = sensorSizeOffset_ + guardRingOffset_; |
| 17 | + offsetPartial_ = guardRingOffset_; |
| 18 | + xmax_ = 0.5 * waferSize_ - offset_; |
| 19 | + ymax_ = xmax_ / sqrt3_; |
| 20 | + c22_ = (v17OrLess_) ? HGCalTypes::c22O : HGCalTypes::c22; |
| 21 | + c27_ = (v17OrLess_) ? HGCalTypes::c27O : HGCalTypes::c27; |
| 22 | +#ifdef EDM_ML_DEBUG |
| 23 | + edm::LogVerbatim("HGCSim") << "Creating HGCGuardRing with wafer size " << waferSize_ << ", Offsets " |
| 24 | + << sensorSizeOffset_ << ":" << guardRingOffset_ << ":" << offset_ << ":" << offsetPartial_ |
| 25 | + << ", mode " << modeUV_ << ", xmax|ymax " << xmax_ << ":" << ymax_ << " and c22:c77 " |
| 26 | + << c22_ << ":" << c77_; |
| 27 | +#endif |
| 28 | +} |
| 29 | + |
| 30 | +bool HGCGuardRing::exclude(G4ThreeVector& point, int zside, int frontBack, int layer, int waferU, int waferV) { |
| 31 | + bool check(false); |
| 32 | + if (hgcons_.waferHexagon8Module()) { |
| 33 | + int index = HGCalWaferIndex::waferIndex(layer, waferU, waferV); |
| 34 | + int partial = HGCalWaferType::getPartial(index, hgcons_.getParameter()->waferInfoMap_); |
| 35 | +#ifdef EDM_ML_DEBUG |
| 36 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::exclude: Layer " << layer << " wafer " << waferU << ":" << waferV |
| 37 | + << " index " << index << " partial " << partial; |
| 38 | +#endif |
| 39 | + if (partial == HGCalTypes::WaferFull) { |
| 40 | + double dx = std::abs(point.x()); |
| 41 | + double dy = std::abs(point.y()); |
| 42 | + if (dx > xmax_) { |
| 43 | + check = true; |
| 44 | + } else if (dy > (2 * ymax_)) { |
| 45 | + check = true; |
| 46 | + } else { |
| 47 | + check = (dx > (sqrt3_ * (2 * ymax_ - dy))); |
| 48 | + } |
| 49 | +#ifdef EDM_ML_DEBUG |
| 50 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::exclude: Point " << point << " zside " << zside << " layer " << layer |
| 51 | + << " wafer " << waferU << ":" << waferV << " partial type " << partial << ":" |
| 52 | + << HGCalTypes::WaferFull << " x " << dx << ":" << xmax_ << " y " << dy << ":" << ymax_ |
| 53 | + << " check " << check; |
| 54 | +#endif |
| 55 | + } else if (partial > 0) { |
| 56 | + int orient = HGCalWaferType::getOrient(index, hgcons_.getParameter()->waferInfoMap_); |
| 57 | +#ifdef EDM_ML_DEBUG |
| 58 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::exclude: Orient " << orient << " Mode " << modeUV_; |
| 59 | +#endif |
| 60 | + if (hgcons_.v16OrLess()) { |
| 61 | + std::vector<std::pair<double, double> > wxy = |
| 62 | + HGCalWaferMask::waferXY(partial, orient, zside, waferSize_, offset_, 0.0, 0.0, v17OrLess_); |
| 63 | + check = !(insidePolygon(point.x(), point.y(), wxy)); |
| 64 | +#ifdef EDM_ML_DEBUG |
| 65 | + std::ostringstream st1; |
| 66 | + st1 << "HGCGuardRing::exclude: Point " << point << " Partial/orient/zside/size/offset " << partial << ":" |
| 67 | + << orient << ":" << zside << ":" << waferSize_ << offset_ << " with " << wxy.size() << " points:"; |
| 68 | + for (unsigned int k = 0; k < wxy.size(); ++k) |
| 69 | + st1 << " (" << wxy[k].first << ", " << wxy[k].second << ")"; |
| 70 | + edm::LogVerbatim("HGCSim") << st1.str(); |
| 71 | +#endif |
| 72 | + } else { |
| 73 | + int placement = HGCalCell::cellPlacementIndex(zside, frontBack, orient); |
| 74 | + std::vector<std::pair<double, double> > wxy = |
| 75 | + HGCalWaferMask::waferXY(partial, placement, waferSize_, offset_, 0.0, 0.0, v17OrLess_); |
| 76 | + check = !(insidePolygon(point.x(), point.y(), wxy)); |
| 77 | +#ifdef EDM_ML_DEBUG |
| 78 | + std::ostringstream st1; |
| 79 | + st1 << "HGCGuardRing::exclude: Point " << point << " Partial/frontback/orient/zside/placeemnt/size/offset " |
| 80 | + << partial << ":" << frontBack << ":" << orient << ":" << zside << ":" << placement << ":" << waferSize_ |
| 81 | + << offset_ << " with " << wxy.size() << " points:"; |
| 82 | + for (unsigned int k = 0; k < wxy.size(); ++k) |
| 83 | + st1 << " (" << wxy[k].first << ", " << wxy[k].second << ")"; |
| 84 | + edm::LogVerbatim("HGCSim") << st1.str(); |
| 85 | +#endif |
| 86 | + } |
| 87 | + } else { |
| 88 | + check = true; |
| 89 | + } |
| 90 | + } |
| 91 | + return check; |
| 92 | +} |
| 93 | + |
| 94 | +bool HGCGuardRing::excludePartial(G4ThreeVector& point, int zside, int frontBack, int layer, int waferU, int waferV) { |
| 95 | + bool check(false); |
| 96 | + if (hgcons_.waferHexagon8Cassette()) { |
| 97 | + int index = HGCalWaferIndex::waferIndex(layer, waferU, waferV); |
| 98 | + int partial = HGCalWaferType::getPartial(index, hgcons_.getParameter()->waferInfoMap_); |
| 99 | + int type = HGCalWaferType::getType(index, hgcons_.getParameter()->waferInfoMap_); |
| 100 | +#ifdef EDM_ML_DEBUG |
| 101 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::excludePartial: Layer " << layer << " wafer " << waferU << ":" |
| 102 | + << waferV << " index " << index << " partial " << partial << " type " << type; |
| 103 | +#endif |
| 104 | + if (partial == HGCalTypes::WaferFull) { |
| 105 | + return (check); |
| 106 | + } else if (partial < 0) { |
| 107 | + return true; |
| 108 | + } else { |
| 109 | + int orient = HGCalWaferType::getOrient(index, hgcons_.getParameter()->waferInfoMap_); |
| 110 | + int placement = HGCalCell::cellPlacementIndex(zside, frontBack, orient); |
| 111 | + double dx = point.x(); |
| 112 | + double dy = point.y(); |
| 113 | +#ifdef EDM_ML_DEBUG |
| 114 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::excludePartial: zside " << zside << " frontBack " << frontBack |
| 115 | + << " orient " << orient << " placement " << placement << " dx " << dx << " dy " << dy; |
| 116 | +#endif |
| 117 | + if (type > 0) { |
| 118 | + for (int ii = HGCalTypes::WaferPartLDOffset; |
| 119 | + ii < (HGCalTypes::WaferPartLDOffset + HGCalTypes::WaferPartLDCount); |
| 120 | + ii++) { |
| 121 | + std::array<double, 4> criterion = |
| 122 | + HGCalWaferMask::maskCut(ii, placement, waferSize_, offsetPartial_, v17OrLess_); |
| 123 | + check |= std::abs(criterion[0] * dy + criterion[1] * dx + criterion[2]) < criterion[3]; |
| 124 | + } |
| 125 | + } else { |
| 126 | + for (int ii = HGCalTypes::WaferPartHDOffset; |
| 127 | + ii < (HGCalTypes::WaferPartHDOffset + HGCalTypes::WaferPartHDCount); |
| 128 | + ii++) { |
| 129 | + std::array<double, 4> criterion = |
| 130 | + HGCalWaferMask::maskCut(ii, placement, waferSize_, offsetPartial_, v17OrLess_); |
| 131 | + check |= std::abs(criterion[0] * dy + criterion[1] * dx + criterion[2]) < criterion[3]; |
| 132 | + } |
| 133 | + } |
| 134 | + } |
| 135 | +#ifdef EDM_ML_DEBUG |
| 136 | + edm::LogVerbatim("HGCSim") << "HGCGuardRing::excludePartial: Point " << point << " zside " << zside << " layer " |
| 137 | + << layer << " wafer " << waferU << ":" << waferV << " partial type " << partial |
| 138 | + << " type " << type << " check " << check; |
| 139 | +#endif |
| 140 | + } |
| 141 | + return check; |
| 142 | +} |
| 143 | + |
| 144 | +bool HGCGuardRing::insidePolygon(double x, double y, const std::vector<std::pair<double, double> >& xyv) { |
| 145 | + int counter(0); |
| 146 | + double x1(xyv[0].first), y1(xyv[0].second); |
| 147 | + for (unsigned i1 = 1; i1 <= xyv.size(); i1++) { |
| 148 | + unsigned i2 = (i1 % xyv.size()); |
| 149 | + double x2(xyv[i2].first), y2(xyv[i2].second); |
| 150 | + if (y > std::min(y1, y2)) { |
| 151 | + if (y <= std::max(y1, y2)) { |
| 152 | + if (x <= std::max(x1, x2)) { |
| 153 | + if (y1 != y2) { |
| 154 | + double xinter = (y - y1) * (x2 - x1) / (y2 - y1) + x1; |
| 155 | + if ((x1 == x2) || (x <= xinter)) |
| 156 | + ++counter; |
| 157 | + } |
| 158 | + } |
| 159 | + } |
| 160 | + } |
| 161 | + x1 = x2; |
| 162 | + y1 = y2; |
| 163 | + } |
| 164 | + |
| 165 | + if (counter % 2 == 0) |
| 166 | + return false; |
| 167 | + else |
| 168 | + return true; |
| 169 | +} |
0 commit comments