We take an additional step towards the optimization of the novel finite-range pseudopotential at constrained Hartree-Fock-Bogolyubov level and implement an optimization procedure within an axial code using harmonic oscillator basis. We perform the optimization using three different numbers of the harmonic oscillator shells. We apply the new parameterizations in the O-Kr part of the nuclear chart and isotopic chain of Sn, and we compare the results with experimental values and those given by a parameterization obtained using a spherical code.