#!/usr/bin/env python
import matplotlib.pyplot as plt
import matplotlib as mpl
import numpy as np

ni_LO = np.array(
[[101.1,	0.2143,	0.1097],
[112.3,	0.2143,	0.0775],
[121.2,	0.1250,	0.0523]]
#[130.7,	0.1667,	0.0359],
#[136.4,	0.1452,	0.0447],
#[146.6, 0.236, 0.041],
#[160.4,	0.1357,	0.0289]]
)

ni_v1 = np.array(
[[71.4, 1.000, 0.034],
[89.2, 0.867, 0.062],
[101.1, 0.822, 0.057],
[112.3,	0.7885,	0.0566],
[121.2,	0.7500,	0.0472],
[130.7,	0.7451,	0.0432],
[136.4,	0.7016,	0.0411],
[146.6, 0.658, 0.037],
[160.4,	0.5317,	0.0349]]
)

pd_LO = np.array(
[[89.2, 0.091, 0.087],
[102.4, 0.227, 0.089],
[111.8, 0.059, 0.040],
[120.0, 0.212, 0.071]]
)

pd_v1 = np.array(
[[89.2, 0.926, 0.050],
[97.4, 0.840, 0.052],
[102.4, 0.922, 0.025],
[111.8, 0.816, 0.036],
[120.0, 0.741, 0.049]]
)

pt_LO = np.array(
[[81.7, 0.261, 0.092],
[89.3, 0.367, 0.088],
[97.4, 0.200, 0.073],
[102.5, 0.116, 0.049],
[111.9, 0.191, 0.057],
[120.1, 0.171, 0.045]]
)

pt_v1 = np.array(
[[60.7, 0.880, 0.065],
[71.4, 0.833, 0.076],
[81.9, 0.606, 0.085],
[92.2, 0.643, 0.074],
[104.6, 0.667, 0.057]]
)

mpl.rc("font", size=10)
fig = plt.figure(figsize=(3.69,3))

plt.errorbar( ni_LO[:,0], ni_LO[:,1], yerr=ni_LO[:,2], marker='o', label='Ni(111)', color='b' )
plt.errorbar( ni_v1[:,0], ni_v1[:,1], yerr=ni_v1[:,2], marker='o', color='b', linestyle='--' )
plt.errorbar( pd_LO[:,0], pd_LO[:,1], yerr=pd_LO[:,2], marker='o', label='Pd(111)', color='k' )
plt.errorbar( pd_v1[:,0], pd_v1[:,1], yerr=pd_v1[:,2], marker='o', color='k', linestyle='--' )
plt.errorbar( pt_LO[:,0], pt_LO[:,1], yerr=pt_LO[:,2], marker='o', label='Pt(111)', color='r' )
plt.errorbar( pt_v1[:,0], pt_v1[:,1], yerr=pt_v1[:,2], marker='o', color='r', linestyle='--' )

plt.legend(loc='best', numpoints=1, fontsize='medium')

plt.xlim(50,170.)
plt.ylim(0,1.)
plt.tick_params(length=6, width=1)
plt.ylabel('Fraction of CH cleavage')
plt.xlabel('Incidence energy (kJ/mol)')

plt.tight_layout()
plt.savefig('cleavage_comparison.pdf')
#plt.show()