md_trajectory = nlread('alumina_trajectory.nc')[-1]

partial_sq = PartialStructureFactor(md_trajectory,
                                    pair_selection=[Aluminum, Aluminum])

# Get the q-values and the partial structure factor.
q_values = partial_sq.qRange().inUnitsOf(Angstrom**-1)
s_q      = partial_sq.data()

# Plot the data using pylab.
import pylab

pylab.plot(q_values, s_q, label='Al-Al structure factor')
pylab.xlabel('q (1/Ang)')
pylab.ylabel('S(q)')
pylab.legend()

pylab.show()