Search for light long-lived neutral particles produced in pp collisions at $$\sqrt{s} = 13~\mathrm {TeV}$$ and decaying into collimated leptons or light hadrons with the ATLAS detector

<jats:title>Abstract</jats:title><jats:p>Several models of physics beyond the Standard Model predict the existence of dark photons, light neutral particles decaying into collimated leptons or light hadrons. This paper presents a search for long-lived dark photons produced from the decay of a Higgs boson or a heavy scalar boson and decaying into displaced collimated Standard Model fermions. The search uses data corresponding to an integrated luminosity of 36.1 fb<jats:inline-formula><jats:alternatives><jats:tex-math>$$^{-1}$$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow/><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:math></jats:alternatives></jats:inline-formula> collected in proton–proton collisions at <jats:inline-formula><jats:alternatives><jats:tex-math>$$\sqrt{s} = 13$$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>13</mml:mn></mml:mrow></mml:math></jats:alternatives></jats:inline-formula> <jats:inline-formula><jats:alternatives><jats:tex-math>$$\text {Te}\text {V}$$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mtext>Te</mml:mtext><mml:mspace/></mml:mrow></mml:math></jats:alternatives></jats:inline-formula> recorded in 2015–2016 with the ATLAS detector at the Large Hadron Collider. The observed number of events is consistent with the expected background, and limits on the production cross section times branching fraction as a function of the proper decay length of the dark photon are reported. A cross section times branching fraction above 4 pb is excluded for a Higgs boson decaying into two dark photons for dark-photon decay lengths between 1.5 mm and 307 mm.</jats:p>

section_content: Several extensions of the Standard Model (SM) predict the existence of a dark sector weakly coupled to the SM [1] [2] [3] [4]. Depending on the structure of the dark sector and its coupling to the SM, some unstable dark states may be produced at colliders, and could decay into SM particles with sizeable branching fractions. In order to avoid a new long-range force, a dark Higgs boson is introduced in such scenarios, to give mass to the dark gauge bosons. The dark Higgs boson may also lead to an exotic decay mode of the Higgs boson, via mixing between the two Higgs sectors, which is one of the favoured production modes that may be probed at the Large Hadron Collider (LHC). This is the mode explored in this search. Branching fractions of up to 10% are currently not excluded for Higgs-boson decays into exotic final states [5, 6]. This paper investigates the case where the two sectors couple via a vector portal, in which a dark photon (γ d ) mixes kinetically with the SM photon and decays into SM leptons and light quarks [7] [8] [9]. The kinetic mixing term ( ), which can vary over a wide range of values, ∼ 10 -11 -10 -2, determines the lifetime of the dark photon. For a small kinetic mixing value, the γ d has a long lifetime, so that it decays at a macroscopic distance from its production point. This analysis focuses on small values of the kinetic mixing term, < 10 -5, and a dark photon mass range between twice the muon mass and the twice the tau mass. Due to their small mass, the dark photons are expected to be produced with large boosts, resulting in collimated groups of leptons and light hadrons in a jet-like structure, referred to hereafter as dark-photon jets (DPJs). 

The search for displaced DPJs presented in this paper uses the dataset collected by the ATLAS detector during 2015-2016 in proton-proton (pp) collisions at a centre-of-mass energy √ s = 13 TeV, corresponding to an integrated luminosity of 36. 1 fb - 1. The analysis exploits multivariate techniques for the suppression of the main multi-jet background, optimised for the different DPJ channels. This technique allows the exploitation of the fully hadronic signature for the first time in ATLAS DPJ searches, resulting in increased sensitivity compared with previous ATLAS results using the data collected in 2011 and 2012 at 7 and 8 TeV respectively [10, 11]. The results are complementary to those from related ATLAS searches for prompt DPJs using 7 and 8 TeV data [12] [13] [14], which probed higher values of, and for displaced dimuon vertices using 13 TeV data [15], which probed higher dark photon mass values. Related searches for dark photons were conducted by the CDF and D0 collaborations at the Tevatron [16, 17] and by the CMS [18] [19] [20] [21] and LHCb [22, 23] collaborations at the LHC. Additional constraints on scenarios with dark photons are extracted from, e. g., beam-dump and fixed-target experiments [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34], e + e -colliders [35] [36] [37] [38] [39] [40] [41] [42] [43], electron and muon anomalous magnetic moment measurements [44] [45] [46] and astrophysical observations [47, 48]. Given the various constraints, a displaced dark photon with a kinetic mixing term < 10 -5 is allowed for γ d masses greater than 100 MeV.
section_name: Introduction

section_content: ATLAS [49] is a multipurpose detector at the LHC, consisting of an inner detector (ID) contained in a superconducting solenoid, which provides a 2 T magnetic field parallel to the beam direction, electromagnetic and hadronic calorimeters (ECAL and HCAL) and a muon spectrometer (MS) that has a system of three large air-core toroid magnets, each composed of eight coils. 

The ID provides measurements of charged-particle momenta in the region of pseudorapidity |η| ≤ 2. 5. 1 The highest spatial resolution is obtained around the vertex region using semiconductor pixel detectors arranged in four barrel layers [50, 51] at average radii of 3. 3 cm, 5. 05 cm, 8. 85 cm, and 12. 25 cm, and three discs on each side, covering radii between 9 cm and 15 cm. The pixel detector is surrounded by four layers of silicon microstrips covering radial distances from 29. 9 cm to 56. 0 cm. These silicon detectors are complemented by a transition radiation tracker (TRT) covering radial distances from 56. 3 cm to 106. 6 cm. 

The ECAL and HCAL calorimeter system covers |η| ≤ 4. 9, and has a total depth of 9. 7 interaction lengths at η = 0, including 22 radiation lengths in the ECAL. The ECAL barrel starts at a radius of 1. 41 m and ends at 1. 96 m with a z extension of ±3. 21 m, covering the |η| ≤ 1. 475 interval. In the 1. 37 ≤ |η| ≤ 3. 2 region, the ECAL endacap starts at z ± 3. 70 m and end at z ± 4. 25 m. The HCAL barrel starts at a radius of 2. 28 m and ends at 4. 25 m with a z extension of ±4. 10 m, covering the |η| ≤ 1. 0 interval. In the endcaps regions up to |η| ≤ 4. 9, the HCAL starts at z ± 4. 3 m and ends at z ± 6. 05 m. 

The MS provides trigger information and momentum measurements for charged particles in the pseudorapidity ranges |η| ≤ 2. 4 and |η| ≤ 2. 7 respectively. It consists of one barrel (|η| ≤ 1. 05) and two endcaps (1. 05 ≤ |η| ≤ 2. 7), each with 16 sectors in φ, equipped with fast detectors for triggering and with chambers for reconstructing the tracks of the outgoing muons with high spatial precision. The MS detectors are arranged in three stations at increasing distances from the IP: inner, middle and outer. Three planes of MS trigger chambers are located in the middle and outer stations. The toroidal magnetic field allows precise reconstruction of charged-particle momenta independent of the ID information. The ATLAS trigger system has two levels [52], level-1 (L1) and the high-level trigger (HLT). The L1 trigger is a hardware-based system using information from the calorimeters and MS. It defines one or more regions-of-interest (RoI), which are geometric regions of the detector identified by (η, φ) coordinates, containing interesting physics objects. The L1 trigger reduces the event rate from the LHC crossing frequency of 40 MHz to a design value of 100 kHz. L1 RoI information provides a seed for the reconstruction of physics objects by the HLT, a software-based system that can access information from all subdetectors. It is implemented in software running on a PC farm that processes the events and reduces the rate of recorded events to 1 kHz.
section_name: The ATLAS detector

section_content: Among the numerous models predicting dark photons, one class particularly interesting for the LHC features a hidden sector communicating with the SM through the Higgs portal for production and through vector portal for decay. The benchmark model used in this analysis is the Falkowski-Ruderman-Volansky-Zupan (FRVZ) model [8, 9], where a pair of dark fermions f d 2 is produced via a Higgs boson (H) decay. Two different cases of this model are considered, involving the production of either two or four dark photons. In the first case, shown in Figure 1 (left), each dark fermion decays into a γ d and a lighter dark fermion assumed to be the hidden lightest stable particle (HLSP). In the second case, shown in Figure 1 (right), each dark fermion decays into an HLSP and a dark scalar s d that in turn decays into a pair of dark photons. 

In general, dark-sector radiation [53] can produce extra dark photons. The number of radiated dark photons is proportional to the size of the dark gauge coupling α d [7]. The dark radiation is not considered in this signal model, which corresponds to an assumed dark coupling α d 0. 01. beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = -ln tan(θ/2). The vector portal communication of the hidden sector with the SM is through kinetic mixing of the dark photon and the standard photon

where B µν and b µν denote the field strengths of the electromagnetic fields for the SM and dark sector respectively, and is the kinetic mixing parameter. A dark photon with a mass m γ d up to a few GeV that mixes kinetically with the SM photon will decay into leptons or light mesons, with branching fractions that depend on its mass [8, 54, 55]. 

The mean lifetime τ, expressed in seconds, of the γ d is related to the kinetic mixing parameter [56] by the relation

Equation ( 1 ) is an approximate expression based on the full relation in Ref. [55].
section_name: Benchmark model

section_content: The analysis presented in this paper uses √ s = 13 TeV pp collision data recorded by the ATLAS detector during the 2015-2016 data-taking periods. Only runs in which all the ATLAS subdetectors were operating normally are selected. The total integrated luminosities are 3. 2 fb -1 and 32. 9 fb -1 for 2015 and 2016 respectively. Data were collected using a set of dedicated triggers that were active during collision bunch crossings as well as during empty and unpaired bunch-crossing slots. The LHC configuration for pp collisions contains 3564 bunch-crossing slots per revolution. An empty bunch-crossing is defined as a slot in which neither beam is filled with protons, and in addition is separated from filled bunches by at least five unfilled bunches on each side. Data collected during empty bunch crossings, referred to as the cosmic dataset, are used for the estimation of the cosmic-ray background. The ratio of the number of filled to empty bunch crossings, F CR = 2. 1, is used to scale the number of events in the cosmic dataset to that in the pp collision data. In unpaired bunch crossings, protons are present in only one of the two beams. Data taken during unpaired bunch crossings are used to study characteristic features of beam-induced backgrounds [57] (BIB) and are referred to as the BIB dataset. 

Monte Carlo (MC) simulation samples were produced for the model considered in this paper and are summarised in Table 1. 

Samples were generated for the Higgs boson mass of 125 GeV, and for a hypothetical beyond-the-SM (BSM) heavy scalar boson with a mass of 800 GeV, considering only the dominant gluon-gluon fusion (ggF) production mechanism. The ggF Higgs boson production cross section in pp collisions at √ s = 13 TeV, estimated at next-to-next-to-leading order (NNLO) [58] [59] [60] [61], is σ SM = 43. 87 pb for m H = 125 GeV. The BSM heavy scalar with a mass of 800 GeV production cross section is conventionally assumed to be σ = 5 pb. 

The mass of the hidden fermion m f d 2 and of the hidden scalar m s d were chosen to be low relative to the Higgs boson mass. Due to the production from a two-body decay of the Higgs boson generated at rest in the transverse plane, events with two back-to-back DPJs are expected. This is also the case leading to four dark photons where each DPJ consist of two collimated dark photons. 

The dark-photon mass was chosen to be 0. 4 GeV, above the pion pair mass threshold, and the γ d decay branching fractions (B) are expected to be B(γ

In the generated samples, the proper decay length cτ of the γ d was chosen such that ∼80% of the decays occur in the volume delimited by the muon trigger chambers (i. e. up to 7 m in radius and 13 m along the z-axis). Since the analysis is sensitive to a wide range of mean proper lifetimes, a weighting method is used to extrapolate the signal efficiency to other mean proper lifetimes. 

All MC samples described above were generated at leading order using M G 5_aMC@NLO 2. 2. 3 [62] interfaced to P 8. 210 [63] for parton shower generation. The A14 set of tuned parameters (tune) for parton showering and hadronisation [64] was used together with the NNPDF2. 3LO parton distribution function (PDF) set [65]. One of the main SM backgrounds in this analysis is multi-jet events. Such events were simulated to perform background studies and to evaluate systematic uncertainties. The MC samples were generated with P 8. 210 using the same tune and PDF as for the signal samples. 

Potential sources of background also include W+jets, Z+jets, t t, single-top-quark, WW, WZ, and ZZ events. Simulation samples are used to study these backgrounds. The W+jets, Z+jets, WW, WZ, and ZZ events were generated using S 2. 2. 2 [66] with the NNPDF 3. 0 NNLO [67] PDF set. Single-top-quark and t t MC samples were generated using P -BOX 1. 2856 [68] [69] [70] [71] and P 6. 428 [72] with the Perugia2012 [73] tune for parton showering and hadronisation, and CT10/CTEQ6L1 [74, 75] PDF sets. 

Data and MC samples of J/ψ → µµ events are used to evaluate systematic uncertainties in muon trigger and reconstruction efficiencies. The MC sample was generated using P 8+P ++ [76] with the A14 tune for parton showering and hadronisation, and the CTEQ6L1 PDF set. The J/ψ → µµ data sample was selected in 2015-2016 pp collisions using the triggers described in Ref. [77]. 

The generated MC events were processed through a full simulation of the ATLAS detector geometry and response [78] using the G 4 [79] toolkit. The simulation included multiple pp interactions per bunch crossing (pile-up), as well as the detector response to interactions in bunch crossings before and after the one producing the hard interaction. To model the effect of pile-up, simulated inelastic pp events were overlaid on each generated signal and background event. The multiple interactions were simulated with P 8. 210 using the A2 tune [80] and the MSTW2008LO PDF set [81]. 

5 Definition of the dark-photon jets
section_name: Data and simulation samples

section_content: Displaced DPJs are reconstructed with criteria that depend on the γ d decay channel. A γ d decaying into a muon pair is searched for by looking for two closely spaced muon tracks in the MS, while a γ d decaying into an electron or pion pair, given the high boost of the γ d, is searched for as an energy deposit in the calorimeters identified as a single narrow jet. MC simulations show that DPJs containing two dark photons both decaying into an electron or pion pair are reconstructed as a single jet. 

Tracks that are reconstructed in the MS and are not matched to any track in the ID are used to identify displaced γ d decays into muons. Since the ID track reconstruction in ATLAS [82] requires at least one hit in one of the two innermost pixel layers, this analysis is sensitive only to displaced γ d decays occurring after the first pixel layers. The search is limited to the pseudorapidity interval |η| < 2. 5, corresponding to the ID coverage, to ensure that selected muons are isolated from ID tracks. Muons with pseudorapidity in the range 1. 0 ≤ |η| ≤ 1. 1 are rejected to avoid the transition region of the MS between barrel and endcap. 

In order to reconstruct γ d decays that occur outside of the innermost layer of muon chambers but before the first MS trigger chamber, muons are required to have at least one hit in two of the three MS tracking station. 

Jets used in this search are reconstructed from clusters [83] of energy deposits in the ECAL and HCAL using the anti-k t algorithm [84, 85] with radius parameter R = 0. 4. The search is limited to γ d decays into electron or hadron pairs in the hadronic calorimeter. Jets produced in the HCAL are expected to be isolated from activity in the ID, with a high ratio of energy deposited in the HCAL (E HCAL ) to energy deposited in the ECAL (E ECAL ), and appear narrower than ordinary jets. The standard jet-cleaning requirements [86] applied in most ATLAS analyses reject jets with high values of E HCAL /E ECAL. A dedicated cleaning algorithm for jets created in the HCAL is applied instead, with no requirements on the ratio E HCAL /E ECAL. 

Jets are required to have transverse momentum p T ≥ 20 GeV and |η| < 2. 5. In addition, the weighted time of the energy deposit in the calorimeter cells is required to be in the range [-4 ns, 4 ns] of the expected arrival time for particles produced at t = 0 (bunch-crossing time) and moving with the speed of light, to reduce cosmic-ray background and BIB jets. 

DPJs are classified according to the number of muons and jets found within a given cone of angular size ∆R ≡ (∆φ) 2 + (∆η) 2 around a muon or jet candidate with the highest transverse momentum. The cone size is fixed to ∆R = 0. 4, since the MC simulations show that this selection retains up to 90% of the dark-photon decay products in the H → 4γ d + X decay channel with m H = 125 GeV. The DPJ classification is summarised as follows:

• muonic-DPJ (µDPJ) -to select DPJs with all constituent dark photons decaying into muons, at least two muons are required and no jets are allowed in the cone. 

• hadronic-DPJ (hDPJ) -to select DPJs with all constituent dark photons decaying into electron or pion pairs in the HCAL, one jet is required and no muons are allowed to be in the cone. The electromagnetic fraction of the jet energy, defined as the ratio of the energy deposited in the ECAL to the total jet energy, is required to be less than 0. 4. This helps reduce the overwhelming background due to multi-jet production. 

Reconstructed DPJs with both muon and jet constituents are not considered in this analysis.
section_name: Dark-photon jet classification

section_content: Muonic-DPJs are reconstructed using a Cambridge-Aachen clustering algorithm [87] that combines all the muons lying within a cone of fixed size in (η, φ) space. The algorithm starts from the highest-p T muon, searching for additional muons within the ∆R = 0. 4 cone around the muon momentum vector. If a second muon is found in the cone, the axis of the cone is rotated to the vector sum of the momenta of the two muons, and the search is repeated until no additional muons are found in the cone. 

Cosmic-ray muons that cross the detector in time coincidence with a pp interaction constitute the main source of background to the muonic-DPJ. The cosmic dataset is used to study this background. A boosted decision tree (BDT) with gradient boosting, implemented in the TMVA framework [88], is trained to discriminate signal DPJs from the DPJ candidates that originate from cosmic-ray background. The BDT uses the following track variables, for each muon in the DPJ, to classify a DPJ as being from signal or background:

-longitudinal impact parameter z 0, defined as the minimum separation in the z-coordinate between the muon track and the primary vertex (PV);2 -arrival times measured by the trigger detectors of the MS; -pseudorapidity η;

-azimuthal angle φ. 

Even if the decay is displaced, signal muons point to the primary vertex because of the high boost of the dark photon, resulting in a narrow z 0 distribution peaking around zero. By contrast, cosmic-ray muons have a broad z 0 distribution. 

Cosmic-ray muons mainly come through the two shafts above the ATLAS detector, resulting in two well-defined peaks in the η and φ distributions. Each hit in the trigger detector of the MS provides a measurement of the time for the muon track, corrected by the time of flight assuming the pp interaction point as the origin of the muon [89]. The difference in time measured by the two layers in the middle station and in the outer station is thus useful for discriminating between cosmic-ray muons and collision muons. Since cosmic-ray muons are downward going, their arrival times in the layers in the upper part of the MS (0 < φ < π) are different from those of collision muons, which are upward-going in this part of the detector. In the lower part of the MS (π < φ < 2π), cosmic-ray muons and collision muons travel downwards, making hit timing less useful for separating between them. 

The cosmic dataset and the signal MC sample H → 2γ d + X with m H = 125 GeV are used for the training of the BDT. The gain in signal significance obtained from dedicated BDT training with the other signal MC samples is found to be negligible. Figure 2 (left) shows the BDT output (µBDT) for the constituent muons of the µDPJs: the distribution provides a clear separation between signal and background muons from cosmic rays. The µBDT output is required to be µBDT > 0. 21; the value is chosen to yield the highest signal significance, S/ √ S + B, where S is the number of signal events and B the number of background events.
section_name: Muonic-DPJ selection

section_content: Signal jets are discriminated from multi-jets using a second classifier also based on a BDT (hBDT). The following variables are used as input to the hBDT:

-jet width, defined as the p T -weighted sum of the ∆R between each energy cluster and the jet axis; -jet vertex tagger (JVT) output [90] ; -E HCAL /E ECAL ; -jet mass, as defined by the jet clustering algorithm [91] ; -jet charge, defined as the momentum-weighted charge sum constructed from tracks associated with the jet; tracks are associated with jets using ghost association [92] ;

-jet timing, defined as the energy-weighted average of the timing for each cell in the jet. 

The JVT is designed to differentiate between pile-up jets and jets originating from the PV. The algorithm uses a multivariate combination of track variables that are sensitive to pile-up. Since jets produced in the hadronic calorimeter have a JVT output distribution similar to that of pile-up jets, the JVT output is used for selection of hadronic-DPJs. Possible pile-up jets contamination is reduced by the analysis selection to a negligible level. 

The signal MC sample H → 2γ d + X with m H = 125 GeV and the simulated multi-jet background events are used for the BDT training. The gain in signal significance obtained from dedicated BDT training with the other signal MC samples is found to be negligible. Figure 2 (right) shows the BDT output for the hDPJs (hBDT). The peak at ∼ -0. 2 in the BDT distributions corresponds to jets with a JVT output that indicates a low pileup-up probability. The hBDT output is required to be hBDT > 0. 91; the value is chosen to yield the highest signal significance.
section_name: Hadronic-DPJ selection

section_content: The standard ATLAS triggers are optimised to select prompt events and are thus usually very inefficient in the selection of displaced objects. This search uses events selected by the logical OR of three dedicated triggers targeting displaced objects: two muon triggers and one calorimeter trigger. 

The L1 muon trigger used in this analysis requires hits in the middle stations to create a low-p T (≥ 6 GeV) muon RoI or hits in both the middle and outer stations for a high-p T (≥ 20 GeV) muon RoI. The muon RoIs have a ∆η × ∆φ spatial extent of 0. 2 × 0. 2 in the barrel and of 0. 1 × 0. 1 in the endcaps. L1 RoI information seeds the reconstruction of muon momenta by the HLT, which uses precision-chamber information to confirm or reject the L1 decision. 

The first muon trigger, the tri-muon MS-only [93], requires at least three L1 muons with p T ≥ 6 GeV in the event, confirmed by the HLT using only MS information. 

The second muon trigger, the muon narrow-scan, is specifically designed to select non-prompt collimated muons originating in the region between the first pixel layer and the first muon trigger plane. It requires an L1 muon with p T ≥ 20 GeV confirmed by the HLT using only MS information. At the HLT a 'scan' is then performed in a cone of ∆R = 0. 5 around this muon, looking for a second muon reconstructed using only MS information. The p T requirement on the second muon was increased from 6 GeV to 15 GeV during the course of the 2015-2016 data taking in order to stay within the allocated trigger-rate limits given the increasing luminosity delivered by the LHC. 

The calorimeter trigger, the CalRatio [93], is designed to select narrow jets produced in the hadronic calorimeter. At L1, the trigger requires a transverse-energy deposit of E T ≥ 60 GeV within a 0. 2 × 0. 2 (∆η × ∆φ) region in the pseudorapidity range |η| ≤ 2. 4. At the HLT, jet reconstruction is then performed with the anti-k t algorithm using a radius parameter of R = 0. 4. Transverse energy E T ≥ 30 GeV and log(E HCAL /E ECAL ) ≥ 1. 2 are required. Jets are required to have no tracks with p T ≥ 2 GeV within ∆R = 0. 2 of the jet axis. Finally, jets are required to pass a BIB removal algorithm that relies on calorimeter cell timing and position. Muons from BIB enter the HCAL and can radiate a bremsstrahlung photon, generating an energy deposit that may be reconstructed as a jet with characteristics similar to the hadronic-DPJ. The algorithm identifies events as containing BIB if the triggering jet has at least four HCAL cells at the same φ and in the same layer with timing consistent with that of a BIB energy deposit. 

Two DPJs satisfying the selection criteria described in Section 5 are required in the events selected by the triggers. If more than two DPJs are reconstructed, the one with the highest transverse momentum, labelled the leading DPJ, and the one farthest in ∆φ from the leading one, labelled the subleading DPJ, are used to classify the event. More than two DPJs are found in 9% of the events in the signal MC sample H → 2γ d + X with m H = 125 GeV. Events are classified as one of the three following channels:

• µDPJ-µDPJ,

• µDPJ-hDPJ,

• hDPJ-hDPJ. 

In the µDPJ-hDPJ channel, either the µDPJ or the hDPJ may be the leading DPJ.
section_name: Trigger and event selection

section_content: A data-driven ABCD method is used to estimate the multi-jet background in each of the three channels. The ABCD method uses two nearly uncorrelated variables defined at the event level to create a two-dimensional plane that is split into four parts: region A, where most signal events are located, and three control regions (B, C, and D) that contain mostly background. The number of background events in A can be predicted from the population of the other three regions: N A = N B × N D /N C, assuming negligible leakage of signal into regions B, C and D. For each channel, the ABCD calculation is performed in two separate regions: one background-dominated validation region (VR) to test the validity of the method, and one signal region (SR). The SRs are defined by the selection criteria described in Sections 5 and 6. These define also the VRs except for the BDT cuts. The VRs BDT cuts for the leading and the subleading DPJs are chosen to have negligible signal contamination, which otherwise can bias the ABCD method validation. SR and VR definitions are summarised in Table 2. 

The two event-level variables used to define the ABCD plane are the isolation of the DPJs relative to tracks in the inner detector and the opening angle between the two DPJs in the transverse plane (|∆φ|). Displaced DPJs are expected to be highly isolated in the ID. The track isolation ( p T ) is defined as the scalar sum of the transverse momenta of the tracks reconstructed in the ID and matched to the PV of the event within a ∆R = 0. 4 cone around the DPJ direction. Matching to the PV helps reduce the dependence of p T on the amount of pile-up. The larger of the two p T values, max(Σp T ), is used as the event-level variable. For signal, the opening angle |∆φ| is expected to be large, due to production of the DPJs in the two-body decay of a Higgs boson generated at rest in the transverse plane. 

The ABCD method relies on there being only one source of background, or multiple sources that have identical distributions in the ABCD plane. Studies using the BIB dataset show that BIB events have a narrow |∆φ| distribution, so the requirement |∆φ| > 0. 1 in the ABCD plane removes BIB events that would otherwise contaminate the method for the hDPJ-hDPJ channel, and has no effect on the signal in each channel. After the final selection, the contribution of BIB events to the signal region is negligible. 

The ABCD plane is defined for all the three channels by 0 ≤ max(Σp T ) ≤ 20 GeV and 0. 1 ≤ |∆φ| ≤ π. The region A is defined by max(Σp T ) < 4. 5 GeV and |∆φ| > 0. 625. Regions B, C, and D are defined by reversing one or both of the requirements: max(Σp T ) > 4. 5 GeV and |∆φ| > 0. 625, max(Σp T ) > 4. 5 GeV and |∆φ| < 0. 625, max(Σp T ) < 4. 5 GeV and |∆φ| < 0. 625 respectively. In order to estimate the residual cosmic-ray background component in the ABCD plane, the event selection is applied to the cosmic dataset, and the resulting event yield is multiplied by F CR. The expected number of cosmic-ray events in the validation regions is: 4 ± 3 in region A and 2 ± 2 in region D for the µDPJ-µDPJ channel; 10 ± 5 in region D for the µDPJ-hDPJ channel. In the signal regions, the expected number of cosmic-ray events is: 8 ± 4 in region A and 2 ± 2 in region D for both the µDPJ-µDPJ and the µDPJ-hDPJ channel; 2 ± 2 in region A for the hDPJ-hDPJ channel. No events are observed in the remaining regions in the cosmic dataset. The estimated cosmic-ray event yields are subtracted from each of the ABCD regions before using the method to estimate the multi-jet background yield. 

Other potential backgrounds to the signal include all the processes that lead to real prompt muons and muons plus jets in the final state, such as the SM production of W+jets, Z+jets, t t, single-top-quark, WW, WZ, and ZZ events. MC samples are used to study these processes. They give no contribution after the trigger selection and the definition of muonic-DPJ and hadronic-DPJ and do not enter in the ABCD plane. 

The linear correlation coefficient between the max(Σp T ) and |∆φ| variables in the VRs is verified to be less than 6%, and the signal contamination is verified to be less than 5% for all channels. Table 3 shows the event counts in each of the four regions of the ABCD plane in the validation regions and the expected number of background events in the validation region A in data. Only statistical uncertainties are shown. The expected contribution from cosmic rays is included in all regions and in the background estimation. 

The observed number of events in the validation region A is in agreement with the number predicted by the ABCD method within the statistical uncertainties. 

As additional validation of the ABCD method, control region D of the SR is divided into four subregions. The subregion with low max(Σp T ) and high |∆φ| is treated as a mock signal region, with the other subregions serving as control regions. Applying the method, the expected and the observed numbers of events in the mock signal region are: 231 ± 58 and 184 for the µDPJ-µDPJ channel, 131 ± 41 and 145 for the µDPJ-hDPJ channel, 402 ± 77 and 479 for the hDPJ-hDPJ channel. These are in agreement within the statistical uncertainties. 

Figure 3 shows the distribution of events in the ABCD plane of the µDPJ-µDPJ channel in the SR for the collision data and the MC signal H → 2γ d + X with m H = 125 GeV, assuming a 10 % Higgs boson decay In order to take into account the small signal contamination in regions B, C and D, a likelihood-based ABCD method is used for the background estimation in the SR. It estimates the background in region A by performing a fit to the background and signal yields in the four regions. A likelihood function is formed from the product of four Poisson functions, one for each of the A, B, C, and D regions, describing signal and background expectations. The likelihood takes the form:

where n A, n B, n C and n D are the four observables that denote the number of events observed in each region in data. The N i are linear combinations of the signal and background expectation in each region, defined as follows:

where s (b) is the signal (background) yield in region A, ε i is the signal contamination derived from MC simulation, and τ B and τ C are the nuisance parameters that describe the ratio of the background expectation in the control region to the background expectation in the signal region. The s, b, τ B and τ C values are allowed to float in the fit to the four data regions. The systematic uncertainty associated with the background estimation reported above is described in Section 8.
section_name: Multi-jet background estimation

section_content: The following effects are considered as possible sources of systematic uncertainty in the signal.
section_name: Systematic uncertainties

section_content: The uncertainty in the combined 2015-2016 integrated luminosity is 2. 1% [94], obtained using the LUCID-2 detector [95] for the primary luminosity measurements.
section_name: Luminosity

section_content: The systematic uncertainty in the narrow-scan trigger efficiency is evaluated using a tag-and-probe method applied to J/ψ → µµ events in data and simulation. The difference between the trigger efficiency in data and that in simulation is evaluated as a function of the opening angle between the two muons. The difference in the region ∆R < 0. 05, corresponding to the ∆R expected for signal, is taken as the uncertainty and is 6%. The systematic uncertainty in the tri-muon MS-only trigger efficiency is 5. 8%, taken from the analysis of 2012 data [11] since the algorithm has not undergone a major change since then. The systematic uncertainty in the CalRatio trigger efficiency is taken from Ref. [96] and is 2%.
section_name: Trigger

section_content: The systematic uncertainty in the MC modelling of the input variables used for the BDT training is evaluated for both the µDPJ and the hDPJ. For the µDPJ, the ratio of the distribution in data to that in simulation is computed for each of the BDT input variables using samples of Z → µµ events. Muons from Z boson decay are reconstructed using information only from the MS. The BDT is retrained using MC signal variables scaled to the data using these ratios, and the fit procedure is repeated. The resulting change in the final signal yield is taken as the systematic uncertainty, and its value is 3%. The same procedure is used for the hDPJ, where the ratios of data to simulated distributions are computed from data and MC samples of dijet events. The resulting uncertainty is 14%.
section_name: BDT shape

section_content: The systematic uncertainty in the single-γ d reconstruction efficiency is evaluated using a tag-and-probe method applied to J/ψ → µµ events in 2015 data and simulation. J/ψ → µµ decays are selected, and the efficiency is evaluated as a function of the opening angle ∆R between the two muons, for both the data and simulated J/ψ decays. For low ∆R values, the efficiency decreases due to the difficulty of reconstructing two tracks with small angular separation in the MS. The difference in J/φ → µµ reconstruction efficiency between simulation and data in the ∆R interval between 0 and 0. 06 (where the DPJ samples are concentrated) amounts to 15%, and is taken as the uncertainty.
section_name: Muon reconstruction

section_content: The jet energy scale and jet energy resolution introduce uncertainties in the signal yield of 1% to 8% and 1% to 5% respectively, depending on the signal process, where the processes with two dark photons are less affected. These uncertainties are calculated using the procedure detailed in Ref. [97]. Since the jets used in this analysis are required to have a low fraction of energy in the electromagnetic calorimeter, additional jet energy uncertainties are derived as a function of electromagnetic energy fraction as well as of pseudorapidity. These additional jet energy uncertainties are found to have an effect of up to 4% on the signal yield, and are taken in quadrature with the regular jet energy uncertainties.
section_name: Jet energy scale and jet energy resolution

section_content: The presence of multiple collisions per bunch crossing affects the efficiency of the ID track isolation criterion quantified in terms of Σp T. The systematic uncertainty is evaluated by comparing Σp T for muons from a sample of reconstructed Z → µµ events in data with that in simulation, as a function of the number of interaction vertices in the event. The systematic uncertainty is evaluated as the maximum difference at the value of the selection requirement on max(Σp T ). It is found to be 5. 1%.
section_name: Effect of pile-up on Σ p T

section_content: The observed numbers of events in the ABCD regions and the expected number of background events in the signal region A are summarised in Table 4. The expected number of background events in region A is estimated using the likelihood-based ABCD method, assuming no signal and not including the observed data in region A. The background estimate includes both the multi-jet and cosmic-ray background, where the former is obtained as described in Section 7, and the latter is estimated from the cosmic dataset. Both sources are included in the expected background given in Table 4. The observed number of events in region A is in agreement with the predicted number of background events. Table 5 shows the expected number of signal events in region A for the FRVZ model parameters of Table 1 and the following assumptions: a value of B(H → f d 2 f d 2 ) = 10% for a Higgs boson with m H = 125 GeV, which is not excluded by the current measurements [5] ; a value of B(H → f d 2 f d 2 ) = 100% and a production cross section of σ = 5 pb for a BSM scalar boson with m H = 800 GeV. 

Upper limits on the production cross section times branching fraction (σ × B) as a function of the γ d proper decay length are derived for the FRVZ H → 2γ d + X and H → 4γ d + X processes using the CL s method [98]. Since each signal sample was generated for a particular proper decay length, it is necessary to extrapolate the signal efficiency to other decay lengths to obtain limits as a function of cτ. This is achieved by applying to the i-th dark photon in the event a weight

where τ ref is the lifetime with which the event sample was simulated, τ new is the lifetime for which it is weighted, and t i is the proper decay time of the i-th dark photon. Each event is weighted by the product of the individual dark-photon weights. The weighted sample is used to evaluate the signal efficiency for τ new. Figure 4 shows the extrapolated signal efficiency for the H → 2γ d + X and H → 4γ d + X processes as a function of cτ of the dark photon in the µDPJ-µDPJ, µDPJ-hDPJ and hDPJ-hDPJ channels. The tri-muon MS-only trigger has a lower efficiency for the H → 2γ d + X process with a m H = 125 GeV Higgs boson than for the other processes, resulting in a lower signal efficiency in the µDPJ-µDPJ channel. The p T requirements of the CalRatio trigger are not optimal for selecting jets produced by γ d decays in the processes with a m H = 125 GeV Higgs boson, resulting in a signal efficiency below 1% in the hDPJ-hDPJ channel. The muon narrow-scan trigger helps to recover some efficiency in the µDPJ-hDPJ channel for these processes. 

The observed 95% CL cross-section upper limits in the µDPJ-µDPJ channel for the H → 2γ d + X and H → 4γ d + X processes are presented in Figure 5 for m H = 125 GeV. The 95% CL exclusion limits in the µDPJ-µDPJ and hDPJ-hDPJ channels for the process H → 2γ d + X are presented in Figure 6 for m H = 800 GeV. The figures also show the expected limits obtained from the likelihood-based ABCD method, using the background estimate derived from the background-only fit using data in the four regions. Excluded cτ ranges are summarised in Table 6, assuming B(H → f boson decay branching fractions into γ d of 10% from a search for displaced dark-photon jets [11] and prompt dark-photon jets [14] at ATLAS. The search of Ref [11], which explored the same region probed by this analysis, is slightly more sensitive in the region of high γ d mass and low. This is due to inclusion of dark-photon jets with both muon and hadron constituents, which are not used in the current analysis. The search of Ref. [14] excluded high values (shorter lifetimes), a region complementary to this analysis. [11] (black line) and prompt [14] (red line) dark-photon jets searches.
section_name: Results and interpretation

section_name: Excluded cτ

section_content: The ATLAS detector at the LHC is used to search for the production of displaced dark-photon jets in a sample of pp collisions at √ s = 13 TeV corresponding to an integrated luminosity of 36. 1 fb -1. No significant excess of events compared with the background expectation is observed, and 95% confidence-level upper limits are set on the production cross section times branching fraction of scalar bosons that decay into dark photons according to the FRVZ model. The upper limits are computed as a function of the proper decay length cτ of the dark photon γ d. In addition to the increase in integrated luminosity and centre-of-mass energy, improvements in background suppression and the exploitation of hadronic γ d decays result in increased sensitivity compared with the ATLAS search using 8 TeV pp data. 

The results for H → 2γ d + X, when H is the Higgs boson, are also interpreted as 90% confidence-level limits on the kinetic mixing parameter as a function of the dark-photon mass. These results improve upon the constraints set in previous LHC searches.
section_name: Conclusions

section_name: The ATLAS Collaboration

section_name: Acknowledgements

section_content: We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia ; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil ; NSERC, NRC and CFI, Canada ; CERN; CONICYT, Chile; CAS, MOST and NSFC, China ; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic ; DNRF and DNSRC, Denmark ; IN2P3-CNRS, CEA-DRF/IRFU, France; SRNSFG, Georgia; BMBF, HGF, and MPG, Germany ; GSRT, Greece; RGC, Hong Kong SAR, China ; ISF and Benoziyo Center, Israel ; INFN, Italy ; MEXT and JSPS, Japan ; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden ; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan ; TAEK, Turkey; STFC, United Kingdom ; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada ; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union ; Investissements d' Avenir Labex and Idex, ANR, France ; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain ; The Royal Society and Leverhulme Trust, United Kingdom. 

The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [99].
section_name: Acknowledgements