A modification of BPSK Costas loop, where low pass filters are applied separately from the loop, is considered. Such modification of the loop can work with non-sinusoidal signals or with some special phase detectors instead of multipliers. In this paper a general approach to the analytical computation of phase detector characteristic of the loop for general periodic signal waveforms is demonstrated and an effective simulation of the loop in signal's phase space is considered.