split_edge.py

import pylab,numpy
from matplotlib.font_manager import fontManager, FontProperties


#set fonts here
font= FontProperties(size='small')
params = {'backend': 'ps',
          'axes.labelsize': 16,
          'xtick.labelsize': 15,
          'ytick.labelsize': 15,
          'legend.fontsize': 16,
          #'font.family': 'sans',
          'font.serif': ['Times','Palatino','serif'],
          'font.size': 16,
          'ps.usedistiller': 'xpdf',  
          }
pylab.rcParams.update(params)


#set figure size
csize=6 #inches = 2.95 used in chemosphere
fig_size=(csize,csize*1.618)
#adjust spacing


#yt = [0, 500, 1500, 2500, 3500]  #y axis ticks

fig = pylab.figure()

fig.subplots_adjust(top=0.95,bottom=0.1,left=0.10,right=0.93,hspace=0.35,wspace=0.2)

# this is plot 1 with Procut orange data


ax = fig.add_subplot(1,1,1)

diameters = [ 1, 2, 3, 4, 5 ]
steps =[ '1', '2', '3', '4', '5' ]
types = [ '$S^{3}/Q$',  '$L(4,1)$' ,'$L(3,1)$' ,'$RP^{3}$' ,'$S^{3}$' ]
markers = ['ws','wo','bx','bs','ro']
msdef = 9
markersizes = [msdef,msdef,12,msdef,msdef]
mewdef = 1
markeredgewidths = [mewdef,mewdef,1.1,mewdef,mewdef]


countd = {}
countd[ '$S^{3}/Q$'] = [ 0, 1, 0, 0, 0]
countd[ '$L(4,1)$' ] = [ 0, 0, 1, 0, 0]
countd[ '$L(3,1)$' ] = [ 0, 1, 1, 0, 0]
countd[ '$RP^{3}$' ] = [ 0,17, 5, 0, 0]
countd[ '$S^{3}$'  ] = [ 3, 574, 3498, 675, 11 ]


#ax.plot(diameters,counts,'k--',alpha=0.5)
for d,c,s in zip(diameters,countd[ '$S^{3}$'  ],steps):
    ax.plot([d,d],[0.7,c],'k--')

for t,m,ms,mew in zip(types,markers,markersizes,markeredgewidths):

    ax.semilogy(diameters,countd[t],m,ms=ms,mew=mew,label=t)

ax.set_xlabel('Edge-Diameter')
#ax.xaxis.set_label_coords(0.5, -0.12)
ax.set_ylabel('Count (log scale)')
ax.set_xticks([0] + diameters + [6])
#ax.set_yticks(yt)
ax.set_xticklabels( ['0','1','2','3','4','5','6'], horizontalalignment='center')
a=ax.axis([0,6,0.7,3498*5])

ax.legend(numpoints=1)


#ax.legend( ('Procut Orange',), loc=1 ) #add trailing comma to make it a tuple so full word is printed



pylab.savefig('./paper/figures/split_edge_diameters.png', dpi=300)
pylab.show()