#!/usr/bin/env python
import matplotlib.pyplot as plt
import matplotlib as mpl
import numpy as np
from math import erf
from sticking_data import *

mpl.rc("text", usetex=True)
plt.style.use('tableau-colorblind10')
Nrows = 1
Ncols = 1
fig, ax = plt.subplots(nrows=Nrows, ncols=Ncols, figsize=(3.69,3.))

ax1 = plt.subplot(Nrows,Ncols,1)

plt.errorbar( AIMD_LO_CHD3_migliorini_2017[:,0], AIMD_LO_CHD3_migliorini_2017[:,1], yerr=AIMD_LO_CHD3_migliorini_2017[:,2], marker='o', label=r'DFT', capsize=4 )
plt.errorbar( QCT_LO_CHD3_migliorini_2017[:,0], QCT_LO_CHD3_migliorini_2017[:,1], yerr=QCT_LO_CHD3_migliorini_2017[:,2], marker='o', label=r'HDNNP', capsize=4 )
plt.errorbar( Exp_LO_CHD3_migliorini_2017[:,0], Exp_LO_CHD3_migliorini_2017[:,1], yerr=Exp_LO_CHD3_migliorini_2017[:,2], marker='d', label=r'Exp.', capsize=4 )

plt.legend()
plt.tick_params(length=6, width=1, direction='in', top=True, right=True)
plt.ylabel('Sticking probability')
plt.xlabel('Incidence energy (kJ/mol)')

plt.tight_layout()
#plt.subplots_adjust(wspace=0, hspace=0)
plt.savefig('stickingprobability_AIMD.pdf')
#plt.show()