#!/usr/bin/env python2.7
# -*- coding: utf8 -*-
"""
Purpose:
File: /home/cz/Program/eclipse/GP/efficiency_factor.py
Author: CHANG Zhi
Date: 2014-12-24
E-mail:
"""
import os
import numpy as np
import matplotlib.pyplot as plt
datafile = 'data/abc.dat'
if(os.path.isfile(datafile)):
cat = np.genfromtxt(datafile,delimiter="")
bars,factor,evts= cat[:,0],cat[:,1],cat[:,2]
else:
print datafile,' do not exist.\n'
exit(0)
#host = plt.subplots()
fig,left_axis=plt.subplots()
#fig.subplots_adjust(right_axis=0,75)
right_axis = left_axis.twinx()
p1, = left_axis.plot(bars, abs(factor), 'b.-')
p2, = right_axis.plot(bars, evts/(5000000*(3*np.sqrt(3)/(2*np.pi))), 'r.-')
left_axis.set_xlim(40,410)
#left_axis.set_xticks(np.arange(0,110,10))
left_axis.set_ylim(0.265,0.355)
left_axis.set_yticks(np.arange(0.27,0.36,0.01))
#right_axis.set_ylim(0,0.26)
#right_axis.set_yticks(np.arange(0,0.26,0.02))
left_axis.set_xlabel('Bar Number in One Edge')
left_axis.set_ylabel('Modulation Factor')
right_axis.set_ylabel('Efficiency')
left_axis.yaxis.label.set_color(p1.get_color())
right_axis.yaxis.label.set_color(p2.get_color())
tkw = dict(size=5, width=1.5)
left_axis.tick_params(axis='y', colors=p1.get_color(), **tkw)
right_axis.tick_params(axis='y', colors=p2.get_color(), **tkw)
left_axis.tick_params(axis='x', **tkw)
plt.savefig('fig/abc.eps')
plt.show()