Skip to content
Snippets Groups Projects
Commit 926fe24b authored by Jonah Kanner's avatar Jonah Kanner :nerd:
Browse files

Try to make star always on plot

git-svn-id: https://svn.ligo.caltech.edu/svn/bayeswave/trunk@28 c56465c9-8126-4a4f-9d7d-ac845eff4865
parent 81fdc2e8
No related branches found
No related tags found
No related merge requests found
......@@ -191,9 +191,9 @@ maxodds=20
# -- Plotting helper function
def limit_odds(odds):
if odds > maxodds:
odds = maxodds
odds = 0.95*maxodds
elif odds < -maxodds:
odds = -maxodds
odds = -0.95*maxodds
return odds
infile = open('evidence.dat', 'r')
......@@ -210,10 +210,16 @@ print sig_noise
print sig_gl
plt.figure()
plt.plot(limit_odds(sig_gl), limit_odds(sig_noise), '*', markersize=15, mfc='yellow')
# plt.plot(limit_odds(sig_gl), limit_odds(sig_noise), '*', markersize=15, mfc='yellow')
plt.plot(sig_gl, sig_noise, '*', markersize=15, mfc='yellow')
plt.ylabel('log( E_signal / E_noise )')
plt.xlabel('log( E_signal / E_glitch )')
plt.axis([-maxodds, maxodds, -maxodds, maxodds])
maxodds = 1.1*np.array( [np.abs(sig_gl), np.abs(sig_noise), 20] ).max()
xmaxodds = 1.1*np.maximum(np.abs(sig_gl), 20)
ymaxodds = 1.1*np.maximum(np.abs(sig_noise), 20)
plt.axis([-xmaxodds, xmaxodds, -ymaxodds, ymaxodds])
plt.fill_between([0,maxodds], [0, 0], [maxodds, maxodds], facecolor='green', interpolate=True, alpha=0.3)
plt.fill_between([-maxodds,0,maxodds], [-maxodds,0,0], [-maxodds, -maxodds, -maxodds], facecolor='blue', interpolate=True, alpha=0.3)
......@@ -221,15 +227,15 @@ plt.fill_between([-maxodds,0], [-maxodds, 0], [maxodds, maxodds], facecolor='red
plt.grid()
plt.text(0.9*maxodds, 0.9*maxodds, 'Signal',
plt.text(0.9*xmaxodds, 0.9*ymaxodds, 'Signal',
horizontalalignment='right',
verticalalignment='top')
plt.text(-0.9*maxodds, 0.9*maxodds, 'Glitch',
plt.text(-0.9*xmaxodds, 0.9*ymaxodds, 'Glitch',
horizontalalignment='left',
verticalalignment='top')
plt.text(0.9*maxodds, -0.9*maxodds, 'Noise',
plt.text(0.9*xmaxodds, -0.9*ymaxodds, 'Noise',
horizontalalignment='right',
verticalalignment='bottom')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment