19#include "BDSExtent.hh"
20#include "BDSGeometryInspector.hh"
21#include "BDSGlobalConstants.hh"
24#include "G4AffineTransform.hh"
26#include "G4CutTubs.hh"
27#include "G4DisplacedSolid.hh"
28#include "G4EllipticalTube.hh"
29#include "G4IntersectionSolid.hh"
30#include "G4SubtractionSolid.hh"
32#include "G4UnionSolid.hh"
33#include "G4Version.hh"
39#if G4VERSION_NUMBER < 1030
40#include "G4VisExtent.hh"
47 G4String className = solid->GetEntityType();
48 if (className ==
"G4DisplacedSolid")
50 else if (className ==
"G4SubtractionSolid")
52 else if (className ==
"G4IntersectionSolid")
54 else if (className ==
"G4UnionSolid")
56 else if (className ==
"G4Box")
58 else if (className ==
"G4Tubs")
60 else if (className ==
"G4CutTubs")
62 else if (className ==
"G4EllipticalTube")
66 G4cout <<
"BDS::DetermineExtents> Solid of type \"" << className <<
"\" is not supported" << G4endl;
73 const G4DisplacedSolid* solid =
dynamic_cast<const G4DisplacedSolid*
>(solidIn);
77 const G4VSolid* constituentSolid = solid->GetConstituentMovedSolid();
80 G4AffineTransform transform = solid->GetTransform();
83 G4ThreeVector translation = transform.NetTranslation();
87 return std::make_pair(outer, inner);
92 const G4SubtractionSolid* solid =
dynamic_cast<const G4SubtractionSolid*
>(solidIn);
96 const G4VSolid* a = solid->GetConstituentSolid(0);
97 const G4VSolid* b = solid->GetConstituentSolid(1);
108 return std::make_pair(outer, inner);
113 const G4IntersectionSolid* solid =
dynamic_cast<const G4IntersectionSolid*
>(solidIn);
117 const G4VSolid* a = solid->GetConstituentSolid(0);
118 const G4VSolid* b = solid->GetConstituentSolid(1);
126 return std::make_pair(outer, inner);
131 const G4UnionSolid* solid =
dynamic_cast<const G4UnionSolid*
>(solidIn);
135 const G4VSolid* a = solid->GetConstituentSolid(0);
136 const G4VSolid* b = solid->GetConstituentSolid(1);
144 return std::make_pair(outer, inner);
149 const G4Box* solid =
dynamic_cast<const G4Box*
>(solidIn);
153 G4double dx = solid->GetXHalfLength();
154 G4double dy = solid->GetYHalfLength();
155 G4double dz = solid->GetZHalfLength();
161 return std::make_pair(outer, inner);
166 const G4Tubs* solid =
dynamic_cast<const G4Tubs*
>(solidIn);
170 G4double innerR = solid->GetInnerRadius();
171 G4double outerR = solid->GetOuterRadius();
172 G4double zHalfL = solid->GetZHalfLength();
175 BDSExtent outer(outerR + lsl, outerR + lsl, zHalfL);
177 return std::make_pair(outer, inner);
182 G4ThreeVector zminV = G4ThreeVector();
183 G4ThreeVector zmaxV = G4ThreeVector();
184#if G4VERSION_NUMBER > 1039
185 solidIn->BoundingLimits(zminV, zmaxV);
186#elif G4VERSION_NUMBER > 1029
187 solidIn->Extent(zminV, zmaxV);
189 G4VisExtent v = solidIn->GetExtent();
190 zmaxV = G4ThreeVector(std::max(std::abs(v.GetXmin()),std::abs(v.GetXmax())),
191 std::max(std::abs(v.GetYmin()),std::abs(v.GetYmax())),
192 std::max(std::abs(v.GetZmin()),std::abs(v.GetZmax())));
195 const G4CutTubs* solid =
dynamic_cast<const G4CutTubs*
>(solidIn);
199 G4double innerR = solid->GetInnerRadius();
200 G4double outerR = solid->GetOuterRadius();
206 -zmaxV.z(), zmaxV.z());
209 -zmaxV.z(), zmaxV.z());
210 return std::make_pair(outer, inner);
215 const G4EllipticalTube* solid =
dynamic_cast<const G4EllipticalTube*
>(solidIn);
219 G4double dX = solid->GetDx();
220 G4double dY = solid->GetDy();
221 G4double dZ = solid->GetDz();
226 return std::make_pair(outer, inner);
Holder for +- extents in 3 dimensions.
BDSExtent Translate(const G4ThreeVector &offset) const
Provide a new copy of this extent with an offset applied.
G4bool TransverselyLessThan(const BDSExtent &r) const
Comparison operator for x,y only. Ignores z (length).
static BDSGlobalConstants * Instance()
Access method.
std::pair< BDSExtent, BDSExtent > InspectIntersectionSolid(const G4VSolid *solidIn)
Inspect a G4IntersectionSolid.
std::pair< BDSExtent, BDSExtent > InspectEllipticalTube(const G4VSolid *solidIn)
Inspect a G4EllipticalTube.
std::pair< BDSExtent, BDSExtent > InspectCutTubs(const G4VSolid *solidIn)
Inspect a G4CutTubs.
std::pair< BDSExtent, BDSExtent > InspectSubtractionSolid(const G4VSolid *solidIn)
std::pair< BDSExtent, BDSExtent > DetermineExtents(const G4VSolid *solid)
std::pair< BDSExtent, BDSExtent > InspectBox(const G4VSolid *solidIn)
Inspect a G4Box.
std::pair< BDSExtent, BDSExtent > InspectTubs(const G4VSolid *solidIn)
Inspect a G4Tubs.
std::pair< BDSExtent, BDSExtent > InspectDisplacedSolid(const G4VSolid *solidIn)
Inspect a G4DisplacedSolid.
std::pair< BDSExtent, BDSExtent > InspectUnionSolid(const G4VSolid *solidIn)
Inspect a G4UnionSolid.