This paper deals with numerical treatment of singularly perturbed parabolic differential equations having delay on the zeroth and first order derivative terms. The solution of the considered problem exhibits boundary layer behaviour as the perturbation parameter tends to zero. The equation is solved using θ-method in temporal discretization and exponentially fitted finite difference method in spatial discretization. The stability of the scheme is proved by using solution bound technique by constructing barrier functions. The parameter uniform convergence analysis of the scheme is carried out and it is shown to be accurate of order O(N −2 N −1 +cε + (∆t) 2) for the case θ = 1 2 , where N is the number of mesh points in spatial discretization and ∆t is the mesh size in temporal discretization. Numerical examples are considered for validating the theoretical analysis of the scheme