import numpy as np
import matplotlib.pyplot as plt
from math import pi, cos
from scipy import loadtxt, optimize
import os

M = 1.41

plt.figure(figsize=(10,7), dpi=80)
ax = plt.axes()

dat = loadtxt("./particles/particles.tsv", skiprows=0, delimiter="\t")
t = dat.transpose()[0]
tracers = dat.transpose()[1:]

r_max = tracers.max()

ax.axis([0, r_max * 1.1 / M, 0, t.max() * 1.0])


for i in range(0, len(tracers)):
    ax.plot(tracers[i] / M, t, label = r"m/M = {}".format((i+1)*0.2))

plot_name = r"Particle world lines"
ax.set_title(plot_name)
ax.set_xlabel("r/M")
ax.set_ylabel("t")
ax.legend()

    
plt.show()



import numpy as np
import matplotlib.pyplot as plt
from math import pi, cos
from scipy import loadtxt, optimize
import os

M = 1.41

plt.figure(figsize=(10,7), dpi=80)
ax = plt.axes()

dat = loadtxt("./particles/particles.tsv", skiprows=0, delimiter="\t")
t = dat.transpose()[0]
tracers = dat.transpose()[1:]

r_max = tracers.max()

ax.axis([0, r_max * 1.1 / M, 0, t.max() * 1.0])


for i in range(0, len(tracers)):
    ax.plot(tracers[i] / M, t, label = r"m/M = {}".format((i+1)*0.2))

plot_name = r"Particle world lines"
ax.set_title(plot_name)
ax.set_xlabel("r/M")
ax.set_ylabel("t")
ax.legend()

    
plt.show()


