Probing lepton flavour violation via neutrinoless τ → 3μ decays with the ATLAS detector

This article presents the sensitivity of the ATLAS experiment to the lepton-ﬂavour-violating decays of τ → 3 μ . A method utilising the production of τ leptons via W → τν decays is used. This method is applied to the sample of 20.3 fb − 1 of pp collision data at a centre-of-mass energy of 8 TeV collected by the ATLAS experiment at the LHC in 2012. No event is observed passing the selection criteria, and the observed (expected) upper limit on the τ lepton branching fraction into three muons, Br (τ → 3 μ) , is 3 . 76 × 10 − 7 (3 . 94 × 10 − 7 ) at 90 % conﬁdence level.


Introduction
The observation of a lepton-flavour-violating (LFV) process involving charged leptons would be a major breakthrough in understanding the matter content of the universe and would support the hypothesis of leptogenesis [1]. In particular, LFV processes involving both a τ lepton and a muon are seen as most promising for such an observation, given the current measurements of neutrino oscillations [2]. In the Standard Model (SM), such processes have a vanishingly small branching fraction, e.g. Br(τ → 3µ) < 10 −14 [3], while a number of models beyond the SM predict it to be of the order of 10 −10 -10 −8 [4][5][6]. The current limits on branching fractions of neutrinoless τ lepton decays are of the order of few times 10 −8 [7][8][9][10], for Z boson LFV decays they are about 10 −5 [2,11,12], and for the LFV decay of a Higgs boson to a τ lepton and a muon they are about 1% [13,14]. The main experimental obstacles to improve the sensitivity with τ leptons are the small number of produced τ leptons world-wide.
In this article, a search for neutrinoless τ lepton decays to three muons is performed with 20.3 fb −1 of pp collision data collected with ATLAS detector in 2012 at 8 TeV centre-of-mass energy. The search is focused on a particular source of τ leptons, namely W → τν decays with subsequent τ → 3µ decay. In such events, τ leptons are produced with a transverse momentum (p T ) mostly in the range of ∼ 25 − 50 GeV. Due to the relativistic boost of the τ lepton, the muons from the τ LFV decay are produced in close geometrical proximity to each other but isolated from other energetic particles in the event. The tau-neutrino from the W boson decay appears as missing transverse momentum (E miss T ) in the detector and together with the transverse momentum of the three muons (p 3µ T ) gives a transverse mass, m T = 2p 3µ T E miss T (1 − cos ∆φ), compatible with the W boson decay, where ∆φ is the angle between the directions of the p 3µ T and the E miss T . The unique signature in the detector is three muons with invariant mass equal to the mass of the τ lepton and with a significant missing transverse momentum that is on average back-to-back with the three muons in the transverse plane. Since no energetic jet is expected in the majority of W boson production events, very small hadronic activity is predicted beyond that from the soft underlying event or multiple simultaneous pp collisions (pile-up). A large fraction of such τ leptons decay sufficiently far from the W production vertex to give a fully reconstructable additional vertex. This allows the selection of three muons originating from a vertex which is displaced from the primary interaction vertex. The background events usually contain one or two muons originating from the decay of hadrons, including decays in flight, while the remaining tracks are hadrons mis-measured as muons, originating from e.g. a pile-up jet or a pion punching through the calorimeter. The dominant background is due to muons originating from decays of b-or c-hadrons (heavy flavour, HF). Although such decays are typically accompanied by jets of particles produced in the direction opposite to the HF jets, in a fraction of the events the associated jet is lost or mis-measured, mimicking the signal E miss T . A small light-flavour multi-jet contribution is also present while the contribution from leptonic decays of vector bosons is negligible.
The analysis strategy is as follows. Events with three muons associated with a common vertex are selected. A loose event selection is applied to collect a high-quality sample of candidate events satisfying |m 3µ − m τ | 1 GeV. The characteristics of the loose sample of events are then analysed with a boosted decision tree (BDT). The BDT input variables are chosen so that the BDT output and the three-muon mass are uncorrelated in the mass range used in the analysis. A tight selection, following an initial cut on the BDT output, is applied to separate the signal from the background. After the optimal cut on the BDT output is found, a search is performed for an excess of events at the τ lepton mass above the expected background level.
The branching fraction is calculated as where N s is the number of observed events above the expected background level in a narrow region around the τ lepton mass, A s × s is the detector acceptance times efficiency for the signal, and N W→τν is the number of τ leptons produced via the W → τν channel (additional contributions to the τ lepton yield are estimated to be less than 3% [21] and the uncertainty in the 7 and 8 TeV luminosities. For the selection applied in the analysis, the contamination from other sources of τ leptons, such as Z → ττ or HF processes, is less than 3% and is therefore neglected. The background is estimated using data as discussed in Section 5.5.

Trigger and reconstruction
To maximise the signal acceptance times efficiency, events are required to pass at least one of seven triggers. These are six multi-muon triggers and one dimuon plus E miss T trigger. The software-based trigger thresholds used for the muons range from 4 to 18 GeV in transverse momentum while the E miss T threshold is 30 GeV. The trigger efficiency for simulated signal events within the muon-trigger acceptance (three generator-level muons with p T > 2.5 GeV and |η| < 2.4) is ∼ 31% for the combination of all triggers used in the analysis. To evaluate the trigger performance in the region where the muons have a small angular separation, as is typical for the signal, a tag-and-probe study is performed using data events containing high-momentum J/ψ → µµ candidates. For this study, the data are collected using a single-muon baseline trigger with a p T threshold of 18 GeV. Single-muon efficiencies are measured separately for the different thresholds which define the six multi-muon triggers. Each multi-muon trigger efficiency is calculated as the product of the single-muon efficiencies. Correction factors are applied to account for the limited performance of the trigger system in identifying a pair of muons as two muon-trigger objects. At small angular separations (∆R 0.2), where most of the signal is expected and where these limitations are most pronounced, these corrections must be taken into account. These factors are measured from the efficiency to identify two independent muon-trigger objects for different ∆R values between the tag-and the probe-muon. The total efficiency of the seven triggers is calculated considering correlations between any of the triggers. The trigger efficiency, measured from the data, is compared to the one measured in simulated J/ψ events for the seven different triggers separately and jointly. Agreement between data and MC simulation was found to be within 11% for all relevant values of ∆R and p T , where the largest difference comes from events where the ∆R separation is smallest. The systematic uncertainty on A s × s due to the trigger is therefore taken to be 11%.
The approach for measuring the muon reconstruction efficiency is similar to that used to measure the trigger efficiency. While the trigger efficiency is measured with respect to muon reconstruction as the baseline, the reconstruction efficiency is measured with respect to ID tracking, which in turn is close to 100% efficient [22]. Small deviations from the assumed value for ID tracking efficiency have a negligible impact on this measurement. The tag-and-probe procedure is performed using muons as tags and ID tracks as probes. The baseline sample for the reconstruction efficiency measurement includes a large number of non-muon tracks, which must be subtracted. This is done in bins of probe-track p T (p trk T ) and bins of the angular separation between the tag-muon and the probe-track, ∆R µ+trk . To describe the J/ψ peak and the background, a small range in tag-muon plus probe-track invariant mass, m µ+trk ∈ [2600, 3500] MeV, is fit to a double Gaussian function plus an exponential function and a second-order polynomial. In each p trk T or ∆R µ+trk bin, the ratio of the J/ψ peak component integral to the full shape (J/ψ plus background) integral is used as a weight to correct the p trk T or ∆R µ+trk shape itself. This is done separately for the probe-track distributions (denominators) and the muon-matched probe-track distributions (numerators). The ratio of the above two weighted distributions is defined as the reconstruction efficiency per p trk T or ∆R µ+trk bin. The efficiency measured with this approach in data is compared with the one from simulation and the difference at small ∆R µ+trk results in an uncertainty of 13.1% per event.

Analysis procedure
The analysis procedure is divided into four steps. First, events containing three high-quality muon objects with a combined invariant mass of less than 2.5 GeV are selected. These muons are required to originate from a common vertex. Second, a loose selection is applied to this sample to obtain a background sample that can be used to train the BDT, which is constructed using the TMVA toolkit [23]. The loose selection cuts (using a number of vertex quantities as well as kinematic quantities) are chosen to obtain a large background sample for training, while rejecting background that is kinematically inconsistent with the signal. Before training the BDT, the data events are divided into three regions based on the three-muon mass value. These are the blinded region (which includes the signal region), a sideband region and a BDT training region as defined in Table 1. Third, a tight selection (tightening the loose selection with a few additional cuts) is applied while simultaneously placing an initial loose cut on the BDT score, denoted by x>x 0 . The x>x 0 cut removes background-like events having a very low BDT score, while the tight selection further reduces the background in the blinded and sideband regions. Fourth, the background rejection as a function of the BDT cut is studied using data events in the sideband region passing the tight+x>x 0 selection. This allows to optimise the final cut on the BDT score, denoted by x>x 1 . The statistical analysis is performed for the tight+x>x 1 selection.  The signal region (SR) is defined as an interval around the τ lepton mass with a half-width corresponding to twice the resolution of the three-muon mass, σ s = 32 MeV, as obtained from the signal MC sample. The analysis was blinded in a slightly wider region to allow variation of the signal region definition. The signal MC sample is divided into two independent samples. One signal sample is used for the BDT training while the second signal sample is used for estimating the A s × s . The background in the signal region is estimated from a fit to the three-muon mass distribution in the sidebands (SB) using the tight+x>x 0 selection. This estimate is then scaled down to the final BDT score cut, x 1 , using a fit to the BDT shape as explained below.

Object selection
Muons are selected to have a transverse momentum greater than 2.5 GeV and are required to pass stringent requirements on the track quality and the associated hits in both the ID and the MS. Only combined ID+MS measurements of track parameters are used. Several matching criteria [22] are imposed to reject non-muon tracks (e.g. tracks from hadron decays in flight). The performance of muon identification is validated in two dedicated dimuon control regions. One region is populated with muons from J/ψ → µµ decay (in 2850 < m 2µ < 3350 MeV), while the second region has an enhanced fraction of non-muon tracks (in events with m 2µ < 750 MeV).
Events with at least three selected muons are considered. All possible three-muon combinations are used as inputs to a vertex fit. The primary vertex (PV) is also refitted after removing the three tracks. Due to the τ lepton lifetime, the three-muon vertex is often separated from the PV. The characteristics of the separation between the three-muon vertex and the PV are therefore used to distinguish signal from background. Particularly, the two projections of the three-muon vertex displacement with respect to the PV in the transverse plane are used; L xy =L T cos θ xy and a 0 xy =L T sin θ xy where L T is the transverse component of the vector connecting the PV and the three-muon vertex and cos θ xy = The threemuon vertex fit probability, p-value, is also used (as calculated from the vertex fit χ 2 and degrees of freedom). After fitting all possible vertices, exactly one three-muon candidate is allowed per event, satisfying m 3µ < 2500 MeV and |Q 3µ | = 1 where Q 3µ is the sum of the charges of the three-muon tracks.
Jets are used to separate the signal from the multi-jet backgrounds (predominantly HF), where more hadronic activity is expected. The jets are reconstructed from topological clusters formed in the calorimeter using the anti-k t algorithm [24] with a radius parameter R = 0.4. The jets are calibrated to the hadronic energy scale using energy-and η-dependent correction factors derived from simulation and with residual corrections from in situ measurements. A detailed description of the jet energy scale measurement and its systematic uncertainties can be found in Ref. [25]. Jets found within a cone of ∆R = 0.2 around a selected three-muon candidate are removed. Jets are required to have p T > 30 GeV and |η| < 2.8; only the leading jet satisfying these criteria is considered. There is no veto of events with more than one jet satisfying these criteria. The leading jet and the three-muon momenta are summed vectorially to define Σ = p jet + p 3µ with Σ T being the magnitude of its transverse component. For events where there are no jets satisfying these criteria (the majority of events for the signal), Σ is simply p 3µ .
The E miss T is calculated as the negative vector sum of the transverse momenta of all high-p T objects reconstructed in the event, as well as a term for other activity in the calorimeter [26]. Clusters associated with electrons, hadronic τ lepton decays and jets are calibrated separately, with other clusters calibrated at the EM energy scale. This E miss T is denoted hereafter by E miss T,cal . In addition, a track-based missing transverse momentum (E miss T,trk ) is calculated as the negative vector sum of the transverse momenta of tracks with |η| < 2.5, p T > 500 MeV and associated with the primary vertex. Both the calorimeter-based and track-based measurements of the E miss T are used.
Several kinematic variables are defined from the reconstructed objects listed above. Two transverse masses are defined using the three-muon transverse momentum (p 3µ T ) as m T = 2p 3µ T E miss T (1 − cos ∆φ 3µ ), where ∆φ 3µ is the angle between the E miss T and p 3µ T directions in the transverse plane. In these definitions, E miss T can be either E miss T,cal or E miss T,trk to obtain m cal T or m trk T respectively. The ∆φ 3µ terms are ∆φ cal 3µ or ∆φ trk 3µ respectively. Similarly, the ∆φ Σ T variable is the angle between the E miss T and Σ T directions in the transverse plane. This adds two additional angles, ∆φ cal Σ T and ∆φ trk Σ T for E miss T,cal and E miss T,trk respectively. These ∆φ Σ T variables provide good separation when a hard jet is found and thus Σ T deviates from p 3µ T in magnitude and direction.

Loose event selection
After the three-muon candidates are formed from the selected muons, a loose event selection is performed, maintaining a signal efficiency of about 80% while rejecting about 95% of the background. This loose selection includes cuts on the displacement of the vertex from the PV, requirements on the three-muon kinematics and on the presence of other tracks (track isolation), and requirements on quantities involving E miss T,cal and E miss T,trk . The loose selection comprises the following requirements: • The L xy significance, S (L xy )=L xy /σ L xy , must satisfy −10<S (L xy )<50, where σ L xy is the uncertainty in the L xy .
• The a 0 xy significance, S (a 0 xy )=a 0 xy /σ a 0 xy , must satisfy S (a 0 xy )<25, where σ a 0 xy is the uncertainty in a 0 xy . • The three-muon track-fit probability product, P trks = p 1 × p 2 × p 3 (where p i is the track fit p-value of track i), must satisfy P trks > 10 −9 .
• The three-muon transverse momentum must satisfy p 3µ T > 10 GeV. • The missing transverse energies, E miss T,cal and E miss T,trk , must both satisfy 10 < E miss T < 250 GeV.
• The transverse masses, m cal T and m trk T , must both satisfy m T > 20 GeV. • The three-muon track isolation is obtained from the sum of the p T of all tracks with p trk T > 500 MeV in a cone of ∆R The loose cuts on the significances, S (L xy ) and S (a 0 xy ), are applied to allow the three-muon vertex to be separated from the PV, while still being compatible with the τ lepton lifetime. The requirement on P trks imposes a goodness-of-fit criterion on the three-muon candidate. This value is based on an examination of signal-like events found in the sideband region in the data. As this is not the only quality requirement imposed on the individual muon objects, it is kept loose in this part of the selection. The efficiency for this cut to select signal events is ∼ 98%, while it is rejecting ∼ 13% of the background events. The kinematic and the isolation variables are very effective in separating the W boson properties of the signal from the HF and the light-flavour multi-jet background, which tend to be non-isolated and with low values of p T , E miss T and m T . The associated cuts remain very loose in this part of the selection to ensure that the sample sizes are large enough for the BDT training.

Multivariate analysis
The events passing the loose selection described above are used as input to the BDT training. There are 6649 events passing the loose selection in the signal MC sample (out of 10 5 ), where 6000 of these events are used for the BDT training and the rest are used for testing the BDT output. Similarly, the number of data events passing the loose selection in the training region is 4672, where 4000 of these events are used for the BDT training. The BDT input variables include kinematic distributions, track and vertex quality discriminants, vertex geometry parameters, and isolation. The following variables (sorted by their importance ranking) are used as inputs to the BDT: T /E miss T,trk − 1. 9. The difference between the three-muon and E miss T,cal directions, ∆φ cal 3µ . 10. Three-muon vertex fit probability, p-value.
11. The three-muon vertex fit a 0 xy significance, S (a 0 xy ). 12. The track fit probability product, P trks . 13. The three-muon transverse momentum, p 3µ T . 14. The number of tracks associated with the PV (after refitting the PV while excluding the three-muon tracks), N PV trk . 15. The three-muon vertex fit L xy significance, S (L xy ). 16. The calorimeter-based missing transverse momentum balance, p 3µ T /E miss T,cal − 1. This configuration was found to give the optimal balance between background rejection and signal efficiency.
The Σ T variable is introduced to avoid vetoing events with at least one jet fulfilling the requirements listed in section 5.1. Although the majority of signal events do not have jets, it is found that keeping such events increases the A s × s by ∼ 15% and also ultimately leads to better rejection power, owing to the significantly larger training and sideband samples. The variables ∆φ cal trk , ∆φ cal 3µ , p 3µ T /E miss T,trk − 1, p 3µ T /E miss T,cal − 1 and N PV trk are complementary to the E miss T -related variables used in the loose selection as well as here. These variables are also very effective in distinguishing the W → τν production of the signal from the HF and light-flavour multi-jet background. The vertex p-value is a variable complementary to the S (L xy ) and S (a 0 xy ) variables used in the loose selection as well as here. The HF and light-flavour multijet backgrounds have mostly random combinations of selected muon objects which do not originate from the same vertex. This variable peaks at very low values for the background while for the signal it is distributed uniformly and thus provides excellent separation.
After training the BDT with data events from the training region and signal MC events from the first signal MC sample, the BDT response is calculated for the data events in the sidebands and for events in the second signal sample. The BDT score, x, ranges between -1 and +1. Events with a very low BDT score, within −1 ≤ x ≤ −0.9 are removed from further consideration, defining x 0 ≡ −0.9.
In order to assess potential modelling problems in the signal MC sample, the BDT input distributions and the BDT response are validated against single-muon data. These data contain mainly W → µν events with a small fraction (< 10%) of background. The single-muon selection is formulated to be as close as possible to the main analysis selection where the differences are mostly driven by the different triggers used (one single-muon trigger with no isolation requirement and with a threshold of 24 GeV is used in the validation) and the exclusion of variables which do not have equivalents in the W → µν case, e.g. the three-muon vertex variables. The training samples used for this validation study, for both data and signal, are the same samples as used in the main analysis, constructed with the same loose selection as described in the previous section. All input variables are used for the training, excluding the p-value, S (L xy ), S (a 0 xy ) and P trks , which cannot be calculated in a single-muon (W → µν) selection. The resulting BDT setup is hereafter referred to as "partial BDT". After training the partial BDT, the response is tested on the second signal sample and on the single-muon data, using the single-muon selection and where the three muon objects in the signal sample are treated as one object (muon). The N PV trk distribution of the signal sample is also modified by subtracting two tracks to reflect the difference with respect to a single-muon selection. The responses in data and simulation are compared and are found to agree within 10% throughout most of the phase-space for all variables. The ratio of the partial BDT responses for the single-muon data and signal MC events is used as an event weight while applying the full selection and calculating the weighted A s × s as described in the next sections. The difference between the weighted and unweighted A s × s is found to be 4% and is taken as a modelling uncertainty.
Any variable which may bias the BDT response by only selecting events very close to the τ lepton mass is not included in the BDT input list. The distribution of the three-muon mass has been examined in several bins of x above x 0 using both the loose and the tight samples, where no hint of potential peaking background around the τ lepton mass has been found. In addition, the shape of the three-muon mass distribution has been found to be insensitive to the BDT cut, as expected given the small correlation coefficient between x and m 3µ , which is found to be about -0.05.

Tight event selection
Additional tight cuts are applied after the BDT training and the application of the x>x 0 cut on the BDT score. The following requirements are tightened or added: • A number of the loose requirements are tightened, namely P trks > 8 × 10 −9 , m cal T > 45 GeV, m trk T > 45 GeV and 1<S (L xy )<50. • Three-muon vertex fit probability must have p−value > 0.2.
• The angle between the Σ T and E miss T,cal (E miss T,trk ) directions is required to be ∆φ cal Σ T > 2 (∆φ trk Σ T > 2).
• The same-charge two-muon mass, m SS , and opposite-charge two-muon mass, m OS1 or m OS2 , satisfy m SS > 300 MeV, m OS1 > 300 MeV and m OS2 > 300 MeV, where m OS1 (m OS2 ) is the mass of the two opposite-charge muon pairs with the highest (second highest) summed scalar p T among the three muons.
• The event is rejected if |m OS − m ω | < 50 MeV or |m OS − m φ | < 50 MeV if either of the p 3µ T , the E miss T,cal or the E miss T,trk is lower than 35 GeV. • The event is rejected if |m OS − m φ | < 50 MeV if |m 3µ − m D s | < 100 MeV.
In the above notation, m OS is m OS1 or m OS2 and m ω , m φ and m D s are the masses of the ω, φ and D s mesons respectively, taken from Ref. [2].
The requirement on the three-muon vertex fit probability is applied in order to ensure a high-quality fit. The cuts on ∆φ cal Σ T (∆φ trk Σ T ) are applied in order to further suppress the HF and multi-jet background where the three-muon candidate is typically produced within or near a jet.
The first two-muon mass requirement is applied to suppress candidates originating from one prompt muon object and two muons from a converted photon. The second requirement on the two-muon masses is applied to prevent the low-mass mesons, ρ/ω and φ, from entering into the region close to the τ lepton mass when combined with an additional track. In the selected three-muon event sample, these resonances appear as two clear peaks in the mass distribution of oppositely charged muon pairs in data. Since the resonances lie in the middle of the signal distribution, the low-p T and E miss T requirement ensures that these can still be distinguished from the signal, and thus it removes the resonances while still maintaining a high enough signal efficiency. Finally, the last requirement is applied to remove a potential D s → π + φ(µµ) contamination from the high-mass sideband. The cuts listed above comprise the tight selection where the tight+x>x 0 selection is used to estimate the background for any cut on x above x 0 . Figure 1 shows the three-muon mass distribution and the BDT response distribution. Figures 2 and 3 show the distributions of the BDT inputs sorted by the separation rank as reported by TMVA during the BDT training. Figure 4 shows the distributions of the complementary variables which are used in the loose or tight selection but not in the BDT.

Background estimation
The events passing the tight+x>x 0 selection are used to estimate the expected number of background events in the signal region for higher cuts on x as described below.
The signal MC and sideband data BDT responses are shown in Fig. 5 after the tight+x>x 0 selection. The distinct shapes illustrate the power of the method in separating the background from the signal. The analytical function also shown in Fig. 5 is a result of a fit to the sideband data, excluding the blinded region, using an unbinned maximum-likelihood estimator. The fit function used is a 0 +a 1 (x +1) a 2 +a 3 (x + 1) a 4 , where a i are the free fit parameters. The parameter a 2 is required to be negative while the other are required to be non-negative. This function can exhibit rising behaviour at both ends of the x distribution (x → ±1) and it is used to scale the quantities measured in x>x 0 to the corresponding quantities in x>x 1 as explained below.
The three-muon mass distribution of the tight+x>x 0 data is fit simultaneously in the two sidebands to a second-order polynomial in m 3µ while excluding the blinded region. This is also done with an unbinned maximum-likelihood estimator. The integral of the resulting fit function in the signal region gives the =8 TeV, 20.3 fb s  Figure 1: The three-muon mass distribution in (a) and the BDT score in (b). The BDT score distribution of the data is shown for the sideband region. The loose data are shown as hollow circles, while the loose signal MC events are shown as light solid grey area. The tight+x>x 0 data are shown as the solid black circles, while the tight+x>x 0 signal MC events are shown as the dark solid grey area. The area of the signal MC shapes is normalised to the area of the loose data shapes and the relative normalisation difference between the loose and the tight+x>x 0 MC signal distributions prior to the normalisation is maintained. For illustration, the signal is not constrained to the SR. expected number of background events, N b (x 0 ) in the signal region before applying the final x 1 cut. The statistical uncertainty of N b (x 0 ) is calculated by scaling the statistical error in the number of events in the sidebands, according to the ratio of analytical integrals in the signal region and sidebands. Figure 6 shows the three-muon mass distribution in the sidebands for the tight+x>x 0 selection as black points together with the fit result. The signal is also shown for reference, scaled up arbitrarily to match the scale of the data.
For any x 1 cut value above x ∼ 0.6, where most of the signal is expected, the estimated N b (x 0 ) in the signal region can be then scaled down according to the ratio of the integrals of the BDT analytical function above and below this cut. This ratio is denoted hereafter by R(x 1 ). The extrapolation procedure can be written as N b (x 1 ) = R(x 1 ) × N b (x 0 ) where N b (x 1 ) is estimated in the signal region for x>x 1 .

Uncertainties and optimisation
The sources of systematic uncertainty associated with the extrapolation procedure in the background estimation are the BDT and sideband fit function choice, the definition of the sideband ranges and the definition of x 0 . To estimate this uncertainty, each of these definitions and choices is varied individually while calculating N b (x 0 ) and R(x 1 ). For each fit function (BDT and sideband), different parameterisations are considered. In addition, to construct the variation of the tight+x>x 0 sample with which the two fits are performed, nine different sideband range variations and ten different x 0 variations are used. The fits, and consequently also the extrapolation procedure, are found to be stable against these variations. The dominant uncertainty component is the impact on R(x 1 ) of varying the sideband ranges definition. The differences from the nominal values of R(x 1 ) and N b (x 0 ) are summed in quadrature and are translated to uncertainties in N b (x 1 ). The systematic uncertainty associated with the extrapolation procedure used to obtain N b (x 1 ) increases with x 1 from ∼ 45% at x 1 = 0.6 to ∼ 80% at x 1 1. The statistical uncertainty of N b (x 1 ) is ∼ 19%, independent of x 1 .
The systematic uncertainty in the signal acceptance times efficiency has contributions from reconstruction (13.1%), trigger (11%) and MC modelling (4%) as discussed in the previous sections. In addition, there is a small (2.1%) contribution due to jet and E miss T calibration. The number of τ leptons produced via the W → τν channel and its uncertainty (3.9%) are estimated as described in Section 3. These uncertainties are independent of x 1 in the range of interest.
The BDT cut is optimised by minimising the expected upper limit on the branching fraction given in Eq. (1), where N s becomes the upper limit on the number of observed events above the expected background level in a narrow region around the τ lepton mass. The procedure is performed by varying x 1 between 0.6 and 1.0 in steps of 0.001 while extracting N b (x 1 ) and its associated errors as explained above.
To obtain the upper limit on N s for each x 1 cut, a single-bin counting experiment is performed using the HistFitter [27] statistical framework, supplied with N b (x 1 ) and its uncertainties. For compatibility with previous searches, the limit on N s and on Br(τ → 3µ) is given at 90% confidence level (CL). In each iteration, A s × s is calculated for the specific x 1 cut using a signal sample that is different from the one used for the BDT training.
During the iterative optimisation process, the extrapolation of the number of events in the sideband region to high x 1 cuts using the BDT shape is tested against a cut-and-count procedure. The two procedures are found to agree very well within the uncertainties, and the extrapolation procedure gives a more conservative result throughout the examined x 1 range. The resulting optimal cut is at x 1 = 0.933. Figure 6 shows the three-muon mass distributions in the full mass range, including the blinded region, for the tight+x>x 1 selection in red squares. Only one event with a three-muon mass of 1860 MeV survives the selection in the full mass range (sideband and blinded regions). This event is found in the range between the signal region and the right sideband region and it does not affect the background estimation or the observation in the signal region.

Results
The event counts entering the different regions at the different steps of the analysis for signal and data are given in  The signal acceptance times efficiency is calculated from the second signal MC sample after applying the full tight+x>x 1 selection. This selection corresponds to A s × s = 0.0231±0.0005 jet/E miss T ±0.0009 modelling ± 0.0025 trig ± 0.0030 reco . With this selection, the expected background yield is N b (x 1 ) = 0.193 ± 0.131 syst ± 0.037 stat . The systematic uncertainty on A s × s is dominated by the uncertainties in the reconstruction and trigger efficiency measurements. The systematic uncertainty on N b (x 1 ) is dominated by the uncertainty in the extrapolation of the background from the tight+x>x 0 selection to the tight+x>x 1 selection. The systematic uncertainties in N b are taken into account when calculating the limit on the number of signal events, N s , via one nuisance parameter. The systematic uncertainties in the product (A s × s ) · N W→τν are summed in quadrature and taken into account as the uncertainty in the signal via one nuisance parameter when calculating the limit. The expected (median) limit on the branching fraction for N o = N b (x 1 ) is 3.94 × 10 −7 at 90% CL. No events are observed in the signal region and the observed limit on the branching fraction is therefore 3.76 × 10 −7 at 90% CL.

Conclusions and outlook
This article presents a search with the ATLAS detector for neutrinoless τ → 3µ decays using 20.3 fb −1 of 2012 LHC pp collision data, utilising τ leptons produced in W → τν decays. No events are observed in the signal region for the final selection while 0.193 ± 0.131 syst ± 0.037 stat background events are expected. This results in an observed (expected) upper limit of 3.76 × 10 −7 (3.94 × 10 −7 ) on Br(τ → 3µ) at 90% CL. Although this limit is not yet competitive with searches performed at B-factories [7,8] and at LHCb [9], it demonstrates the potential of LHC data collected by ATLAS as a probe of lepton flavour violation in τ lepton decays. This analysis utilises single τ lepton production in an environment very different from B-factories, which rely on τ lepton pair production in e + e − collisions. The method and sample presented here were used to improve the ATLAS muon trigger and reconstruction of low-p T , collimated muons relevant to the τ → 3µ search. The analysis is limited by the number of W → τν decays and by the systematic uncertainty, which depends on the size of the data sample. With the much larger data sets anticipated at Run 2 of the LHC, the sensitivity of ATLAS to lepton-flavour-violating decays will be improved significantly.       T /E miss T,cal − 1 in (h). The loose data in the sidebands are shown as hollow circles, while the loose signal MC events are shown as light solid grey area. The tight+x>x 0 data in the sidebands are shown as the solid black circles, while the tight+x>x 0 signal MC events are shown as the dark solid grey area. The area of the signal MC shapes is normalised to the area of the loose data shapes and the relative normalisation difference between the loose and the tight+x>x 0 MC signal distributions prior to the normalisation is maintained. For illustration, the signal is not constrained to the SR.  The loose data in the sidebands are shown as hollow circles, while the loose signal MC events are shown as light solid grey area. The tight+x>x 0 data in the sidebands are shown as the solid black circles, while the tight+x>x 0 signal MC events are shown as the dark solid grey area. The area of the signal MC shapes is normalised to the area of the loose data shapes and the relative normalisation difference between the loose and the tight+x>x 0 MC signal distributions prior to the normalisation is maintained. For illustration, the signal is not constrained to the SR.