# $CVSHeader: pal8e/dat/stat/plot-Qbond,v 1.2 2005/02/15 01:10:24 indraneel Rel $

## Run using GNUPlot

set data style lines
set mxtics 5
show mxtics
set mytics 10
show mytics
#set key below box
set key left

#a == amplitude
#s == standard deviation
#c == center

#fs(x) = as * exp (-0.5 * (((x-cs)/ss)**2))

#fa1(x) = aa1 * exp (-0.5 * (((x-ca1)/sa)**2))
#fa2(x) = aa2 * exp (-0.5 * (((x-ca2)/sa)**2))
#fa(x) = fa1(x) + fa2(x)

#as = 0.170
#cs = 4.810
#ss = 0.227

#aa1 = 0.073
#aa2 = 0.077
#ca1 = 4.463
#ca2 = 5.245
#sa = 0.261

#fit [3.2:7.2] fs(x) "QdistSyn" using 1:2 via as, ss, cs
#fit [3.2:7.2] fa(x) "QdistAnti" using 1:2 via aa1, aa2, ca1, ca2, sa

set xrange [0:90]
#set yrange [0:0.2]
set xlabel "Angle between lines for Quadruplet parameter 2 ({/Symbol \260})"
set ylabel "Normalised Count"
set title "C{/Symbol a} Angles for Quadruplet Parameter 2 from DSSP"

set arrow from 87.0,graph 0 to 87.0,graph 1 nohead lt 1 lw 4
set label "c1" at 83,graph 0.96
set arrow from 82.2,graph 0 to 82.2,graph 1 nohead lt 1 lw 4
set label "c2" at 78,graph 0.96

set terminal x11
plot "QbondSyn" with lines lt 2 lw 4 title "Parallel strands", \
	"QbondAnti" with lines lt 3 lw 4 title "Antiparallel strands"

#set terminal postscript enhanced color blacktext solid "Times-Roman" 12
set terminal postscript eps enhanced "Times-Roman" 18
set output "Qbond.eps"
replot


