Proton Therapy Treatment Simulations with the Bergen DTC Prototype for Range Verification

<p>This dataset contains&nbsp;proton therapy Monte Carlo simulations using GATE (Jan et al. 2004)&nbsp;version 9.2 and Geant4&nbsp;(Agostinelli et al. 2003; Allison et al. 2006;&nbsp;Allison et al. 2016)&nbsp;version 11.0.0. The treatment targets are&nbsp;the&nbsp;pediatric head phantom 715-HN&nbsp;by CIRS Inc.&nbsp;(Norfolk, VA, United States), digitized&nbsp;by Giacometti et al. (2017), and The Visible Human Female (VHF) head, courtesy of the U.S. National Library of Medicine, resampled to 1 mm voxels and scaled down&nbsp;to 80%. Distal to the phantom, is a simplified geometry of the digital tracking calorimeter designed by the&nbsp;<a href="https://www.uib.no/en/ift/142356/medical-physics-bergen-pct-project">Bergen pCT collaboration</a>&nbsp;for proton computed tomography (Alme et al. 2020).</p> <ul> <li><strong>treatment_simulations.tar.gz:</strong>&nbsp;36258 simulations with the 715-HN phantom</li> <li><strong>treatment_simulations_vhf.tar.gz:</strong>&nbsp;35673 simulations with the VHF phantom</li> </ul> <p>Each simulation is a single pencil beam&nbsp;with <span class="math-tex">\(10^7\)</span>&nbsp;primary protons. To determine valid spots to use,&nbsp;probing simulations with&nbsp;<span class="math-tex">\(10^5\)</span>&nbsp;primaries were conducted, where the beam position, energy, and phantom rotation were varied:&nbsp;10 mm interval in the xy-plane for the beam position,&nbsp;3 mm water range interval for the beam energy, and 30&deg; rotations of the phantom. If no primary was found in the detecter after the probing simulation, it is a valid spot for the dataset.</p> <p>Resulting hit files, containing simulated detector readout, can be found in the archive files for the respective phantom. Each simulation is accompanied by a metafile in JSON format, specifying all parameters used to produce the output. E.g., beam position and energy can be found in &quot;parameters&quot; -&gt; &quot;beam_spot_x&quot;, &quot;beam_spot_y&quot;, &quot;beam_energy&quot; and the phantom rotation can be found in &quot;parameters&quot; -&gt; &quot;phantom_rotation_angle&quot;.</p> <p>Simulation data has been pre-processed with two steps. First, hits for the same track in the same layer are binned together into a single entry by averaging the positions (posX, posY, posZ)&nbsp;and summing the energy depositions (edep). Second, the energy deposition values are discretized into a cluster size and back into MeV, to represent the energy resolution of the ALPIDE chip (Alme et al. 2020). Clusters of size 0 are removed.&nbsp;Additionally, some of the unused GATE output columns are removed and a column specifying the&nbsp;detector layer is created through the formula layer = 2*volumeID[2] + volumeID[3]. Details about the&nbsp;meaning of the other output columns can be found in the <a href="https://opengate.readthedocs.io/en/latest/data_output_management.html">GATE documentation</a>.</p>