Search for pair-produced higgsinos decaying via Higgs or Z bosons to final states containing a pair of photons and a pair of b-jets with the ATLAS detector

A search is presented for the pair production of higgsinos ̃𝜒 in gauge-mediated supersymmetry models, where the lightest neutralinos ̃𝜒 01 decay into a light gravitino ̃𝐺 in association with either a Higgs ℎ or a 𝑍 boson. The search is performed with the ATLAS detector at the Large Hadron Collider using 139 fb −1 of proton–proton collisions at a centre-of-mass energy of √ 𝑠 = 13 TeV . It targets ﬁnal states in which a Higgs boson decays into a photon pair, while the other Higgs or 𝑍 boson decays into a 𝑏̄𝑏 pair, with missing transverse momentum associated with the two gravitinos. Search regions dependent on the amount of missing transverse momentum are deﬁned by the requirements that the diphoton mass should be consistent with the mass of the Higgs boson, and the 𝑏̄𝑏 mass with the mass of the Higgs or 𝑍 boson. The main backgrounds are estimated with data-driven methods using the sidebands of the diphoton mass distribution. No excesses beyond Standard Model expectations are observed and higgsinos with masses up to 320 GeV are excluded, assuming a branching fraction of 100% for ̃𝜒 01 → ℎ ̃𝐺 . This analysis excludes higgsinos with masses of 130 GeV for branching fractions to ℎ ̃𝐺 as low as 36%, thus providing complementarity to previous ATLAS searches in ﬁnal states with multiple leptons

In the Minimal Supersymmetric Standard Model [1,2] two neutral higgsinos, the proposed spin-half supersymmetric (SUSY) [3] partners of Higgs bosons, mix together with the spin-half partners of the neutral SU(2)  and U(1)  gauge bosons, forming four neutralino states χ0 1,2,3,4 , while the charged higgsino mixes with the spin-half partners of the  ± to form two pairs of chargino states χ± 1,2 , where indices are ordered in increasing mass.When the mass terms for the spin-half partners of the gauge bosons are large, the most accessible electroweak states at colliders are the lightest, higgsino-dominated, neutralinos χ0 1 , χ0 2 and charginos χ± 1 .These particles are expected to have a mass around the electroweak scale following naturalness arguments [4,5], with mass splittings of 1 GeV or less.Such mass spectra are predicted, for example, in general gauge mediation (GGM) models [6][7][8][9][10] or in gauge-mediated supersymmetry breaking (GMSB) models [6][7][8] in which the lightest neutralino χ0 1 is higgsino dominated and decays into a light gravitino ( G) [11], the proposed spin- 3  2 partner of the graviton, in association with a Higgs boson ℎ or  boson.The heavier χ0 2 and χ± 1 decay into the lighter χ0 1 in association with quarks and leptons.
These processes represent a promising supersymmetry search target at the LHC, particularly if other SUSY particles are sufficiently heavy as to be kinematically inaccessible or to have suppressed cross-sections.The resulting gravitinos would not be detected by collider experiments, but their presence could be inferred from momentum imbalance in the direction perpendicular to the beams ( miss T ).The branching fractions B ( χ0 1 → ℎ G) and B ( χ0 1 →  G) will depend on the components of the neutralino mixing matrix.In the following, it is assumed that there are no additional decay modes of the lightest neutralino, i.e., B ( The subsequent decay products of the ℎ and the  can be detected in a variety of different decay modes.Decays of the lightest neutralino into photons are also possible, but result in significantly different final states, which are considered in other searches, e.g.Ref. [12].This paper presents a search that targets final states in which a Higgs boson decays into a photon pair, while the other Higgs boson decays into a  b pair, with missing transverse momentum associated with the two gravitinos.This final state also has sensitivity to mixed decays in which one neutralino decays into ℎ and one into  where the  decays into a pair of -quarks, as shown in Figure 1.The ATLAS Collaboration has previously performed searches using events with at least three -jets [13], in ℎ →  events without requiring the -jet decay mode of the other boson [14], events with four or more leptons [15], events containing two leptons, jets and missing transverse momentum [16], and in events with two boosted bosons, missing transverse momentum and no leptons [17].The CMS collaboration has performed several searches in a variety of decay channels [18][19][20] including  [20].For the searches involving two Higgs bosons in the final state, decays into  combine the advantages of excellent mass resolution in ℎ →  with the large Higgs boson branching fraction to  b. The analysis is based on data collected by the ATLAS detector [21] using proton-proton ( ) collisions delivered by the LHC [22] between 2015 and 2018, at a centre-of-mass energy of √  = 13 TeV with a bunch-crossing interval of 25 ns.The average number of   interactions per bunch crossing ranged from about 13 in 2015 to about 39 in 2018.After applying trigger and data-quality requirements [23] the integrated luminosity of the data sample, measured with the LUCID-2 detector [24] using the methods described in Ref. [25], corresponds to 139.0 ± 2.4 fb −1 .
The ATLAS detector at the LHC is a multipurpose particle detector with a forward-backward symmetric cylindrical geometry and a near 4 coverage in solid angle. 1 It consists of an inner tracking detector (ID) 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the -axis along the beam pipe.The -axis points from the IP to the centre of the LHC ring, and the -axis points upwards.Polar coordinates (, ) are used in the transverse plane,  being the azimuthal angle around the -axis.The pseudorapidity is defined in terms of the polar angle  as  = − ln tan(/2).Angular distance is measured in units of Δ ≡ √︁ (Δ) 2 + (Δ) 2 .surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic and hadron calorimeters, and a muon spectrometer (MS).The inner tracking detector covers the pseudorapidity range || < 2.5.It consists of silicon pixel, silicon microstrip, and transition radiation tracking detectors.Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic (EM) energy measurements with high granularity.A steel/scintillator-tile hadron calorimeter covers the central pseudorapidity range || < 1.7.The endcap and forward regions are instrumented with LAr calorimeters for both the EM and hadronic energy measurements up to || = 4.9.The muon spectrometer surrounds the calorimeters and is based on three large superconducting air-core toroidal magnets with eight coils each.The field integral of the toroids ranges between 2.0 and 6.0 Tm across most of the detector.The muon spectrometer includes a system of precision tracking chambers and fast detectors for triggering.A two-level trigger system is used to select events.The first-level trigger is implemented in hardware and uses a subset of the detector information to accept events at a rate below 100 kHz.This is followed by a software-based trigger that reduces the accepted event rate to 1 kHz on average, depending on the data-taking conditions.An extensive software suite [26] is used in data simulation, in the reconstruction and analysis of real and simulated data, in detector operations, and in the trigger and data acquisition systems of the experiment.
Monte Carlo (MC) simulations are used to design and optimise the search strategy, for calculating the expected yields of signal events, and as part of the background estimation process.The simulated backgrounds include prompt diphoton production () as well as resonant Higgs production via  tℎ, vector boson fusion (VBF),  ℎ, gluon-gluon fusion (ggF), ℎ, ℎ and  bℎ.The main details regarding the generator settings for such processes are summarised in Table 1.In addition, for the prompt diphoton production the next-to-leading-order (NLO) matrix elements for up to one parton, and leading-order (LO) matrix elements for up to three partons were calculated with the Comix [27] and OpenLoops [28][29][30] libraries.They were matched to the Sherpa parton shower [31] using the MEPS@NLO prescription [32][33][34][35] with a dynamic merging cut [36] of 10 GeV.Photons were required to be isolated according to a smooth isolation-cone criterion [37].The ggF process is simulated with next-to-next-to-leading-order (NNLO) accuracy for arbitrary inclusive  → ℎ observables by reweighting the Higgs boson rapidity spectrum in Hj-MiNLO [38][39][40] to that of HNNLO [41].In the ℎ simulated sample, the overlap with the  tℎ production was removed using the diagram removal scheme [42,43].The decays of bottom and charm hadrons were simulated using the EvtGen 1.6.0program [44] for all Higgs boson samples.The normalisation of these samples accounts for the decay branching ratio calculated with HDECAY [45][46][47] and Prophecy4f [48][49][50].All the MC samples including Higgs bosons are simulated assuming  ℎ = 125 GeV.
In the following, the backgrounds arising from resonant Higgs processes will be categorised as  tℎ and ℎ (other), summarising events from VBF,  ℎ, ggF, ℎ, ℎ and  bℎ production, which account for only a minor contribution to the total background in the signal regions (SR).
The signal MC simulation includes the production of pairs of higgsinos through initial production of χ+ , where  ,  ′ are soft fermions that remain undetected.The mass of the gravitino was set to 1 MeV, and the decays χ0 1 → ℎ G and χ0 1 →  G were similarly assumed to be prompt.Events were generated from LO matrix elements with up to two extra partons, using the MadGraph5_aMC@NLO 2.9.3 [73] generator interfaced to Pythia 8.245.Events were matched using the CKKW-L prescription [75], with a matching scale set to one quarter of the mass of the pair-produced SUSY particles.Signal cross-sections were calculated to NLO in the strong coupling constant, adding the resummation of soft gluon emission at next-to-leading-logarithm accuracy (NLO+NLL) [76][77][78][79][80][81][82][83].The cross-section for ( χ0 1 ) = 150 GeV amounts to 3.83 ± 0.16 pb.
The generation of all simulated signal samples includes the effect of multiple interactions in the same and neighbouring bunch crossings (pile-up), modelled by overlaying the simulated hard-scattering event with inelastic   events generated with Pythia 8.186 using the NNPDF2.3loset of parton distribution functions [52] and the A3 set of tuned parameters [84].
The various background processes contributing to the signal regions are categorised into two classes: resonant and non-resonant backgrounds.The resonant backgrounds are those that contain Higgs boson decays into , and so are expected to have a peak in   near  ℎ .They are subdominant, and are determined from MC simulations.
For this analysis, photons are reconstructed [85] with pseudorapidity || < 2.37, while vetoing those photons in the EM calorimeter transition region 1.37 < || < 1.52 between the calorimeter barrel and endcap, and are required to have  T > 25 GeV.Signal photons are further required to satisfy 'Tight' selection and isolation requirements [85].Photons with a looser ('LoosePrime3' [85]) identification and without any isolation requirement, defined as baseline photons, are used to aid the data-driven estimation of some of the backgrounds, as described below.
Jets are reconstructed using the particle-flow algorithm, which combines both calorimeter and tracking information [86].The anti-  algorithm [87, 88] with a radius parameter of  = 0.4 is used, with a four-momentum recombination scheme.From the selected jets, -jets are identified using the DL1r tagger at the 77% efficiency working point, as determined from simulated  t events [89].
Electrons are reconstructed from isolated EM calorimeter energy deposits matched to ID tracks and are required to have || < 2.47, a transverse momentum  T > 4.5 GeV, and to satisfy the 'LooseAndBLayer' requirement defined in Ref.
[85], which is based on a likelihood using measurements of shower shapes in the calorimeter and track properties in the ID as input variables.Muons are reconstructed in the region || < 2.7 from MS tracks matching ID tracks.They are required to have  T > 4 GeV and satisfy the 'Medium' identification requirements [90], based on the number of hits in the different ID and MS subsystems, and on the significance of the charge to momentum ratio, /.
To avoid reconstruction ambiguities and double counting of analysis objects, electrons, muons and jets are removed from events if they are within Δ = 0.4 of any photon passing 'Loose' identification criteria [85].While this selection favours photons originating from electron bremsstrahlung over the parent electron itself, such photon candidates do not satisfy the photon isolation requirements and therefore do not enter the signal regions.Electrons and muons within Δ = 0.4 of a jet are also removed.
The missing transverse momentum,  miss T , is constructed [91] from the negative sum of all transverse momenta −  ( ì  T )  of all reconstructed objects  in the event.These consist of photons, jets, muons and electrons, as well as a track soft-term that includes the transverse momenta of all tracks associated with the primary vertex2 that are not reconstructed as one of the particles listed.Candidate events were selected using a combination of diphoton triggers [92, 93], and are required to have exactly two signal photons.To ensure good trigger and reconstruction efficiency the highest- T (leading) photon is required to have  T > 35 GeV, and the sub-leading photon to have  T > 25 GeV.Additionally, the diphoton invariant mass is required to lie in the range of 95 GeV <   < 160 GeV, and each photon is required to satisfy  T /  > 0.2.Events are required to have exactly two -tagged jets, and any event containing one or more reconstructed electrons or muons is rejected.
Three different signal regions are then defined (Table 2), to gain sensitivity to different χ0 1 mass hypotheses and to the different decay modes.The three regions are designed to be non-overlapping, allowing a statistical combination that further enhances sensitivity.The SR selections are optimised based on studies of simplified GMSB models with varying χ0 1 mass, and of differing values of the branching fraction Each of the three SRs requires a diphoton invariant mass   within 5 GeV of the observed mass of the Higgs boson of 125 GeV [94,95].The signal regions differ in their requirements on the invariant mass   of the -tagged jets, and in the requirement on the missing transverse momentum.SR1h targets those events that contain two Higgs bosons, and so has an   requirement that is also broadly consistent with  ℎ .A second region, SR1Z, is designed to select events where   is instead consistent with   , thus gaining sensitivity to events in which one χ0 1 decays into ℎ G and the other into  G.These two regions, SR1h and SR1Z, both target the lower end of the allowed χ0 1 mass range, a region in which  miss T is relatively small < 100 GeV.Both regions include an additional selection on the transverse momentum of the diphoton system   T , as well as on   T /  , to improve the signal-to-background ratio.SR2 targets higher-mass χ0 1 decays, which consequently demands high missing transverse momentum,  miss T > 100 GeV.The expected yield for these higher-mass models is small, so to enhance sensitivity SR2 has a rather loose selection on   , which can admit both Higgs and  boson decays.No selection on   T is included to retain more signal events.
The dominant backgrounds are from non-resonant events, those for which no peak at   ≈  ℎ is expected.They include prompt diphoton events and also those with fewer than two prompt photons that can contribute when other detected particles -most often those associated with jets -are falsely identified as photons.These prompt and fake photon non-resonant contributions to the signal regions are modelled using data-driven techniques, since there is a variety of different ways in which mismeasurements can contribute, not all of which are expected to be well-modelled in MC simulations.
The data-driven method used to determine the non-resonant backgrounds is known as the '2×2D sideband method', a technique developed in the context of ℎ →  measurements [96], and based on similar techniques used at the D0 [97] and CDF [98] experiments.The method makes extensive supporting measurements in data to improve the accuracy of the predictions from MC simulations.An outline of the method is given below, while further details may be found in Ref. [96].
Control regions are defined as outlined below, selecting events in the diphoton sidebands (  < 120 GeV and   > 130 GeV) using baseline photons.The looser requirements provide a large yield of prompt and fake photons, while remaining unbiased by the photon identification requirements of the triggers.These events are then partitioned into categories, based on the photon identification status (tight or not-tight) and the photon isolation status (isolated or non-isolated).The partitioning is performed separately for the leading and sub-leading photon candidates, thus providing 16 different observable categories of events.
The expected yields of , -jet ( ), jet- (  ) and jet-jet (  ) events in each region can be expressed separately using the yields obtained by selecting baseline photons, the prompt photon identification and isolation efficiencies, the fake-photon mis-identification and mis-isolation probabilities (fake-factors) and the relevant correlations between the fake-factors.It is assumed, based on previous studies [96], that any photon-photon efficiency correlations, photon-jet efficiency correlations or fake-photon mis-identification correlations are negligible.
A simultaneous likelihood fit constrained by the yields in each of the 16 observable categories is then performed to determine the parameters above, where the prompt photon identification and isolation efficiencies are determined using the diphoton MC sample.From the fitted parameters the fractions of ,  ,   and   events in the control regions are calculated and applied to the corresponding signal and validation regions (VR) defined below.The key assumption, that the fraction of the events of each of these four types is the same between the control region and the corresponding signal region, is justified by the smooth behaviour in   of the non-resonant backgrounds across and between the sidebands, and by the fact that the photon identification and isolation efficiencies do not depend strongly on   [99].
Since the background composition is expected to change in regions with different  miss T requirements, backgrounds are estimated separately for the low- miss T signal regions SR1h and SR1Z, and for the high- miss T signal region SR2.All control regions select events in the   sidebands (  < 120 GeV and   > 130 GeV).Additionally, the low- miss signal region, additional supporting measurements in the   sidebands in regions defined by inverting the  miss T and/or the -jet multiplicity requirements, are also used in conjunction with those closer to SR2, and the final background prediction is obtained by extrapolating from such regions with the aid of the ABCD method [100].
To test the accuracy of the background prediction, two VRs are defined.These regions have similar kinematic selections to the signal regions, except that they remove the requirement on   while they require that   be outside the selection band of the corresponding signal region, resulting in negligible signal contamination (≲ 3 %).The two VRs are designed to reflect the same  miss T spectra of the SRs, with VR1 requiring  miss T < 100 GeV and VR2 requiring  miss T > 100 GeV.The background estimation method described above is applied to the VRs in the same manner as for the SRs, where correction factors for the predicted numbers of ,  ,   and   events are determined in the respective CR and applied in both the VR and the SR.A schematic overview of the signal, control-and validation-region strategy is shown in Figure 2.
The overall yields and the shapes of the backgrounds in both VRs are found to be well predicted.Distributions of two of the most important kinematic variables and the corresponding agreement between data and the Standard Model (SM) prediction are illustrated in Figure 3.  Table 3: Observed and expected numbers of events in the three signal regions.The background category 'ℎ (other)' includes events originating from VBF,  ℎ, ggF, ℎ, ℎ and  bℎ, all subdominant in this signature.The table also includes model-independent 95% CL upper limits on the visible number of BSM events ( 95 obs ), the number of BSM events given the expected number of background events ( 95 exp ) and the visible BSM cross-section (⟨ ⟩ 95 obs ), all calculated from pseudo-experiments.The discovery -value ( 0 ) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.Systematic uncertainties are calculated for three classes of sources: from data-driven background methods, from the modelling of the detector response and from theoretical calculations of different processes that are used to normalise the corresponding MC predictions.The total size of the systematic uncertainty in the prediction for the background yields ranges from 15% to 39%, depending on the signal region.The dominant sources of systematic uncertainty are from the weights calculated by the 2×2D data-driven background method, and from the corresponding statistical uncertainty of the associated control regions.Systematic uncertainties associated with the detector response are found to be negligible.The theoretical uncertainties are found to be at the level of 3% to 6% depending on the SR.The statistical uncertainties are in the range of 40% to 77%, depending on the signal region, and so are the dominant source of uncertainty for the available data sample.

Channel
The   distributions associated with the three signal regions are shown in Figure 4.The observed yields in each SR, reported in Table 3, are found to be consistent with the corresponding SM predictions, and inconsistent with substantial additional contributions from the supersymmetric signal models under test.Table 3 also shows, for each SR, the upper limits at 95% confidence level (CL) on the possible number of beyond the Standard Model (BSM) events  95 , and on the visible BSM cross-section ⟨ ⟩ 95 obs =  95 obs / ∫ Ld, defined as the product of the ATLAS experimental acceptance and efficiency and the production crosssection.
A statistical combination of the three signal regions is performed, and expected and observed limits on the cross-section for higgsino pair-production are again calculated following the   s formalism [101] implemented in HistFitter [102], while using the profile likelihood ratio as a test statistic [103].The results of this statistical combination are shown in Figure 5 for the assumption B ( χ0 1 → ℎ G) = 100%.Under the same assumption, cross-sections above 1 pb are excluded at 95% CL for ( χ0 1 ) < 150 GeV, and the Events / 10 GeV    =13 TeV, 139 fb s

ATLAS
All limits at 95% CL -1 = 13 TeV, 139 fb s ) assuming it decays via either χ0 1 → ℎ G or χ0 1 →  G. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2.The ± 1 variation on the expected limit is shown.The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 .Values of B ( χ0 1 → ℎ G) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.theoretical prediction for the pure higgsino cross-section is excluded at 95% CL for neutralino masses < 320 GeV.

ATLAS
The branching fraction assumption can be relaxed in such a way that the higgsino is allowed to decay either via χ0 . Limits for this scenario are set as a function of higgsino mass and branching fraction, as shown in Figure 6.As expected, the sensitivity is strongest in the region with smaller higgsino mass, and larger B ( χ0 1 → ℎ G).Nevertheless, the exclusion at 95% CL extends down to B ( χ0 1 → ℎ G) = 36% for neutralino masses of 130 GeV.The sensitivity in the low-mass regime (( χ0 1 ) ≲ 200 GeV) is improved over the search in final states with multiple b-jets presented in Ref. [13].
In conclusion, a search is performed with the ATLAS detector using 139 fb −1 of proton-proton collisions at a centre-of-mass energy of √  = 13 TeV for excess events containing either pairs of Higgs bosons, or one Higgs boson with one  boson, in a final state containing a photon pair, a  b pair and missing transverse momentum.The observed data is in agreement with the SM prediction in all three tested signal regions and the results are used to derive upper limits on the pair production of higgsinos in gauge-mediated supersymmetry models.The observed 95% CL exclusion is extended up to ( χ0 1 ) = 320 GeV for B (

Figure 1 :
Figure 1: The diagram shows the topology considered by this search, the production of two neutralinos χ0 1 , which decay into a gravitino G and either a Higgs boson ℎ or a  boson.The search aims to select events in which a pair of photons is produced by a Higgs boson from one χ0 1 , while a pair of -quarks is produced by either the decay of a Higgs boson or a  boson from the other.

Figure 2 :
Figure 2: Schematic overview of the event categorisation.Events with (a) low- miss T and (b) high- miss Tare treated separately, while the invariant masses of the photons   and -jets   are used to categorise events.The resulting control regions (CR) are used to estimate backgrounds.This estimation is applied to both the validation regions (VR) and the signal regions (SR).

Figure 3 :
Figure 3: Distribution of the diphoton invariant mass in both validation regions: (a) VR1 and (c) VR2 as well as the missing transverse momentum in both validation regions: (b) VR1 and (d) VR2.The solid histograms are stacked to show the SM expectations after the 2×2D background estimation technique is applied.Background and signal predictions are normalised to the luminosity.The background category 'ℎ (other)' includes events originating from VBF,  ℎ, ggF, ℎ, ℎ and  bℎ, all subdominant in this signature.Statistical and systematic uncertainties are indicated by the shaded area.The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin.The first and last bins include the underflows and overflows respectively.

Figure 4 :
Figure 4: Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on   itself, for the three signal regions: (a) SR1h, (b) SR1Z and (c) SR2.The background estimation techniques described in the text are applied.The different backgrounds are stacked to add up to the total SM prediction in each bin.The predicted yields for signal benchmark models of varying χ0 1 mass are also overlaid (not stacked) assuming B ( χ0 1 → ℎ G) to equal (a) 100%, (b) 50% and (c) 100%.Background and signal predictions are normalised to the luminosity.The background category 'ℎ (other)' includes events originating from VBF,  ℎ, ggF, ℎ, ℎ and  bℎ, all subdominant in this signature.The sizes of the statistical and systematic uncertainties are indicated by the shaded areas.The lower panels show the ratio of the data to the SM prediction.Arrows indicate the borders of the signal region (|  − 125 GeV| < 5 GeV).The first and last bins include the underflows and overflows respectively.

Figure 5 :
Figure 5: Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B ( χ0 1 → ℎ G) = 100% for different χ0 1 masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2.The inner and outer bands indicate the 1 and 2 variation on the expected limit respectively.

Figure 6 :
Figure 6: Observed and expected 95% CL limits on the pure-higgsino branching fraction to B ( χ0 1 → ℎ G) as a function of the higgsino mass ( χ0 1) assuming it decays via either χ0 1 → ℎ G or χ0 1 →  G. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2.The ± 1 variation on the expected limit is shown.The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 .Values of B ( χ0 1 → ℎ G) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Table 1 :
Background MC simulated samples used in this analysis.The matrix element and parton shower generators, underlying-event tune, PDF set and cross-section order are provided.

Table 2 :
Definition of the three different signal region selections, where SR1h is aimed at final states containing two Higgs bosons, SR1Z is aimed at final states containing one Higgs boson and one  boson, and SR2 is aimed at higher-mass neutralinos decaying through either channel.