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
|
/// \file BVHNavigator.h
/// \author Guilherme Amadio
#ifndef VECGEOM_NAVIGATION_BVHNAVIGATOR_H_
#define VECGEOM_NAVIGATION_BVHNAVIGATOR_H_
#include "VecGeom/management/BVHManager.h"
#include "VecGeom/navigation/BVHSafetyEstimator.h"
#include "VecGeom/navigation/VNavigator.h"
#include "VecGeom/volumes/LogicalVolume.h"
#include "VecGeom/volumes/PlacedVolume.h"
namespace vecgeom {
inline namespace VECGEOM_IMPL_NAMESPACE {
/**
* @brief Navigator class using the bounding volume hierarchy of each logical volume for acceleration.
*/
template <bool MotherIsConvex = false>
class BVHNavigator : public VNavigatorHelper<BVHNavigator<MotherIsConvex>, MotherIsConvex> {
private:
/** Constructor. Private since this is a singleton class accessed only via the @c Instance() static method. */
VECCORE_ATT_DEVICE
BVHNavigator()
: VNavigatorHelper<BVHNavigator<MotherIsConvex>, MotherIsConvex>(BVHSafetyEstimator::Instance())
{
}
public:
using SafetyEstimator_t = BVHSafetyEstimator;
using Base = VNavigatorHelper<BVHNavigator<MotherIsConvex>, MotherIsConvex>;
using Base::CheckDaughterIntersections;
static constexpr const char *gClassNameString = "BVHNavigator";
#ifndef VECCORE_CUDA
/** Returns the instance of this singleton class. */
static VNavigator *Instance()
{
static BVHNavigator instance;
return &instance;
}
#else
// If used on device, this needs to be implemented in a .cu file rather than in this header
// This hack is used also by NewSimpleNavigator, implemented in LogicalVolume.cpp
// This is now implemented in BVHManager.cu
VECCORE_ATT_DEVICE
static VNavigator *Instance();
#endif
/**
* Checks for intersections against child volumes of logical volume @p lvol, using the BVH
* associated with it.
* @param[in] lvol Logical volume being checked.
* @param[in] localpoint Point in the local coordinates of the logical volume.
* @param[in] localdir Direction in the local coordinates of the logical volume.
* @param[in] in_state Incoming navigation state.
* @param[in] out_state Outgoing navigation state (not used by this method).
* @param[in] step Maximum step size. Volumes beyond this distance are ignored.
* @param[out] hitcandidate
* @returns Whether @p out_state has been modified or not. Always false for this method.
*/
VECCORE_ATT_HOST_DEVICE
bool CheckDaughterIntersections(LogicalVolume const *lvol, Vector3D<Precision> const &localpoint,
Vector3D<Precision> const &localdir, NavigationState const *in_state,
NavigationState * /* out_state */, Precision &step,
VPlacedVolume const *&hitcandidate) const final
{
if (auto bvh = BVHManager::GetBVH(lvol)) {
VPlacedVolume const *last = in_state ? in_state->GetLastExited() : nullptr;
bvh->CheckDaughterIntersections(localpoint, localdir, step, last, hitcandidate);
}
return false; /* return value indicates whether out_state has been modified */
}
};
} // namespace VECGEOM_IMPL_NAMESPACE
} // namespace vecgeom
#endif
|
