10–12 Sept 2014
University of Pisa
Europe/Rome timezone

Sampling secondary particles in high energy physics simulation on the GPU

12 Sept 2014, 09:00
30m
University of Pisa

University of Pisa

<a target="_blank" href=https://www.google.com/maps/place/Dipartimento+di+Fisica/@43.720239,10.407985,17z/data=!3m1!4b1!4m2!3m1!1s0x12d591bb7d8c8ec9:0xbf91ddd442e32978>Polo Fibonacci</a> Largo Bruno Pontecorvo, 3 I-56127 Pisa <em>phone +39 050 2214 327</em>

Speaker

Soon Yung Jun (Fermilab)

Description

We present a massively parallel application for sampling secondary particles in high energy physics (HEP) simulation on a Graphics Processing Unit (GPU). HEP experiments primarily uses the Geant4 toolkit (Geant4) to simulate the passage of particles through a general-purpose detector, which requires intensive computing resources due to the complexity of the geometry as well as physics processes applied to particles copiously produced by primary collisions and secondary interactions. The combined composition and rejection methods to sample secondary particles often used in Geant4 may not be suitable for optimal performance of HEP events simulation using recent hardware architectures of accelerated or many-core processors owing to the stochastic nature of the Monte Carlo technique. An alternative approach based on a discrete inverse cumulative probability distribution is explored to minimize the divergence in thread level parallelism as well as to vectorize physics processes for spatial locality and instruction throughput. The inverse cumulative distribution of the differential cross section associated with each electromagnetic physics process is tabulated based on algorithms excerpted from Geant4 and a simple random sampling technique with a linear interpolation is implemented for GPU. Validation and performance evaluation with the alternative technique compared to the conventional composition and rejection method both on GPU and CPU are presented.

Primary author

Soon Yung Jun (Fermilab)

Co-authors

Andrei GHEATA (CERN) Federico CARMINATI (CERN) Georgios BITZES (University of Athens (GR)) Guilherme LIMA (Fermi National Accelerator Lab. (US)) Johannes Christof DE FINE LICHT (University of Copenhagen (DK)) John APOSTOLAKIS (CERN) Laurent DUHEM (INTEL) Marilena BANDIERAMONTE (Universita e INFN (IT)) Mihaly NOVAK (CERN) Oksana SHADURA (National Technical Univ. of Ukraine "Kyiv Polytechnic Institute) Philippe CANAL (Fermi National Accelerator Lab. (US)) Raman SEHGAL (Bhabha Atomic Research Centre (IN)) Rene BRUN (CERN) Sandro Christian WENZEL (CERN) Victor Daniel ELVIRA (Fermi National Accelerator Lab. (US))

Presentation materials