The LIGO/Virgo collaboration detected a compact object in the GW190814 event on August 14, 2019, with an estimated mass of \(\sim 2.59^{+0.08}_{-0.09}\) M \(_{\odot }\) , alongside a 23 M \(_{\odot }\) black hole. This observation ignited a scientific discussion due to the object’s unique position in the mass gap between neutron stars and black holes. In our study, we investigate the upper limit of the central magnetic field of anisotropic magnetized neutron stars, contextualized by the GW190814 event. We employ the density-dependent relativistic mean-field theory equation of state and assume a specific density dependence for the magnetic field. Our findings indicate that the magnetic field strength, anisotropy, and the magnetic field’s orientation significantly affect the physical properties of neutron stars. We examine two orientations: radial (where local magnetic fields point radially) and transverse (where they are perpendicular to the radial direction). Notably, stars with a transverse magnetic field orientation exhibit increased mass with higher anisotropy and magnetic field strength. Additionally, we show that the magnetic field, its orientation, and anisotropy substantially influence the tidal deformability of neutron stars. We aim to constrain the central magnetic field of neutron stars based on the tidal deformability measurement for a canonical \(1.4\) M \(_{\odot }\) neutron star, \(\Lambda _{1.4}=616^{+273}_{-158}\) , obtained from GW190814, which suggests a preference for a stiffer equation of state.