CMS 3D CMS Logo

List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes | Static Private Attributes
HGCalDDDConstants Class Reference

#include <HGCalDDDConstants.h>

Public Member Functions

std::pair< int, int > assignCell (float x, float y, int lay, int subSec, bool reco) const
 
std::array< int, 5 > assignCellHex (float x, float y, int zside, int lay, bool reco, bool extend, bool debug) const
 
std::array< int, 3 > assignCellTrap (float x, float y, float z, int lay, bool reco) const
 
bool cassetteShiftScintillator (int layer, int iphi) const
 
bool cassetteShiftSilicon (int layer, int waferU, int waferV) const
 
int cassetteTile (int iphi) const
 
std::pair< double, double > cellEtaPhiTrap (int type, int irad) const
 
bool cellInLayer (int waferU, int waferV, int cellU, int cellV, int lay, int zside, bool reco) const
 
double cellSizeHex (int type) const
 
std::pair< double, double > cellSizeTrap (int type, int irad) const
 
double cellThickness (int layer, int waferU, int waferV) const
 
int32_t cellType (int type, int waferU, int waferV, int iz, int fwdBack, int orient) const
 
double distFromEdgeHex (double x, double y, double z) const
 
double distFromEdgeTrap (double x, double y, double z) const
 
void etaPhiFromPosition (const double x, const double y, const double z, const int layer, int &ieta, int &iphi, int &type, double &wt) const
 
int firstLayer () const
 
HGCalGeometryMode::GeometryMode geomMode () const
 
int getLayer (double z, bool reco) const
 
int getLayerOffset () const
 
HGCalParameters::hgtrap getModule (unsigned int k, bool hexType, bool reco) const
 
std::vector< HGCalParameters::hgtrapgetModules () const
 
const HGCalParametersgetParameter () const
 
int getPhiBins (int lay) const
 
const std::vector< double > & getRadiusLayer (int layer) const
 
std::pair< double, double > getRangeR (int, bool reco) const
 
std::pair< int, int > getREtaRange (int lay) const
 
HGCalParameters::hgtrform getTrForm (unsigned int k) const
 
unsigned int getTrFormN () const
 
std::vector< HGCalParameters::hgtrformgetTrForms () const
 
int getTypeHex (int layer, int waferU, int waferV) const
 
int getTypeTrap (int layer) const
 
int getUVMax (int type) const
 
std::pair< double, double > getXY (int layer, double x, double y, bool forwd) const
 
 HGCalDDDConstants (const HGCalParameters *hp, const std::string &name)
 
bool isHalfCell (int waferType, int cell) const
 
bool isValidHex (int lay, int mod, int cell, bool reco) const
 
bool isValidHex8 (int lay, int waferU, int waferV, bool fullAndPart) const
 
bool isValidHex8 (int lay, int modU, int modV, int cellU, int cellV, bool fullAndPart) const
 
bool isValidTrap (int zside, int lay, int ieta, int iphi) const
 
int lastLayer (bool reco) const
 
int layerIndex (int lay, bool reco) const
 
unsigned int layers (bool reco) const
 
unsigned int layersInit (bool reco) const
 
int layerType (int lay) const
 
int levelTop (int ind=0) const
 
std::pair< float, float > localToGlobal8 (int zside, int lay, int waferU, int waferV, double localX, double localY, bool reco, bool debug) const
 
std::pair< float, float > locateCell (int cell, int lay, int type, bool reco) const
 
std::pair< float, float > locateCell (int zside, int lay, int waferU, int waferV, int cellU, int cellV, bool reco, bool all, bool norot, bool debug) const
 
std::pair< float, float > locateCell (const HGCSiliconDetId &, bool debug) const
 
std::pair< float, float > locateCell (const HGCScintillatorDetId &, bool debug) const
 
std::pair< float, float > locateCellHex (int cell, int wafer, bool reco) const
 
std::pair< float, float > locateCellTrap (int zside, int lay, int ieta, int iphi, bool reco, bool debug) const
 
bool maskCell (const DetId &id, int corners) const
 
int maxCells (bool reco) const
 
int maxCells (int lay, bool reco) const
 
int maxCellUV () const
 
int maxModules () const
 
int maxModulesPerLayer () const
 
int maxRows (int lay, bool reco) const
 
double minSlope () const
 
int modifyUV (int uv, int type1, int type2) const
 
int modules (int lay, bool reco) const
 
int modulesInit (int lay, bool reco) const
 
double mouseBite (bool reco) const
 
int numberCells (bool reco) const
 
std::vector< int > numberCells (int lay, bool reco) const
 
int numberCellsHexagon (int wafer) const
 
int numberCellsHexagon (int lay, int waferU, int waferV, bool flag) const
 
std::pair< double, double > rangeR (double z, bool reco) const
 
std::pair< double, double > rangeRLayer (int lay, bool reco) const
 
std::pair< double, double > rangeZ (bool reco) const
 
std::pair< int, int > rowColumnWafer (const int wafer) const
 
int sectors () const
 
std::pair< int, int > simToReco (int cell, int layer, int mod, bool half) const
 
int tileCount (int layer, int ring) const
 
bool tileExist (int zside, int layer, int ring, int phi) const
 
HGCalParameters::tileInfo tileInfo (int zside, int layer, int ring) const
 
bool tilePhiEdge (double phi, int layer, int iphi) const
 
bool tileRingEdge (double rho, int layer, int ring) const
 
std::pair< int, int > tileRings (int layer) const
 
int tileSiPM (int sipm) const
 
bool tileTrapezoid () const
 
std::pair< int, int > tileType (int layer, int ring, int phi) const
 
bool trapezoidFile () const
 
unsigned int volumes () const
 
int waferCount (const int type) const
 
int waferFileIndex (unsigned int kk) const
 
std::tuple< int, int, int > waferFileInfo (unsigned int kk) const
 
bool waferFileInfoExist (int kk) const
 
std::tuple< int, int, int > waferFileInfoFromIndex (int kk) const
 
unsigned int waferFileSize () const
 
int waferFromCopy (int copy) const
 
void waferFromPosition (const double x, const double y, int &wafer, int &icell, int &celltyp) const
 
void waferFromPosition (const double x, const double y, const int zside, const int layer, int &waferU, int &waferV, int &cellU, int &cellV, int &celltype, double &wt, bool extend, bool debug) const
 
bool waferFullInLayer (int wafer, int lay, bool reco) const
 
bool waferHexagon6 () const
 
bool waferHexagon8 () const
 
bool waferHexagon8File () const
 
HGCalParameters::waferInfo waferInfo (int lay, int waferU, int waferV) const
 
bool waferInLayer (int wafer, int lay, bool reco) const
 
GlobalPoint waferLocal2Global (HepGeom::Point3D< float > &loc, const DetId &id, bool useWafer, bool reco, bool debug) const
 
int waferMax () const
 
int waferMin () const
 
std::pair< double, double > waferParameters (bool reco) const
 
std::pair< double, double > waferPosition (int wafer, bool reco) const
 
std::pair< double, double > waferPosition (int lay, int waferU, int waferV, bool reco, bool debug) const
 
int wafers () const
 
int wafers (int layer, int type) const
 
double waferSepar (bool reco) const
 
double waferSize (bool reco) const
 
int waferToCopy (int wafer) const
 
int waferType (DetId const &id, bool fromFile) const
 
int waferType (int layer, int waferU, int waferV, bool fromFile) const
 
std::tuple< int, int, int > waferType (HGCSiliconDetId const &id, bool fromFile) const
 
int waferTypeL (int wafer) const
 
std::pair< int, int > waferTypeRotation (int layer, int waferU, int waferV, bool fromFile, bool debug) const
 
int waferTypeT (int wafer) const
 
int waferUVMax () const
 
bool waferVirtual (int layer, int waferU, int waferV) const
 
double waferZ (int layer, bool reco) const
 
 ~HGCalDDDConstants ()=default
 

Private Types

typedef std::array< int, 3 > HGCWaferParam
 
typedef std::array< std::vector< int32_t >, 2 > Simrecovecs
 

Private Member Functions

int cellHex (double xx, double yy, const double &cellR, const std::vector< double > &posX, const std::vector< double > &posY) const
 
void cellHex (double xloc, double yloc, int cellType, int place, int part, int &cellU, int &cellV, bool extend, bool debug) const
 
std::pair< int, float > getIndex (int lay, bool reco) const
 
bool isValidCell (int layindex, int wafer, int cell) const
 
bool isValidCell8 (int lay, int waferU, int waferV, int cellU, int cellV, int type) const
 
int layerFromIndex (int index, bool reco) const
 
int32_t waferIndex (int wafer, int index) const
 
bool waferInLayerTest (int wafer, int lay, bool full) const
 
std::pair< double, double > waferPosition (int waferU, int waferV, bool reco) const
 
std::pair< double, double > waferPositionNoRot (int lay, int waferU, int waferV, bool reco, bool debug) const
 

Private Attributes

const bool fullAndPart_
 
HGCalGeomTools geomTools_
 
double hexside_
 
double hexsideT_
 
HGCalCassette hgcassette_
 
std::unique_ptr< HGCalCellhgcell_
 
std::unique_ptr< HGCalCellUVhgcellUV_
 
const HGCalParametershgpar_
 
Simrecovecs max_modules_layer_
 
int32_t maxWafersPerLayer_
 
const HGCalGeometryMode::GeometryMode mode_
 
int32_t modHalf_
 
double rmax_
 
double rmaxT_
 
const double sqrt3_
 
std::array< uint32_t, 2 > tot_layers_
 
int32_t tot_wafers_
 
std::unordered_map< int32_t, bool > waferIn_
 
std::map< int, HGCWaferParamwaferLayer_
 
std::array< int, 4 > waferMax_
 

Static Private Attributes

static constexpr float dPhiMin = 0.02
 
static constexpr double k_horizontalShift = 1.0
 
static constexpr double tan30deg_ = 0.5773502693
 
static constexpr double tol_ = 0.001
 

Detailed Description

this class reads the constant section of the numbering xml-files of the high granulairy calorimeter

Date
2014/03/20 00:06:50
Author
Sunanda Banerjee, SINP sunan.nosp@m.da.b.nosp@m.anerj.nosp@m.ee@c.nosp@m.ern.c.nosp@m.h

Definition at line 35 of file HGCalDDDConstants.h.

Member Typedef Documentation

◆ HGCWaferParam

typedef std::array<int, 3> HGCalDDDConstants::HGCWaferParam
private

Definition at line 246 of file HGCalDDDConstants.h.

◆ Simrecovecs

typedef std::array<std::vector<int32_t>, 2> HGCalDDDConstants::Simrecovecs
private

Definition at line 245 of file HGCalDDDConstants.h.

Constructor & Destructor Documentation

◆ HGCalDDDConstants()

HGCalDDDConstants::HGCalDDDConstants ( const HGCalParameters hp,
const std::string &  name 
)

Definition at line 25 of file HGCalDDDConstants.cc.

References HGCalParameters::cassettes_, HGCalParameters::cassetteShift_, HGCalParameters::cellSize_, funct::cos(), fullAndPart_, getIndex(), getTrForm(), getTrFormN(), HGCalGeometryMode::Hexagon8Cassette, hexside_, hexsideT_, hgcassette_, hgcell_, hgcellUV_, hgpar_, mps_fire::i, HGCalParameters::k_ScaleFromDDD, HGCalParameters::hgtrform::lay, pixelTopology::layer, layers(), layersInit(), SiStripPI::max, max_modules_layer_, maxCells(), maxWafersPerLayer_, mode_, modHalf_, modulesInit(), Skims_PA_cff::name, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, rmax_, rmaxT_, sectors(), HGCalParameters::sensorSeparation_, HGCalCassette::setParameter(), tan30deg_, tot_layers_, tot_wafers_, HGCalGeometryMode::TrapezoidCassette, HGCalParameters::waferCopy_, waferHexagon6(), waferHexagon8(), waferIn_, HGCalWaferIndex::waferIndex(), waferInLayer(), waferLayer_, waferMax_, HGCalParameters::waferR_, wafers(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), cms::cuda::wmax, and cms::cuda::wmin.

27 #ifdef EDM_ML_DEBUG
28  edm::LogVerbatim("HGCalGeom") << "Mode " << mode_ << " FullAndPart " << fullAndPart_;
29 #endif
30  if (waferHexagon6() || waferHexagon8()) {
33  hexside_ = 2.0 * rmax_ * tan30deg_;
34  hexsideT_ = 2.0 * rmaxT_ * tan30deg_;
35  hgcell_ = std::make_unique<HGCalCell>(2.0 * rmaxT_, hgpar_->nCellsFine_, hgpar_->nCellsCoarse_);
36  hgcellUV_ = std::make_unique<HGCalCellUV>(
38 #ifdef EDM_ML_DEBUG
39  edm::LogVerbatim("HGCalGeom") << "rmax_ " << rmax_ << ":" << rmaxT_ << ":" << hexside_ << ":" << hexsideT_
40  << " CellSize " << 0.5 * HGCalParameters::k_ScaleFromDDD * hgpar_->cellSize_[0] << ":"
42 #endif
43  }
46 #ifdef EDM_ML_DEBUG
47  edm::LogVerbatim("HGCalGeom") << "Setup HGCalCassette for " << hgpar_->cassettes_ << " cassettes";
48 #endif
49  }
50  // init maps and constants
51  modHalf_ = 0;
53  for (int simreco = 0; simreco < 2; ++simreco) {
54  tot_layers_[simreco] = layersInit((bool)simreco);
55  max_modules_layer_[simreco].resize(tot_layers_[simreco] + 1);
56  for (unsigned int layer = 1; layer <= tot_layers_[simreco]; ++layer) {
57  max_modules_layer_[simreco][layer] = modulesInit(layer, (bool)simreco);
58  if (simreco == 1) {
59  modHalf_ += max_modules_layer_[simreco][layer];
61 #ifdef EDM_ML_DEBUG
62  edm::LogVerbatim("HGCalGeom") << "Layer " << layer << " with " << max_modules_layer_[simreco][layer] << ":"
63  << modHalf_ << " modules in RECO";
64  } else {
65  edm::LogVerbatim("HGCalGeom") << "Layer " << layer << " with " << max_modules_layer_[simreco][layer]
66  << " modules in SIM";
67 #endif
68  }
69  }
70 #ifdef EDM_ML_DEBUG
71  edm::LogVerbatim("HGCalGeom") << "SimReco " << simreco << " with " << tot_layers_[simreco] << " Layers";
72 #endif
73  }
74  tot_wafers_ = wafers();
75 
76 #ifdef EDM_ML_DEBUG
77  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants initialized for " << name << " with " << layers(false) << ":"
78  << layers(true) << " layers, " << wafers() << ":" << 2 * modHalf_
79  << " wafers with maximum " << maxWafersPerLayer_ << " per layer and "
80  << "maximum of " << maxCells(false) << ":" << maxCells(true) << " cells";
81 #endif
82  if (waferHexagon6() || waferHexagon8()) {
83  int wminT(9999999), wmaxT(-9999999), kount1(0), kount2(0);
84  for (unsigned int i = 0; i < getTrFormN(); ++i) {
85  int lay0 = getTrForm(i).lay;
86  int wmin(9999999), wmax(-9999999), kount(0);
87  for (int wafer = 0; wafer < sectors(); ++wafer) {
88  bool waferIn = waferInLayer(wafer, lay0, true);
89  if (waferHexagon8()) {
90  int kndx = HGCalWaferIndex::waferIndex(lay0,
93  waferIn_[kndx] = waferIn;
94  }
95  if (waferIn) {
96  int waferU = ((waferHexagon6()) ? wafer : HGCalWaferIndex::waferU(hgpar_->waferCopy_[wafer]));
97  if (waferU < wmin)
98  wmin = waferU;
99  if (waferU > wmax)
100  wmax = waferU;
101  ++kount;
102  }
103  }
104  if (wminT > wmin)
105  wminT = wmin;
106  if (wmaxT < wmax)
107  wmaxT = wmax;
108  if (kount1 < kount)
109  kount1 = kount;
110  kount2 += kount;
111 #ifdef EDM_ML_DEBUG
112  int lay1 = getIndex(lay0, true).first;
113  edm::LogVerbatim("HGCalGeom") << "Index " << i << " Layer " << lay0 << ":" << lay1 << " Wafer " << wmin << ":"
114  << wmax << ":" << kount;
115 #endif
116  HGCWaferParam a1{{wmin, wmax, kount}};
117  waferLayer_[lay0] = a1;
118  }
119  waferMax_ = std::array<int, 4>{{wminT, wmaxT, kount1, kount2}};
120 #ifdef EDM_ML_DEBUG
121  edm::LogVerbatim("HGCalGeom") << "Overall wafer statistics: " << wminT << ":" << wmaxT << ":" << kount1 << ":"
122  << kount2;
123 #endif
124  }
125 }
Log< level::Info, true > LogVerbatim
unsigned int layersInit(bool reco) const
HGCalCassette hgcassette_
std::unique_ptr< HGCalCellUV > hgcellUV_
Simrecovecs max_modules_layer_
int32_t waferU(const int32_t index)
bool waferHexagon8() const
int modulesInit(int lay, bool reco) const
std::unordered_map< int32_t, bool > waferIn_
bool waferHexagon6() const
HGCalGeometryMode::GeometryMode mode_
bool waferInLayer(int wafer, int lay, bool reco) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::map< int, HGCWaferParam > waferLayer_
std::unique_ptr< HGCalCell > hgcell_
int maxCells(bool reco) const
std::vector< double > cellSize_
T sqrt(T t)
Definition: SSEVec.h:19
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
std::pair< int, float > getIndex(int lay, bool reco) const
unsigned int layers(bool reco) const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
std::array< int, 4 > waferMax_
unsigned int getTrFormN() const
HGCalParameters::hgtrform getTrForm(unsigned int k) const
bool waferHexagon8File() const
const HGCalGeometryMode::GeometryMode mode_
static constexpr double k_ScaleFromDDD
std::array< int, 3 > HGCWaferParam
int32_t waferV(const int32_t index)
std::vector< int > waferCopy_
std::vector< double > cassetteShift_
static constexpr double tan30deg_
std::array< uint32_t, 2 > tot_layers_
__host__ __device__ V wmin
const HGCalParameters * hgpar_
void setParameter(int cassette, const std::vector< double > &shifts)
Definition: HGCalCassette.cc:8
__host__ __device__ V V wmax

◆ ~HGCalDDDConstants()

HGCalDDDConstants::~HGCalDDDConstants ( )
default

Member Function Documentation

◆ assignCell()

std::pair< int, int > HGCalDDDConstants::assignCell ( float  x,
float  y,
int  lay,
int  subSec,
bool  reco 
) const

Definition at line 127 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, cellHex(), HGCalParameters::cellSize_, getIndex(), hgpar_, HGCalParameters::k_ScaleFromDDD, dt_dqm_sourceclient_common_cff::reco, rmax_, findQualityFiles::size, waferHexagon6(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::waferTypeT_, x, geometryCSVtoXML::xx, y, and geometryCSVtoXML::yy.

Referenced by HGCNumberingScheme::assignCell(), and HGCalGeometry::getClosestCell().

127  {
128  const auto& index = getIndex(lay, reco);
129  if (index.first < 0)
130  return std::make_pair(-1, -1);
131  if (waferHexagon6()) {
132  float xx = (reco) ? x : HGCalParameters::k_ScaleFromDDD * x;
133  float yy = (reco) ? y : HGCalParameters::k_ScaleFromDDD * y;
134 
135  // First the wafer
136  int wafer = cellHex(xx, yy, rmax_, hgpar_->waferPosX_, hgpar_->waferPosY_);
137  if (wafer < 0 || wafer >= static_cast<int>(hgpar_->waferTypeT_.size())) {
138  edm::LogWarning("HGCalGeom") << "Wafer no. out of bound for " << wafer << ":" << (hgpar_->waferTypeT_).size()
139  << ":" << (hgpar_->waferPosX_).size() << ":" << (hgpar_->waferPosY_).size()
140  << " ***** ERROR *****";
141  return std::make_pair(-1, -1);
142  } else {
143  // Now the cell
144  xx -= hgpar_->waferPosX_[wafer];
145  yy -= hgpar_->waferPosY_[wafer];
146  if (hgpar_->waferTypeT_[wafer] == 1)
147  return std::make_pair(wafer,
148  cellHex(xx,
149  yy,
152  hgpar_->cellFineY_));
153  else
154  return std::make_pair(wafer,
155  cellHex(xx,
156  yy,
159  hgpar_->cellCoarseY_));
160  }
161  } else {
162  return std::make_pair(-1, -1);
163  }
164 }
size
Write out results.
std::vector< double > waferPosY_
std::vector< double > cellFineY_
bool waferHexagon6() const
std::vector< double > cellCoarseX_
std::vector< double > cellSize_
std::pair< int, float > getIndex(int lay, bool reco) const
std::vector< double > cellFineX_
static constexpr double k_ScaleFromDDD
fixed size matrix
std::vector< int > waferTypeT_
std::vector< double > cellCoarseY_
Log< level::Warning, false > LogWarning
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
int cellHex(double xx, double yy, const double &cellR, const std::vector< double > &posX, const std::vector< double > &posY) const

◆ assignCellHex()

std::array< int, 5 > HGCalDDDConstants::assignCellHex ( float  x,
float  y,
int  zside,
int  lay,
bool  reco,
bool  extend,
bool  debug 
) const

Definition at line 166 of file HGCalDDDConstants.cc.

References debug, HGCalParameters::k_ScaleToDDD, dt_dqm_sourceclient_common_cff::reco, waferFromPosition(), waferHexagon8(), waferType(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, geometryCSVtoXML::xx, y, geometryCSVtoXML::yy, and ecaldqm::zside().

Referenced by HGCalGeometry::getClosestCell(), and HGCalGeometry::getClosestCellHex().

167  {
168  int waferU(0), waferV(0), waferType(-1), cellU(0), cellV(0);
169  if (waferHexagon8()) {
170  double xx = (reco) ? HGCalParameters::k_ScaleToDDD * x : x;
171  double yy = (reco) ? HGCalParameters::k_ScaleToDDD * y : y;
172  double wt(1.0);
173 #ifdef EDM_ML_DEBUG
174  edm::LogVerbatim("HGCalGeom") << "assignCellHex x " << x << ":" << xx << " y " << y << ":" << yy << " Lay " << lay;
175 #endif
176  waferFromPosition(xx, yy, zside, lay, waferU, waferV, cellU, cellV, waferType, wt, extend, debug);
177  }
178  return std::array<int, 5>{{waferU, waferV, waferType, cellU, cellV}};
179 }
Log< level::Info, true > LogVerbatim
void waferFromPosition(const double x, const double y, int &wafer, int &icell, int &celltyp) const
int32_t waferU(const int32_t index)
bool waferHexagon8() const
int zside(DetId const &)
static constexpr double k_ScaleToDDD
#define debug
Definition: HDRShower.cc:19
int waferType(DetId const &id, bool fromFile) const
int32_t waferV(const int32_t index)

◆ assignCellTrap()

std::array< int, 3 > HGCalDDDConstants::assignCellTrap ( float  x,
float  y,
float  z,
int  lay,
bool  reco 
) const

Definition at line 181 of file HGCalDDDConstants.cc.

References funct::abs(), HGCalParameters::cassettes_, angle_units::operators::convertRadToDeg(), HGCalParameters::firstLayer_, getIndex(), HGCalCassette::getShift(), hgcassette_, hgpar_, LEDCalibrationChannels::iphi, HGCalParameters::iradMaxBH_, HGCalParameters::iradMinBH_, HGCalParameters::k_ScaleFromDDD, pixelTopology::layer, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, M_PI, mode_, HGCalParameters::nphiCassette_, phi, HGCalParameters::phiOffset_, alignCSCRings::r, HGCalParameters::radiusLayer_, dt_dqm_sourceclient_common_cff::reco, HGCalParameters::scintCells(), HGCalParameters::scintCellSize(), HGCalParameters::scintType(), mathSSE::sqrt(), HGCalTileIndex::tileCassette(), tileExist(), tilePhiEdge(), tileRingEdge(), tol_, HGCalGeometryMode::TrapezoidCassette, x, HGCalParameters::xLayerHex_, geometryCSVtoXML::xx, y, HGCalParameters::yLayerHex_, geometryCSVtoXML::yy, z, and ecaldqm::zside().

Referenced by HGCalGeometry::getClosestCell(), and HGCalNumberingScheme::getUnitID().

181  {
182  int irad(-1), iphi(-1), type(-1);
183  const auto& indx = getIndex(layer, reco);
184  if (indx.first < 0)
185  return std::array<int, 3>{{irad, iphi, type}};
186  int zside = (z > 0) ? 1 : -1;
187  double xx = (reco) ? (zside * x) : (zside * HGCalParameters::k_ScaleFromDDD * x);
188  double yy = (reco) ? y : HGCalParameters::k_ScaleFromDDD * y;
189  int ll = layer - hgpar_->firstLayer_;
190  xx -= hgpar_->xLayerHex_[ll];
191  yy -= hgpar_->yLayerHex_[ll];
194  auto cshift = hgcassette_.getShift(layer, zside, cassette);
195 #ifdef EDM_ML_DEBUG
196  edm::LogVerbatim("HGCalGeomT") << "Cassette " << cassette << " Shift " << -(zside * cshift.first) << ":"
197  << cshift.second;
198 #endif
199  xx += (zside * cshift.first);
200  yy -= cshift.second;
201  }
202  double phi = (((yy == 0.0) && (xx == 0.0)) ? 0. : std::atan2(yy, xx));
203  if (phi < 0)
204  phi += (2.0 * M_PI);
205  if (indx.second != 0)
206  iphi = (1 + static_cast<int>(phi / indx.second)) % hgpar_->scintCells(layer);
207  if (iphi == 0)
210  double r = std::sqrt(xx * xx + yy * yy);
211  auto ir = std::lower_bound(hgpar_->radiusLayer_[type].begin(), hgpar_->radiusLayer_[type].end(), r);
212  irad = static_cast<int>(ir - hgpar_->radiusLayer_[type].begin());
213  irad = std::clamp(irad, hgpar_->iradMinBH_[indx.first], hgpar_->iradMaxBH_[indx.first]);
214 #ifdef EDM_ML_DEBUG
215  edm::LogVerbatim("HGCalGeomT") << "assignCellTrap Input " << x << ":" << y << ":" << z << ":" << layer << ":" << reco
216  << " x|y|r " << xx << ":" << yy << ":" << r << " phi " << phi << ":"
217  << convertRadToDeg(phi) << " o/p " << irad << ":" << iphi << ":" << type;
218 #endif
219  if (!tileExist(zside, layer, irad, iphi)) {
220  if (tileRingEdge(r, layer, irad)) {
221  if (std::abs(r - hgpar_->radiusLayer_[type][irad - 1]) < tol_) {
222  --irad;
223  if (irad <= hgpar_->iradMinBH_[indx.first])
224  irad = hgpar_->iradMinBH_[indx.first];
225  } else {
226  ++irad;
227  if (irad > hgpar_->iradMaxBH_[indx.first])
228  irad = hgpar_->iradMaxBH_[indx.first];
229  }
230 #ifdef EDM_ML_DEBUG
231  edm::LogVerbatim("HGCalGeomT") << "assignCellTrap: ring # modified to " << irad << ":"
232  << hgpar_->iradMinBH_[indx.first] << ":" << hgpar_->iradMaxBH_[indx.first];
233  ;
234 #endif
235  } else if (tilePhiEdge(phi, layer, iphi)) {
236  if (std::abs(phi - hgpar_->scintCellSize(layer) * (iphi - 1)) < tol_) {
237  --iphi;
238  if (iphi <= 0)
239  iphi = 1;
240  } else {
241  ++iphi;
242  if (iphi > hgpar_->scintCells(layer))
243  iphi = 1;
244  }
245 #ifdef EDM_ML_DEBUG
246  edm::LogVerbatim("HGCalGeomT") << "assignCellTrap: iphi # modified to " << iphi << ":"
247  << hgpar_->scintCells(layer);
248 #endif
249  }
250  }
251  return std::array<int, 3>{{irad, iphi, type}};
252 }
std::vector< int > iradMaxBH_
Log< level::Info, true > LogVerbatim
static constexpr double tol_
HGCalCassette hgcassette_
int scintType(const int layer) const
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
bool tileRingEdge(double rho, int layer, int ring) const
std::pair< double, double > getShift(int layer, int zside, int cassette) const
int zside(DetId const &)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
double scintCellSize(const int layer) const
std::vector< double > yLayerHex_
T sqrt(T t)
Definition: SSEVec.h:19
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::pair< int, float > getIndex(int lay, bool reco) const
bool tilePhiEdge(double phi, int layer, int iphi) const
int32_t tileCassette(int32_t, int32_t, int32_t, int32_t)
#define M_PI
bool tileExist(int zside, int layer, int ring, int phi) const
std::vector< int > iradMinBH_
const HGCalGeometryMode::GeometryMode mode_
static constexpr double k_ScaleFromDDD
fixed size matrix
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_
int scintCells(const int layer) const
std::vector< double > xLayerHex_

◆ cassetteShiftScintillator()

bool HGCalDDDConstants::cassetteShiftScintillator ( int  layer,
int  iphi 
) const

Definition at line 268 of file HGCalDDDConstants.cc.

References cassetteTile(), HGCalCassette::getShift(), hgcassette_, LEDCalibrationChannels::iphi, pixelTopology::layer, mode_, edm::shift, and HGCalGeometryMode::TrapezoidCassette.

Referenced by HGCScintSD::setDetUnitId().

268  {
269  bool shift(false);
271  auto cshift = hgcassette_.getShift(layer, 1, cassetteTile(iphi));
272  if ((cshift.first != 0) || (cshift.second != 0))
273  shift = true;
274  }
275  return shift;
276 }
HGCalCassette hgcassette_
std::pair< double, double > getShift(int layer, int zside, int cassette) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
const HGCalGeometryMode::GeometryMode mode_
static unsigned int const shift
int cassetteTile(int iphi) const

◆ cassetteShiftSilicon()

bool HGCalDDDConstants::cassetteShiftSilicon ( int  layer,
int  waferU,
int  waferV 
) const

Definition at line 254 of file HGCalDDDConstants.cc.

References HGCalCassette::getShift(), HGCalGeometryMode::Hexagon8Cassette, hgcassette_, hgpar_, pixelTopology::layer, mode_, edm::shift, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

Referenced by HGCalSD::setDetUnitId().

254  {
255  bool shift(false);
258  auto ktr = hgpar_->waferInfoMap_.find(indx);
259  if (ktr != hgpar_->waferInfoMap_.end()) {
260  auto cshift = hgcassette_.getShift(layer, 1, (ktr->second).cassette);
261  if ((cshift.first != 0) || (cshift.second != 0))
262  shift = true;
263  }
264  }
265  return shift;
266 }
HGCalCassette hgcassette_
int32_t waferU(const int32_t index)
std::pair< double, double > getShift(int layer, int zside, int cassette) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
const HGCalGeometryMode::GeometryMode mode_
int32_t waferV(const int32_t index)
static unsigned int const shift
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ cassetteTile()

int HGCalDDDConstants::cassetteTile ( int  iphi) const
inline

◆ cellEtaPhiTrap()

std::pair< double, double > HGCalDDDConstants::cellEtaPhiTrap ( int  type,
int  irad 
) const

Definition at line 278 of file HGCalDDDConstants.cc.

References HGCalParameters::cellSize_, hgcalPerformanceValidation::df, l1ctLayer1_cff::dr, hgpar_, alignCSCRings::r, HGCalParameters::radiusLayer_, and tileTrapezoid().

278  {
279  double dr(0), df(0);
280  if (tileTrapezoid()) {
281  double r = 0.5 * ((hgpar_->radiusLayer_[type][irad - 1] + hgpar_->radiusLayer_[type][irad]));
282  dr = (hgpar_->radiusLayer_[type][irad] - hgpar_->radiusLayer_[type][irad - 1]);
283  df = r * hgpar_->cellSize_[type];
284  }
285  return std::make_pair(dr, df);
286 }
std::vector< double > cellSize_
bool tileTrapezoid() const
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_

◆ cellHex() [1/2]

int HGCalDDDConstants::cellHex ( double  xx,
double  yy,
const double &  cellR,
const std::vector< double > &  posX,
const std::vector< double > &  posY 
) const
private

Definition at line 1825 of file HGCalDDDConstants.cc.

References funct::abs(), PVValHelper::dx, PVValHelper::dy, dqmdumpme::k, EgammaValidation_cff::num, RecoTauValidation_cfi::posX, RecoTauValidation_cfi::posY, tan30deg_, TrackerOfflineValidation_Dqm_cff::xmax, geometryCSVtoXML::xx, and geometryCSVtoXML::yy.

Referenced by assignCell(), and waferFromPosition().

1826  {
1827  int num(0);
1828  const double tol(0.00001);
1829  double cellY = 2.0 * cellR * tan30deg_;
1830  for (unsigned int k = 0; k < posX.size(); ++k) {
1831  double dx = std::abs(xx - posX[k]);
1832  double dy = std::abs(yy - posY[k]);
1833  if (dx <= (cellR + tol) && dy <= (cellY + tol)) {
1834  double xmax = (dy <= 0.5 * cellY) ? cellR : (cellR - (dy - 0.5 * cellY) / tan30deg_);
1835  if (dx <= (xmax + tol)) {
1836  num = k;
1837  break;
1838  }
1839  }
1840  }
1841  return num;
1842 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
static constexpr double tan30deg_

◆ cellHex() [2/2]

void HGCalDDDConstants::cellHex ( double  xloc,
double  yloc,
int  cellType,
int  place,
int  part,
int &  cellU,
int &  cellV,
bool  extend,
bool  debug 
) const
private

Definition at line 1844 of file HGCalDDDConstants.cc.

References funct::abs(), cellType(), debug, newFWLiteAna::found, HGCalGeometryMode::Hexagon8Cassette, hgcellUV_, hgpar_, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, SiStripPI::max, SiStripPI::min, mode_, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, rmax_, rmaxT_, edm::shift, sqrt3_, V0Monitor_cff::v0, HGCalTypes::WaferFull, and waferHexagon8File().

1845  {
1847  auto uv = (part == HGCalTypes::WaferFull)
1848  ? hgcellUV_->cellUVFromXY3(xloc, yloc, place, cellType, true, debug)
1849  : hgcellUV_->cellUVFromXY1(xloc, yloc, place, cellType, part, true, debug);
1850  cellU = uv.first;
1851  cellV = uv.second;
1852  } else if (waferHexagon8File()) {
1853  auto uv = hgcellUV_->cellUVFromXY3(xloc, yloc, place, cellType, extend, debug);
1854  cellU = uv.first;
1855  cellV = uv.second;
1856  } else {
1857  int ncell = (cellType == 0) ? hgpar_->nCellsFine_ : hgpar_->nCellsCoarse_;
1858  double delY = 2 * rmax_ / (3 * ncell);
1859  double delX = 0.5 * delY * sqrt3_;
1860  double delYT = (extend) ? (2 * rmaxT_ / (3 * ncell)) : delY;
1861  double delXT = 0.5 * delYT * sqrt3_;
1862  double v0 = ((xloc / delY - 1.0) / 1.5);
1863  int cv0 = (v0 > 0) ? (ncell + static_cast<int>(v0 + 0.5)) : (ncell - static_cast<int>(-v0 + 0.5));
1864  double u0 = (0.5 * yloc / delX + 0.5 * cv0);
1865  int cu0 = (u0 > 0) ? (ncell / 2 + static_cast<int>(u0 + 0.5)) : (ncell / 2 - static_cast<int>(-u0 + 0.5));
1866  cu0 = std::max(0, std::min(cu0, 2 * ncell - 1));
1867  cv0 = std::max(0, std::min(cv0, 2 * ncell - 1));
1868  if (cv0 - cu0 >= ncell)
1869  cv0 = cu0 + ncell - 1;
1870  if (debug)
1871  edm::LogVerbatim("HGCalGeom") << "cellHex: input " << xloc << ":" << yloc << ":" << cellType << " parameter "
1872  << delX << ":" << delY << " u0 " << u0 << ":" << cu0 << " v0 " << v0 << ":" << cv0;
1873  bool found(false);
1874  static constexpr int shift[3] = {0, 1, -1};
1875  for (int i1 = 0; i1 < 3; ++i1) {
1876  cellU = cu0 + shift[i1];
1877  for (int i2 = 0; i2 < 3; ++i2) {
1878  cellV = cv0 + shift[i2];
1879  if (((cellV - cellU) < ncell) && ((cellU - cellV) <= ncell) && (cellU >= 0) && (cellV >= 0) &&
1880  (cellU < 2 * ncell) && (cellV < 2 * ncell)) {
1881  double xc = (1.5 * (cellV - ncell) + 1.0) * delY;
1882  double yc = (2 * cellU - cellV - ncell) * delX;
1883  if ((std::abs(yloc - yc) <= delXT) && (std::abs(xloc - xc) <= delYT) &&
1884  ((std::abs(xloc - xc) <= 0.5 * delYT) ||
1885  (std::abs(yloc - yc) <= sqrt3_ * (delYT - std::abs(xloc - xc))))) {
1886  if (debug)
1887  edm::LogVerbatim("HGCalGeom")
1888  << "cellHex: local " << xc << ":" << yc << " difference " << std::abs(xloc - xc) << ":"
1889  << std::abs(yloc - yc) << ":" << sqrt3_ * (delY - std::abs(yloc - yc)) << " comparator " << delX
1890  << ":" << delY << " (u,v) = (" << cellU << "," << cellV << ")";
1891  found = true;
1892  break;
1893  }
1894  }
1895  }
1896  if (found)
1897  break;
1898  }
1899  if (!found) {
1900  cellU = cu0;
1901  cellV = cv0;
1902  }
1903  }
1904 }
Log< level::Info, true > LogVerbatim
std::unique_ptr< HGCalCellUV > hgcellUV_
int32_t cellType(int type, int waferU, int waferV, int iz, int fwdBack, int orient) const
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
#define debug
Definition: HDRShower.cc:19
part
Definition: HCALResponse.h:20
bool waferHexagon8File() const
const HGCalGeometryMode::GeometryMode mode_
static unsigned int const shift
const HGCalParameters * hgpar_

◆ cellInLayer()

bool HGCalDDDConstants::cellInLayer ( int  waferU,
int  waferV,
int  cellU,
int  cellV,
int  lay,
int  zside,
bool  reco 
) const

Definition at line 288 of file HGCalDDDConstants.cc.

References getIndex(), HGCalWaferMask::goodCell(), HGCalGeometryMode::Hexagon8Cassette, HGCalGeometryMode::Hexagon8Module, hgpar_, locateCell(), mode_, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, hgcalPerformanceValidation::orient, HGCalParameters::rMinLayHex_, mathSSE::sqrt(), HGCalTypes::WaferFineThin, HGCalTypes::WaferFull, waferHexagon6(), waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalTypes::WaferOrient0, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), geometryCSVtoXML::xy, and ecaldqm::zside().

Referenced by HGCalGeometry::newCell().

288  {
289  const auto& indx = getIndex(lay, true);
290  if (indx.first >= 0) {
292  int indx = HGCalWaferIndex::waferIndex(lay, waferU, waferV);
293  auto ktr = hgpar_->waferInfoMap_.find(indx);
294  int part = (ktr != hgpar_->waferInfoMap_.end()) ? (ktr->second).part : HGCalTypes::WaferFull;
295  return HGCalWaferMask::goodCell(cellU, cellV, part);
297  int indx = HGCalWaferIndex::waferIndex(lay, waferU, waferV);
298  auto ktr = hgpar_->waferInfoMap_.find(indx);
300  if (ktr != hgpar_->waferInfoMap_.end()) {
301  thck = (ktr->second).type;
302  part = (ktr->second).part;
303  rotn = (ktr->second).orient;
304  }
306  return HGCalWaferMask::goodCell(cellU, cellV, ncell, part, rotn);
307  } else if (waferHexagon8() || waferHexagon6()) {
308  const auto& xy =
309  ((waferHexagon8()) ? locateCell(zside, lay, waferU, waferV, cellU, cellV, reco, true, false, false)
310  : locateCell(cellU, lay, waferU, reco));
311  double rpos = sqrt(xy.first * xy.first + xy.second * xy.second);
312  return ((rpos >= hgpar_->rMinLayHex_[indx.first]) && (rpos <= hgpar_->rMaxLayHex_[indx.first]));
313  } else {
314  return true;
315  }
316  } else {
317  return false;
318  }
319 }
int32_t waferU(const int32_t index)
static bool goodCell(int u, int v, int N, int type, int rotn)
bool waferHexagon8() const
int zside(DetId const &)
bool waferHexagon6() const
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
T sqrt(T t)
Definition: SSEVec.h:19
static constexpr int32_t WaferOrient0
Definition: HGCalTypes.h:61
std::pair< int, float > getIndex(int lay, bool reco) const
std::vector< double > rMinLayHex_
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
part
Definition: HCALResponse.h:20
std::pair< float, float > locateCell(int cell, int lay, int type, bool reco) const
const HGCalGeometryMode::GeometryMode mode_
fixed size matrix
int32_t waferV(const int32_t index)
static constexpr int32_t WaferFineThin
Definition: HGCalTypes.h:30
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ cellSizeHex()

double HGCalDDDConstants::cellSizeHex ( int  type) const

Definition at line 334 of file HGCalDDDConstants.cc.

References HGCalParameters::cellSize_, hgpar_, HGCalParameters::k_ScaleFromDDD, tileTrapezoid(), and waferHexagon8().

Referenced by HGCalGeometryLoader::build(), and hgcal::RecHitTools::getRadiusToSide().

334  {
335  int indx = ((waferHexagon8()) ? ((type >= 1) ? 1 : 0) : ((type == 1) ? 1 : 0));
336  double cell = (tileTrapezoid() ? 0.5 * hgpar_->cellSize_[indx]
338  return cell;
339 }
bool waferHexagon8() const
std::vector< double > cellSize_
bool tileTrapezoid() const
static constexpr double k_ScaleFromDDD
const HGCalParameters * hgpar_

◆ cellSizeTrap()

std::pair<double, double> HGCalDDDConstants::cellSizeTrap ( int  type,
int  irad 
) const
inline

Definition at line 51 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::radiusLayer_.

Referenced by HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), and HGCalGeometry::getNewCorners().

51  {
52  return std::make_pair(hgpar_->radiusLayer_[type][irad - 1], hgpar_->radiusLayer_[type][irad]);
53  }
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_

◆ cellThickness()

double HGCalDDDConstants::cellThickness ( int  layer,
int  waferU,
int  waferV 
) const

Definition at line 321 of file HGCalDDDConstants.cc.

References HGCalParameters::cellThickness_, hgpar_, pixelTopology::layer, waferHexagon6(), waferHexagon8(), waferType(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

Referenced by HGCalTBAnalyzer::analyzeSimHits(), and hgcal::RecHitTools::getSiThickness().

321  {
322  double thick(-1);
323  int type = waferType(layer, waferU, waferV, false);
324  if (type >= 0) {
325  if (waferHexagon8()) {
326  thick = 10000.0 * hgpar_->cellThickness_[type]; // cm to micron
327  } else if (waferHexagon6()) {
328  thick = 100.0 * (type + 1); // type = 1,2,3 for 100,200,300 micron
329  }
330  }
331  return thick;
332 }
int32_t waferU(const int32_t index)
bool waferHexagon8() const
bool waferHexagon6() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::vector< double > cellThickness_
int waferType(DetId const &id, bool fromFile) const
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_

◆ cellType()

int32_t HGCalDDDConstants::cellType ( int  type,
int  waferU,
int  waferV,
int  iz,
int  fwdBack,
int  orient 
) const

Definition at line 341 of file HGCalDDDConstants.cc.

References HGCalCell::cellPlacementIndex(), HGCalCell::cellPlacementOld, HGCalCell::cellType(), hgpar_, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, and hgcalPerformanceValidation::orient.

Referenced by cellHex(), and HGCalTopology::neighbors().

341  {
342  int placement = (orient < 0) ? HGCalCell::cellPlacementOld : HGCalCell::cellPlacementIndex(iz, fwdBack, orient);
343  int ncell = (type == 0) ? hgpar_->nCellsFine_ : hgpar_->nCellsCoarse_;
344  auto cellType = HGCalCell::cellType(cellU, cellV, ncell, placement);
345  return cellType.first;
346 }
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:233
static constexpr int32_t cellPlacementOld
Definition: HGCalCell.h:25
static std::pair< int32_t, int32_t > cellType(int32_t u, int32_t v, int32_t ncell, int32_t placementIndex)
Definition: HGCalCell.cc:245
int32_t cellType(int type, int waferU, int waferV, int iz, int fwdBack, int orient) const
const HGCalParameters * hgpar_

◆ distFromEdgeHex()

double HGCalDDDConstants::distFromEdgeHex ( double  x,
double  y,
double  z 
) const

Definition at line 348 of file HGCalDDDConstants.cc.

References funct::abs(), PVValHelper::dx, PVValHelper::dy, getLayer(), hexside_, hgpar_, dqmdumpme::k, HGCalParameters::k_ScaleFromDDD, HGCalParameters::k_ScaleToDDD, layerIndex(), rmax_, sqrt3_, tan30deg_, waferHexagon8(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, x, HGCalParameters::xLayerHex_, geometryCSVtoXML::xx, y, HGCalParameters::yLayerHex_, geometryCSVtoXML::yy, and z.

Referenced by HFNoseSD::isItinFidVolume(), and HGCalSD::isItinFidVolume().

348  {
349  // Assming the point is within a hexagonal plane of the wafer, calculate
350  // the shortest distance from the edge
351  if (z < 0)
352  x = -x;
353  double dist(0);
354  // Input x, y in Geant4 unit and transformed to CMSSW standard
357  if (waferHexagon8()) {
358  int ll = layerIndex(getLayer(z, false), false);
359  xx -= hgpar_->xLayerHex_[ll];
360  yy -= hgpar_->yLayerHex_[ll];
361  }
362  int sizew = static_cast<int>(hgpar_->waferPosX_.size());
363  int wafer = sizew;
364  // Transform to the local coordinate frame of the wafer first
365  for (int k = 0; k < sizew; ++k) {
366  double dx = std::abs(xx - hgpar_->waferPosX_[k]);
367  double dy = std::abs(yy - hgpar_->waferPosY_[k]);
368  if ((dx <= rmax_) && (dy <= hexside_) && ((dy <= 0.5 * hexside_) || (dx * tan30deg_ <= (hexside_ - dy)))) {
369  wafer = k;
370  xx -= hgpar_->waferPosX_[k];
371  yy -= hgpar_->waferPosY_[k];
372  break;
373  }
374  }
375  // Look at only one quarter (both x,y are positive)
376  if (wafer < sizew) {
377  if (std::abs(yy) < 0.5 * hexside_) {
378  dist = rmax_ - std::abs(xx);
379  } else {
380  dist = 0.5 * ((rmax_ - std::abs(xx)) - sqrt3_ * (std::abs(yy) - 0.5 * hexside_));
381  }
382  } else {
383  dist = 0;
384  }
386 #ifdef EDM_ML_DEBUG
387  edm::LogVerbatim("HGCalGeom") << "DistFromEdgeHex: Local " << xx << ":" << yy << " wafer " << wafer << " flag "
388  << (wafer < sizew) << " Distance " << rmax_ << ":" << (rmax_ - std::abs(xx)) << ":"
389  << (std::abs(yy) - 0.5 * hexside_) << ":" << 0.5 * hexside_ << ":" << dist;
390 #endif
391  return dist;
392 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
bool waferHexagon8() const
std::vector< double > yLayerHex_
int layerIndex(int lay, bool reco) const
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
static constexpr double k_ScaleToDDD
static constexpr double k_ScaleFromDDD
int getLayer(double z, bool reco) const
static constexpr double tan30deg_
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
std::vector< double > xLayerHex_

◆ distFromEdgeTrap()

double HGCalDDDConstants::distFromEdgeTrap ( double  x,
double  y,
double  z 
) const

Definition at line 394 of file HGCalDDDConstants.cc.

References funct::abs(), getLayer(), hgpar_, LEDCalibrationChannels::iphi, HGCalParameters::iradMaxBH_, HGCalParameters::iradMinBH_, HGCalParameters::k_ScaleFromDDD, HGCalParameters::k_ScaleToDDD, layerIndex(), pfDeepBoostedJetPreprocessParams_cfi::lower_bound, M_PI, SiStripPI::max, SiStripPI::min, phi, alignCSCRings::r, HGCalParameters::radiusLayer_, HGCalParameters::scintCellSize(), HGCalParameters::scintType(), mathSSE::sqrt(), x, geometryCSVtoXML::xx, y, and z.

Referenced by HGCScintSD::isItinFidVolume().

394  {
395  // Assming the point is within the eta-phi plane of the scintillator tile,
396  // calculate the shortest distance from the edge
397  int lay = getLayer(z, false);
398  double xx = (z < 0) ? -x : x;
399  int indx = layerIndex(lay, false);
400  double r = HGCalParameters::k_ScaleFromDDD * std::sqrt(x * x + y * y);
401  double phi = (r == 0. ? 0. : std::atan2(y, xx));
402  if (phi < 0)
403  phi += (2.0 * M_PI);
404  int type = hgpar_->scintType(lay);
405  double cell = hgpar_->scintCellSize(lay);
406  // Compare with the center of the tile find distances along R and also phi
407  // Take the smaller value
408  auto ir = std::lower_bound(hgpar_->radiusLayer_[type].begin(), hgpar_->radiusLayer_[type].end(), r);
409  int irad = static_cast<int>(ir - hgpar_->radiusLayer_[type].begin());
410  irad = std::clamp(irad, hgpar_->iradMinBH_[indx], hgpar_->iradMaxBH_[indx]);
411  int iphi = 1 + static_cast<int>(phi / cell);
412  double dphi = std::max(0.0, (0.5 * cell - std::abs(phi - (iphi - 0.5) * cell)));
413  double dist = std::min((r - hgpar_->radiusLayer_[type][irad - 1]), (hgpar_->radiusLayer_[type][irad] - r));
414 #ifdef EDM_ML_DEBUG
415  edm::LogVerbatim("HGCalGeom") << "DistFromEdgeTrap: Global " << x << ":" << y << ":" << z << " Layer " << lay
416  << " Index " << indx << ":" << type << " xx " << xx << " R " << r << ":" << irad << ":"
417  << hgpar_->radiusLayer_[type][irad - 1] << ":" << hgpar_->radiusLayer_[type][irad]
418  << " Phi " << phi << ":" << iphi << ":" << (iphi - 0.5) * cell << " cell " << cell
419  << " Dphi " << dphi << " Dist " << dist << ":" << r * dphi;
420 #endif
421  return HGCalParameters::k_ScaleToDDD * std::min(r * dphi, dist);
422 }
std::vector< int > iradMaxBH_
Log< level::Info, true > LogVerbatim
int scintType(const int layer) const
double scintCellSize(const int layer) const
T sqrt(T t)
Definition: SSEVec.h:19
int layerIndex(int lay, bool reco) const
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
#define M_PI
static constexpr double k_ScaleToDDD
std::vector< int > iradMinBH_
static constexpr double k_ScaleFromDDD
int getLayer(double z, bool reco) const
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_

◆ etaPhiFromPosition()

void HGCalDDDConstants::etaPhiFromPosition ( const double  x,
const double  y,
const double  z,
const int  layer,
int &  ieta,
int &  iphi,
int &  type,
double &  wt 
) const

◆ firstLayer()

int HGCalDDDConstants::firstLayer ( ) const
inline

◆ geomMode()

HGCalGeometryMode::GeometryMode HGCalDDDConstants::geomMode ( ) const
inline

◆ getIndex()

std::pair< int, float > HGCalDDDConstants::getIndex ( int  lay,
bool  reco 
) const
private

Definition at line 1906 of file HGCalDDDConstants.cc.

References hgpar_, layerIndex(), HGCalParameters::moduleCellR_, HGCalParameters::moduleCellS_, HGCalParameters::scintCellSize(), waferHexagon6(), and waferHexagon8().

Referenced by assignCell(), assignCellTrap(), cellInLayer(), HGCalDDDConstants(), isValidTrap(), locateCell(), locateCellTrap(), maxCells(), maxRows(), modules(), modulesInit(), numberCells(), rangeRLayer(), simToReco(), waferFullInLayer(), waferInLayer(), wafers(), and waferZ().

1906  {
1907  int indx = layerIndex(lay, reco);
1908  if (indx < 0)
1909  return std::make_pair(-1, 0);
1910  float cell(0);
1911  if (waferHexagon6()) {
1912  cell = (reco ? hgpar_->moduleCellR_[0] : hgpar_->moduleCellS_[0]);
1913  } else {
1914  if (waferHexagon8()) {
1915  cell = (reco ? hgpar_->moduleCellR_[0] : hgpar_->moduleCellS_[0]);
1916  } else {
1917  cell = hgpar_->scintCellSize(lay);
1918  }
1919  }
1920  return std::make_pair(indx, cell);
1921 }
std::vector< double > moduleCellR_
bool waferHexagon8() const
bool waferHexagon6() const
double scintCellSize(const int layer) const
int layerIndex(int lay, bool reco) const
std::vector< double > moduleCellS_
fixed size matrix
const HGCalParameters * hgpar_

◆ getLayer()

int HGCalDDDConstants::getLayer ( double  z,
bool  reco 
) const

Definition at line 424 of file HGCalDDDConstants.cc.

References funct::abs(), HGCalParameters::firstLayer_, hgpar_, dqmdumpme::k, HGCalParameters::k_ScaleFromDDD, HGCalParameters::layerGroupO_, layerIndex(), dt_dqm_sourceclient_common_cff::reco, waferHexagon6(), z, HGCalParameters::zLayerHex_, and geometryCSVtoXML::zz.

Referenced by distFromEdgeHex(), distFromEdgeTrap(), HGCMouseBite::exclude(), HGCalGeometry::getClosestCell(), and HGCalGeometry::getClosestCellHex().

424  {
425  // Get the layer # from the gloabl z coordinate
426  unsigned int k = 0;
427  double zz = (reco ? std::abs(z) : HGCalParameters::k_ScaleFromDDD * std::abs(z));
428  const auto& zLayerHex = hgpar_->zLayerHex_;
429  auto itr = std::find_if(zLayerHex.begin() + 1, zLayerHex.end(), [&k, &zz, &zLayerHex](double zLayer) {
430  ++k;
431  return zz < 0.5 * (zLayerHex[k - 1] + zLayerHex[k]);
432  });
433  int lay = (itr == zLayerHex.end()) ? static_cast<int>(zLayerHex.size()) : k;
434  if (waferHexagon6() && reco) {
435  int indx = layerIndex(lay, false);
436  if (indx >= 0)
437  lay = hgpar_->layerGroupO_[indx];
438  } else {
439  lay += (hgpar_->firstLayer_ - 1);
440  }
441  return lay;
442 }
bool waferHexagon6() const
int layerIndex(int lay, bool reco) const
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< double > zLayerHex_
std::vector< int > layerGroupO_
fixed size matrix
const HGCalParameters * hgpar_

◆ getLayerOffset()

int HGCalDDDConstants::getLayerOffset ( ) const
inline

Definition at line 69 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::layerOffset_.

Referenced by HGCalNumberingScheme::HGCalNumberingScheme().

69 { return hgpar_->layerOffset_; }
const HGCalParameters * hgpar_

◆ getModule()

HGCalParameters::hgtrap HGCalDDDConstants::getModule ( unsigned int  k,
bool  hexType,
bool  reco 
) const

Definition at line 444 of file HGCalDDDConstants.cc.

References HGCalParameters::getModule(), HGCSiliconDetId::HGCalCoarseThick, hgpar_, findQualityFiles::size, HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, and HGCalParameters::waferTypeL_.

Referenced by HGCalGeometryLoader::build(), and HGCalGeometry::getSummary().

444  {
446  if (hexType) {
447  if (indx >= hgpar_->waferTypeL_.size())
448  edm::LogWarning("HGCalGeom") << "Wafer no. out bound for index " << indx << ":" << (hgpar_->waferTypeL_).size()
449  << ":" << (hgpar_->waferPosX_).size() << ":" << (hgpar_->waferPosY_).size()
450  << " ***** ERROR *****";
451  unsigned int type =
452  ((indx < hgpar_->waferTypeL_.size()) ? hgpar_->waferTypeL_[indx] - 1 : HGCSiliconDetId::HGCalCoarseThick);
453  mytr = hgpar_->getModule(type, reco);
454  } else {
455  mytr = hgpar_->getModule(indx, reco);
456  }
457  return mytr;
458 }
size
Write out results.
std::vector< double > waferPosY_
hgtrap getModule(unsigned int k, bool reco) const
fixed size matrix
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
std::vector< int > waferTypeL_

◆ getModules()

std::vector< HGCalParameters::hgtrap > HGCalDDDConstants::getModules ( ) const

Definition at line 460 of file HGCalDDDConstants.cc.

References HGCalParameters::getModule(), hgpar_, dqmdumpme::k, and HGCalParameters::moduleLayR_.

460  {
461  std::vector<HGCalParameters::hgtrap> mytrs;
462  for (unsigned int k = 0; k < hgpar_->moduleLayR_.size(); ++k)
463  mytrs.emplace_back(hgpar_->getModule(k, true));
464  return mytrs;
465 }
hgtrap getModule(unsigned int k, bool reco) const
std::vector< int > moduleLayR_
const HGCalParameters * hgpar_

◆ getParameter()

const HGCalParameters* HGCalDDDConstants::getParameter ( ) const
inline

◆ getPhiBins()

int HGCalDDDConstants::getPhiBins ( int  lay) const

Definition at line 467 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::scintCells(), and tileTrapezoid().

467 { return (tileTrapezoid() ? hgpar_->scintCells(lay) : 0); }
bool tileTrapezoid() const
const HGCalParameters * hgpar_
int scintCells(const int layer) const

◆ getRadiusLayer()

const std::vector<double>& HGCalDDDConstants::getRadiusLayer ( int  layer) const
inline

Definition at line 76 of file HGCalDDDConstants.h.

References hgpar_, pixelTopology::layer, HGCalParameters::radiusLayer_, HGCalParameters::scintType(), and tileTrapezoid().

76  {
78  }
int scintType(const int layer) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
bool tileTrapezoid() const
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_

◆ getRangeR()

std::pair< double, double > HGCalDDDConstants::getRangeR ( int  lay,
bool  reco 
) const

Definition at line 469 of file HGCalDDDConstants.cc.

References hgpar_, layerIndex(), HGCalParameters::rMaxLayHex_, and HGCalParameters::rMinLayHex_.

469  {
470  int indx = layerIndex(lay, false);
471  if ((indx >= 0) && (indx < static_cast<int>(hgpar_->rMinLayHex_.size())))
472  return std::make_pair(hgpar_->rMinLayHex_[indx], hgpar_->rMaxLayHex_[indx]);
473  else
474  return std::make_pair(0, -1.);
475 }
int layerIndex(int lay, bool reco) const
std::vector< double > rMinLayHex_
std::vector< double > rMaxLayHex_
const HGCalParameters * hgpar_

◆ getREtaRange()

std::pair< int, int > HGCalDDDConstants::getREtaRange ( int  lay) const

Definition at line 477 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::iradMaxBH_, HGCalParameters::iradMinBH_, layerIndex(), and tileTrapezoid().

Referenced by HGCalTriggerGeometryV9Imp2::getModuleFromTriggerCell(), and HGCalTriggerGeometryV9Imp2::getTriggerCellsFromModule().

477  {
478  int irmin(0), irmax(0);
479  if (tileTrapezoid()) {
480  int indx = layerIndex(lay, false);
481  if ((indx >= 0) && (indx < static_cast<int>(hgpar_->iradMinBH_.size()))) {
482  irmin = hgpar_->iradMinBH_[indx];
483  irmax = hgpar_->iradMaxBH_[indx];
484  }
485  }
486  return std::make_pair(irmin, irmax);
487 }
std::vector< int > iradMaxBH_
int layerIndex(int lay, bool reco) const
bool tileTrapezoid() const
std::vector< int > iradMinBH_
const HGCalParameters * hgpar_

◆ getTrForm()

HGCalParameters::hgtrform HGCalDDDConstants::getTrForm ( unsigned int  k) const
inline

Definition at line 79 of file HGCalDDDConstants.h.

References HGCalParameters::getTrForm(), hgpar_, and dqmdumpme::k.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getSummary(), and HGCalDDDConstants().

79 { return hgpar_->getTrForm(k); }
hgtrform getTrForm(unsigned int k) const
const HGCalParameters * hgpar_

◆ getTrFormN()

unsigned int HGCalDDDConstants::getTrFormN ( ) const
inline

Definition at line 80 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::trformIndex_.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getSummary(), and HGCalDDDConstants().

80 { return hgpar_->trformIndex_.size(); }
std::vector< uint32_t > trformIndex_
const HGCalParameters * hgpar_

◆ getTrForms()

std::vector< HGCalParameters::hgtrform > HGCalDDDConstants::getTrForms ( ) const

Definition at line 489 of file HGCalDDDConstants.cc.

References HGCalParameters::getTrForm(), hgpar_, dqmdumpme::k, and HGCalParameters::trformIndex_.

Referenced by HGCalTrackCollectionProducer::beginLuminosityBlock().

489  {
490  std::vector<HGCalParameters::hgtrform> mytrs;
491  for (unsigned int k = 0; k < hgpar_->trformIndex_.size(); ++k)
492  mytrs.emplace_back(hgpar_->getTrForm(k));
493  return mytrs;
494 }
std::vector< uint32_t > trformIndex_
hgtrform getTrForm(unsigned int k) const
const HGCalParameters * hgpar_

◆ getTypeHex()

int HGCalDDDConstants::getTypeHex ( int  layer,
int  waferU,
int  waferV 
) const

Definition at line 505 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, HGCalParameters::typesInLayers_, waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferTypeL_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

Referenced by HGCalGeometryLoader::build(), HGCalTriggerGeometryV9Imp2::detIdWaferType(), HGCalTriggerGeometryV9Imp3::detIdWaferType(), and HGCalTopology::neighbors().

505  {
506  // Get the module type for a silicon wafer
507  if (waferHexagon8()) {
509  return ((itr == hgpar_->typesInLayers_.end() ? 2 : hgpar_->waferTypeL_[itr->second]));
510  } else {
511  return -1;
512  }
513 }
int32_t waferU(const int32_t index)
bool waferHexagon8() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
wafer_map typesInLayers_
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_

◆ getTypeTrap()

int HGCalDDDConstants::getTypeTrap ( int  layer) const

Definition at line 496 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, HGCalParameters::scintType(), and tileTrapezoid().

Referenced by HGCalTriggerGeometryV9Imp3::detIdWaferType().

496  {
497  // Get the module type for scinitllator
498  if (tileTrapezoid()) {
499  return hgpar_->scintType(layer);
500  } else {
501  return -1;
502  }
503 }
int scintType(const int layer) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
bool tileTrapezoid() const
const HGCalParameters * hgpar_

◆ getUVMax()

int HGCalDDDConstants::getUVMax ( int  type) const
inline

◆ getXY()

std::pair< double, double > HGCalDDDConstants::getXY ( int  layer,
double  x,
double  y,
bool  forwd 
) const

Definition at line 515 of file HGCalDDDConstants.cc.

References HGCalParameters::firstLayer_, hgpar_, pixelTopology::layer, HGCalParameters::layerRotV_, HGCalParameters::layerType_, x, testProducerWithPsetDescEmpty_cfi::x1, y, and testProducerWithPsetDescEmpty_cfi::y1.

Referenced by localToGlobal8(), locateCell(), waferFromPosition(), waferLocal2Global(), and waferPosition().

515  {
516  int ll = layer - hgpar_->firstLayer_;
517  double x0(x), y0(y);
518  if ((!hgpar_->layerType_.empty()) && (ll < static_cast<int>(hgpar_->layerRotV_.size()))) {
519  if (forwd) {
520  x0 = x * hgpar_->layerRotV_[ll].first + y * hgpar_->layerRotV_[ll].second;
521  y0 = y * hgpar_->layerRotV_[ll].first - x * hgpar_->layerRotV_[ll].second;
522  } else {
523  x0 = x * hgpar_->layerRotV_[ll].first - y * hgpar_->layerRotV_[ll].second;
524  y0 = y * hgpar_->layerRotV_[ll].first + x * hgpar_->layerRotV_[ll].second;
525  }
526  }
527 #ifdef EDM_ML_DEBUG
528  double x1(x0), y1(y0);
529  if (ll < static_cast<int>(hgpar_->layerRotV_.size())) {
530  if (forwd) {
531  x1 = x0 * hgpar_->layerRotV_[ll].first - y0 * hgpar_->layerRotV_[ll].second;
532  y1 = y0 * hgpar_->layerRotV_[ll].first + x0 * hgpar_->layerRotV_[ll].second;
533  } else {
534  x1 = x0 * hgpar_->layerRotV_[ll].first + y0 * hgpar_->layerRotV_[ll].second;
535  y1 = y0 * hgpar_->layerRotV_[ll].first - x0 * hgpar_->layerRotV_[ll].second;
536  }
537  }
538  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants: Layer " << layer << ":" << ll << " mode " << forwd << " x " << x
539  << ":" << x0 << ":" << x1 << " y " << y << ":" << y0 << ":" << y1;
540 #endif
541  return std::make_pair(x0, y0);
542 }
Log< level::Info, true > LogVerbatim
std::vector< std::pair< double, double > > layerRotV_
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::vector< int > layerType_
const HGCalParameters * hgpar_

◆ isHalfCell()

bool HGCalDDDConstants::isHalfCell ( int  waferType,
int  cell 
) const

Definition at line 544 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseHalf_, HGCalParameters::cellFineHalf_, hgpar_, and waferType().

Referenced by hgcal::RecHitTools::isHalfCell().

544  {
545  if (waferType < 1 || cell < 0)
546  return false;
547  return waferType == 2 ? hgpar_->cellCoarseHalf_[cell] : hgpar_->cellFineHalf_[cell];
548 }
std::vector< bool > cellCoarseHalf_
std::vector< bool > cellFineHalf_
int waferType(DetId const &id, bool fromFile) const
const HGCalParameters * hgpar_

◆ isValidCell()

bool HGCalDDDConstants::isValidCell ( int  layindex,
int  wafer,
int  cell 
) const
private

Definition at line 1943 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, HGCSiliconDetId::HGCalFine, hgpar_, mps_fire::result, HGCalParameters::rMaxLayHex_, HGCalParameters::rMinLayHex_, findQualityFiles::rr, mathSSE::sqrt(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::waferTypeT_, x, and y.

Referenced by isValidHex().

1943  {
1944  // Calculate the position of the cell
1945  // Works for options HGCalHexagon/HGCalHexagonFull
1946  double x = hgpar_->waferPosX_[wafer];
1947  double y = hgpar_->waferPosY_[wafer];
1948  if (hgpar_->waferTypeT_[wafer] - 1 == HGCSiliconDetId::HGCalFine) {
1949  x += hgpar_->cellFineX_[cell];
1950  y += hgpar_->cellFineY_[cell];
1951  } else {
1952  x += hgpar_->cellCoarseX_[cell];
1953  y += hgpar_->cellCoarseY_[cell];
1954  }
1955  double rr = sqrt(x * x + y * y);
1956  bool result = ((rr >= hgpar_->rMinLayHex_[lay - 1]) && (rr <= hgpar_->rMaxLayHex_[lay - 1]) &&
1957  (wafer < static_cast<int>(hgpar_->waferPosX_.size())));
1958 #ifdef EDM_ML_DEBUG
1959  if (!result)
1960  edm::LogVerbatim("HGCalGeom") << "Input " << lay << ":" << wafer << ":" << cell << " Position " << x << ":" << y
1961  << ":" << rr << " Compare Limits " << hgpar_->rMinLayHex_[lay - 1] << ":"
1962  << hgpar_->rMaxLayHex_[lay - 1] << " Flag " << result;
1963 #endif
1964  return result;
1965 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
std::vector< double > cellFineY_
std::vector< double > cellCoarseX_
T sqrt(T t)
Definition: SSEVec.h:19
std::vector< double > rMinLayHex_
std::vector< double > rMaxLayHex_
std::vector< double > cellFineX_
std::vector< int > waferTypeT_
std::vector< double > cellCoarseY_
const HGCalParameters * hgpar_
std::vector< double > waferPosX_

◆ isValidCell8()

bool HGCalDDDConstants::isValidCell8 ( int  lay,
int  waferU,
int  waferV,
int  cellU,
int  cellV,
int  type 
) const
private

Definition at line 1967 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseIndex_, HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineIndex_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, HGCalParameters::firstLayer_, HGCalWaferMask::goodCell(), HGCalGeometryMode::Hexagon8Cassette, hgpar_, mode_, N, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, mps_fire::result, HGCalParameters::rMaxLayHex_, HGCalParameters::rMinLayHex_, findQualityFiles::rr, mathSSE::sqrt(), waferHexagon8File(), waferPositionNoRot(), waferTypeRotation(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, geometryCSVtoXML::xy, and y.

Referenced by isValidHex8().

1967  {
1968  bool result(false);
1969  auto partn = waferTypeRotation(lay, waferU, waferV, false, false);
1970 #ifdef EDM_ML_DEBUG
1971  edm::LogVerbatim("HGCalGeom") << "waferHexagon8 " << waferHexagon8File() << ":" << mode_ << ":"
1972  << HGCalGeometryMode::Hexagon8Cassette << " part " << partn.first << ":"
1973  << partn.second;
1974 #endif
1976  result = HGCalWaferMask::goodCell(cellU, cellV, partn.first);
1977 #ifdef EDM_ML_DEBUG
1978  edm::LogVerbatim("HGCalGeom") << "Input " << lay << ":" << waferU << ":" << waferV << ":" << cellU << ":" << cellV
1979  << " Result " << result << " from goodCell";
1980 #endif
1981  } else {
1982  float x(0), y(0);
1983  int kndx = cellV * 100 + cellU;
1984  if (type == 0) {
1985  auto ktr = hgpar_->cellFineIndex_.find(kndx);
1986  if (ktr != hgpar_->cellFineIndex_.end()) {
1987  x = hgpar_->cellFineX_[ktr->second];
1988  y = hgpar_->cellFineY_[ktr->second];
1989  }
1990 #ifdef EDM_ML_DEBUG
1991  edm::LogVerbatim("HGCalGeom") << "Fine " << cellU << ":" << cellV << ":" << kndx << ":" << x << ":" << y << ":"
1992  << (ktr != hgpar_->cellFineIndex_.end());
1993 #endif
1994  } else {
1995  auto ktr = hgpar_->cellCoarseIndex_.find(kndx);
1996  if (ktr != hgpar_->cellCoarseIndex_.end()) {
1997  x = hgpar_->cellCoarseX_[ktr->second];
1998  y = hgpar_->cellCoarseY_[ktr->second];
1999  }
2000 #ifdef EDM_ML_DEBUG
2001  edm::LogVerbatim("HGCalGeom") << "Coarse " << cellU << ":" << cellV << ":" << kndx << ":" << x << ":" << y << ":"
2002  << (ktr != hgpar_->cellCoarseIndex_.end());
2003 #endif
2004  }
2005  const auto& xy = waferPositionNoRot(lay, waferU, waferV, true, false);
2006  x += xy.first;
2007  y += xy.second;
2008 #ifdef EDM_ML_DEBUG
2009  edm::LogVerbatim("HGCalGeom") << "With wafer (" << waferU << "," << waferV << ") " << x << ":" << y;
2010 #endif
2011  double rr = sqrt(x * x + y * y);
2012  int ll = lay - hgpar_->firstLayer_;
2013  double tol = waferHexagon8File() ? 0.5 : 0.0;
2014  result = (((rr + tol) >= hgpar_->rMinLayHex_[ll]) && (rr <= hgpar_->rMaxLayHex_[ll]));
2015 #ifdef EDM_ML_DEBUG
2016  edm::LogVerbatim("HGCalGeom") << "Input " << lay << ":" << ll << ":" << waferU << ":" << waferV << ":" << cellU
2017  << ":" << cellV << " Position " << x << ":" << y << ":" << rr << " Compare Limits "
2018  << hgpar_->rMinLayHex_[ll] << ":" << hgpar_->rMaxLayHex_[ll] << " Flag " << result
2019  << " from Radius Limits";
2020 #endif
2021  if (result && waferHexagon8File()) {
2022  int N = (type == 0) ? hgpar_->nCellsFine_ : hgpar_->nCellsCoarse_;
2023  result = HGCalWaferMask::goodCell(cellU, cellV, N, partn.first, partn.second);
2024 #ifdef EDM_ML_DEBUG
2025  edm::LogVerbatim("HGCalGeom") << "Input " << lay << ":" << waferU << ":" << waferV << ":" << cellU << ":" << cellV
2026  << " N " << N << " part " << partn.first << ":" << partn.second << " Result "
2027  << result << " from goodCell";
2028 #endif
2029  }
2030  }
2031  return result;
2032 }
Log< level::Info, true > LogVerbatim
wafer_map cellFineIndex_
int32_t waferU(const int32_t index)
static bool goodCell(int u, int v, int N, int type, int rotn)
std::vector< double > cellFineY_
std::pair< int, int > waferTypeRotation(int layer, int waferU, int waferV, bool fromFile, bool debug) const
std::vector< double > cellCoarseX_
T sqrt(T t)
Definition: SSEVec.h:19
std::vector< double > rMinLayHex_
std::vector< double > rMaxLayHex_
#define N
Definition: blowfish.cc:9
wafer_map cellCoarseIndex_
bool waferHexagon8File() const
std::vector< double > cellFineX_
const HGCalGeometryMode::GeometryMode mode_
int32_t waferV(const int32_t index)
std::vector< double > cellCoarseY_
const HGCalParameters * hgpar_
std::pair< double, double > waferPositionNoRot(int lay, int waferU, int waferV, bool reco, bool debug) const

◆ isValidHex()

bool HGCalDDDConstants::isValidHex ( int  lay,
int  mod,
int  cell,
bool  reco 
) const

Definition at line 550 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellFineX_, HGCalParameters::copiesInLayers_, HGCSiliconDetId::HGCalFine, hgpar_, isValidCell(), layers(), maxCells(), mod(), mps_fire::result, findQualityFiles::size, HGCalParameters::waferCopy_, waferHexagon6(), and HGCalParameters::waferTypeT_.

Referenced by HGCNumberingScheme::getUnitID(), and HGCalTopology::valid().

550  {
551  // Check validity for a layer|wafer|cell of pre-TDR version
552  bool result(false), resultMod(false);
553  int cellmax(0);
554  if (waferHexagon6()) {
555  int32_t copyNumber = hgpar_->waferCopy_[mod];
556  result = ((lay > 0 && lay <= static_cast<int>(layers(reco))));
557  if (result) {
558  const int32_t lay_idx = reco ? (lay - 1) * 3 + 1 : lay;
559  const auto& the_modules = hgpar_->copiesInLayers_[lay_idx];
560  auto moditr = the_modules.find(copyNumber);
561  result = resultMod = (moditr != the_modules.end());
562 #ifdef EDM_ML_DEBUG
563  if (!result)
564  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants: Layer " << lay << ":" << lay_idx << " Copy " << copyNumber
565  << ":" << mod << " Flag " << result;
566 #endif
567  if (result) {
568  if (moditr->second >= 0) {
569  if (mod >= static_cast<int>(hgpar_->waferTypeT_.size()))
570  edm::LogWarning("HGCalGeom") << "Module no. out of bound for " << mod << " to be compared with "
571  << (hgpar_->waferTypeT_).size() << " ***** ERROR *****";
572  cellmax = ((hgpar_->waferTypeT_[mod] - 1 == HGCSiliconDetId::HGCalFine)
573  ? static_cast<int>(hgpar_->cellFineX_.size())
574  : static_cast<int>(hgpar_->cellCoarseX_.size()));
575  result = (cell >= 0 && cell <= cellmax);
576  } else {
577  result = isValidCell(lay_idx, mod, cell);
578  }
579  }
580  }
581  }
582 
583 #ifdef EDM_ML_DEBUG
584  if (!result)
585  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants: Layer " << lay << ":"
586  << (lay > 0 && (lay <= static_cast<int>(layers(reco)))) << " Module " << mod << ":"
587  << resultMod << " Cell " << cell << ":" << cellmax << ":"
588  << (cell >= 0 && cell <= cellmax) << ":" << maxCells(reco);
589 #endif
590  return result;
591 }
size
Write out results.
Log< level::Info, true > LogVerbatim
bool isValidCell(int layindex, int wafer, int cell) const
layer_map copiesInLayers_
bool waferHexagon6() const
std::vector< double > cellCoarseX_
int maxCells(bool reco) const
unsigned int layers(bool reco) const
std::vector< double > cellFineX_
fixed size matrix
std::vector< int > waferCopy_
std::vector< int > waferTypeT_
Log< level::Warning, false > LogWarning
const HGCalParameters * hgpar_
T mod(const T &a, const T &b)
Definition: ecalDccMap.h:4

◆ isValidHex8() [1/2]

bool HGCalDDDConstants::isValidHex8 ( int  lay,
int  waferU,
int  waferV,
bool  fullAndPart 
) const

Definition at line 593 of file HGCalDDDConstants.cc.

References fullAndPart_, hgpar_, pixelTopology::layer, HGCalParameters::typesInLayers_, HGCalTypes::WaferCornerMin, waferIn_, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalParameters::waferMaskMode_, HGCalTypes::WaferOut, and HGCalParameters::waferTypes_.

Referenced by HGCalTopology::addHGCSiliconId(), isValidHex8(), HGCalSD::setDetUnitId(), HGCalTopology::valid(), and HGCalTopology::validModule().

593  {
594  // Check validity for a layer|wafer|cell of post-TDR version
595  int indx = HGCalWaferIndex::waferIndex(layer, modU, modV);
596  auto itr = hgpar_->typesInLayers_.find(indx);
597 #ifdef EDM_ML_DEBUG
598  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants::isValidHex8:WaferType " << layer << ":" << modU << ":" << modV
599  << ":" << indx << " Test " << (itr != hgpar_->typesInLayers_.end());
600 #endif
601  if (itr == hgpar_->typesInLayers_.end()) {
602 #ifdef EDM_ML_DEBUG
603  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot find " << layer << ":" << modU << ":" << modV
604  << " in wadferIndex";
605 #endif
606  return false;
607  }
608 
609  if (fullAndPart_) {
610  auto ktr = hgpar_->waferInfoMap_.find(indx);
611 #ifdef EDM_ML_DEBUG
612  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants::isValidHex8:WaferInfoMap " << layer << ":" << modU << ":"
613  << modV << ":" << indx << " Test " << (ktr != hgpar_->waferInfoMap_.end());
614 #endif
615  if (ktr == hgpar_->waferInfoMap_.end()) {
616 #ifdef EDM_ML_DEBUG
617  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot find " << layer << ":" << modU << ":" << modV
618  << " in wadferInfoMap";
619 #endif
620  return false;
621  }
622  } else {
623  auto jtr = waferIn_.find(indx);
624 #ifdef EDM_ML_DEBUG
625  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants::isValidHex8:WaferIn " << jtr->first << ":" << jtr->second;
626 #endif
627  if (!(jtr->second)) {
628 #ifdef EDM_ML_DEBUG
629  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot find " << layer << ":" << modU << ":" << modV
630  << " in wadferIn";
631 #endif
632  return false;
633  }
634  }
635 
636  if (fullAndPart || fullAndPart_) {
637  auto ktr = hgpar_->waferTypes_.find(indx);
638  if (ktr != hgpar_->waferTypes_.end()) {
639  if (hgpar_->waferMaskMode_ > 0) {
640  if (ktr->second.first == HGCalTypes::WaferOut) {
641 #ifdef EDM_ML_DEBUG
642  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot find " << layer << ":" << modU << ":" << modV
643  << " due to WaferOut";
644 #endif
645  return false;
646  }
647  } else {
648  if (ktr->second.first < HGCalTypes::WaferCornerMin) {
649 #ifdef EDM_ML_DEBUG
650  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot find " << layer << ":" << modU << ":" << modV
651  << " due to WaferCornerMin";
652 #endif
653  return false;
654  }
655  }
656  }
657  }
658  return true;
659 }
Log< level::Info, true > LogVerbatim
static constexpr int32_t WaferOut
Definition: HGCalTypes.h:56
std::unordered_map< int32_t, bool > waferIn_
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
waferT_map waferTypes_
wafer_map typesInLayers_
static constexpr int32_t WaferCornerMin
Definition: HGCalTypes.h:74
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ isValidHex8() [2/2]

bool HGCalDDDConstants::isValidHex8 ( int  lay,
int  modU,
int  modV,
int  cellU,
int  cellV,
bool  fullAndPart 
) const

Definition at line 661 of file HGCalDDDConstants.cc.

References hgpar_, isValidCell8(), isValidHex8(), pixelTopology::layer, N, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, HGCalParameters::typesInLayers_, HGCalWaferIndex::waferIndex(), and HGCalParameters::waferTypeL_.

661  {
662  // First check validity for a layer|wafer| of post TDR version
663 #ifdef EDM_ML_DEBUG
664  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: layer|wafer " << layer << ":" << modU << ":" << modV << ":"
665  << fullAndPart << " Valid " << isValidHex8(layer, modU, modV, fullAndPart);
666 #endif
667  if (!isValidHex8(layer, modU, modV, fullAndPart))
668  return false;
669  int indx = HGCalWaferIndex::waferIndex(layer, modU, modV);
670  auto itr = hgpar_->typesInLayers_.find(indx);
671  int type = hgpar_->waferTypeL_[itr->second];
672  int N = ((hgpar_->waferTypeL_[itr->second] == 0) ? hgpar_->nCellsFine_ : hgpar_->nCellsCoarse_);
673 #ifdef EDM_ML_DEBUG
674  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants::isValidHex8:Cell " << cellU << ":" << cellV << ":" << N
675  << " Tests " << (cellU >= 0) << ":" << (cellU < 2 * N) << ":" << (cellV >= 0) << ":"
676  << (cellV < 2 * N) << ":" << ((cellV - cellU) < N) << ":" << ((cellU - cellV) <= N);
677 #endif
678  if ((cellU < 0) || (cellU >= 2 * N) || (cellV < 0) || (cellV >= 2 * N)) {
679 #ifdef EDM_ML_DEBUG
680  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot statisfy Cell 1 condition " << cellU << ":" << cellV
681  << ":" << N;
682 #endif
683  return false;
684  }
685  if (((cellV - cellU) >= N) || ((cellU - cellV) > N)) {
686 #ifdef EDM_ML_DEBUG
687  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:: Cannot statisfy Cell 2 condition " << cellU << ":" << cellV
688  << ":" << N;
689 #endif
690  return false;
691  }
692  return isValidCell8(layer, modU, modV, cellU, cellV, type);
693 }
Log< level::Info, true > LogVerbatim
bool isValidHex8(int lay, int waferU, int waferV, bool fullAndPart) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
bool isValidCell8(int lay, int waferU, int waferV, int cellU, int cellV, int type) const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
#define N
Definition: blowfish.cc:9
wafer_map typesInLayers_
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_

◆ isValidTrap()

bool HGCalDDDConstants::isValidTrap ( int  zside,
int  lay,
int  ieta,
int  iphi 
) const

Definition at line 695 of file HGCalDDDConstants.cc.

References getIndex(), hgpar_, LEDCalibrationChannels::iphi, HGCalParameters::iradMaxBH_, HGCalParameters::iradMinBH_, pixelTopology::layer, convertSQLiteXML::ok, HGCalParameters::scintCells(), tileExist(), trapezoidFile(), validateGeometry_cfg::valid, and ecaldqm::zside().

Referenced by HGCalTopology::addHGCSCintillatorId(), HGCScintSD::setDetUnitId(), and HGCalTopology::valid().

695  {
696  // Check validity for a layer|eta|phi of scintillator
697  const auto& indx = getIndex(layer, true);
698  if (indx.first < 0)
699  return false;
700  bool ok = ((irad >= hgpar_->iradMinBH_[indx.first]) && (irad <= (hgpar_->iradMaxBH_[indx.first] + 1)) && (iphi > 0) &&
701  (iphi <= hgpar_->scintCells(layer)));
702  bool valid = ((ok && trapezoidFile()) ? tileExist(zside, layer, irad, iphi) : ok);
703 #ifdef EDM_ML_DEBUG
704  bool tileEx = trapezoidFile() ? tileExist(zside, layer, irad, iphi) : true;
705  edm::LogVerbatim("HGCalGeomT") << "HGCalDDDConstants::isValidityTrap: Input " << zside << ":" << layer << ":" << irad
706  << ":" << iphi << " Range on Ring " << hgpar_->iradMinBH_[indx.first] << ":"
707  << (hgpar_->iradMaxBH_[indx.first] + 1)
708  << " Range on phi 0:" << hgpar_->scintCells(layer) << " tileExist " << tileEx
709  << " Valid " << ok << ":" << valid;
710 #endif
711  return valid;
712 }
std::vector< int > iradMaxBH_
Log< level::Info, true > LogVerbatim
bool trapezoidFile() const
int zside(DetId const &)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::pair< int, float > getIndex(int lay, bool reco) const
bool tileExist(int zside, int layer, int ring, int phi) const
std::vector< int > iradMinBH_
const HGCalParameters * hgpar_
int scintCells(const int layer) const

◆ lastLayer()

int HGCalDDDConstants::lastLayer ( bool  reco) const

Definition at line 714 of file HGCalDDDConstants.cc.

References HGCalParameters::firstLayer_, hgpar_, and tot_layers_.

Referenced by HGCalGeometry::neighborZ(), HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_(), and tileCount().

714  {
715  return (hgpar_->firstLayer_ + tot_layers_[static_cast<int>(reco)] - 1);
716 }
fixed size matrix
std::array< uint32_t, 2 > tot_layers_
const HGCalParameters * hgpar_

◆ layerFromIndex()

int HGCalDDDConstants::layerFromIndex ( int  index,
bool  reco 
) const
private

Definition at line 1923 of file HGCalDDDConstants.cc.

References HGCalParameters::depthLayerF_, spr::find(), HGCalParameters::firstLayer_, hgpar_, HGCalParameters::layerIndex_, and waferHexagon6().

Referenced by waferIndex().

1923  {
1924  int ll(-1);
1925  if (waferHexagon6() && reco) {
1926  ll = static_cast<int>(std::find(hgpar_->depthLayerF_.begin(), hgpar_->depthLayerF_.end(), index) -
1927  hgpar_->depthLayerF_.begin());
1928  if (ll == static_cast<int>(hgpar_->depthLayerF_.size()))
1929  ll = -1;
1930  } else {
1931  ll = static_cast<int>(std::find(hgpar_->layerIndex_.begin(), hgpar_->layerIndex_.end(), index) -
1932  hgpar_->layerIndex_.begin());
1933  if (ll == static_cast<int>(hgpar_->layerIndex_.size()))
1934  ll = -1;
1935  }
1936 #ifdef EDM_ML_DEBUG
1937  edm::LogVerbatim("HGCalGeom") << "LayerFromIndex for " << index << ":" << reco << ":" << waferHexagon6() << " is"
1938  << ll << ":" << (ll + hgpar_->firstLayer_);
1939 #endif
1940  return ((ll < 0) ? ll : (ll + hgpar_->firstLayer_));
1941 }
Log< level::Info, true > LogVerbatim
std::vector< int > depthLayerF_
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
bool waferHexagon6() const
std::vector< int > layerIndex_
fixed size matrix
const HGCalParameters * hgpar_

◆ layerIndex()

int HGCalDDDConstants::layerIndex ( int  lay,
bool  reco 
) const

Definition at line 720 of file HGCalDDDConstants.cc.

References HGCalParameters::depthIndex_, HGCalParameters::depthLayerF_, HGCalParameters::firstLayer_, hgpar_, HGCalParameters::layerIndex_, and waferHexagon6().

Referenced by HGCalGeometryLoader::build(), distFromEdgeHex(), distFromEdgeTrap(), getIndex(), getLayer(), getRangeR(), getREtaRange(), and HGCalGeometry::getSummary().

720  {
721  int ll = lay - hgpar_->firstLayer_;
722  if (ll < 0 || ll >= static_cast<int>(hgpar_->layerIndex_.size()))
723  return -1;
724  if (waferHexagon6()) {
725  if (reco && ll >= static_cast<int>(hgpar_->depthIndex_.size()))
726  return -1;
727  return (reco ? hgpar_->depthLayerF_[ll] : hgpar_->layerIndex_[ll]);
728  } else {
729  return (hgpar_->layerIndex_[ll]);
730  }
731 }
std::vector< int > depthLayerF_
bool waferHexagon6() const
std::vector< int > layerIndex_
fixed size matrix
std::vector< int > depthIndex_
const HGCalParameters * hgpar_

◆ layers()

unsigned int HGCalDDDConstants::layers ( bool  reco) const

◆ layersInit()

unsigned int HGCalDDDConstants::layersInit ( bool  reco) const

Definition at line 733 of file HGCalDDDConstants.cc.

References HGCalParameters::depthIndex_, hgpar_, and HGCalParameters::layerIndex_.

Referenced by HGCalDDDConstants().

733  {
734  return (reco ? hgpar_->depthIndex_.size() : hgpar_->layerIndex_.size());
735 }
std::vector< int > layerIndex_
fixed size matrix
std::vector< int > depthIndex_
const HGCalParameters * hgpar_

◆ layerType()

int HGCalDDDConstants::layerType ( int  lay) const
inline

Definition at line 95 of file HGCalDDDConstants.h.

References HGCalParameters::firstLayer_, hgpar_, HGCalParameters::layerType_, and HGCalTypes::WaferCenter.

Referenced by locateCell(), and waferFromPosition().

95  {
97  }
static constexpr int32_t WaferCenter
Definition: HGCalTypes.h:21
std::vector< int > layerType_
const HGCalParameters * hgpar_

◆ levelTop()

int HGCalDDDConstants::levelTop ( int  ind = 0) const
inline

Definition at line 108 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::levelT_.

Referenced by SimG4HGCalValidation::beginRun(), HGCScintSD::update(), HGCalSD::update(), HFNoseSD::update(), and HGCSD::update().

108 { return hgpar_->levelT_[ind]; }
std::vector< int > levelT_
const HGCalParameters * hgpar_

◆ localToGlobal8()

std::pair< float, float > HGCalDDDConstants::localToGlobal8 ( int  zside,
int  lay,
int  waferU,
int  waferV,
double  localX,
double  localY,
bool  reco,
bool  debug 
) const

Definition at line 737 of file HGCalDDDConstants.cc.

References debug, HGCalParameters::firstLayer_, HGCalCassette::getShift(), getXY(), HGCalGeometryMode::Hexagon8Cassette, hgcassette_, hgpar_, HGCalParameters::k_ScaleToDDD, HGCalParameters::layerType_, mode_, HGCalTypes::WaferCenterR, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, waferPositionNoRot(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, geometryCSVtoXML::xy, y, and ecaldqm::zside().

Referenced by HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), and HGCalGeometry::getNewCorners().

738  {
739  double x(localX), y(localY);
740  bool rotx =
742  if (debug)
743  edm::LogVerbatim("HGCalGeom") << "LocalToGlobal8 " << lay << ":" << (lay - hgpar_->firstLayer_) << ":" << rotx
744  << " Local (" << x << ":" << y << ") Reco " << reco;
745  if (!reco) {
748  }
749  const auto& xy = waferPositionNoRot(lay, waferU, waferV, reco, debug);
750  x += xy.first;
751  y += xy.second;
752  int indx = HGCalWaferIndex::waferIndex(lay, waferU, waferV);
753  auto ktr = hgpar_->waferInfoMap_.find(indx);
754  if ((mode_ == HGCalGeometryMode::Hexagon8Cassette) && (ktr != hgpar_->waferInfoMap_.end())) {
755  auto cshift = hgcassette_.getShift(lay, -zside, (ktr->second).cassette);
756  if (debug)
757  edm::LogVerbatim("HGCalGeom") << "Cassette " << (ktr->second).cassette << " Shift " << -(zside * cshift.first)
758  << ":" << cshift.second;
759  if (!reco) {
760  x -= ((HGCalParameters::k_ScaleToDDD)*zside * cshift.first);
761  y += ((HGCalParameters::k_ScaleToDDD)*cshift.second);
762  } else {
763  x -= (zside * cshift.first);
764  y += cshift.second;
765  }
766  }
767  if (debug)
768  edm::LogVerbatim("HGCalGeom") << "With wafer " << x << ":" << y << " by adding " << xy.first << ":" << xy.second;
769  return (rotx ? getXY(lay, x, y, false) : std::make_pair(x, y));
770 }
Log< level::Info, true > LogVerbatim
HGCalCassette hgcassette_
int32_t waferU(const int32_t index)
std::pair< double, double > getShift(int layer, int zside, int cassette) const
std::pair< double, double > getXY(int layer, double x, double y, bool forwd) const
int zside(DetId const &)
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
std::vector< int > layerType_
static constexpr int32_t WaferCenterR
Definition: HGCalTypes.h:27
static constexpr double k_ScaleToDDD
#define debug
Definition: HDRShower.cc:19
const HGCalGeometryMode::GeometryMode mode_
fixed size matrix
int32_t waferV(const int32_t index)
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
std::pair< double, double > waferPositionNoRot(int lay, int waferU, int waferV, bool reco, bool debug) const

◆ locateCell() [1/4]

std::pair< float, float > HGCalDDDConstants::locateCell ( int  cell,
int  lay,
int  type,
bool  reco 
) const

Definition at line 772 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, getIndex(), HGCSiliconDetId::HGCalFine, hgpar_, mps_fire::i, HGCalParameters::k_ScaleToDDD, waferHexagon6(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::waferTypeT_, x, and y.

Referenced by HGCalSimHitValidation::analyzeHits(), HGCalTBAnalyzer::analyzeSimHits(), cellInLayer(), HFNoseNumberingScheme::checkPosition(), HGCalNumberingScheme::checkPosition(), HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), HGCNumberingScheme::getLocalCoords(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), locateCell(), HGCScintSD::setDetUnitId(), and HGCalSD::setDetUnitId().

772  {
773  // type refers to wafer # for hexagon cell
774  float x(999999.), y(999999.);
775  const auto& index = getIndex(lay, reco);
776  int i = index.first;
777  if (i < 0)
778  return std::make_pair(x, y);
779  if (waferHexagon6()) {
780  x = hgpar_->waferPosX_[type];
781  y = hgpar_->waferPosY_[type];
782 #ifdef EDM_ML_DEBUG
783  float x0(x), y0(y);
784 #endif
786  x += hgpar_->cellFineX_[cell];
787  y += hgpar_->cellFineY_[cell];
788  } else {
789  x += hgpar_->cellCoarseX_[cell];
790  y += hgpar_->cellCoarseY_[cell];
791  }
792 #ifdef EDM_ML_DEBUG
793  edm::LogVerbatim("HGCalGeom") << "LocateCell (Wafer) " << x0 << ":" << y0 << " Final " << x << ":" << y;
794 #endif
795  if (!reco) {
798  }
799  }
800  return std::make_pair(x, y);
801 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
std::vector< double > cellFineY_
bool waferHexagon6() const
std::vector< double > cellCoarseX_
std::pair< int, float > getIndex(int lay, bool reco) const
static constexpr double k_ScaleToDDD
std::vector< double > cellFineX_
fixed size matrix
std::vector< int > waferTypeT_
std::vector< double > cellCoarseY_
const HGCalParameters * hgpar_
std::vector< double > waferPosX_

◆ locateCell() [2/4]

std::pair< float, float > HGCalDDDConstants::locateCell ( int  zside,
int  lay,
int  waferU,
int  waferV,
int  cellU,
int  cellV,
bool  reco,
bool  all,
bool  norot,
bool  debug 
) const

Definition at line 803 of file HGCalDDDConstants.cc.

References python.cmstools::all(), HGCalParameters::cellCoarseIndex_, HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineIndex_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, HGCalCell::cellPlacementIndex(), HGCalCell::cellPlacementOld, debug, HGCalParameters::firstLayer_, HGCalCassette::getShift(), getXY(), HGCalGeometryMode::Hexagon8Cassette, hgcassette_, hgcell_, hgpar_, HGCalParameters::k_ScaleToDDD, HGCalTypes::layerFrontBack(), layerType(), mode_, hgcalPerformanceValidation::orient, HGCalParameters::typesInLayers_, HGCalTypes::WaferCenterR, waferHexagon8File(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, waferPositionNoRot(), HGCalParameters::waferTypeL_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, geometryCSVtoXML::xy, y, and ecaldqm::zside().

805  {
806  double x(0), y(0);
807  int indx = HGCalWaferIndex::waferIndex(lay, waferU, waferV);
808  auto itr = hgpar_->typesInLayers_.find(indx);
809  int type = ((itr == hgpar_->typesInLayers_.end()) ? 2 : hgpar_->waferTypeL_[itr->second]);
810  int layertype = layerType(lay);
811  bool rotx = (norot) ? false : (layertype == HGCalTypes::WaferCenterR);
812  if (debug) {
813  edm::LogVerbatim("HGCalGeom") << "LocateCell " << lay << ":" << (lay - hgpar_->firstLayer_) << ":" << layertype
814  << ":" << rotx << ":" << waferU << ":" << waferV << ":" << indx << ":"
815  << (itr == hgpar_->typesInLayers_.end()) << ":" << type << " Flags " << reco << ":"
816  << all;
817  }
818  auto ktr = hgpar_->waferInfoMap_.end();
819  int place(HGCalCell::cellPlacementOld);
820  if (waferHexagon8File()) {
822  ktr = hgpar_->waferInfoMap_.find(indx);
823  if (ktr != hgpar_->waferInfoMap_.end())
824  place = HGCalCell::cellPlacementIndex(1, HGCalTypes::layerFrontBack(layertype), (ktr->second).orient);
825  }
826  auto xy = hgcell_->cellUV2XY2(cellU, cellV, place, type);
827  x = xy.first;
828  y = xy.second;
829  if (debug)
830  edm::LogVerbatim("HGCalGeom") << "Type " << type << " Place " << place << " Cell " << cellU << ":" << cellV
831  << " Position " << x << ":" << y;
832  } else {
833  int kndx = cellV * 100 + cellU;
834  if (type == 0) {
835  auto jtr = hgpar_->cellFineIndex_.find(kndx);
836  if (jtr != hgpar_->cellFineIndex_.end()) {
837  x = hgpar_->cellFineX_[jtr->second];
838  y = hgpar_->cellFineY_[jtr->second];
839  }
840  if (debug)
841  edm::LogVerbatim("HGCalGeom") << "Fine " << cellU << ":" << cellV << ":" << kndx << ":" << x << ":" << y << ":"
842  << (jtr != hgpar_->cellFineIndex_.end());
843  } else {
844  auto jtr = hgpar_->cellCoarseIndex_.find(kndx);
845  if (jtr != hgpar_->cellCoarseIndex_.end()) {
846  x = hgpar_->cellCoarseX_[jtr->second];
847  y = hgpar_->cellCoarseY_[jtr->second];
848  }
849  if (debug)
850  edm::LogVerbatim("HGCalGeom") << "Coarse " << cellU << ":" << cellV << ":" << kndx << ":" << x << ":" << y
851  << ":" << (jtr != hgpar_->cellCoarseIndex_.end());
852  }
853  }
854  if (!reco) {
857  }
858  if (all) {
859  const auto& xy = waferPositionNoRot(lay, waferU, waferV, reco, debug);
860  x += xy.first;
861  y += xy.second;
862  if ((mode_ == HGCalGeometryMode::Hexagon8Cassette) && (ktr != hgpar_->waferInfoMap_.end())) {
863  auto cshift = hgcassette_.getShift(lay, -zside, (ktr->second).cassette);
864  if (debug)
865  edm::LogVerbatim("HGCalGeom") << "Cassette " << (ktr->second).cassette << " Shift " << -(zside * cshift.first)
866  << ":" << cshift.second;
867  if (!reco) {
868  x -= ((HGCalParameters::k_ScaleToDDD)*zside * cshift.first);
869  y += ((HGCalParameters::k_ScaleToDDD)*cshift.second);
870  } else {
871  x -= (zside * cshift.first);
872  y += cshift.second;
873  }
874  }
875  if (debug)
876  edm::LogVerbatim("HGCalGeom") << "With wafer " << x << ":" << y << " by adding " << xy.first << ":" << xy.second;
877  }
878  return (rotx ? getXY(lay, x, y, false) : std::make_pair(x, y));
879 }
Log< level::Info, true > LogVerbatim
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:233
static constexpr int32_t cellPlacementOld
Definition: HGCalCell.h:25
wafer_map cellFineIndex_
HGCalCassette hgcassette_
def all(container)
workaround iterator generators for ROOT classes
Definition: cmstools.py:25
int32_t waferU(const int32_t index)
std::pair< double, double > getShift(int layer, int zside, int cassette) const
std::vector< double > cellFineY_
std::pair< double, double > getXY(int layer, double x, double y, bool forwd) const
int zside(DetId const &)
std::vector< double > cellCoarseX_
std::unique_ptr< HGCalCell > hgcell_
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
static constexpr int32_t WaferCenterR
Definition: HGCalTypes.h:27
static constexpr double k_ScaleToDDD
#define debug
Definition: HDRShower.cc:19
wafer_map cellCoarseIndex_
bool waferHexagon8File() const
std::vector< double > cellFineX_
wafer_map typesInLayers_
const HGCalGeometryMode::GeometryMode mode_
fixed size matrix
int32_t waferV(const int32_t index)
std::vector< double > cellCoarseY_
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_
std::pair< double, double > waferPositionNoRot(int lay, int waferU, int waferV, bool reco, bool debug) const
static int32_t layerFrontBack(int32_t layerOrient)
Definition: HGCalTypes.h:125
int layerType(int lay) const

◆ locateCell() [3/4]

std::pair< float, float > HGCalDDDConstants::locateCell ( const HGCSiliconDetId id,
bool  debug 
) const

Definition at line 881 of file HGCalDDDConstants.cc.

References debug, pixelTopology::layer, locateCell(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and ecaldqm::zside().

881  {
882  return locateCell(id.zside(), id.layer(), id.waferU(), id.waferV(), id.cellU(), id.cellV(), true, true, false, debug);
883 }
int32_t waferU(const int32_t index)
int zside(DetId const &)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
#define debug
Definition: HDRShower.cc:19
std::pair< float, float > locateCell(int cell, int lay, int type, bool reco) const
int32_t waferV(const int32_t index)

◆ locateCell() [4/4]

std::pair< float, float > HGCalDDDConstants::locateCell ( const HGCScintillatorDetId id,
bool  debug 
) const

Definition at line 885 of file HGCalDDDConstants.cc.

References debug, LEDCalibrationChannels::iphi, pixelTopology::layer, locateCellTrap(), and ecaldqm::zside().

885  {
886  return locateCellTrap(id.zside(), id.layer(), id.iradius(), id.iphi(), true, debug);
887 }
int zside(DetId const &)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::pair< float, float > locateCellTrap(int zside, int lay, int ieta, int iphi, bool reco, bool debug) const
#define debug
Definition: HDRShower.cc:19

◆ locateCellHex()

std::pair< float, float > HGCalDDDConstants::locateCellHex ( int  cell,
int  wafer,
bool  reco 
) const

Definition at line 889 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, HGCSiliconDetId::HGCalFine, hgpar_, HGCalParameters::k_ScaleToDDD, HGCalParameters::waferTypeT_, x, and y.

Referenced by HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), HGCalGeometry::getNewCorners(), and HGCalGeometry::getPosition().

889  {
890  float x(0), y(0);
891  if (hgpar_->waferTypeT_[wafer] - 1 == HGCSiliconDetId::HGCalFine) {
892  x = hgpar_->cellFineX_[cell];
893  y = hgpar_->cellFineY_[cell];
894  } else {
895  x = hgpar_->cellCoarseX_[cell];
896  y = hgpar_->cellCoarseY_[cell];
897  }
898  if (!reco) {
901  }
902  return std::make_pair(x, y);
903 }
std::vector< double > cellFineY_
std::vector< double > cellCoarseX_
static constexpr double k_ScaleToDDD
std::vector< double > cellFineX_
fixed size matrix
std::vector< int > waferTypeT_
std::vector< double > cellCoarseY_
const HGCalParameters * hgpar_

◆ locateCellTrap()

std::pair< float, float > HGCalDDDConstants::locateCellTrap ( int  zside,
int  lay,
int  ieta,
int  iphi,
bool  reco,
bool  debug 
) const

Definition at line 905 of file HGCalDDDConstants.cc.

References funct::abs(), HGCalParameters::cassettes_, angle_units::operators::convertRadToDeg(), funct::cos(), debug, HGCalParameters::firstLayer_, getIndex(), HGCalCassette::getShift(), hgcassette_, hgpar_, LEDCalibrationChannels::iphi, HGCalParameters::iradMaxBH_, HGCalParameters::iradMinBH_, HGCalParameters::k_ScaleToDDD, SiStripPI::max, SiStripPI::min, mode_, HGCalParameters::nphiCassette_, phi, HGCalParameters::phiOffset_, alignCSCRings::r, HGCalParameters::radiusLayer_, FastTimerService_cff::range, rangeR(), HGCalParameters::scintType(), funct::sin(), HGCalTileIndex::tileCassette(), HGCalGeometryMode::TrapezoidCassette, HGCalGeometryMode::TrapezoidFile, HGCalGeometryMode::TrapezoidModule, x, HGCalParameters::xLayerHex_, y, HGCalParameters::yLayerHex_, z, HGCalParameters::zLayerHex_, and ecaldqm::zside().

Referenced by HGCalSimHitValidation::analyzeHits(), HGCalGeometryLoader::build(), HGCalNumberingScheme::checkPosition(), and locateCell().

906  {
907  float x(0), y(0);
908  const auto& indx = getIndex(lay, reco);
909  if (indx.first >= 0) {
910  int ir = std::abs(irad);
911  int type = hgpar_->scintType(lay);
912  double phi = (iphi - 0.5) * indx.second;
913  double z = hgpar_->zLayerHex_[indx.first];
914  double r = 0.5 * (hgpar_->radiusLayer_[type][ir - 1] + hgpar_->radiusLayer_[type][ir]);
915  std::pair<double, double> range = rangeR(z, true);
916  if (debug)
917  edm::LogVerbatim("HGCalGeom") << "locateCellTrap:: Input " << lay << ":" << irad << ":" << iphi << ":" << reco
918  << " IR " << ir << ":" << hgpar_->iradMinBH_[indx.first] << ":"
919  << hgpar_->iradMaxBH_[indx.first] << " Type " << type << " Z " << indx.first << ":"
920  << z << " phi " << phi << ":" << convertRadToDeg(phi) << " R " << r << ":"
921  << range.first << ":" << range.second;
924  r = std::max(range.first, std::min(r, range.second));
925  x = r * std::cos(phi);
926  y = r * std::sin(phi);
927  int ll = lay - hgpar_->firstLayer_;
928  x += hgpar_->xLayerHex_[ll];
929  y += hgpar_->yLayerHex_[ll];
932  auto cshift = hgcassette_.getShift(lay, zside, cassette);
933  if (debug)
934  edm::LogVerbatim("HGCalGeom") << "Cassette " << cassette << " Shift " << -(zside * cshift.first) << ":"
935  << cshift.second;
936  x -= (zside * cshift.first);
937  y += cshift.second;
938  }
939  if (irad < 0)
940  x = -x;
941  }
942  if (!reco) {
945  }
946  return std::make_pair(x, y);
947 }
std::vector< int > iradMaxBH_
Log< level::Info, true > LogVerbatim
HGCalCassette hgcassette_
int scintType(const int layer) const
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
std::pair< double, double > getShift(int layer, int zside, int cassette) const
int zside(DetId const &)
std::vector< double > yLayerHex_
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::pair< int, float > getIndex(int lay, bool reco) const
int32_t tileCassette(int32_t, int32_t, int32_t, int32_t)
std::vector< double > zLayerHex_
std::pair< double, double > rangeR(double z, bool reco) const
static constexpr double k_ScaleToDDD
#define debug
Definition: HDRShower.cc:19
std::vector< int > iradMinBH_
const HGCalGeometryMode::GeometryMode mode_
fixed size matrix
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_
std::vector< double > xLayerHex_

◆ maskCell()

bool HGCalDDDConstants::maskCell ( const DetId id,
int  corners 
) const

Definition at line 949 of file HGCalDDDConstants.cc.

References DetId::det(), DetId::Forward, getUVMax(), HGCalWaferMask::goodCell(), HGCalGeometryMode::Hexagon8Cassette, hgpar_, l1ctLayer2EG_cff::id, HGCalParameters::k_CornerSize, HGCalTypes::k_OffsetRotation, pixelTopology::layer, gpuClustering::pixelStatus::mask, HGCalWaferMask::maskCell(), mode_, N, edm::second(), findQualityFiles::v, HGCalTypes::WaferFull, waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalParameters::waferTypes_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and LaserClient_cfi::wl.

Referenced by hgcal::RecHitTools::maskCell(), HGCalTopology::maskCell(), HFNoseSD::setDetUnitId(), HGCalSD::setDetUnitId(), and HGCalTopology::valid().

949  {
950  bool mask(false);
951  if (corners > 2 && corners <= static_cast<int>(HGCalParameters::k_CornerSize)) {
952  if (waferHexagon8()) {
953  int N(0), layer(0), waferU(0), waferV(0), u(0), v(0);
954  if (detId.det() == DetId::Forward) {
955  HFNoseDetId id(detId);
956  N = getUVMax(id.type());
957  layer = id.layer();
958  waferU = id.waferU();
959  waferV = id.waferV();
960  u = id.cellU();
961  v = id.cellV();
962  } else {
963  HGCSiliconDetId id(detId);
964  N = getUVMax(id.type());
965  layer = id.layer();
966  waferU = id.waferU();
967  waferV = id.waferV();
968  u = id.cellU();
969  v = id.cellV();
970  }
972  auto itr = hgpar_->waferTypes_.find(wl);
973  auto ktr = hgpar_->waferInfoMap_.find(wl);
974 #ifdef EDM_ML_DEBUG
975  edm::LogVerbatim("HGCalGeom") << "MaskCell: Layer " << layer << " Wafer " << waferU << ":" << waferV << " Index "
976  << wl << ":" << (itr != hgpar_->waferTypes_.end()) << ":"
977  << (ktr != hgpar_->waferInfoMap_.end());
978 #endif
980  int part = (ktr != hgpar_->waferInfoMap_.end()) ? (ktr->second).part : HGCalTypes::WaferFull;
982  } else if (itr != hgpar_->waferTypes_.end()) {
983  if ((itr->second).second <= HGCalTypes::k_OffsetRotation)
984  mask = HGCalWaferMask::maskCell(u, v, N, (itr->second).first, (itr->second).second, corners);
985  else
987  u, v, N, (itr->second).first, ((itr->second).second - HGCalTypes::k_OffsetRotation)));
988  }
989  }
990  }
991  return mask;
992 }
Log< level::Info, true > LogVerbatim
static constexpr int k_OffsetRotation
Definition: HGCalTypes.h:91
int32_t waferU(const int32_t index)
static bool goodCell(int u, int v, int N, int type, int rotn)
bool waferHexagon8() const
constexpr uint32_t mask
Definition: gpuClustering.h:26
static constexpr uint32_t k_CornerSize
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
U second(std::pair< T, U > const &p)
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
int getUVMax(int type) const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
waferT_map waferTypes_
#define N
Definition: blowfish.cc:9
part
Definition: HCALResponse.h:20
const HGCalGeometryMode::GeometryMode mode_
int32_t waferV(const int32_t index)
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
static bool maskCell(int u, int v, int N, int ncor, int fcor, int corners)

◆ maxCells() [1/2]

int HGCalDDDConstants::maxCells ( bool  reco) const

Definition at line 994 of file HGCalDDDConstants.cc.

References hgcalTBTopologyTester_cfi::cells, HGCalParameters::depth_, hgpar_, mps_fire::i, HGCalParameters::layer_, and layers().

Referenced by HGCalDDDConstants(), HGCalTopology::HGCalTopology(), and isValidHex().

994  {
995  int cells(0);
996  for (unsigned int i = 0; i < layers(reco); ++i) {
997  int lay = reco ? hgpar_->depth_[i] : hgpar_->layer_[i];
998  if (cells < maxCells(lay, reco))
999  cells = maxCells(lay, reco);
1000  }
1001  return cells;
1002 }
std::vector< int > layer_
std::vector< int > depth_
int maxCells(bool reco) const
unsigned int layers(bool reco) const
fixed size matrix
const HGCalParameters * hgpar_

◆ maxCells() [2/2]

int HGCalDDDConstants::maxCells ( int  lay,
bool  reco 
) const

Definition at line 1004 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellFineX_, hgcalTBTopologyTester_cfi::cells, HGCalParameters::defineFull_, HGCalParameters::firstLayer_, getIndex(), HGCSiliconDetId::HGCalCoarseThick, HGCSiliconDetId::HGCalFine, hgpar_, dqmdumpme::k, SiStripPI::max, N, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, HGCalParameters::scintCells(), tileTrapezoid(), HGCalParameters::typesInLayers_, HGCalParameters::waferCopy_, waferHexagon6(), waferHexagon8(), HGCalWaferIndex::waferIndex(), waferInLayerTest(), HGCalParameters::waferTypeL_, HGCalParameters::waferTypeT_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

1004  {
1005  const auto& index = getIndex(lay, reco);
1006  if (index.first < 0)
1007  return 0;
1008  if (waferHexagon6()) {
1009  unsigned int cells(0);
1010  for (unsigned int k = 0; k < hgpar_->waferTypeT_.size(); ++k) {
1011  if (waferInLayerTest(k, index.first, hgpar_->defineFull_)) {
1012  unsigned int cell = (hgpar_->waferTypeT_[k] - 1 == HGCSiliconDetId::HGCalFine) ? (hgpar_->cellFineX_.size())
1013  : (hgpar_->cellCoarseX_.size());
1014  if (cell > cells)
1015  cells = cell;
1016  }
1017  }
1018  return static_cast<int>(cells);
1019  } else if (waferHexagon8()) {
1020  int cells(0);
1021  for (unsigned int k = 0; k < hgpar_->waferCopy_.size(); ++k) {
1022  if (waferInLayerTest(k, index.first, hgpar_->defineFull_)) {
1026  : hgpar_->waferTypeL_[itr->second]);
1028  cells = std::max(cells, 3 * N * N);
1029  }
1030  }
1031  return cells;
1032  } else if (tileTrapezoid()) {
1033  return hgpar_->scintCells(index.first + hgpar_->firstLayer_);
1034  } else {
1035  return 0;
1036  }
1037 }
bool waferInLayerTest(int wafer, int lay, bool full) const
int32_t waferU(const int32_t index)
bool waferHexagon8() const
bool waferHexagon6() const
std::vector< double > cellCoarseX_
std::pair< int, float > getIndex(int lay, bool reco) const
bool tileTrapezoid() const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
#define N
Definition: blowfish.cc:9
std::vector< double > cellFineX_
wafer_map typesInLayers_
fixed size matrix
int32_t waferV(const int32_t index)
std::vector< int > waferCopy_
std::vector< int > waferTypeT_
const HGCalParameters * hgpar_
int scintCells(const int layer) const
std::vector< int > waferTypeL_

◆ maxCellUV()

int HGCalDDDConstants::maxCellUV ( ) const
inline

Definition at line 110 of file HGCalDDDConstants.h.

References hgpar_, HGCalParameters::nCellsFine_, and tileTrapezoid().

Referenced by HGCalTopology::HGCalTopology().

110 { return (tileTrapezoid() ? hgpar_->nCellsFine_ : 2 * hgpar_->nCellsFine_); }
bool tileTrapezoid() const
const HGCalParameters * hgpar_

◆ maxModules()

int HGCalDDDConstants::maxModules ( ) const
inline

Definition at line 113 of file HGCalDDDConstants.h.

References modHalf_.

113 { return modHalf_; }

◆ maxModulesPerLayer()

int HGCalDDDConstants::maxModulesPerLayer ( ) const
inline

Definition at line 114 of file HGCalDDDConstants.h.

References maxWafersPerLayer_.

114 { return maxWafersPerLayer_; }

◆ maxRows()

int HGCalDDDConstants::maxRows ( int  lay,
bool  reco 
) const

Definition at line 1039 of file HGCalDDDConstants.cc.

References HGCalParameters::defineFull_, getIndex(), hgpar_, mps_fire::i, dqmdumpme::k, HGCalParameters::waferCopy_, waferHexagon6(), waferHexagon8(), waferInLayerTest(), and HGCalParameters::waferUVMaxLayer_.

1039  {
1040  int kymax(0);
1041  const auto& index = getIndex(lay, reco);
1042  int i = index.first;
1043  if (i < 0)
1044  return kymax;
1045  if (waferHexagon6()) {
1046  for (unsigned int k = 0; k < hgpar_->waferCopy_.size(); ++k) {
1048  int ky = ((hgpar_->waferCopy_[k]) / 100) % 100;
1049  if (ky > kymax)
1050  kymax = ky;
1051  }
1052  }
1053  } else if (waferHexagon8()) {
1054  kymax = 1 + 2 * hgpar_->waferUVMaxLayer_[index.first];
1055  }
1056  return kymax;
1057 }
bool waferInLayerTest(int wafer, int lay, bool full) const
bool waferHexagon8() const
bool waferHexagon6() const
std::vector< int > waferUVMaxLayer_
std::pair< int, float > getIndex(int lay, bool reco) const
fixed size matrix
std::vector< int > waferCopy_
const HGCalParameters * hgpar_

◆ minSlope()

double HGCalDDDConstants::minSlope ( ) const
inline

Definition at line 116 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::slopeMin_.

Referenced by HGCScintSD::update(), HGCalSD::update(), HFNoseSD::update(), and HGCSD::update().

116 { return hgpar_->slopeMin_[0]; }
std::vector< double > slopeMin_
const HGCalParameters * hgpar_

◆ modifyUV()

int HGCalDDDConstants::modifyUV ( int  uv,
int  type1,
int  type2 
) const

Definition at line 1059 of file HGCalDDDConstants.cc.

Referenced by HGCalTopology::neighbors().

1059  {
1060  // Modify u/v for transition of type1 to type2
1061  return (((type1 == type2) || (type1 * type2 != 0)) ? uv : ((type1 == 0) ? (2 * uv + 1) / 3 : (3 * uv) / 2));
1062 }

◆ modules()

int HGCalDDDConstants::modules ( int  lay,
bool  reco 
) const

Definition at line 1064 of file HGCalDDDConstants.cc.

References dqmdumpme::first, getIndex(), max_modules_layer_, and dt_dqm_sourceclient_common_cff::reco.

Referenced by wafers().

1064  {
1065  if (getIndex(lay, reco).first < 0)
1066  return 0;
1067  else
1068  return max_modules_layer_[static_cast<int>(reco)][lay];
1069 }
Simrecovecs max_modules_layer_
std::pair< int, float > getIndex(int lay, bool reco) const
fixed size matrix

◆ modulesInit()

int HGCalDDDConstants::modulesInit ( int  lay,
bool  reco 
) const

Definition at line 1071 of file HGCalDDDConstants.cc.

References HGCalParameters::defineFull_, HGCalParameters::firstModule_, getIndex(), hgpar_, dqmdumpme::k, HGCalParameters::lastModule_, tileTrapezoid(), waferInLayerTest(), and HGCalParameters::waferPosX_.

Referenced by HGCalDDDConstants().

1071  {
1072  int nmod(0);
1073  const auto& index = getIndex(lay, reco);
1074  if (index.first < 0)
1075  return nmod;
1076  if (!tileTrapezoid()) {
1077  for (unsigned int k = 0; k < hgpar_->waferPosX_.size(); ++k) {
1078  if (waferInLayerTest(k, index.first, hgpar_->defineFull_))
1079  ++nmod;
1080  }
1081  } else {
1082  nmod = 1 + hgpar_->lastModule_[index.first] - hgpar_->firstModule_[index.first];
1083  }
1084  return nmod;
1085 }
bool waferInLayerTest(int wafer, int lay, bool full) const
std::vector< int > firstModule_
std::pair< int, float > getIndex(int lay, bool reco) const
bool tileTrapezoid() const
std::vector< int > lastModule_
fixed size matrix
const HGCalParameters * hgpar_
std::vector< double > waferPosX_

◆ mouseBite()

double HGCalDDDConstants::mouseBite ( bool  reco) const

Definition at line 1087 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::k_ScaleToDDD, and HGCalParameters::mouseBite_.

Referenced by HFNoseSD::update(), and HGCalSD::update().

1087  {
1089 }
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

◆ numberCells() [1/2]

int HGCalDDDConstants::numberCells ( bool  reco) const

Definition at line 1091 of file HGCalDDDConstants.cc.

References hgcalTBTopologyTester_cfi::cells, HGCalParameters::depth_, hgpar_, dqmdumpme::k, HGCalParameters::layer_, dt_dqm_sourceclient_common_cff::reco, HGCalGeomParameters::scintillatorCassette, HGCalGeomParameters::scintillatorFile, tileCount(), tileTrapezoid(), and HGCalParameters::waferMaskMode_.

Referenced by HGCalTopology::allGeomModules().

1091  {
1094  ? tileCount(0, -1)
1095  : 0;
1096  if (cells == 0) {
1097  unsigned int nlayer = (reco) ? hgpar_->depth_.size() : hgpar_->layer_.size();
1098  for (unsigned k = 0; k < nlayer; ++k) {
1099  std::vector<int> ncells = numberCells(((reco) ? hgpar_->depth_[k] : hgpar_->layer_[k]), reco);
1100  cells = std::accumulate(ncells.begin(), ncells.end(), cells);
1101  }
1102  }
1103  return cells;
1104 }
std::vector< int > layer_
static constexpr int scintillatorCassette
std::vector< int > depth_
static constexpr int scintillatorFile
int tileCount(int layer, int ring) const
bool tileTrapezoid() const
int numberCells(bool reco) const
fixed size matrix
const HGCalParameters * hgpar_

◆ numberCells() [2/2]

std::vector< int > HGCalDDDConstants::numberCells ( int  lay,
bool  reco 
) const

Definition at line 1106 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellFineX_, HGCalParameters::defineFull_, HGCalParameters::firstModule_, getIndex(), HGCSiliconDetId::HGCalFine, hgpar_, mps_fire::i, dqmdumpme::k, nphi, numberCellsHexagon(), HGCalParameters::scintCells(), tileTrapezoid(), HGCalParameters::waferCopy_, waferHexagon6(), waferInLayerTest(), HGCalParameters::waferTypeT_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

1106  {
1107  const auto& index = getIndex(lay, reco);
1108  int i = index.first;
1109  std::vector<int> ncell;
1110  if (i >= 0) {
1111  if (waferHexagon6()) {
1112  for (unsigned int k = 0; k < hgpar_->waferTypeT_.size(); ++k) {
1114  unsigned int cell = (hgpar_->waferTypeT_[k] - 1 == HGCSiliconDetId::HGCalFine)
1115  ? (hgpar_->cellFineX_.size())
1116  : (hgpar_->cellCoarseX_.size());
1117  ncell.emplace_back(static_cast<int>(cell));
1118  }
1119  }
1120  } else if (tileTrapezoid()) {
1121  int nphi = hgpar_->scintCells(lay);
1122  for (int k = hgpar_->firstModule_[i]; k <= hgpar_->lastModule_[i]; ++k)
1123  ncell.emplace_back(nphi);
1124  } else {
1125  for (unsigned int k = 0; k < hgpar_->waferCopy_.size(); ++k) {
1126  if (waferInLayerTest(k, index.first, hgpar_->defineFull_)) {
1127  int cell = numberCellsHexagon(lay,
1130  true);
1131  ncell.emplace_back(cell);
1132  }
1133  }
1134  }
1135  }
1136  return ncell;
1137 }
const int nphi
bool waferInLayerTest(int wafer, int lay, bool full) const
int32_t waferU(const int32_t index)
bool waferHexagon6() const
std::vector< double > cellCoarseX_
std::vector< int > firstModule_
std::pair< int, float > getIndex(int lay, bool reco) const
bool tileTrapezoid() const
int numberCellsHexagon(int wafer) const
std::vector< double > cellFineX_
fixed size matrix
int32_t waferV(const int32_t index)
std::vector< int > waferCopy_
std::vector< int > waferTypeT_
const HGCalParameters * hgpar_
int scintCells(const int layer) const

◆ numberCellsHexagon() [1/2]

int HGCalDDDConstants::numberCellsHexagon ( int  wafer) const

Definition at line 1139 of file HGCalDDDConstants.cc.

References HGCalParameters::cellCoarseX_, HGCalParameters::cellFineX_, HGCSiliconDetId::HGCalFine, hgpar_, and HGCalParameters::waferTypeT_.

Referenced by HGCalGeometry::newCell(), numberCells(), and HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_().

1139  {
1140  if (wafer >= 0 && wafer < static_cast<int>(hgpar_->waferTypeT_.size())) {
1141  if (hgpar_->waferTypeT_[wafer] - 1 == HGCSiliconDetId::HGCalFine)
1142  return static_cast<int>(hgpar_->cellFineX_.size());
1143  else
1144  return static_cast<int>(hgpar_->cellCoarseX_.size());
1145  } else {
1146  return 0;
1147  }
1148 }
std::vector< double > cellCoarseX_
std::vector< double > cellFineX_
std::vector< int > waferTypeT_
const HGCalParameters * hgpar_

◆ numberCellsHexagon() [2/2]

int HGCalDDDConstants::numberCellsHexagon ( int  lay,
int  waferU,
int  waferV,
bool  flag 
) const

Definition at line 1150 of file HGCalDDDConstants.cc.

References RemoveAddSevLevel::flag, HGCSiliconDetId::HGCalCoarseThick, hgpar_, caHitNtupletGeneratorKernels::if(), N, HGCalParameters::nCellsCoarse_, HGCalParameters::nCellsFine_, HGCalParameters::typesInLayers_, HGCalWaferIndex::waferIndex(), HGCalParameters::waferTypeL_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

1150  {
1152  int type =
1153  ((itr == hgpar_->typesInLayers_.end()) ? HGCSiliconDetId::HGCalCoarseThick : hgpar_->waferTypeL_[itr->second]);
1154  int N = (type == 0) ? hgpar_->nCellsFine_ : hgpar_->nCellsCoarse_;
1155  if (flag)
1156  return (3 * N * N);
1157  else
1158  return N;
1159 }
int32_t waferU(const int32_t index)
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
#define N
Definition: blowfish.cc:9
wafer_map typesInLayers_
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_

◆ rangeR()

std::pair< double, double > HGCalDDDConstants::rangeR ( double  z,
bool  reco 
) const

Definition at line 1161 of file HGCalDDDConstants.cc.

References funct::abs(), HGCalParameters::detectorType_, HGCalParameters::firstLayer_, HGCalParameters::firstMixedLayer_, hgpar_, HGCalParameters::k_ScaleFromDDD, HGCalParameters::k_ScaleToDDD, HGCalGeomTools::radius(), HGCalParameters::radiusMixBoundary_, HGCalParameters::rMaxFront_, HGCalParameters::rMinFront_, HGCalParameters::slopeMin_, HGCalParameters::slopeTop_, z, HGCalParameters::zFrontMin_, HGCalParameters::zFrontTop_, HGCalParameters::zLayerHex_, and geometryCSVtoXML::zz.

Referenced by ticl::SeedingRegionByTracks::buildFirstLayers(), ticl::LinkingAlgoByDirectionGeometric::buildLayers(), HFNoseNumberingScheme::checkPosition(), HGCalNumberingScheme::checkPosition(), HGCDigitizer::checkPosition(), locateCellTrap(), and HGCalGeometry::neighborZ().

1161  {
1162  double rmin(0), rmax(0), zz(0);
1163  if (hgpar_->detectorType_ > 0) {
1164  zz = (reco ? std::abs(z) : HGCalParameters::k_ScaleFromDDD * std::abs(z));
1165  if (hgpar_->detectorType_ <= 2) {
1167  } else {
1168  rmin = HGCalGeomTools::radius(
1170  }
1171  if ((hgpar_->detectorType_ == 2) && (zz >= hgpar_->zLayerHex_[hgpar_->firstMixedLayer_ - 1])) {
1172  rmax = HGCalGeomTools::radius(
1174  } else {
1176  }
1177  }
1178  if (!reco) {
1181  }
1182 #ifdef EDM_ML_DEBUG
1183  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:rangeR: " << z << ":" << zz << " R " << rmin << ":" << rmax;
1184 #endif
1185  return std::make_pair(rmin, rmax);
1186 }
Log< level::Info, true > LogVerbatim
std::vector< double > zFrontMin_
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)
std::vector< double > rMaxFront_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< double > slopeTop_
std::vector< double > zLayerHex_
std::vector< double > slopeMin_
static constexpr double k_ScaleToDDD
std::vector< double > radiusMixBoundary_
std::vector< double > rMinFront_
fixed size matrix
std::vector< double > zFrontTop_
const HGCalParameters * hgpar_

◆ rangeRLayer()

std::pair< double, double > HGCalDDDConstants::rangeRLayer ( int  lay,
bool  reco 
) const

Definition at line 1188 of file HGCalDDDConstants.cc.

References getIndex(), hgpar_, HGCalParameters::k_ScaleToDDD, HGCalParameters::rMaxLayHex_, and HGCalParameters::rMinLayHex_.

1188  {
1189  double rmin(0), rmax(0);
1190  const auto& index = getIndex(lay, reco);
1191  if (index.first >= 0 && index.first < static_cast<int>(hgpar_->rMinLayHex_.size())) {
1192  rmin = hgpar_->rMinLayHex_[index.first];
1193  rmax = hgpar_->rMaxLayHex_[index.first];
1194  }
1195  if (!reco) {
1198  }
1199 #ifdef EDM_ML_DEBUG
1200  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:rangeR: " << lay << ":" << index.first << " R " << rmin << ":"
1201  << rmax;
1202 #endif
1203  return std::make_pair(rmin, rmax);
1204 }
Log< level::Info, true > LogVerbatim
std::pair< int, float > getIndex(int lay, bool reco) const
std::vector< double > rMinLayHex_
std::vector< double > rMaxLayHex_
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

◆ rangeZ()

std::pair< double, double > HGCalDDDConstants::rangeZ ( bool  reco) const

Definition at line 1206 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::k_ScaleToDDD, HGCalParameters::waferThick_, HGCalParameters::zLayerHex_, SiStripMonitorCluster_cfi::zmax, and SiStripMonitorCluster_cfi::zmin.

Referenced by HFNoseNumberingScheme::checkPosition(), HGCalNumberingScheme::checkPosition(), and HGCDigitizer::checkPosition().

1206  {
1207  double zmin = (hgpar_->zLayerHex_[0] - hgpar_->waferThick_);
1208  double zmax = (hgpar_->zLayerHex_[hgpar_->zLayerHex_.size() - 1] + hgpar_->waferThick_);
1209 #ifdef EDM_ML_DEBUG
1210  edm::LogVerbatim("HGCalGeom") << "HGCalDDDConstants:rangeZ: " << zmin << ":" << zmax << ":" << hgpar_->waferThick_;
1211 #endif
1212  if (!reco) {
1215  }
1216  return std::make_pair(zmin, zmax);
1217 }
Log< level::Info, true > LogVerbatim
std::vector< double > zLayerHex_
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

◆ rowColumnWafer()

std::pair< int, int > HGCalDDDConstants::rowColumnWafer ( const int  wafer) const

Definition at line 1219 of file HGCalDDDConstants.cc.

References cuy::col, filterCSVwithJSON::copy, HGCalTypes::getUnpackedU(), HGCalTypes::getUnpackedV(), hgpar_, and HGCalParameters::waferCopy_.

1219  {
1220  int row(0), col(0);
1221  if (wafer < static_cast<int>(hgpar_->waferCopy_.size())) {
1222  int copy = hgpar_->waferCopy_[wafer];
1225  ;
1226  }
1227  return std::make_pair(row, col);
1228 }
static int32_t getUnpackedU(int id)
Definition: HGCalTypes.cc:16
static int32_t getUnpackedV(int id)
Definition: HGCalTypes.cc:22
std::vector< int > waferCopy_
col
Definition: cuy.py:1009
const HGCalParameters * hgpar_

◆ sectors()

int HGCalDDDConstants::sectors ( ) const
inline

◆ simToReco()

std::pair< int, int > HGCalDDDConstants::simToReco ( int  cell,
int  layer,
int  mod,
bool  half 
) const

Definition at line 1230 of file HGCalDDDConstants.cc.

References LEDCalibrationChannels::depth, getIndex(), hgpar_, mps_fire::i, HGCalParameters::layerGroup_, HGCalParameters::layerGroupM_, HGCalParameters::layerGroupO_, mod(), findQualityFiles::size, waferHexagon6(), and HGCalParameters::waferTypeL_.

Referenced by HGCalTimingAnalyzer::analyzeSimHits(), HGCalTBAnalyzer::analyzeSimHits(), CaloParticleDebugger::fillSimHits(), and CaloTruthAccumulator::fillSimHits().

1230  {
1231  if (!waferHexagon6()) {
1232  return std::make_pair(cell, lay);
1233  } else {
1234  const auto& index = getIndex(lay, false);
1235  int i = index.first;
1236  if (i < 0) {
1237  edm::LogWarning("HGCalGeom") << "Wrong Layer # " << lay << " not in the list ***** ERROR *****";
1238  return std::make_pair(-1, -1);
1239  }
1240  if (mod >= static_cast<int>(hgpar_->waferTypeL_.size())) {
1241  edm::LogWarning("HGCalGeom") << "Invalid Wafer # " << mod << "should be < " << (hgpar_->waferTypeL_).size()
1242  << " ***** ERROR *****";
1243  return std::make_pair(-1, -1);
1244  }
1245  int depth(-1);
1246  int kx = cell;
1247  int type = hgpar_->waferTypeL_[mod];
1248  if (type == 1) {
1249  depth = hgpar_->layerGroup_[i];
1250  } else if (type == 2) {
1251  depth = hgpar_->layerGroupM_[i];
1252  } else {
1253  depth = hgpar_->layerGroupO_[i];
1254  }
1255  return std::make_pair(kx, depth);
1256  }
1257 }
size
Write out results.
std::vector< int > layerGroupM_
bool waferHexagon6() const
std::pair< int, float > getIndex(int lay, bool reco) const
std::vector< int > layerGroup_
std::vector< int > layerGroupO_
Log< level::Warning, false > LogWarning
const HGCalParameters * hgpar_
T mod(const T &a, const T &b)
Definition: ecalDccMap.h:4
std::vector< int > waferTypeL_

◆ tileCount()

int HGCalDDDConstants::tileCount ( int  layer,
int  ring 
) const

Definition at line 1259 of file HGCalDDDConstants.cc.

References b, HGCalParameters::firstLayer_, hgpar_, dqmdumpme::k, lastLayer(), pixelTopology::layer, relativeConstraints::ring, HGCalTileIndex::tileIndex(), HGCalParameters::tileInfoMap_, and HGCalParameters::tileRingRange_.

Referenced by numberCells().

1259  {
1260  int laymin(layer), laymax(layer), ringmin(ring), ringmax(ring), kount(0);
1261  if (layer == 0) {
1262  laymin = hgpar_->firstLayer_;
1263  laymax = lastLayer(true);
1264  }
1265 #ifdef EDM_ML_DEBUG
1266  edm::LogVerbatim("HGCalGeom") << "tileCount: layer " << layer << " ring " << ring << " layerMin/Max " << laymin << ":"
1267  << laymax;
1268 #endif
1269  for (int lay = laymin; lay <= laymax; ++lay) {
1270  if (ring < 0) {
1271  int ll = lay - hgpar_->firstLayer_;
1272  ringmin = hgpar_->tileRingRange_[ll].first;
1273  ringmax = hgpar_->tileRingRange_[ll].second;
1274  }
1275 #ifdef EDM_ML_DEBUG
1276  edm::LogVerbatim("HGCalGeom") << "tileCount: lay " << lay << ":" << (lay - hgpar_->firstLayer_) << " rings "
1277  << ringmin << ":" << ringmax;
1278 #endif
1279  for (int rin = ringmin; rin <= ringmax; ++rin) {
1280  int indx = HGCalTileIndex::tileIndex(lay, rin + 1, 0);
1281  auto itr = hgpar_->tileInfoMap_.find(indx);
1282 #ifdef EDM_ML_DEBUG
1283  edm::LogVerbatim("HGCalGeom") << "tileCount: rin " << rin << " indx " << indx << " itr "
1284  << (itr != hgpar_->tileInfoMap_.end());
1285 #endif
1286  if (itr != hgpar_->tileInfoMap_.end()) {
1287  for (int k = 0; k < 4; ++k) {
1288  std::bitset<24> b(itr->second.hex[k]);
1289  kount += b.count();
1290  }
1291  }
1292 #ifdef EDM_ML_DEBUG
1293  edm::LogVerbatim("HGCalGeom") << "tileCount: lay|rin " << lay << ":" << rin << " kount " << kount;
1294 #endif
1295  }
1296  }
1297  return (3 * kount);
1298 }
Log< level::Info, true > LogVerbatim
int lastLayer(bool reco) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t tileIndex(int32_t layer, int32_t ring, int32_t phi)
std::vector< std::pair< int, int > > tileRingRange_
double b
Definition: hdecay.h:118
const HGCalParameters * hgpar_
tileInfo_map tileInfoMap_

◆ tileExist()

bool HGCalDDDConstants::tileExist ( int  zside,
int  layer,
int  ring,
int  phi 
) const

Definition at line 1300 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, convertSQLiteXML::ok, phi, relativeConstraints::ring, HGCalTileIndex::tileExist(), HGCalTileIndex::tileIndex(), HGCalParameters::tileInfoMap_, and ecaldqm::zside().

Referenced by assignCellTrap(), HGCalGeometryLoader::build(), and isValidTrap().

1300  {
1301  int indx = HGCalTileIndex::tileIndex(layer, ring, 0);
1302  auto itr = hgpar_->tileInfoMap_.find(indx);
1303  bool ok = (itr == hgpar_->tileInfoMap_.end()) ? false : HGCalTileIndex::tileExist(itr->second.hex, zside, phi);
1304  return ok;
1305 }
int zside(DetId const &)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t tileIndex(int32_t layer, int32_t ring, int32_t phi)
bool tileExist(const int32_t *hex, int32_t zside, int32_t phi)
const HGCalParameters * hgpar_
tileInfo_map tileInfoMap_

◆ tileInfo()

HGCalParameters::tileInfo HGCalDDDConstants::tileInfo ( int  zside,
int  layer,
int  ring 
) const

Definition at line 1307 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, relativeConstraints::ring, HGCalTileIndex::tileIndex(), and HGCalParameters::tileInfoMap_.

1307  {
1308  int indx = HGCalTileIndex::tileIndex(layer, ring, 0);
1309  auto itr = hgpar_->tileInfoMap_.find(indx);
1310  return ((itr == hgpar_->tileInfoMap_.end()) ? HGCalParameters::tileInfo() : itr->second);
1311 }
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t tileIndex(int32_t layer, int32_t ring, int32_t phi)
const HGCalParameters * hgpar_
tileInfo_map tileInfoMap_

◆ tilePhiEdge()

bool HGCalDDDConstants::tilePhiEdge ( double  phi,
int  layer,
int  iphi 
) const

Definition at line 1313 of file HGCalDDDConstants.cc.

References funct::abs(), cuy::dif1, cuy::dif2, hgpar_, LEDCalibrationChannels::iphi, pixelTopology::layer, phi, HGCalParameters::scintCellSize(), and tol_.

Referenced by assignCellTrap().

1313  {
1314  double dif1 = std::abs(phi - hgpar_->scintCellSize(layer) * (iphi - 1));
1315  double dif2 = std::abs(phi - hgpar_->scintCellSize(layer) * iphi);
1316 #ifdef EDM_ML_DEBUG
1317  edm::LogVerbatim("HGCalGeomT") << "HGCalDDDConstants::tilePhiEdge:: input: " << phi << ":" << layer << ":" << iphi
1318  << " Differences " << dif1 << ":" << dif2;
1319 #endif
1320  return ((dif1 < tol_) || (dif2 < tol_));
1321 }
Log< level::Info, true > LogVerbatim
static constexpr double tol_
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
double scintCellSize(const int layer) const
dif1
Definition: cuy.py:893
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
dif2
Definition: cuy.py:895
const HGCalParameters * hgpar_

◆ tileRingEdge()

bool HGCalDDDConstants::tileRingEdge ( double  rho,
int  layer,
int  ring 
) const

Definition at line 1323 of file HGCalDDDConstants.cc.

References funct::abs(), cuy::dif1, cuy::dif2, hgpar_, pixelTopology::layer, alignCSCRings::r, HGCalParameters::radiusLayer_, relativeConstraints::ring, HGCalParameters::scintType(), and tol_.

Referenced by assignCellTrap().

1323  {
1324  int type = hgpar_->scintType(layer);
1325  double dif1 = std::abs(r - hgpar_->radiusLayer_[type][ring - 1]);
1326  double dif2 = std::abs(r - hgpar_->radiusLayer_[type][ring]);
1327 #ifdef EDM_ML_DEBUG
1328  edm::LogVerbatim("HGCalGeomT") << "HGCalDDDConstants::tileRingEdge:: input: " << r << ":" << layer << ":" << ring
1329  << " Differences " << dif1 << ":" << dif2;
1330 #endif
1331  return ((dif1 < tol_) || (dif2 < tol_));
1332 }
Log< level::Info, true > LogVerbatim
static constexpr double tol_
int scintType(const int layer) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
dif1
Definition: cuy.py:893
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
dif2
Definition: cuy.py:895
std::vector< double > radiusLayer_[2]
const HGCalParameters * hgpar_

◆ tileRings()

std::pair< int, int > HGCalDDDConstants::tileRings ( int  layer) const

Definition at line 1333 of file HGCalDDDConstants.cc.

References HGCalParameters::firstLayer_, hgpar_, pixelTopology::layer, mode_, HGCalParameters::tileRingRange_, HGCalGeometryMode::TrapezoidCassette, HGCalGeometryMode::TrapezoidFile, and HGCalGeometryMode::TrapezoidModule.

1333  {
1336  int ll = layer - hgpar_->firstLayer_;
1337  if (ll >= 0 && ll < static_cast<int>(hgpar_->tileRingRange_.size()))
1338  return hgpar_->tileRingRange_[ll];
1339  }
1340  return std::make_pair(0, 0);
1341 }
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::vector< std::pair< int, int > > tileRingRange_
const HGCalGeometryMode::GeometryMode mode_
const HGCalParameters * hgpar_

◆ tileSiPM()

int HGCalDDDConstants::tileSiPM ( int  sipm) const
inline

◆ tileTrapezoid()

bool HGCalDDDConstants::tileTrapezoid ( ) const
inline

◆ tileType()

std::pair< int, int > HGCalDDDConstants::tileType ( int  layer,
int  ring,
int  phi 
) const

Definition at line 1343 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, phi, relativeConstraints::ring, HGCalTypes::SiPMLarge, HGCalTileIndex::tileIndex(), and HGCalParameters::tileInfoMap_.

Referenced by HGCalGeometryLoader::build(), HGCalTopology::encode(), HGCalNumberingScheme::getUnitID(), and HGCalGeometry::newCell().

1343  {
1344  int indx = HGCalTileIndex::tileIndex(layer, ring, phi);
1345  int type(-1), sipm(-1);
1346  auto itr = hgpar_->tileInfoMap_.find(indx);
1347  if (itr != hgpar_->tileInfoMap_.end()) {
1348  type = 1 + (itr->second).type;
1349  sipm = ((itr->second).sipm == HGCalTypes::SiPMLarge) ? 0 : 1;
1350  }
1351  return std::make_pair(type, sipm);
1352 }
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t tileIndex(int32_t layer, int32_t ring, int32_t phi)
const HGCalParameters * hgpar_
tileInfo_map tileInfoMap_

◆ trapezoidFile()

bool HGCalDDDConstants::trapezoidFile ( ) const
inline

◆ volumes()

unsigned int HGCalDDDConstants::volumes ( ) const
inline

Definition at line 147 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::moduleLayR_.

Referenced by HGCalGeometryLoader::build().

147 { return hgpar_->moduleLayR_.size(); }
std::vector< int > moduleLayR_
const HGCalParameters * hgpar_

◆ waferCount()

int HGCalDDDConstants::waferCount ( const int  type) const
inline

Definition at line 176 of file HGCalDDDConstants.h.

References waferMax_.

176 { return ((type == 0) ? waferMax_[2] : waferMax_[3]); }
std::array< int, 4 > waferMax_

◆ waferFileIndex()

int HGCalDDDConstants::waferFileIndex ( unsigned int  kk) const

Definition at line 1603 of file HGCalDDDConstants.cc.

References hgpar_, GetRecoTauVFromDQM_MC_cff::kk, and HGCalParameters::waferInfoMap_.

1603  {
1604  if (kk < hgpar_->waferInfoMap_.size()) {
1605  auto itr = hgpar_->waferInfoMap_.begin();
1606  std::advance(itr, kk);
1607  return itr->first;
1608  } else
1609  return 0;
1610 }
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferFileInfo()

std::tuple< int, int, int > HGCalDDDConstants::waferFileInfo ( unsigned int  kk) const

Definition at line 1612 of file HGCalDDDConstants.cc.

References hgpar_, GetRecoTauVFromDQM_MC_cff::kk, and HGCalParameters::waferInfoMap_.

1612  {
1613  if (kk < hgpar_->waferInfoMap_.size()) {
1614  auto itr = hgpar_->waferInfoMap_.begin();
1615  std::advance(itr, kk);
1616  return std::make_tuple(itr->second.type, itr->second.part, itr->second.orient);
1617  } else
1618  return std::make_tuple(0, 0, 0);
1619 }
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferFileInfoExist()

bool HGCalDDDConstants::waferFileInfoExist ( int  kk) const
inline

Definition at line 187 of file HGCalDDDConstants.h.

References hgpar_, GetRecoTauVFromDQM_MC_cff::kk, and HGCalParameters::waferInfoMap_.

187  {
188  return (hgpar_->waferInfoMap_.find(kk) != hgpar_->waferInfoMap_.end());
189  }
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferFileInfoFromIndex()

std::tuple< int, int, int > HGCalDDDConstants::waferFileInfoFromIndex ( int  kk) const

Definition at line 1621 of file HGCalDDDConstants.cc.

References hgpar_, GetRecoTauVFromDQM_MC_cff::kk, and HGCalParameters::waferInfoMap_.

1621  {
1622  auto itr = hgpar_->waferInfoMap_.find(kk);
1623  if (itr != hgpar_->waferInfoMap_.end()) {
1624  return std::make_tuple(itr->second.type, itr->second.part, itr->second.orient);
1625  } else
1626  return std::make_tuple(0, 0, 0);
1627 }
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferFileSize()

unsigned int HGCalDDDConstants::waferFileSize ( ) const
inline

Definition at line 183 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::waferInfoMap_.

Referenced by waferType(), and waferTypeRotation().

183 { return hgpar_->waferInfoMap_.size(); }
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferFromCopy()

int HGCalDDDConstants::waferFromCopy ( int  copy) const

Definition at line 1354 of file HGCalDDDConstants.cc.

References filterCSVwithJSON::copy, hgpar_, dqmdumpme::k, mps_fire::result, and HGCalParameters::waferCopy_.

Referenced by HGCNumberingScheme::getUnitID().

