# $CVSHeader: pal8e/dat/stat/plot-QtorTS,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 5
show mytics
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 [-5:3]
#set yrange [0:0.2]
set xlabel "Times Sigma"
set ylabel "Probability"
set title "Probability for Quadruplet Torsion Angles (Parameter 3) from DSSP"

set terminal x11
plot "QtorSynTS" with lines lt 2 lw 4 title "Parallel strands", \
	"QtorAntiTS" 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 "QtorTS.eps"
replot


