New results on νμ → ντ appearance with the OPERA experiment in the CNGS beam

The OPERA neutrino experiment is designed to perform the first observation of neutrino oscillations in direct appearance mode in the νμ → ντ channel, via the detection of the τ-leptons created in charged current ντ interactions. The detector, located in the underground Gran Sasso Laboratory, consists of an emulsion/lead target with an average mass of about 1.2 kt, complemented by electronic detectors. It is exposed to the CERN Neutrinos to Gran Sasso beam, with a baseline of 730 km and a mean energy of 17 GeV. The observation of the first ντ candidate event and the analysis of the 2008-2009 neutrino sample have been reported in previous publications. This work describes substantial improvements in the analysis and in the evaluation of the detection efficiencies and backgrounds using new simulation tools. The analysis is extended to a sub-sample of 2010 and 2011 data, resulting from an electronic detector-based pre-selection, in which an additional ντ candidate has been observed. The significance of the two events in terms of a νμ → ντ oscillation signal is of 2.40σ.


Introduction
Flavour transitions between neutrino species were predicted nearly 50 years ago [1][2][3]. In 1998 the Super-Kamiokande experiment observed a strong deficit of atmospheric muon neutrinos in the data, and interpreted it as a result of such transitions [4][5][6]. A few years later, after being investigated for decades by real-time and geochemical experiments [7][8][9][10][11][12], the so-called solar neutrino problem could also be interpreted (after the SNO results [13])

JHEP11(2013)036
in terms of neutrino oscillations in presence of the MSW effect [14] in the Sun's matter. The disappearance of ν µ in atmospheric neutrinos [4][5][6][15][16][17] was confirmed at accelerator-based long baseline experiments [18][19][20]. The observation of an appearance of neutrinos from oscillations consistent with the disappearance results is still missing. The OPERA experiment [21][22][23][24] has the capability of detecting the appearance of a small ν τ component in a ν µ beam. The fact that the ν µ → ν e oscillation could not be the dominant reason of the deficit of ν µ was also established by the measurement of atmospheric ν e rates and confirmed by nuclear reactor experiments at short baselines [25,26] in the late nineties. Only recently, the probability amplitude of the ν µ → ν e transition, governed by the θ 13 mixing angle, has been measured at a long baseline neutrino beam (T2K [27]). The first evidence of the disappearance ofν e has come from reactor experiments at short baselines and θ 13 has been measured (Daya Bay [28], Double Chooz [29], RENO [30]). Super-Kamiokande also recently reported a statistical evidence of ν µ → ν τ transitions [31] in its atmospheric neutrino sample but with a modest signal to background noise due to the difficulty in cleanly reconstructing ν τ events in a water-Cherenkov detector. A positive evidence from OPERA would definitely prove that the ν µ → ν τ transition is the reason of ν µ deficit and it is the dominant mechanism at the atmospheric scale, providing an essential constraint for the establishment of the 3-flavor mixing scheme. Furthermore the measurement can constrain phenomenological models such as neutrino Non Standard Interactions (NSI) [32,33] or sterile neutrinos, which, for certain choices of the parameters, result in modifications to the expected ν τ appearance rate.
After two years of technical and low-intensity runs (2006)(2007) the OPERA detector in the Gran Sasso underground laboratory (LNGS) has been exposed from 2008 to 2012 to the high-energy CERN Neutrinos to Gran Sasso (CNGS) beam [34][35][36]. The detection of the first CNGS neutrino interactions in OPERA was reported in [37,38] and the observation of a first ν τ candidate event was presented in [39], using a sample corresponding to 1.9 × 10 19 protons on target (pot). An update based on the data collected in years 2008-09 (4.88 × 10 19 pot) was presented in [40]. This paper presents the progress of the ν τ appearance search, through a re-evaluation of the efficiencies and of the backgrounds based on a full simulation of the complete analysis chain. Moreover a larger data sample has been exploited which allowed the observation of a second ν τ candidate event.
After an introduction explaining the operation of the experiment (section 2), the data selection (section 3) and the simulation (section 4), the results (section 5) will be described, giving particular emphasis to the second ν τ candidate. The significance of the observation of two candidates given the present statistics and the estimated level of background is addressed in section 6 in terms of exclusion of the null-hypothesis.

Detector and neutrino beam
The OPERA detector [41] is designed to tackle a challenging task: achieving micro-metric tracking accuracy over a very large detector volume spanning about (6.5 × 6.5 × 8) m 3 . The scale of the required granularity is set by the flight length of τ leptons, which for the CNGS beam has a roughly exponential distribution with a mean of about 600 µm. This JHEP11(2013)036 Figure 1. A picture of the OPERA detector. CNGS neutrinos travel from left to right. The (right-handed) reference frame is oriented such that: the y-axis is perpendicular to the hall floor and pointing up; the z-axis is orthogonal to the brick walls and oriented as the incoming neutrinos. The angle between the neutrino direction and the z-axis projected into the yz plane is of 58 mrad. challenge was addressed by using nuclear emulsion based trackers. Another important constraint is related to the practical impossibility to analyze the full emulsion surface O(0.1 km 2 ), even with the state-of-the art automatic scanning technology. This, together with other constraints, resulted in a highly modular target made of units based on the Emulsion Cloud Chamber (ECC) technique, hereafter called bricks, interspersed with pair of planes of horizontal and vertical scintillator strips (called Target Tracker or TT [42]) that allow locating with a centimetric resolution the unit in which the neutrino interaction occurred. A brick is composed of 57 emulsion films interleaved with 56, 1 mm thick, lead plates for a mass of 8.3 kg. Its thickness along the beam direction corresponds to about 10 radiation lengths and its transverse size is 128 × 102 mm 2 . A film consists of two 44 µm layers deposited on each side of a 205 µm plastic base. Another key ingredient for the experiment are the Changeable Sheet (CS) doublets [43], attached to the downstream face of each brick. This is a pair of films having received in the underground laboratory a special treatment (refreshing) aiming at erasing most of the cosmic background accumulated since their fabrication. Their scanning allows a relatively fast feedback on the prediction of the electronic detectors (ED) and provides a prediction of the event position in the brick at the O(10)µm level, thus greatly helping the vertex location. Finally a magnetic spectrometer system instrumented with Resistive Plate Chamber (RPC) detectors and high-precision Drift Tubes (DT), is used for the task of identifying muons and measuring their charge and momentum. A good muon identification capability is essential to reduce the background to τ decays from charmed particles produced in charged current (CC) ν µ interactions. The detector (figure 1) is divided into two identical units called Super Modules (SM), each consisting of a target and a spectrometer section. The average number of bricks has been about 140000 for a target mass of about 1.2 kt. OPERA was exposed to the CNGS ν µ beam [34][35][36] at a long-baseline, 730 km away from the source. The neutrino beam, produced by 400 GeV-protons accelerated in the SPS, has an average energy of about 17 GeV, optimised for the observation of ν τ CC interactions JHEP11(2013)036 in the OPERA detector. In terms of interactions, theν µ contamination is 2.1%, the ν e andν e contaminations are together below 1%, while the intrinsic ν τ component (from D s decays in the CNGS target and beam-dump) is of O(10 −6 ), hence negligible.

Data samples and event selection
CNGS completed its operation on December 3, 2012. A sample corresponding to 17.97 × 10 19 pot has been registered by the detector since the beginning of the program in 2008. The first task of the electronic detectors DAQ [44] is the time-tagging of the hits allowing the selection of events in coincidence with the two 10.5 µs-wide CNGS spills separated by 50 ms (the so-called "on-time" events).
An on-line filtering (TT local trigger) is then applied in order to remove background from random noise in the detector: hits are required in the horizontal and vertical views of at least two planes or the presence of at least 4 hits in a single plane is required with the sum of their photomultiplier ADC signals exceeding 500 counts (corresponding to about ten photo-electrons). 1 The event is then classified as being CC-like (hereafter "1µ") or neutral current (NC) like ("0µ") using the data of the target tracker and the spectrometers. Recorded hits are processed by a pattern recognition algorithm and sub-samples of hits in both views are grouped into three dimensional (3D) tracks. A 3D-track is tagged as a muon if the product of its length and the density along its path is larger than 660 g/cm 2 . An event is classified as 1µ if either it contains at least one 3D-track tagged as a muon [45] or the total number of TT and RPC planes having at least one hit is larger than 19. The complementary sample is defined as 0µ. About 19% of the NC events are classified as 1µ while only 6% of CC events are classified as 0µ. The momentum of 3D-tracks is calculated from their bending in the spectrometer magnetic field and/or their range with a Kalman filter-based reconstruction algorithm.
A pure sample of 106422 CNGS on-time neutrino interactions is selected (table 1, 3 rd column). About 60% of this sample results from neutrino interactions in the rock in front of the detector typically producing long passing-through muon tracks, while the rest is from interactions occurring, in about equal fractions, in the target (contained events) and in the spectrometer. A classifier algorithm, OpCarac [46], is applied to select the contained events yielding a total of 19505 interactions (table 1, 4 th column). This number corresponds to an average rate of about 18 contained neutrino interactions per day which have been stably recorded with an overall dead-time of the data acquisition and detectors of 1.5%.
Contained events are processed by a brick-finding algorithm [47]. The topology and the energy deposition in the TT scintillator strips, as well as the muon track information (when available) are used to define a three-dimensional probability density map for the vertex position. This probability is integrated over the volumes of the bricks and these are ranked in order of decreasing probability, for extraction and analysis. In the following, the highest-probability brick will be denoted as HPB.

Brick pre-selection and choice of the analysed sample
The analysis of the 2008-2009 inclusive sample (30% of the overall number of pot) was reported in [40]. In that sample of 2738 fully analysed events one event was recognised as a ν τ candidate decaying to a single charged hadron (noted as τ → 1h) [39]. More precisely the measurement of the final state topology and kinematics strongly favors this interpretation: τ − → ν τ ρ − with ρ − → π − π 0 and π 0 → γγ.
Since then the scanning strategy has been modified to accelerate the finding of a significant signal for ν µ → ν τ oscillation. Priority was given to the scanning of a preselected sample of bricks: • 2008-2009 data: analysis of up to two bricks per event without event pre-selection.
• 2010 data: analysis of the HPBs, for all 0µ events and for 1µ events with a muon momentum p µ < 15 GeV/c.
• 2011 data: analysis of the HPBs for all 0µ events.
This pre-selection is temporary and the analysis of lower-priority bricks will occur in the near future. 0µ events form a preferred sample to search for τ decays since these essentially contain signals of the electron final state, with a branching ratio of (17.85 ± 0.05)%, and of the 1-prong, (49.52 ± 0.07)% and 3-prong, (15.19 ± 0.08)% hadronic channels.
In the 1µ-decay channel, with a branching ratio of (17.36 ± 0.05)%, most events occur at low muon momentum. This is due to the fact that the atmospheric-scale oscillation, at this baseline, mainly affects the low-energy region. Moreover, a large fraction of the τ momentum is transferred to the two final-state neutrinos. This is illustrated in the left plot of figure 2 which also shows the distributions of the reconstructed signed muonmomentum (q × p µ ) for 1µ-events for data and Monte Carlo (MC). The shape of the data is well described by MC, indicating that the pre-selection of events is well understood. The momentum cut at 15 GeV/c (dashed lines) results in a loss in the sample of muonic τ decays of only 4%. Instead, this cut reduces by 33% the amount of ν CC µ events to be analysed and, more important, by 28% the size of the charm sample. The momentum cut The right plot in figure 2 shows the distributions of the brick probability for 0µ and 1µ events of the 2010 and 2011 samples, separately for the first two bricks in the probability ranking. Again, the shape of the brick probability is well described by MC. Details on the MC simulation will be given in section 4.
The number of events selected at the level of electronic detectors (N ED ) is 11149 (listed for each year in the first row of table 2). Part of this sample is discarded by the brick finding algorithm based on the topology of the measured hits in the TT. The rejected sample is mainly accounted for by punch-through of external interactions. The numbers of events with a good brick prediction, which are sent to the scanning laboratories (N pred , 9267 in total), are given in the 2 nd row in table 2.
Selected bricks are routinely extracted by a robotic Brick Manipulator System (BMS) capable of keeping up with the average weekly rate of neutrino interactions and thus allowing for the emulsion-based analysis to potentially proceed in parallel to the neutrino interactions data taking. Removed bricks with a negative CS result are re-inserted in the target after replacement of the CS while dismantled bricks are not replaced. Empty spaces are filled by a rearrangement of bricks. This approach aims at keeping the target homogeneous and minimise the occurrence of events with irregularities in the energy flow containment. The position of each brick is registered in a dynamic database to enable the event location.   Table 2. Summary of the analysed data samples. 2B stands for two-brick analysis while 1B denotes the highest-probability brick analysis. N ED (section 3.1) are the pre-selected events, N pred (section 3.1) are the events with a brick prediction, N brick (section 3.3) the events located in the bricks and N DS (section 3.4) the number of events for which the decay-search procedure was applied.
the DAQ system, red curve). The loss in target mass at the end of the running amounts to 8.5% of its maximum value. The pre-selected events (N ED ) represent a fraction of 78% and 33% of the 2010 and 2011 contained-event samples, respectively.

Analysis of the changeable sheets
After the extraction of the selected brick from the detector, the scanning of the CS doublets is performed in order to validate or disprove the brick-finding result. The pattern of tracks reconstructed in the CS doublets can confirm the prediction of the electronic detector, or act as veto and thus trigger the extraction of neighbouring bricks. CS doublets are inspected by automatic optical scanning microscopes [49][50][51][52][53][54][55][56][57][58] in two specialised scanning stations at the LNGS laboratory and at Nagoya University. A rectangular shape is defined centered on the prediction of the electronic detectors with a resulting JHEP11(2013)036 average scanning area of 20 cm 2 for 1µ events and 35 cm 2 for 0µ events, the pointing accuracy here being poorer due to the absence of a muon track.
The tracking efficiency of the CS doublets has been measured in neutrino data and test beam exposures [43,59]. A charged particle crossing the CS will produce 4 track segments called micro-tracks. The micro-track detection efficiency for minimum ionising tracks is 95% [59] with the automatic scanning system. The distributions of the residuals in position and angle between the ED and CS tracks have a standard deviation of about 8 mm and 15 mrad respectively.
The analysis and scanning strategy have been improved during the years. In the analysis of the 2008 and 2009 runs, a brick was sent to development as soon as at least one CS track had a possible match with some ED hits. Subsequently the strategy ("CS trigger") has been refined for the 2010 and 2011 runs by requiring the presence of at least: • for 1µ events, a CS track compatible with the ED muon track within 60 mrad • for 0µ events, a CS track matching an isolated ED track • 2 or more CS tracks possibly converging towards a common origin in the brick.
If none of these conditions are fulfilled, the brick is put back in the target with a new CS doublet and the next brick in the probability map is extracted (in case of multi-brick analysis). This method allows saving scanning and analysis time and minimising the target mass loss. In case of a positive outcome, the brick is exposed to cosmic rays in a dedicated pit in the surface LNGS laboratory for 14 hours for high-precision film-to-film alignment, later to be dismantled in the dark room where the emulsion films are developed. The films are finally dispatched to the scanning laboratories of the Collaboration for the "vertex location" and "decay search" analysis.

Vertex location in the brick
CS tracks are projected to the most downstream emulsion film in the brick (through a distance of about 4.5 mm mainly filled with the plastic and Aluminium film of the CS and brick boxes) where they are searched for. The residuals in position are at the level of 50-60 µm such that the tracks are typically found within the predicted microscope view (400 × 300 µm). They are then followed upstream film by film ("scan-back" [60,61]) adjusting the predictions in angle and position at each step to cope with multiple Coulomb scattering (MCS) in the lead plates and measurement errors effects.
The scan-back procedure is stopped when no track candidate is found in three consecutive films; the lead plate just upstream of the last detected segment is defined as the candidate vertex plate. At this stage a volume is defined with a transverse area of 1 cm 2 for 5 films upstream and 10 films downstream of the stopping point (or less if the stopping point lies too upstream or downstream in the brick) and tracks within an angular acceptance tan θ < 0.6 (θ being the angle of the track with the z axis, figure 1) are searched for in this volume (general scan).
The scan-back procedure is modified for events in which the CS-trigger is produced only by a shower-like topology (mainly due to γ conversions from π 0 decays). In such cases JHEP11(2013)036 an area of 1 cm 2 centered on the shower axis is analysed for 20 films (corresponding to about 3.5 X 0 ) starting from the last plate.
All track segments collected in the scanned volumes are analysed by offline algorithms which perform precise alignment, tracking and vertexing. The alignment of the films (with an accuracy of a few µm) is assured by the cosmic ray tracks (section 3.2).
In some cases the scan-back track is not associated to any other track inside the volume due to a low-multiplicity event or because of the angular acceptance limitation or tracking inefficiencies. In this condition the vertex is consider as detected. A further signature of the neutrino interaction can then be obtained from the analysis of nuclear fragments which might be visible in the forward and/or backward hemisphere especially if the interaction happened close to the emulsion layers. Dedicated image analysis tools have been developed [62,63] to detect highly ionising large-angle tracks in the proximity of the stopping point.
Interactions are tagged as "dead material" when the extrapolated vertex does not lie in the lead plates or in the films (i.e. it is in the scintillator strips, in the brick supporting structure) and no event-related tracks are found in the upstream-brick CS-doublet.
For the 2008-2011 analysed sample, the events located in the brick lead plates (N brick ) amount to 4135 (table 2).

Data-MC comparison for the location efficiency
An important check of the understanding of the complex task of event location is done at this stage. We consider the data collected in the year 2008 and 2009 and the two-brick analysis strategy.
For the definition of the initial control sample of 0µ events (N 0µ,ED ) we adopt a stricter selection on the fiducial volume and a minimal requirement on the energy released in the TT (E T T > 70 MeV). This allows getting the contamination of external events to the level of 5% in the 0µ sample. This residual contamination has been statistically subtracted. The location is then defined, f DM being the fraction of interactions in the dead material which is estimated from the MC simulation (f DM = 7.8%) and f BQ being the fraction of films which could not be analysed due to their bad quality (f BQ = 6%).
The dependence of the location efficiency loc on E T T is shown in figure 4 for data and MC for 1µ and for 0µ events. The two quantities are highly correlated in the 0µ sample where the hadronic activity plays a crucial role in the location. The bands in the MC prediction are representing the systematic uncertainty which is 10% for the 1µ sample and the 0µ sample with E T T > 200 MeV and 20% for the 0µ sample with E T T < 200 MeV. The systematic uncertainties account for residual differences in the MC description of the data: the implementation of the CS trigger and brick scanning strategy, the definition of the dead material contribution. For the 0µ sample an additional uncertainty arises from the subtraction of the external background component.
The agreement between data and MC is good in shape and normalisation for the 1µ sample while for the 0µ sample the agreement is good in shape, but data tends to be below the MC by about 15%, independently of E T T . It must be noted, however, that this  residual difference in the location efficiency in data and Monte Carlo for the 0µ sample has no effect on the predicted signal or background events which is normalised to the number of localised events in the data (see section 4.1).

Decay search (topological selection)
The decay search (DS) procedure is aimed at detecting the decay topologies of τ leptons produced in ν CC τ interactions once a vertex has been identified in the volume scan data. The decay is defined as "short" if it happens in the same lead plate where the neutrino interaction occurred or in the first downstream emulsion layer and as "long" if it happens further downstream such that at least one complete micro-track is produced by the τlepton. About 46% of the τ decays are expected to be short (43% in the first lead plate and 3% in the first emulsion layer) while the remaining 54% are long decays happening in the first plastic base (11%), in the second emulsion layer (2%), in the second lead plate (25%) or further downstream (16%).
Candidate daughter tracks from short-lived particles decays are selected by requiring their impact parameter (IP ) with respect to the reconstructed neutrino interaction vertex to be larger than 10 µm if the depth in lead (λ) is lower than 500 µm or by loosening this requirement to IP > (5 + 0.01 × λ) µm for deeper vertices. The next requirement is that the momentum of the candidate daughter track measured from its MCS (p mcs ) [64] is larger than 1 GeV/c. If the number of planes in the scanned volume is not sufficient to estimate p mcs the angular spread S θ of the available segments is evaluated (see [65] for a detailed definition). If S θ > 15 mrad in both views, the track is discarded, otherwise the MCS measurement is extended to more plates.
A search is then performed for additional tracks (called extra-tracks) which are not directly connected to the primary vertex: they must be detected in at least 3 films and the distance along the z-axis between the most upstream segment and the neutrino vertex JHEP11(2013)036 (∆z = z up − z vtx ) is required to be positive and less than 3.6 mm. In addition we require IP < 300 µm if ∆z < 1 mm and IP < 500 µm for ∆z > 1 mm.
In order to deal with topologies in which the first-reconstructed vertex is the decay vertex and thus downstream of the primary neutrino vertex (e.g. in a τ → 3h decay) the search for extra tracks is also extended upstream, requiring ∆z > −2 mm and IP < 500 µm. For the events located on the basis of a single track (which could be due to a highermultiplicity decay with tracks outside of the scanning angular acceptance), extra-tracks are searched for within a ∆z of ± 3 mm and IP < 500 µm. In this case the z position of the vertex is assumed to be at the center of the lead plate.
A search for charged short-lived parents is then applied to both upstream and downstream extra-tracks by looking for tracks or single segments with IP < 10 µm and a distance of closest approach d ca < 20 µm with respect to the possible daughter tracks.
Finally the presence of a significant kink (larger than 20 mrad) is checked on the muon track for 1µ events and on all tracks for 0µ events. The four most upstream segments of the track are used to evaluate the ratio (R) between the maximum angular difference between pairs of segments and the overall angular spread [65]. If R > 5 then a "kink trigger" is issued and further analysis is performed.
All these operations are supplemented by manual eye-inspection at the microscope in order to improve the information provided by automatic scanning in terms of angular resolution (the automatic system accuracy can sometimes be spoiled by the association of uncorrelated background grains) and efficiency (by recovering the information for inefficient planes). Eye-inspection is also employed to improve the purity by discarding fake tracks (most important for large angles and high-fog 2 films), and furthermore allows excluding passing-through cosmic ray tracks which can mimic extra-tracks due to inefficiencies of the automatic scanning in the five upstream films.
The number of events with a completed decay-search (N DS ) for the present analysis amounts to 3969 (table 2). This number differs from N brick by about 4% due to the events which fall at the edges of the brick (either in the longitudinal plane or in the transverse direction).

The charm control sample
Given the similarity in mass and decay topologies, the detection of charmed particles constitutes not only a background but also an important tool to verify the understanding of the τ detection efficiency up to the topological selection (decay-search) level.
The charm sample is selected using the same DS analysis as used for the signal search with the exception that the kinematic selection of ν τ candidate events (which will be described in the following section) is not applied. 3 The charm yield has been predicted using the latest analysis of the CHORUS data [66] at the WANF neutrino beam ( E ν 27 GeV) and re-weighting it for the different neutrino 2 AgBr crystals can be activated in the emulsion due to thermal or mechanical excitations resulting in the presence of randomly distributed grains which are referred to as "fog".  figure 5 for data and MC. Not only the absolute yields but also the shapes of the distributions are in very good agreement, which indicates that the systematic error on the estimated efficiency of the full analysis chain cannot exceed 20%. A more extensive discussion of the charm sample will be presented in [67].

Kinematic selection
Several kinematic quantities of the neutrino interaction are accessible at brick-level via momentum reconstruction using MCS. The energy of electrons and photons is also measured by employing calorimetric techniques [68,69]. Kinematic criteria can be defined to improve the signal-to-background ratio. To improve the acceptance for electromagnetic showers and reduce the error on the track momentum measurement the standard volume considered for the location (section 3.3) is enlarged and tracks are followed downstream, eventually in other bricks (section 3.6).
After denoting the charged tracks emerging from the neutrino interaction vertex or the decay vertex as "primaries" and "daughters" respectively, and the short-lived primary track as "parent", we define the following relevant variables:

JHEP11(2013)036
• z dec : the z-coordinate of the decay vertex with respect to the downstream face of the lead plate containing the primary vertex (z dec < 44 µm for short decays).
• p 2ry T : the transverse momentum of the daughter with respect to the parent direction, for 1-prong decays.
• p miss T : the magnitude of the vectorial sum of the transverse momenta of primaries (except the parent) and daughters with respect to the neutrino beam direction.
• p 2ry : the scalar sum of the momenta of the daughters.
• θ kink : the average 3D angle between the parent and its daughters (kink angle).
• m: the invariant mass of the daughters (calculated attributing the π mass).
• φ lH : the angle between the parent and the vectorial sum of the other primaries calculated in the plane perpendicular to the CNGS axis.
In the calculation of φ lH , if the primary multiplicity (including the τ track candidate) is larger than two, the primaries are defined after removal of the track with the largest difference in φ with respect to the τ candidate (TLD) unless it is identified as a hadron with high probability by the track follow-down procedure (section 3.6). This is intended to remove events where the TLD is actually a genuine muon which has not been identified as such by the electronic detectors (for these events, if the TLD would not be excluded the angle would be large and thus the kinematic cut on this variable would not be effective). The applied selections are summarised for each channel in table 3. The choice of cut parameters has been studied and defined in [23,24] for the 1h, electron and muon channels and in [71] for the 3h channel.

Track follow-down
The track kinematics of interesting events fulfilling the selection are further studied in neighbouring bricks. Primary tracks are followed (with an angular acceptance extending up to tan θ = 1) until either a stopping point, an interaction or a muon decay topology is found (track follow-down, TFD). Thus, by the study of momentum-range correlations, track length, energy loss in proximity of the stopping point, and (if possible) the tagging of interactions or muon decays, a muon/hadron separation exceeding the limitations of the electronic detector reconstruction is obtained.
This procedure results in a significant reduction of the background in all channels from ν CC µ charm production where the primary µ is not identified in the ED. The TFD technique also provides an important reduction of the background in the τ → µ channel from ν CC µ events in which the primary µ is wrongly associated to a secondary vertex with a kink topology from a hadronic re-interaction. This also holds for ν N C µ events with a fake muon. The improved muon/hadron separation provided by TFD is also effective in reducing the

JHEP11(2013)036 4 Monte Carlo simulation of signal and backgrounds
Compared to the results published in [40] several improvements have been adopted in the simulation of the detector. Efficiencies and backgrounds are based on a new software framework which combines the information of the electronic and ECC detectors. All particle trajectories are digitised at the micro-track level in a volume consisting of 3×3×3 bricks, with the brick containing the neutrino interaction at the center. For dedicated studies (i.e. for the TFD simulation), also larger volumes have been considered. The efficiency and the resolution (in angle and position) of the scanning microscopes are simulated using parametrisations obtained from real data. The framework allows reproducing the analysis flow from the prioritisation based on the brick probability, the CS-trigger, the CS-to-brick connection, the scan-back, the alignment and the vertexing with full access to the reconstruction parameters. Using this tool the decay search procedure has been included in the efficiency evaluation. Furthermore the detection efficiencies have been re-evaluated for all channels and the simulation of the brick extraction strategies for the different years of data taking has also been updated.

Expected signal event rates
The neutrino fluxes used in the calculation of the expected signal rates are based on a FLUKA [72,73] simulation of the CNGS beam-line (2005 revision [74]). The neutrino interactions in the detector are generated using the NEGN generator after a tuning of the parameters based on the high-statistics data sample of the NOMAD experiment [75] . The energy-dependence of the ν τ cross section that has been used is the default implementantion contained in the GENIE v2.6 simulation program [76]. The ν µ → ν τ oscillation probability was evaluated using ∆m 2 23 = 2.32 × 10 −3 eV 2 [77] and sin 2 2θ 23 = 1. With these ingredients, the expected rate of ν CC τ on lead at true level is equal to 3.32 events/(10 19 pot kt). Taking into account the time evolution of the pot delivery rate and the effective mass of the detector throughout the running period (figure 3) a total number of 66.4 ν τ CC interactions in lead are expected (with the non-DIS component being about one third [76] of the total).
The expected rate of ν CC τ events in the 0µ sample can be expressed in terms of the measured event rate in the same category (n 0µ from both ν CC µ and ν N C µ ) as: denoting with r (k) the efficiency for reconstructing the process k (ν CC µ , ν N C µ ) in the reconstructed category r (0µ, 1µ) after convolution with the CNGS flux φ; N (ν CC µ ) being the ν CC µ interactions at true level and α = σ(ν N C µ )/σ(ν CC µ ) ⊗ φ. In the present analysis the actual number of n 0µ is taken from the data while other quantities are MC driven. The expected numbers obtained with the above-mentioned procedure are insensitive to systematic effects on the efficiencies up to the location level being common to ν τ and ν µ events. The same approach is followed in the prediction of the ν τ component in the 1µ sample. In this particular case, in addition to the standard 1µ classification (section 3),

JHEP11(2013)036
the presence of a muon tagged 3D track is also required due to the need of matching this ED track to the ECC secondary µ candidate track for the sake of background reduction.
Finally it must be noted that in the present approach the signal efficiency is not taking into account possible migrations between different channels and could therefore be slightly underestimated.

Expected background event rates
Three sources of backgrounds are giving significant contributions to the final sample: charmed particles decays, hadronic interactions and large-angle muon scattering (LAS). The last one only affects the τ → µ channel while the other classes differently affect all decay channels. These backgrounds have been discussed extensively in [23,24]. The basic ingredients for the more recent evaluation of these components will be discussed in the following.

Charmed particle decays
The most effective tool for the rejection of the background from charmed particles is an efficient identification of the primary muon in ν CC µ interactions. The fraction of ν CC µ events having an associated charm quark at the CNGS energies is estimated from the CHORUS measurement [66] as (σ(ν CC µ + c)/σ(ν CC µ )) ⊗ φ = (4.38 ± 0.26)%. The tagging of the primary muon is achieved first, at the level of the ED, via the reconstruction of penetrating 3D-tracks and the vetoing of events involving many TT planes (to complement the tracking) as described in section 3. Requiring an event to be classified as 0µ reduces the yield down to 6% of the initial charm sample. At the emulsion detector level, the TFD procedure gives a further suppression by 60%. Taking into account the effects of the full analysis chain and considering only the HPB, an overall suppression of 3.3 × 10 −5 is obtained in the 1h channel and 3.4 × 10 −4 in the 3h channel. The two-brick analysis gives an increase by a factor 20-25% resulting in an expectation of 0.027 ± 0.005 events and 0.12 ± 0.02 events in the analysed sample for the 1h and 3h channels respectively (figure 7, green-coloured histogram and table 4).
The background from charmed particles in the τ → e decay channel (where, unlike in the muon case, the positive charge of the decay electron cannot be measured) gives a contribution in the analysed sample of 0.020 ± 0.004 events. This is the only significant background in this channel.
The background from charmed particles in the τ → µ decay channel is relatively small since it only arises in two cases i.e. 1) when the primary µ is not identified, a muonic decay occurs (f (C + → µ + h 0 ) = (5.3 ± 2.8)%), and the positive charge of the secondary µ is not measured or mis-identified; 2) the secondary µ is not identified and the primary µ is wrongly matched to the decay daughter. This results in a background contribution of 0.011 ± 0.005 events ( figure 7, green-coloured histogram and table 4).
Finally, the double charm production in NC interactions process has to be considered. The yield has been measured by the CHORUS collaboration to be σ(ccν)/σ DIS N C = (3.62 +2.9 −2.4 (stat.) ± 0.54(sys.)) × 10 −3 [78]. This component cannot profit of the suppression given by the primary muon identification but, on the other hand, it can be removed by JHEP11(2013)036 measuring the presence of the primary and secondary vertices. This contribution, which is expected to be a second order correction, is not included in the present estimates.
The uncertainty on the charm background component is estimated within the 20% level based on the measured sample of charm events (section 3.4.1).

Hadronic re-interactions.
This background has been estimated by giving hadronic tracks from located 0µ events (from an initial sample of 9×10 6 ν N C µ interactions) as an input to a FLUKA-based MC simulation. With respect to a similar analysis reported in [40] this approach implements the FLUKA hadronic models to fully reconstructed events thus including the biases introduced by the selection chain, up to the event-location level.
The fraction of located 0µ events with at least one track mimicking the topology of a single-prong long decay and θ kink > 20 mrad is found to be 2.1% while the fraction with at least a track producing three visible prongs in the scanning acceptance is 0.18% (short and long). The kinematic selection (in particular the p 2ry and p 2ry T cuts) strongly reduces the "topological" sample to (0.27 ± 0.02(stat.))% (1h) and (1.2 ± 0.1(stat.))% (3h), due to the typical low transverse momentum of the secondary products. At this stage an additional reduction of 30% is obtained by requiring the absence of nuclear fragments (either in the backward or forward hemisphere) or minimum ionising particles up to tan θ = 3. Finally the estimated rate of background events amounts to (3.9 ± 0.2(stat.))×10 −5 of the located 0µ events for the 1h channel and (1.5 ± 0.2(stat.))×10 −5 for the 3h channel. The expected interaction rate per unit length in the plastic base is about 4.5 times smaller than in the lead plate due to the combined effect of the density and the mean atomic number.
Several data-driven checks of the FLUKA description of hadronic interactions in the OPERA bricks were performed [40]. In general a good agreement between data and simulation is observed for different data-sets: a sample of hadronic interactions of pions of 2, 4 and 10 GeV/c momenta from a CERN-based test beam; a sample of hadronic tracks from CNGS neutrino interactions measured in the emulsions (total length 19 m) and a sample of hadronic nuclear fragments from test-beam pion interactions for which the yields and angular distributions have been studied. With respect to [40], analyses based on larger data samples have been developed [79][80][81][82] allowing estimating an accuracy of the predictions for the hadronic background at the 30% level. The uncertainty in the rate of high-angle nuclear fragments emitted in hadronic re-interactions is estimated in [82] at the 10% level. In conclusion, considering the 1h and 3h channels together, in the analysed sample a total of 0.018 ± 0.005 events from hadronic re-interactions is expected (figure 7, yellow-coloured histogram and table 4).

