ATLAS Offline Software
BevelledCylinderVolumeBounds.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3 */
4 
6 // BevelledCylinderVolumeBounds.cxx, (c) ATLAS Detector software
8 
9 // Trk
10 #include "TrkVolumes/BevelledCylinderVolumeBounds.h"
11 // TrkSurfaces
12 #include "TrkSurfaces/CylinderSurface.h"
13 #include "TrkSurfaces/DiscSurface.h"
14 #include "TrkSurfaces/PlaneSurface.h"
15 //#include "TrkGeometrySurfaces/BevelledCylinderBounds.h"
16 #include "TrkGeometrySurfaces/SubtractedDiscSurface.h"
17 #include "TrkSurfaces/DiscBounds.h"
18 #include "TrkSurfaces/EllipseBounds.h"
19 #include "TrkSurfaces/RectangleBounds.h"
20 #include "TrkSurfaces/TrapezoidBounds.h"
21 #include "TrkVolumes/CombinedVolumeBounds.h"
22 #include "TrkVolumes/CuboidVolumeBounds.h"
23 #include "TrkVolumes/VolumeExcluder.h"
24 // Gaudi
25 #include "GaudiKernel/MsgStream.h"
26 #include "GaudiKernel/SystemOfUnits.h"
27 // STD
28 #include <cmath>
29 #include <iostream>
30 
31 const double Trk::BevelledCylinderVolumeBounds::s_numericalStable =
32  10e-2 * Gaudi::Units::millimeter;
33 
34 Trk::BevelledCylinderVolumeBounds::BevelledCylinderVolumeBounds()
35  : VolumeBounds()
36  , m_innerRadius(0.)
37  , m_outerRadius(0.)
38  , m_halfPhiSector(0.)
39  , m_halfZ(0.)
40  , m_thetaMinus(0.)
41  , m_thetaPlus(0.)
42  , m_type(-1)
43  , m_boundaryAccessors()
44  , m_subtractedVolume(nullptr)
45 {}
46 
47 Trk::BevelledCylinderVolumeBounds::BevelledCylinderVolumeBounds(
48  double rinner,
49  double router,
50  double haphi,
51  double halez,
52  int type)
53  : VolumeBounds()
54  , m_innerRadius(std::abs(rinner))
55  , m_outerRadius(std::abs(router))
56  , m_halfPhiSector(haphi)
57  , m_halfZ(std::abs(halez))
58  , m_thetaMinus(0.)
59  , m_thetaPlus(0.)
60  , m_type(type)
61  , m_boundaryAccessors()
62  , m_subtractedVolume(nullptr)
63 {}
64 
65 Trk::BevelledCylinderVolumeBounds::BevelledCylinderVolumeBounds(
66  const Trk::BevelledCylinderVolumeBounds& cylbo)
67  : VolumeBounds()
68  , m_innerRadius(cylbo.m_innerRadius)
69  , m_outerRadius(cylbo.m_outerRadius)
70  , m_halfPhiSector(cylbo.m_halfPhiSector)
71  , m_halfZ(cylbo.m_halfZ)
72  , m_thetaMinus(cylbo.m_thetaMinus)
73  , m_thetaPlus(cylbo.m_thetaPlus)
74  , m_type(cylbo.m_type)
75  , m_boundaryAccessors()
76  , m_subtractedVolume(nullptr)
77 {}
78 
79 Trk::BevelledCylinderVolumeBounds::~BevelledCylinderVolumeBounds() = default;
80 
81 Trk::BevelledCylinderVolumeBounds&
82 Trk::BevelledCylinderVolumeBounds::operator=(
83  const Trk::BevelledCylinderVolumeBounds& cylbo)
84 {
85  if (this != &cylbo) {
86  m_innerRadius = cylbo.m_innerRadius;
87  m_outerRadius = cylbo.m_outerRadius;
88  m_halfPhiSector = cylbo.m_halfPhiSector;
89  m_halfZ = cylbo.m_halfZ;
90  m_thetaMinus = cylbo.m_thetaMinus;
91  m_thetaPlus = cylbo.m_thetaPlus;
92  m_type = cylbo.m_type;
93  m_subtractedVolume.store(nullptr);
94  }
95  return *this;
96 }
97 
98 std::vector<std::unique_ptr<Trk::Surface>>
99  Trk::BevelledCylinderVolumeBounds::decomposeToSurfaces
100  (const Amg::Transform3D& transform)
101 {
102  auto retsf = std::vector<std::unique_ptr<Trk::Surface>>();
103  // memory optimisation (reserve a save number of 20)
104  retsf.reserve(6);
105 
106  Amg::RotationMatrix3D discRot(transform.rotation());
107  Amg::Vector3D cylCenter(transform.translation());
108 
109  // Lazy init m_m_subtractedVolume
110  if (m_type > -1 && !m_subtractedVolume) {
111  m_subtractedVolume.set(std::unique_ptr<Trk::Volume>(subtractedVolume()));
112  }
113 
114  // bottom Ellipse/Disc (negative z)
115  if (m_type < 0)
116  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
117  Amg::Transform3D(
118  (discRot *
119  Amg::AngleAxis3D(-m_thetaMinus + M_PI, Amg::Vector3D(0., 1., 0.))) *
120  Amg::Translation3D(
121  cylCenter - (halflengthZ() - m_outerRadius * tan(m_thetaMinus)) *
122  discRot.col(2))),
123  bottomEllipseBounds()));
124  else {
125  if (m_subtractedVolume) {
126  auto subtrVol = std::make_unique<Trk::Volume>(*m_subtractedVolume);
127  Trk::DiscSurface bottomDisc(
128  Amg::Transform3D(
129  transform * Amg::AngleAxis3D(M_PI, Amg::Vector3D(1., 0., 0.)) *
130  Amg::Translation3D(Amg::Vector3D(0., 0., halflengthZ()))),
131  discBounds());
132  retsf.push_back(std::make_unique<Trk::SubtractedDiscSurface>(
133  bottomDisc, std::make_shared<Trk::VolumeExcluder>(std::move(subtrVol)), false));
134  } else
135  retsf.push_back(std::make_unique<Trk::DiscSurface>(
136  Amg::Transform3D(
137  transform * Amg::AngleAxis3D(M_PI, Amg::Vector3D(1., 0., 0.)) *
138  Amg::Translation3D(Amg::Vector3D(0., 0., halflengthZ()))),
139  discBounds()));
140  }
141 
142  // top Ellipse/Disc (positive z)
143  if (m_type < 0)
144  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
145  Amg::Transform3D(
146  discRot * Amg::AngleAxis3D(m_thetaPlus, Amg::Vector3D(0., 1., 0.)) *
147  Amg::Translation3D(
148  cylCenter +
149  (halflengthZ() - m_outerRadius * tan(m_thetaPlus)) * discRot.col(2))),
150  topEllipseBounds()));
151  else {
152  if (m_subtractedVolume) {
153  auto subtrVol = std::make_unique<Trk::Volume>(*m_subtractedVolume);
154  Trk::DiscSurface topDisc(
155  Amg::Transform3D(
156  transform * Amg::Translation3D(Amg::Vector3D(0., 0., halflengthZ()))),
157  discBounds());
158  retsf.push_back(std::make_unique<Trk::SubtractedDiscSurface>(
159  topDisc, std::make_shared<Trk::VolumeExcluder>(std::move(subtrVol)), false));
160  } else
161  retsf.push_back(std::make_unique<Trk::DiscSurface>(
162  Amg::Transform3D(
163  transform * Amg::Translation3D(Amg::Vector3D(0., 0., halflengthZ()))),
164  discBounds()));
165  }
166 
167  // outer BevelledCylinder/Plane
168  if (m_type < 0)
169  retsf.push_back(std::make_unique<Trk::CylinderSurface>(
170  Amg::Transform3D(transform), outerBevelledCylinderBounds()));
171  else if (m_type < 2)
172  retsf.push_back(std::make_unique<Trk::CylinderSurface>(
173  Amg::Transform3D(transform), outerCylinderBounds()));
174  else
175  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
176  Amg::Transform3D(
177  transform *
178  Amg::Translation3D(Amg::Vector3D(this->outerRadius(), 0., 0.)) *
179  Amg::AngleAxis3D(-90 * Gaudi::Units::deg, Amg::Vector3D(0., 0., 1.)) *
180  Amg::AngleAxis3D(-90. * Gaudi::Units::deg, Amg::Vector3D(1., 0., 0.))),
181  outerBevelledPlaneBounds()));
182 
183  // inner BevelledCylinder/Plane
184  if (innerRadius() > s_numericalStable) {
185  if (m_type < 1)
186  retsf.push_back(std::make_unique<Trk::CylinderSurface>(
187  Amg::Transform3D(transform), innerBevelledCylinderBounds()));
188  else if (m_type == 2)
189  retsf.push_back(std::make_unique<Trk::CylinderSurface>(
190  Amg::Transform3D(transform), innerCylinderBounds()));
191  else
192  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
193  Amg::Transform3D(
194  transform *
195  Amg::Translation3D(Amg::Vector3D(this->innerRadius(), 0., 0.)) *
196  Amg::AngleAxis3D(+90 * Gaudi::Units::deg, Amg::Vector3D(0., 0., 1.)) *
197  Amg::AngleAxis3D(
198  -90. * Gaudi::Units::deg, Amg::Vector3D(1., 0., 0.))),
199  innerBevelledPlaneBounds()));
200  }
201 
202  if (std::abs(halfPhiSector() - M_PI) > s_numericalStable) {
203  if (m_type < 0) {
204  // sectorPlane 1 (negative phi)
205  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
206  Amg::Transform3D(
207  transform *
208  Amg::AngleAxis3D(-halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
209  Amg::Translation3D(Amg::Vector3D(mediumRadius(), 0., 0.)) *
210  Amg::AngleAxis3D(M_PI / 2, Amg::Vector3D(1., 0., 0.))),
211  sectorTrdBounds()));
212  // sectorPlane 2 (positive phi)
213  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
214  Amg::Transform3D(
215  transform *
216  Amg::AngleAxis3D(halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
217  Amg::Translation3D(Amg::Vector3D(mediumRadius(), 0., 0.)) *
218  Amg::AngleAxis3D(-M_PI / 2, Amg::Vector3D(1., 0., 0.))),
219  sectorTrdBounds()));
220  } else {
221  // sectorPlane 1 (negative phi)
222  double ri = innerRadius();
223  double ro = outerRadius();
224  if (m_type == 1 || m_type == 3)
225  ri *= 1. / cos(halfPhiSector());
226  if (m_type > 1)
227  ro *= 1. / cos(halfPhiSector());
228  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
229  Amg::Transform3D(
230  transform *
231  Amg::AngleAxis3D(-halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
232  Amg::Translation3D(Amg::Vector3D(0.5 * (ri + ro), 0., 0.)) *
233  Amg::AngleAxis3D(M_PI / 2, Amg::Vector3D(1., 0., 0.))),
234  sectorPlaneBounds()));
235  // sectorPlane 2 (positive phi)
236  retsf.push_back(std::make_unique<Trk::PlaneSurface>(
237  Amg::Transform3D(
238  transform *
239  Amg::AngleAxis3D(halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
240  Amg::Translation3D(Amg::Vector3D(0.5 * (ri + ro), 0., 0.)) *
241  Amg::AngleAxis3D(-M_PI / 2, Amg::Vector3D(1., 0., 0.))),
242  sectorPlaneBounds()));
243  }
244  }
245  return retsf;
246 }
247 
248 Trk::ObjectAccessor
249 Trk::BevelledCylinderVolumeBounds::boundarySurfaceAccessor(
250  const Amg::Vector3D& gp,
251  const Amg::Vector3D& dir,
252  bool) const
253 {
254  // the tube case - most likely
255  if (m_innerRadius != 0. && std::abs(m_halfPhiSector - M_PI) < s_numericalStable) {
256  // prepare the data
257  double posZ = gp.z();
258  double posR = gp.perp();
259  // difference Vector
260  Amg::Vector3D diff(gp + dir);
261  // differences
262  double deltaZ = diff.z() - posZ;
263  double deltaR = diff.perp() - posR;
264 
265  // the isOnSurface cases + switchers (that change with the cases)
266  bool isOnFace = false;
267  bool intersectionIndicator = (deltaR > 0.);
268  bool choiceIndicator = false;
269 
270  // on surface or look-a-likes
271  // on zMin or slightly outside
272  if (
273  std::abs(posZ + m_halfZ) < s_numericalStable ||
274  (posZ < (-m_halfZ) && deltaZ > 0.)) {
275  isOnFace = true;
276  } else if (
277  std::abs(posZ - m_halfZ) < s_numericalStable ||
278  (posZ > m_halfZ && deltaZ < 0.)) {
279  // on zMax or slighly outside
280  isOnFace = true;
281  } else if (
282  std::abs(posR - m_innerRadius) < s_numericalStable ||
283  (posR < m_innerRadius && deltaR > 0.)) {
284  // on innerRadius or slightly outside
285  isOnFace = true;
286  choiceIndicator =
287  (deltaZ > 0.); // the side choice indicator <========== WRONG ==========
288  } else if (
289  std::abs(posR - m_outerRadius) < s_numericalStable ||
290  (posR > m_outerRadius && deltaR < 0.)) {
291  // on outRadius or slightly outside
292  isOnFace = true;
293  choiceIndicator =
294  (deltaZ > 0.); // the side choice indicator <========== WRONG ==========
295  }
296 
297  // the onSurface case
298  // =================================================================================
299  if (isOnFace) {
300  if (intersectionIndicator) {
301  // intersect the Rmax
302  Trk::BevelledBoundaryIntersector intersectRmax(
303  posR, posZ, deltaZ / deltaR, m_outerRadius);
304  double zOfIntersect = intersectRmax.yOfX;
305  // now check if the intersect is inside m_halfZ
306  if (std::abs(zOfIntersect) <= m_halfZ)
307  return {(choiceIndicator || zOfIntersect > 0.)
308  ? m_boundaryAccessors.bevelledtubeAccessor(
309  Trk::BevelledTubeRincreaseZincrease)
310  : m_boundaryAccessors.bevelledtubeAccessor(
311  Trk::BevelledTubeRincreaseZdecrease)};
312  // if the intersect is outside
313  return {(choiceIndicator || zOfIntersect > 0.)
314  ? m_boundaryAccessors.bevelledtubeAccessor(
315  Trk::BevelledTubeZincreaseRincrease)
316  : m_boundaryAccessors.bevelledtubeAccessor(
317  Trk::BevelledTubeZdecreaseRincrease)};
318  }
319  // intersect the Rmin
320  Trk::BevelledBoundaryIntersector intersectRmin(
321  posR, posZ, deltaZ / deltaR, m_innerRadius);
322  double zOfIntersect = intersectRmin.yOfX;
323  // now check if the intersect is inside m_halfZ
324  if (std::abs(zOfIntersect) <= m_halfZ)
325  return {(choiceIndicator || zOfIntersect > 0.)
326  ? m_boundaryAccessors.bevelledtubeAccessor(
327  Trk::BevelledTubeRdecreaseZincrease)
328  : m_boundaryAccessors.bevelledtubeAccessor(
329  Trk::BevelledTubeRdecreaseZdecrease)};
330  // if the intersect is outside
331  return {(choiceIndicator || zOfIntersect > 0.)
332  ? m_boundaryAccessors.bevelledtubeAccessor(
333  Trk::BevelledTubeZincreaseRdecrease)
334  : m_boundaryAccessors.bevelledtubeAccessor(
335  Trk::BevelledTubeZdecreaseRdecrease)};
336  }
337  // =================================================================================================
338 
339  // ======================= the inside/outside part remains
340  // =========================================
341  // (a) outside cases
342  if (posR < m_innerRadius && deltaR < 0.)
343  return {m_boundaryAccessors.bevelledtubeAccessor(
344  Trk::BevelledTubeOutsideRminRdecrease)};
345  if (posR > m_outerRadius && deltaR > 0.)
346  return {m_boundaryAccessors.bevelledtubeAccessor(
347  Trk::BevelledTubeOutsideRmaxRincrease)};
348  if (posZ < -m_halfZ && deltaZ < 0.)
349  return {m_boundaryAccessors.bevelledtubeAccessor(
350  Trk::BevelledTubeOutsideZminZdecrease)};
351  if (posZ > m_halfZ && deltaZ > 0.)
352  return {m_boundaryAccessors.bevelledtubeAccessor(
353  Trk::BevelledTubeOutsideZmaxZincrease)};
354  // (b) inside cases
355  // the increase R case
356  if (deltaR > 0.) {
357  // solve the linear equation for the outer Radius
358  Trk::BevelledBoundaryIntersector intersectRmax(
359  posR, posZ, deltaZ / deltaR, m_outerRadius);
360  double zOfIntersect = intersectRmax.yOfX;
361 
362  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect > 0.)
363  return {m_boundaryAccessors.bevelledtubeAccessor(
364  Trk::BevelledTubeRincreaseZincrease)};
365  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect < 0.)
366  return {m_boundaryAccessors.bevelledtubeAccessor(
367  Trk::BevelledTubeRincreaseZdecrease)};
368  if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect < 0.)
369  return {m_boundaryAccessors.bevelledtubeAccessor(
370  Trk::BevelledTubeZdecreaseRincrease)};
371  if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect > 0.)
372  return {m_boundaryAccessors.bevelledtubeAccessor(
373  Trk::BevelledTubeZincreaseRincrease)};
374 
375  } else {
376  // solve the linear equation for the inner Radius
377  Trk::BevelledBoundaryIntersector intersectRmin(
378  posR, posZ, deltaZ / deltaR, m_innerRadius);
379  double zOfIntersect = intersectRmin.yOfX;
380 
381  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect > 0.)
382  return {m_boundaryAccessors.bevelledtubeAccessor(
383  Trk::BevelledTubeRdecreaseZincrease)};
384  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect < 0.)
385  return {m_boundaryAccessors.bevelledtubeAccessor(
386  Trk::BevelledTubeRdecreaseZdecrease)};
387  if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect > 0.)
388  return {m_boundaryAccessors.bevelledtubeAccessor(
389  Trk::BevelledTubeZincreaseRdecrease)};
390  return {m_boundaryAccessors.bevelledtubeAccessor(
391  Trk::BevelledTubeZdecreaseRdecrease)};
392  }
393  }
394  // the cylinder case
395  if (m_innerRadius == 0. && std::abs(m_halfPhiSector - M_PI) < 10e-3) {
396 
397  double posZ = gp.z();
398  double posR = gp.perp();
399  // difference Vector
400  Amg::Vector3D diff(gp + dir);
401  // differences
402  double deltaZ = diff.z() - posZ;
403  double deltaR = diff.perp() - posR;
404 
405  // solve the linear equation for the cylinder Radius
406  Trk::BevelledBoundaryIntersector intersectR(
407  posR, posZ, deltaZ / deltaR, m_outerRadius);
408  // the intersection point
409  double zOfIntersect = intersectR.yOfX;
410 
411  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect > 0.)
412  return {m_boundaryAccessors.bevelledcylinderAccessor(
413  Trk::BevelledCylinderZincrease)};
414  if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect < 0.)
415  return {m_boundaryAccessors.bevelledcylinderAccessor(
416  Trk::BevelledCylinderZdecrease)};
417  if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect > 0.)
418  return {m_boundaryAccessors.bevelledcylinderAccessor(
419  Trk::BevelledCylinderPositiveFace)};
420  return {m_boundaryAccessors.bevelledcylinderAccessor(
421  Trk::BevelledCylinderNegativeFace)};
422  }
423  // the sectoral cylinder case
424  if (m_innerRadius != 0. && std::abs(m_halfPhiSector - M_PI) > 10e-3)
425  return {m_boundaryAccessors.sectoralBevelledCylinderAccessor(
426  Trk::StandardSectoralBevelledCylinder)};
427  // it remains the sectoral tube case
428  return {m_boundaryAccessors.sectoralBevelledTubeAccessor(
429  Trk::StandardSectoralBevelledTube)};
430 }
431 
432 std::shared_ptr<Trk::CylinderBounds>
433 Trk::BevelledCylinderVolumeBounds::innerBevelledCylinderBounds() const
434 {
435  return std::make_shared<Trk::CylinderBounds>(m_innerRadius, m_halfPhiSector, m_halfZ);
436 }
437 
438 std::shared_ptr<Trk::CylinderBounds>
439 Trk::BevelledCylinderVolumeBounds::outerBevelledCylinderBounds() const
440 {
441  return std::make_shared<Trk::CylinderBounds>(m_outerRadius, m_halfPhiSector, m_halfZ);
442 }
443 
444 std::shared_ptr<Trk::RectangleBounds>
445 Trk::BevelledCylinderVolumeBounds::outerBevelledPlaneBounds() const
446 {
447  return std::make_shared<Trk::RectangleBounds>(m_outerRadius * tan(m_halfPhiSector), m_halfZ);
448 }
449 
450 std::shared_ptr<Trk::RectangleBounds>
451 Trk::BevelledCylinderVolumeBounds::innerBevelledPlaneBounds() const
452 {
453  return std::make_shared<Trk::RectangleBounds>(m_innerRadius * tan(m_halfPhiSector), m_halfZ);
454 }
455 
456 std::shared_ptr<Trk::EllipseBounds>
457 Trk::BevelledCylinderVolumeBounds::bottomEllipseBounds() const
458 {
459  return std::make_shared<Trk::EllipseBounds>(
460  m_innerRadius / cos(m_thetaMinus),
461  m_innerRadius,
462  m_outerRadius / cos(m_thetaMinus),
463  m_outerRadius,
464  m_halfPhiSector);
465 }
466 
467 std::shared_ptr<Trk::EllipseBounds>
468 Trk::BevelledCylinderVolumeBounds::topEllipseBounds() const
469 {
470  return std::make_shared<Trk::EllipseBounds>(
471  m_innerRadius / cos(m_thetaPlus),
472  m_innerRadius,
473  m_outerRadius / cos(m_thetaPlus),
474  m_outerRadius,
475  m_halfPhiSector);
476 }
477 
478 std::shared_ptr<Trk::CylinderBounds>
479 Trk::BevelledCylinderVolumeBounds::innerCylinderBounds() const
480 {
481  return std::make_shared<Trk::CylinderBounds>(m_innerRadius, m_halfPhiSector, m_halfZ);
482 }
483 
484 std::shared_ptr<Trk::CylinderBounds>
485 Trk::BevelledCylinderVolumeBounds::outerCylinderBounds() const
486 {
487  return std::make_shared<Trk::CylinderBounds>(m_outerRadius, m_halfPhiSector, m_halfZ);
488 }
489 
490 std::shared_ptr<Trk::DiscBounds>
491 Trk::BevelledCylinderVolumeBounds::discBounds() const
492 {
493  // adjust radius to make sure all surface covered
494  double outerRadius =
495  (m_type > 1) ? m_outerRadius / cos(m_halfPhiSector) : m_outerRadius;
496  return std::make_shared<Trk::DiscBounds>(m_innerRadius, outerRadius, m_halfPhiSector);
497 }
498 
499 std::shared_ptr<Trk::TrapezoidBounds>
500 Trk::BevelledCylinderVolumeBounds::sectorTrdBounds() const
501 {
502  return std::make_shared<Trk::TrapezoidBounds>(
503  0.5 * (outerRadius() - innerRadius()), m_halfZ, m_thetaMinus,
504  m_thetaPlus);
505 }
506 
507 std::shared_ptr<Trk::RectangleBounds>
508 Trk::BevelledCylinderVolumeBounds::sectorPlaneBounds() const
509 {
510  double ri = innerRadius();
511  double ro = outerRadius();
512  if (m_type == 1 || m_type == 3)
513  ri *= 1. / cos(halfPhiSector());
514  if (m_type > 1)
515  ro *= 1. / cos(halfPhiSector());
516  return std::make_shared<Trk::RectangleBounds>(0.5 * (ro - ri), m_halfZ);
517 }
518 
519 Trk::Volume*
520 Trk::BevelledCylinderVolumeBounds::subtractedVolume() const
521 {
522  if (m_type < 1)
523  return nullptr;
524 
525  double tp = tan(m_halfPhiSector);
526  std::unique_ptr<Trk::Volume> volIn;
527  std::unique_ptr<Trk::Volume> volOut;
528  if (m_type == 1 || m_type == 3) { // cut inner cylinder
529  volIn = std::make_unique<Trk::Volume>(
530  nullptr,
531  std::make_shared<Trk::CuboidVolumeBounds>(m_innerRadius, m_innerRadius * tp + 0.1, m_halfZ + 0.1));
532  }
533  if (m_type > 1) {
534  double hz = m_outerRadius * (1. / cos(m_halfPhiSector) - 1.);
535  volOut = std::make_unique<Trk::Volume>(
536  std::make_unique<Amg::Transform3D>(Amg::Translation3D(Amg::Vector3D(m_outerRadius + hz, 0., 0.))),
537  std::make_shared<Trk::CuboidVolumeBounds>(hz, m_outerRadius * tp + 0.1, m_halfZ + 0.1));
538  }
539 
540  if (!volIn){
541  return volOut.release();
542  }
543  else if (!volOut){
544  return volIn.release();
545  }
546  return new Trk::Volume(
547  nullptr, std::make_shared<Trk::CombinedVolumeBounds>(std::move(volIn), std::move(volOut), false));
548 }
549 
550 // ostream operator overload
551 
552 MsgStream&
553 Trk::BevelledCylinderVolumeBounds::dump(MsgStream& sl) const
554 {
555  std::stringstream sl_temp;
556  sl_temp << std::setiosflags(std::ios::fixed);
557  sl_temp << std::setprecision(7);
558  sl_temp << "Trk::BevelledCylinderVolumeBounds: (innerR, outerR, "
559  "halfPhiSector, halflengthInZ, thetaMinus, thetaPlus) = ";
560  sl_temp << "(" << m_innerRadius << ", " << m_outerRadius << ", "
561  << m_halfPhiSector << ", " << m_halfZ << ", " << m_thetaMinus << ", "
562  << m_thetaPlus << ")";
563  sl << sl_temp.str();
564  return sl;
565 }
566 
567 std::ostream&
568 Trk::BevelledCylinderVolumeBounds::dump(std::ostream& sl) const
569 {
570  std::stringstream sl_temp;
571  sl_temp << std::setiosflags(std::ios::fixed);
572  sl_temp << std::setprecision(7);
573  sl_temp << "Trk::BevelledCylinderVolumeBounds: (innerR, outerR, "
574  "halfPhiSector, halflengthInZ, thetaMinus, thetaPlus) = ";
575  sl_temp << "(" << m_innerRadius << ", " << m_outerRadius << ", "
576  << m_halfPhiSector << ", " << m_halfZ << m_thetaMinus << ", "
577  << m_thetaPlus << ")";
578  sl << sl_temp.str();
579  return sl;
580 }
581 
Trk::BevelledCylinderVolumeBounds::discBounds
std::shared_ptr< DiscBounds > discBounds() const
This method returns the associated DiscBounds for the bottom/top DiscSurface.
Definition: BevelledCylinderVolumeBounds.cxx:491
TrapezoidBounds.h
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
EllipseBounds.h
Trk::BevelledCylinderVolumeBounds::decomposeToSurfaces
virtual std::vector< std::unique_ptr< Trk::Surface > > decomposeToSurfaces(const Amg::Transform3D &transform) override
Method to decompose the Bounds into boundarySurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:100
Trk::BevelledTubeRincreaseZincrease
@ BevelledTubeRincreaseZincrease
Accessor type [ 2,1,0,3 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:29
Trk::BevelledCylinderVolumeBounds::m_thetaMinus
double m_thetaMinus
Definition: BevelledCylinderVolumeBounds.h:199
DiscBounds.h
Trk::BevelledCylinderVolumeBounds
Definition: BevelledCylinderVolumeBounds.h:100
RectangleBounds.h
Trk::BevelledTubeRdecreaseZincrease
@ BevelledTubeRdecreaseZincrease
Accessor type [ 3,1,0,2 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:33
Trk::BevelledCylinderPositiveFace
@ BevelledCylinderPositiveFace
Pos.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:21
Trk::BevelledCylinderNegativeFace
@ BevelledCylinderNegativeFace
Neg.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:22
Trk::BevelledCylinderVolumeBounds::s_numericalStable
static const double s_numericalStable
numerical stability
Definition: BevelledCylinderVolumeBounds.h:212
Trk::BevelledCylinderVolumeBounds::m_type
int m_type
Definition: BevelledCylinderVolumeBounds.h:201
Trk::BevelledCylinderVolumeBounds::m_outerRadius
double m_outerRadius
Definition: BevelledCylinderVolumeBounds.h:196
M_PI
#define M_PI
Definition: ActiveFraction.h:11
deg
#define deg
Definition: SbPolyhedron.cxx:17
mc.diff
diff
Definition: mc.SFGenPy8_MuMu_DD.py:14
Trk::BevelledCylinderVolumeBounds::BevelledCylinderVolumeBounds
BevelledCylinderVolumeBounds()
Default Constructor.
Definition: BevelledCylinderVolumeBounds.cxx:34
ParticleTest.tp
tp
Definition: ParticleTest.py:25
Trk::BevelledTubeZincreaseRincrease
@ BevelledTubeZincreaseRincrease
Accessor type [ 1,2,0,3 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:31
Trk::BevelledTubeOutsideRminRdecrease
@ BevelledTubeOutsideRminRdecrease
Accessor type [ 3,1,0,2] - inverse case.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:37
Trk::BevelledCylinderVolumeBounds::sectorPlaneBounds
std::shared_ptr< RectangleBounds > sectorPlaneBounds() const
Definition: BevelledCylinderVolumeBounds.cxx:508
Trk::BevelledTubeOutsideRmaxRincrease
@ BevelledTubeOutsideRmaxRincrease
Accessor type [ 2,1,0,3 ] - inverse case.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:38
Trk::DiscSurface
Definition: DiscSurface.h:54
Trk::BevelledTubeOutsideZminZdecrease
@ BevelledTubeOutsideZminZdecrease
Accessor type [ 0,3,2,1 ] - inverse case.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:39
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
Trk::BevelledCylinderVolumeBounds::m_innerRadius
double m_innerRadius
Definition: BevelledCylinderVolumeBounds.h:195
Trk::BevelledCylinderVolumeBounds::outerBevelledCylinderBounds
std::shared_ptr< CylinderBounds > outerBevelledCylinderBounds() const
This method returns the associated BevelledCylinderBounds of the outer BevelledCylinderSurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:439
Trk::BevelledCylinderVolumeBounds::outerCylinderBounds
std::shared_ptr< CylinderBounds > outerCylinderBounds() const
This method returns the associated CylinderBounds of the outer CylinderSurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:485
python.SystemOfUnits.millimeter
float millimeter
Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration # This software is distribute...
Definition: SystemOfUnits.py:62
Trk::BevelledCylinderVolumeBounds::bottomEllipseBounds
std::shared_ptr< EllipseBounds > bottomEllipseBounds() const
This method returns the associated EllipseBounds for the bottom/top EllipseSurface.
Definition: BevelledCylinderVolumeBounds.cxx:457
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
Trk::BevelledTubeZdecreaseRincrease
@ BevelledTubeZdecreaseRincrease
Accessor type [ 0,2,1,3 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:32
Trk::BevelledCylinderVolumeBounds::m_halfZ
double m_halfZ
Definition: BevelledCylinderVolumeBounds.h:198
Trk::StandardSectoralBevelledCylinder
@ StandardSectoralBevelledCylinder
Definition: BevelledCylinderVolumeBoundaryAccessors.h:48
Trk::BevelledCylinderVolumeBounds::m_thetaPlus
double m_thetaPlus
Definition: BevelledCylinderVolumeBounds.h:200
m_type
TokenType m_type
the type
Definition: TProperty.cxx:44
VolumeExcluder.h
Trk::VolumeBounds
Definition: VolumeBounds.h:45
Trk::BevelledTubeRdecreaseZdecrease
@ BevelledTubeRdecreaseZdecrease
Accessor type [ 3,0,1,2 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:34
Trk::BevelledCylinderVolumeBounds::dump
MsgStream & dump(MsgStream &sl) const override
Output Method for MsgStream.
Definition: BevelledCylinderVolumeBounds.cxx:553
Trk::BevelledCylinderVolumeBounds::boundarySurfaceAccessor
ObjectAccessor boundarySurfaceAccessor(const Amg::Vector3D &gp, const Amg::Vector3D &dir, bool forceInside=false) const override
Provide accessor for BoundarySurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:249
beamspotman.posZ
posZ
Definition: beamspotman.py:1622
BevelledCylinderVolumeBounds.h
Amg::Transform3D
Eigen::Affine3d Transform3D
Definition: GeoPrimitives.h:46
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
CylinderSurface.h
Trk::BevelledTubeRincreaseZdecrease
@ BevelledTubeRincreaseZdecrease
Accessor type [ 2,0,1,3 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:30
Trk::BevelledCylinderVolumeBounds::operator=
BevelledCylinderVolumeBounds & operator=(const BevelledCylinderVolumeBounds &cylbo)
Assignment operator.
Definition: BevelledCylinderVolumeBounds.cxx:82
drawFromPickle.tan
tan
Definition: drawFromPickle.py:36
Trk::BevelledTubeZdecreaseRdecrease
@ BevelledTubeZdecreaseRdecrease
Accessor type [ 0,3,1,2 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:36
Trk::BevelledTubeZincreaseRdecrease
@ BevelledTubeZincreaseRdecrease
Accessor type [ 1,3,0,2 ].
Definition: BevelledCylinderVolumeBoundaryAccessors.h:35
beamspotman.dir
string dir
Definition: beamspotman.py:621
Trk::BevelledCylinderVolumeBounds::topEllipseBounds
std::shared_ptr< EllipseBounds > topEllipseBounds() const
This method returns the associated EllipseBounds for the bottom/top EllipseSurface.
Definition: BevelledCylinderVolumeBounds.cxx:468
Trk::BevelledBoundaryIntersector::yOfX
double yOfX
the result of x
Definition: BevelledCylinderVolumeBounds.h:29
Trk::BevelledCylinderZincrease
@ BevelledCylinderZincrease
Cylinder hit, then pos.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:19
Trk::BevelledCylinderVolumeBounds::innerBevelledPlaneBounds
std::shared_ptr< RectangleBounds > innerBevelledPlaneBounds() const
This method returns the associated plane surface bounds of the inner bevelled surface.
Definition: BevelledCylinderVolumeBounds.cxx:451
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
Trk::ObjectAccessor
Definition: ObjectAccessor.h:15
Trk::BevelledCylinderVolumeBounds::sectorTrdBounds
std::shared_ptr< TrapezoidBounds > sectorTrdBounds() const
This method returns the associated PlaneBounds limiting a sectoral BevelledCylinderVolume.
Definition: BevelledCylinderVolumeBounds.cxx:500
Trk::BevelledBoundaryIntersector
Definition: BevelledCylinderVolumeBounds.h:27
Amg::RotationMatrix3D
Eigen::Matrix< double, 3, 3 > RotationMatrix3D
Definition: GeoPrimitives.h:49
PlaneSurface.h
Trk::BevelledCylinderVolumeBounds::innerCylinderBounds
std::shared_ptr< CylinderBounds > innerCylinderBounds() const
This method returns the associated CylinderBounds of the inner CylinderSurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:479
Trk::BevelledCylinderVolumeBounds::innerBevelledCylinderBounds
std::shared_ptr< CylinderBounds > innerBevelledCylinderBounds() const
This method returns the associated BevelledCylinderBounds of the inner BevelledCylinderSurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:433
Amg::Translation3D
Eigen::Translation< double, 3 > Translation3D
Definition: GeoPrimitives.h:44
Trk::BevelledTubeOutsideZmaxZincrease
@ BevelledTubeOutsideZmaxZincrease
Accessor type [ 1,3,2,0 ] - inverse case.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:40
CombinedVolumeBounds.h
DiscSurface.h
Trk::BevelledCylinderVolumeBounds::~BevelledCylinderVolumeBounds
virtual ~BevelledCylinderVolumeBounds()
Destructor.
Trk::BevelledCylinderVolumeBounds::subtractedVolume
Volume * subtractedVolume() const
This method returns the bevelled area volume.
Definition: BevelledCylinderVolumeBounds.cxx:520
Amg::AngleAxis3D
Eigen::AngleAxisd AngleAxis3D
Definition: GeoPrimitives.h:45
SubtractedDiscSurface.h
fitman.hz
def hz
Definition: fitman.py:516
makeComparison.deltaZ
int deltaZ
Definition: makeComparison.py:46
Trk::BevelledCylinderZdecrease
@ BevelledCylinderZdecrease
Cylinder hit, the neg.
Definition: BevelledCylinderVolumeBoundaryAccessors.h:20
Trk::Volume
Definition: Volume.h:36
Trk::BevelledCylinderVolumeBounds::outerBevelledPlaneBounds
std::shared_ptr< RectangleBounds > outerBevelledPlaneBounds() const
This method returns the associated BevelledCylinderBounds of the outer BevelledCylinderSurfaces.
Definition: BevelledCylinderVolumeBounds.cxx:445
makeComparison.deltaR
float deltaR
Definition: makeComparison.py:36
Trk::StandardSectoralBevelledTube
@ StandardSectoralBevelledTube
Definition: BevelledCylinderVolumeBoundaryAccessors.h:55
CuboidVolumeBounds.h
Trk::BevelledCylinderVolumeBounds::m_halfPhiSector
double m_halfPhiSector
Definition: BevelledCylinderVolumeBounds.h:197
Generated on Mon Sep 1 2025 21:06:58 for ATLAS Offline Software by   doxygen 1.8.18