Using analytic results obtained in a meson effective theory, that includes all infrared sensitive logs, we build a dedicated Monte Carlo framework to describe QED corrections in $$ \overline{B} $$
B
¯
→$$ \overline{K} $$
K
¯
ℓ+ℓ− for a generic form factor. For the neutral mode $$ {\overline{B}}^0 $$
B
¯
0
→$$ {\overline{K}}^0 $$
K
¯
0
ℓ+ℓ−, we perform a detailed numerical comparison versus those obtained with the general-purpose photon-shower tool PHOTOS. The comparison indicates a good agreement, at the few per-mil level, when focusing on the rare mode only. In addition, our framework allows us to investigate the impact of the charmonium resonances. Interference effects, not described by PHOTOS in the experimental analysis, are found to be small in the dilepton invariant mass region q2< 6GeV2, which is used to determine $$ {R}_{K^{\left(\ast \right)}} $$
R
K
∗
. Using a semi-analytic framework we assess the full, rare and resonant, mode. Based thereupon, we discuss strategies to check the subtraction of the resonant mode, which has a sizeable impact at q2 ≈ 6GeV2 in the electron mode.