The first observation of the $D^0 \to \pi^+\pi^\mu^+\mu^$ and $D^0 \to K^+K^\mu^+\mu^$ decays is reported using a sample of protonproton collisions collected by LHCb at a centerofmass energy of 8$ $TeV, and corresponding to 2$ $fb$^{1}$ of integrated luminosity. The corresponding branching fractions are measured using as normalization the decay $D^0 \to K^ \pi^+[\mu^+\mu^]_{\rho^0/\omega}$, where the two muons are consistent with coming from the decay of a $\rho^0$ or $\omega$ meson. The results are $\mathcal{B}(D^0 \to \pi^+\pi^\mu^+\mu^)=(9.64\pm0.48\pm0.51\pm0.97)\times10^{7}$ and $\mathcal{B}(D^0 \to K^+K^\mu^+\mu^)=( 1.54\pm0.27\pm0.09\pm0.16)\times10^{7}$, where the uncertainties are statistical, systematic, and due to the limited knowledge of the normalization branching fraction. The dependence of the branching fraction on the dimuon mass is also investigated.
Example diagrams describing the (left) short and (right) longdistance contributions to $ D ^0 \rightarrow h^+h^\mu^+\mu^$ decays, where $q=d,s$ and $h=\pi,K$. 
Fig1.pdf [84 KiB] HiDef png [112 KiB] Thumbnail [76 KiB] *.C file 

Distributions of $ m( D ^0 )$ for the $ D ^0 \rightarrow \pi ^+ \pi ^ \mu^+\mu^$ candidates in the low $ m(\mu ^+ \mu ^ )$ , $\eta$, $\rho ^0 /\omega $, $\phi$ and high $ m(\mu ^+ \mu ^ )$ regions. Fit projections are overlaid. 
Fig2.pdf [39 KiB] HiDef png [571 KiB] Thumbnail [443 KiB] *.C file 

Distributions of $ m( D ^0 )$ for the $ D ^0 \rightarrow K ^+ K ^ \mu^+\mu^$ candidates in the low $ m(\mu ^+ \mu ^ )$ , $\eta$ and $\rho ^0 /\omega $ regions. Fit projections are overlaid. No fit is performed in the $\eta$ region, where only two candidates are observed. 
Fig3.pdf [22 KiB] HiDef png [281 KiB] Thumbnail [246 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2017019.gif Thumbnail 
Yields of (top) $ D ^0 \rightarrow \pi ^+ \pi ^ \mu^+\mu^$ and (bottom) $ D ^0 \rightarrow K ^+ K ^ \mu^+\mu^$ signal decays, their significance with respect to the backgroundonly hypothesis, and ratio of efficiencies between signal and normalization decays ($R_\epsilon^i$) for each dimuonmass region. The yield and the significance ($\mathcal{S}$) are not reported for the $\eta$ region of $ D ^0 \rightarrow K ^+ K ^ \mu^+\mu^$ , where only two candidates are observed. 
Table_1.pdf [56 KiB] HiDef png [102 KiB] Thumbnail [51 KiB] tex code 

Branching fractions of (top) $ D ^0 \rightarrow \pi ^+ \pi ^ \mu^+\mu^$ and (bottom) $ D ^0 \rightarrow K ^+ K ^ \mu^+\mu^$ decays in different ranges of dimuon mass, where the uncertainties are statistical, systematic and due to the limited knowledge of the normalization branching fraction. The reported upper limits correspond to 90% (95%) CL. The correlations between the various dimuonmass ranges are reported in the supplemental material [23]. 
Table_2.pdf [48 KiB] HiDef png [112 KiB] Thumbnail [56 KiB] tex code 

Correlation coefficients between the $ D ^0 \rightarrow \pi ^+ \pi ^ \mu^+\mu^$ branching fractions in the dimuonmass ranges. 
Table_3.pdf [43 KiB] HiDef png [54 KiB] Thumbnail [25 KiB] tex code 

Correlation coefficients between the $ D ^0 \rightarrow K ^+ K ^ \mu^+\mu^$ branching fractions in the dimuonmass ranges. 
Table_4.pdf [42 KiB] HiDef png [64 KiB] Thumbnail [28 KiB] tex code 
This ZIP file contains supplemetary material for the publication LHCbPAPER2017019. The files are: SupplementaryLHCbPAPER2017019.pdf : An overview of the extra figures *.pdf, *.png, *.eps, *.C : The figures in variuous formats 
Fig4.pdf [24 KiB] HiDef png [273 KiB] Thumbnail [205 KiB] *C file 

Fig5a.pdf [21 KiB] HiDef png [146 KiB] Thumbnail [80 KiB] *C file 

Fig5b.pdf [15 KiB] HiDef png [117 KiB] Thumbnail [64 KiB] *C file 

Fig6.pdf [24 KiB] HiDef png [106 KiB] Thumbnail [59 KiB] *C file 

Fig7a.pdf [14 KiB] HiDef png [84 KiB] Thumbnail [53 KiB] *C file 

Fig7b.pdf [14 KiB] HiDef png [90 KiB] Thumbnail [56 KiB] *C file 

LHCbP[..].pdf [184 KiB] HiDef png [429 KiB] Thumbnail [283 KiB] *C file 
Created on 21 November 2020.