1354  {
1355  const int ncopies = hgpar_->waferCopy_.size();
1356  int wafer(ncopies);
1357  bool result(false);
1358  for (int k = 0; k < ncopies; ++k) {
1359  if (copy == hgpar_->waferCopy_[k]) {
1360  wafer = k;
1361  result = true;
1362  break;
1363  }
1364  }
1365  if (!result) {
1366  wafer = -1;
1367 #ifdef EDM_ML_DEBUG
1368  edm::LogVerbatim("HGCalGeom") << "Cannot find " << copy << " in a list of " << ncopies << " members";
1369  for (int k = 0; k < ncopies; ++k)
1370  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << hgpar_->waferCopy_[k];
1371 #endif
1372  }
1373 #ifdef EDM_ML_DEBUG
1374  edm::LogVerbatim("HGCalGeom") << "WaferFromCopy " << copy << ":" << wafer << ":" << result;
1375 #endif
1376  return wafer;
1377 }
Log< level::Info, true > LogVerbatim
std::vector< int > waferCopy_
const HGCalParameters * hgpar_

◆ waferFromPosition() [1/2]

void HGCalDDDConstants::waferFromPosition ( const double  x,
const double  y,
int &  wafer,
int &  icell,
int &  celltyp 
) const

Definition at line 1379 of file HGCalDDDConstants.cc.

References funct::abs(), HGCalParameters::cellCoarseX_, HGCalParameters::cellCoarseY_, HGCalParameters::cellFineX_, HGCalParameters::cellFineY_, cellHex(), HGCalParameters::cellSize_, PVValHelper::dx, PVValHelper::dy, hexside_, HGCSiliconDetId::HGCalFine, hgpar_, dqmdumpme::k, HGCalParameters::k_ScaleFromDDD, rmax_, tan30deg_, HGCalParameters::waferCopy_, HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::waferTypeT_, x, geometryCSVtoXML::xx, y, and geometryCSVtoXML::yy.

Referenced by assignCellHex(), HFNoseNumberingScheme::checkPosition(), HGCalNumberingScheme::checkPosition(), HFNoseNumberingScheme::getUnitID(), HGCNumberingScheme::getUnitID(), and HGCalNumberingScheme::getUnitID().

1379  {
1380  // Input x, y in Geant4 unit and transformed to CMSSW standard
1383  int size_ = static_cast<int>(hgpar_->waferCopy_.size());
1384  wafer = size_;
1385  for (int k = 0; k < size_; ++k) {
1386  double dx = std::abs(xx - hgpar_->waferPosX_[k]);
1387  double dy = std::abs(yy - hgpar_->waferPosY_[k]);
1388  if (dx <= rmax_ && dy <= hexside_) {
1389  if ((dy <= 0.5 * hexside_) || (dx * tan30deg_ <= (hexside_ - dy))) {
1390  wafer = k;
1391  celltyp = hgpar_->waferTypeT_[k];
1392  xx -= hgpar_->waferPosX_[k];
1393  yy -= hgpar_->waferPosY_[k];
1394  break;
1395  }
1396  }
1397  }
1398  if (wafer < size_) {
1399  if (celltyp - 1 == HGCSiliconDetId::HGCalFine)
1400  icell = cellHex(
1402  else
1403  icell = cellHex(xx,
1404  yy,
1407  hgpar_->cellCoarseY_);
1408  } else {
1409  wafer = -1;
1410 #ifdef EDM_ML_DEBUG
1411  edm::LogWarning("HGCalGeom") << "Cannot get wafer type corresponding to " << x << ":" << y << " " << xx << ":"
1412  << yy;
1413 #endif
1414  }
1415 #ifdef EDM_ML_DEBUG
1416  edm::LogVerbatim("HGCalGeom") << "Position " << x << ":" << y << " Wafer " << wafer << ":" << size_ << " XX " << xx
1417  << ":" << yy << " Cell " << icell << " Type " << celltyp;
1418 #endif
1419 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
std::vector< double > cellFineY_
std::vector< double > cellCoarseX_
std::vector< double > cellSize_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< double > cellFineX_
static constexpr double k_ScaleFromDDD
std::vector< int > waferCopy_
static constexpr double tan30deg_
std::vector< int > waferTypeT_
std::vector< double > cellCoarseY_
Log< level::Warning, false > LogWarning
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
int cellHex(double xx, double yy, const double &cellR, const std::vector< double > &posX, const std::vector< double > &posY) const

◆ waferFromPosition() [2/2]

void HGCalDDDConstants::waferFromPosition ( const double  x,
const double  y,
const int  zside,
const int  layer,
int &  waferU,
int &  waferV,
int &  cellU,
int &  cellV,
int &  celltype,
double &  wt,
bool  extend,
bool  debug 
) const

Definition at line 1421 of file HGCalDDDConstants.cc.

References funct::abs(), cellHex(), HGCalCell::cellPlacementIndex(), HGCalCell::cellPlacementOld, HGCalParameters::cellThickness_, debug, PVValHelper::dx, PVValHelper::dy, HGCalParameters::firstLayer_, HGCalWaferType::getCassette(), HGCalWaferType::getOrient(), HGCalWaferType::getPartial(), HGCalCassette::getShift(), HGCalWaferType::getType(), getXY(), HGCalGeometryMode::Hexagon8Cassette, hexside_, hexsideT_, HGCSiliconDetId::HGCalCoarseThick, hgcassette_, hgpar_, dqmdumpme::k, HGCalParameters::k_ScaleFromDDD, pixelTopology::layer, HGCalTypes::layerFrontBack(), layerType(), HGCalParameters::layerType_, mode_, HGCalParameters::nCellsFine_, hgcalPerformanceValidation::orient, rmax_, rmaxT_, tan30deg_, HGCalParameters::typesInLayers_, HGCalParameters::useSimWt_, HGCalTypes::WaferCenterR, HGCalParameters::waferCopy_, HGCalTypes::WaferFull, waferHexagon8File(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::waferThick_, HGCalParameters::waferTypeL_, HGCalWaferIndex::waferU(), HGCalParameters::waferUVMax_, HGCalWaferIndex::waferV(), x, testProducerWithPsetDescEmpty_cfi::x1, HGCalParameters::xLayerHex_, geometryCSVtoXML::xx, geometryCSVtoXML::xy, y, testProducerWithPsetDescEmpty_cfi::y1, HGCalParameters::yLayerHex_, geometryCSVtoXML::yy, and ecaldqm::zside().

1432  {
1433  // Expect x, y as in SIM step
1434  waferU = waferV = 1 + hgpar_->waferUVMax_;
1435  cellU = cellV = celltype = 0;
1436  if ((hgpar_->xLayerHex_.empty()) || (hgpar_->yLayerHex_.empty()))
1437  return;
1438  int ll = layer - hgpar_->firstLayer_;
1439  int layertype = layerType(layer);
1440  bool rotx = ((!hgpar_->layerType_.empty()) && (layertype == HGCalTypes::WaferCenterR));
1441  double xx(0), yy(0);
1442  if (rotx) {
1443  std::pair<double, double> xy =
1445  xx = xy.first - hgpar_->xLayerHex_[ll];
1446  yy = xy.second - hgpar_->yLayerHex_[ll];
1447  } else {
1450  }
1451  if (debug)
1452  edm::LogVerbatim("HGCalGeom") << "waferFromPosition:: Layer " << layer << ":" << ll << " Rot " << rotx << " X " << x
1453  << ":" << xx << " Y " << y << ":" << yy;
1454  double rmax = extend ? rmaxT_ : rmax_;
1455  double hexside = extend ? hexsideT_ : hexside_;
1456  for (unsigned int k = 0; k < hgpar_->waferPosX_.size(); ++k) {
1457  double dx0(0), dy0(0);
1462  auto ktr = hgpar_->waferInfoMap_.find(indx);
1463  if (ktr != hgpar_->waferInfoMap_.end()) {
1464  auto cshift = hgcassette_.getShift(layer, -zside, (ktr->second).cassette);
1465  if (debug)
1466  edm::LogVerbatim("HGCalGeom") << "Cassette " << (ktr->second).cassette << " Shift " << -(zside * cshift.first)
1467  << ":" << cshift.second;
1468  dx0 = -(zside * cshift.first);
1469  dy0 = cshift.second;
1470  }
1471  }
1472  double dx = std::abs(xx - dx0 - hgpar_->waferPosX_[k]);
1473  double dy = std::abs(yy - dy0 - hgpar_->waferPosY_[k]);
1474  if (dx <= rmax && dy <= hexside) {
1475  if ((dy <= 0.5 * hexside) || (dx * tan30deg_ <= (hexside - dy))) {
1476  if (waferHexagon8File()) {
1479  if (debug)
1480  edm::LogVerbatim("HGCalGeom") << "Position (" << x << ", " << y << ") Wafer type:partial:orient:cassette "
1481  << celltype << ":" << HGCalWaferType::getPartial(index, hgpar_->waferInfoMap_)
1484  } else {
1486  celltype = ((itr == hgpar_->typesInLayers_.end()) ? HGCSiliconDetId::HGCalCoarseThick
1487  : hgpar_->waferTypeL_[itr->second]);
1488  }
1489  if (debug)
1490  edm::LogVerbatim("HGCalGeom") << "WaferFromPosition:: Input " << layer << ":" << ll << ":"
1491  << hgpar_->firstLayer_ << ":" << rotx << ":" << x << ":" << y << ":"
1492  << hgpar_->xLayerHex_[ll] << ":" << hgpar_->yLayerHex_[ll] << ":" << xx << ":"
1493  << yy << " compared with " << hgpar_->waferPosX_[k] << ":"
1494  << hgpar_->waferPosY_[k] << " difference " << dx << ":" << dy << ":"
1495  << dx * tan30deg_ << ":" << (hexside_ - dy) << " comparator " << rmax_ << ":"
1496  << rmaxT_ << ":" << hexside_ << ":" << hexsideT_ << " wafer " << waferU << ":"
1497  << waferV << ":" << celltype;
1498  xx -= (dx0 + hgpar_->waferPosX_[k]);
1499  yy -= (dy0 + hgpar_->waferPosY_[k]);
1500  break;
1501  }
1502  }
1503  }
1504  if ((std::abs(waferU) <= hgpar_->waferUVMax_) && (celltype >= 0)) {
1508  auto ktr = hgpar_->waferInfoMap_.find(indx);
1509  if (ktr != hgpar_->waferInfoMap_.end()) {
1510  place = HGCalCell::cellPlacementIndex(1, HGCalTypes::layerFrontBack(layertype), (ktr->second).orient);
1511  part = (ktr->second).part;
1512  if (debug)
1513  edm::LogVerbatim("HGCalGeom") << "waferFromPosition: frontback " << layertype << ":"
1514  << HGCalTypes::layerFrontBack(layertype) << " Orient " << (ktr->second).orient
1515  << " place " << place << " part " << part;
1516  }
1517  }
1518  cellHex(xx, yy, celltype, place, part, cellU, cellV, extend, debug);
1519  wt = (((celltype < 2) && (hgpar_->useSimWt_ > 0)) ? (hgpar_->cellThickness_[celltype] / hgpar_->waferThick_) : 1.0);
1520  } else {
1521  cellU = cellV = 2 * hgpar_->nCellsFine_;
1522  wt = 1.0;
1523  celltype = -1;
1524  }
1525  if ((celltype < 0) && debug) {
1526  double x1(xx);
1527  double y1(yy);
1528  edm::LogVerbatim("HGCalGeom") << "waferfFromPosition: Bad type for X " << x << ":" << x1 << ":" << xx << " Y " << y
1529  << ":" << y1 << ":" << yy << " Wafer " << waferU << ":" << waferV << " Cell " << cellU
1530  << ":" << cellV;
1531  for (unsigned int k = 0; k < hgpar_->waferPosX_.size(); ++k) {
1532  double dx = std::abs(x1 - hgpar_->waferPosX_[k]);
1533  double dy = std::abs(y1 - hgpar_->waferPosY_[k]);
1534  edm::LogVerbatim("HGCalGeom") << "Wafer [" << k << "] Position (" << hgpar_->waferPosX_[k] << ", "
1535  << hgpar_->waferPosY_[k] << ") difference " << dx << ":" << dy << ":"
1536  << dx * tan30deg_ << ":" << hexside - dy << " Paramerers " << rmax << ":"
1537  << hexside;
1538  }
1539  }
1540  edm::LogVerbatim("HGCalGeomX") << "Input x:y:layer " << x << ":" << y << ":" << layer << " Wafer " << waferU << ":"
1541  << waferV << " Cell " << cellU << ":" << cellV << ":" << celltype << " wt " << wt;
1542 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:233
static constexpr int32_t cellPlacementOld
Definition: HGCalCell.h:25
static int getType(int index, const HGCalParameters::waferInfo_map &wafers)
HGCalCassette hgcassette_
int32_t waferU(const int32_t index)
std::pair< double, double > getShift(int layer, int zside, int cassette) const
static int getPartial(int index, const HGCalParameters::waferInfo_map &wafers)
std::pair< double, double > getXY(int layer, double x, double y, bool forwd) const
int zside(DetId const &)
static int getOrient(int index, const HGCalParameters::waferInfo_map &wafers)
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
std::vector< double > yLayerHex_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
static int getCassette(int index, const HGCalParameters::waferInfo_map &wafers)
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
std::vector< int > layerType_
static constexpr int32_t WaferCenterR
Definition: HGCalTypes.h:27
#define debug
Definition: HDRShower.cc:19
std::vector< double > cellThickness_
part
Definition: HCALResponse.h:20
bool waferHexagon8File() const
wafer_map typesInLayers_
const HGCalGeometryMode::GeometryMode mode_
static constexpr double k_ScaleFromDDD
int32_t waferV(const int32_t index)
std::vector< int > waferCopy_
static constexpr double tan30deg_
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
int cellHex(double xx, double yy, const double &cellR, const std::vector< double > &posX, const std::vector< double > &posY) const
std::vector< int > waferTypeL_
std::vector< double > xLayerHex_
static int32_t layerFrontBack(int32_t layerOrient)
Definition: HGCalTypes.h:125
int layerType(int lay) const

◆ waferFullInLayer()

bool HGCalDDDConstants::waferFullInLayer ( int  wafer,
int  lay,
bool  reco 
) const

Definition at line 1551 of file HGCalDDDConstants.cc.

References getIndex(), and waferInLayerTest().

1551  {
1552  const auto& indx = getIndex(lay, reco);
1553  if (indx.first < 0)
1554  return false;
1555  return waferInLayerTest(wafer, indx.first, false);
1556 }
bool waferInLayerTest(int wafer, int lay, bool full) const
std::pair< int, float > getIndex(int lay, bool reco) const
fixed size matrix

◆ waferHexagon6()

bool HGCalDDDConstants::waferHexagon6 ( ) const
inline

◆ waferHexagon8()

bool HGCalDDDConstants::waferHexagon8 ( ) const
inline

◆ waferHexagon8File()

bool HGCalDDDConstants::waferHexagon8File ( ) const
inline

◆ waferIndex()

int32_t HGCalDDDConstants::waferIndex ( int  wafer,
int  index 
) const
private

Definition at line 2034 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, layerFromIndex(), HGCalParameters::waferCopy_, HGCalWaferIndex::waferIndex(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

Referenced by waferInLayerTest().

2034  {
2035  int layer = layerFromIndex(index, true);
2039 #ifdef EDM_ML_DEBUG
2040  edm::LogVerbatim("HGCalGeom") << "WaferIndex for " << wafer << ":" << index << " (" << layer << ":" << waferU << ":"
2041  << waferV << ") " << indx;
2042 #endif
2043  return indx;
2044 }
Log< level::Info, true > LogVerbatim
int32_t waferU(const int32_t index)
int layerFromIndex(int index, bool reco) const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
int32_t waferV(const int32_t index)
std::vector< int > waferCopy_
const HGCalParameters * hgpar_

◆ waferInfo()

HGCalParameters::waferInfo HGCalDDDConstants::waferInfo ( int  lay,
int  waferU,
int  waferV 
) const

Definition at line 1558 of file HGCalDDDConstants.cc.

References hgpar_, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

Referenced by HGCalNumberingScheme::getUnitID().

1558  {
1559  int indx = HGCalWaferIndex::waferIndex(lay, waferU, waferV);
1560  auto itr = hgpar_->waferInfoMap_.find(indx);
1561  return ((itr == hgpar_->waferInfoMap_.end()) ? HGCalParameters::waferInfo() : itr->second);
1562 }
int32_t waferU(const int32_t index)
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
int32_t waferV(const int32_t index)
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferInLayer()

bool HGCalDDDConstants::waferInLayer ( int  wafer,
int  lay,
bool  reco 
) const

Definition at line 1544 of file HGCalDDDConstants.cc.

References HGCalParameters::defineFull_, getIndex(), hgpar_, and waferInLayerTest().

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getSummary(), and HGCalDDDConstants().

1544  {
1545  const auto& indx = getIndex(lay, reco);
1546  if (indx.first < 0)
1547  return false;
1548  return waferInLayerTest(wafer, indx.first, hgpar_->defineFull_);
1549 }
bool waferInLayerTest(int wafer, int lay, bool full) const
std::pair< int, float > getIndex(int lay, bool reco) const
fixed size matrix
const HGCalParameters * hgpar_

◆ waferInLayerTest()

bool HGCalDDDConstants::waferInLayerTest ( int  wafer,
int  lay,
bool  full 
) const
private

Definition at line 2046 of file HGCalDDDConstants.cc.

References TCMET_cfi::corner, full, fullAndPart_, hexside_, hgpar_, recoMuon::in, HGCalParameters::k_CornerSize, rmax_, HGCalParameters::rMaxLayHex_, HGCalParameters::rMinLayHex_, HGCalGeomTools::waferCorner(), waferHexagon6(), waferIndex(), HGCalParameters::waferInfoMap_, HGCalWaferIndex::waferLayer(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), HGCalParameters::xLayerHex_, and HGCalParameters::yLayerHex_.

Referenced by maxCells(), maxRows(), modulesInit(), numberCells(), waferFullInLayer(), and waferInLayer().

2046  {
2047  bool in = (waferHexagon6()) ? true : false;
2048  if (!in) {
2049  double xpos = hgpar_->waferPosX_[wafer] + hgpar_->xLayerHex_[lay];
2050  double ypos = hgpar_->waferPosY_[wafer] + hgpar_->yLayerHex_[lay];
2051  std::pair<int, int> corner = HGCalGeomTools::waferCorner(
2052  xpos, ypos, rmax_, hexside_, hgpar_->rMinLayHex_[lay], hgpar_->rMaxLayHex_[lay], in);
2053  in = (full ? (corner.first > 0) : (corner.first == static_cast<int>(HGCalParameters::k_CornerSize)));
2054  if (in && fullAndPart_) {
2055  int indx = waferIndex(wafer, lay);
2056  in = (hgpar_->waferInfoMap_.find(indx) != hgpar_->waferInfoMap_.end());
2057 #ifdef EDM_ML_DEBUG
2058  if (!in)
2059  edm::LogVerbatim("HGCalGeom") << "WaferInLayerTest: Layer " << lay << " wafer " << wafer << " index " << indx
2060  << "( " << HGCalWaferIndex::waferLayer(indx) << ", "
2061  << HGCalWaferIndex::waferU(indx) << ", " << HGCalWaferIndex::waferV(indx)
2062  << ") in " << in;
2063 #endif
2064  }
2065 #ifdef EDM_ML_DEBUG
2066  edm::LogVerbatim("HGCalGeom") << "WaferInLayerTest: Layer " << lay << " wafer " << wafer << " R-limits "
2067  << hgpar_->rMinLayHex_[lay] << ":" << hgpar_->rMaxLayHex_[lay] << " Corners "
2068  << corner.first << ":" << corner.second << " In " << in;
2069 #endif
2070  }
2071  return in;
2072 }
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
bool waferHexagon6() const
static constexpr uint32_t k_CornerSize
static std::pair< int32_t, int32_t > waferCorner(double xpos, double ypos, double r, double R, double rMin, double rMax, bool oldBug=false)
std::vector< double > yLayerHex_
Definition: GenABIO.cc:168
std::vector< double > rMinLayHex_
std::vector< double > rMaxLayHex_
int32_t waferIndex(int wafer, int index) const
int32_t waferV(const int32_t index)
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
std::vector< double > waferPosX_
std::vector< double > xLayerHex_

◆ waferLocal2Global()

GlobalPoint HGCalDDDConstants::waferLocal2Global ( HepGeom::Point3D< float > &  loc,
const DetId id,
bool  useWafer,
bool  reco,
bool  debug 
) const

Definition at line 1629 of file HGCalDDDConstants.cc.

References debug, getXY(), HGCSiliconDetId::layer(), dt_dqm_sourceclient_common_cff::reco, waferPositionNoRot(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferV(), waferZ(), x, geometryCSVtoXML::xx, geometryCSVtoXML::xy, y, HGCSiliconDetId::zside(), and geometryCSVtoXML::zz.

1630  {
1631  HGCSiliconDetId detid(id);
1632  double x(0), y(0);
1633  if (useWafer) {
1634  auto xyw = waferPositionNoRot(detid.layer(), detid.waferU(), detid.waferV(), reco, debug);
1635  x = xyw.first;
1636  y = xyw.second;
1637  }
1638  auto xy = getXY(detid.layer(), (x + loc.x()), (y + loc.y()), false);
1639  double zz = (detid.zside() < 0) ? -(loc.z() + waferZ(detid.layer(), reco)) : (loc.z() + waferZ(detid.layer(), reco));
1640  double xx = (detid.zside() < 0) ? -xy.first : xy.first;
1641  return GlobalPoint(xx, xy.second, zz);
1642 }
double waferZ(int layer, bool reco) const
std::pair< double, double > getXY(int layer, double x, double y, bool forwd) const
#define debug
Definition: HDRShower.cc:19
std::pair< double, double > waferPositionNoRot(int lay, int waferU, int waferV, bool reco, bool debug) const

◆ waferMax()

int HGCalDDDConstants::waferMax ( ) const
inline

Definition at line 178 of file HGCalDDDConstants.h.

References waferMax_.

Referenced by HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_().

178 { return waferMax_[1]; }
std::array< int, 4 > waferMax_

◆ waferMin()

int HGCalDDDConstants::waferMin ( ) const
inline

Definition at line 179 of file HGCalDDDConstants.h.

References waferMax_.

Referenced by HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_().

179 { return waferMax_[0]; }
std::array< int, 4 > waferMax_

◆ waferParameters()

std::pair< double, double > HGCalDDDConstants::waferParameters ( bool  reco) const

Definition at line 1564 of file HGCalDDDConstants.cc.

References hexside_, HGCalParameters::k_ScaleToDDD, and rmax_.

1564  {
1565  if (reco)
1566  return std::make_pair(rmax_, hexside_);
1567  else
1569 }
static constexpr double k_ScaleToDDD
fixed size matrix

◆ waferPosition() [1/3]

std::pair< double, double > HGCalDDDConstants::waferPosition ( int  wafer,
bool  reco 
) const

Definition at line 1571 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::k_ScaleToDDD, HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, geometryCSVtoXML::xx, and geometryCSVtoXML::yy.

Referenced by HGCalGeometryLoader::build(), HGCMouseBite::exclude(), and waferPositionNoRot().

1571  {
1572  double xx(0), yy(0);
1573  if (wafer >= 0 && wafer < static_cast<int>(hgpar_->waferPosX_.size())) {
1574  xx = hgpar_->waferPosX_[wafer];
1575  yy = hgpar_->waferPosY_[wafer];
1576  }
1577  if (!reco) {
1580  }
1581  return std::make_pair(xx, yy);
1582 }
std::vector< double > waferPosY_
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_
std::vector< double > waferPosX_

◆ waferPosition() [2/3]

std::pair< double, double > HGCalDDDConstants::waferPosition ( int  lay,
int  waferU,
int  waferV,
bool  reco,
bool  debug 
) const

Definition at line 1584 of file HGCalDDDConstants.cc.

References debug, HGCalParameters::firstLayer_, getXY(), hgpar_, HGCalParameters::layerType_, HGCalTypes::WaferCenterR, waferPositionNoRot(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and geometryCSVtoXML::xy.

1585  {
1586  int ll = lay - hgpar_->firstLayer_;
1587  bool rotx = ((!hgpar_->layerType_.empty()) && (hgpar_->layerType_[ll] == HGCalTypes::WaferCenterR));
1588 #ifdef EDM_ML_DEBUG
1589  if (debug)
1590  edm::LogVerbatim("HGCalGeom") << "Layer " << lay << ":" << ll << " Rotation " << rotx << " U:V " << waferU << ":"
1591  << waferV;
1592 #endif
1593  auto xy = waferPositionNoRot(lay, waferU, waferV, reco, debug);
1594  std::pair<double, double> xy0 = (rotx) ? getXY(lay, xy.first, xy.second, false) : xy;
1595 #ifdef EDM_ML_DEBUG
1596  if (debug)
1597  edm::LogVerbatim("HGCalGeom") << "Without and with rotation " << xy.first << ":" << xy.second << ":" << xy0.first
1598  << ":" << xy0.second;
1599 #endif
1600  return xy0;
1601 }
Log< level::Info, true > LogVerbatim
int32_t waferU(const int32_t index)
std::pair< double, double > getXY(int layer, double x, double y, bool forwd) const
std::vector< int > layerType_
static constexpr int32_t WaferCenterR
Definition: HGCalTypes.h:27
#define debug
Definition: HDRShower.cc:19
fixed size matrix
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_
std::pair< double, double > waferPositionNoRot(int lay, int waferU, int waferV, bool reco, bool debug) const

◆ waferPosition() [3/3]

std::pair< double, double > HGCalDDDConstants::waferPosition ( int  waferU,
int  waferV,
bool  reco 
) const
private

Definition at line 2098 of file HGCalDDDConstants.cc.

References hgpar_, HGCalParameters::k_ScaleToDDD, HGCalWaferIndex::waferIndex(), HGCalParameters::waferPosX_, HGCalParameters::waferPosY_, HGCalParameters::wafersInLayers_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), geometryCSVtoXML::xx, and geometryCSVtoXML::yy.

2098  {
2099  double xx(0), yy(0);
2100  int indx = HGCalWaferIndex::waferIndex(0, waferU, waferV);
2101  auto itr = hgpar_->wafersInLayers_.find(indx);
2102  if (itr != hgpar_->wafersInLayers_.end()) {
2103  xx = hgpar_->waferPosX_[itr->second];
2104  yy = hgpar_->waferPosY_[itr->second];
2105  }
2106  if (!reco) {
2109  }
2110  return std::make_pair(xx, yy);
2111 }
std::vector< double > waferPosY_
int32_t waferU(const int32_t index)
wafer_map wafersInLayers_
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
static constexpr double k_ScaleToDDD
fixed size matrix
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_
std::vector< double > waferPosX_

◆ waferPositionNoRot()

std::pair< double, double > HGCalDDDConstants::waferPositionNoRot ( int  lay,
int  waferU,
int  waferV,
bool  reco,
bool  debug 
) const
private

Definition at line 2074 of file HGCalDDDConstants.cc.

References debug, HGCalParameters::firstLayer_, hgpar_, HGCalParameters::k_ScaleToDDD, waferPosition(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, HGCalParameters::xLayerHex_, geometryCSVtoXML::xy, y, and HGCalParameters::yLayerHex_.

Referenced by isValidCell8(), localToGlobal8(), locateCell(), waferLocal2Global(), and waferPosition().

2075  {
2076  int ll = lay - hgpar_->firstLayer_;
2077  double x = hgpar_->xLayerHex_[ll];
2078  double y = hgpar_->yLayerHex_[ll];
2079 #ifdef EDM_ML_DEBUG
2080  if (debug)
2081  edm::LogVerbatim("HGCalGeom") << "Layer " << lay << ":" << ll << " Shift " << hgpar_->xLayerHex_[ll] << ":"
2082  << hgpar_->yLayerHex_[ll] << " U:V " << waferU << ":" << waferV;
2083 #endif
2084  if (!reco) {
2087  }
2088  const auto& xy = waferPosition(waferU, waferV, reco);
2089  x += xy.first;
2090  y += xy.second;
2091 #ifdef EDM_ML_DEBUG
2092  if (debug)
2093  edm::LogVerbatim("HGCalGeom") << "With wafer " << x << ":" << y << ":" << xy.first << ":" << xy.second;
2094 #endif
2095  return std::make_pair(x, y);
2096 }
Log< level::Info, true > LogVerbatim
int32_t waferU(const int32_t index)
std::vector< double > yLayerHex_
std::pair< double, double > waferPosition(int wafer, bool reco) const
static constexpr double k_ScaleToDDD
#define debug
Definition: HDRShower.cc:19
fixed size matrix
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_
std::vector< double > xLayerHex_

◆ wafers() [1/2]

int HGCalDDDConstants::wafers ( ) const

Definition at line 1644 of file HGCalDDDConstants.cc.

References HGCalParameters::depth_, hgpar_, mps_fire::i, layers(), HGCalParameters::moduleLayR_, modules(), and tileTrapezoid().

Referenced by HGCalTopology::allGeomModules(), and HGCalDDDConstants().

1644  {
1645  int wafer(0);
1646  if (!tileTrapezoid()) {
1647  for (unsigned int i = 0; i < layers(true); ++i) {
1648  int lay = hgpar_->depth_[i];
1649  wafer += modules(lay, true);
1650  }
1651  } else {
1652  wafer = static_cast<int>(hgpar_->moduleLayR_.size());
1653  }
1654  return wafer;
1655 }
std::vector< int > depth_
std::vector< int > moduleLayR_
unsigned int layers(bool reco) const
bool tileTrapezoid() const
int modules(int lay, bool reco) const
const HGCalParameters * hgpar_

◆ wafers() [2/2]

int HGCalDDDConstants::wafers ( int  layer,
int  type 
) const

Definition at line 1657 of file HGCalDDDConstants.cc.

References HGCalParameters::firstModule_, getIndex(), hgpar_, HGCalParameters::lastModule_, pixelTopology::layer, tileTrapezoid(), and waferLayer_.

1657  {
1658  int wafer(0);
1659  if (!tileTrapezoid()) {
1660  auto itr = waferLayer_.find(layer);
1661  if (itr != waferLayer_.end()) {
1662  unsigned ity = (type > 0 && type <= 2) ? type : 0;
1663  wafer = (itr->second)[ity];
1664  }
1665  } else {
1666  const auto& index = getIndex(layer, true);
1667  wafer = 1 + hgpar_->lastModule_[index.first] - hgpar_->firstModule_[index.first];
1668  }
1669  return wafer;
1670 }
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::map< int, HGCWaferParam > waferLayer_
std::vector< int > firstModule_
std::pair< int, float > getIndex(int lay, bool reco) const
bool tileTrapezoid() const
std::vector< int > lastModule_
const HGCalParameters * hgpar_

◆ waferSepar()

double HGCalDDDConstants::waferSepar ( bool  reco) const
inline

Definition at line 192 of file HGCalDDDConstants.h.

References hgpar_, HGCalParameters::k_ScaleToDDD, and HGCalParameters::sensorSeparation_.

192  {
194  }
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

◆ waferSize()

double HGCalDDDConstants::waferSize ( bool  reco) const
inline

Definition at line 195 of file HGCalDDDConstants.h.

References hgpar_, HGCalParameters::k_ScaleToDDD, and HGCalParameters::waferSize_.

Referenced by HFNoseSD::update(), and HGCalSD::update().

195  {
197  }
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

◆ waferToCopy()

int HGCalDDDConstants::waferToCopy ( int  wafer) const
inline

Definition at line 200 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::waferCopy_.

200  {
201  return ((wafer >= 0) && (wafer < static_cast<int>(hgpar_->waferCopy_.size())))
202  ? hgpar_->waferCopy_[wafer]
203  : static_cast<int>(hgpar_->waferCopy_.size());
204  }
std::vector< int > waferCopy_
const HGCalParameters * hgpar_

◆ waferType() [1/3]

int HGCalDDDConstants::waferType ( DetId const &  id,
bool  fromFile 
) const

Definition at line 1672 of file HGCalDDDConstants.cc.

References DetId::Forward, hgpar_, HFNoseDetId::layer(), HGCSiliconDetId::layer(), pixelTopology::layer, HGCSiliconDetId::type(), HFNoseDetId::type(), waferFileSize(), waferHexagon6(), waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, waferTypeL(), HGCalWaferIndex::waferU(), HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HGCalWaferIndex::waferV(), HGCSiliconDetId::waferV(), and HFNoseDetId::waferV().

Referenced by HGCalSimHitValidation::analyzeHits(), assignCellHex(), cellThickness(), HGCalNumberingScheme::getUnitID(), isHalfCell(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_().

1672  {
1673  int type(1);
1674  if (waferHexagon8()) {
1675  if (fromFile && (waferFileSize() > 0)) {
1676  int layer(0), waferU(0), waferV(0);
1677  if (id.det() != DetId::Forward) {
1678  HGCSiliconDetId hid(id);
1679  layer = hid.layer();
1680  waferU = hid.waferU();
1681  waferV = hid.waferV();
1682  } else {
1683  HFNoseDetId hid(id);
1684  layer = hid.layer();
1685  waferU = hid.waferU();
1686  waferV = hid.waferV();
1687  }
1689  if (itr != hgpar_->waferInfoMap_.end())
1690  type = (itr->second).type;
1691  } else {
1692  type = ((id.det() != DetId::Forward) ? HGCSiliconDetId(id).type() : HFNoseDetId(id).type());
1693  }
1694  } else if (waferHexagon6()) {
1695  type = waferTypeL(HGCalDetId(id).wafer()) - 1;
1696  }
1697  return type;
1698 }
int32_t waferU(const int32_t index)
bool waferHexagon8() const
bool waferHexagon6() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
unsigned int waferFileSize() const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
int type() const
get the type
Definition: HFNoseDetId.h:51
int32_t waferV(const int32_t index)
int waferTypeL(int wafer) const
int type() const
get the type
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_

◆ waferType() [2/3]

int HGCalDDDConstants::waferType ( int  layer,
int  waferU,
int  waferV,
bool  fromFile 
) const

Definition at line 1700 of file HGCalDDDConstants.cc.

References HGCSiliconDetId::HGCalCoarseThick, hgpar_, pixelTopology::layer, HGCalParameters::typesInLayers_, waferFileSize(), waferHexagon6(), waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalParameters::waferTypeL_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

1700  {
1702  if (waferHexagon8()) {
1703  if (fromFile && (waferFileSize() > 0)) {
1705  if (itr != hgpar_->waferInfoMap_.end())
1706  type = (itr->second).type;
1707  } else {
1709  if (itr != hgpar_->typesInLayers_.end())
1710  type = hgpar_->waferTypeL_[itr->second];
1711  }
1712  } else if (waferHexagon6()) {
1713  if ((waferU >= 0) && (waferU < static_cast<int>(hgpar_->waferTypeL_.size())))
1714  type = (hgpar_->waferTypeL_[waferU] - 1);
1715  }
1716  return type;
1717 }
int32_t waferU(const int32_t index)
bool waferHexagon8() const
bool waferHexagon6() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
unsigned int waferFileSize() const
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
wafer_map typesInLayers_
int32_t waferV(const int32_t index)
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_

◆ waferType() [3/3]

std::tuple< int, int, int > HGCalDDDConstants::waferType ( HGCSiliconDetId const &  id,
bool  fromFile 
) const

Definition at line 1719 of file HGCalDDDConstants.cc.

References dqmdumpme::first, hgpar_, HGCalGeomTools::k_allCorners, HGCalGeomTools::k_fiveCorners, HGCalGeomTools::k_fourCorners, HGCalTypes::k_OffsetRotation, HGCalGeomTools::k_threeCorners, pixelTopology::layer, hgcalPerformanceValidation::orient, edm::second(), HGCalParameters::typesInLayers_, waferFileSize(), HGCalTypes::WaferFive, HGCalTypes::WaferFull, HGCalTypes::WaferHalf, HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalTypes::WaferThree, HGCalParameters::waferTypeL_, HGCalParameters::waferTypes_, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

1719  {
1720  const auto& index = HGCalWaferIndex::waferIndex(id.layer(), id.waferU(), id.waferV());
1721  int type(-1), part(-1), orient(-1);
1722  if (fromFile && (waferFileSize() > 0)) {
1723  auto itr = hgpar_->waferInfoMap_.find(index);
1724  if (itr != hgpar_->waferInfoMap_.end()) {
1725  type = (itr->second).type;
1726  part = (itr->second).part;
1727  orient = (itr->second).orient;
1728  }
1729  } else {
1730  auto ktr = hgpar_->typesInLayers_.find(index);
1731  if (ktr != hgpar_->typesInLayers_.end())
1732  type = hgpar_->waferTypeL_[ktr->second];
1733  auto itr = hgpar_->waferTypes_.find(index);
1734  if (itr != hgpar_->waferTypes_.end()) {
1735  if ((itr->second).second < HGCalTypes::k_OffsetRotation) {
1736  orient = (itr->second).second;
1737  if ((itr->second).first == HGCalGeomTools::k_allCorners) {
1739  } else if ((itr->second).first == HGCalGeomTools::k_fiveCorners) {
1741  } else if ((itr->second).first == HGCalGeomTools::k_fourCorners) {
1743  } else if ((itr->second).first == HGCalGeomTools::k_threeCorners) {
1745  }
1746  } else {
1747  part = (itr->second).first;
1748  orient = ((itr->second).second - HGCalTypes::k_OffsetRotation);
1749  }
1750  } else {
1752  orient = 0;
1753  }
1754  }
1755  return std::make_tuple(type, part, orient);
1756 }
static constexpr int k_OffsetRotation
Definition: HGCalTypes.h:91
static constexpr int k_fourCorners
int32_t waferU(const int32_t index)
static constexpr int32_t WaferThree
Definition: HGCalTypes.h:42
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
U second(std::pair< T, U > const &p)
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
unsigned int waferFileSize() const
static constexpr int32_t WaferHalf
Definition: HGCalTypes.h:39
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
waferT_map waferTypes_
static constexpr int k_fiveCorners
part
Definition: HCALResponse.h:20
wafer_map typesInLayers_
static constexpr int k_threeCorners
int32_t waferV(const int32_t index)
static constexpr int32_t WaferFive
Definition: HGCalTypes.h:36
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
static constexpr int k_allCorners
std::vector< int > waferTypeL_

◆ waferTypeL()

int HGCalDDDConstants::waferTypeL ( int  wafer) const
inline

Definition at line 210 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::waferTypeL_.

Referenced by HGCalTBAnalyzer::analyzeRecHits(), PFRecHitQTestThresholdInThicknessNormalizedMIPs::test(), and waferType().

210  {
211  return ((wafer >= 0) && (wafer < static_cast<int>(hgpar_->waferTypeL_.size()))) ? hgpar_->waferTypeL_[wafer] : 0;
212  }
const HGCalParameters * hgpar_
std::vector< int > waferTypeL_

◆ waferTypeRotation()

std::pair< int, int > HGCalDDDConstants::waferTypeRotation ( int  layer,
int  waferU,
int  waferV,
bool  fromFile,
bool  debug 
) const

Definition at line 1758 of file HGCalDDDConstants.cc.

References debug, TauDecayModes::dec, dqmdumpme::first, hgpar_, HGCalGeomTools::k_allCorners, HGCalGeomTools::k_fiveCorners, HGCalGeomTools::k_fourCorners, HGCalTypes::k_OffsetRotation, HGCalGeomTools::k_threeCorners, pixelTopology::layer, hgcalPerformanceValidation::orient, edm::second(), HGCalTypes::WaferCorner0, waferFileSize(), HGCalTypes::WaferFive, HGCalTypes::WaferFull, HGCalTypes::WaferHalf, waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferInfoMap_, HGCalTypes::WaferOut, HGCalTypes::WaferThree, HGCalParameters::waferTypes_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and LaserClient_cfi::wl.

Referenced by isValidCell8(), and HGCalSD::setDetUnitId().

1759  {
1760  int type(HGCalTypes::WaferOut), rotn(0);
1762  bool withinList(true);
1763  if (fromFile && (waferFileSize() > 0)) {
1764  auto itr = hgpar_->waferInfoMap_.find(wl);
1765  withinList = (itr != hgpar_->waferInfoMap_.end());
1766  if (withinList) {
1767  type = (itr->second).part;
1768  rotn = (itr->second).orient;
1769  }
1770  } else {
1771  auto itr = hgpar_->waferTypes_.find(wl);
1772  if (waferHexagon8()) {
1773  withinList = (itr != hgpar_->waferTypes_.end());
1774  if (withinList) {
1775  if ((itr->second).second < HGCalTypes::k_OffsetRotation) {
1776  rotn = (itr->second).second;
1777  if ((itr->second).first == HGCalGeomTools::k_allCorners) {
1779  } else if ((itr->second).first == HGCalGeomTools::k_fiveCorners) {
1781  } else if ((itr->second).first == HGCalGeomTools::k_fourCorners) {
1783  } else if ((itr->second).first == HGCalGeomTools::k_threeCorners) {
1785  }
1786  } else {
1787  type = (itr->second).first;
1788  rotn = ((itr->second).second - HGCalTypes::k_OffsetRotation);
1789  }
1790  } else {
1792  rotn = HGCalTypes::WaferCorner0;
1793  }
1794  }
1795  }
1796 #ifdef EDM_ML_DEBUG
1797  if (debug)
1798  edm::LogVerbatim("HGCalGeom") << "waferTypeRotation: Layer " << layer << " Wafer " << waferU << ":" << waferV
1799  << " Index " << std::hex << wl << std::dec << ":" << withinList << " Type " << type
1800  << " Rotation " << rotn;
1801 #endif
1802  return std::make_pair(type, rotn);
1803 }
Log< level::Info, true > LogVerbatim
static constexpr int k_OffsetRotation
Definition: HGCalTypes.h:91
static constexpr int k_fourCorners
int32_t waferU(const int32_t index)
static constexpr int32_t WaferOut
Definition: HGCalTypes.h:56
static constexpr int32_t WaferThree
Definition: HGCalTypes.h:42
bool waferHexagon8() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
U second(std::pair< T, U > const &p)
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
unsigned int waferFileSize() const
static constexpr int32_t WaferHalf
Definition: HGCalTypes.h:39
static constexpr int32_t WaferCorner0
Definition: HGCalTypes.h:13
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
waferT_map waferTypes_
#define debug
Definition: HDRShower.cc:19
static constexpr int k_fiveCorners
part
Definition: HCALResponse.h:20
static constexpr int k_threeCorners
int32_t waferV(const int32_t index)
static constexpr int32_t WaferFive
Definition: HGCalTypes.h:36
waferInfo_map waferInfoMap_
const HGCalParameters * hgpar_
static constexpr int k_allCorners

◆ waferTypeT()

int HGCalDDDConstants::waferTypeT ( int  wafer) const
inline

Definition at line 206 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::waferTypeT_.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getClosestCell(), and hgcal::RecHitTools::isHalfCell().

206  {
207  return ((wafer >= 0) && (wafer < static_cast<int>(hgpar_->waferTypeT_.size()))) ? hgpar_->waferTypeT_[wafer] : 0;
208  }
std::vector< int > waferTypeT_
const HGCalParameters * hgpar_

◆ waferUVMax()

int HGCalDDDConstants::waferUVMax ( ) const
inline

Definition at line 217 of file HGCalDDDConstants.h.

References hgpar_, and HGCalParameters::waferUVMax_.

Referenced by HGCalTopology::HGCalTopology().

217 { return hgpar_->waferUVMax_; }
const HGCalParameters * hgpar_

◆ waferVirtual()

bool HGCalDDDConstants::waferVirtual ( int  layer,
int  waferU,
int  waferV 
) const

Definition at line 1805 of file HGCalDDDConstants.cc.

References hgpar_, pixelTopology::layer, waferHexagon6(), waferHexagon8(), HGCalWaferIndex::waferIndex(), HGCalParameters::waferTypes_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and LaserClient_cfi::wl.

1805  {
1806  bool type(false);
1807  if (waferHexagon8()) {
1809  type = (hgpar_->waferTypes_.find(wl) != hgpar_->waferTypes_.end());
1810  } else if (waferHexagon6()) {
1811  int wl = HGCalWaferIndex::waferIndex(layer, waferU, 0, true);
1812  type = (hgpar_->waferTypes_.find(wl) != hgpar_->waferTypes_.end());
1813  }
1814  return type;
1815 }
int32_t waferU(const int32_t index)
bool waferHexagon8() const
bool waferHexagon6() const
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
waferT_map waferTypes_
int32_t waferV(const int32_t index)
const HGCalParameters * hgpar_

◆ waferZ()

double HGCalDDDConstants::waferZ ( int  layer,
bool  reco 
) const

Definition at line 1817 of file HGCalDDDConstants.cc.

References getIndex(), hgpar_, HGCalParameters::k_ScaleToDDD, and HGCalParameters::zLayerHex_.

Referenced by HGCalSimHitValidation::analyzeHits(), HGCalTBAnalyzer::analyzeRecHits(), HGCalTBAnalyzer::analyzeSimHits(), ticl::SeedingRegionByTracks::buildFirstLayers(), ticl::LinkingAlgoByDirectionGeometric::buildLayers(), HFNoseNumberingScheme::checkPosition(), HGCalNumberingScheme::checkPosition(), HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), HGCalTriggerTools::getLayerZ(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), hgcal::RecHitTools::getPositionLayer(), ticl::LinkingAlgoByDirectionGeometric::linkTracksters(), HGCalGeometry::neighborZ(), HeterogeneousHGCalHEFCellPositionsFiller::set_conditions_(), and waferLocal2Global().

1817  {
1818  const auto& index = getIndex(lay, reco);
1819  if (index.first < 0)
1820  return 0;
1821  else
1823 }
std::pair< int, float > getIndex(int lay, bool reco) const
std::vector< double > zLayerHex_
static constexpr double k_ScaleToDDD
fixed size matrix
const HGCalParameters * hgpar_

Member Data Documentation

◆ dPhiMin

constexpr float HGCalDDDConstants::dPhiMin = 0.02
staticprivate

Definition at line 244 of file HGCalDDDConstants.h.

◆ fullAndPart_

const bool HGCalDDDConstants::fullAndPart_
private

Definition at line 252 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), isValidHex8(), and waferInLayerTest().

◆ geomTools_

HGCalGeomTools HGCalDDDConstants::geomTools_
private

Definition at line 242 of file HGCalDDDConstants.h.

◆ hexside_

double HGCalDDDConstants::hexside_
private

◆ hexsideT_

double HGCalDDDConstants::hexsideT_
private

Definition at line 254 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and waferFromPosition().

◆ hgcassette_

HGCalCassette HGCalDDDConstants::hgcassette_
private

◆ hgcell_

std::unique_ptr<HGCalCell> HGCalDDDConstants::hgcell_
private

Definition at line 240 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and locateCell().

◆ hgcellUV_

std::unique_ptr<HGCalCellUV> HGCalDDDConstants::hgcellUV_
private

Definition at line 241 of file HGCalDDDConstants.h.

Referenced by cellHex(), and HGCalDDDConstants().

◆ hgpar_

const HGCalParameters* HGCalDDDConstants::hgpar_
private

◆ k_horizontalShift

constexpr double HGCalDDDConstants::k_horizontalShift = 1.0
staticprivate

Definition at line 243 of file HGCalDDDConstants.h.

◆ max_modules_layer_

Simrecovecs HGCalDDDConstants::max_modules_layer_
private

Definition at line 257 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and modules().

◆ maxWafersPerLayer_

int32_t HGCalDDDConstants::maxWafersPerLayer_
private

Definition at line 258 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and maxModulesPerLayer().

◆ mode_

const HGCalGeometryMode::GeometryMode HGCalDDDConstants::mode_
private

◆ modHalf_

int32_t HGCalDDDConstants::modHalf_
private

Definition at line 255 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and maxModules().

◆ rmax_

double HGCalDDDConstants::rmax_
private

◆ rmaxT_

double HGCalDDDConstants::rmaxT_
private

Definition at line 254 of file HGCalDDDConstants.h.

Referenced by cellHex(), HGCalDDDConstants(), and waferFromPosition().

◆ sqrt3_

const double HGCalDDDConstants::sqrt3_
private

Definition at line 250 of file HGCalDDDConstants.h.

Referenced by cellHex(), and distFromEdgeHex().

◆ tan30deg_

constexpr double HGCalDDDConstants::tan30deg_ = 0.5773502693
staticprivate

◆ tol_

constexpr double HGCalDDDConstants::tol_ = 0.001
staticprivate

Definition at line 249 of file HGCalDDDConstants.h.

Referenced by assignCellTrap(), tilePhiEdge(), and tileRingEdge().

◆ tot_layers_

std::array<uint32_t, 2> HGCalDDDConstants::tot_layers_
private

Definition at line 256 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), lastLayer(), and layers().

◆ tot_wafers_

int32_t HGCalDDDConstants::tot_wafers_
private

Definition at line 255 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants().

◆ waferIn_

std::unordered_map<int32_t, bool> HGCalDDDConstants::waferIn_
private

Definition at line 261 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and isValidHex8().

◆ waferLayer_

std::map<int, HGCWaferParam> HGCalDDDConstants::waferLayer_
private

Definition at line 259 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), and wafers().

◆ waferMax_

std::array<int, 4> HGCalDDDConstants::waferMax_
private

Definition at line 260 of file HGCalDDDConstants.h.

Referenced by HGCalDDDConstants(), waferCount(), waferMax(), and waferMin().