Terahertz (THz) communication signals are susceptible to severe degradation because of the molecular interaction with the atmosphere in the form of subsequent absorption and re-radiation. Recently, reconfigurable intelligent surface (RIS) has emerged as a potential technology to assist in THz communications by boosting signal power or providing virtual line-of-sight (LOS) paths. However, the re-radiated energy has either been modeled as a scattering component or as additive Gaussian noise in the literature. Since the precise characterization is still a work in progress, this paper presents the first comparative investigation of the performance of an RIS-aided THz system under these two extreme re-radiation models. In particular, we first develop a novel parametric channel model that encompasses both models of the re-radiation through a simple parameter change, and then utilize that to design a robust block-coordinate descent (BCD) algorithmic framework which maximizes a lower bound on channel capacity while accounting for imperfect channel state information (CSI). In this framework, the original problem is split into two sub-problems: a) receive beamformer optimization, and b) RIS phase-shift optimization. As the latter sub-problem (unlike the former) has no analytical solution, we propose three approaches for it: a) semi-definite relaxation (SDR) (high complexity), b) signal alignment (SA) (low complexity), and c) gradient descent (GD) (low complexity). The time complexities associated with the proposed approaches are explicitly derived. We analytically demonstrate the limited interference suppression capability of a passive RIS by deriving the stationary points of signal-to-interference and noise ratio (SINR) of a one-element RIS system with one interferer. Our numerical results also demonstrate that slightly better throughput is achieved when the re-radiation manifests as scattering.