CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
HcalDDDSimConstants Class Reference

#include <HcalDDDSimConstants.h>

Public Member Functions

HcalCellType::HcalCell cell (const int &det, const int &zside, const int &depth, const int &etaR, const int &iphi) const
 
int findDepth (const int &det, const int &eta, const int &phi, const int &zside, const int &lay) const
 
unsigned int findLayer (const int &layer, const std::vector< HcalParameters::LayerItem > &layerGroup) const
 
std::vector< std::pair< double, double > > getConstHBHE (const int &type) const
 
int getDepthEta16 (const int &det, const int &phi, const int &zside) const
 
int getDepthEta16M (const int &det) const
 
int getDepthEta29 (const int &phi, const int &zside, const int &i) const
 
int getDepthEta29M (const int &i, const bool &planOne) const
 
std::pair< int, double > getDetEta (const double &eta, const int &depth) const
 
int getEta (const int &det, const int &lay, const double &hetaR) const
 
std::pair< int, int > getEtaDepth (const int &det, int etaR, const int &phi, const int &zside, int depth, const int &lay) const
 
double getEtaHO (const double &etaR, const double &x, const double &y, const double &z) const
 
const std::vector< double > & getEtaTableHF () const
 
int getFrontLayer (const int &det, const int &eta) const
 
const std::vector< double > & getGparHF () const
 
const std::vector< HcalDetId > & getIdHF2QIE () const
 
std::pair< int, int > getiEtaRange (const int &i) const
 
int getLastLayer (const int &det, const int &eta) const
 
double getLayer0Wt (const int &det, const int &phi, const int &zside) const
 
int getLayerBack (const int &det, const int &eta, const int &phi, const int &zside, const int &depth) const
 
int getLayerFront (const int &det, const int &eta, const int &phi, const int &zside, const int &depth) const
 
int getLayerMax (const int &eta, const int &depth) const
 
int getMaxDepth (const int &det, const int &eta, const int &phi, const int &zside, const bool &partialOnly) const
 
int getMaxDepth (const int &type) const
 
std::pair< int, int > getMaxDepthDet (const int &i) const
 
int getMinDepth (const int &det, const int &eta, const int &phi, const int &zside, const bool &partialOnly) const
 
std::pair< int, int > getModHalfHBHE (const int &type) const
 
std::pair< double, double > getPhiCons (const int &det, const int &ieta) const
 
std::vector< std::pair< int, double > > getPhis (const int &subdet, const int &ieta) const
 
const std::vector< double > & getPhiTableHF () const
 
const std::vector< double > & getRTableHF () const
 
std::vector< HcalCellType > HcalCellTypes () const
 
std::vector< HcalCellType > HcalCellTypes (const HcalSubdetector &, int ieta=-1, int depth=-1) const
 
 HcalDDDSimConstants (const HcalParameters *hp)
 
bool isBH () const
 
const HcalLayerDepthMap * ldMap () const
 
int maxHFDepth (const int &ieta, const int &iphi) const
 
unsigned int numberOfCells (const HcalSubdetector &) const
 
const HcalParameters * parameter () const
 
int phiNumber (const int &phi, const int &unit) const
 
void printTiles () const
 
int unitPhi (const double &dphi) const
 
int unitPhi (const int &det, const int &etaR) const
 
 ~HcalDDDSimConstants ()
 

Private Member Functions

double deltaEta (const int &det, const int &eta, const int &depth) const
 
double getEta (const double &r, const double &z) const
 
double getEta (const int &det, const int &etaR, const int &zside, int depth=1) const
 
double getGain (const HcalSubdetector &subdet, const int &depth) const
 
int getShift (const HcalSubdetector &subdet, const int &depth) const
 
void initialize ()
 
unsigned int layerGroup (int det, int eta, int phi, int zside, int i) const
 
unsigned int layerGroup (int eta, int i) const
 
unsigned int layerGroupSize (int eta) const
 
void printTileHB (const int &eta, const int &phi, const int &zside, const int &depth) const
 
void printTileHE (const int &eta, const int &phi, const int &zside, const int &depth) const
 

Private Attributes

int depthEta16 [2]
 
int depthEta29 [2]
 
std::pair< int, int > depthMaxDf_
 
std::pair< int, int > depthMaxSp_
 
std::vector< int > depths [nDepthMax]
 
double dlShort
 
double dzVcal
 
double etaHO [4]
 
const HcalParameters * hpar
 
std::vector< HcalDetId > idHF2QIE
 
bool isBH_
 
int layBHB [3]
 
int layBHE [4]
 
int layFHB [2]
 
int layFHE [3]
 
HcalLayerDepthMap ldmap_
 
std::vector< int > maxDepth
 
int nDepth
 
int nEta
 
int nmodHB
 
int nmodHE
 
int nPhiF
 
int nR
 
int nzHB
 
int nzHE
 
double rminHO
 
double zVcal
 

Static Private Attributes

static const int maxLayer_ = 18
 
static const int maxLayerHB_ = 16
 
static const int nDepthMax = 9
 

Detailed Description

this class reads the constant section of the hcal-sim-numbering xml-file

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 24 of file HcalDDDSimConstants.h.

Constructor & Destructor Documentation

◆ HcalDDDSimConstants()

HcalDDDSimConstants::HcalDDDSimConstants ( const HcalParameters *  hp)

Definition at line 10 of file HcalDDDSimConstants.cc.

10  : hpar(hp) {
11 #ifdef EDM_ML_DEBUG
12  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::HcalDDDSimConstants (const HcalParameters* hp) constructor\n";
13 #endif
14 
15  initialize();
16 #ifdef EDM_ML_DEBUG
17  std::vector<HcalCellType> cellTypes = HcalCellTypes();
18  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal (All)\n";
19 #endif
20 }

References HcalCellTypes(), and initialize().

◆ ~HcalDDDSimConstants()

HcalDDDSimConstants::~HcalDDDSimConstants ( )

Definition at line 22 of file HcalDDDSimConstants.cc.

22  {
23 #ifdef EDM_ML_DEBUG
24  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::destructed!!!\n";
25 #endif
26 }

Member Function Documentation

◆ cell()

HcalCellType::HcalCell HcalDDDSimConstants::cell ( const int &  det,
const int &  zside,
const int &  depth,
const int &  etaR,
const int &  iphi 
) const

Definition at line 28 of file HcalDDDSimConstants.cc.

29  {
30  double etaMn = hpar->etaMin[0];
31  double etaMx = hpar->etaMax[0];
32  if (idet == static_cast<int>(HcalEndcap)) {
33  etaMn = hpar->etaMin[1];
34  etaMx = hpar->etaMax[1];
35  } else if (idet == static_cast<int>(HcalForward)) {
36  etaMn = hpar->etaMin[2];
37  etaMx = hpar->etaMax[2];
38  }
39  double eta = 0, deta = 0, phi = 0, dphi = 0, rz = 0, drz = 0;
40  bool ok = false, flagrz = true;
41  if ((idet == static_cast<int>(HcalBarrel) || idet == static_cast<int>(HcalEndcap) ||
42  idet == static_cast<int>(HcalOuter) || idet == static_cast<int>(HcalForward)) &&
43  etaR >= etaMn && etaR <= etaMx && depth > 0)
44  ok = true;
45  if (idet == static_cast<int>(HcalEndcap) && depth > (int)(hpar->zHE.size()))
46  ok = false;
47  else if (idet == static_cast<int>(HcalBarrel) && depth > maxLayerHB_ + 1)
48  ok = false;
49  else if (idet == static_cast<int>(HcalOuter) && depth != 4)
50  ok = false;
51  else if (idet == static_cast<int>(HcalForward) && depth > maxDepth[2])
52  ok = false;
53  if (ok) {
54  eta = getEta(idet, etaR, zside, depth);
55  deta = deltaEta(idet, etaR, depth);
56  double fibin, fioff;
57  if (idet == static_cast<int>(HcalBarrel) || idet == static_cast<int>(HcalOuter)) {
58  fioff = hpar->phioff[0];
59  fibin = hpar->phibin[etaR - 1];
60  } else if (idet == static_cast<int>(HcalEndcap)) {
61  fioff = hpar->phioff[1];
62  fibin = hpar->phibin[etaR - 1];
63  } else {
64  fioff = hpar->phioff[2];
65  fibin = hpar->phitable[etaR - hpar->etaMin[2]];
66  if (unitPhi(fibin) > 2)
67  fioff = hpar->phioff[4];
68  }
69  phi = -fioff + (iphi - 0.5) * fibin;
70  dphi = 0.5 * fibin;
71  if (idet == static_cast<int>(HcalForward)) {
72  int ir = nR + hpar->etaMin[2] - etaR - 1;
73  if (ir > 0 && ir < nR) {
74  rz = 0.5 * (hpar->rTable[ir] + hpar->rTable[ir - 1]);
75  drz = 0.5 * (hpar->rTable[ir] - hpar->rTable[ir - 1]);
76  } else {
77  ok = false;
78 #ifdef EDM_ML_DEBUG
79  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: wrong eta " << etaR << " (" << ir << "/" << nR
80  << ") Detector " << idet << std::endl;
81 #endif
82  }
83  } else if (etaR <= nEta) {
84  int laymin(depth), laymax(depth);
85  if (idet == static_cast<int>(HcalOuter)) {
86  laymin = (etaR > hpar->noff[2]) ? ((int)(hpar->zHE.size())) : ((int)(hpar->zHE.size())) - 1;
87  laymax = ((int)(hpar->zHE.size()));
88  }
89  double d1 = 0, d2 = 0;
90  if (idet == static_cast<int>(HcalEndcap)) {
91  flagrz = false;
92  d1 = hpar->zHE[laymin - 1] - hpar->dzHE[laymin - 1];
93  d2 = hpar->zHE[laymax - 1] + hpar->dzHE[laymax - 1];
94  } else {
95  d1 = hpar->rHB[laymin - 1] - hpar->drHB[laymin - 1];
96  d2 = hpar->rHB[laymax - 1] + hpar->drHB[laymax - 1];
97  }
98  rz = 0.5 * (d2 + d1);
99  drz = 0.5 * (d2 - d1);
100  } else {
101  ok = false;
102  edm::LogWarning("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth << " or etaR " << etaR
103  << " for detector " << idet;
104  }
105  } else {
106  ok = false;
107  edm::LogWarning("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth << " det " << idet;
108  }
109  HcalCellType::HcalCell tmp(ok, eta, deta, phi, dphi, rz, drz, flagrz);
110 
111 #ifdef EDM_ML_DEBUG
112  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants: det/side/depth/etaR/"
113  << "phi " << idet << "/" << zside << "/" << depth << "/" << etaR << "/" << iphi
114  << " Cell Flag " << tmp.ok << " " << tmp.eta << " " << tmp.deta << " phi " << tmp.phi
115  << " " << tmp.dphi << " r(z) " << tmp.rz << " " << tmp.drz << " " << tmp.flagrz;
116 #endif
117  return tmp;
118 }

References d1, deltaEta(), LEDCalibrationChannels::depth, HcalParameters::drHB, HcalParameters::dzHE, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, getEta(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, createfilelist::int, LEDCalibrationChannels::iphi, maxDepth, maxLayerHB_, nEta, HcalParameters::noff, nR, convertSQLiteXML::ok, phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHB, HcalParameters::rTable, createJobs::tmp, unitPhi(), HcalParameters::zHE, and ecaldqm::zside().

Referenced by HcalCellTypes(), and HcalTestAnalysis::update().

◆ deltaEta()

double HcalDDDSimConstants::deltaEta ( const int &  det,
const int &  eta,
const int &  depth 
) const
private

Definition at line 979 of file HcalDDDSimConstants.cc.

979  {
980  double tmp = 0;
981  if (det == static_cast<int>(HcalForward)) {
982  int ir = nR + hpar->etaMin[2] - etaR - 1;
983  if (ir > 0 && ir < nR) {
984  double z = zVcal;
985  if (depth % 2 != 1)
986  z += dlShort;
987  tmp = 0.5 * (getEta(hpar->rTable[ir - 1], z) - getEta(hpar->rTable[ir], z));
988  }
989  } else {
990  if (etaR > 0 && etaR < nEta) {
991  if (etaR == hpar->noff[1] - 1 && depth > depthEta29[0]) {
992  tmp = 0.5 * (hpar->etaTable[etaR + 1] - hpar->etaTable[etaR - 1]);
993  } else if (det == static_cast<int>(HcalOuter)) {
994  if (etaR == hpar->noff[2]) {
995  tmp = 0.5 * (etaHO[0] - hpar->etaTable[etaR - 1]);
996  } else if (etaR == hpar->noff[2] + 1) {
997  tmp = 0.5 * (hpar->etaTable[etaR] - etaHO[1]);
998  } else if (etaR == hpar->noff[3]) {
999  tmp = 0.5 * (etaHO[2] - hpar->etaTable[etaR - 1]);
1000  } else if (etaR == hpar->noff[3] + 1) {
1001  tmp = 0.5 * (hpar->etaTable[etaR] - etaHO[3]);
1002  } else {
1003  tmp = 0.5 * (hpar->etaTable[etaR] - hpar->etaTable[etaR - 1]);
1004  }
1005  } else {
1006  tmp = 0.5 * (hpar->etaTable[etaR] - hpar->etaTable[etaR - 1]);
1007  }
1008  }
1009  }
1010 #ifdef EDM_ML_DEBUG
1011  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants::deltaEta " << etaR << " " << depth << " ==> " << tmp;
1012 #endif
1013  return tmp;
1014 }

References LEDCalibrationChannels::depth, depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, createJobs::tmp, z, and zVcal.

Referenced by cell().

◆ findDepth()

int HcalDDDSimConstants::findDepth ( const int &  det,
const int &  eta,
const int &  phi,
const int &  zside,
const int &  lay 
) const

Definition at line 120 of file HcalDDDSimConstants.cc.

121  {
122  int depth = (ldmap_.isValid(det, phi, zside)) ? ldmap_.getDepth(det, eta, phi, zside, lay) : -1;
123  return depth;
124 }

References LEDCalibrationChannels::depth, PVValHelper::eta, HcalLayerDepthMap::getDepth(), HcalLayerDepthMap::isValid(), ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getHCID(), and layerGroup().

◆ findLayer()

unsigned int HcalDDDSimConstants::findLayer ( const int &  layer,
const std::vector< HcalParameters::LayerItem > &  layerGroup 
) const

Definition at line 126 of file HcalDDDSimConstants.cc.

127  {
128  unsigned int id = layerGroup.size();
129  for (unsigned int i = 0; i < layerGroup.size(); i++) {
130  if (layer == (int)(layerGroup[i].layer)) {
131  id = i;
132  break;
133  }
134  }
135  return id;
136 }

References mps_fire::i, triggerObjects_cff::id, phase1PixelTopology::layer, and layerGroup().

◆ getConstHBHE()

std::vector< std::pair< double, double > > HcalDDDSimConstants::getConstHBHE ( const int &  type) const

Definition at line 138 of file HcalDDDSimConstants.cc.

138  {
139  std::vector<std::pair<double, double> > gcons;
140  if (type == 0) {
141  for (unsigned int i = 0; i < hpar->rHB.size(); ++i) {
142  gcons.emplace_back(std::pair<double, double>(hpar->rHB[i], hpar->drHB[i]));
143  }
144  } else {
145  for (unsigned int i = 0; i < hpar->zHE.size(); ++i) {
146  gcons.emplace_back(std::pair<double, double>(hpar->zHE[i], hpar->dzHE[i]));
147  }
148  }
149  return gcons;
150 }

References HcalParameters::drHB, HcalParameters::dzHE, hpar, mps_fire::i, HcalParameters::rHB, and HcalParameters::zHE.

◆ getDepthEta16()

int HcalDDDSimConstants::getDepthEta16 ( const int &  det,
const int &  phi,
const int &  zside 
) const

Definition at line 152 of file HcalDDDSimConstants.cc.

152  {
153  int depth = ldmap_.getDepth16(det, phi, zside);
154  if (depth < 0)
155  depth = (det == 2) ? depthEta16[1] : depthEta16[0];
156 #ifdef EDM_ML_DEBUG
157  edm::LogVerbatim("HcalGeom") << "getDepthEta16: " << det << ":" << depth;
158 #endif
159  return depth;
160 }

References LEDCalibrationChannels::depth, depthEta16, HcalLayerDepthMap::getDepth16(), ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getDepthEta16(), getEtaDepth(), HcalDDDRecConstants::getHCID(), getMaxDepth(), getMinDepth(), HcalDDDRecConstants::getOneEtaBin(), HcalDDDRecConstants::initialize(), and printTiles().

◆ getDepthEta16M()

int HcalDDDSimConstants::getDepthEta16M ( const int &  det) const

Definition at line 162 of file HcalDDDSimConstants.cc.

162  {
163  int depth = (det == 2) ? depthEta16[1] : depthEta16[0];
164  std::vector<int> phis;
165  int detsp = ldmap_.validDet(phis);
166  if (detsp == det) {
167  int zside = (phis[0] > 0) ? 1 : -1;
168  int iphi = (phis[0] > 0) ? phis[0] : -phis[0];
169  int depthsp = ldmap_.getDepth16(det, iphi, zside);
170  if (det == 1 && depthsp > depth)
171  depth = depthsp;
172  if (det == 2 && depthsp < depth)
173  depth = depthsp;
174  }
175  return depth;
176 }

References LEDCalibrationChannels::depth, depthEta16, HcalLayerDepthMap::getDepth16(), LEDCalibrationChannels::iphi, ldmap_, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getMinDepth(), HcalDDDRecConstants::getOneEtaBin(), and HcalDDDRecConstants::initialize().

◆ getDepthEta29()

int HcalDDDSimConstants::getDepthEta29 ( const int &  phi,
const int &  zside,
const int &  i 
) const

Definition at line 178 of file HcalDDDSimConstants.cc.

178  {
179  int depth = (i == 0) ? ldmap_.getMaxDepthLastHE(2, phi, zside) : -1;
180  if (depth < 0)
181  depth = (i == 1) ? depthEta29[1] : depthEta29[0];
182  return depth;
183 }

References LEDCalibrationChannels::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), mps_fire::i, ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getDepthEta29(), getEtaDepth(), HcalDDDRecConstants::getHCID(), HcalDDDRecConstants::getOneEtaBin(), and HcalDDDRecConstants::mergedDepthList29().

◆ getDepthEta29M()

int HcalDDDSimConstants::getDepthEta29M ( const int &  i,
const bool &  planOne 
) const

Definition at line 185 of file HcalDDDSimConstants.cc.

185  {
186  int depth = (i == 1) ? depthEta29[1] : depthEta29[0];
187  if (i == 0 && planOne) {
188  std::vector<int> phis;
189  int detsp = ldmap_.validDet(phis);
190  if (detsp == 2) {
191  int zside = (phis[0] > 0) ? 1 : -1;
192  int iphi = (phis[0] > 0) ? phis[0] : -phis[0];
193  int depthsp = ldmap_.getMaxDepthLastHE(2, iphi, zside);
194  if (depthsp > depth)
195  depth = depthsp;
196  }
197  }
198  return depth;
199 }

References LEDCalibrationChannels::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), mps_fire::i, LEDCalibrationChannels::iphi, ldmap_, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getOneEtaBin(), and HcalDDDRecConstants::initialize().

◆ getDetEta()

std::pair< int, double > HcalDDDSimConstants::getDetEta ( const double &  eta,
const int &  depth 
) const

Definition at line 201 of file HcalDDDSimConstants.cc.

201  {
202  int hsubdet(0), ieta(0);
203  double etaR(0);
204  double heta = fabs(eta);
205  for (int i = 0; i < nEta; i++)
206  if (heta > hpar->etaTable[i])
207  ieta = i + 1;
208  if (heta <= hpar->etaRange[1]) {
209  if (((ieta == hpar->etaMin[1] && depth == depthEta16[1]) || (ieta > hpar->etaMax[0])) &&
210  (ieta <= hpar->etaMax[1])) {
211  hsubdet = static_cast<int>(HcalEndcap);
212  } else {
213  hsubdet = static_cast<int>(HcalBarrel);
214  }
215  etaR = eta;
216  } else {
217  hsubdet = static_cast<int>(HcalForward);
218  double theta = 2. * atan(exp(-heta));
219  double hR = zVcal * tan(theta);
220  etaR = (eta >= 0. ? hR : -hR);
221  }
222  return std::pair<int, double>(hsubdet, etaR);
223 }

References LEDCalibrationChannels::depth, depthEta16, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, generateTowerEtThresholdLUT::etaRange, HcalParameters::etaTable, JetChargeProducer_cfi::exp, HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, LEDCalibrationChannels::ieta, nEta, funct::tan(), theta(), and zVcal.

Referenced by HcalNumberingFromDDD::unitID().

◆ getEta() [1/3]

double HcalDDDSimConstants::getEta ( const double &  r,
const double &  z 
) const
private

Definition at line 1056 of file HcalDDDSimConstants.cc.

1056  {
1057  double tmp = 0;
1058  if (z != 0)
1059  tmp = -log(tan(0.5 * atan(r / z)));
1060 #ifdef EDM_ML_DEBUG
1061  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants::getEta " << r << " " << z << " ==> " << tmp;
1062 #endif
1063  return tmp;
1064 }

References dqm-mbProfile::log, alignCSCRings::r, funct::tan(), createJobs::tmp, and z.

◆ getEta() [2/3]

double HcalDDDSimConstants::getEta ( const int &  det,
const int &  etaR,
const int &  zside,
int  depth = 1 
) const
private

Definition at line 1016 of file HcalDDDSimConstants.cc.

1016  {
1017  double tmp = 0;
1018  if (det == static_cast<int>(HcalForward)) {
1019  int ir = nR + hpar->etaMin[2] - etaR - 1;
1020  if (ir > 0 && ir < nR) {
1021  double z = zVcal;
1022  if (depth % 2 != 1)
1023  z += dlShort;
1024  tmp = 0.5 * (getEta(hpar->rTable[ir - 1], z) + getEta(hpar->rTable[ir], z));
1025  }
1026  } else {
1027  if (etaR > 0 && etaR < nEta) {
1028  if (etaR == hpar->noff[1] - 1 && depth > depthEta29[0]) {
1029  tmp = 0.5 * (hpar->etaTable[etaR + 1] + hpar->etaTable[etaR - 1]);
1030  } else if (det == static_cast<int>(HcalOuter)) {
1031  if (etaR == hpar->noff[2]) {
1032  tmp = 0.5 * (etaHO[0] + hpar->etaTable[etaR - 1]);
1033  } else if (etaR == hpar->noff[2] + 1) {
1034  tmp = 0.5 * (hpar->etaTable[etaR] + etaHO[1]);
1035  } else if (etaR == hpar->noff[3]) {
1036  tmp = 0.5 * (etaHO[2] + hpar->etaTable[etaR - 1]);
1037  } else if (etaR == hpar->noff[3] + 1) {
1038  tmp = 0.5 * (hpar->etaTable[etaR] + etaHO[3]);
1039  } else {
1040  tmp = 0.5 * (hpar->etaTable[etaR] + hpar->etaTable[etaR - 1]);
1041  }
1042  } else {
1043  tmp = 0.5 * (hpar->etaTable[etaR] + hpar->etaTable[etaR - 1]);
1044  }
1045  }
1046  }
1047  if (zside == 0)
1048  tmp = -tmp;
1049 #ifdef EDM_ML_DEBUG
1050  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants::getEta " << etaR << " " << zside << " " << depth << " ==> "
1051  << tmp;
1052 #endif
1053  return tmp;
1054 }

References LEDCalibrationChannels::depth, depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, createJobs::tmp, z, ecaldqm::zside(), and zVcal.

◆ getEta() [3/3]

int HcalDDDSimConstants::getEta ( const int &  det,
const int &  lay,
const double &  hetaR 
) const

Definition at line 225 of file HcalDDDSimConstants.cc.

225  {
226  int ieta(0);
227  if (det == static_cast<int>(HcalForward)) { // Forward HCal
228  ieta = hpar->etaMax[2];
229  for (int i = nR - 1; i > 0; i--)
230  if (hetaR < hpar->rTable[i])
231  ieta = hpar->etaMin[2] + nR - i - 1;
232  } else { // Barrel or Endcap
233  ieta = 1;
234  for (int i = 0; i < nEta - 1; i++)
235  if (hetaR > hpar->etaTable[i])
236  ieta = i + 1;
237  if (det == static_cast<int>(HcalBarrel)) {
238  if (ieta > hpar->etaMax[0])
239  ieta = hpar->etaMax[0];
240  if (lay == maxLayer_) {
241  if (hetaR > etaHO[1] && ieta == hpar->noff[2])
242  ieta++;
243  }
244  } else if (det == static_cast<int>(HcalEndcap)) {
245  if (ieta <= hpar->etaMin[1])
246  ieta = hpar->etaMin[1];
247  }
248  }
249  return ieta;
250 }

References etaHO, HcalParameters::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, HcalParameters::etaMin, HcalParameters::etaTable, HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, LEDCalibrationChannels::ieta, maxLayer_, nEta, HcalParameters::noff, and nR.

Referenced by cell(), deltaEta(), getEta(), and HcalNumberingFromDDD::unitID().

◆ getEtaDepth()

std::pair< int, int > HcalDDDSimConstants::getEtaDepth ( const int &  det,
int  etaR,
const int &  phi,
const int &  zside,
int  depth,
const int &  lay 
) const

Definition at line 252 of file HcalDDDSimConstants.cc.

253  {
254 #ifdef EDM_ML_DEBUG
255  edm::LogVerbatim("HcalGeom") << "HcalDDDEsimConstants:getEtaDepth: I/P " << det << ":" << etaR << ":" << phi << ":"
256  << zside << ":" << depth << ":" << lay;
257 #endif
258  //Modify the depth index
259  if ((det == static_cast<int>(HcalEndcap)) && (etaR == 17) && (lay == 1))
260  etaR = 18;
261  if (det == static_cast<int>(HcalForward)) { // Forward HCal
262  } else if (det == static_cast<int>(HcalOuter)) {
263  depth = 4;
264  } else {
265  if (lay >= 0) {
266  depth = layerGroup(det, etaR, phi, zside, lay - 1);
267  if (etaR == hpar->noff[0] && lay > 1) {
268  int kphi = phi + int((hpar->phioff[3] + 0.1) / hpar->phibin[etaR - 1]);
269  kphi = (kphi - 1) % 4 + 1;
270  if (kphi == 2 || kphi == 3)
271  depth = layerGroup(det, etaR, phi, zside, lay - 2);
272  }
273  } else if (det == static_cast<int>(HcalBarrel)) {
274  if (depth > getMaxDepth(det, etaR, phi, zside, false))
275  depth = getMaxDepth(det, etaR, phi, zside, false);
276  }
277  if (etaR >= hpar->noff[1] && depth > getDepthEta29(phi, zside, 0)) {
278  etaR -= getDepthEta29(phi, zside, 1);
279  } else if (etaR == hpar->etaMin[1]) {
280  if (det == static_cast<int>(HcalBarrel)) {
281  if (depth > getDepthEta16(det, phi, zside))
282  depth = getDepthEta16(det, phi, zside);
283  } else {
284  if (depth < getDepthEta16(det, phi, zside))
285  depth = getDepthEta16(det, phi, zside);
286  }
287  }
288  }
289 #ifdef EDM_ML_DEBUG
290  edm::LogVerbatim("HcalGeom") << "HcalDDDEsimConstants:getEtaDepth: O/P " << etaR << ":" << depth;
291 #endif
292  return std::pair<int, int>(etaR, depth);
293 }

References LEDCalibrationChannels::depth, HcalParameters::etaMin, getDepthEta16(), getDepthEta29(), getMaxDepth(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, createfilelist::int, layerGroup(), HcalParameters::noff, phi, HcalParameters::phibin, HcalParameters::phioff, and ecaldqm::zside().

Referenced by HcalTestNumberingTester::analyze(), and HcalNumberingFromDDD::unitID().

◆ getEtaHO()

double HcalDDDSimConstants::getEtaHO ( const double &  etaR,
const double &  x,
const double &  y,
const double &  z 
) const

Definition at line 295 of file HcalDDDSimConstants.cc.

295  {
296  if (hpar->zHO.size() > 4) {
297  double eta = fabs(etaR);
298  double r = std::sqrt(x * x + y * y);
299  if (r > rminHO) {
300  double zz = fabs(z);
301  if (zz > hpar->zHO[3]) {
302  if (eta <= hpar->etaTable[10])
303  eta = hpar->etaTable[10] + 0.001;
304  } else if (zz > hpar->zHO[1]) {
305  if (eta <= hpar->etaTable[4])
306  eta = hpar->etaTable[4] + 0.001;
307  }
308  }
309  eta = (z >= 0. ? eta : -eta);
310 #ifdef EDM_ML_DEBUG
311  edm::LogVerbatim("HcalGeom") << "R " << r << " Z " << z << " eta " << etaR << ":" << eta;
312  if (eta != etaR)
313  edm::LogInfo("HCalGeom") << "**** Check *****";
314 #endif
315  return eta;
316  } else {
317  return etaR;
318  }
319 }

References PVValHelper::eta, HcalParameters::etaTable, hpar, alignCSCRings::r, rminHO, mathSSE::sqrt(), x, y, z, HcalParameters::zHO, and geometryCSVtoXML::zz.

Referenced by HcalNumberingFromDDD::unitID().

◆ getEtaTableHF()

const std::vector<double>& HcalDDDSimConstants::getEtaTableHF ( ) const
inline

Definition at line 44 of file HcalDDDSimConstants.h.

44 { return hpar->etaTableHF; }

References HcalParameters::etaTableHF, and hpar.

◆ getFrontLayer()

int HcalDDDSimConstants::getFrontLayer ( const int &  det,
const int &  eta 
) const

Definition at line 321 of file HcalDDDSimConstants.cc.

321  {
322  int lay = 0;
323  if (det == 1) {
324  if (std::abs(eta) == 16)
325  lay = layFHB[1];
326  else
327  lay = layFHB[0];
328  } else {
329  if (std::abs(eta) == 16)
330  lay = layFHE[1];
331  else if (std::abs(eta) == 18)
332  lay = layFHE[2];
333  else
334  lay = layFHE[0];
335  }
336  return lay;
337 }

References funct::abs(), PVValHelper::eta, layFHB, and layFHE.

Referenced by HcalDDDRecConstants::getLayerFront().

◆ getGain()

double HcalDDDSimConstants::getGain ( const HcalSubdetector &  subdet,
const int &  depth 
) const
private

Definition at line 1085 of file HcalDDDSimConstants.cc.

1085  {
1086  double gain;
1087  switch (subdet) {
1088  case HcalEndcap:
1089  gain = hpar->HEGains[0];
1090  break;
1091  case HcalForward:
1092  gain = hpar->HFGains[(depth - 1) % 2];
1093  break;
1094  case HcalOuter:
1095  gain = hpar->HBGains[3];
1096  break;
1097  default:
1098  gain = hpar->HBGains[0];
1099  break;
1100  }
1101  return gain;
1102 }

References LEDCalibrationChannels::depth, PedestalClient_cfi::gain, HcalParameters::HBGains, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEGains, HcalParameters::HFGains, and hpar.

Referenced by HcalCellTypes().

◆ getGparHF()

const std::vector<double>& HcalDDDSimConstants::getGparHF ( ) const
inline

Definition at line 45 of file HcalDDDSimConstants.h.

45 { return hpar->gparHF; }

References HcalParameters::gparHF, and hpar.

Referenced by HCalSD::HCalSD(), HFFibre::HFFibre(), HFShower::HFShower(), HFShowerLibrary::HFShowerLibrary(), and HFShowerParam::HFShowerParam().

◆ getIdHF2QIE()

const std::vector<HcalDetId>& HcalDDDSimConstants::getIdHF2QIE ( ) const
inline

Definition at line 46 of file HcalDDDSimConstants.h.

46 { return idHF2QIE; }

References idHF2QIE.

Referenced by HcalDDDRecConstants::getHFCellParameters().

◆ getiEtaRange()

std::pair<int, int> HcalDDDSimConstants::getiEtaRange ( const int &  i) const
inline

Definition at line 43 of file HcalDDDSimConstants.h.

43 { return std::pair<int, int>(hpar->etaMin[i], hpar->etaMax[i]); }

References HcalParameters::etaMax, HcalParameters::etaMin, hpar, and mps_fire::i.

◆ getLastLayer()

int HcalDDDSimConstants::getLastLayer ( const int &  det,
const int &  eta 
) const

Definition at line 339 of file HcalDDDSimConstants.cc.

339  {
340  int lay = 0;
341  if (det == 1) {
342  if (std::abs(eta) == 15)
343  lay = layBHB[1];
344  else if (std::abs(eta) == 16)
345  lay = layBHB[2];
346  else
347  lay = layBHB[0];
348  } else {
349  if (std::abs(eta) == 16)
350  lay = layBHE[1];
351  else if (std::abs(eta) == 17)
352  lay = layBHE[2];
353  else if (std::abs(eta) == 18)
354  lay = layBHE[3];
355  else
356  lay = layBHE[0];
357  }
358  return lay;
359 }

References funct::abs(), PVValHelper::eta, layBHB, and layBHE.

Referenced by HcalDDDRecConstants::getLayerBack(), and HcalDDDRecConstants::getThickActive().

◆ getLayer0Wt()

double HcalDDDSimConstants::getLayer0Wt ( const int &  det,
const int &  phi,
const int &  zside 
) const

Definition at line 361 of file HcalDDDSimConstants.cc.

361  {
362  double wt = ldmap_.getLayer0Wt(det, phi, zside);
363  if (wt < 0)
364  wt = (det == 2) ? hpar->Layer0Wt[1] : hpar->Layer0Wt[0];
365  return wt;
366 }

References HcalLayerDepthMap::getLayer0Wt(), hpar, HcalParameters::Layer0Wt, ldmap_, phi, and ecaldqm::zside().

Referenced by HCalSD::getEnergyDeposit(), and HcalDDDRecConstants::getLayer0Wt().

◆ getLayerBack()

int HcalDDDSimConstants::getLayerBack ( const int &  det,
const int &  eta,
const int &  phi,
const int &  zside,
const int &  depth 
) const

Definition at line 387 of file HcalDDDSimConstants.cc.

388  {
389  int layer = ldmap_.getLayerBack(det, eta, phi, zside, depth);
390  if (layer < 0) {
391  if (det == 1 || det == 2) {
392  layer = depths[depth - 1][eta - 1];
393  } else {
394  layer = maxLayer_;
395  }
396  }
397  return layer;
398 }

References LEDCalibrationChannels::depth, depths, PVValHelper::eta, HcalLayerDepthMap::getLayerBack(), phase1PixelTopology::layer, ldmap_, maxLayer_, phi, and ecaldqm::zside().

Referenced by printTileHB(), and printTileHE().

◆ getLayerFront()

int HcalDDDSimConstants::getLayerFront ( const int &  det,
const int &  eta,
const int &  phi,
const int &  zside,
const int &  depth 
) const

Definition at line 368 of file HcalDDDSimConstants.cc.

369  {
370  int layer = ldmap_.getLayerFront(det, eta, phi, zside, depth);
371  if (layer < 0) {
372  if (det == 1 || det == 2) {
373  layer = 1;
374  for (int l = 0; l < getLayerMax(eta, depth); ++l) {
375  if ((int)(layerGroup(eta - 1, l)) == depth) {
376  layer = l + 1;
377  break;
378  }
379  }
380  } else {
381  layer = (eta > hpar->noff[2]) ? maxLayerHB_ + 1 : maxLayer_;
382  }
383  }
384  return layer;
385 }

References LEDCalibrationChannels::depth, PVValHelper::eta, HcalLayerDepthMap::getLayerFront(), getLayerMax(), hpar, cmsLHEtoEOSManager::l, phase1PixelTopology::layer, layerGroup(), ldmap_, maxLayer_, maxLayerHB_, HcalParameters::noff, phi, and ecaldqm::zside().

Referenced by printTileHB(), and printTileHE().

◆ getLayerMax()

int HcalDDDSimConstants::getLayerMax ( const int &  eta,
const int &  depth 
) const

Definition at line 400 of file HcalDDDSimConstants.cc.

400  {
401  int layermx = ((eta < hpar->etaMin[1]) && depth - 1 < maxDepth[0]) ? maxLayerHB_ + 1 : (int)layerGroupSize(eta - 1);
402  return layermx;
403 }

References LEDCalibrationChannels::depth, PVValHelper::eta, HcalParameters::etaMin, hpar, createfilelist::int, layerGroupSize(), maxDepth, and maxLayerHB_.

Referenced by getLayerFront(), and initialize().

◆ getMaxDepth() [1/2]

int HcalDDDSimConstants::getMaxDepth ( const int &  det,
const int &  eta,
const int &  phi,
const int &  zside,
const bool &  partialOnly 
) const

Definition at line 405 of file HcalDDDSimConstants.cc.

406  {
407  int dmax(-1);
408  if (partialDetOnly) {
409  if (ldmap_.isValid(det, phi, zside)) {
410  dmax = ldmap_.getDepths(eta).second;
411  }
412  } else if (det == 1 || det == 2) {
413  if (ldmap_.isValid(det, phi, zside))
414  dmax = ldmap_.getDepths(eta).second;
415  else if (det == 2)
416  dmax = (maxDepth[1] > 0) ? layerGroup(eta - 1, maxLayer_) : 0;
417  else if (eta == hpar->etaMax[0])
418  dmax = getDepthEta16(det, phi, zside);
419  else
420  dmax = layerGroup(eta - 1, maxLayerHB_);
421  } else if (det == 3) { // HF
422  dmax = maxHFDepth(zside * eta, phi);
423  } else if (det == 4) { // HO
424  dmax = maxDepth[3];
425  } else {
426  dmax = -1;
427  }
428  return dmax;
429 }

References PVValHelper::eta, HcalParameters::etaMax, getDepthEta16(), HcalLayerDepthMap::getDepths(), hpar, HcalLayerDepthMap::isValid(), layerGroup(), ldmap_, maxDepth, maxHFDepth(), maxLayer_, maxLayerHB_, phi, and ecaldqm::zside().

◆ getMaxDepth() [2/2]

int HcalDDDSimConstants::getMaxDepth ( const int &  type) const
inline

Definition at line 53 of file HcalDDDSimConstants.h.

53 { return maxDepth[type]; }

References maxDepth.

Referenced by getEtaDepth(), HcalDDDRecConstants::getMaxDepth(), HCalSD::HCalSD(), and printTiles().

◆ getMaxDepthDet()

std::pair<int, int> HcalDDDSimConstants::getMaxDepthDet ( const int &  i) const
inline

Definition at line 55 of file HcalDDDSimConstants.h.

55 { return ((i == 1) ? depthMaxDf_ : depthMaxSp_); }

References depthMaxDf_, depthMaxSp_, and mps_fire::i.

Referenced by HcalDDDRecConstants::initialize().

◆ getMinDepth()

int HcalDDDSimConstants::getMinDepth ( const int &  det,
const int &  eta,
const int &  phi,
const int &  zside,
const bool &  partialOnly 
) const

Definition at line 431 of file HcalDDDSimConstants.cc.

432  {
433  int lmin(-1);
434  if (partialDetOnly) {
435  if (ldmap_.isValid(det, phi, zside)) {
436  lmin = ldmap_.getDepths(eta).first;
437  }
438  } else if (det == 3) { // HF
439  lmin = 1;
440  } else if (det == 4) { // HO
441  lmin = maxDepth[3];
442  } else {
443  if (ldmap_.isValid(det, phi, zside)) {
444  lmin = ldmap_.getDepths(eta).first;
445  } else if (layerGroupSize(eta - 1) > 0) {
446  lmin = (int)(layerGroup(eta - 1, 0));
447  unsigned int type = (det == 1) ? 0 : 1;
448  if (type == 1 && eta == hpar->etaMin[1])
449  lmin = getDepthEta16(det, phi, zside);
450  } else {
451  lmin = 1;
452  }
453  }
454  return lmin;
455 }

References PVValHelper::eta, HcalParameters::etaMin, getDepthEta16(), HcalLayerDepthMap::getDepths(), hpar, createfilelist::int, HcalLayerDepthMap::isValid(), layerGroup(), layerGroupSize(), ldmap_, maxDepth, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getMinDepth().

◆ getModHalfHBHE()

std::pair< int, int > HcalDDDSimConstants::getModHalfHBHE ( const int &  type) const

Definition at line 457 of file HcalDDDSimConstants.cc.

457  {
458  if (type == 0) {
459  return std::pair<int, int>(nmodHB, nzHB);
460  } else {
461  return std::pair<int, int>(nmodHE, nzHE);
462  }
463 }

References nmodHB, nmodHE, nzHB, and nzHE.

◆ getPhiCons()

std::pair< double, double > HcalDDDSimConstants::getPhiCons ( const int &  det,
const int &  ieta 
) const

Definition at line 465 of file HcalDDDSimConstants.cc.

465  {
466  double fioff(0), fibin(0);
467  if (det == static_cast<int>(HcalForward)) { // Forward HCal
468  fioff = hpar->phioff[2];
469  fibin = hpar->phitable[ieta - hpar->etaMin[2]];
470  if (unitPhi(fibin) > 2) { // HF double-phi
471  fioff = hpar->phioff[4];
472  }
473  } else { // Barrel or Endcap
474  if (det == static_cast<int>(HcalEndcap)) {
475  fioff = hpar->phioff[1];
476  } else {
477  fioff = hpar->phioff[0];
478  }
479  fibin = hpar->phibin[ieta - 1];
480  }
481  return std::pair<double, double>(fioff, fibin);
482 }

References HcalParameters::etaMin, HcalEndcap, HcalForward, hpar, LEDCalibrationChannels::ieta, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, and unitPhi().

Referenced by getPhis(), HcalDDDRecConstants::getPhis(), and HcalNumberingFromDDD::unitID().

◆ getPhis()

std::vector< std::pair< int, double > > HcalDDDSimConstants::getPhis ( const int &  subdet,
const int &  ieta 
) const

Definition at line 484 of file HcalDDDSimConstants.cc.

484  {
485  std::vector<std::pair<int, double> > phis;
486  int ietaAbs = (ieta > 0) ? ieta : -ieta;
487  std::pair<double, double> ficons = getPhiCons(subdet, ietaAbs);
488  int nphi = int((2._pi + 0.1 * ficons.second) / ficons.second);
489  int units = unitPhi(subdet, ietaAbs);
490  for (int ifi = 0; ifi < nphi; ++ifi) {
491  double phi = -ficons.first + (ifi + 0.5) * ficons.second;
492  int iphi = phiNumber(ifi + 1, units);
493  phis.emplace_back(std::pair<int, double>(iphi, phi));
494  }
495 #ifdef EDM_ML_DEBUG
496  edm::LogVerbatim("HcalGeom") << "getPhis: subdet|ieta|iphi " << subdet << "|" << ieta << " with " << phis.size()
497  << " phi bins";
498  for (unsigned int k = 0; k < phis.size(); ++k)
499  edm::LogVerbatim("HcalGeom") << "[" << k << "] iphi " << phis[k].first << " phi "
500  << convertRadToDeg(phis[k].second);
501 #endif
502  return phis;
503 }

References angle_units::operators::convertRadToDeg(), getPhiCons(), LEDCalibrationChannels::ieta, createfilelist::int, LEDCalibrationChannels::iphi, dqmdumpme::k, nphi, phi, phiNumber(), edm::second(), unitPhi(), and units().

Referenced by HcalCellTypes().

◆ getPhiTableHF()

const std::vector<double>& HcalDDDSimConstants::getPhiTableHF ( ) const
inline

Definition at line 60 of file HcalDDDSimConstants.h.

60 { return hpar->phitable; }

References hpar, and HcalParameters::phitable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), and HFShowerLibrary::HFShowerLibrary().

◆ getRTableHF()

const std::vector<double>& HcalDDDSimConstants::getRTableHF ( ) const
inline

Definition at line 61 of file HcalDDDSimConstants.h.

61 { return hpar->rTable; }

References hpar, and HcalParameters::rTable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), HFFibre::HFFibre(), HFShowerFibreBundle::HFShowerFibreBundle(), HFShowerLibrary::HFShowerLibrary(), and HFShowerPMT::HFShowerPMT().

◆ getShift()

int HcalDDDSimConstants::getShift ( const HcalSubdetector &  subdet,
const int &  depth 
) const
private

Definition at line 1066 of file HcalDDDSimConstants.cc.

1066  {
1067  int shift;
1068  switch (subdet) {
1069  case HcalEndcap:
1070  shift = hpar->HEShift[0];
1071  break;
1072  case HcalForward:
1073  shift = hpar->HFShift[(depth - 1) % 2];
1074  break;
1075  case HcalOuter:
1076  shift = hpar->HBShift[3];
1077  break;
1078  default:
1079  shift = hpar->HBShift[0];
1080  break;
1081  }
1082  return shift;
1083 }

References LEDCalibrationChannels::depth, HcalParameters::HBShift, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEShift, HcalParameters::HFShift, hpar, and edm::shift.

Referenced by HcalCellTypes().

◆ HcalCellTypes() [1/2]

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes ( ) const

Definition at line 505 of file HcalDDDSimConstants.cc.

505  {
506  std::vector<HcalCellType> cellTypes = HcalCellTypes(HcalBarrel);
507 #ifdef EDM_ML_DEBUG
508  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal Barrel\n";
509  for (unsigned int i = 0; i < cellTypes.size(); i++)
510  edm::LogInfo("HCalGeom") << "Cell " << i << " " << cellTypes[i] << "\n";
511 #endif
512 
513  std::vector<HcalCellType> hoCells = HcalCellTypes(HcalOuter);
514 #ifdef EDM_ML_DEBUG
515  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << hoCells.size() << " cells of type HCal Outer\n";
516  for (unsigned int i = 0; i < hoCells.size(); i++)
517  edm::LogInfo("HCalGeom") << "Cell " << i << " " << hoCells[i] << "\n";
518 #endif
519  cellTypes.insert(cellTypes.end(), hoCells.begin(), hoCells.end());
520 
521  std::vector<HcalCellType> heCells = HcalCellTypes(HcalEndcap);
522 #ifdef EDM_ML_DEBUG
523  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << heCells.size() << " cells of type HCal Endcap\n";
524  for (unsigned int i = 0; i < heCells.size(); i++)
525  edm::LogInfo("HCalGeom") << "Cell " << i << " " << heCells[i] << "\n";
526 #endif
527  cellTypes.insert(cellTypes.end(), heCells.begin(), heCells.end());
528 
529  std::vector<HcalCellType> hfCells = HcalCellTypes(HcalForward);
530 #ifdef EDM_ML_DEBUG
531  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << hfCells.size() << " cells of type HCal Forward\n";
532  for (unsigned int i = 0; i < hfCells.size(); i++)
533  edm::LogInfo("HCalGeom") << "Cell " << i << " " << hfCells[i] << "\n";
534 #endif
535  cellTypes.insert(cellTypes.end(), hfCells.begin(), hfCells.end());
536 
537  return cellTypes;
538 }

References HcalBarrel, HcalEndcap, HcalForward, HcalOuter, and mps_fire::i.

Referenced by HcalDDDRecConstants::HcalCellTypes(), HcalDDDSimConstants(), and numberOfCells().

◆ HcalCellTypes() [2/2]

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes ( const HcalSubdetector &  subdet,
int  ieta = -1,
int  depth = -1 
) const

Definition at line 540 of file HcalDDDSimConstants.cc.

542  {
543  std::vector<HcalCellType> cellTypes;
544  if (subdet == HcalForward) {
545  if (dzVcal < 0)
546  return cellTypes;
547  }
548 
549  int dmin, dmax, indx, nz;
550  double hsize = 0;
551  switch (subdet) {
552  case HcalEndcap:
553  dmin = 1;
554  dmax = (maxDepth[1] > 0) ? maxLayer_ + 1 : 0;
555  indx = 1;
556  nz = nzHE;
557  break;
558  case HcalForward:
559  dmin = 1;
560  dmax = (!idHF2QIE.empty()) ? 2 : maxDepth[2];
561  indx = 2;
562  nz = 2;
563  break;
564  case HcalOuter:
565  dmin = 4;
566  dmax = 4;
567  indx = 0;
568  nz = nzHB;
569  break;
570  default:
571  dmin = 1;
572  dmax = maxLayerHB_ + 1;
573  indx = 0;
574  nz = nzHB;
575  break;
576  }
577  if (depthl > 0)
578  dmin = dmax = depthl;
579  int ietamin = (ieta > 0) ? ieta : hpar->etaMin[indx];
580  int ietamax = (ieta > 0) ? ieta : hpar->etaMax[indx];
581  int phi = (indx == 2) ? 3 : 1;
582 
583  // Get the Cells
584  int subdet0 = static_cast<int>(subdet);
585  for (int depth = dmin; depth <= dmax; depth++) {
586  int shift = getShift(subdet, depth);
587  double gain = getGain(subdet, depth);
588  if (subdet == HcalForward) {
589  if (depth % 2 == 1)
590  hsize = dzVcal;
591  else
592  hsize = dzVcal - 0.5 * dlShort;
593  }
594  for (int eta = ietamin; eta <= ietamax; eta++) {
595  for (int iz = 0; iz < nz; ++iz) {
596  int zside = -2 * iz + 1;
597  HcalCellType::HcalCell temp1 = cell(subdet0, zside, depth, eta, phi);
598  if (temp1.ok) {
599  std::vector<std::pair<int, double> > phis = getPhis(subdet0, eta);
600  HcalCellType temp2(subdet, eta, zside, depth, temp1, shift, gain, hsize);
601  double dphi, fioff;
602  std::vector<int> phiMiss2;
603  if ((subdet == HcalBarrel) || (subdet == HcalOuter)) {
604  fioff = hpar->phioff[0];
605  dphi = hpar->phibin[eta - 1];
606  if (subdet == HcalOuter) {
607  if (eta == hpar->noff[4]) {
608  int kk = (iz == 0) ? 7 : (7 + hpar->noff[5]);
609  for (int miss = 0; miss < hpar->noff[5 + iz]; miss++) {
610  phiMiss2.emplace_back(hpar->noff[kk]);
611  kk++;
612  }
613  }
614  }
615  } else if (subdet == HcalEndcap) {
616  fioff = hpar->phioff[1];
617  dphi = hpar->phibin[eta - 1];
618  } else {
619  fioff = hpar->phioff[2];
620  dphi = hpar->phitable[eta - hpar->etaMin[2]];
621  if (unitPhi(dphi) > 2)
622  fioff = hpar->phioff[4];
623  }
624  int unit = unitPhi(dphi);
625  temp2.setPhi(phis, phiMiss2, fioff, dphi, unit);
626  cellTypes.emplace_back(temp2);
627  // For HF look at extra cells
628  if ((subdet == HcalForward) && (!idHF2QIE.empty())) {
629  HcalCellType temp3(subdet, eta, zside + 2, depth, temp1, shift, gain, hsize);
630  std::vector<int> phiMiss3;
631  for (auto& phi : phis) {
632  bool ok(false);
633  for (auto l : idHF2QIE) {
634  if ((eta * zside == l.ieta()) && (phi.first == l.iphi())) {
635  ok = true;
636  break;
637  }
638  }
639  if (!ok)
640  phiMiss3.emplace_back(phi.first);
641  }
642  dphi = hpar->phitable[eta - hpar->etaMin[2]];
643  unit = unitPhi(dphi);
644  fioff = (unit > 2) ? hpar->phioff[4] : hpar->phioff[2];
645  temp3.setPhi(phis, phiMiss2, fioff, dphi, unit);
646  cellTypes.emplace_back(temp3);
647  }
648  }
649  }
650  }
651  }
652  return cellTypes;
653 }

References cell(), LEDCalibrationChannels::depth, dlShort, dzVcal, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, PedestalClient_cfi::gain, getGain(), getPhis(), getShift(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, LEDCalibrationChannels::ieta, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, maxDepth, maxLayer_, maxLayerHB_, HcalParameters::noff, nzHB, nzHE, HcalCellType::HcalCell::ok, convertSQLiteXML::ok, phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalCellType::setPhi(), edm::shift, unit(), unitPhi(), and ecaldqm::zside().

◆ initialize()

void HcalDDDSimConstants::initialize ( void  )
private

Definition at line 743 of file HcalDDDSimConstants.cc.

743  {
744  nEta = hpar->etaTable.size();
745  nR = hpar->rTable.size();
746  nPhiF = nR - 1;
747  isBH_ = false;
748 
749 #ifdef EDM_ML_DEBUG
750  for (int i = 0; i < nEta - 1; ++i) {
751  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants:Read LayerGroup" << i << ":";
752  for (unsigned int k = 0; k < layerGroupSize(i); k++)
753  edm::LogVerbatim("HcalGeom") << " [" << k << "] = " << layerGroup(i, k);
754  }
755 #endif
756 
757  // Geometry parameters for HF
758  dlShort = hpar->gparHF[0];
759  zVcal = hpar->gparHF[4];
760  dzVcal = hpar->dzVcal;
761 #ifdef EDM_ML_DEBUG
762  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants: dlShort " << dlShort << " zVcal " << zVcal << " and dzVcal "
763  << dzVcal;
764 #endif
765 
766  //Transform some of the parameters
767  maxDepth = hpar->maxDepth;
768  maxDepth[0] = maxDepth[1] = 0;
769  for (int i = 0; i < nEta - 1; ++i) {
770  unsigned int imx = layerGroupSize(i);
771  int laymax = (imx > 0) ? layerGroup(i, imx - 1) : 0;
772  if (i < hpar->etaMax[0]) {
773  int laymax0 = (imx > 16) ? layerGroup(i, 16) : laymax;
774  if (i + 1 == hpar->etaMax[0] && laymax0 > 2)
775  laymax0 = 2;
776  if (maxDepth[0] < laymax0)
777  maxDepth[0] = laymax0;
778  }
779  if (i >= hpar->etaMin[1] - 1 && i < hpar->etaMax[1]) {
780  if (maxDepth[1] < laymax)
781  maxDepth[1] = laymax;
782  }
783  }
784 #ifdef EDM_ML_DEBUG
785  for (int i = 0; i < 4; ++i)
786  edm::LogVerbatim("HcalGeom") << "Detector Type [" << i << "] iEta " << hpar->etaMin[i] << ":" << hpar->etaMax[i]
787  << " MaxDepth " << maxDepth[i];
788 #endif
789 
790  int maxdepth = (maxDepth[1] > maxDepth[0]) ? maxDepth[1] : maxDepth[0];
791  for (int i = 0; i < maxdepth; ++i) {
792  for (int k = 0; k < nEta - 1; ++k) {
793  int layermx = getLayerMax(k + 1, i + 1);
794  int ll = layermx;
795  for (int l = layermx - 1; l >= 0; --l) {
796  if ((int)layerGroup(k, l) == i + 1) {
797  ll = l + 1;
798  break;
799  }
800  }
801  depths[i].emplace_back(ll);
802  }
803 
804 #ifdef EDM_ML_DEBUG
805  edm::LogVerbatim("HcalGeom") << "Depth " << i + 1 << " with " << depths[i].size() << " etas:";
806  for (int k = 0; k < nEta - 1; ++k)
807  edm::LogVerbatim("HcalGeom") << " [" << k << "] " << depths[i][k];
808 #endif
809  }
810 
811  nzHB = hpar->modHB[1];
812  nmodHB = hpar->modHB[0];
813  nzHE = hpar->modHE[1];
814  nmodHE = hpar->modHE[0];
815 #ifdef EDM_ML_DEBUG
816  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants:: " << nzHB << ":" << nmodHB << " barrel and " << nzHE << ":"
817  << nmodHE << " endcap half-sectors";
818 #endif
819 
820  if (hpar->rHB.size() > maxLayerHB_ + 1 && hpar->zHO.size() > 4) {
821  rminHO = hpar->rHO[0];
822  for (int k = 0; k < 4; ++k)
823  etaHO[k] = hpar->rHO[k + 1];
824  } else {
825  rminHO = -1.0;
826  etaHO[0] = hpar->etaTable[4];
827  etaHO[1] = hpar->etaTable[4];
828  etaHO[2] = hpar->etaTable[10];
829  etaHO[3] = hpar->etaTable[10];
830  }
831 #ifdef EDM_ML_DEBUG
832  edm::LogVerbatim("HcalGeom") << "HO Eta boundaries " << etaHO[0] << " " << etaHO[1] << " " << etaHO[2] << " "
833  << etaHO[3];
834  edm::LogVerbatim("HcalGeom") << "HO Parameters " << rminHO << " " << hpar->zHO.size();
835  for (int i = 0; i < 4; ++i)
836  edm::LogVerbatim("HcalGeom") << " eta[" << i << "] = " << etaHO[i];
837  for (unsigned int i = 0; i < hpar->zHO.size(); ++i)
838  edm::LogVerbatim("HcalGeom") << " zHO[" << i << "] = " << hpar->zHO[i];
839 #endif
840 
841  int noffsize = 7 + hpar->noff[5] + hpar->noff[6];
842  int noffl(noffsize + 5);
843  if ((int)(hpar->noff.size()) > (noffsize + 3)) {
844  depthEta16[0] = hpar->noff[noffsize];
845  depthEta16[1] = hpar->noff[noffsize + 1];
846  depthEta29[0] = hpar->noff[noffsize + 2];
847  depthEta29[1] = hpar->noff[noffsize + 3];
848  if ((int)(hpar->noff.size()) > (noffsize + 4)) {
849  noffl += (2 * hpar->noff[noffsize + 4]);
850  if ((int)(hpar->noff.size()) > noffl)
851  isBH_ = (hpar->noff[noffl] > 0);
852  }
853  } else {
854  depthEta16[0] = 2;
855  depthEta16[1] = 3;
856  depthEta29[0] = 2;
857  depthEta29[1] = 1;
858  }
859 #ifdef EDM_ML_DEBUG
860  edm::LogVerbatim("HcalGeom") << "isBH_ " << hpar->noff.size() << ":" << noffsize << ":" << noffl << ":" << isBH_;
861  edm::LogVerbatim("HcalGeom") << "Depth index at ieta = 16 for HB (max) " << depthEta16[0] << " HE (min) "
862  << depthEta16[1] << "; max depth for itea = 29 : (" << depthEta29[0] << ":"
863  << depthEta29[1] << ")";
864 #endif
865 
866  if ((int)(hpar->noff.size()) > (noffsize + 4)) {
867  int npair = hpar->noff[noffsize + 4];
868  int kk = noffsize + 4;
869  for (int k = 0; k < npair; ++k) {
870  idHF2QIE.emplace_back(HcalDetId(HcalForward, hpar->noff[kk + 1], hpar->noff[kk + 2], 1));
871  kk += 2;
872  }
873  }
874 #ifdef EDM_ML_DEBUG
875  edm::LogVerbatim("HcalGeom") << idHF2QIE.size() << " detector channels having 2 QIE cards:";
876  for (unsigned int k = 0; k < idHF2QIE.size(); ++k)
877  edm::LogVerbatim("HcalGeom") << " [" << k << "] " << idHF2QIE[k];
878 #endif
879 
880  layFHB[0] = 0;
881  layFHB[1] = 1;
882  layBHB[0] = 16;
883  layBHB[1] = 15;
884  layBHB[2] = 8;
885  if (maxDepth[1] == 0) {
886  layFHE[0] = layFHE[1] = layFHE[2] = 0;
887  layBHE[0] = layBHE[1] = layBHE[2] = layBHE[3] = 0;
888  } else {
889  layFHE[0] = 1;
890  layFHE[1] = 4;
891  layFHE[2] = 0;
892  layBHE[0] = 18;
893  layBHE[1] = 9;
894  layBHE[2] = 14;
895  layBHE[3] = 16;
896  }
897  depthMaxSp_ = std::pair<int, int>(0, 0);
898  int noffk(noffsize + 5);
899  if ((int)(hpar->noff.size()) > (noffsize + 5)) {
900  noffk += (2 * hpar->noff[noffsize + 4]);
901  if ((int)(hpar->noff.size()) >= noffk + 7) {
902  int dtype = hpar->noff[noffk + 1];
903  int nphi = hpar->noff[noffk + 2];
904  int ndeps = hpar->noff[noffk + 3];
905  int ndp16 = hpar->noff[noffk + 4];
906  int ndp29 = hpar->noff[noffk + 5];
907  double wt = 0.1 * (hpar->noff[noffk + 6]);
908  if ((int)(hpar->noff.size()) >= (noffk + 7 + nphi + 3 * ndeps)) {
909  if (dtype == 1 || dtype == 2) {
910  std::vector<int> ifi, iet, ily, idp;
911  for (int i = 0; i < nphi; ++i)
912  ifi.emplace_back(hpar->noff[noffk + 7 + i]);
913  for (int i = 0; i < ndeps; ++i) {
914  iet.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i]);
915  ily.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i + 1]);
916  idp.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i + 2]);
917  }
918 #ifdef EDM_ML_DEBUG
919  edm::LogVerbatim("HcalGeom") << "Initialize HcalLayerDepthMap for "
920  << "Detector " << dtype << " etaMax " << hpar->etaMax[dtype] << " with " << nphi
921  << " sectors";
922  for (int i = 0; i < nphi; ++i)
923  edm::LogVerbatim("HcalGeom") << " [" << i << "] " << ifi[i];
924  edm::LogVerbatim("HcalGeom") << "And " << ndeps << " depth sections";
925  for (int i = 0; i < ndeps; ++i)
926  edm::LogVerbatim("HcalGeom") << " [" << i << "] " << iet[i] << " " << ily[i] << " " << idp[i];
927  edm::LogVerbatim("HcalGeom") << "Maximum depth for last HE Eta tower " << depthEta29[0] << ":" << ndp16 << ":"
928  << ndp29 << " L0 Wt " << hpar->Layer0Wt[dtype - 1] << ":" << wt;
929 #endif
930  ldmap_.initialize(dtype, hpar->etaMax[dtype - 1], ndp16, ndp29, wt, ifi, iet, ily, idp);
931  int zside = (ifi[0] > 0) ? 1 : -1;
932  int iphi = (ifi[0] > 0) ? ifi[0] : -ifi[0];
933  depthMaxSp_ = std::pair<int, int>(dtype, ldmap_.getDepthMax(dtype, iphi, zside));
934  }
935  }
936  int noffm = (noffk + 7 + nphi + 3 * ndeps);
937  if ((int)(hpar->noff.size()) > noffm) {
938  int ndnext = hpar->noff[noffm];
939  if (ndnext > 4 && (int)(hpar->noff.size()) >= noffm + ndnext) {
940  for (int i = 0; i < 2; ++i)
941  layFHB[i] = hpar->noff[noffm + i + 1];
942  for (int i = 0; i < 3; ++i)
943  layFHE[i] = hpar->noff[noffm + i + 3];
944  }
945  if (ndnext > 11 && (int)(hpar->noff.size()) >= noffm + ndnext) {
946  for (int i = 0; i < 3; ++i)
947  layBHB[i] = hpar->noff[noffm + i + 6];
948  for (int i = 0; i < 4; ++i)
949  layBHE[i] = hpar->noff[noffm + i + 9];
950  }
951  }
952  }
953  }
954 #ifdef EDM_ML_DEBUG
955  edm::LogVerbatim("HcalGeom") << "Front Layer Definition for HB: " << layFHB[0] << ":" << layFHB[1]
956  << " and for HE: " << layFHE[0] << ":" << layFHE[1] << ":" << layFHE[2];
957  edm::LogVerbatim("HcalGeom") << "Last Layer Definition for HB: " << layBHB[0] << ":" << layBHB[1] << ":" << layBHB[2]
958  << " and for HE: " << layBHE[0] << ":" << layBHE[1] << ":" << layBHE[2] << ":"
959  << layBHE[3];
960 #endif
961  if (depthMaxSp_.first == 0) {
962  depthMaxSp_ = depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
963  } else if (depthMaxSp_.first == 1) {
964  depthMaxDf_ = std::pair<int, int>(1, maxDepth[0]);
965  if (depthMaxSp_.second > maxDepth[0])
966  maxDepth[0] = depthMaxSp_.second;
967  } else {
968  depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
969  if (depthMaxSp_.second > maxDepth[1])
970  maxDepth[1] = depthMaxSp_.second;
971  }
972 #ifdef EDM_ML_DEBUG
973  edm::LogVerbatim("HcalGeom") << "Detector type and maximum depth for all RBX " << depthMaxDf_.first << ":"
974  << depthMaxDf_.second << " and for special RBX " << depthMaxSp_.first << ":"
975  << depthMaxSp_.second;
976 #endif
977 }

References depthEta16, depthEta29, depthMaxDf_, depthMaxSp_, depths, dlShort, runTauDisplay::dtype, HcalParameters::dzVcal, dzVcal, etaHO, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaTable, cms::cuda::for(), HcalLayerDepthMap::getDepthMax(), getLayerMax(), HcalParameters::gparHF, HcalForward, hpar, mps_fire::i, idHF2QIE, HcalLayerDepthMap::initialize(), LEDCalibrationChannels::iphi, isBH_, dqmdumpme::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, layBHB, layBHE, HcalParameters::Layer0Wt, layerGroup(), layerGroupSize(), layFHB, layFHE, ldmap_, HcalParameters::maxDepth, maxDepth, maxLayerHB_, HcalParameters::modHB, HcalParameters::modHE, nEta, nmodHB, nmodHE, HcalParameters::noff, nphi, nPhiF, nR, nzHB, nzHE, HcalParameters::rHB, HcalParameters::rHO, rminHO, HcalParameters::rTable, HcalParameters::zHO, ecaldqm::zside(), and zVcal.

Referenced by HcalDDDSimConstants().

◆ isBH()

bool HcalDDDSimConstants::isBH ( ) const
inline

Definition at line 64 of file HcalDDDSimConstants.h.

64 { return isBH_; }

References isBH_.

Referenced by HcalDDDRecConstants::isBH().

◆ layerGroup() [1/2]

unsigned int HcalDDDSimConstants::layerGroup ( int  det,
int  eta,
int  phi,
int  zside,
int  i 
) const
private

Definition at line 1229 of file HcalDDDSimConstants.cc.

1229  {
1230  int depth0 = findDepth(det, eta, phi, zside, lay);
1231  unsigned int depth = (depth0 > 0) ? (unsigned int)(depth0) : layerGroup(eta - 1, lay);
1232  return depth;
1233 }

References LEDCalibrationChannels::depth, PVValHelper::eta, findDepth(), layerGroup(), phi, and ecaldqm::zside().

◆ layerGroup() [2/2]

unsigned int HcalDDDSimConstants::layerGroup ( int  eta,
int  i 
) const
private

Definition at line 1216 of file HcalDDDSimConstants.cc.

1216  {
1217  unsigned int k = 0;
1218  for (auto const& it : hpar->layerGroupEtaSim) {
1219  if (it.layer == (unsigned int)(eta + 1)) {
1220  return it.layerGroup.at(i);
1221  }
1222  if (it.layer > (unsigned int)(eta + 1))
1223  break;
1224  k = it.layerGroup.at(i);
1225  }
1226  return k;
1227 }

References PVValHelper::eta, hpar, mps_fire::i, dqmdumpme::k, and HcalParameters::layerGroupEtaSim.

Referenced by findLayer(), getEtaDepth(), getLayerFront(), getMaxDepth(), getMinDepth(), initialize(), and layerGroup().

◆ layerGroupSize()

unsigned int HcalDDDSimConstants::layerGroupSize ( int  eta) const
private

Definition at line 1203 of file HcalDDDSimConstants.cc.

1203  {
1204  unsigned int k = 0;
1205  for (auto const& it : hpar->layerGroupEtaSim) {
1206  if (it.layer == (unsigned int)(eta + 1)) {
1207  return it.layerGroup.size();
1208  }
1209  if (it.layer > (unsigned int)(eta + 1))
1210  break;
1211  k = it.layerGroup.size();
1212  }
1213  return k;
1214 }

References PVValHelper::eta, hpar, dqmdumpme::k, and HcalParameters::layerGroupEtaSim.

Referenced by getLayerMax(), getMinDepth(), and initialize().

◆ ldMap()

const HcalLayerDepthMap* HcalDDDSimConstants::ldMap ( ) const
inline

Definition at line 65 of file HcalDDDSimConstants.h.

65 { return &ldmap_; }

References ldmap_.

Referenced by HcalDDDRecConstants::getDepth(), HcalDDDRecConstants::getEtaBins(), HcalDDDRecConstants::getLayerBack(), HcalDDDRecConstants::getLayerFront(), HcalDDDRecConstants::getPhiZOne(), HcalDDDRecConstants::initialize(), and HcalDDDRecConstants::withSpecialRBXHBHE().

◆ maxHFDepth()

int HcalDDDSimConstants::maxHFDepth ( const int &  ieta,
const int &  iphi 
) const

Definition at line 655 of file HcalDDDSimConstants.cc.

655  {
656  int mxdepth = maxDepth[2];
657  if (!idHF2QIE.empty()) {
658  bool ok(false);
659  for (auto k : idHF2QIE) {
660  if ((eta == k.ieta()) && (iphi == k.iphi())) {
661  ok = true;
662  break;
663  }
664  }
665  if (!ok)
666  mxdepth = 2;
667  }
668  return mxdepth;
669 }

References PVValHelper::eta, idHF2QIE, LEDCalibrationChannels::iphi, dqmdumpme::k, maxDepth, and convertSQLiteXML::ok.

Referenced by getMaxDepth(), HcalDDDRecConstants::maxHFDepth(), FastHFShowerLibrary::modifyDepth(), and HCalSD::modifyDepth().

◆ numberOfCells()

unsigned int HcalDDDSimConstants::numberOfCells ( const HcalSubdetector &  subdet) const

Definition at line 671 of file HcalDDDSimConstants.cc.

671  {
672  unsigned int num = 0;
673  std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
674  for (auto& cellType : cellTypes) {
675  num += (unsigned int)(cellType.nPhiBins());
676  }
677 #ifdef EDM_ML_DEBUG
678  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants:numberOfCells " << cellTypes.size() << " " << num
679  << " for subdetector " << subdet;
680 #endif
681  return num;
682 }

References HcalCellTypes(), createfilelist::int, and EgammaValidation_cff::num.

Referenced by HcalDDDRecConstants::nCells(), and HcalDDDRecConstants::numberOfCells().

◆ parameter()

const HcalParameters* HcalDDDSimConstants::parameter ( ) const
inline

Definition at line 68 of file HcalDDDSimConstants.h.

68 { return hpar; }

References hpar.

◆ phiNumber()

int HcalDDDSimConstants::phiNumber ( const int &  phi,
const int &  unit 
) const

Definition at line 684 of file HcalDDDSimConstants.cc.

684  {
685  int iphi_skip = phi;
686  if (units == 2)
687  iphi_skip = (phi - 1) * 2 + 1;
688  else if (units == 4)
689  iphi_skip = (phi - 1) * 4 - 1;
690  if (iphi_skip < 0)
691  iphi_skip += 72;
692  return iphi_skip;
693 }

References phi, and units().

Referenced by HcalDDDRecConstants::getHCID(), getPhis(), HcalDDDRecConstants::getPhis(), and HcalNumberingFromDDD::unitID().

◆ printTileHB()

void HcalDDDSimConstants::printTileHB ( const int &  eta,
const int &  phi,
const int &  zside,
const int &  depth 
) const
private

Definition at line 1104 of file HcalDDDSimConstants.cc.

1104  {
1105  edm::LogVerbatim("HcalGeom") << "HcalDDDSimConstants::printTileHB for eta " << eta << " and depth " << depth;
1106 
1107  double etaL = hpar->etaTable.at(eta - 1);
1108  double thetaL = 2. * atan(exp(-etaL));
1109  double etaH = hpar->etaTable.at(eta);
1110  double thetaH = 2. * atan(exp(-etaH));
1111  int layL = getLayerFront(1, eta, phi, zside, depth);
1112  int layH = getLayerBack(1, eta, phi, zside, depth);
1113  edm::LogVerbatim("HcalGeom") << "\ntileHB:: eta|depth " << zside * eta << "|" << depth << " theta "
1114  << convertRadToDeg(thetaH) << ":" << convertRadToDeg(thetaL) << " Layer " << layL - 1
1115  << ":" << layH - 1;
1116  for (int lay = layL - 1; lay < layH; ++lay) {
1117  std::vector<double> area(2, 0);
1118  int kk(0);
1119  double mean(0);
1120  for (unsigned int k = 0; k < hpar->layHB.size(); ++k) {
1121  if (lay == hpar->layHB[k]) {
1122  double zmin = hpar->rhoxHB[k] * std::cos(thetaL) / std::sin(thetaL);
1123  double zmax = hpar->rhoxHB[k] * std::cos(thetaH) / std::sin(thetaH);
1124  double dz = (std::min(zmax, hpar->dxHB[k]) - zmin);
1125  if (dz > 0) {
1126  area[kk] = dz * hpar->dyHB[k];
1127  mean += area[kk];
1128  kk++;
1129  }
1130  }
1131  }
1132  if (area[0] > 0) {
1133  mean /= (kk * 100);
1134  edm::LogVerbatim("HcalGeom") << std::setw(2) << lay << " Area " << std::setw(8) << area[0] << " " << std::setw(8)
1135  << area[1] << " Mean " << mean;
1136  }
1137  }
1138 }

References custom_jme_cff::area, angle_units::operators::convertRadToDeg(), funct::cos(), LEDCalibrationChannels::depth, HcalParameters::dxHB, HcalParameters::dyHB, PVValHelper::dz, PVValHelper::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, dqmdumpme::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, SiStripPI::mean, min(), phi, HcalParameters::rhoxHB, funct::sin(), SiStripMonitorCluster_cfi::zmax, SiStripMonitorCluster_cfi::zmin, and ecaldqm::zside().

Referenced by printTiles().

◆ printTileHE()

void HcalDDDSimConstants::printTileHE ( const int &  eta,
const int &  phi,
const int &  zside,
const int &  depth 
) const
private

Definition at line 1140 of file HcalDDDSimConstants.cc.

1140  {
1141  double etaL = hpar->etaTable[eta - 1];
1142  double thetaL = 2. * atan(exp(-etaL));
1143  double etaH = hpar->etaTable[eta];
1144  double thetaH = 2. * atan(exp(-etaH));
1145  int layL = getLayerFront(2, eta, phi, zside, depth);
1146  int layH = getLayerBack(2, eta, phi, zside, depth);
1147  double phib = hpar->phibin[eta - 1];
1148  int nphi = 2;
1149  if (phib > 6._deg)
1150  nphi = 1;
1151  edm::LogVerbatim("HcalGeom") << "\ntileHE:: Eta/depth " << zside * eta << "|" << depth << " theta "
1152  << convertRadToDeg(thetaH) << ":" << convertRadToDeg(thetaL) << " Layer " << layL - 1
1153  << ":" << layH - 1 << " phi " << nphi;
1154  for (int lay = layL - 1; lay < layH; ++lay) {
1155  std::vector<double> area(4, 0);
1156  int kk(0);
1157  double mean(0);
1158  for (unsigned int k = 0; k < hpar->layHE.size(); ++k) {
1159  if (lay == hpar->layHE[k]) {
1160  double rmin = hpar->zxHE[k] * std::tan(thetaH);
1161  double rmax = hpar->zxHE[k] * std::tan(thetaL);
1162  if ((lay != 0 || eta == 18) &&
1163  (lay != 1 || (eta == 18 && hpar->rhoxHE[k] - hpar->dyHE[k] > 1000) ||
1164  (eta != 18 && hpar->rhoxHE[k] - hpar->dyHE[k] < 1000)) &&
1165  rmin + 30 < hpar->rhoxHE[k] + hpar->dyHE[k] && rmax > hpar->rhoxHE[k] - hpar->dyHE[k]) {
1166  rmin = std::max(rmin, hpar->rhoxHE[k] - hpar->dyHE[k]);
1167  rmax = std::min(rmax, hpar->rhoxHE[k] + hpar->dyHE[k]);
1168  double dx1 = rmin * std::tan(phib);
1169  double dx2 = rmax * std::tan(phib);
1170  double ar1 = 0, ar2 = 0;
1171  if (nphi == 1) {
1172  ar1 = 0.5 * (rmax - rmin) * (dx1 + dx2 - 4. * hpar->dx1HE[k]);
1173  mean += ar1;
1174  } else {
1175  ar1 = 0.5 * (rmax - rmin) * (dx1 + dx2 - 2. * hpar->dx1HE[k]);
1176  ar2 = 0.5 * (rmax - rmin) * ((rmax + rmin) * tan(10._deg) - 4 * hpar->dx1HE[k]) - ar1;
1177  mean += (ar1 + ar2);
1178  }
1179  area[kk] = ar1;
1180  area[kk + 2] = ar2;
1181  kk++;
1182  }
1183  }
1184  }
1185  if (area[0] > 0 && area[1] > 0) {
1186  int lay0 = lay - 1;
1187  if (eta == 18)
1188  lay0++;
1189  if (nphi == 1) {
1190  mean /= (kk * 100);
1191  edm::LogVerbatim("HcalGeom") << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " "
1192  << std::setw(8) << area[1] << " Mean " << mean;
1193  } else {
1194  mean /= (kk * 200);
1195  edm::LogVerbatim("HcalGeom") << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " "
1196  << std::setw(8) << area[1] << ":" << std::setw(8) << area[2] << " " << std::setw(8)
1197  << area[3] << " Mean " << mean;
1198  }
1199  }
1200  }
1201 }

References custom_jme_cff::area, angle_units::operators::convertRadToDeg(), LEDCalibrationChannels::depth, HcalParameters::dx1HE, HcalParameters::dyHE, PVValHelper::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, dqmdumpme::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, SiStripPI::max, SiStripPI::mean, min(), nphi, phi, HcalParameters::phibin, HcalParameters::rhoxHE, funct::tan(), ecaldqm::zside(), and HcalParameters::zxHE.

Referenced by printTiles().

◆ printTiles()

void HcalDDDSimConstants::printTiles ( ) const

Definition at line 695 of file HcalDDDSimConstants.cc.

695  {
696  std::vector<int> phis;
697  int detsp = ldmap_.validDet(phis);
698  int kphi = (detsp > 0) ? phis[0] : 1;
699  int zside = (kphi > 0) ? 1 : -1;
700  int iphi = (kphi > 0) ? kphi : -kphi;
701  edm::LogVerbatim("HcalGeom") << "Tile Information for HB from " << hpar->etaMin[0] << " to " << hpar->etaMax[0]
702  << "\n";
703  for (int eta = hpar->etaMin[0]; eta <= hpar->etaMax[0]; eta++) {
704  int dmax = getMaxDepth(1, eta, iphi, -zside, false);
705  for (int depth = 1; depth <= dmax; depth++)
706  printTileHB(eta, iphi, -zside, depth);
707  if (detsp == 1) {
708  int dmax = getMaxDepth(1, eta, iphi, zside, false);
709  for (int depth = 1; depth <= dmax; depth++)
710  printTileHB(eta, iphi, zside, depth);
711  }
712  }
713 
714  edm::LogVerbatim("HcalGeom") << "\nTile Information for HE from " << hpar->etaMin[1] << " to " << hpar->etaMax[1]
715  << "\n";
716  for (int eta = hpar->etaMin[1]; eta <= hpar->etaMax[1]; eta++) {
717  int dmin = (eta == hpar->etaMin[1]) ? getDepthEta16(2, iphi, -zside) : 1;
718  int dmax = getMaxDepth(2, eta, iphi, -zside, false);
719  for (int depth = dmin; depth <= dmax; depth++)
720  printTileHE(eta, iphi, -zside, depth);
721  if (detsp == 2) {
722  int dmax = getMaxDepth(2, eta, iphi, zside, false);
723  for (int depth = 1; depth <= dmax; depth++)
724  printTileHE(eta, iphi, zside, depth);
725  }
726  }
727 }

References LEDCalibrationChannels::depth, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, getDepthEta16(), getMaxDepth(), hpar, LEDCalibrationChannels::iphi, ldmap_, printTileHB(), printTileHE(), HcalLayerDepthMap::validDet(), and ecaldqm::zside().

◆ unitPhi() [1/2]

int HcalDDDSimConstants::unitPhi ( const double &  dphi) const

Definition at line 735 of file HcalDDDSimConstants.cc.

735  {
736  const double fiveDegInRad = 2 * M_PI / 72;
737  int units = int(dphi / fiveDegInRad + 0.5);
738  if (units < 1)
739  units = 1;
740  return units;
741 }

References createfilelist::int, M_PI, and units().

◆ unitPhi() [2/2]

int HcalDDDSimConstants::unitPhi ( const int &  det,
const int &  etaR 
) const

Definition at line 729 of file HcalDDDSimConstants.cc.

729  {
730  double dphi =
731  (det == static_cast<int>(HcalForward)) ? hpar->phitable[etaR - hpar->etaMin[2]] : hpar->phibin[etaR - 1];
732  return unitPhi(dphi);
733 }

References HcalParameters::etaMin, HcalForward, hpar, HcalParameters::phibin, and HcalParameters::phitable.

Referenced by cell(), HcalDDDRecConstants::getEtaPhi(), HcalDDDRecConstants::getHCID(), getPhiCons(), getPhis(), HcalDDDRecConstants::getPhis(), HcalCellTypes(), HcalDDDRecConstants::HcalCellTypes(), HcalDDDRecConstants::initialize(), and HcalNumberingFromDDD::unitID().

Member Data Documentation

◆ depthEta16

int HcalDDDSimConstants::depthEta16[2]
private

Definition at line 106 of file HcalDDDSimConstants.h.

Referenced by getDepthEta16(), getDepthEta16M(), getDetEta(), and initialize().

◆ depthEta29

int HcalDDDSimConstants::depthEta29[2]
private

Definition at line 107 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getDepthEta29(), getDepthEta29M(), getEta(), and initialize().

◆ depthMaxDf_

std::pair<int, int> HcalDDDSimConstants::depthMaxDf_
private

Definition at line 114 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

◆ depthMaxSp_

std::pair<int, int> HcalDDDSimConstants::depthMaxSp_
private

Definition at line 114 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

◆ depths

std::vector<int> HcalDDDSimConstants::depths[nDepthMax]
private

Definition at line 98 of file HcalDDDSimConstants.h.

Referenced by getLayerBack(), and initialize().

◆ dlShort

double HcalDDDSimConstants::dlShort
private

Definition at line 105 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), HcalCellTypes(), and initialize().

◆ dzVcal

double HcalDDDSimConstants::dzVcal
private

Definition at line 104 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

◆ etaHO

double HcalDDDSimConstants::etaHO[4]
private

Definition at line 102 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), and initialize().

◆ hpar

const HcalParameters* HcalDDDSimConstants::hpar
private

Definition at line 91 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getConstHBHE(), getDetEta(), getEta(), getEtaDepth(), getEtaHO(), getEtaTableHF(), getGain(), getGparHF(), getiEtaRange(), getLayer0Wt(), getLayerFront(), getLayerMax(), getMaxDepth(), getMinDepth(), getPhiCons(), getPhiTableHF(), getRTableHF(), getShift(), HcalCellTypes(), initialize(), layerGroup(), layerGroupSize(), parameter(), printTileHB(), printTileHE(), printTiles(), and unitPhi().

◆ idHF2QIE

std::vector<HcalDetId> HcalDDDSimConstants::idHF2QIE
private

Definition at line 113 of file HcalDDDSimConstants.h.

Referenced by getIdHF2QIE(), HcalCellTypes(), initialize(), and maxHFDepth().

◆ isBH_

bool HcalDDDSimConstants::isBH_
private

Definition at line 112 of file HcalDDDSimConstants.h.

Referenced by initialize(), and isBH().

◆ layBHB

int HcalDDDSimConstants::layBHB[3]
private

Definition at line 110 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

◆ layBHE

int HcalDDDSimConstants::layBHE[4]
private

Definition at line 111 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

◆ layFHB

int HcalDDDSimConstants::layFHB[2]
private

Definition at line 108 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

◆ layFHE

int HcalDDDSimConstants::layFHE[3]
private

Definition at line 109 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

◆ ldmap_

HcalLayerDepthMap HcalDDDSimConstants::ldmap_
private

Definition at line 92 of file HcalDDDSimConstants.h.

Referenced by findDepth(), getDepthEta16(), getDepthEta16M(), getDepthEta29(), getDepthEta29M(), getLayer0Wt(), getLayerBack(), getLayerFront(), getMaxDepth(), getMinDepth(), initialize(), ldMap(), and printTiles().

◆ maxDepth

std::vector<int> HcalDDDSimConstants::maxDepth
private

Definition at line 94 of file HcalDDDSimConstants.h.

Referenced by cell(), getLayerMax(), getMaxDepth(), getMinDepth(), HcalCellTypes(), initialize(), and maxHFDepth().

◆ maxLayer_

const int HcalDDDSimConstants::maxLayer_ = 18
staticprivate

Definition at line 76 of file HcalDDDSimConstants.h.

Referenced by getEta(), getLayerBack(), getLayerFront(), getMaxDepth(), and HcalCellTypes().

◆ maxLayerHB_

const int HcalDDDSimConstants::maxLayerHB_ = 16
staticprivate

Definition at line 77 of file HcalDDDSimConstants.h.

Referenced by cell(), getLayerFront(), getLayerMax(), getMaxDepth(), HcalCellTypes(), and initialize().

◆ nDepth

int HcalDDDSimConstants::nDepth
private

Definition at line 99 of file HcalDDDSimConstants.h.

◆ nDepthMax

const int HcalDDDSimConstants::nDepthMax = 9
staticprivate

Definition at line 75 of file HcalDDDSimConstants.h.

◆ nEta

int HcalDDDSimConstants::nEta
private

Definition at line 95 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getDetEta(), getEta(), and initialize().

◆ nmodHB

int HcalDDDSimConstants::nmodHB
private

Definition at line 100 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

◆ nmodHE

int HcalDDDSimConstants::nmodHE
private

Definition at line 101 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

◆ nPhiF

int HcalDDDSimConstants::nPhiF
private

Definition at line 97 of file HcalDDDSimConstants.h.

Referenced by initialize().

◆ nR

int HcalDDDSimConstants::nR
private

Definition at line 96 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getEta(), and initialize().

◆ nzHB

int HcalDDDSimConstants::nzHB
private

Definition at line 100 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

◆ nzHE

int HcalDDDSimConstants::nzHE
private

Definition at line 101 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), HcalCellTypes(), and initialize().

◆ rminHO

double HcalDDDSimConstants::rminHO
private

Definition at line 102 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

◆ zVcal

double HcalDDDSimConstants::zVcal
private

Definition at line 103 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getDetEta(), getEta(), and initialize().

HcalDDDSimConstants::nmodHB
int nmodHB
Definition: HcalDDDSimConstants.h:100
HcalDDDSimConstants::idHF2QIE
std::vector< HcalDetId > idHF2QIE
Definition: HcalDDDSimConstants.h:113
DDAxes::y
HcalParameters::HFGains
std::vector< double > HFGains
Definition: HcalParameters.h:44
HcalDDDSimConstants::layBHB
int layBHB[3]
Definition: HcalDDDSimConstants.h:110
ALCARECOTkAlBeamHalo_cff.etaMin
etaMin
GeV.
Definition: ALCARECOTkAlBeamHalo_cff.py:32
mps_fire.i
i
Definition: mps_fire.py:428
HcalParameters::maxDepth
std::vector< int > maxDepth
Definition: HcalParameters.h:48
geometryCSVtoXML.zz
zz
Definition: geometryCSVtoXML.py:19
HcalDDDSimConstants::getDepthEta16
int getDepthEta16(const int &det, const int &phi, const int &zside) const
Definition: HcalDDDSimConstants.cc:152
HcalLayerDepthMap::getLayer0Wt
double getLayer0Wt(const int subdet, const int iphi, const int zside) const
Definition: HcalLayerDepthMap.cc:242
HcalParameters::layHB
std::vector< int > layHB
Definition: HcalParameters.h:51
SiStripPI::mean
Definition: SiStripPayloadInspectorHelper.h:169
generateTowerEtThresholdLUT.etaRange
etaRange
Definition: generateTowerEtThresholdLUT.py:55
HcalParameters::etaTable
std::vector< double > etaTable
Definition: HcalParameters.h:35
HcalDDDSimConstants::getEta
int getEta(const int &det, const int &lay, const double &hetaR) const
Definition: HcalDDDSimConstants.cc:225
min
T min(T a, T b)
Definition: MathUtil.h:58
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:189
angle_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
HcalLayerDepthMap::getLayerFront
int getLayerFront(const int subdet, const int ieta, const int iphi, const int zside, const int depth) const
Definition: HcalLayerDepthMap.cc:168
HcalParameters::HBGains
std::vector< double > HBGains
Definition: HcalParameters.h:42
HcalDDDSimConstants::nEta
int nEta
Definition: HcalDDDSimConstants.h:95
HcalLayerDepthMap::isValid
bool isValid(const int det, const int phi, const int zside) const
Definition: HcalLayerDepthMap.cc:250
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
HcalLayerDepthMap::getDepth
int getDepth(const int subdet, const int ieta, const int iphi, const int zside, const int layer) const
Definition: HcalLayerDepthMap.cc:111
HcalParameters::etaTableHF
std::vector< double > etaTableHF
Definition: HcalParameters.h:45
HcalDDDSimConstants::maxDepth
std::vector< int > maxDepth
Definition: HcalDDDSimConstants.h:94
HcalDDDSimConstants::rminHO
double rminHO
Definition: HcalDDDSimConstants.h:102
DDAxes::x
HcalParameters::zHO
std::vector< double > zHO
Definition: HcalParameters.h:21
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
HcalBarrel
Definition: HcalAssistant.h:33
HcalDDDSimConstants::HcalCellTypes
std::vector< HcalCellType > HcalCellTypes() const
Definition: HcalDDDSimConstants.cc:505
createJobs.tmp
tmp
align.sh
Definition: createJobs.py:716
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
SiStripMonitorCluster_cfi.zmin
zmin
Definition: SiStripMonitorCluster_cfi.py:200
HcalDDDSimConstants::ldmap_
HcalLayerDepthMap ldmap_
Definition: HcalDDDSimConstants.h:92
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HcalLayerDepthMap::initialize
void initialize(const int subdet, const int ietaMax, const int dep16C, const int dep29C, const double wtl0C, std::vector< int > const &iphi, std::vector< int > const &ieta, std::vector< int > const &layer, std::vector< int > const &depth)
Definition: HcalLayerDepthMap.cc:20
HcalParameters::layerGroupEtaSim
std::vector< LayerItem > layerGroupEtaSim
Definition: HcalParameters.h:63
HcalDDDSimConstants::cell
HcalCellType::HcalCell cell(const int &det, const int &zside, const int &depth, const int &etaR, const int &iphi) const
Definition: HcalDDDSimConstants.cc:28
HcalDDDSimConstants::findDepth
int findDepth(const int &det, const int &eta, const int &phi, const int &zside, const int &lay) const
Definition: HcalDDDSimConstants.cc:120
HcalLayerDepthMap::getMaxDepthLastHE
int getMaxDepthLastHE(const int subdet, const int iphi, const int zside) const
Definition: HcalLayerDepthMap.cc:234
HcalParameters::etaMax
std::vector< int > etaMax
Definition: HcalParameters.h:55
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
HcalDDDSimConstants::depthEta29
int depthEta29[2]
Definition: HcalDDDSimConstants.h:107
units
TString units(TString variable, Char_t axis)
HcalDDDSimConstants::nzHB
int nzHB
Definition: HcalDDDSimConstants.h:100
LEDCalibrationChannels.iphi
iphi
Definition: LEDCalibrationChannels.py:64
HcalParameters::layHE
std::vector< int > layHE
Definition: HcalParameters.h:52
nphi
const int nphi
Definition: CMTRawAnalyzer.h:424
HcalDDDSimConstants::layFHE
int layFHE[3]
Definition: HcalDDDSimConstants.h:109
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
trackingPlots.hp
hp
Definition: trackingPlots.py:1231
HcalDDDSimConstants::dlShort
double dlShort
Definition: HcalDDDSimConstants.h:105
HcalParameters::dx1HE
std::vector< double > dx1HE
Definition: HcalParameters.h:30
HcalParameters::rhoxHE
std::vector< double > rhoxHE
Definition: HcalParameters.h:27
SiStripMonitorCluster_cfi.zmax
zmax
Definition: SiStripMonitorCluster_cfi.py:201
HcalDDDSimConstants::getGain
double getGain(const HcalSubdetector &subdet, const int &depth) const
Definition: HcalDDDSimConstants.cc:1085
HcalParameters::rHO
std::vector< double > rHO
Definition: HcalParameters.h:32
PVValHelper::eta
Definition: PVValidationHelpers.h:70
HcalParameters::dxHB
std::vector< double > dxHB
Definition: HcalParameters.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
DDAxes::z
HcalDDDSimConstants::maxLayerHB_
static const int maxLayerHB_
Definition: HcalDDDSimConstants.h:77
HcalParameters::phitable
std::vector< double > phitable
Definition: HcalParameters.h:38
HcalOuter
Definition: HcalAssistant.h:35
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:150
GetRecoTauVFromDQM_MC_cff.kk
kk
Definition: GetRecoTauVFromDQM_MC_cff.py:84
HcalParameters::dzHE
std::vector< double > dzHE
Definition: HcalParameters.h:20
dqmdumpme.k
k
Definition: dqmdumpme.py:60
HcalDDDSimConstants::maxHFDepth
int maxHFDepth(const int &ieta, const int &iphi) const
Definition: HcalDDDSimConstants.cc:655
HcalParameters::rHB
std::vector< double > rHB
Definition: HcalParameters.h:17
LEDCalibrationChannels.depth
depth
Definition: LEDCalibrationChannels.py:65
HcalDDDSimConstants::unitPhi
int unitPhi(const int &det, const int &etaR) const
Definition: HcalDDDSimConstants.cc:729
HcalParameters::dzVcal
double dzVcal
Definition: HcalParameters.h:46
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
HcalDDDSimConstants::layerGroupSize
unsigned int layerGroupSize(int eta) const
Definition: HcalDDDSimConstants.cc:1203
HcalParameters::modHB
std::vector< int > modHB
Definition: HcalParameters.h:49
HcalLayerDepthMap::getDepths
std::pair< int, int > getDepths(const int eta) const
Definition: HcalLayerDepthMap.cc:160
HcalParameters::modHE
std::vector< int > modHE
Definition: HcalParameters.h:50
HcalParameters::rTable
std::vector< double > rTable
Definition: HcalParameters.h:36
LEDCalibrationChannels.ieta
ieta
Definition: LEDCalibrationChannels.py:63
HcalDDDSimConstants::layerGroup
unsigned int layerGroup(int eta, int i) const
Definition: HcalDDDSimConstants.cc:1216
HcalDDDSimConstants::layFHB
int layFHB[2]
Definition: HcalDDDSimConstants.h:108
HcalDDDSimConstants::depthMaxSp_
std::pair< int, int > depthMaxSp_
Definition: HcalDDDSimConstants.h:114
HcalLayerDepthMap::getLayerBack
int getLayerBack(const int subdet, const int ieta, const int iphi, const int zside, const int depth) const
Definition: HcalLayerDepthMap.cc:183
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
HcalParameters::zxHE
std::vector< double > zxHE
Definition: HcalParameters.h:28
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:40
HcalDDDSimConstants::dzVcal
double dzVcal
Definition: HcalDDDSimConstants.h:104
runTauDisplay.dtype
dtype
Definition: runTauDisplay.py:121
HcalDDDSimConstants::depths
std::vector< int > depths[nDepthMax]
Definition: HcalDDDSimConstants.h:98
HcalDDDSimConstants::hpar
const HcalParameters * hpar
Definition: HcalDDDSimConstants.h:91
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
HcalParameters::zHE
std::vector< double > zHE
Definition: HcalParameters.h:19
HcalDetId
Definition: HcalDetId.h:12
createfilelist.int
int
Definition: createfilelist.py:10
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
HcalDDDSimConstants::layBHE
int layBHE[4]
Definition: HcalDDDSimConstants.h:111
HcalDDDSimConstants::phiNumber
int phiNumber(const int &phi, const int &unit) const
Definition: HcalDDDSimConstants.cc:684
HcalCellType::HcalCell
Definition: HcalCellType.h:15
HcalParameters::phioff
std::vector< double > phioff
Definition: HcalParameters.h:34
HcalParameters::dyHE
std::vector< double > dyHE
Definition: HcalParameters.h:29
HcalDDDSimConstants::etaHO
double etaHO[4]
Definition: HcalDDDSimConstants.h:102
cms::cuda::for
for(int i=first, nt=offsets[nh];i< nt;i+=gridDim.x *blockDim.x)
Definition: HistoContainer.h:15
EgammaValidation_cff.num
num
Definition: EgammaValidation_cff.py:33
HcalDDDSimConstants::getDepthEta29
int getDepthEta29(const int &phi, const int &zside, const int &i) const
Definition: HcalDDDSimConstants.cc:178
cmsLHEtoEOSManager.l
l
Definition: cmsLHEtoEOSManager.py:204
unit
Basic3DVector unit() const
Definition: Basic3DVectorLD.h:162
HcalDDDSimConstants::getLayerBack
int getLayerBack(const int &det, const int &eta, const int &phi, const int &zside, const int &depth) const
Definition: HcalDDDSimConstants.cc:387
alignCSCRings.r
r
Definition: alignCSCRings.py:93
PedestalClient_cfi.gain
gain
Definition: PedestalClient_cfi.py:37
HcalCellType
Definition: HcalCellType.h:13
HcalForward
Definition: HcalAssistant.h:36
HcalLayerDepthMap::getDepth16
int getDepth16(const int subdet, const int iphi, const int zside) const
Definition: HcalLayerDepthMap.cc:126
DDAxes::phi
HcalParameters::dyHB
std::vector< double > dyHB
Definition: HcalParameters.h:25
HcalParameters::HFShift
std::vector< int > HFShift
Definition: HcalParameters.h:59
HcalDDDSimConstants::depthMaxDf_
std::pair< int, int > depthMaxDf_
Definition: HcalDDDSimConstants.h:114
HcalDDDSimConstants::getLayerFront
int getLayerFront(const int &det, const int &eta, const int &phi, const int &zside, const int &depth) const
Definition: HcalDDDSimConstants.cc:368
HcalDDDSimConstants::zVcal
double zVcal
Definition: HcalDDDSimConstants.h:103
HcalDDDSimConstants::printTileHE
void printTileHE(const int &eta, const int &phi, const int &zside, const int &depth) const
Definition: HcalDDDSimConstants.cc:1140
edm::shift
static unsigned const int shift
Definition: LuminosityBlockID.cc:7
HcalParameters::Layer0Wt
std::vector< double > Layer0Wt
Definition: HcalParameters.h:41
HcalEndcap
Definition: HcalAssistant.h:34
HcalDDDSimConstants::depthEta16
int depthEta16[2]
Definition: HcalDDDSimConstants.h:106
PVValHelper::dz
Definition: PVValidationHelpers.h:51
HcalParameters::noff
std::vector< int > noff
Definition: HcalParameters.h:56
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
HcalDDDSimConstants::getPhiCons
std::pair< double, double > getPhiCons(const int &det, const int &ieta) const
Definition: HcalDDDSimConstants.cc:465
HcalLayerDepthMap::validDet
int validDet(std::vector< int > &phis) const
Definition: HcalLayerDepthMap.cc:258
HcalParameters::drHB
std::vector< double > drHB
Definition: HcalParameters.h:18
HcalParameters::HBShift
std::vector< int > HBShift
Definition: HcalParameters.h:57
ALCARECOTkAlBeamHalo_cff.etaMax
etaMax
Definition: ALCARECOTkAlBeamHalo_cff.py:33
HcalDDDSimConstants::printTileHB
void printTileHB(const int &eta, const int &phi, const int &zside, const int &depth) const
Definition: HcalDDDSimConstants.cc:1104
HcalDDDSimConstants::getPhis
std::vector< std::pair< int, double > > getPhis(const int &subdet, const int &ieta) const
Definition: HcalDDDSimConstants.cc:484
HcalDDDSimConstants::getLayerMax
int getLayerMax(const int &eta, const int &depth) const
Definition: HcalDDDSimConstants.cc:400
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
HcalDDDSimConstants::getMaxDepth
int getMaxDepth(const int &type) const
Definition: HcalDDDSimConstants.h:53
HcalDDDSimConstants::getShift
int getShift(const HcalSubdetector &subdet, const int &depth) const
Definition: HcalDDDSimConstants.cc:1066
HcalDDDSimConstants::nPhiF
int nPhiF
Definition: HcalDDDSimConstants.h:97
HcalParameters::HEShift
std::vector< int > HEShift
Definition: HcalParameters.h:58
HcalDDDSimConstants::nzHE
int nzHE
Definition: HcalDDDSimConstants.h:101
HcalParameters::phibin
std::vector< double > phibin
Definition: HcalParameters.h:37
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HcalParameters::HEGains
std::vector< double > HEGains
Definition: HcalParameters.h:43
HcalParameters::etaMin
std::vector< int > etaMin
Definition: HcalParameters.h:54
HcalCellType::HcalCell::ok
bool ok
Definition: HcalCellType.h:16
HcalParameters::gparHF
std::vector< double > gparHF
Definition: HcalParameters.h:40
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
HcalDDDSimConstants::isBH_
bool isBH_
Definition: HcalDDDSimConstants.h:112
HcalDDDSimConstants::maxLayer_
static const int maxLayer_
Definition: HcalDDDSimConstants.h:76
edm::Log
Definition: MessageLogger.h:70
HcalDDDSimConstants::initialize
void initialize()
Definition: HcalDDDSimConstants.cc:743
HcalDDDSimConstants::deltaEta
double deltaEta(const int &det, const int &eta, const int &depth) const
Definition: HcalDDDSimConstants.cc:979
HcalDDDSimConstants::nR
int nR
Definition: HcalDDDSimConstants.h:96
custom_jme_cff.area
area
Definition: custom_jme_cff.py:140
HcalLayerDepthMap::getDepthMax
int getDepthMax(const int subdet, const int iphi, const int zside) const
Definition: HcalLayerDepthMap.cc:144
d1
static constexpr float d1
Definition: L1EGammaCrystalsEmulatorProducer.cc:85
HcalParameters::rhoxHB
std::vector< double > rhoxHB
Definition: HcalParameters.h:23
HcalDDDSimConstants::nmodHE
int nmodHE
Definition: HcalDDDSimConstants.h:101
Generated for CMSSW Reference Manual by doxygen 1.8.16