def LISA_ESA():
LISA_ESA = detector.SpaceBased(
"LISA_ESA",
T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IFO,
f_IFO_break,
Background=Background,
T_type=T_type,
)
LISA_ESA.T_obs = [T_obs, T_obs_min, T_obs_max]
LISA_ESA.L = [L, L_min, L_max]
LISA_ESA.A_acc = [A_acc, A_acc_min, A_acc_max]
LISA_ESA.f_acc_break_low = [
f_acc_break_low,
f_acc_break_low_min,
f_acc_break_low_max,
]
LISA_ESA.f_acc_break_high = [
f_acc_break_high,
f_acc_break_high_min,
f_acc_break_high_max,
]
LISA_ESA.A_IFO = [A_IFO, A_IFO_min, A_IFO_max]
LISA_ESA.f_IFO_break = [f_IFO_break, f_IFO_break_min, f_IFO_break_max]
return LISA_ESA
def test_LISA_prop1(source_space_based):
T_obs = 4 * u.yr # Observing time in years
L = 2.5e9 * u.m # armlength in meters
A_acc = 3e-15 * u.m / u.s / u.s
f_acc_break_low = 0.4 * u.mHz.to("Hz") * u.Hz
f_acc_break_high = 8.0 * u.mHz.to("Hz") * u.Hz
f_IMS_break = 2.0 * u.mHz.to("Hz") * u.Hz
A_IMS = 10e-12 * u.m
Background = False
T_type = "N"
LISA_prop1 = detector.SpaceBased(
"LISA_ESA",
T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
T_type=T_type,
)
[lisa_sample_x, lisa_sample_y, lisa_SNR] = snr.Get_SNR_Matrix(
source_space_based, LISA_prop1, var_x, sampleRate_x, var_y, sampleRate_y
)
def test_LISA_ex3():
LISA_Other_filedirectory = load_directory + 'LISA_Other/StrainFiles/'
LISA_ex3_filename = 'scg_6981.dat'
LISA_ex3_filelocation = LISA_Other_filedirectory + LISA_ex3_filename
#`I_type` should be Amplitude Spectral Density
LISA_ex3 = detector.SpaceBased('LISA Example 3',
SpaceBased_T_obs,
load_location=LISA_ex3_filelocation,
I_type='A')
def test_LISA_ex2():
LISA_ex2_filedirectory = load_directory + 'LISA_Other/StrainFiles/'
LISA_ex2_filename = 'LISA_sensitivity.txt'
LISA_ex2_filelocation = LISA_ex2_filedirectory + LISA_ex2_filename
#`I_type` should be Effective Noise Spectral Density
LISA_ex2 = detector.SpaceBased('LISA Example 2',
SpaceBased_T_obs,
load_location=LISA_ex2_filelocation,
I_type='E')
def test_LISA_ex1():
LISA_Other_filedirectory = load_directory + 'LISA_Other/StrainFiles/'
LISA_ex1_filename = 'LISA_Allocation_S_h_tot.txt'
LISA_ex1_filelocation = LISA_Other_filedirectory + LISA_ex1_filename
#`I_type` should be Effective Noise Spectral Density
LISA_ex1 = detector.SpaceBased('LISA Example 1',
SpaceBased_T_obs,
load_location=LISA_ex1_filelocation,
I_type='E')
def test_LISA_prop2():
LISA_prop2 = detector.SpaceBased('LISA Approximate',
LISA_T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
T_type='A')
def LISA_prop1():
LISA_prop1 = detector.SpaceBased('LISA',
LISA_T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background)
return LISA_prop1
def test_LISA_ex3():
LISA_Other_filedirectory = load_directory + "LISA_Other/"
LISA_ex3_filename = "scg_6981.dat"
LISA_ex3_filelocation = LISA_Other_filedirectory + LISA_ex3_filename
# `I_type` should be Amplitude Spectral Density
LISA_ex3 = detector.SpaceBased(
"LISA Example 3",
SpaceBased_T_obs,
load_location=LISA_ex3_filelocation,
I_type="A",
)
LISA_ex3.fT
def test_LISA_ex2():
LISA_ex2_filedirectory = load_directory + "LISA_Other/"
LISA_ex2_filename = "LISA_sensitivity.txt"
LISA_ex2_filelocation = LISA_ex2_filedirectory + LISA_ex2_filename
# `I_type` should be Effective Noise Spectral Density
LISA_ex2 = detector.SpaceBased(
"LISA Example 2",
SpaceBased_T_obs,
load_location=LISA_ex2_filelocation,
I_type="E",
)
LISA_ex2.S_n_f
def test_LISA_ex1():
LISA_Other_filedirectory = load_directory + "LISA_Other/"
LISA_ex1_filename = "LISA_Allocation_S_h_tot.txt"
LISA_ex1_filelocation = LISA_Other_filedirectory + LISA_ex1_filename
# `I_type` should be Effective Noise Spectral Density
LISA_ex1 = detector.SpaceBased(
"LISA Example 1",
SpaceBased_T_obs,
load_location=LISA_ex1_filelocation,
I_type="E",
)
LISA_ex1.h_n_f
def LISA_prop2():
LISA_prop2 = detector.SpaceBased(
"LISA Approximate",
LISA_T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
Background_model=1,
T_type="A",
)
return LISA_prop2
def test_LISA_prop2():
LISA_prop2 = detector.SpaceBased(
"LISA Approximate",
LISA_T_obs_2,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
Background_model=1,
openingangle=60,
T_type="A",
)
LISA_prop2.S_n_f
def Get_Instrument(model):
if model == 0: #Einstein Telescope
load_name = 'ET_D_data.txt'
load_location = load_directory + 'EinsteinTelescope/StrainFiles/' + load_name
T_obs = 4 * u.yr #Observing time in years
T_obs_min = 1 * u.yr
T_obs_max = 10 * u.yr
instrument = detector.GroundBased('ET',
T_obs,
load_location=load_location,
I_type='A')
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
elif model == 1: #aLIGO
load_name = 'aLIGODesign.txt'
load_location = load_directory + 'aLIGO/StrainFiles/' + load_name
T_obs = 4 * u.yr #Observing time in years
T_obs_min = 1 * u.yr
T_obs_max = 10 * u.yr
instrument = detector.GroundBased('aLIGO',
T_obs,
load_location=load_location,
I_type='A')
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
elif model == 2: #NANOGrav 15 yr
###############################################
#NANOGrav calculation using 11.5yr parameters https://arxiv.org/abs/1801.01837
T_obs = 15 * u.yr #Observing time in years
T_obs_min = 10 * u.yr
T_obs_max = 30 * u.yr
sigma = 100 * u.ns.to('s') * u.s #rms timing residuals in seconds
N_p = 18 #Number of pulsars
cadence = 1 / (
2 * u.wk.to('yr') * u.yr
) #Avg observation cadence of 1 every 2 weeks in num/year
instrument = detector.PTA('NANOGrav', T_obs, N_p, sigma, cadence)
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
elif model == 3: #SKA (2030s)
###############################################
#SKA calculation using parameters and methods from arXiv:0804.4476 section 7.1
T_obs = 15 * u.yr #Observing time (years)
T_obs_min = 10 * u.yr
T_obs_max = 30 * u.yr
sigma = 10 * u.ns.to('s') * u.s #rms timing residuals in nanoseconds
N_p = 20 #Number of pulsars
cadence = 1 / (u.wk.to('yr') * u.yr
) #Avg observation cadence of 1 every week in num/year
instrument = detector.PTA('SKA', T_obs, N_p, sigma, cadence)
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
elif model == 4: #Robson,Cornish,and Liu 2018, LISA (https://arxiv.org/pdf/1803.01944.pdf)
T_obs = 4 * u.yr #Observing time in years
T_obs_min = 1 * u.yr
T_obs_max = 10 * u.yr
L = 2.5e9 * u.m #armlength in meters
L_min = 1.0e7 * u.m
L_max = 1.0e11 * u.m
A_acc = 3e-15 * u.m / u.s / u.s #M/s**2
A_IMS = 1.5e-11 * u.m
f_IMS_break = 2. * u.mHz.to('Hz') * u.Hz
f_acc_break_low = .4 * u.mHz.to('Hz') * u.Hz
f_acc_break_high = 8. * u.mHz.to('Hz') * u.Hz
Background = False
instrument = detector.SpaceBased('LISA_Alt',
T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
T_type='A')
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
instrument.L = [L, L_min, L_max]
else: #L3 proposal
#Default Params!
T_obs = 4 * u.yr #Observing time in years
T_obs_min = 1 * u.yr
T_obs_max = 10 * u.yr
L = 2.5e9 * u.m #armlength in meters
L_min = 1.0e7 * u.m
L_max = 1.0e11 * u.m
f_acc_break_low = .4 * u.mHz.to('Hz') * u.Hz
f_acc_break_high = 8. * u.mHz.to('Hz') * u.Hz
f_IMS_break = 2. * u.mHz.to('Hz') * u.Hz
A_acc = 3e-15 * u.m / u.s / u.s
A_IMS = 10e-12 * u.m
Background = False
instrument = detector.SpaceBased('LISA_ESA',
T_obs,
L,
A_acc,
f_acc_break_low,
f_acc_break_high,
A_IMS,
f_IMS_break,
Background=Background,
T_type='N')
instrument.T_obs = [T_obs, T_obs_min, T_obs_max]
instrument.L = [L, L_min, L_max]
return instrument
