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 HL-LHC 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.

We study the existence of $S^1$-equivariant characteristic classes on certain natural infinite rank bundles over the loop space $LM$ of a manifold $M$. We discuss the different $S^1$-equivariant cohomology theories in the literature and clarify their relationships. We attempt to use $S^1$-equivariant Chern-Weil techniques to construct $S^1$-equivariant characteristic classes. The main result is the construction of a sequence of $S^1$-equivariant characteristic classes on the total space of the bundles, but these classes do not descend to the base $LM$. Nevertheless, we conclude by identifying a class of bundles for which the $S^1$-equivariant first Chern class does descend to $LM$.

We study the Maslov index as a tool to analyze stability of steady state solutions to a reaction-diffusion equation in one spatial dimension. We show that the path of unstable subspaces associated to this equation is governed by a matrix Riccati equation whose solution $S$ develops singularities when changes in the Maslov index occur. Our main result proves that at these singularities the change in Maslov index equals the number of eigenvalues of $S$ that increase to $+\infty$ minus the number of eigenvalues that decrease to $-\infty$.

We have created 3D models of the CMS detector and particle collision events in SketchUp, a 3D modelling program. SketchUp provides a Ruby API which we use to interface with the CMS Detector Description to create 3D models of the CMS detector. With the Ruby API, we also have created an interface to the JSON-based event format used for the iSpy event display to create 3D models of CMS events. These models have many applications related to 3D representation of the CMS detector and events. Figures produced based on these models were used in conference presentations, journal publications, technical design reports for the detector upgrades, art projects, outreach programs, and other presentations.

We present and compare the predictions of various cosmic-ray Monte Carlo models for the energy (dE/deta) and particle (dN/deta) flows in p-p, p-Pb and Pb-Pb collisions at sqrt(s) = 14, 8.8, and 5.5 TeV respectively, in the range covered by forward LHC detectors like CASTOR or TOTEM (5.2<|eta|<6.6) and ZDC or LHCf (|eta|>8.1 for neutrals).

In the course of its operation, the EGRET experiment detected high-energy gamma ray sources at energies above 100 MeV over the whole sky. In this communication, we search for large-scale anisotropy patterns among the catalogued EGRET sources using an expansion in spherical harmonics, accounting for EGRET's highly non-uniform exposure. We find significant excess in the quadrupole and octopole moments. This is consistent with the hypothesis that, in addition to the galactic plane, a second mid-latitude (5^{\circ} < |b| < 30^{\circ}) population, perhaps associated with the Gould belt, contributes to the gamma ray flux above 100 MeV.

The properties of cosmic rays with energies above 10**6 GeV have to be deduced from the spacetime structure and particle content of the air showers which they initiate. In this review we summarize the phenomenology of these giant air showers. We describe the hadronic interaction models used to extrapolate results from collider data to ultra high energies, and discuss the prospects for insights into forward physics at the LHC. We also describe the main electromagnetic processes that govern the longitudinal shower evolution, as well as the lateral spread of particles. Armed with these two principal shower ingredients and motivation from the underlying physics, we provide an overview of some of the different methods proposed to distinguish primary species. The properties of neutrino interactions and the potential of forthcoming experiments to isolate deeply penetrating showers from baryonic cascades are also discussed. We finally venture into a terra incognita endowed with TeV-scale gravity and explore anomalous neutrino-induced showers.

Using data from the SUGAR and the AGASA experiments taken during a 10 yr period with nearly uniform exposure to the entire sky, we search for large-scale anisotropy patterns in the arrival directions of cosmic rays with energies > 10^{19.6} eV. We determine the angular power spectrum from an expansion in spherical harmonics for modes out to \ell=5. Based on available statistics, we find no significant deviation from isotropy. We compare the rather modest results which can be extracted from existing data samples with the results that should be forthcoming as new full-sky observatories begin operation.

One of the possible consequences of the existence of extra degrees of freedom beyond the electroweak scale is the increase of neutrino-nucleon cross sections ($\sigma_{\nu N}$) beyond Standard Model predictions. At ultra-high energies this may allow the existence of neutrino-initiated extensive air showers. In this paper, we examine the most relevant observables of such showers. Our analysis indicates that the future Pierre Auger Observatory could be potentially powerful in probing models with large compact dimensions.

We describe a scalable parallelization of Geant4 using commodity hardware in a collaborative effort between the College of Computer Science and the Department of Physics at Northeastern University. The system consists of a Beowulf cluster of 32 Pentium II processors with 128 MBytes of memory each, connected via ATM and fast Ethernet. The bulk of the parallelization is done using TOP-C (Task Oriented Parallel C), software widely used in the computational algebra community. TOP-C provides a flexible and powerful framework for parallel algorithm development, is easy to learn, and is available at no cost. Its task oriented nature allows one to parallelize legacy code while hiding the details of interprocess communications. Applications include fast interactive simulation of computationally intensive processes such as electromagnetic showers. General results motivate wider applications of TOP-C to other simulation problems as well as to pattern recognition in high energy physics.