In [1]:
from odin import xray
In [3]:
# load up the shotset from disk
ss = xray.Shotset.load('shotset.shot')
print ss
In [4]:
ip = ss.intensity_profile()
figure()
plot(ip[:,0], ip[:,1], lw=2)
show()
In [8]:
q_values = array([1.9]) # not gold but you get the idea...
# ---- THIS IS THE IMPORTANT ONE ----
rings = ss.to_rings(q_values, num_phi=360)
# -----------------------------------
corr = rings.correlate_intra(1.9, 1.9, mean_only=True)
figure()
plot(rings.phi_values, corr)
xlabel(r'$\phi$')
ylabel('Correlation')
show()
In [ ]: