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