Search for dark matter in events with heavy quarks and missing transverse momentum in pp collisions with the ATLAS detector

This article reports on a search for dark matter pair production in association with bottom or top quarks in 20 . 3 fb − 1 of pp collisions collected at √ s = 8 TeV by the ATLAS detector at the LHC. Events with large missing transverse momentum are selected when produced in association with high-momentum jets of which one or more are identiﬁed as jets containing b -quarks. Final states with top quarks are selected by requiring a high jet multiplicity and in some cases a single lepton. The data are found to be consistent with the Standard Model expectations and limits are set on the mass scale of effective ﬁeld theories that describe scalar and tensor interactions between dark matter and Standard Model particles. Limits on the dark-matter–nucleon cross-section for spin-independent and spin-dependent interactions are also provided. These limits are particularly strong for low-mass dark matter. Using a simpliﬁed model, constraints are set on the mass of dark matter and of a coloured mediator suitable to explain a possible signal of annihilating dark matter.


Introduction
The existence of dark matter (DM) in the Universe is highly motivated by many astrophysical and cosmological observations [1][2][3][4].However, its nature remains a mystery.One of the best motivated candidates for a DM particle is a weakly interacting massive particle (WIMP) [5].At the Large Hadron Collider (LHC), one can search for DM particles (χ) that are pair produced in pp collisions.These studies are sensitive to low DM masses (m χ ≤ 10 GeV), and therefore provide information complementary to direct DM searches, which are most sensitive to larger DM masses [6][7][8][9].
If the particles that mediate the interactions between DM and Standard Model (SM) particles are too heavy to be produced directly in the experiment, their interactions can be described by contact operators in the framework of an effective field theory [10][11][12].For each operator considered, the reach is expressed in terms of the effective mass scale of the interaction, M * , and of the χ-nucleon cross-section, σ χ-N , as a function of m χ .
Since DM particles do not interact in the detector, the main signature of DM pair production at colliders is large missing transverse momentum.Initial-state radiation (ISR) of jets, photons, Z, or W bosons, was used to tag DM pair production at colliders in several searches at the Tevatron [13] and the LHC [14][15][16][17][18][19][20][21].
A new search for DM pair production in association with one b-quark or a pair of heavy quarks (b or t) was proposed in Ref. [22].The dominant Feynman diagrams for these processes are shown in Fig. 1.To search for these processes, dedicated selections are defined to reconstruct the various production and decay modes of these heavy-quark final states.For final states containing a semileptonic decay of a top quark, the results of the search for a supersymmetric partner of the top quark are used [23].
The analysis presented in this article is particularly sensitive to effective scalar interactions between DM and quarks described by the operator [12] O scalar = q m q M N * qq χχ, (1) where N = 3 for Dirac DM (D1 operator) and N = 2 for complex scalar DM (C1 operator).The quark and DM fields are denoted by q and χ, respectively.The scalar operators are normalized by m q , which mitigates contributions to flavour-changing processes, strongly constrained by flavour physics observables [24,25], through the framework of minimal flavour violation (MFV).The dependence on the quark mass makes final states with bottom and top quarks the most sensitive to these operators.This search is also sensitive to tensor couplings between DM and quarks.The tensor operator (D9), which describes a magnetic moment coupling, is parameterized as [12]: MFV suggests that the D9 operator should have a mass dependence from Yukawa couplings although canonically this is not parametrised as such.
The results are also interpreted in light of a bottom-Flavoured Dark Matter model (b-FDM) [26].The b-FDM model was proposed to explain the excess of gamma rays from the galactic centre, recently observed by the Fermi Gamma-ray Space Telescope, and interpreted as a signal for DM annihilation [27].This analysis of the data recorded by the Fermi-LAT collaboration favours DM with a mass of approximately 35 GeV annihilating into b-quarks via a coloured mediator.In this model, a new scalar field, φ, mediates the interactions between DM and quarks as shown in Fig. 2. DM is assumed to be a Dirac fermion that couples to right-handed, down-type quarks.The lightest DM particle, which constitutes cosmic DM, preferentially couples to b-quarks.The collider signature of this model is b-quarks produced in association with missing transverse momentum.This analysis sets constraints on the mass of the mediator and DM particle in the framework of the b-FDM model.

Detector description and physics objects
The ATLAS detector [33] at the LHC covers the pseudorapidity1 range of |η| < 4.9 and is hermetic in azimuth φ.It consists of an inner tracking detector surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters, and an external muon spectrometer incorporating large superconducting toroidal magnets.A three-level trigger system is used to select events for subsequent offline analysis.The data set used in this analysis consists of 20.3 fb −1 of pp collision data recorded at a centre-of-mass energy of √ s = 8 TeV with stable beam conditions [34] during the 2012 LHC run.All subsystems listed above were required to be operational.
This analysis requires the reconstruction of muons, electrons, jets, and missing transverse momentum.Muon candidates are identified from tracks that are well reconstructed inside both the inner detector and the muon spectrometer [35].To reject cosmic-ray muons, muon candidates are required to be consistent with production at the primary vertex, defined as the vertex with the highest Σ(p track T ) 2 , where p track T refers to the transverse momentum of each track.Electrons are identified as tracks that are matched to a well-reconstructed cluster in the electromagnetic calorimeter.Electron candidates must satisfy the tight electron shower shape and track selection criteria of Ref. [36].Both electrons and muons are required to have transverse momenta p T > 20 GeV and |η| < 2.5.Potential ambiguities between overlapping candidate objects are resolved based on their angular separation.If an electron candidate and a jet overlap within ∆R < 0.2, then the object is considered to be an electron and the jet is discarded.If an electron candidate and any jet overlap within 0.2 < ∆R < 0.4, or if an electron candidate and a T with ATLAS 3 Table 1 Selections for signal regions 1-4.Variables p j i T (p b i T ) represent the transverse momentum of the i-th jet (b-tagged jet).The asymmetric transverse mass am T2 [28][29][30], topness [31], m jjj and Razor R [32] GeV b-tagged jet overlap within ∆R < 0.2 of each other, then the electron is discarded and the jet is retained.Photon candidates must satisfy the tight quality criteria and |η| < 2.37 [37].
Jet candidates are reconstructed using the anti-k t clustering algorithm [38] with a radius parameter of 0.4.The inputs to this algorithm are three-dimensional topological clusters [39].The four-momentum of the jet is defined as the sum of the four-momenta of the topological clusters, assuming that each cluster originates from a particle defined to be massless and to come from the interaction point.
To calibrate the reconstructed energy, jets are corrected for the effects of calorimeter response and inhomogenities using energy-and η-dependent calibration factors based on simulation and validated with extensive test-beam and collision-data studies [40].In the simulation, this procedure calibrates the jet energies to those of the corresponding jets constructed from stable simulated particles.In-situ measurements are used to further correct the data to match the energy scale in simulated events.Effects due to additional pp interactions in the same and preceding bunch crossings (pile-up effects) are corrected [41].
Jets containing particles from the hadronisation of a b-quark (b-jets) are tagged using a multivariate algorithm [42,43].The b-tagging algorithm combines the measurement of several quantities distinguishing heavy quarks from light quarks based on their longer lifetime and heavier mass.These quantities include the distance of closest approach of tracks in the jet to the primary event vertex, the number and position of secondary vertices formed by tracks within the jet, as well as the invariant mass associated with such vertices.The algorithm is trained on Monte Carlo (MC) simulations and its performance is calibrated using data.To optimize the sensitivity of this analysis, a requirement on the output of the b-tagging algorithm which provides a 60% (70%) b-jet efficiency operating point is used in signal regions (SR) 1 and 2 (3 and 4) defined below.The corresponding misidentification probability is 15% (20%) for c-jets, and less than 1% for light-quark jets.The aforementioned b-tagging efficiencies and misidentification probabilities were derived in a simulated t t sample with jet transverse momenta of p T > 20 GeV and |η| < 2.5.
The missing transverse momentum, with magnitude E miss T , is defined as the negative vector sum of the transverse momenta of jets, muons, electrons, photons, and topological clusters not assigned to any reconstructed objects [44].

Event selection
Candidate signal events containing at least one high-p T jet and large E miss T are assigned to one of four orthogonal signal regions.The first two signal regions focus on events with DM produced in conjunction with one (SR1) or two (SR2) b-quarks in the final state.SR3 and SR4 target events in which DM is produced in conjunction with a tt pair, where either both top quarks decay hadronically (SR3) or one top quark decays hadronically and the other semileptonically (SR4).SR4 was developed for a top squark search by the ATLAS Collaboration and coincides with the "tNbC mix" signal region described in Ref. [23].
Table 2 The quoted uncertainties include all statistical and systematic effects added in quadrature.The effective mass scale, M * , is set to be 100/40/600 GeV for the D1/C1/D9 operators, approximately corresponding to the expected limit.The probabilities of the background-only hypothesis, p-values, are also given.The last two lines show the observed and expected 95% CL upper limits on the number of beyond-the-SM events.Events assigned to SR1 and SR2 are required to pass a calorimeter-based E miss T trigger with a threshold of 80 GeV.To enrich the sample in pp → χ χ + b( b), events are required to have a low jet multiplicity (n jets <5), E miss T > 300 GeV, and the most energetic b-tagged jet must have a p T > 100 GeV.The azimuthal separation between the directions of the jets and the missing transverse momentum is required to be more than 1.0 radian.Events with at least one identified muon or electron are discarded to reject leptonic decays of W and Z bosons.Events satisfying these selection criteria are assigned to SR1 provided that the jet multiplicity does not exceed two.Events are assigned to SR2 when at least three jets are reconstructed in the event and the second most energetic jet has p T > 100 GeV.If there is a second b-tagged jet it has to satisfy p T > 60 GeV.
Events assigned to SR3 are required to pass triggers specifically designed to select hadronic decays of top quark pairs.Such triggers require either five jets with p T ≥ 55 GeV each or four jets with p T ≥ 45 GeV, of which one is tagged as a b-jet.To select pp → χ χ + t t events, at least five reconstructed jets are required, of which at least two are b-tagged, and E miss T > 200 GeV.Furthermore, the azimuthal separation between the most energetic b-jet and the missing transverse momentum is required to be at least 1.6 radians.To reduce W/Z leptonic decays and leptonic top quark decays, events with at least one identified muon or electron are discarded.To maximize the rejection of the abundant t t background, the Razor variable R [32] is used.This variable utilizes both transverse and longitudinal information about the event to fully exploit the kinematics of the decay.To separate signal and background, R > 0.75 is required.
To enrich the sample in pp → χ χ + t t with one semileptonic decay of the t quark, events assigned to SR4 use single-lepton or E miss

√
GeV, with H 4j T = 4 i=1 p T (jet i ) and where the jets are ordered by decreasing p T .The azimuthal angle between the missing transverse momentum and the two most energetic jets is required to be greater than 0.6 radians.
Special variables, such as the asymmetric transverse mass am T2 [28][29][30] and the topness variable [31], are used to reject the dileptonic t t component of the background.Details can be found in Ref. [23].The diboson background is suppressed by a requirement on the three-jet invariant mass (m jjj < 360 GeV) [23].A τ veto rejects t t events with hadronically decaying τ leptons in the final state.Additional selection criteria [23] on the angles between the lepton and the various jets 2 Since the longitudinal component of the momentum of the neutrinos is not measured, the measured properties of the W boson candidates are limited to their transverse momentum and transverse mass, defined as where E miss T is the magnitude of the missing transverse momentum vector, p T is the transverse momentum of the lepton and p x and p y (E miss x and E miss y ) are the magnitude of the x and y components of the lepton momentum (missing transverse momentum) respectively.are imposed to further reduce the t t background.Table 1 provides an overview of the selections applied in all four signal regions.
The product of the detector acceptance A and the reconstruction efficiency for the selections described above varies between 0.1% and 8% depending on the signal region, operator, and specific channel considered.SR1 and SR2 have the highest efficiencies (A × > 2%) for the D9 operator, while SR3 and SR4 are most efficient for the D1 and C1 operators (A × > 1%).
The dominant background for SR1 and SR2 is due to Z → νν events produced in conjunction with one or more jets.This irreducible background is estimated from data using two control regions (CRs).The first CR exploits Z+jets events with Z → µ + µ − , while the second uses γ+jets events for which the production at high transverse momentum (p γ T > M Z ) mimics that of Z+jets [45].The γ+jets control region substantially increases the number of events at large missing transverse momentum.The transverse momentum of the dimuon pair or photon is added vectorially to the E miss T of the event to simulate the Z → ν ν background.Corrections to compensate for the differences in efficiency and acceptance between the Z(νν)+jets and Z(µ + µ − )+jets or γ+jets are derived from data using control regions without b-tagged jets before applying any requirements on the missing transverse momentum.Remaining kinematic selections correspond to the ones described in Tab. 1.A muon control region is chosen because the energy loss of muons in the detector is comparatively small.The systematic uncertainties introduced by this data-driven procedure on the Z(ν ν)+jets background are approximately 10%, mainly from the flavour composition of background processes, kinematic differences between the control and signal regions and relative normalizations of backgrounds.
Production of W/Z+jets with subsequent leptonic decays of W and to a much smaller degree Z is also a substantial source of background for SR1 and SR2 when the resulting charged leptons fail to be identified or if the W or Z bosons decay to τ leptons.These contributions are estimated from Z( + − )+jets and W ( ν)+jets MC samples generated using ALPGEN2.3[46] with the CTEQ6L1 [47] parton distribution function (PDF) set.These samples are generated with up to five light partons (u, d, s) and one c quark or two heavy quarks (c, b) per event.W + b production is highly suppressed and therefore negligible.A control region enriched in W ( ν)+jets events is selected by adding a lepton requirement to the selection and is used to validate the estimate of this background.The purity of W ( ν)+jets in the control region for SR1 (SR2) is 67% (47%).After full selection the contribution of b(c)-quarks to the dominant W ( ν)+jets background is approximately 39% (38%) for SR1 and 52% (37%) for SR2.The systematic uncertainty on this background is approximately 20%.Finally, the small contribution from t t is estimated using MC samples and validated in data control regions before applying signal selection requirements.The t t process is selected with very high purity by requiring events with one lepton and large jet multiplicities.
The dominant source of background for SR3 and SR4 is t t events.In SR3, this contribution is estimated from data using a control region not overlapping with SR4 and largely dominated by t t events with one of the two top quarks decaying semileptonically.The fivejets requirement is relaxed to three jets.Additionally, the event is required to contain exactly one lepton with p e(µ) T > 30 (25) GeV and must fullfil E miss T + m T > 25 (30) GeV for the electron (muon) channel.The potential signal contribution to this selection is less than 0.1%.The uncertainties are small because the SR3 data control region uses a kinematic region similar to the signal region with the lepton veto and jet multiplicity being the main difference.These effects were studied and considered as systematic uncertainties.Dominant uncertainties are related to jets and the top quark momentum distribution.Corrections to compensate for the differences in efficiency and acceptance between hadronic and semileptonic top decays are derived from MC samples generated using the POWHEG BOX generator [48] interfaced with JIMMY4.31[49] with the next-to-leadingorder (NLO) PDF set CT10 [50].The systematic uncertainty on the t t background in SR3 of approximately 7% is derived by studying corrections for the top quark momentum distribution, and shower modelling by interfacing the same generator with PYTHIA6 [51,52].
In SR4, the t t background is estimated from data using a control region obtained by requiring 60 GeV < m T < 90 GeV and loosening the selection criteria on E miss T , am T2 , and E miss T / H 4j T .A similar selection, but applying an inverted b-tagging requirement, is used to estimate the W ( ν)+jets background.The uncertainty on the t t background is estimated to be approximately 20% [23], which is larger than the uncertainty in SR3 due to the tight selection, while the uncertainty on the W ( ν)+jets background is approximately 7%.These uncertainties are evaluated by varying the renormalisation and factorisation scale of the simulations, comparing alternative PDF sets, and studying the effects of different shower generators and of ISR and final-state radiation.Additional sources of background, which include singletop, t t + Z/W , and diboson production, are estimated in all signal regions using simulations.The single-top (s-channel) and W t background is generated using the POWHEG generator.The single-top t-channel is generated with ACERMC3.8 [53] interfaced with PYTHIA6.Associated production of t t and a vector boson (W , Z) are generated with MADGRAPH5 [54] with up to two additional partons interfaced with PYTHIA6.The cross-sections for t t production in association with a W (Z) boson are determined using the MSTW2008 NLO (CTEQ6.6M)PDF sets.The diboson samples are generated using HERWIG6.520[55,56] and JIMMY4.31 with the CTEQ6L1 PDF set.The multijet background is estimated using data-driven methods [57] and is found to be negligible in all signal regions after full selection.
Object reconstruction efficiencies in simulated events are corrected to reproduce the performance measured in data.The systematic uncertainty of the background estimates derived from simulation combines the uncertainties on the efficiency of the b-tagging algorithm, the uncertainties on the determination of the energy scale and resolution of the jet energy and E miss T , the theoretical uncertainty on the various cross-sections, changes in the shapes of distributions used to extrapolate event counts from control regions to the signal region, data driven corrections and the PDF uncertainties.Overall, the systematic uncertainty on the background estimated from simulation is calculated to be between 12% and 18%, depending on the signal region.
The simulation of the signal samples of pp → χχ + b(b), pp → χχ + tt, and b-FDM employs the MAD-GRAPH5 generator interfaced with PYTHIA6 using the CTEQ6L1 PDF.Samples are generated for operators D1, C1, and D9, assuming M * = 1 TeV and m χ between 10 GeV and 1300 GeV.Samples for the b-FDM model are generated for m χ values between 1 GeV and 1300 GeV and mediator masses, m φ , between 5 GeV and 3000 GeV.The instrumental uncertainties on the simulated signal yields for D1, C1, and D9 operators are between 11% and 15%, depending on the signal region.The equivalent uncertainties for the b-FDM model range between 6%-16% depending on m χ and the mediator mass.The uncertainties from the PDF are computed by comparing the rates obtained with the default PDF set (CTEQ6L1) with those obtained with two alternative sets (MSTW2008LO and NNPDF21LO [58,59]).The uncertainties on the signal acceptance from PDF and scale variations are estimated to be approximately 10% for the D1, C1, and D9 operators for m χ =10 GeV and approximately 6% for b-FDM models.
The validity of the effective field theory assumption depends on the momentum transfer of the process modelled, which should be below the energy scale of the underlying interactions [60].To account for this, the momentum transfer m(χχ) = Q tr in the events is required to be less than the energy scale probed.Specifically, Q tr must be smaller than the mass M of the heavy mediator.For an ultraviolet completion this implies M * = M/ √ g q g χ .Along with perturbativity of the couplings g q g χ < 4π this leads to the following validity requirements on MC truth level:

