The Gamow-Teller (GT) transition is inarguably one of the most important nuclear weak transitions of the spin-isosopin $\sigma\tau$ type. It has many applications in nuclear and astrophysics. These include, but are not limited to, r-process $\beta$-decays, stellar electron captures, neutrino cooling rates, neutrino absorption and inelastic scattering on nuclei. The quasiparticle random phase approximation (QRPA) is an efficient way to generate GT strength distribution. In order to better understand both theoretical systematics and uncertainties, we compare the GT strength distributions, centroid and width calculations for $^{40-60}$Ti isotopes, using the pn-QRPA, Pyatov method (PM) and the Schematic model (SM). The pn-QRPA and SM are further sub-divided into three categories in order to highlight the role of particle-particle (pp) force and deformation of the nucleus in the GT strength calculations. In PM, we study only the influence of the pp force in the calculation. We also compare with experimental results and other calculations where available. We found that the inclusion of pp force and deformation significantly improves the performance of SM and pn-QRPA models. Incorporation of pp force leads to pinning down the centroid value in the PM. The calculated GT strength functions using the pn-QRPA (C) and SM (C) models are in reasonable agreement with measured data.