1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
/*
* PlacedPolycone.h
*
* Created on: May 14, 2014
* Author: swenzel
*/
#ifndef VECGEOM_VOLUMES_PLACEDPOLYCONE_H_
#define VECGEOM_VOLUMES_PLACEDPOLYCONE_H_
#include "VecGeom/base/Global.h"
#ifndef VECCORE_CUDA
#include "VecGeom/base/RNG.h"
#include <cmath>
#endif
#include "VecGeom/volumes/PlacedVolume.h"
#include "VecGeom/volumes/UnplacedVolume.h"
#include "VecGeom/volumes/kernel/PolyconeImplementation.h"
#include "VecGeom/volumes/PlacedVolImplHelper.h"
#include "VecGeom/volumes/UnplacedPolycone.h"
namespace vecgeom {
VECGEOM_DEVICE_FORWARD_DECLARE(class PlacedPolycone;);
VECGEOM_DEVICE_DECLARE_CONV(class, PlacedPolycone);
inline namespace VECGEOM_IMPL_NAMESPACE {
class PlacedPolycone : public VPlacedVolume {
public:
#ifndef VECCORE_CUDA
PlacedPolycone(char const *const label, LogicalVolume const *const logicalVolume,
Transformation3D const *const transformation)
: VPlacedVolume(label, logicalVolume, transformation)
{
}
PlacedPolycone(LogicalVolume const *const logicalVolume, Transformation3D const *const transformation)
: PlacedPolycone("", logicalVolume, transformation)
{
}
#else
VECCORE_ATT_DEVICE PlacedPolycone(LogicalVolume const *const logicalVolume,
Transformation3D const *const transformation, const int id, const int copy_no,
const int child_id)
: VPlacedVolume(logicalVolume, transformation, id, copy_no, child_id)
{
}
#endif
VECCORE_ATT_HOST_DEVICE
virtual ~PlacedPolycone() {}
VECCORE_ATT_HOST_DEVICE
UnplacedPolycone const *GetUnplacedVolume() const
{
return static_cast<UnplacedPolycone const *>(GetLogicalVolume()->GetUnplacedVolume());
}
VECCORE_ATT_HOST_DEVICE
virtual void PrintType() const override;
virtual void PrintType(std::ostream &s) const override;
bool IsOpen() const { return (GetUnplacedVolume()->GetDeltaPhi() < kTwoPi); }
Precision GetStartPhi() const { return GetUnplacedVolume()->GetStartPhi(); }
Precision GetEndPhi() const { return GetUnplacedVolume()->GetEndPhi(); }
int GetNumRZCorner() const { return 2 * (int)(GetUnplacedVolume()->GetNz()); } // in USolids nCorners = 2*nPlanes
PolyconeHistorical *GetOriginalParameters() const { return GetUnplacedVolume()->GetOriginalParameters(); }
void Reset() { const_cast<UnplacedPolycone *>(GetUnplacedVolume())->Reset(); }
#ifndef VECCORE_CUDA
virtual VPlacedVolume const *ConvertToUnspecialized() const override;
#ifdef VECGEOM_ROOT
virtual TGeoShape const *ConvertToRoot() const override;
#endif
#ifdef VECGEOM_GEANT4
virtual G4VSolid const *ConvertToGeant4() const override;
#endif
#endif // VECGEOM_TEST_BENCHMARK
#if !defined(VECCORE_CUDA)
virtual Precision SurfaceArea() const override { return GetUnplacedVolume()->SurfaceArea(); }
#endif
}; // end class
template <typename UnplacedPolycone_t>
class SPlacedPolycone : public PlacedVolumeImplHelper<UnplacedPolycone_t, PlacedPolycone> {
using Base = PlacedVolumeImplHelper<UnplacedPolycone_t, PlacedPolycone>;
public:
typedef UnplacedPolycone UnplacedShape_t;
using Base::Base;
};
} // namespace VECGEOM_IMPL_NAMESPACE
} // namespace vecgeom
#endif // VECGEOM_VOLUMES_PLACEDPOLYCONE_H_
|
