TURTLE: A C library for an optimistic stepping through a topography

Abstract TURTLE is a C library providing utilities allowing to navigate through a topography described by a Digital Elevation Model (DEM). The library has been primarily designed for the Monte Carlo transport of particles scattering over medium to long ranges, e.g. atmospheric muons. But, it can also efficiently handle ray tracing problems with very large DEMs ( 1 0 9 nodes or more), e.g. for neutrino simulations. The TURTLE library was built on an optimistic ray tracing algorithm, detailed in the present paper. This algorithm proceeds by trials and errors, approximating the topography within the modelling uncertainties of the DEM data. This allows to traverse a topography in constant time, i.e. independently of the number of grid nodes, and with no added memory. Detailed performance studies are provided by comparison to other ray tracing algorithms and as an application to muon transport in a Monte Carlo simulation. Program summary Program Title: The TURTLE library Program Files doi: http://dx.doi.org/10.17632/pbnsctnbk8.1 Licensing provisions: LGPL-3.0 Programming language: C99 Nature of problem: Transport long range Monte Carlo particles through a topography described by a high resolution Digital Elevation Model (DEM). Solution method: Mapping the ground surface on the fly using a heuristic requiring only the original DEM data.

[1]  A. Dell'Acqua,et al.  Geant4 - A simulation toolkit , 2003 .

[2]  R. D. Garg,et al.  Evaluation of vertical accuracy of open source Digital Elevation Model (DEM) , 2013, Int. J. Appl. Earth Obs. Geoinformation.

[3]  Josef Weidendorfer,et al.  A Tool Suite for Simulation Based Analysis of Memory Access Behavior , 2004, International Conference on Computational Science.

[4]  Valentin Niess,et al.  Backward Monte-Carlo applied to muon transport , 2017, Comput. Phys. Commun..

[5]  D.K. Olson,et al.  Converting Earth-centered, Earth-fixed coordinates to geodetic coordinates , 1996, IEEE Transactions on Aerospace and Electronic Systems.

[6]  D. Caruge,et al.  Recent Developments in Geant4 , 2014, ICS 2014.

[7]  G. Cosmo,et al.  The Geant4 geometry modeler , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[8]  Y. Zheng,et al.  First detection of extensive air showers by the TREND self-triggering radio experiment , 2010, 1007.4359.

[9]  Christopher M Poole,et al.  Fast Tessellated Solid Navigation in GEANT4 , 2012, IEEE Transactions on Nuclear Science.

[10]  Melanie Johnston-Hollitt,et al.  RADIO SOURCES IN THE NCP REGION OBSERVED WITH THE 21 CENTIMETER ARRAY , 2016, 1602.06624.

[11]  Nicholas Nethercote,et al.  Valgrind: a framework for heavyweight dynamic binary instrumentation , 2007, PLDI '07.

[12]  R. T. Williams.  Lambert and Mercator map projections in geology and geophysics , 1995 .

[13]  John Apostolakis,et al.  A vectorization approach for multifaceted solids in VecGeom , 2019, EPJ Web of Conferences.

[14]  John S. Hendricks,et al.  MCNPX 2.7.0 extensions , 2011 .

[15]  Shanghai,et al.  Giant Radio Array for Neutrino Detection: Science and Design , 2018 .

[16]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[17]  S. Incerti,et al.  Geant4 developments and applications , 2006, IEEE Transactions on Nuclear Science.

[18]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[19]  Pierre Alliez,et al.  Computational geometry algorithms library , 2008, SIGGRAPH '08.

[20]  David E. Smith,et al.  The Development of the NASA GSFC and NIMA Joint Geopotential Model , 1997 .

[21]  F. Garufi,et al.  Joint measurement of the atmospheric muon flux through the Puy de Dôme volcano with plastic scintillators and Resistive Plate Chambers detectors , 2015 .