CMS 3D CMS Logo
Functions
DDHGCalWaferFullRotated.cc File Reference
#include "DD4hep/DetFactoryHelper.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "DetectorDescription/DDCMS/interface/DDutils.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/HGCalCommonData/interface/HGCalTypes.h"
#include "Geometry/HGCalCommonData/interface/HGCalCell.h"
#include <string>
#include <vector>
#include <sstream>

Go to the source code of this file.

Functions

static long algorithm (dd4hep::Detector &, cms::DDParsingContext &ctxt, xml_h e)
 
 DD4HEP_OPEN_PLUGIN (dd4hep, ddcms_det_element_DDCMS_hgcal_DDHGCalWaferFullRotated)
 

Function Documentation

◆ algorithm()

static long algorithm ( dd4hep::Detector &  ,
cms::DDParsingContext &  ctxt,
xml_h  e 
)
static

Definition at line 19 of file DDHGCalWaferFullRotated.cc.

References funct::abs(), cms::DDNamespace::addSolidNS(), cms::DDNamespace::addVolumeNS(), writedatasetfile::args, hgcalTopologyTester_cfi::cell1, HGCalCell::cellPlacementIndex(), HGCalCell::cellUV2Cell(), HGCalCell::cellUV2XY1(), filterCSVwithJSON::copy, MillePedeFileConverter_cfg::e, HGCalGeometryMode::ExtrudedPolygon, mps_fire::i, cuy::ii, dqmdumpme::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, phase1PixelTopology::layerName, hgcalTopologyTester_cfi::layers, cms::DDNamespace::material(), g4SimHits_cfi::Material, cms::DDNamespace::name(), nCells, hgcalPerformanceValidation::orient, HGCalTypes::packCellTypeUV(), PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), diffTwoXMLs::r2, cms::s_executed, L1EGammaClusterEmuProducer_cfi::scale, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, findQualityFiles::v, cms::DDNamespace::volume(), HGCalTypes::waferFrontBack(), HGCalProperty::waferThick(), and zw().

19  {
20 #ifdef EDM_ML_DEBUG
21  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: Creating an instance";
22 #endif
23  cms::DDNamespace ns(ctxt, e, true);
24  cms::DDAlgoArguments args(ctxt, e);
25 
26  const auto& material = args.value<std::string>("ModuleMaterial");
27  const auto& thick = args.value<double>("ModuleThickness");
28  const auto& waferSize = args.value<double>("WaferSize");
29  const auto& waferSepar = args.value<double>("SensorSeparation");
30  const auto& waferThick = args.value<double>("WaferThickness");
31  const auto& waferTag = args.value<std::string>("WaferTag");
32 #ifdef EDM_ML_DEBUG
33  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: Module " << args.parentName() << " made of " << material
34  << " T " << thick << " Wafer 2r " << waferSize << " Half Separation " << waferSepar
35  << " T " << waferThick;
36 #endif
37  const auto& orient = args.value<std::vector<int> >("WaferOrient");
38  const auto& face = args.value<std::vector<int> >("WaferFace");
39  const auto& tag = args.value<std::vector<std::string> >("WaferPlacementIndex");
40  const auto& layerNames = args.value<std::vector<std::string> >("LayerNames");
41  const auto& materials = args.value<std::vector<std::string> >("LayerMaterials");
42  const auto& layerThick = args.value<std::vector<double> >("LayerThickness");
43  const auto& layerType = args.value<std::vector<int> >("LayerTypes");
44  std::vector<int> copyNumber;
45  copyNumber.resize(materials.size(), 1);
46 #ifdef EDM_ML_DEBUG
47  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: " << layerNames.size() << " types of volumes";
48  for (unsigned int i = 0; i < layerNames.size(); ++i)
49  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << layerNames[i] << " of thickness " << layerThick[i]
50  << " filled with " << materials[i] << " type " << layerType[i];
51 #endif
52  const auto& layers = args.value<std::vector<int> >("Layers");
53 #ifdef EDM_ML_DEBUG
54  std::ostringstream st1;
55  for (unsigned int i = 0; i < layers.size(); ++i)
56  st1 << " [" << i << "] " << layers[i];
57  edm::LogVerbatim("HGCalGeom") << "There are " << layers.size() << " blocks" << st1.str();
58 #endif
59  const auto& nCells = args.value<int>("NCells");
60  int cellType(-1);
61  std::vector<std::string> cellNames;
62  std::vector<int> cellOffset;
63  if (nCells > 0) {
64  cellType = args.value<int>("CellType");
65  cellNames = args.value<std::vector<std::string> >("CellNames");
66  cellOffset = args.value<std::vector<int> >("CellOffset");
67  }
68 #ifdef EDM_ML_DEBUG
69  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: Cells/Wafer " << nCells << " Cell Type " << cellType
70  << " NameSpace " << ns.name() << " # of cells " << cellNames.size();
71  std::ostringstream st2;
72  for (unsigned int i = 0; i < cellOffset.size(); ++i)
73  st2 << " [" << i << "] " << cellOffset[i];
74  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: " << cellOffset.size() << " types of cells with offsets "
75  << st2.str();
76  for (unsigned int k = 0; k < cellNames.size(); ++k)
77  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: Cell[" << k << "] " << cellNames[k];
78 
79  edm::LogVerbatim("HGCalGeom") << "==>> Executing DDHGCalWaferFullRotated...";
80 #endif
81 
82  static constexpr double tol = 0.00001;
83  static const double sqrt3 = std::sqrt(3.0);
84  double rM = 0.5 * waferSize;
85  double RM2 = rM / sqrt3;
86  double r2 = 0.5 * waferSize;
87  double R2 = r2 / sqrt3;
88  const int nFine(nCells), nCoarse(nCells);
89  HGCalCell wafer((waferSize + waferSepar), nFine, nCoarse);
90  for (unsigned int k = 0; k < tag.size(); ++k) {
91  // First the mother
92  std::vector<double> xM = {rM, 0, -rM, -rM, 0, rM};
93  std::vector<double> yM = {RM2, 2 * RM2, RM2, -RM2, -2 * RM2, -RM2};
94  std::vector<double> zw = {-0.5 * thick, 0.5 * thick};
95  std::vector<double> zx(2, 0), zy(2, 0), scale(2, 1.0);
96  std::string parentName = args.parentName() + tag[k] + waferTag;
97  dd4hep::Material matter = ns.material(material);
98  dd4hep::Solid solid = dd4hep::ExtrudedPolygon(xM, yM, zw, zx, zy, scale);
99  ns.addSolidNS(ns.prepend(parentName), solid);
100  dd4hep::Volume glogM = dd4hep::Volume(solid.name(), solid, matter);
101  ns.addVolumeNS(glogM);
102 #ifdef EDM_ML_DEBUG
103  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: " << solid.name() << " extruded polygon made of "
104  << material << " z|x|y|s (0) " << zw[0] << ":" << zx[0] << ":" << zy[0] << ":"
105  << scale[0] << " z|x|y|s (1) " << zw[1] << ":" << zx[1] << ":" << zy[1] << ":"
106  << scale[1] << " and " << xM.size() << " edges";
107  for (unsigned int kk = 0; kk < xM.size(); ++kk)
108  edm::LogVerbatim("HGCalGeom") << "[" << kk << "] " << xM[kk] << ":" << yM[kk];
109 #endif
110 
111  // Then the layers
112  std::vector<double> xL = {r2, 0, -r2, -r2, 0, r2};
113  std::vector<double> yL = {R2, 2 * R2, R2, -R2, -2 * R2, -R2};
114  std::vector<dd4hep::Volume> glogs(materials.size());
115  for (unsigned int ii = 0; ii < copyNumber.size(); ii++) {
116  copyNumber[ii] = 1;
117  }
118  double zi(-0.5 * thick), thickTot(0.0);
119  for (unsigned int l = 0; l < layers.size(); l++) {
120  unsigned int i = layers[l];
121  if (copyNumber[i] == 1) {
122  if (layerType[i] > 0) {
123  zw[0] = -0.5 * waferThick;
124  zw[1] = 0.5 * waferThick;
125  } else {
126  zw[0] = -0.5 * layerThick[i];
127  zw[1] = 0.5 * layerThick[i];
128  }
129  std::string layerName = layerNames[i] + tag[k] + waferTag;
130  solid = dd4hep::ExtrudedPolygon(xL, yL, zw, zx, zy, scale);
131  ns.addSolidNS(ns.prepend(layerName), solid);
132  matter = ns.material(materials[i]);
133  glogs[i] = dd4hep::Volume(solid.name(), solid, matter);
134  ns.addVolumeNS(glogs[i]);
135 #ifdef EDM_ML_DEBUG
136  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: " << solid.name() << " extruded polygon made of "
137  << materials[i] << " z|x|y|s (0) " << zw[0] << ":" << zx[0] << ":" << zy[0] << ":"
138  << scale[0] << " z|x|y|s (1) " << zw[1] << ":" << zx[1] << ":" << zy[1] << ":"
139  << scale[1] << " and " << xL.size() << " edges";
140  for (unsigned int kk = 0; kk < xL.size(); ++kk)
141  edm::LogVerbatim("HGCalGeom") << "[" << kk << "] " << xL[kk] << ":" << yL[kk];
142 #endif
143  }
144  dd4hep::Position tran0(0, 0, (zi + 0.5 * layerThick[i]));
145  glogM.placeVolume(glogs[i], copyNumber[i], tran0);
146 #ifdef EDM_ML_DEBUG
147  edm::LogVerbatim("HGCalGeom") << "DDHGCalWaferFullRotated: " << glogs[i].name() << " number " << copyNumber[i]
148  << " positioned in " << glogM.name() << " at " << tran0 << " with no rotation";
149 #endif
150  ++copyNumber[i];
151  zi += layerThick[i];
152  thickTot += layerThick[i];
153  if ((layerType[i] > 0) && (nCells > 0)) {
154  for (int u = 0; u < 2 * nCells; ++u) {
155  for (int v = 0; v < 2 * nCells; ++v) {
156  if (((v - u) < nCells) && ((u - v) <= nCells)) {
157  int placeIndex = wafer.cellPlacementIndex(1, HGCalTypes::waferFrontBack(face[k]), orient[k]);
158  std::pair<double, double> xy1 = wafer.cellUV2XY1(u, v, placeIndex, cellType);
159  double yp = xy1.second;
160  double xp = xy1.first;
161  int cell(0);
162  std::pair<int, int> cell1 = wafer.cellUV2Cell(u, v, placeIndex, cellType);
163  cell = cell1.first + cellOffset[cell1.second];
164  dd4hep::Position tran(xp, yp, 0);
165  int copy = HGCalTypes::packCellTypeUV(cellType, u, v);
166  glogs[i].placeVolume(ns.volume(cellNames[cell]), copy, tran);
167 #ifdef EDM_ML_DEBUG
168  edm::LogVerbatim("HGCalGeom")
169  << "DDHGCalWaferFullRotated: " << cellNames[cell] << " number " << copy << " positioned in "
170  << glogs[i].name() << " at " << tran << " with no rotation";
171 #endif
172  }
173  }
174  }
175  }
176  }
177  if (std::abs(thickTot - thick) >= tol) {
178  if (thickTot > thick) {
179  edm::LogError("HGCalGeom") << "Thickness of the partition " << thick << " is smaller than " << thickTot
180  << ": thickness of all its components **** ERROR ****";
181  } else {
182  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << thick << " does not match with " << thickTot
183  << " of the components";
184  }
185  }
186  }
187  return cms::s_executed;
188 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
mps_fire.i
i
Definition: mps_fire.py:429
HGCalTypes::packCellTypeUV
static int32_t packCellTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:28
phase1PixelTopology::layerName
constexpr char const * layerName[numberOfLayers]
Definition: SimplePixelTopology.h:86
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
nCells
caConstants::TupleMultiplicity const CAHitNtupletGeneratorKernelsGPU::HitToTuple const cms::cuda::AtomicPairCounter GPUCACell const *__restrict__ uint32_t const *__restrict__ nCells
Definition: CAHitNtupletGeneratorKernelsImpl.h:43
g4SimHits_cfi.Material
Material
Definition: g4SimHits_cfi.py:603
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
cms::s_executed
static constexpr long s_executed
Definition: DDAlgoArguments.h:16
makeGlobalPositionRcd_cfg.tag
tag
Definition: makeGlobalPositionRcd_cfg.py:6
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HGCalTypes::waferFrontBack
static int32_t waferFrontBack(int32_t index)
Definition: HGCalTypes.h:126
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:13
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
cuy.ii
ii
Definition: cuy.py:589
L1EGammaClusterEmuProducer_cfi.scale
scale
Definition: L1EGammaClusterEmuProducer_cfi.py:15
GetRecoTauVFromDQM_MC_cff.kk
kk
Definition: GetRecoTauVFromDQM_MC_cff.py:84
zw
auto zw(V v) -> Vec2< typename std::remove_reference< decltype(v[0])>::type >
Definition: ExtVec.h:75
hgcalTopologyTester_cfi.cell1
cell1
Definition: hgcalTopologyTester_cfi.py:11
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
hgcalTopologyTester_cfi.layers
layers
Definition: hgcalTopologyTester_cfi.py:8

◆ DD4HEP_OPEN_PLUGIN()

DD4HEP_OPEN_PLUGIN ( dd4hep  ,
ddcms_det_element_DDCMS_hgcal_DDHGCalWaferFullRotated   
)

Definition at line 190 of file DDHGCalWaferFullRotated.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.14