Commit 51dba30f authored by hannah's avatar hannah

updating files paths

parent 1cae0af6
......@@ -6,24 +6,10 @@ from astropy.cosmology import Planck18 as cosmo
def getAGNPost(nFreqBins):
if nFreqBins==9:
path0 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq9/posterior.dat'.format(nFreqBins)
path1 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq9_1/posterior.dat'.format(nFreqBins)
path2 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq9_2/posterior.dat'.format(nFreqBins)
post0 = np.genfromtxt(path0,names=True)
post1 = np.genfromtxt(path1,names=True)
post2 = np.genfromtxt(path2,names=True)
elif nFreqBins==5:
path0 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq5/posterior.dat'.format(nFreqBins)
path1 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq5_1/posterior.dat'.format(nFreqBins)
path2 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq5_2/posterior.dat'.format(nFreqBins)
post0 = np.genfromtxt(path0,names=True)
post1 = np.genfromtxt(path1,names=True)
post2 = np.genfromtxt(path2,names=True)
elif nFreqBins==3:
path0 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq3/posterior.dat'.format(nFreqBins)
path1 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq3_1/posterior.dat'.format(nFreqBins)
path2 = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{}/nFreq3_2/posterior.dat'.format(nFreqBins)
path0 = '../results/agnosticPosterior/nFreq{}/posterior_0.dat'.format(nFreqBins)
path1 = '../results/agnosticPosterior/nFreq{}/posterior_1.dat'.format(nFreqBins)
path2 = '../results/agnosticPosterior/nFreq{}/posterior_2.dat'.format(nFreqBins)
post0 = np.genfromtxt(path0,names=True)
post1 = np.genfromtxt(path1,names=True)
post2 = np.genfromtxt(path2,names=True)
......@@ -50,7 +36,7 @@ def getAGNPost(nFreqBins):
'beta': betas,
'gamma': gammas}
saveLinesHere = '/home/ADF/middlehr/Documents/PTA/EPTAInterp/results/agn_model/f{0}/f{0}_plotting_things/'.format(nFreqBins)
saveLinesHere = 'lines/nFreq{}/'.format(nFreqBins)
return posteriorCombined, saveLinesHere
......
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