pterm = True
if args.snr is not None:
print 'Scaling distance to give SNR = {0}'.format(args.snr)
snr2 = 0
args.gwdist = 1
for ct, p in enumerate(pp):
inducedRes = (PALutils.createResiduals(psr[ct], np.pi/2-args.gwdec, args.gwra, \
args.gwchirpmass, args.gwdist, args.gwfreq, args.gwphase, \
args.gwpolarization, args.gwinc, psrTerm=pterm))
# compute snr
snr2 += PALutils.calculateMatchedFilterSNR(psr[ct], inducedRes, inducedRes)**2
# get total snr
snr = np.sqrt(snr2)
# scale distance appropiately
args.gwdist = snr/args.snr
print 'Scaled GW distance = {0} for SNR = {1}'.format(args.gwdist, args.snr)
# make residuals
for ct, p in enumerate(pp):
inducedRes = (PALutils.createResiduals(psr[ct], np.pi/2-args.gwdec, args.gwra, \
args.gwchirpmass, args.gwdist, args.gwfreq, args.gwphase, \
args.gwpolarization, args.gwinc, psrTerm=pterm))
