diff --git a/config.txt b/config.txt index e54a589..e9b1d4b 100644 --- a/config.txt +++ b/config.txt @@ -1,7 +1,7 @@ # Simulation part num_walker=20000 # Motion model part -delta=126e3 +delta=161e3 eta=0.0 # Distribution part tau=1e-3 @@ -15,9 +15,9 @@ techo_start=0e-6 techo_stop=40e-6 techo_steps=5 # STE part -tevo_start=5e-6 -tevo_stop=60e-6 -tevo_steps=1 +tevo_start=2e-6 +tevo_stop=50e-6 +tevo_steps=49 tmix_start=1e-5 tmix_stop=1e1 tmix_steps=31 diff --git a/main.py b/main.py index af0b289..7bcc978 100644 --- a/main.py +++ b/main.py @@ -8,7 +8,8 @@ motion = 'IsotropicAngle' distribution = 'Delta' # parameter = {} parameter = { - "angle": np.linspace(3, 150, num=148), + "angle": np.linspace(1, 180, num=180), + "eta": 0.1 # "sigma": 0.1, # "tau": np.logspace(-1, 1, num=3) } @@ -69,7 +70,7 @@ ax.plot(angles, tau_cc_angles, 'o', label='CC (4.9mus)') ax.plot(angles, tau_2_angles, 'o', label='F2') ax.legend() -np.savetxt('angle.dat', np.c_[angles, tau_cc_angles, tau_ss_angles, tau_2_angles], header='#x\tcc\tss\tf2') +np.savetxt('angle_eta.dat', np.c_[angles, tau_cc_angles, tau_ss_angles, tau_2_angles], header='#x\tcc\tss\tf2') fig2, ax2 = plt.subplots() ax2.plot(angles, np.array(tau_cc_angles)/np.array(tau_2_angles), 'o')