The direct $C P$ asymmetries of the decays $B^0 \rightarrow K^{*0} \mu^+ \mu^$ and $B^+ \rightarrow K^{+} \mu^+ \mu^$ are measured using $pp$ collision data corresponding to an integrated luminosity of 3.0$ fb^{1}$ collected with the LHCb detector. The respective control modes $B^0 \rightarrow J/\psi K^{*0}$ and $B^+ \rightarrow J/\psi K^{+}$ are used to account for detection and production asymmetries. The measurements are made in several intervals of $\mu^+ \mu^$ invariant mass squared, with the $\phi(1020)$ and charmonium resonance regions excluded. Under the hypothesis of zero $C P$ asymmetry in the control modes, the average values of the asymmetries are \begin{align} {\cal A}_{C P}(B^0 \rightarrow K^{*0} \mu^+ \mu^) &= 0.035 \pm 0.024 \pm 0.003, \cr {\cal A}_{C P}(B^+ \rightarrow K^{+} \mu^+ \mu^) &= \phantom{}0.012 \pm 0.017 \pm 0.001, \end{align} where the first uncertainties are statistical and the second are due to systematic effects. Both measurements are consistent with the Standard Model prediction of small $C P$ asymmetry in these decays.
Unbinned maximumlikelihood fits to the $ K ^+ \pi ^ \mu ^+\mu ^ $ mass distributions of the integrated data set for (a) $ { B ^0 \rightarrow K^{*0} \mu^+ \mu^}$ and (b) $\overline{ B }{} ^0 \rightarrow \overline{ K }{} ^{*0} \mu^+ \mu^$ decays for one magnet polarity, and (c) $ { B ^0 \rightarrow K^{*0} \mu^+ \mu^}$ and (d) $\overline{ B }{} ^0 \rightarrow \overline{ K }{} ^{*0} \mu^+ \mu^$ for the other. The blue, solid line represents the total fit, the red, shortdashed line represents the signal component and the grey, longdashed line represents the combinatorial background. 
Fig1.pdf [41 KiB] HiDef png [380 KiB] Thumbnail [283 KiB] *.C file 

Unbinned maximumlikelihood fits to the $ K ^+ \mu ^+\mu ^ $ mass distributions of the integrated data set for (a) $ { B ^+ \rightarrow K^{+} \mu^+ \mu^}$ and (b) $ B ^ \rightarrow K^{} \mu^+ \mu^$ decays for one magnet polarity, and (c) $ { B ^+ \rightarrow K^{+} \mu^+ \mu^}$ and (d) $ B ^ \rightarrow K^{} \mu^+ \mu^$ for the other. The blue, solid line represents the total fit, the red, shortdashed line represents the signal component and the grey, longdashed line represents the combinatorial background. 
Fig2.pdf [42 KiB] HiDef png [391 KiB] Thumbnail [290 KiB] *.C file 

Values of $ {\cal A}_{ C P }$ for $ { B ^0 \rightarrow K^{*0} \mu^+ \mu^}$ decays in each of the the 14 $ q^2$ bins used in the analysis. The error bars are the sum of the statistical and systematic uncertainties in quadrature. The dashed line represents the weighted average value, and the grey band indicates $\pm1\sigma$. The vertical red lines show the $\phi(1020)$, $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ , and $\psi {(2S)}$ regions, which are vetoed. 
Fig3.pdf [15 KiB] HiDef png [114 KiB] Thumbnail [113 KiB] *.C file 

Values of $ {\cal A}_{ C P }$ for $ { B ^+ \rightarrow K^{+} \mu^+ \mu^}$ decays in each of the the 17 $ q^2$ bins used in the analysis. The error bars are the sum of the statistical and systematic uncertainties in quadrature. The dashed line represents the weighted average value, and the grey band indicates $\pm1\sigma$. The vertical red lines show the $\phi(1020)$, $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ , and $\psi {(2S)}$ regions, which are vetoed. 
Fig4.pdf [16 KiB] HiDef png [129 KiB] Thumbnail [132 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2014032.gif Thumbnail 
Summary of the sources of systematic uncertainty for the measurements of $ {\cal A}_{ C P } ( {B \rightarrow K^{(*)} \mu^+ \mu^} )$. The ranges shown in parentheses indicate the minimum and maximum values of the systematic uncertainties in different $ q^2$ bins, while the numbers outside the parentheses are the values averaged over $ q^2$ . These may be outside the ranges as the uncertainties are determined by methods affected by statistical fluctuations. There is no systematic uncertainty due to duplicate candidates in the $ { B ^+ \rightarrow K^{+} \mu^+ \mu^}$ decay. 
Table_1.pdf [46 KiB] HiDef png [81 KiB] Thumbnail [37 KiB] tex code 

Values of $ {\cal A}_{ C P }$ in $ { B ^0 \rightarrow K^{*0} \mu^+ \mu^}$ decays in each of the 14 $ q^2$ bins used in the analysis. The first uncertainties are statistical and the second are systematic. 
Table_2.pdf [43 KiB] HiDef png [209 KiB] Thumbnail [110 KiB] tex code 

Values of $ {\cal A}_{ C P }$ in $ { B ^+ \rightarrow K^{+} \mu^+ \mu^}$ decays in each of the 17 $ q^2$ bins used in the analysis. The first uncertainties are statistical and the second are systematic. 
Table_3.pdf [43 KiB] HiDef png [260 KiB] Thumbnail [133 KiB] tex code 
Supplementary material full pdf 
Supple[..].pdf [77 KiB] 

This ZIP file contains supplemetary material for the publication LHCbPAPER2014032. The files are: Supplementary.pdf : An overview of the extra figures *.pdf, *.png, *.eps : The figures in various formats 
Fig1.pdf [15 KiB] HiDef png [120 KiB] Thumbnail [118 KiB] *C file 

Fig2.pdf [16 KiB] HiDef png [135 KiB] Thumbnail [134 KiB] *C file 
Created on 17 August 2019.Citation count from INSPIRE on 20 August 2019.