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 [185 KiB] Thumbnail [35 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 [937 KiB] Thumbnail [180 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 [502 KiB] Thumbnail [106 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 [100 KiB] Thumbnail [17 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 [109 KiB] Thumbnail [19 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 [52 KiB] Thumbnail [9 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 [60 KiB] Thumbnail [11 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 [459 KiB] Thumbnail [81 KiB] *C file 

Fig5a.pdf [21 KiB] HiDef png [171 KiB] Thumbnail [37 KiB] *C file 

Fig5b.pdf [15 KiB] HiDef png [130 KiB] Thumbnail [32 KiB] *C file 

Fig6.pdf [24 KiB] HiDef png [124 KiB] Thumbnail [27 KiB] *C file 

Fig7a.pdf [14 KiB] HiDef png [101 KiB] Thumbnail [24 KiB] *C file 

Fig7b.pdf [14 KiB] HiDef png [108 KiB] Thumbnail [25 KiB] *C file 

LHCbP[..].pdf [184 KiB] HiDef png [657 KiB] Thumbnail [114 KiB] *C file 
Created on 11 December 2018.Citation count from INSPIRE on 11 December 2018.