#!/usr/bin/python

import sys
import numpy as np
import matplotlib.pyplot as plt

# python ./compare_ITS_Pelliccioni.py ./dose_vs_E_kapton_150um.txt ./Pelliccioni2000-Table_A2.2_ElectronDoseEquivPerFluence.txt

ITS_dose_vs_E_filename             = sys.argv[1]
Pelliccioni2000_dose_vs_E_filename = sys.argv[2]

E_ITS, dose_ITS = np.loadtxt(ITS_dose_vs_E_filename, unpack=True) 
E_P2000, dose_P2000 = np.loadtxt(Pelliccioni2000_dose_vs_E_filename, unpack=True)


plt.loglog(E_ITS*1E6, dose_ITS, marker='.', label='ITS', markersize=15)
plt.loglog(E_P2000, dose_P2000, marker='^', label='Pelliccioni 2000 Table A2.2')
plt.xlabel('Energy (eV)')
plt.ylabel(r'Dose per Fluence (Sv / (#/cm$^2$))')
plt.grid(b=True, which='both') # https://stackoverflow.com/questions/9127434/how-to-create-major-and-minor-gridlines-with-different-linestyles-in-python
plt.legend()
plt.savefig("comparison_dose_vs_energy.png", dpi=1000)
plt.show()
3464cc7a880e342bc1b9c9f650a7729c1f5e6d22
<mask token>

<mask token>
plt.loglog(E_ITS * 1000000.0, dose_ITS, marker='.', label='ITS', markersize=15)
plt.loglog(E_P2000, dose_P2000, marker='^', label='Pelliccioni 2000 Table A2.2'
    )
plt.xlabel('Energy (eV)')
plt.ylabel('Dose per Fluence (Sv / (#/cm$^2$))')
plt.grid(b=True, which='both')
plt.legend()
plt.savefig('comparison_dose_vs_energy.png', dpi=1000)
plt.show()

<mask token>
ITS_dose_vs_E_filename = sys.argv[1]
Pelliccioni2000_dose_vs_E_filename = sys.argv[2]
E_ITS, dose_ITS = np.loadtxt(ITS_dose_vs_E_filename, unpack=True)
E_P2000, dose_P2000 = np.loadtxt(Pelliccioni2000_dose_vs_E_filename, unpack
    =True)
plt.loglog(E_ITS * 1000000.0, dose_ITS, marker='.', label='ITS', markersize=15)
plt.loglog(E_P2000, dose_P2000, marker='^', label='Pelliccioni 2000 Table A2.2'
    )
plt.xlabel('Energy (eV)')
plt.ylabel('Dose per Fluence (Sv / (#/cm$^2$))')
plt.grid(b=True, which='both')
plt.legend()
plt.savefig('comparison_dose_vs_energy.png', dpi=1000)
plt.show()

import sys
import numpy as np
import matplotlib.pyplot as plt
ITS_dose_vs_E_filename = sys.argv[1]
Pelliccioni2000_dose_vs_E_filename = sys.argv[2]
E_ITS, dose_ITS = np.loadtxt(ITS_dose_vs_E_filename, unpack=True)
E_P2000, dose_P2000 = np.loadtxt(Pelliccioni2000_dose_vs_E_filename, unpack
    =True)
plt.loglog(E_ITS * 1000000.0, dose_ITS, marker='.', label='ITS', markersize=15)
plt.loglog(E_P2000, dose_P2000, marker='^', label='Pelliccioni 2000 Table A2.2'
    )
plt.xlabel('Energy (eV)')
plt.ylabel('Dose per Fluence (Sv / (#/cm$^2$))')
plt.grid(b=True, which='both')
plt.legend()
plt.savefig('comparison_dose_vs_energy.png', dpi=1000)
plt.show()

#!/usr/bin/python

import sys
import numpy as np
import matplotlib.pyplot as plt

# python ./compare_ITS_Pelliccioni.py ./dose_vs_E_kapton_150um.txt ./Pelliccioni2000-Table_A2.2_ElectronDoseEquivPerFluence.txt

ITS_dose_vs_E_filename             = sys.argv[1]
Pelliccioni2000_dose_vs_E_filename = sys.argv[2]

E_ITS, dose_ITS = np.loadtxt(ITS_dose_vs_E_filename, unpack=True) 
E_P2000, dose_P2000 = np.loadtxt(Pelliccioni2000_dose_vs_E_filename, unpack=True)


plt.loglog(E_ITS*1E6, dose_ITS, marker='.', label='ITS', markersize=15)
plt.loglog(E_P2000, dose_P2000, marker='^', label='Pelliccioni 2000 Table A2.2')
plt.xlabel('Energy (eV)')
plt.ylabel(r'Dose per Fluence (Sv / (#/cm$^2$))')
plt.grid(b=True, which='both') # https://stackoverflow.com/questions/9127434/how-to-create-major-and-minor-gridlines-with-different-linestyles-in-python
plt.legend()
plt.savefig("comparison_dose_vs_energy.png", dpi=1000)
plt.show()