I’m working on optimizing geometric navigation in radiation therapy simulations using VecGeom with CT-based geometries. CT images contain regions of varying tissue densities (e.g., bone ≈1.8 g/cm³, soft tissue ≈1.0 g/cm³, lungs ≈0.3 g/cm³). To balance efficiency and accuracy, I want to dynamically adjust BVH parameters: finer BVH structures for high-density regions (e.g., bone) and coarser structures for low-density regions (e.g., lungs).
However, I’m unsure how to implement this density-aware BVH parameter tuning in practice. VecGeom’s documentation doesn’t explicitly mention a strategy class like vecgeom::navigation::BVHParameterStrategy.
My questions are:
- Which existing classes or APIs in VecGeom control BVH parameters (e.g., max leaf size, split strategy) that I can extend or modify?
- How can I integrate density-based logic to adjust these parameters—should this be done within VecGeom’s source code, or via the Geant4-VecGeom interface layer?
- Are there code examples or best practices for dynamic BVH parameter tuning based on material properties in VecGeom?