CMS-HIG-13-010 ; CERN-PH-EP-2015-116 | ||
A search for pair production of new light bosons decaying into muons | ||
CMS Collaboration | ||
31 May 2015 | ||
Phys. Lett. B 752 (2016) 146 | ||
Abstract: A search for the pair production of new light bosons, each decaying into a pair of muons, is performed with the CMS experiment at the LHC, using a dataset corresponding to an integrated luminosity of 20.7 fb−1 collected in proton-proton collisions at center-of-mass energy of √s= 8 TeV. No excess is observed in the data relative to standard model background expectation and a model independent upper limit on the product of the cross section, branching fraction, and acceptance is derived. The results are compared with two benchmark models, the first one in the context of the next-to-minimal supersymmetric standard model, and the second one in scenarios containing a hidden sector, including those predicting a non-negligible light boson lifetime. | ||
Links: e-print arXiv:1506.00424 [hep-ex] (PDF) ; CDS record ; inSPIRE record ; CADI line (restricted) ; |
Figures & Tables | Summary | Additional Figures | CMS Publications |
---|
Figures | |
![]() png pdf |
Figure 1:
Distribution of the invariant masses m1μμ vs. m2μμ for the isolated dimuon events following the application of all constraints except the m1μμ≃m2μμ requirement of compatibility within the detector resolution. The compatible diagonal signal region (outlined with dashed lines) contains one data event (triangle) at m1μμ= 0.33 GeV and m2μμ= 0.22 GeV. There are also nine data events (white circles) which fail the m1μμ≃m2μμ compatibility requirement. The color scale indicates the expected SM background in range 2mμ<m1μμ,m2μμ<2mτ. |
![]() png pdf |
Figure 2-a:
(a) for benchmark model 1: 95%CL upper limits from this search for the NMSSM scenarios with ma1= 0.25 GeV (dashed curve), ma1= 2 GeV (dash-dotted curve) and ma1= 3.55 GeV (dotted curve) on σ(pp→h1/2→2a1)B2(a1→2μ) as a function of mh1 in the range 86 <mh1< 125 GeV and of mh2 for mh2> 125 GeV. As an illustration, the limits are compared to the predicted rate (solid curve) obtained using a simplified scenario with σ(pp→hi→2a1)= 0.008 σSM, which yields predictions for the rates of dimuon pair events comparable to the obtained experimental limits, and B(a1→2μ)= 7.7%. The chosen B(a1→2μ) is taken from [17, Dermisek and Gunion, PRD 81 (2010) 075003] for ma1=2 GeV and tanβ= 20. (b) for benchmark model 2: 95% CL upper limits (black solid curves) from this search on σ(pp→h→2γD+X)B(h→2γD+X) (with mn1= 10 GeV, mnD= 1 GeV) in the plane of two of the parameters (ε and mγD) for the dark SUSY scenarios, along with constraints from other experiments. The colored contours represent different values of B(h→2γD+X) in the range 0.1--40%. |
![]() png pdf |
Figure 2-b:
(a) for benchmark model 1: 95%CL upper limits from this search for the NMSSM scenarios with ma1= 0.25 GeV (dashed curve), ma1= 2 GeV (dash-dotted curve) and ma1= 3.55 GeV (dotted curve) on σ(pp→h1/2→2a1)B2(a1→2μ) as a function of mh1 in the range 86 <mh1< 125 GeV and of mh2 for mh2> 125 GeV. As an illustration, the limits are compared to the predicted rate (solid curve) obtained using a simplified scenario with σ(pp→hi→2a1)= 0.008 σSM, which yields predictions for the rates of dimuon pair events comparable to the obtained experimental limits, and B(a1→2μ)= 7.7%. The chosen B(a1→2μ) is taken from [17, Dermisek and Gunion, PRD 81 (2010) 075003] for ma1=2 GeV and tanβ= 20. (b) for benchmark model 2: 95% CL upper limits (black solid curves) from this search on σ(pp→h→2γD+X)B(h→2γD+X) (with mn1= 10 GeV, mnD= 1 GeV) in the plane of two of the parameters (ε and mγD) for the dark SUSY scenarios, along with constraints from other experiments. The colored contours represent different values of B(h→2γD+X) in the range 0.1--40%. |
Tables | |
![]() png pdf |
Table 1:
Event selection efficiencies ϵsim(mh1,ma1) and ϵsim(mγD,cτγD), as obtained from simulation, the geometric and kinematic acceptances αgen(mh1,ma1) and αgen(mγD,cτγD), calculated using only generator-level information, and their ratios (with statistical uncertainties), for a few representative NMSSM and dark SUSY benchmark samples. The experimental data-to- simulation scale factor (ϵdata/ϵsim, described later) is not applied. |
Summary |
A search for pairs of new light bosons produced in the decay of a Higgs boson, that subse- quently decay to pairs of oppositely charged muons, is presented. One event is observed in the signal region, with 2.2 ± 0.7 events expected from the SM backgrounds. A model independent upper limit at 95% CL on the product of the cross section, branching fraction, and acceptance is obtained. This limit is valid for light boson masses in the range 2mμ<ma<2mτ. The obtained results allow a straightforward interpretation within a broad range of physics models that predict the same type of signature. The results are compared with two benchmark models in the context of the NMSSM and dark SUSY, including scenarios predicting a non-negligible light boson lifetime. |
Additional Figures | |
![]() png pdf |
Additional Figure 1:
Separation between dimuon vertexes. The distance between vertexes measured for two dimuons (Δz=zμμ1−zμμ2) for a NMSSM simulated sample with mh1= 125 GeV and ma1= 2 GeV and for dimuon and ``orphan'' muon (Δz=zμμ−zμ) in events from bˉb-background-enriched part of data sample. |
![]() png pdf |
Additional Figure 2:
Dimuon isolation. The dimuon isolation distribution max(Isoμμ1,Isoμμ2) in a NMSSM simulated sample with mh1= 125 GeV and ma1= 2 GeV and in events from bˉb-background-enriched part of data sample. |
![]() png pdf |
Additional Figure 3:
bˉb background 1D template. The dimuon invariant mass distribution for a bˉb-background-enriched data sample (solid circles) in which at least one of the muons is "high pT" (pT> 17 GeV,|η|< 0.9), and the template (solid line) used to describe the low mass resonances, where the η,ω and ϕ are modeled using a gaussian distribution and the J/ψ is parametrized with a Crystal Ball function, this distribution is commonly referred as B17+8. |
![]() png pdf |
Additional Figure 4:
bˉb background 1D template. The dimuon invariant mass distribution for a bˉb-background-enriched data sample (solid circles) in which both muons satisfy the minimum pt and eta requirements (pT> 8 GeV,|η|<2.4), and the template (solid line) used to describe the low mass resonances, where the η,ω and ϕ are modeled using a gaussian distribution and the J/ψ is parametrized with a Crystal Ball function, this distribution is commonly referred as B8+8. |
![]() png pdf |
Additional Figure 5:
bˉb background 3D template. The analytical template for the distribution of the dimuon masses mμμ1 and mμμ2 obtained using a bˉb-background-enriched data sample. |
![]() png pdf |
Additional Figure 6:
Prompt double J/ψ control sample (Invariant mass selection). Mass distribution for the first (mμμ1) and second (mμμ2) dimuon in data after applying initial double J/ψ control sample selection (i.e trigger, muon pt, muon eta), the long tails correspond to J/ψ's coming from b-decays or non-resonant dimuon pairs from semileptonic decays of b-quark, only events in which the dimuon invariant mass is between 2.8 and 3.3 GeV are selected (events inside the box). |
![]() png pdf |
Additional Figure 7:
Prompt double J/ψ control sample ("non-prompt" subtraction 1). Isolation distribution for the first (Isoμμ1) and second (Isoμμ2) dimuon in data (including the invariant mass cut 2.8 <mμμ< 3.3 GeV), most of the "prompt" events are located inside region A, while the contribution from "non-prompt" events in this region is subtracted using non-isolated sidebands (B,C,D) in data. |
![]() png pdf |
Additional Figure 8:
Prompt double J/ψ control sample ("non-prompt" subtraction 2). The Lxy distribution for the second dimuon (μμ2) after applying the control sample selection including the invariant mass cut 2.8 <mμμ< 3.3 GeV and dimuon isolation Isoμμ< 2 GeV (black dots), the combined ("prompt"+"non-prompt") model is obtained by fitting the two templates to data, the "prompt" template was obtained from Monte Carlo signal sample and the "non-prompt" from a non-isolated sideband in data. |
![]() png pdf |
Additional Figure 9:
Prompt double J/ψ control sample. The two dimuon invariant mass distribution (Mμμ1+μμ2) in data following the control sample selection with "non-prompt" contribution subtracted (black dots), and the corresponding DPS and SPS templates from Monte Carlo simulation (normalized to the number of events in data). |
![]() png pdf |
Additional Figure 10:
Model Independent Limit. Model independent 95% CL upper limit on the product of the cross section times branching fraction times acceptance: σ(pp→2a+X)×B2(a→2μ)×αgen , where αgen is the generator level acceptance including the fiducial cut (Lxy< 4.4 cm and |Lz|< 34.5 cm). |
![]() png pdf |
Additional Figure 11:
Signal candidate event display. Display of the one event in the data which survived all selection requirements. Muons are shown as red curves with parameters labeled. Muon chambers are shown as red shapes. Inner silicon tracker is outlined with charged tracks shown as green curves, the image shows a 2D cross section of the CMS detector. |
![]() png pdf |
Additional Figure 12:
Signal candidate event display. Display of the one event in the data which survived all selection requirements. Muons are shown as red curves with parameters labeled. Muon chambers are shown as red shapes. Inner silicon tracker is outlined with charged tracks shown as green curves, the image shows a 3D representation of the CMS detector. |
![]() png pdf |
Additional Figure 13:
Interpolated acceptance ϵrec of the analysis selection requirements, evaluated using the dark SUSY simulated MC samples for several sets of mγD and cτγD (solid bullets) and extrapolated to other points in the (mγD, cτγD) plane. The plot shows the interpolation of the acceptance as a function of cτγD (solid line). |
![]() |
Compact Muon Solenoid LHC, CERN |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |