#!/usr/bin/env python

import numpy as np
import matplotlib.pyplot as plt

#----------------------------------------------------------------------------------------#

data  = np.loadtxt("hexagomal_pristine_crystal_GOM.txt", skiprows=18)

#----------------------------------------------------------------------------------------#

fig = plt.figure()

plt.plot(data[:,0],data[:,1],'k-')
plt.yscale('log')

plt.ylim(10**(-2.0),10**7)
plt.grid(True)

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_1.png', bbox_inches='tight')
plt.close()

#----------------------------------------------------------------------------------------#

plt.plot(data[:,0],data[:,2],'k-')

plt.ylim(-0.5,0.5)
plt.grid(True,which="both")

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$-P_{12}/P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_2.png', bbox_inches='tight')

plt.close()

#----------------------------------------------------------------------------------------#

plt.plot(data[:,0],data[:,3],'k-')

plt.ylim(-1.0,1.0)
plt.grid(True,which="both")

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$P_{22}/P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_3.png', bbox_inches='tight')

plt.close()

#----------------------------------------------------------------------------------------#

plt.plot(data[:,0],data[:,4],'k-')

plt.ylim(-1.0,1.0)
plt.grid(True,which="both")

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$P_{33}/P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_4.png', bbox_inches='tight')

plt.close()

#----------------------------------------------------------------------------------------#

plt.plot(data[:,0],data[:,5],'k-')

plt.ylim(-0.5,0.5)
plt.grid(True,which="both")

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$-P_{43}/P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_5.png', bbox_inches='tight')

plt.close()

#----------------------------------------------------------------------------------------#

plt.plot(data[:,0],data[:,6],'k-')

plt.ylim(0.0,1.0)
plt.grid(True,which="both")

plt.xlabel(r"Scattering Angle $\Theta$ ($^\circ$)")
plt.title(r"$P_{44}/P_{11}$")

plt.savefig('geometric_optic_method_hexagonal_pristine_crystals_6.png', bbox_inches='tight')

plt.close()

#----------------------------------------------------------------------------------------#

About / Au sujet de

Ben
Created:
25 avril 2017 22:30:06

Updated:
25 avril 2017 22:30:06

License / Licence

MIT License

Abstract / Résumé