Results
Table 2 shows the expected background from various sources in the four signal regions as well as the observed yields in data.The expected signal yields for the operators D1, C1, and D9, as well as for the b-FDM model are also shown.The probabilities of the backgroundonly hypothesis, p-values, for the signal regions SR1, SR2, SR3, and SR4 are 0.09, 0.29, 0.24, and 0.18, respectively.As no significant excess is observed, limits on the signal yield are set using a profile likelihood ratio test following the CL s prescription [61].Also given is the 95% confidence level (CL) upper limit on the number of beyond-the-SM events.The yields for the b-FDM model are obtained assuming m χ =10 GeV and a mediator mass m φ = 600 GeV.The limit on M * for a given assumption on m χ is determined by varying M * and scaling the number of signal events predicted by the corresponding sample generated with M * = 1 TeV until it is equal to the observed upper limit on beyond-the-SM events.The corresponding production cross-section for DM produced via the D1 operator in association with b(t)-quarks and m χ = 10 GeV is 38 (221) fb.The cross-section for b-FDM models with m φ = 600 GeV and m χ = 10 GeV is 134 fb.The signal efficiency is independent of M * .Figure 3 shows the E miss T distributions for (a) SR1, (b) SR2, and (d) SR4 and (c) the R variable for SR3.
Figure 4 shows the 90% CL exclusion curves for the effective mass scale M * as a function of m χ .The results for the operators D1, C1, and D9 are presented individually for all four signal regions.The best limits on the D1 and C1 operators are obtained using SR4, while SR1 provides the best limits on the D9 operator, as shown in Fig. 4.These limits are then converted into limits on the χ-nucleon cross-section [12].Figures 5  and 6 show the corresponding 90% CL exclusion curves for the spin-independent and spin-dependent χ-nucleon cross-section for the scalar (D1) and tensor (D9) operators as a function of m χ for the strongest results obtained in any signal region.The most stringent limits set by direct detection experiments [6][7][8][9] are also shown.Only m χ where more then 90% of the events fulfill the effective field theory validity constraints are shown in Figs. 5 and 6.

LUX (2013)
Fig. 5 Upper limits at 90% CL on the spin-independent χnucleon cross-section (σ SI χ-N ) for the scalar operator D1 (red) as a function of mχ.The yellow and green curves represent the exclusion limits recently set by the LUX and Super-CDMS collaborations [6,7,62].The coupling is assumed to be gqgχ = g = 4π.
sults reported in this article represent the first ATLAS limits on the scalar operator C1 and they significantly improve the sensitivity to χ-nucleon interactions mediated by the scalar operator D1 compared to previous ATLAS results [14,16,18,19].
Figure 7 shows the exclusion curves observed and expected for the b-FDM model as a function of the mediator and DM masses.For each point in (m χ , m φ ), the signal region with the best expected sensitivity is used.For a DM particle of approximately 35 GeV, as suggested by the interpretation of data recorded by the Fermi-LAT collaboration, mediator masses between approximately 300 GeV and 500 GeV are excluded at 95% CL. [GeV]

PICASSO (2012)
Fig. 6 Upper limits at 90% CL on the spin-dependent χnucleon cross-section (σ SD χ-N ) for the tensor operator D9 (red) as a function of mχ.The yellow and green curves represent the exclusion limits recently set by the COUPP and PICASSO collaborations [8,9,62].The coupling is assumed to be gqgχ = g = 4π.7 Exclusion contour at 95% CL for the b-FDM model from combined results of SR1 and SR2.The expected limit is given by the dashed line, and the yellow band indicates the ±1 σ uncertainty.The observed limit is given by the solid red line.The region beneath the curve indicating the observed limit is excluded.

Conclusions
In summary, this article reports a search for dark-matter pair production in association with bottom or top quarks.The analysis is performed using 20.3 fb −1 of pp collisions collected at √ s = 8 TeV by the ATLAS de-tector at the LHC.The results are interpreted in the framework of an effective field theory to set stringent limits on scalar and tensor interactions between Standard Model and DM particles.The data are found to be consistent with the Standard Model expectations, and limits are set on the mass scale of effective field theories that describe scalar and tensor interactions between DM and Standard Model particles.The exclusion limits are strongest at low DM masses.The limit on the χ-nucleon cross-section mediated by the D1 operator is improved significantly with respect to previously published ATLAS results by obtaining sensitivities of approximately σ SI χ-N = 10 −42 cm 2 for m χ = 10 GeV.Constraints on b-Flavoured Dark Matter models, suitable to explain a possible signal of annihilating DM, are also presented.The excluded regions depend on m χ and m φ .For m χ = 35 GeV, mediator particles with m φ = 300-500 GeV are excluded.

Fig. 1 Fig. 2
Fig. 1 Dominant Feynman diagrams for DM production in conjunction with (a) a single b-quark and (b) a heavy quark (bottom or top) pair using an effective field theory approach.

T
triggers, and require exactly one isolated lepton (electron or muon) with p T > 25 GeV, at least four high-p T jets, where one jet is b-tagged with p T > 60 GeV.Events with E miss T > 270 GeV are selected when the transverse mass 2 formed by the lepton and E miss T , m T ( , E miss T ), exceeds 130 GeV and E miss T / H 4j T > 9

Fig. 3
Fig. 3 Comparison between data and expected SM background.(a), (b): E miss T variable for SR1 and SR2 and for an example signal with the operator D9. (c): R variable for SR3 excluding the selection on R and for an example signal with the operator D1. (d): E miss T variable for SR4 excluding the selection on E miss T and for an example signal with the operator D1.Other backgrounds are composed of diboson and multijet production.The expected signal for χ χ + b( b) (SR1, 2) and for χ χ + t t (SR3, 4) production for mχ = 10 GeV is given by the red line.The final selection requirements are indicated by an arrow.The error bars represent the statistical uncertainty.The dashed area shows the systematic uncertainty on the background estimation.Events with values exceeding the range presented are included in the highest bin.

Fig. 4
Fig.4 Lower limits on M * at 90% CL for the SR1 (red), SR2 (black), SR3 (green), and SR4 (blue) as a function of mχ for the operators (a) D1, (b) C1, and (c) D9.Solid lines and markers indicate the validity range of the effective field theory assuming couplings gqgχ < 4π, the dashed lines and hollow makers represent the full collider constraints.
have recently claimed possible observations of DM.The re-

Fig.
Fig.7Exclusion contour at 95% CL for the b-FDM model from combined results of SR1 and SR2.The expected limit is given by the dashed line, and the yellow band indicates the ±1 σ uncertainty.The observed limit is given by the solid red line.The region beneath the curve indicating the observed limit is excluded.