files/plot_saxs.py

   1 #! /usr/bin/env python
   2
   3 import matplotlib
   4 #matplotlib.use('GTK')
   5 matplotlib.use('Agg')
   6 import matplotlib.pyplot as plt
   7 import matplotlib.cm as cm
   8 import numpy as np
   9 import sys
  10
  11 x,y= np.loadtxt('saxs.data',usecols=(0,1),unpack=True)
  12 x1,y1= np.loadtxt('MODEL1_saxs.data',usecols=(0,2),unpack=True)
  13 x2,y2= np.loadtxt('MODEL2_saxs.data',usecols=(0,2),unpack=True)
  14 x3,y3= np.loadtxt('MODEL3_saxs.data',usecols=(0,2),unpack=True)
  15 x4,y4= np.loadtxt('MODEL4_saxs.data',usecols=(0,2),unpack=True)
  16 x5,y5= np.loadtxt('MODEL5_saxs.data',usecols=(0,2),unpack=True)
  17
  18
  19 plt.clf()
  20 plt.xlabel('r')
  21 plt.ylabel('P(r)')
  22 plt.xlim(x[0], x[-1])
  23 plt.plot(x,y,'-',color='black')
  24 plt.plot(x1,y1,'-',color='red')
  25 plt.plot(x2,y2,'-',color='lime')
  26 plt.plot(x3,y3,'-',color='blue')
  27 plt.plot(x4,y4,'-',color='magenta')
  28 plt.plot(x5,y5,'-',color='cyan')
  29 plt.legend(['reference','MODEL1','MODEL2','MODEL3','MODEL4','MODEL5'])
  30 plt.savefig('saxs.png')