  We investigate how well a simple leading order perturbation theory model of
the bispectrum can fit the BAO feature in the measured bispectrum monopole of
galaxies. Previous works showed that perturbative models of galaxy bispectrum
start failing at the wavenumbers of k ~ 0.1 Mpc/h. We show that when the BAO
feature in the bispectrum is separated it can be successfully modeled up to
much higher wavenumbers. We validate our modeling on GLAM simulations that were
run with and without the BAO feature in the initial conditions. We also
quantify the amount of systematic error due to BAO template being offset from
the true cosmology. We find that the systematic errors do not exceed 0.3 per
cent for reasonable deviations from the true cosmology.
