A discovery of the rare decay $B^{+} \rightarrow \pi^{+} \mu^{+} \mu^{}$ is presented. This decay is observed for the first time, with 5.2 $\sigma$ significance. The observation is made using $pp$ collision data, corresponding to an integrated luminosity of 1.0 fb$^{1}$, collected with the LHCb detector. The measured branching fraction is (2.3 $\pm$ 0.6 (stat.) $\pm$ 0.1 (syst.))$\times 10^{8}$, and the ratio of the $B^{+} \rightarrow \pi^{+} \mu^{+} \mu^{}$ and $B^{+} \rightarrow \K^{+} \mu^{+} \mu^{}$ branching fractions is measured to be 0.053 $\pm$ 0.014 (stat.) $\pm$ 0.001 (syst.).
BDT output distribution for simulated $ B ^+ \rightarrow \pi ^+ \mu ^+ \mu ^ $ events (black solid line) and candidates taken from the mass sidebands in the data (red dotted line). Both distributions are normalised to unit area. The vertical line indicates the chosen cut value of 0.325. 
BDTout.pdf [14 KiB] HiDef png [131 KiB] Thumbnail [114 KiB] *.C file 

Invariant mass distribution for $ B ^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ $ candidates under the (a) $ K ^+ \mu ^+ \mu ^ $ and (b) $\pi ^+ \mu ^+ \mu ^ $ mass hypotheses with the fit projections overlaid. In the legend, "part. reco" refers to partially reconstructed background. The fit models are described in the text. 
pionfit.pdf [19 KiB] HiDef png [196 KiB] Thumbnail [165 KiB] *.C file 

kaonfit.pdf [18 KiB] HiDef png [184 KiB] Thumbnail [156 KiB] *.C file 

Invariant mass distribution of $ B ^+ \rightarrow \pi ^+ \mu ^+ \mu ^ $ candidates with the fit projection overlaid (a) in the full mass range and (b) in the region around the $ B $ mass. In the legend, "part. reco." and "combinatorial" refer to partially reconstructed and combinatorial backgrounds respectively. The discontinuity at 5500 $ {\mathrm{ Me V /}c^2}$ is due to the removal of data used for training the BDT. 
pimumufit.pdf [22 KiB] HiDef png [267 KiB] Thumbnail [219 KiB] *.C file 

pimumu[..].pdf [15 KiB] HiDef png [191 KiB] Thumbnail [144 KiB] *.C file 

Invariant mass distribution of $ B ^+ \rightarrow K ^+ \mu ^+ \mu ^ $ candidates with the fit projection overlaid (a) in the full mass range and (b) in the region around the $ B $ mass. In the legend, "combinatorial" refers to the combinatorial background. 
kmumufit.pdf [33 KiB] HiDef png [183 KiB] Thumbnail [145 KiB] *.C file 

kmumuf[..].pdf [18 KiB] HiDef png [203 KiB] Thumbnail [159 KiB] *.C file 

Invariant mass distribution of ${B^+\rightarrow J / \psi \pi ^+ }$ candidates with the fit projection overlaid. In the legend, "part. reco." and "combinatorial" refer to partially reconstructed and combinatorial backgrounds respectively. The fit model is described in the text. 
jpsipifit.pdf [32 KiB] HiDef png [227 KiB] Thumbnail [186 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2012020.gif Thumbnail 
Summary of systematic uncertainties. 
Table_1.pdf [70 KiB] HiDef png [45 KiB] Thumbnail [21 KiB] tex code 
Created on 11 July 2020.