but the shapes and relative weights of the cross sections are well evaluated with the previous formalism. The wave function for the $f_1(1285)$ is given by $$\label{Eq:24}
|f_1(1285)\rangle=-\frac{1}{2}(K^{*+}K^-+K^{*0}\bar{K}^{0}-K^{*-}K^+-\bar{K}^{*0}K^0)$$ and then everything is identical to what have been done so far, except that one has four diagrams. The evaluation of the widths for $f_1(1285)\to \pi a_0(980)$ and $f_1(1285)\to \pi K \bar K$ at the peak of the $f_1(1285)$ are done in Refs. [@Aceti:2015zva] and [@Aceti:2015pma], respectively. All that must be done is to perform the same evaluation as a function of $M_{\rm inv}$ (for the $f_1(1285)$ at higher energies), and implement the $f_1(1285)$ propagator in Eq. (\[Eq:5\]). The results are shown in Figs. \[Fig:8\], \[Fig:9\] and \[Fig:9bis\], which should be compared with Figs. \[Fig:6\], \[Fig:4\] and \[Fig:4bis\], respectively. As we can see, not only the shapes but the absolute numbers are about the same as with the simplified wave function. Moreover, in Fig. \[Fig:9\] a peak develops for the $f_1(1285)\to \pi K\bar K$ decay around 1400 MeV (see Fig. \[Fig:9bis\] for more details). The peak and width of the distribution around this energy are in fair agreement with the mass of 1420 MeV and width of about 55 MeV quoted in 