Dune cdr the Single-Phase Protodune

Chapter 5: Software and Computing
5–132
provided [28] for supported versions of Linux and Mac OS X. Tools for compiling the framework
and applications are also provided, along with all of the required dependencies, including the gnu
C++ compiler. The versions of the software and its dependencies are managed by the UNIX
Product Support (UPS) system, which allows easy version selection, setup and configuration of
the LArSoft environment on computers with previous versions already installed, such as those at
Fermilab.
LArSoft is under rapid development by both the core LArSoft team and by contributors from
participating experiments: ArgoNeuT, LArIAT, MicroBooNE, SBND, and DUNE. Within DUNE,
the LArTPC near-detector option, both ProtoDUNE detectors, and the Far Detector are clients
and contributors of LArSoft.
5.6
Event simulation
Three main sources of particles are simulated in protoDUNE-SP: beam particles coming through
the beamline, and beam halo particles from both our own line (H4) and from upstream sources
and cosmic ray particles. Each particle type has its dedicated generating process and the mixing
of the right relative and absolute amounts in the relevant time frames is of crucial importance.
Beam particles are generated using the G4Beamline event generator which is a standard tool used
at CERN for simulation of beamlines.
The upstream beam instrumentation devices – wire chambers, Cherenkov counters, and time-of-
flight-counters – are fully simulated as described in Section 6.3. Cosmic-ray events are simulated
either with the CRY [29] event generator or CORSIKA [30]. Neutrino scattering events are sim-
ulated using GENIE [31]; while these events will be very rare in the detector, the extrapolation
of the performance of the ProtoDUNE-SP detector to the FD will require simulating them. A
dedicated generator in LArSoft simulates radionuclide decay products which can be overlaid on
other events.
The detector geometry is coded in GDML files [27] that are generated by the GeGeDe [32] geometry
system. These files contain the locations, sizes, shapes, and material content of the detector
components, the active liquid argon volume, and the surrounding materials, such as the field
cage, the beam windows, the cryostat the supporting structure, and the experimental hall. These
external features will impact the distributions of cosmic-ray particles impinging on the active
detector. The channels and volumes are numbered and named in the GDML files, with conventions
followed by the LArSoft simulation code.
The active volume of the detector is divided into cubes 300 µm on a side, called voxels. GEANT4
tracks particles through the argon with each step ending on a voxel boundary. This allows the
simulation of small-scale physics processes, such as delta-ray emission and showering, at a level of
detail smaller than the intrinsic resolution of the detector. While GEANT4 calculates the energy
deposited by each particle for each step, the simulation of ionization and scintillation-photon
emission is performed using one of two algorithms in LArSoft: a dedicated parameterization that
depends on the electric field in the liquid argon and the ionization density [33], or NEST [34], which
ProtoDUNE Single-Phase Technical Design Report

---

Chapter 5: Software and Computing
5–133
is tuned to previous noble-liquid experimental results and introduces an anti-correlation between
the photon yield and the ionization electron yield for each step.
The availability of alternate interaction models in MC simulation may be helpful for estimating
systematic uncertainty. FLUKA [35, 36, 37] is a MC package with detailed modeling of particles
traversing the detector. The option of creating an interface to LArSoft will be considered if
resources are available.
The average specific energy loss for a minimum-ionizing particle (MIP) is approximately 2.12 MeV/cm.
The W -value for ionization is 23.6 eV per electron-ion pair, and the W -value for scintillation is
19.5 eV per photon, resulting in tens of thousands of drifting electrons and photons per cm of
charged-particle track in the detector. It is impractical to simulate the paths of each of these
electrons and photons using GEANT4, and therefore computational techniques are incorporated
into LArSoft to achieve a high simulation speed while preserving accuracy. The electrons are prop-
agated by LArSoft-specific tools, including a tool that integrates over the distributions created by
longitudinal and transverse diffusion, and a geometry tool that indicates the wire locations that
record charge, assuming uniform wire spacing.
The effect of charge loss due to attachment of electrons to impurities (i.e., the effect of the electron
lifetime), is implemented in this step. Effects due to space charge are simulated using a smoothly
parameterized map of distortions in the apparent position components (x, y, z) of the charge de-
posits as functions of the true position. This map can be made using SPaCE [38], a program that
traces particle trajectories in LAr based on the electric field calculated using Poisson’s equation, a
given space-charge density map, and the boundary conditions provided by the cathodes, anodes,
and field cages. It is anticipated that the space-charge distortions in ProtoDUNE-SP may be as
large as 20 cm [38].
Photon propagation is simulated using a library that contains the probabilities of observing a
photon emitted at a particular point in space by a particular photon detector. Here, the space is
divided into cubical voxels 6 cm on a side, and the library is indexed by photon detector element.
The simulated arrival times and charge amounts on each wire, and of each set of photons arriving
at each photon detector, along with the identity of the particles generating them, are stored in the
simulation output file for use in determining the performance of the downstream reconstruction
algorithms. These charge depositions and photons are inputs to the detector response functions –
the field response and the electronics response are convoluted with the true arrival times to make
simulated waveforms. The detector field response functions are simulated using GARFIELD [39],
but they will be validated with real data, as the simulation contains oversimplifications, such as
inadequate modeling of induction signals.
The electronics gain is applied so that the simulated signals match the expected responses. Simu-
lated noise is then added, and the result is quantized to reproduce the behavior of a 12-bit ADC,
including realistic pedestals and saturation. A similar process is followed to simulate the response
of the photon detectors, given the arrival times of the photons. Functionality exists within LArSoft
to overlay MC-simulated particles with raw digits in the data in order to simulate pileup of cosmics
and other beam interaction particles. The simulated raw digits are then written to compressed
ROOT files for further analysis.
ProtoDUNE Single-Phase Technical Design Report

---