
Particle physics has an ambitious and broad experimental programme for the
coming decades. This programme requires large investments in detector hardware,
either to build new facilities and experiments, or to upgrade existing ones.
Similarly, it requires commensurate investment in the R&D of software to
acquire, manage, process, and analyse the shear amounts of data to be recorded.
In planning for the HLLHC in particular, it is critical that all of the
collaborating stakeholders agree on the software goals and priorities, and that
the efforts complement each other. In this spirit, this white paper describes
the R&D activities required to prepare for this software upgrade.

The CoGeNT and CRESST WIMP direct detection experiments have recently
observed excesses of nuclear recoil events, while the DAMA/LIBRA experiment has
a long standing annual modulation signal. It has been suggested that these
excesses may be due to light mass, m_chi ~ 510 GeV, WIMPs. The Earth's motion
with respect to the Galactic rest frame leads to a directional dependence in
the WIMP scattering rate, providing a powerful signal of the Galactic origin of
any recoil excess. We investigate whether direct detection experiments with
directional sensitivity have the potential to observe this anisotropic
scattering rate with the elastically scattering light WIMPs proposed to explain
the observed excesses. We find that the number of recoils required to detect an
anisotropic signal from light WIMPs at 5 sigma significance varies from 7 to
more than 190 over the set of target nuclei and energy thresholds expected for
directional detectors. Smaller numbers arise from configurations where the
detector is only sensitive to recoils from the highest speed, and hence most
anisotropic, WIMPs. However, the event rate above threshold is very small in
these cases, leading to the need for large experimental exposures to accumulate
even a small number of events. To account for this sensitivity to the tail of
the WIMP velocity distribution, whose shape is not well known, we consider two
exemplar halo models spanning the range of possibilities. We also note that for
an accurate calculation the Earth's orbital speed must be averaged over. We
find that the exposures required to detect 10 GeV WIMPs at a WIMPproton
crosssection of 10^4 pb are of order 10^3 kg day for a 20 keV energy
threshold, within reach of planned directional detectors. Lower WIMP masses
require higher exposures and/or lower energy thresholds for detection.

This paper presents an investigation of the use of interest point detection
algorithms from image processing applied to reconstruction of interactions in
high granularity tracking detectors. Their purpose is to extract keypoints from
the data as input to higher level reconstruction algorithms, replacing the role
of human operators in event selection and reconstruction guidance. Simulations
of nu_mu charged current events in a small liquid argon time projection chamber
are used as a concrete example of a modern high granularity tracking detector.
Data from the simulations are used to characterize the localization of interest
points to physical features and the efficiency of finding interest points
associated with the primary vertex and track ends is measured. A high degree of
localization is found, with 93% of detected interest points found within 5mm of
a physical feature. Working in two 2D projections, the primary vertex and both
track ends are found in both projections in 85% of events. It is also shown
that delta electrons can be detected.

Direct detection experiments have reached the sensitivity required to detect
dark matter WIMPs. Demonstrating that a putative signal is due to WIMPs, and
not backgrounds, is a major challenge however. The direction dependence of the
WIMP scattering rate provides a potential WIMP `smoking gun'. If the WIMP
distribution is predominantly smooth, the Galactic recoil distribution is
peaked in the direction opposite to the direction of Solar motion. Previous
studies have found that, for an ideal detector, of order 10 WIMP events would
be sufficient to reject isotropy, and rule out an isotropic background. We
examine how the median recoil direction could be used to confirm the WIMP
origin of an anisotropic recoil signal. Specifically we determine the number of
events required to reject the null hypothesis that the median direction is
random (corresponding to an isotropic Galactic recoil distribution) at 95%
confidence. We find that for zero background 31 events are required, a factor
of roughly 2 more than are required to simply reject isotropy. We also
investigate the effect of a nonzero isotropic background. As the background
rate is increased the number of events required increases, initially fairly
gradually and then more rapidly, once the signal becomes subdominant. We also
discuss the effect of features in the speed distribution at large speeds, as
found in recent high resolution simulations, on the median recoil direction.

We study the consequences of limited recoil sense reconstruction on the
number of events required to reject isotropy and detect a WIMP signal using a
directional detector. For a constant probability of determining the sense
correctly, 3d readout and zero background, we find that as the probability is
decreased from 1.0 to 0.75 the number of events required to reject isotropy
using the mean angle statistic is increased by a factor of a few. As the
probability is decreased further the number of events required using this
statistic increases sharply, and in fact isotropy can be rejected more easily
by discarding the sense information and using axial statistics. This however
requires an order of magnitude more events than vectorial data with perfect
sense determination. We also consider energy dependent probabilities of
correctly measuring the sense, 2d readout and nonzero background. Our main
conclusion regarding the sense determination is that correctly determining the
sense of the abundant, but less anisotropic, low energy recoils is most
important for minimising the number of events required.

We study the dependence of the exposure required to directly detect a WIMP
directional recoil signal on the capabilities of a directional detector.
Specifically we consider variations in the nuclear recoil energy threshold, the
background rate, whether the detector measures the recoil momentum vector in 2
or 3 dimensions and whether or not the sense of the momentum vector can be
determined. We find that the property with the biggest effect on the required
exposure is the measurement of the momentum vector sense. If the detector
cannot determine the recoil sense, the exposure required is increased by an
order of magnitude for 3d readout and two orders of magnitude for 2d
readout. For 2d readout the required exposure, in particular if the senses
can not be measured, can be significantly reduced by analyzing the reduced
angles with the, time dependent, projected direction of solar motion
subtracted. The background rate effectively places a lower limit on the WIMP
crosssection to which the detector is sensitive; it will be very difficult to
detect WIMPs with a signal rate more than an order of magnitude below the
background rate. Lowering the energy threshold also reduces the required
exposure, but only for thresholds above 20 keV.

The direction dependence of the WIMP direct detection rate provides a
powerful tool for distinguishing a WIMP signal from possible backgrounds. We
study the the number of events required to discriminate a WIMP signal from an
isotropic background for a detector with 2d readout using nonparametric
circular statistics. We also examine the number of events needed to i) detect a
deviation from rotational symmetry, due to flattening of the Milky Way halo and
ii) detect a deviation in the mean direction due to a tidal stream. If the
senses of the recoils are measured then of order 2070 events (depending on the
plane of the 2d read out) will be sufficient to reject isotropy of the raw
recoil angles at 90% confidence. If the senses can not be measured these number
increase by roughly two orders of magnitude (compared with an increase of one
order of magnitude for the case of full 3d readout). The distributions of the
reduced angles, with the (time dependent) direction of solar motion subtracted,
are far more anisotropic, however, and if the isotropy tests are applied to
these angles then the numbers of events required are similar to the case of 3d
readout. A deviation from rotational symmetry will only be detectable if the
Milky Way halo is significantly flattened. The deviation in the mean direction
due to a tidal stream is potentially detectable, however, depending on the
density and direction of the stream.

The direction dependence of the event rate in WIMP direct detection
experiments provides a powerful tool for distinguishing WIMP events from
potential backgrounds. We use a variety of (nonparametric) statistical tests
to examine the number of events required to distinguish a WIMP signal from an
isotropic background when the uncertainty in the reconstruction of the nuclear
recoil direction is included in the calculation of the expected signal. We
consider a range of models for the Milky Way halo, and also study rotational
symmetry tests aimed at detecting nonsphericity/isotropy of the Milky Way
halo. Finally we examine ways of detecting tidal streams of WIMPs. We find that
if the senses of the recoils are known then of order ten events will be
sufficient to distinguish a WIMP signal from an isotropic background for all of
the halo models considered, with the uncertainties in reconstructing the recoil
direction only mildly increasing the required number of events. If the senses
of the recoils are not known the number of events required is an order of
magnitude larger, with a large variation between halo models, and the recoil
resolution is now an important factor. The rotational symmetry tests require of
order a thousand events to distinguish between spherical and significantly
triaxial halos, however a deviation of the peak recoil direction from the
direction of the solar motion due to a tidal stream could be detected with of
order a hundred events, regardless of whether the sense of the recoils is
known.