Large-angle muon scattering.
The occurrence of large-angle scattering of muons in thin (O(0.1)X 0 ) lead plates is, at present, not well constrained by measurements. Upper limits from extrapolations of measurements on copper or nuclear emulsions have been reported in [23,24]: the rates in the signal region θ kink > 20 mrad and p µ T > 0.25 GeV/c for muon tracks with a realistic angular and momentum spectrum, are, at 90% C.L., < 2.3 × 10 −5 and < 4.1 × 10 −5 , respectively.   Table 4. Signal and backgrounds expectations for the analysed sample. The numbers correspond to the filled histograms in figure 7.
A dedicated measurement of large angle scattering of 9 GeV/c muons from thin lead plates is also reported in [23,24] giving 0.6 +0.7 −0.6 × 10 −5 . A GEANT 3.21 [83] based simulation, modified to take into account data on the nuclear form factor of lead and inelastic interactions, predicts a rate of 2 × 10 −6 . More experimental activities to determine this process are in progress. In the present work the same contribution used for the experiment proposal [23,24] (1 × 10 −5 ) is assumed corresponding to 0.009 ± 0.005 events in the analysed sample (figure 7, blue-coloured histogram and table 4).

Summary of the expected event numbers.
The estimated signal and background events for the sample considered in the present analysis are summarised in a graphical form in figure 7, for each decay channel and in tabular form in table 4. The background from charmed particles is the dominant one in the τ → e and τ → 3h channels while LAS in the τ → µ channel and hadronic re-interactions in the τ → 1h channel are the largest backgrounds. The systematic uncertainty on the signal and on the charm background is estimated to be 20% (section 3.4.1) and the one on JHEP11(2013)036 the hadronic background to be 30%. The LAS background uncertainty is assumed to be of the order of 50% (grey bars in figure 7).
The expectation for the signal events using the two-brick analysis in the final data sample (2008-2012) is of 3.2±0.3 observed ν τ events (dashed line in figure 7, left). This number is calculated using ∆m 2 23 = 2.32 × 10 −3 eV 2 as a central value, 4 the accumulated statistics of 17.97 × 10 19 pot and a dynamic target mass (figure 3) corrected for the dead material and the fraction of bad-quality films (section 3.3). Furthermore this estimate is based on a realistic simulation of all the steps of the analysis chain and in particular the decay search phase (section 3.4) which had not been simulated in full detail in previous analyses.
The current estimate does not include some factors which are expected to increase the efficiency. The extension of the analysis to the 3 rd and 4 th bricks in the probability map is foreseen. Smaller effects due to the migration of ν τ events from one channel to another are not taken into account at present. These effects increase only the expected number of ν τ events, while keeping the background at the same level. Finally an optimisation of the selection is ongoing within the new simulation framework with the goal of maximising the expected sensitivity.

Analysed sample and results
The number of events from the present analysed sample surviving the complete selection chain described above amounts to two τ -candidates (hereafter called A and B). Given that the candidate which occurred first (A, τ → 1h) has already been extensively documented in [39,40], this paper will focus on the description of the second candidate.

Description of the second ν τ candidate event
This neutrino interaction occurred on 23 April 2011 at 7.15 UTC time. The pattern of hit scintillator strips in the TT is shown in figure 8. The event, of which the estimated hadronic energy is (22.0 ± 6.2) GeV, is classified as 0µ and the bulk of the activity in the electronic detector is contained within about 6 to 8 brick walls (more than 60 X 0 and about 2.5 pion interaction lengths). The interaction took place in the target of the upstream Super Module and lies well within the brick-filled target region. The neutrino vertex brick (having a probability of 63%) was located in the longitudinal direction (z), in the second most upstream brick layer (called W 0 ), in the horizontal direction (x), in the 3 rd brick layer from the closest side (left side looking towards the CNGS) and, in the vertical direction (y), in the 19 th brick layer from the bottom. In figure 8  z order) while in the transverse plane it lies at 6.8 mm and 41.4 mm from the closest sides in x and y.
The topology of the primary vertex (V 0 in figure 9) consists of two tracks, the τ lepton candidate and another track (called p 0 ). The distance of closest approach of the p 0 and τ tracks equals to 0.2 µm and the vertex lies close to the downstream emulsion film, at a depth in lead of only 120 µm. A forward-going nuclear fragment associated to the primary vertex has also been detected at a large angle with slopes of: (1.15, −0.28). 5 The flight length of the τ lepton candidate is (1466 ± 10) µm and its decay occurs in the plastic base allowing excluding with a high efficiency (above 99.8% at 90% C.L. up to tan θ = 3, [62]) the emission of highly ionising nuclear fragments. The secondary vertex (V 1 , figure 9) consists of three tracks (called d 1 , d 2 , d 3 ). A display of the reconstructed grains in the emulsion layers is presented in figure 10. The impact parameters of the decay products with respect to the reconstructed secondary vertex are 1.3, 1.2 and 0.3 µm for d 1 , d 2 and d 3 respectively. After eye-inspection the background from instrumental fake tracks or tracks due to Compton electrons is negligible.
In order to strongly constrain the hypothesis that the secondary vertex could be a hadronic interaction, a search for nuclear fragments has been performed both upstream and downstream of the vertex with an automatic scanning up to tan θ = 3.5 as well as by visual inspection. No such fragment was found.
A search for γ conversions has been performed up to tan θ = 1 for the 35 films (about 6 X 0 ) downstream of the vertex yielding no candidates.  The analysis of each track to determine its nature and momentum is described below: • track p 0 has a measured momentum p p 0 = (2.8 +0.7 −0.7 ) GeV/c. This track points towards the centre of the detector in the top-view such that its signature as a hadron is already well constrained using the target tracker information only. Nevertheless the track has been followed in the downstream wall (W 1 ) where it is found to exit the brick to the side. Compatible tracks have neither been found in the adjacent brick in W 1 nor in downstream bricks (and CS doublets) up to W 3 . Gamma rays from a possible hadronic interaction were also searched for in two bricks in W 2 , however, no γ was found. This track is then assumed to interact in the dead material in-between two bricks. A muon with a momentum of the magnitude measured by MCS would be expected to travel from 26 to 44 brick layers before stopping. This makes the muon hypothesis very unlikely: D TFD is 0.05 for this track. Its projected slopes are (0.155, −0.365), well inside the angular acceptance of the scanning. Even considering the possibility of having missed a muon track crossing the downstream emulsion detectors, the pattern in the scintillators does not allow the existence of such a track for more than about 7-8 brick walls.
• track d 1 has slopes of (-0.056, 0.101) and a measured momentum p d 1 = (6.6 +2.0 −1.4 ) GeV/c. A hadronic interaction is detected in the emulsions of the brick in wall W 4 (see figure 8) producing two charged tracks with slopes of (0.234, 0.489) and (0.034, -0.305). The signature of the interaction is also indicated by the target tracker scintillators.
• track d 2 has a slope of (-0.041, 0.260) and a measured momentum p d 2 = (1.3 +0.2 −0.2 ) GeV/c. It has not been found in W 2 or in the following walls corresponding to a rangemomentum correlation parameter D TFD = 0.25.
• track d 3 has a slope of (-0.134, 0.220) and a measured momentum p d 3 = (2.0 +0.9 −0.6 ) GeV/c. This track interacts in the brick containing the neutrino vertex, after 11 lead plates i.e. about 1.3 cm downstream (V 2 in figure 9). The interaction occurs inside the emulsion resulting in a very clear signature. The final state is composed of two charged tracks (d 3−1 and d 3−2 ) and four back-scattered nuclear fragments.
The scalar sum of the momenta of all the measured charged particles in the event is 12.7 +2.3 −1.7 GeV/c.  Table 5. Selection criteria for ν τ candidate events and corresponding measured values.

Summary of the two ν τ candidate events
From the arguments developed in [39,40] and in section 5.1, it is concluded that the two events are candidates of τ decays into the 1-charged hadron (1h) and 3-charged hadrons (3h) channels, respectively. table 5 summarises the values taken by the kinematic variables for the two candidates together with the cuts applied in the analysis. The distributions of φ lH , for the expected signal and background components in the analysed sample are presented in figure 12 for the 1h and 3h channels separately. The measured values in the data are indicated by vertical lines.

Significance of the observation
The significance of the observation of two candidates with the present estimate of the backgrounds is addressed by considering the confidence in the exclusion of the null hypothesis (i.e. having observed no ν µ → ν τ oscillation signal). The probability that the two events JHEP11(2013)036 might be due to an upward fluctuation of the background is defined as the integral of the Poisson distribution evaluated for n ≥ 2: p = ∞ n=2 µ n e −µ n! . Using for µ the sum of the background in the four channels (see figure 7), the value p = 1.36 × 10 −2 is obtained, corresponding to a 2.20σ exclusion (adopting the one-sided definition).
In a similar manner the p-values of the two channels can be defined: p 1h = 2.6 × 10 −2 (1.94σ ) and p 3h = 1.10×10 −1 (1.23σ ). To combine them, an approach based on generating a large number of pseudo-experiments has then been followed. For each of the four τ decay channels an integer, n i=1,...,4 , is extracted from the Poisson distribution of background. The p i values are then calculated as above for each extraction and the estimator of the results p is obtained by taking their product p = p 1 p 2 p 3 p 4 (different choices of estimator exist). The counting of the fraction of pseudo experiments for which p ≤ p 1h p 3h yields a significance of 2.40σ. More elaborated estimates of the significance could be obtained using a log-likelihood ratio analysis or other multi-variate techniques. This approach will be addressed in future works.

Conclusions and prospects
The results of the ν τ appearance analysis on a pre-selected sample of the neutrino interactions collected by the OPERA experiment in the years 2008 to 2011 are reported. A two-brick analysis was performed on all 2008 and 2009 predictions while for 2010 and 2011 the analysis was restricted to the most probable brick of all 0µ events and of 1µ events with p µ < 15 GeV/c of the 2010 sample. In terms of expected signal events the considered sample corresponds to about 50% of the final expectation. A total of two ν τ candidates has been observed, one in the τ → 1h channel and one in the τ → 3h channel. This result is compatible with the expected amount of 1.53 ± 0.16 signal events and 0.18 ± 0.02 background events in all channels. The good agreement between data and MC, both for the location efficiencies of ν µ events and for the detection of charmed particles, indicates that the overall τ finding efficiencies are well understood. Using the presently analysed sample the absence of a signal from ν µ → ν τ oscillations (null hypothesis) is excluded at 2.40σ.
In the near future, the analysis of the 2011 sample will be completed using the same event selection as for 2010 as well as that of the 2012 sample currently in progress. Finally the search for events not found in the HPB will be extended to second-priority bricks as it was done for the 2008 and 2009 samples. The inclusion of 3 rd and 4 th bricks, which is also foreseen as a further step, will bring an additional increase of the location efficiency. The current significance can be improved with the ongoing increase of the analysed sample.