In this paper, we have considered the problem of recovering the time dependent source term for the fourth order time fractional pseudoparabolic equation theoretically and numerically, for the first time, by considering an additional measurement at the left boundary of the space interval. This is very challenging and interesting inverse problem with many important applications in various fields of physics and mechanics. The existence of unique solution to the problem has been discussed by means of the contraction principle and Banach Fixed-point theorem on a small time interval and the unique solvability theorem is proved. The stability results for the inverse problem have also been presented. However, since the governing equation is yet ill-posed (very slight errors in the additional input may cause relatively significant errors in the output force), the regularization of the solution is needed. Therefore, to get a stable solution, a regularized objective function is to be minimized for retrieval of the unknown coefficient of the forcing term. The proposed problem is discretized using the quintic B-spline (QnB-spline) collocation technique and has been reshaped as a non-linear least-squares optimization of the Tikhonov regularization function. The stability analysis of the direct numerical scheme has also been presented. The MATLAB subroutine lsqnonlin tool has been used to expedite the numerical computations. Both analytical and perturbed data are inverted and the numerical outcomes for two benchmark test examples are reported and discussed.