An angular analysis of $B^0\to J/\psi K^+\pi^$ decays is performed, using protonproton collision data corresponding to an integrated luminosity of 3 ${\rm fb}^{1}$ collected with the LHCb detector. The $m(K^+\pi^)$ spectrum is divided into fine bins. In each $m(K^+\pi^)$ bin, the hypothesis that the threedimensional angular distribution can be described by structures induced only by $K^\ast$ resonances is examined, making minimal assumptions about the $K^+\pi^$ system. The data reject the $K^\ast$only hypothesis with a large significance, implying the observation of exotic contributions in a modelindependent fashion. Inspection of the $m(J/\psi\pi^)$ versus $m(K^+\pi^)$ plane suggests structures near $m(J/\psi\pi^)=4200$ MeV and 4600 MeV.
Spectrum of $K^\ast_J$ resonances from Ref. [7], with the vertical span of the boxes indicating $\pm\Gamma_0$, where $\Gamma_0$ is the width of each resonance. The horizontal dashed lines mark the $ m( K ^+ \pi ^ )$ physical region for ${B^0\rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K^+\pi^}$ decays, while the dotdashed lines mark the specific region, ${ m( K ^+ \pi ^ ) \in [1085,1445]}$ $\text{\,MeV}$ , employed for determining the significance of exotic contributions. 
Fig1.pdf [581 KiB] HiDef png [497 KiB] Thumbnail [91 KiB] *.C file 

Comparison of $m(J\psi \pi^)$ in the $ m( K ^+ \pi ^ ) \in[1085,1265]$ $\text{\,MeV}$ region between the backgroundsubtracted data and simulated events weighted by moments models with $J^k_{\rm max}=2$ and $J^k_{\rm max}=15$. 
Fig2.pdf [505 KiB] HiDef png [348 KiB] Thumbnail [81 KiB] *.C file 

Likelihoodratio test for exotic significance. The data shows a $10\sigma$ deviation from the pseudoexperiments generated according to the null hypothesis ($K^\ast_J$only contributions). 
Fig3.pdf [538 KiB] HiDef png [393 KiB] Thumbnail [91 KiB] *.C file 

Backgroundsubtracted 2D distribution of $m( { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \pi^)$ versus $ m( K ^+ \pi ^ ) $ in the region ${ m( K ^+ \pi ^ ) \in[745,1545]}$ $\text{\,MeV}$ . The intensity ($z$axis) scale has been highly truncated to limit the strong $K^\ast(892)^0$ contribution. 
Fig4.pdf [2 MiB] HiDef png [1 MiB] Thumbnail [190 KiB] *.C file 

Angle conventions as described in Ref. [18] for (a) $\overline{ B } {}^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} (\rightarrow \mu ^ \mu ^+ ) K ^ \pi ^+ $ and (b) ${ B ^0 \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} (\rightarrow \mu ^+ \mu ^ ) K ^+ \pi ^ }$. The leptonic ($\mu ^+ \mu ^ $) and hadronic ($K^+\pi^$) frames are backtoback with a common $\hat{y}$ axis. 
Fig5a.pdf [14 KiB] HiDef png [280 KiB] Thumbnail [59 KiB] *.C file 

Fig5b.pdf [14 KiB] HiDef png [280 KiB] Thumbnail [59 KiB] *.C file 

Fit to the invariant mass $m( { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K^+\pi^)$ in the $ m( K ^+ \pi ^ ) \in [1235,1265]$ $\text{\,MeV}$ bin. 
Fig6.pdf [558 KiB] HiDef png [438 KiB] Thumbnail [90 KiB] *.C file 

Comparison of $m(J/\psi\pi^)$ in the $m(K^+\pi^)\in[880,910]$ MeV region between the backgroundsubtracted data and simulation data weighted by moments models with $J^k_{\rm max}=2$ and $J^k_{\rm max}=15$. 
Fig7.pdf [479 KiB] HiDef png [304 KiB] Thumbnail [70 KiB] *.C file 

Comparison of $m(J/\psi\pi^)$ in the $m(K^+\pi^)\in[1265,1445]$ MeV region between the backgroundsubtracted data and simulation data weighted by moments models with $J^k_{\rm max}=2$ and $J^k_{\rm max}=15$. 
Fig8.pdf [500 KiB] HiDef png [333 KiB] Thumbnail [77 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2018043.gif Thumbnail 
The transversitybasis moments of the 28 orthonormal angular functions $f_i(\Omega)$ in Eq. 2 till spin2 in the $ K ^+ \pi ^ $ system. 
Table_1.pdf [66 KiB] HiDef png [260 KiB] Thumbnail [44 KiB] tex code 

The 12 angular terms for each additional spin$J$ wave in the $ K ^+ \pi ^ $ system, for $J\geq3$. 
Table_2.pdf [49 KiB] HiDef png [380 KiB] Thumbnail [55 KiB] tex code 
Created on 16 February 2019.Citation count from INSPIRE on 21 February 2019.