Skip to content
Snippets Groups Projects
Commit 709cbaee authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Updated the plotting scripts.

parent 7b0d18f5
No related branches found
No related tags found
No related merge requests found
...@@ -34,6 +34,8 @@ rc('font', family='serif') ...@@ -34,6 +34,8 @@ rc('font', family='serif')
import sys import sys
import os import os
print "Plotting..."
# Read Quickshed accelerations # Read Quickshed accelerations
data=loadtxt("particle_dump.dat") data=loadtxt("particle_dump.dat")
id = data[:,0] id = data[:,0]
...@@ -142,10 +144,10 @@ grid() ...@@ -142,10 +144,10 @@ grid()
subplot(313, title="Acceleration along Z") subplot(313, title="Acceleration along Z")
#plot(id, errz_g , 'gs') #plot(id, errz_g , 'gs')
plot(id, errz_bh , 'rx') plot(id, errz_new , 'b.', label="QuickShed")
plot(id, errz_new , 'b.') plot(id, errz_bh , 'rx', label="Legacy")
text(id[-1], 0.18, "B-H: $%5.3f\\pm%5.3f$\n QuickShed: $%5.3f\\pm%5.3f$"%(meanz_bh, stdz_bh, meanz_new, stdz_new), backgroundcolor="w", va="top", ha="right" ) #text(id[-1], 0.18, "B-H: $%5.3f\\pm%5.3f$\n QuickShed: $%5.3f\\pm%5.3f$"%(meanz_bh, stdz_bh, meanz_new, stdz_new), backgroundcolor="w", va="top", ha="right" )
legend(loc="upper right")
ylim(-0.2, 0.2) ylim(-0.2, 0.2)
xlim(0,id[-1]) xlim(0,id[-1])
...@@ -177,6 +179,6 @@ subplot(313, title="Acceleration along Z") ...@@ -177,6 +179,6 @@ subplot(313, title="Acceleration along Z")
#hist(errz_g, bins=bins, normed=1, histtype='step', rwidth=0.01, color='g') #hist(errz_g, bins=bins, normed=1, histtype='step', rwidth=0.01, color='g')
hist(errz_bh, bins=bins, normed=1, histtype='step', rwidth=0.01, color='r') hist(errz_bh, bins=bins, normed=1, histtype='step', rwidth=0.01, color='r')
hist(errz_new, bins=bins, normed=1, histtype='step', rwidth=0.01, color='b') hist(errz_new, bins=bins, normed=1, histtype='step', rwidth=0.01, color='b')
xlim(-0.03, 0.03) #xlim(-0.03, 0.03)
savefig("histogram.png") savefig("histogram.png")
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment