We extracted the constraint on BR(B^+\rightarrow K^+ \varphi)\timesBR(\varphi \rightarrow \mu \bar{\mu }) as a function of the scalar massm_\varphiand lifetime\tau _\varphifrom Fig. 4 in [20] which is model-independent. In Fig.2(left) we use this constraint to set a limit ...