def set_gaia(self):
beta_star1 = 5; r_DM = 1000
if gp.case == 1:
gamma_star1=0.1; r_star1=100; r_a1=100; gamma_DM=1; rho0=0.064
elif gp.case == 2:
gamma_star1=0.1; r_star1=250; r_a1=250; gamma_DM=0; rho0=0.400
elif gp.case == 3:
gamma_star1=0.1; r_star1=250; r_a1=np.inf; gamma_DM=1; rho0=0.064
elif gp.case == 4:
gamma_star1=0.1; r_star1=1000; r_a1=np.inf; gamma_DM=0; rho0=0.400
elif gp.case == 5:
gamma_star1=1.0; r_star1=100; r_a1=100; gamma_DM=1; rho0=0.064
elif gp.case == 6:
gamma_star1=1.0; r_star1=250; r_a1=250; gamma_DM=0; rho0=0.400
elif gp.case == 7:
gamma_star1=1.0; r_star1=250; r_a1=np.inf; gamma_DM=1; rho0=0.064
elif gp.case == 8:
gamma_star1=1.0; r_star1=1000; r_a1=np.inf; gamma_DM=0; rho0=0.400
self.params = [beta_star1, r_DM, gamma_star1, r_star1, r_a1, gamma_DM, rho0]
AAA = gh.myfill(100*gamma_star1) # 100
BBB = gh.myfill(10*beta_star1) # 050
CCC = gh.myfill(100*r_star1/r_DM) # 100
DDD = gh.myfill(100*r_a1/r_star1) # 100
EEEE = {0: "core",
1: "cusp"
}[gamma_DM] # core
FFFF = gh.myfill(1000*rho0, 4)
self.longdir = 'gs'+AAA+'_bs'+BBB+'_rcrs'+CCC+'_rarc'+DDD+'_'+EEEE+'_'+FFFF+'mpc3_df/'
self.dir = self.modedir + self.longdir
## new variable to hold the .dat input file
self.datafile = self.dir + 'dat'
self.massfile = self.dir+'enclosedmass/enclosedmass_0.txt'
self.nufiles.append(self.dir+'nu/nunotnorm_0.txt') # all comp.
self.sigfiles.append(self.dir+'siglos/siglos_0.txt')
self.kappafiles.append(self.dir+'kappalos/kappalos_0.txt')
self.nufiles.append(self.dir+'nu/nunotnorm_1.txt') # first and only
self.sigfiles.append(self.dir+'siglos/siglos_1.txt')
self.kappafiles.append(self.dir+'kappalos/kappalos_1.txt')
def set_walk(self):
self.dir = self.machine + 'DTwalk/'
if gp.case == 0:
gamma_star1 = 0.1; gamma_star2 = 1.0 # 1. or 0.1
beta_star1 = 5.0; beta_star2 = 5.0 # fixed to 5
r_star1 = 1000.; r_star2 = 1000. # 500 or 1000
r_a1 = 1.0; r_a2 = 1.0
gamma_DM = 0 # 0 or 1
elif gp.case == 1:
gamma_star1 = 1.0; gamma_star2 = 1.0 # 1. or 0.1
beta_star1 = 5.0; beta_star2 = 5.0 # fixed to 5
r_star1 = 500.; r_star2 = 1000. # 500 or 1000
r_a1 = 1.0; r_a2 = 1.0
gamma_DM = 0 # core
elif gp.case == 2:
gamma_star1 = 1.0; gamma_star2 = 1.0 # 1. or 0.1
beta_star1 = 5.0; beta_star2 = 5.0 # fixed to 5
r_star1 = 500.; r_star2 = 1000. # 500 or 1000
r_a1 = 1.0; r_a2 = 1.0
gamma_DM = 1 # cusp
elif gp.case == 4:
gamma_star1 = 0.1; gamma_star2 = 0.1 # 1. or 0.1
beta_star1 = 5.0; beta_star2 = 5.0 # fixed to 5
r_star1 = 100.; r_star2 = 500. # 500 or 1000
r_a1 = 1.0; r_a2 = 1.0
gamma_DM = 1 # cusp
elif gp.case == 5:
gamma_star1 = 0.1; gamma_star2 = 0.1 # 1. or 0.1
beta_star1 = 5.0; beta_star2 = 5.0 # fixed to 5
r_star1 = 100.; r_star2 = 500. # 500 or 1000
r_a1 = 2.5; r_a2 = 2.5
gamma_DM = 0 # core
alpha_DM = 1; beta_DM = 3;
r_DM = 1000 # fixed to 1000pc
AAA = gh.myfill(100*gamma_star1) # 100
BBB = gh.myfill(10*beta_star1) # 050
CCC = gh.myfill(r_star1/10) # 100
DDD = gh.myfill(100*r_a1) # 100
EEE = {0: "core",
1: "cusp"
}[gamma_DM] # core
FFF = gh.myfill(100*gamma_star2) # 010
GGG = gh.myfill(10*beta_star2) # 050
HHH = gh.myfill(r_star2/10) # 100
III = gh.myfill(100*r_a2) # 100
JJJ = EEE # core
NNN = gh.myfill(3) # 003 # realization (1..10)
self.longdir = "c1_"+AAA+"_"+BBB+"_"+CCC+"_"+DDD+"_"+EEE+"_c2_"+FFF+"_"+GGG+"_"+HHH+"_"+III+"_"+JJJ+"_"+NNN+"_6d/"
self.dir = self.dir + self.longdir
self.massfile = self.dir + 'enclosedmass/enclosedmass_0.txt'
self.analytic = self.dir + 'samplepars'
# print('analytic set to ', self.analytic)
LINE = np.loadtxt(self.analytic, skiprows=0, unpack=False)
rho0 = LINE[19] # read from the corresp. samplepars file
self.params = [beta_star1, r_DM, gamma_star1, r_star1, r_a1, gamma_DM, rho0]
for i in np.arange(gp.pops+1): # 0, 1, 2 for gp.pops=2
self.nufiles.append(self.dir+'nu/nunotnorm_'+str(i)+'.txt')
self.sigfiles.append(self.dir+'siglos/siglos_'+str(i)+'.txt')
self.kappafiles.append(self.dir+'kappalos/kappalos_'+str(i)+'.txt')
return
