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 [31 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 [597 KiB] Thumbnail [172 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 [285 KiB] Thumbnail [99 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 [123 KiB] Thumbnail [20 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 [137 KiB] Thumbnail [23 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 [64 KiB] Thumbnail [10 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 [78 KiB] Thumbnail [12 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 [278 KiB] Thumbnail [77 KiB] *C file 

Fig5a.pdf [21 KiB] HiDef png [183 KiB] Thumbnail [76 KiB] *C file 

Fig5b.pdf [15 KiB] HiDef png [152 KiB] Thumbnail [65 KiB] *C file 

Fig6.pdf [24 KiB] HiDef png [132 KiB] Thumbnail [56 KiB] *C file 

Fig7a.pdf [14 KiB] HiDef png [107 KiB] Thumbnail [49 KiB] *C file 

Fig7b.pdf [14 KiB] HiDef png [114 KiB] Thumbnail [52 KiB] *C file 

LHCbP[..].pdf [184 KiB] HiDef png [425 KiB] Thumbnail [109 KiB] *C file 
Created on 16 February 2019.Citation count from INSPIRE on 21 February 2019.