def test_ET(source_ground_based):
ET_filedirectory = load_directory + "/EinsteinTelescope/"
ET_filename = "ET_D_data.txt"
ET_filelocation = ET_filedirectory + ET_filename
T_obs = 4 * u.yr # Observing time in years
ET = detector.GroundBased("ET", T_obs, load_location=ET_filelocation, I_type="A")
[et_sample_x, et_sample_y, et_SNR] = snr.Get_SNR_Matrix(
source_ground_based, ET, var_x, sampleRate_x, var_y, sampleRate_y
)
def aLIGO():
aLIGO_filedirectory = load_directory + 'aLIGO/StrainFiles/'
aLIGO_filename = 'aLIGODesign.txt'
aLIGO_filelocation = aLIGO_filedirectory + aLIGO_filename
aLIGO = detector.GroundBased('aLIGO',
Ground_T_obs,
load_location=aLIGO_filelocation,
I_type='A')
return aLIGO
def ET():
ET_filedirectory = load_directory + 'EinsteinTelescope/StrainFiles/'
ET_filename = 'ET_B_data.txt'
ET_filelocation = ET_filedirectory + ET_filename
ET = detector.GroundBased('ET',
Ground_T_obs,
load_location=ET_filelocation,
I_type='A')
return ET
def ET():
ET_filedirectory = load_directory + "EinsteinTelescope/"
ET_filename = "ET_B_data.txt"
ET_filelocation = ET_filedirectory + ET_filename
ET = detector.GroundBased("ET",
Ground_T_obs,
load_location=ET_filelocation,
I_type="A")
return ET
def test_aLIGO(source_ground_based):
aLIGO_filedirectory = load_directory + "/aLIGO/"
aLIGO_filename = "aLIGODesign.txt"
aLIGO_filelocation = aLIGO_filedirectory + aLIGO_filename
T_obs = 4 * u.yr # Observing time in years
aLIGO = detector.GroundBased(
"aLIGO", T_obs, load_location=aLIGO_filelocation, I_type="A"
)
[aLIGO_sample_x, aLIGO_sample_y, aLIGO_SNR] = snr.Get_SNR_Matrix(
source_ground_based, aLIGO, var_x, sampleRate_x, var_y, sampleRate_y
)
def aLIGO_gwinc():
aLIGO_gwinc = detector.GroundBased(
"aLIGO gwinc",
T_obs,
noise_dict=noise_dict,
f_min=1.0,
f_max=1e4,
nfreqs=int(1e3),
)
aLIGO_gwinc.T_obs = [T_obs, T_obs_min, T_obs_max]
return aLIGO_gwinc
def aLIGO():
aLIGO_filedirectory = load_directory + "aLIGO/"
aLIGO_filename = "aLIGODesign.txt"
aLIGO_filelocation = aLIGO_filedirectory + aLIGO_filename
aLIGO = detector.GroundBased("aLIGO",
Ground_T_obs,
load_location=aLIGO_filelocation,
I_type="A")
aLIGO.S_n_f
return aLIGO
def aLIGO_gwinc():
aLIGO_gwinc = detector.GroundBased("aLIGO gwinc",
Ground_T_obs,
noise_dict=gwinc_dict)
return aLIGO_gwinc
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