File indexing completed on 2025-04-30 22:24:20
0001 #include "Geometry/HGCalCommonData/interface/HGCalParameters.h"
0002 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0003 #include "Geometry/HGCalCommonData/interface/HGCalWaferIndex.h"
0004
0005 #include <array>
0006
0007
0008
0009 HGCalParameters::HGCalParameters(const std::string& nam) : name_(nam), nCells_(0), waferMaskMode_(0) {
0010 #ifdef EDM_ML_DEBUG
0011 edm::LogVerbatim("HGCalGeom") << "Construct HGCalParameters for " << name_;
0012 #endif
0013 }
0014
0015 void HGCalParameters::fillModule(const HGCalParameters::hgtrap& mytr, bool reco) {
0016 if (reco) {
0017 moduleLayR_.emplace_back(mytr.lay);
0018 moduleBlR_.emplace_back(mytr.bl);
0019 moduleTlR_.emplace_back(mytr.tl);
0020 moduleHR_.emplace_back(mytr.h);
0021 moduleDzR_.emplace_back(mytr.dz);
0022 moduleAlphaR_.emplace_back(mytr.alpha);
0023 moduleCellR_.emplace_back(mytr.cellSize);
0024 } else {
0025 moduleLayS_.emplace_back(mytr.lay);
0026 moduleBlS_.emplace_back(mytr.bl);
0027 moduleTlS_.emplace_back(mytr.tl);
0028 moduleHS_.emplace_back(mytr.h);
0029 moduleDzS_.emplace_back(mytr.dz);
0030 moduleAlphaS_.emplace_back(mytr.alpha);
0031 moduleCellS_.emplace_back(mytr.cellSize);
0032 }
0033 }
0034
0035 HGCalParameters::hgtrap HGCalParameters::getModule(unsigned int k, bool reco) const {
0036 HGCalParameters::hgtrap mytr;
0037 if (reco) {
0038 if (k < moduleLayR_.size()) {
0039 mytr.lay = moduleLayR_[k];
0040 mytr.bl = moduleBlR_[k];
0041 mytr.tl = moduleTlR_[k];
0042 mytr.h = moduleHR_[k];
0043 mytr.dz = moduleDzR_[k];
0044 mytr.alpha = moduleAlphaR_[k];
0045 mytr.cellSize = moduleCellR_[k];
0046 } else {
0047 mytr.lay = -1;
0048 mytr.bl = mytr.tl = mytr.h = mytr.dz = mytr.alpha = mytr.cellSize = 0;
0049 }
0050 } else {
0051 if (k < moduleLayS_.size()) {
0052 mytr.lay = moduleLayS_[k];
0053 mytr.bl = moduleBlS_[k];
0054 mytr.tl = moduleTlS_[k];
0055 mytr.h = moduleHS_[k];
0056 mytr.dz = moduleDzS_[k];
0057 mytr.alpha = moduleAlphaS_[k];
0058 mytr.cellSize = moduleCellS_[k];
0059 } else {
0060 mytr.lay = -1;
0061 mytr.bl = mytr.tl = mytr.h = mytr.dz = mytr.alpha = mytr.cellSize = 0;
0062 }
0063 }
0064 return mytr;
0065 }
0066
0067 void HGCalParameters::fillTrForm(const HGCalParameters::hgtrform& mytr) {
0068 int zp = (mytr.zp == 1) ? 1 : 0;
0069 uint32_t indx = ((zp & kMaskZside) << kShiftZside);
0070 indx |= ((mytr.lay & kMaskLayer) << kShiftLayer);
0071 indx |= ((mytr.sec & kMaskSector) << kShiftSector);
0072 indx |= ((mytr.subsec & kMaskSubSec) << kShiftSubSec);
0073 #ifdef EDM_ML_DEBUG
0074 edm::LogVerbatim("HGCalGeom") << "ZP " << zp << ":" << kMaskZside << ":" << kShiftZside
0075 << ((zp & kMaskZside) << kShiftZside) << " Lay " << mytr.lay << ":" << kMaskLayer << ":"
0076 << kShiftLayer << ":" << ((mytr.lay & kMaskLayer) << kShiftLayer) << " Sector "
0077 << mytr.sec << ":" << kMaskSector << ":" << kShiftSector << ":"
0078 << ((mytr.sec & kMaskSector) << kShiftSector) << " SubSec " << mytr.subsec << ":"
0079 << kMaskSubSec << ":" << kShiftSubSec << ":"
0080 << ((mytr.subsec & kMaskSubSec) << kShiftSubSec) << " Index " << std::hex << indx
0081 << std::dec;
0082 #endif
0083 trformIndex_.emplace_back(indx);
0084 trformTranX_.emplace_back(mytr.h3v.x());
0085 trformTranY_.emplace_back(mytr.h3v.y());
0086 trformTranZ_.emplace_back(mytr.h3v.z());
0087 trformRotXX_.emplace_back((std::abs(mytr.hr.xx()) > tol) ? mytr.hr.xx() : 0);
0088 trformRotYX_.emplace_back((std::abs(mytr.hr.yx()) > tol) ? mytr.hr.yx() : 0);
0089 trformRotZX_.emplace_back((std::abs(mytr.hr.zx()) > tol) ? mytr.hr.zx() : 0);
0090 trformRotXY_.emplace_back((std::abs(mytr.hr.xy()) > tol) ? mytr.hr.xy() : 0);
0091 trformRotYY_.emplace_back((std::abs(mytr.hr.yy()) > tol) ? mytr.hr.yy() : 0);
0092 trformRotZY_.emplace_back((std::abs(mytr.hr.zy()) > tol) ? mytr.hr.zy() : 0);
0093 trformRotXZ_.emplace_back((std::abs(mytr.hr.xz()) > tol) ? mytr.hr.xz() : 0);
0094 trformRotYZ_.emplace_back((std::abs(mytr.hr.yz()) > tol) ? mytr.hr.yz() : 0);
0095 trformRotZZ_.emplace_back((std::abs(mytr.hr.zz()) > tol) ? mytr.hr.zz() : 0);
0096 #ifdef EDM_ML_DEBUG
0097 unsigned int k = trformIndex_.size() - 1;
0098 edm::LogVerbatim("HGCalGeom") << "HGCalParameters[" << k << "] Index " << std::hex << trformIndex_[k] << std::dec
0099 << " (" << mytr.zp << ", " << mytr.lay << ", " << mytr.sec << ", " << mytr.subsec
0100 << ") Translation (" << trformTranX_[k] << ", " << trformTranY_[k] << ", "
0101 << trformTranZ_[k] << ") Rotation (" << trformRotXX_[k] << ", " << trformRotYX_[k]
0102 << ", " << trformRotZX_[k] << ", " << trformRotXY_[k] << ", " << trformRotYY_[k] << ", "
0103 << trformRotZY_[k] << ", " << trformRotXZ_[k] << ", " << trformRotYZ_[k] << ", "
0104 << trformRotZZ_[k];
0105 #endif
0106 }
0107
0108 HGCalParameters::hgtrform HGCalParameters::getTrForm(unsigned int k) const {
0109 HGCalParameters::hgtrform mytr;
0110 if (k < trformIndex_.size()) {
0111 const auto& id = getID(k);
0112 mytr.zp = id[0];
0113 mytr.lay = id[1];
0114 mytr.sec = id[2];
0115 mytr.subsec = id[3];
0116 mytr.h3v = CLHEP::Hep3Vector(trformTranX_[k], trformTranY_[k], trformTranZ_[k]);
0117 const CLHEP::HepRep3x3 rotation(trformRotXX_[k],
0118 trformRotXY_[k],
0119 trformRotXZ_[k],
0120 trformRotYX_[k],
0121 trformRotYY_[k],
0122 trformRotYZ_[k],
0123 trformRotZX_[k],
0124 trformRotZY_[k],
0125 trformRotZZ_[k]);
0126 mytr.hr = CLHEP::HepRotation(rotation);
0127 } else {
0128 mytr.zp = mytr.lay = mytr.sec = mytr.subsec = 0;
0129 }
0130 #ifdef EDM_ML_DEBUG
0131 edm::LogVerbatim("HGCalGeom") << "HGCalParameters[" << k << "] Index " << std::hex << trformIndex_[k] << std::dec
0132 << " (" << mytr.zp << ", " << mytr.lay << ", " << mytr.sec << ", " << mytr.subsec
0133 << ") Translation (" << mytr.h3v.x() << ", " << mytr.h3v.y() << ", " << mytr.h3v.z()
0134 << ") Rotation (" << mytr.hr.xx() << ", " << mytr.hr.yx() << ", " << mytr.hr.zx()
0135 << ", " << mytr.hr.xy() << ", " << mytr.hr.yy() << ", " << mytr.hr.zy() << ", "
0136 << mytr.hr.xz() << ", " << mytr.hr.yz() << ", " << mytr.hr.zz();
0137 #endif
0138 return mytr;
0139 }
0140
0141 void HGCalParameters::addTrForm(const CLHEP::Hep3Vector& h3v) {
0142 unsigned int k = trformTranX_.size();
0143 if (k > 0) {
0144 trformTranX_[k - 1] += h3v.x();
0145 trformTranY_[k - 1] += h3v.y();
0146 trformTranZ_[k - 1] += h3v.z();
0147 }
0148 }
0149
0150 void HGCalParameters::scaleTrForm(double scale) {
0151 unsigned int k = trformTranX_.size();
0152 if (k > 0) {
0153 trformTranX_[k - 1] *= scale;
0154 trformTranY_[k - 1] *= scale;
0155 trformTranZ_[k - 1] *= scale;
0156 }
0157 }
0158
0159 std::array<int, 4> HGCalParameters::getID(unsigned int k) const {
0160 int zp = ((trformIndex_[k] >> kShiftZside) & kMaskZside);
0161 if (zp != 1)
0162 zp = -1;
0163 int lay = ((trformIndex_[k] >> kShiftLayer) & kMaskLayer);
0164 int sec = ((trformIndex_[k] >> kShiftSector) & kMaskSector);
0165 int subsec = ((trformIndex_[k] >> kShiftSubSec) & kMaskSubSec);
0166 return std::array<int, 4>{{zp, lay, sec, subsec}};
0167 }
0168
0169 double HGCalParameters::scintRing(int indx, int irad) const {
0170 double r(0);
0171 if (scintValidRing(indx, irad)) {
0172 r = 0.5 * (radiusLayer_[scintType(indx + 1)][irad - 1] + radiusLayer_[scintType(indx + 1)][irad]);
0173 }
0174 return r;
0175 }
0176
0177 #include "FWCore/Utilities/interface/typelookup.h"
0178
0179 TYPELOOKUP_DATA_REG(HGCalParameters);
