1. 1. Introduction
  2. 2. Neutrinos & Neutrino Oscillations
  3. 3. Neutrino Interactions
  4. 4. Coherent Pion Production
  5. 5. The T2K Experiment
    1. 5.1. The J‑PARC Neutrino Beam
    2. 5.2. Super‑Kamiokande
    3. 5.3. INGRID
    4. 5.4. ND280
      1. 5.4.1. The Pi-Zero Detector
      2. 5.4.2. The Fine Grained Detectors
      3. 5.4.3. The Time Projection Chambers
      4. 5.4.4. The Electromagnetic Calorimeters
      5. 5.4.5. The Magnet & SMRD
    5. 5.5. Common Components
    6. 5.6. Simulations
    7. 5.7. ND280 Software & Reconstruction
  6. 6. A Search for Coherent Pion Production at T2K
  7. 7. Conclusions

1. Introduction

2. Neutrinos & Neutrino Oscillations

3. Neutrino Interactions

4. Coherent Pion Production

5. The T2K Experiment

T2K [1] is a long-baseline neutrino oscillation experiment designed to make precision measurements of θ13, θ23 and Δm232. It does this by firing a beam of mostly νμ along a baseline of 295 km and measuring the resulting νμ and νe energy spectra.

The νμ beam is generated at the J‑PARC facility in Tokai, on the East coast of Japan, and detected by the famous Super‑Kamiokande water-Cherenkov detector near Kamioka. The Super‑Kamiokande detector is well suited for this task. Its 50 kT water volume ensures a sufficient number of the neutrinos which reach it will interact, and reconstructing a primary lepton's Cherenkov ring allows it to determine the timing, direction, energy and flavour of the incident neutrino. Timing and direction are essential to identify the neutrino's origin as being the T2K beam, while flavour and energy are required to measure the energy spectra for νμ and νe separately. T2K is also notable as the world's first off-axis neutrino experiment with the J‑PARC beam directed such that the far-detector is 2.5 ° off the beam's axis.

In addition to the production and detection of neutrinos, the experiment requires several other studies.

The beam's direction, normalisation and energy spectrum all need to be known accurately to correctly predict the un-oscillated flux at Super‑Kamiokande. For example, a shift of the beam direction by 1 mrad would result in a ∼20 MeV change in the peak of the neutrino energy spectrum 1. In addition, the initial beam is not composed purely of νμ but also contains small amounts of νμ, νe and νe. It is important for T2K to know this flavour composition if it is to separate those in the initial beam from those that appear due to oscillations. Since the beam direction is determined by the state of the hardware and, as will be discussed briefly in Section 5.1, neutrino production is difficult to simulate, all of these need to be determined experimentally.

Additionally, as discussed in previous chapters, T2K will need a good knowledge of neutrino-nucleus interactions if it is to successfully understand selection efficiencies, reconstruction of neutrino energies and backgrounds. Since there is so little data over T2K's neutrino energy range this also needs to be done experimentally. Finally, because photons can mimic an electron at Super‑Kamiokande, NC π0 production is a particularly dangerous background for νe detection, so the cross-section for this process in water needs to be accurately determined.

In order to satisfy these requirements T2K built a near-detector facility 280 m downstream of the beam target which contains two separate detectors (Figure 5.1). The on-axis “INGRID” is primarily responsible for measuring the beam's normalisation and direction, while the off-axis “ND280” 2 measures the beam's flavour composition, and studies neutrino-nucleus interactions including NC π0 production.

Figure 5.1: The near detector complex. On the top floor stands the off-axis ND280, shown with the magent open to expose the inner detectors. The floor below holds the horizontal arm of INGRID, in-front of which stands the vertical arm. [1]

5.1. The J‑PARC Neutrino Beam

Figure 5.2: An overview of the J‑PARC accelerator chain []

The Japanese Proton Accelerator Research Complex (J‑PARC) in Tokai on the East coast of Japan is a new facility built on the site where Japan developed their first nuclear reactors (Figure 5.2). The accelerator chain begins with a linear accelerator (LINAC) which boosts H- ions up to a kinetic energy of 181 MeV before passing through charge-stripping foils. The resulting protons are then injected into a 3 GeV Rapid-Cycling Synchrotron (RCS) which in turn feeds a 30 GeV Main Ring (MR) from which the protons destined for the neutrino beamline are extracted.

Fast kicker magnets extract protons from the MR so that they can be bent inside the ring and along the neutrino beamline. As they approach the target the protons are focused and directed downwards by 3.6 ° with respect to the horizontal. Although the direction to Super‑Kamiokande is actually only 1.3 ° downwards, this additional vertical offset forms part of the 2.5 ° off-axis angle of Super‑Kamiokande with respect to the beam.

Figure 5.3: Schematic of the neutrino-beamline from the target onwards.

The protons then collide with the target which is a graphite rod 914 mm long and 26 mm in diameter, with a density of 1.8 g cm-3. It is cooled with helium gas and is designed to operate at a beam power of 750 kW, though to date the maximum delivered has been 200 kW.

The target is located inside the first of three magnetic horns which select and focus the secondary particles (mostly π+) produced in the target (Figure 5.3). The first horn is primarily responsible for “collecting” the pions, while the second and third horns further downstream focus them down a 96 m decay volume.

Positive pions are the most common neutrino-producing particles which enter the decay volume though there are also a significant number of positive kaons. In addition, a small number of wrong-sign pions and kaons are unavoidable and it is also possible for μ+ produced from π+/K+ decays to also decay and contribute to νμ and νe contaminations (Table 5.1).

Neutrino Decay Fraction (%)
νμ π+ → μ+ νμ 73
K+ → μ+ νμ 13
K+ → π0 μ+ νμ 12
K0 → π0 μ+ νμ 2
μ-e- νe νμ 0.02
νμ K-μ- νμ 42
K- → π0 μ- νμ 21
π-μ- νμ 14
K0 → π0 μ+ νμ 13
μ+e+ νe νμ 9
νe K+ → π0 e+ νe 50
π+e+ νe 38
K0 → π- e+ νe 8
μ+e+ νe νμ 4
νe K-e- νe 57
K0 → π+ e- νe 33
π-e- νe 9
μ-e- νe νμ 2
Table 5.1: The fraction of each decay mode which contributes to neutrinos of each neutrino flavour under the flux peak (0.4 GeV ≤ Eν ≤ 1.0 GeV)

The decay volume ends 117.5 m from the target with two muon monitoring detectors which provide valuable data for validating the beam simulation.

Figure 5.4: Simulated neutrino flux at ND280 as a function of neutrino energy.

One of the notable features of the T2K neutrino beam is that it is the world's first off-axis neutrino beam. As was mentioned earlier, the beam is directed such that the Super‑Kamiokande detector is situated 2.5 ° off the beam's primary axis. Directly along the beam's axis, a detector is exposed to a broad range of decay kinematics giving a correspondingly broad neutrino energy spectrum. At fixed angles off-axis however, only a more select range of decay kinematics contribute resulting in a beam with a much narrower energy spectrum (Figure 5.5). There are two good reasons for doing this. The first is that the peak of an off-axis spectrum is lower than for an on-axis beam, allowing T2K to push the peak energy down towards the oscillation maximum energy of ∼600 MeV - giving higher statistics where it matters most. The second advantage comes from the removal of a large flux of higher energy neutrinos. These neutrinos do not contribute to the oscillation signal at Super‑Kamiokande but will generate background interactions - in particular NC interactions - which can be reconstructed in the signal region. For oscillation analyses then, an off-axis beam delivers higher statistics with lower backgrounds.

Figure 5.5: The effect of moving off-axis on the νμ energy spectrum. [2]

At the end of T2K Run 4, in May 2013, the beam had delivered a total of 2.57 × 1020 POT (Table 5.2). This is delivered in ∼5 µs long “spills” of 8 (6 in T2K Run 1) “bunches” each of which is approximately 15.0 ns wide (Figure 5.6).

T2K Run Bunches per Spill Horn Current (kA) Run POT (1020) Total POT (1020)
1 6 250 0.31 0.31
2 8 250 1.12 1.43
3b 8 205 0.22 1.65
3c 8 250 1.37 3.02
4 8 250 3.37 6.40
Table 5.2: Conditions and performance of the T2K beam up to May 2013.
Figure 5.6: (left) The POT delivered by the beam up to the end of T2K Run 4 in May 2013, and the amount of that recorded by ND280. (right) The 8 bunch structure of a beam spill with respect to the beam trigger can be clearly seen in Super‑Kamiokande's initial νe selection. [T2K official plots]

Although designed to take a currents of 320 kA, thus far all T2K physics data has been taken while operating the horns at 250 kA, with the exception of T2K Run 3b when they were operated at 205 kA. The effect of this reduction in horn current is mainly a reduction of around 10 % in the total flux at the peak (Figure 5.7

Figure 5.7: The effect of changing the horn current on the νμ flux at ND280.

5.2. Super‑Kamiokande

Super‑Kamiokande is perhaps the most famous experiment in the history of neutrino physics and played an important part in establishing evidence for both solar and atmospheric neutrinos oscillating. It also holds the world's leading limit on proton decay [3]. Now in it's fourth running configuration, Super‑Kamiokande is being re-utilised as the far detector for T2K.

Figure 5.8: Schematic of the Super‑Kamiokande detector. []

The device itself is a large water-Cherenkov detector within a vertical cylindrical tank. A stainless steel frame forming a smaller cylinder within the tank separates the volume into two regions known as the inner and outer detectors. The entire volume is filled with 50 kT of ultra-pure water though, in practice, the fiducial mass is somewhat smaller.

The inner detector is the target region, and it is instrumented on all sides by 11129 inward-facing photo-multiplier tubes (PMTs) which cover ∼40 % of the surface area. Particles produced in the inner detector which are above the Cherenkov threshold produce rings of light which can be detected by these PMTs and from which the event can be reconstructed. The position of the interaction vertex can be determined from the position, shape and diameter of the ring. The thickness of the light ring is a measure of the distance the particle travelled, with exiting particles generating a completely filled circle. Finally, primitive particle identification for separating showering and non-showering particles (effectively electron-muon separation) can be done using the diffuseness of the outer part of the ring (Figure 5.9). It is the fact that this particle identification cannot discriminate between a showering electron or photon that makes NC π0 production such an important background for νe selection.

Figure 5.9: Event displays from T2K beam interactions in Super‑Kamiokande. The outer part of the Cherenkov ring resulting from an electron (left) is much “fuzzier” than one from a muon (right). [T2K official plots]

The outer detector was designed as a veto for incoming (mostly cosmic-ray) particles and is essentially a 2 m thick region surrounding the inner detector. The outwards facing wall is also instrumented with PMTs albeit much more sparsely - the tubes are fewer and smaller. In T2K however, because the direction and timing requirements are already quite effective at selecting beam events, interactions in the outer detector can also contribute to the oscillation analysis.


The Interactive Neutrino Grid (INGRID) [4] is the on-axis near detector primarily designed to constrain the beam normalisation and direction. It achieves this by placing identical modules at a variety of off-axis angles in a plane perpendicular to the beam. The direction can then be determined from the beam's profile across this plane, and the normalisation from the rate of interactions.

The baseline design consists of 16 identical modules (Figure 5.10). 14 of these are placed in a 10 m × 10 m cross-shape, centred on the beam's axis. This is formed of a vertical stack of 7 modules, which sits in front of a horizontal row of a further 7 modules. In between these two sets are two modules located above the horizontal row and on either side of the vertical stack. These off-axis modules assist in characterising any asymmetries in the beam shape - which is important for an off-axis beam.

Figure 5.10: The layout of the INGRID modules. 14 modules are arranged in a cross-shape: 7 vertical modules in front of 7 horizontal modules. 2 additional modules are placed in off-axis positions. The proton module is located in the centre of the cross, between the vertical and horizontal arms. [1]

Each standard module is primarily constructed of a central tracking region surrounded by veto planes (Figure 5.11). The tracking region is approximately 1.2 m square and is comprised of 9 × 65 mm thick iron plates, separating the first 10 of 11 planes of plastic scintillator. Each of these planes contains 2 layers of 10 mm × 50 mm bars arranged in an x-y configuration. The central tracking region is surrounded above, below and on the sides by veto planes of 22 plastic scintillator bars aligned perpendicular to the beam.

Figure 5.11: An exploded view of one of the standard INGRID modules. The left diagram shows the central tracking region consisting of 9 iron plates and 11 x-y scintillator planes. The right diagram shows the addition of the surrounding veto planes. [1]

Because of the large amount of iron dead-material and low tracking granularity, INGRID modules are unable to reconstruct much more than a primary μ-, though this is sufficient for its task of profiling the beam. After the construction of INGRID, additional components were assembled in to a scintillator-only module which is referred to as the “proton module” - so named because its totally active design allows it, for example, to also reconstruct a proton from a CC QE interaction. This module was placed centrally and between the vertical and horizontal arms (inline with the off-axis modules).

During T2K Runs 1 and 2 INGRID consistently recorded a mean event rate of 1.5 / 1014POT, agreeing with the predicted rate from simulation within uncertainty: data / simulation = 1.06 ± 0.04. Fitting the number of events in each module as a function of off-axis angle allowed it to measured the direction of the T2K beam to an accuracy of around 0.4 mrad, well below the target resolution of 1.0 mrad (Figure 5.12).

Figure 5.12: (left) An example INGRID beam profile measurement showing number of events in Beam Run 32 as a function of horizontal distance from the centre. (right) The beam direction is measured and stable well within the experiment's target of ±1.0 mrad. [4]

5.4. ND280

ND280 is the experiment's off-axis near-detector, located in the near-detector facility 280m downstream of the target, and slightly downstream of INGRID.

As discussed earlier, the primary requirements of the detector are:

The first goal requires a target of water on which neutrinos can interact to produce π0s. Almost immediately after creation these pions will decay via π0γγ requiring the presence of high density material to encourage the photons to shower if they are to be detected.

Meanwhile, the study of CC νμ interactions requires precise measurements of the primary lepton's angle and momentum, in addition to the ability to observe secondary particles emerging from the nucleus. The former drives the design towards lower-densities and the latter requires an active target.

In an ideal world a neutrino detector should have a single target region but, unfortunately for ND280, the dual requirements of efficiently converting photons while maintaining good position resolution cannot be simultaneously satisfied. The same is true for the requirements of having a water target, and an active target. As a result the detector is a compromise of multiple sub-systems.

Another significant influence on the design of ND280 is the need for a magnet. Due to their penetrating nature, the momentum of muons cannot practically be measured calorimetrically in such a confined space but must instead be done via curvature in a magnetic field. Such a field also brings the added capability to distinguish μ- from μ+ - valuable when νμ make up ∼6% of the total flux and the majority of cosmic-ray muons are μ+.

Figure 5.13: An exploded view of the ND280. The central 'basket' region contains the P0D, TPCs, FGDs and Ds ECal. They are surrounded by the Barrel ECal and P0D ECal, and the whole detector is enclosed by the UA1 magnet. [1]

The design conceived to meet all these requirements (Figure 5.13) placed the primary detectors in a central region, known as the “basket”, and surrounded them on all but the upstream-side with electro-magnetic calorimeters. Upstream-most in the basket is the Pi-Zero Detector (P0D) which combines water, metal foils and plastic scintillator to provide the target, photon conversion and reconstruction required for the study of NC π0s. Downstream of the P0D is the “tracker”, in which three Time Projection Chambers (TPCs) sandwich two plastic scintillator detectors referred to as Fine Grained Detectors (FGDs). The FGDs act as the target for ND280's interaction studies - offering sufficient mass for a reasonable interaction rate in addition to an active target to access some of the activity near an interaction vertex. The argon-gas TPCs give excellent tracking, for direction and momentum reconstruction, and particle identification from dE/dx. They serve both the FGDs and the P0D, tracking any forward-escaping particles.

The surrounding electromagnetic calorimeters (ECals) consist of two types with very different remits. The more substantial Tracker ECals provide particle identification and energy measurements for final-state particles leaving the tracker. The much more limited P0D ECals which surround the P0D are there only to catch muons and photons escaping at high angles, and to veto incoming backgrounds.

All these systems are enclosed within a magnet which generates a 0.2 T magnetic field and houses the Side Muon Ranging Detector (SMRD) - planes of scintillator between the plates of the iron yoke - that acts primarily as a cosmic trigger. However encasing the detector within the former UA1 and NOMAD magnet which was available substantially restricts the available space to a box approximately 7 m long × 3.6 m wide × 3.5 m high.

5.4.1. The Pi-Zero Detector

The P0D [5] is the upstream-most detector in the ND280 basket with a central water target region sandwiched by two ECal regions (Figure 5.14). The water targets alternate with plastic scintillator modules and brass foils to encourage the photons to convert. The ECal regions meanwhile alternate lead sheets between the plastic scintillator modules to maximise photon conversion and containment. As a final feature, the water targets of the P0D can be filled and drained during running to aid in separating the water contribution to interactions in the fiducial volume. In its most recent technical note [6] the P0D reported a π0 reconstruction efficiency of 3.6 % (though updated analyses are reporting improvements upon this).

Figure 5.14: A schematic diagram showing a side-on view of the P0D's design. [1]

5.4.2. The Fine Grained Detectors

The two Fine Grained Detectors (FGDs) [7] are the most important target masses for studies of neutrino-nucleus interactions. The upstream FGD, referred to as FGD 1 consists only of plastic scintillator modules, while the downstream FGD, FGD 2, includes six 25 mm thick water targets in addition to plastic scintillator bars.

The FGD plastic scintillator bars are sufficiently dense to provide a good interaction rate (∼1.6 / 1016POT), but the bars' 9.61 mm square cross-section also give reasonably good position resolution. This allows FGD 1 in particular to give good vertex position resolution and to reconstruct final-state particles in addition to the primary lepton. For particles which are contained, particle identification is also possible by comparing the range and energy deposited. While predominantly composed of carbon-12 the scintillator, coatings, fibres and glue in FGD 1's active region also bring along small fractions of oxygen, hydrogen, titanium, silicon and nitrogen, which collectively average to give an Aeffective = 12.10 [8].

FGD 1 has 30 plastic scintillator layers, each with 192 × 1864.3 mm long bars, alternating between horizontal and vertical orientations. There are 14 identical layers in FGD 2, two of which lie either side of each water target. Because of its water targets FGD 2 does not have quite the same vertex resolution, but it does provide a comparable environment to FGD 1. This will enable studies of neutrino-carbon interactions from FGD 1 to be compared with neutrino-water interactions in FGD 2 for a better understanding of interactions in Super‑Kamiokande.

5.4.3. The Time Projection Chambers

The three Time Projection Chambers (TPCs) [9] are placed immediately downstream of each of the three target sub-detectors (P0D, FGD 1 and FGD 2) to provide tracking of the particles which emerge. For all but the highest-angle tracks they achieve position resolutions < 1.0 mm, allowing them to accurately determine the charge and momentum of particles from their curvature in the magnetic field (with transverse momentum resolution σp/p2 ∼ 0.1 %).

They also provide particle identification by vertically segmenting tracks, and calculating the mean rate of energy loss (dE/dx) from the path-length corrected energy deposited in each segment (after truncating off the highest 30 %). This is then compared with the expected value from simulation to identify the particle (Figure 5.15). The resolution on the mean dE/dx is better than 8 % for MIPs, making the TPCs the ND280's primary particle identification tool [9].

Figure 5.15: The truncated mean dE/dx for particles in T2K Run 1 data, with the fitted curves for simulated muons, protons, pions and electrons overlaid. [1]

The TPCs are all constructed identically (Figure 5.16). A central cathode separates two drift chambers containing an argon-based gas, and readout electronics are placed on the lateral sides. Each drift chamber has twelve 342 mm × 359 mm “micromegas” readout tiles, each of which provides a grid of 48 × 36 pads.

Figure 5.16: Diagram showing the design of a TPC. [1]

5.4.4. The Electromagnetic Calorimeters

The ND280 ECals are really two different sub-systems: the Tracker ECal, comprising the Downstream ECal (Ds ECal) and Barrel ECals, surrounds the tracker (FGDs and TPCs), while the P0D ECal surrounds the P0D.

The Tracker ECal was designed to convert photons generated in the FGDs, provide particle identification (particularly for muon-pion separation), make energy measurements of showering particles, and track anything escaping the FGDs at high-angle. Because this requires a full reconstruction, the Tracker ECal modules have many thin lead sheets separated by layers of plastic scintillator bars which alternate between two views perpendicular to the tracker.

The P0D ECal is significantly less sophisticated since the P0D is capable of converting most of the photons generated within it, and has the tracker downstream for more penetrating particles. Instead the P0D ECal serves only to catch muons and photons escaping the P0D at very high-angles, and to veto incoming activity. These tasks can be performed by just 6 layers of uni-directional plastic scintillator, separated by much thicker lead sheets.

One thing in common between the Barrel ECal and P0D ECal is the need for 6 modules to cover the 4 sides of the basket. Each of the lateral sides can be covered by one module, but the top and bottom sides require 2 modules each to allow the magnet to open.

Module(s) Layer Bars per
Bar length
Ds ECal 17 Horizontal 50 2.04 2 1.75
17 Vertical 50 2.04 2
Barrel ECal
16 Perpendicular 96 1.52 1
15 Parallel 38 3.84 2
Barrel ECal
16 Perpendicular 96 2.36 1
15 Parallel 57 3.84 2
P0D ECal
6 Parallel 38 2.34 1 4.0
P0D ECal
6 Parallel 96 2.34 1
Table 5.3: Specifications of the ECal modules. Layer names relate to the direction along which the longest dimension of the bars lie with respect to to the beam or the world.

5.4.5. The Magnet & SMRD

A muon's energy cannot practically be measured calorimetrically in such a confined space, since they are so highly penetrating, so this must be done from curvature in a magnetic field. A magnetic field also gives ND280 the ability to measure the charge of a particle track - a valuable capability for distinguishing μ- from μ+ and π+ from π-. For this reason T2K acquired a 0.2 T magnet from CERN which had previously been used in the NOMAD and UA1 experiments.

The magnet's coil runs vertically and parallel to the beam to generate a magnetic field which, from the view of an observer standing upstream of the detector, would point horizontally to the right. This results in negatively charged particles which are initially travelling downstream being bent downwards. The coil is surrounded by a yoke assembled from plates of iron.

Finally, some of the gaps between the iron plates are instrumented by plastic scintillator planes, with optical fibres running through them. These planes form the Side Muon Ranging Detector (SMRD). In concert with the P0D and Ds ECal, the SMRD primarily acts as the detector's cosmic trigger and veto, though its hits from escaping muons can also be associated with a reconstructed track from the inner detectors.

5.5. Common Components

The P0D, FGDs, ECals, SMRD and INGRID all contain plastic scintillator as their active components. The extruded plastic scintillator was produced at Fermilab and is composed of polystyrene doped with wavelength shifting scintillators PPO (1 %) and POPOP (0.03 %) which results in blue light emission (420 nm). The plastic was co-extruded with a 0.25 mm thick reflective layer of TiO2 to keep scintillation light in and external light out.

This scintillation light is then collected and transported by 1 mm diameter optical fibres which are strung down 3 mm diameter extruded holes in the centre of the bars (or laid in a groove in the SMRD planes). These fibres, manufactured by Kuraray, are also wavelength shifting and convert the blue light spectrum to a green (476 nm) that is well matched to the light sensors attached to the ends.

These light sensors are Multi-Pixel Photon Counters (MPPC) [10] manufactured by Hamamatsu which, unlike photomultiplier tubes, can still be operated in a magnetic field. They consist of 667 pixels arranged accross a 1.3 mm × 1.3 mm active area, and boast a photon detection efficiency of ∼20 % for the light received from the optical fibres.

5.6. Simulations

As with all complex experiments it is necessary for T2K to run simulations in order to understand how the measurements recorded by its detectors correspond to the physics which occurred.

The J‑PARC neutrino beam is simulated for the purpose of providing flux predictions [2] to the near and far detectors using software collectively known as “JNUBEAM”. The simulated beamline geometry includes the target, horns, decay volume and beam monitors. 30 GeV protons are directed onto the graphite target, and the resulting particles are tracked until they decay to neutrinos or stop. The simulation of primary interactions inside the target are based on data from the NA61/SHINE [11] experiment, secondary interactions are simulated with FLUKA [12], and interactions outside the target use GEANT3/GCALOR [13].

The simulation of the near detector geometries is done in GEANT4 [14] as is the passage of final-state particles emerging from neutrino-nucleus interactions. The neutrino interactions themselves can be simulated with either the GENIE or NEUT generators (see Chapter 3), which propagate neutrinos from the beam flux simulation through the geometry exported from the GEANT4 simulation. After the particles have been propagated through the near detectors a custom software package, called “ELECSIM”, simulates the response of the detector and its electronics to the energy deposited. In the scintillator sub-detectors this includes simulating the light emitted in response to energy deposition, the light's transport through the bar and down the optical fibres, the response of the MPPCs to that light and the electronics chain thereafter. For the TPCs ELECSIM simulates the electron drift, response of the micromegas and, again, the electronics chain which follows.

Super‑Kamiokande's “SKDETSIM” is a fortran-based software package which is responsible for the geometry and final-state particle simulations. Again, the neutrino interactions can be generated with GENIE or NEUT using the neutrino flux from the beam simulation.

5.7. ND280 Software & Reconstruction

The first processing stage in the ND280 software chain is calibration, and is the last stage at which simulation and data are treated differently. For data the goal is to correct for the various effects introduced by the scintillator bars, optical fibres and readout electronics. For simulation the goal is to replicate the output from this process. In both cases the result is a collection of calibrated “hits” with a position, time and information on the energy deposited (often referred to as “charge”).

The ND280 reconstruction software is essentially a two-stage process as a result of the very different time/space/charge information provided by the different sub-detector systems, which makes a unified reconstruction algorithm unsuitable. Instead, each sub-detector has its own dedicated reconstruction algorithm which groups associated hits together, fits the resulting objects, and calculates properties for use by analysers (such as the TPC particle identification, or FGD vertex activity). The last of these to be run is the TPC reconstruction, which needs to find matching objects in neighbouring detectors in order to determine the time at which the track was created.

The results from the individual sub-detector reconstructions are then passed to a “global” reconstruction, which matches objects with compatible start/end points and directions between detectors. These combined objects are re-fitted using a Kalman filter from the Recpack [15] toolkit, utilising geometric and magnetic field models, to create final “reconstructed objects” with position, direction and momentum measurements..

6. A Search for Coherent Pion Production at T2K

7. Conclusions


  1. This can be determined from Figure 5.5, where the fluxes at 2.5 ° and 2.0 ° peak at approximately 0.6 GeV and 0.75 GeV respectively.
  2. T2K's naming convention for “ND280” is ambiguous, and is used to refer both to the entire near detector complex and specifically to the off-axis near detector. In this thesis I refer to the former as the “near detector complex”, reserving “ND280” exclusively for the off-axis near detector.


  1. The T2K Collaboration
    The T2K Experiment
    NIM A, Volume 659, p106-135, 2011
  2. The T2K Collaboration
    T2K neutrino flux prediction
    Physical Review D, Volume 87, p012001, 2013
  3. The Super-Kamiokande Collaboration
    Search for Proton Decay via p → e+π0 and p → μ+π0 in a Large Water Cherenkov Detector
    Physical Review Letters, Volume 102, p141801, 2009
  4. The T2K Collaboration
    Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector
    NIM A, Volume 694, p211-223, 2012
  5. The T2K ND280 P0D Collaboration
    The T2K ND280 Off-Axis Pi-Zero Detector
    NIM A, Volume 686, p48-63, 2012
  6. G. Lopez et al.
    Measurement of NC π0 Production with the T2K π0 Detector (P0D)
    T2K Technical Note 056, Version 6
  7. The T2K ND280 FGD Collaboration
    The T2K fine-grained detectors
    NIM A, Volume 696, p1-31, 2012
  8. S. Oser
    Elemental Composition of the FGD XY Modules
    T2K Technical Note 091, Version 1.2
  9. The T2K ND280 TPC collaboration
    Time Projection Chambers for the T2K Near Detectors
    NIM A, Volume 637, p25-46, 2011
  10. A. Vacheret et al.
    Characterization and Simulation of the Response of Multi Pixel Photon Counters to Low Light Levels
    NIM A, Volume 656, p69-83, 2011
  11. The NA61/SHINE Collaboration
    Pion emission from the T2K replica target: Method, results and application
    NIM A, Volume 701, p99-114, 2013
  12. G. Battistoni et al.
    The FLUKA code: description and benchmarking
    AIP Conference Proceedings, Volume 896, p31-49, 2007
  13. R. Brun, F. Carminati, and S. Giani
    GEANT: Detector description and simulation tool
    CERN-W5013, 1994
  14. The Geant4 Collaboration
    Geant4 - a simulation toolkit
    NIM A, Volume 506, p250-303, 2003
  15. A. Cervera-Villanueva et al.
    “RecPack” a reconstruction toolkit
    NIM A, Volume 534, p180-183, 2004