Commit d50c5bc9 authored by Karsten Wiesner's avatar Karsten Wiesner 💬
Browse files

added a simple test script using synthlisa

parent 78dcc18e
import matplotlib as mpl
mpl.use('Agg')
import matplotlib.pyplot as plt
plt.rcParams['agg.path.chunksize'] = 10000
import numpy as np
from math import *
import numpy.fft as FFT
import random
import scipy
from numpy import transpose
from decimal import Decimal
from synthlisa import *
c = 2.99792458e8
patches = 256
OneYear = 2**25
NYear = 1.0
stime = 4.0
samples = int(NYear*OneYear/stime)
lisa = EccentricInclined(2.5e9/2.99792458e8,0.0,0.,-1,0)
Sl = 0 # no laser noise since it would be canceled with first-generation TDI
Sop = ((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
#Spm = 0
# Sop = 0
Spm = 3.0e-15**2/(4.0*math.pi*math.pi*c**2) ## 2.5e-48
Sop = 2*((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
noiseTDI = TDInoise(lisa, stime,Spm, stime,Sop, stime,0.0)
[tn, Xmn] = np.transpose(getobs(samples,stime,[noiseTDI.t, noiseTDI.Xm]))
specXmn = spect(Xmn,stime,patches=256)
freqmn, spXmn = specXmn[1:,0], specXmn[1:,1]
fig = plt.figure()
ax = fig.add_subplot(111)
#ax.plot(range(10))
ax.loglog(freqmn, spXmn, ls='-', color='green')
fig.savefig('/vagrant/swetas_precious_results.png')
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment