# This script needs gnuplot3.6beta or later. set term post port enh "Helvetica" 17 #set term post enh "Helvetica" 15 #set term post eps enh "Helvetica" 15 set output 'fig2.ps' # Multiplicity Distibution # Generated by CASC # Bimp f(b) Nc_i Nimf_i Nc_f Nimf_f # No Fil W.Fil No Fil W.Fil No Fil W.Fil No Fil W.Fil # 0.50 78.54 97.50 91.00 6.00 5.00 100.00 93.50 5.00 4.00 #ymx=11.5 ymx=12 xsft=12.0 xsft0=xsft*0 xsft1=xsft*1 xsft2=xsft*2 set xrange [0:xsft*3] set yrange [0:ymx*2] set trange [0:1.0] #set fixedb set nokey #set xlabel "b (fm)" set label 1 "b (fm)" at 6+xsft0,-1.5 center set label 2 "b (fm)" at 6+xsft1,-1.5 center set label 3 "b (fm)" at 6+xsft2,-1.5 center set label 4 "E/A = 100 MeV" at 6+xsft0,24.5 center set label 5 "E/A = 250 MeV" at 6+xsft1,24.5 center set label 6 "E/A = 400 MeV" at 6+xsft2,24.5 center set label 7 "Before SDM" at 3.5+xsft2,11+ymx set label 8 "After SDM" at 3.5+xsft2,11 set ylabel "" 1,0 set pointsize 1.4 set bar small set ticscale 0.5 xrat=0.2 yrat=0.3 #set linewidth 0.5 #set linewidth 0.2 set linestyle 1 lt 1 lw 3.0 set linestyle 3 lt 3 lw 3.0 # set size 1.0,0.75 set origin 0.0,0.0 #set title "E/A = 100 MeV" set ytics (\ '0' 0, '2' 2, '4' 4\ ,'6' 6, '8' 8, '10' 10\ ,'0' 0+ymx, '2' 2+ymx, '4' 4+ymx\ ,'6' 6+ymx, '8' 8+ymx, '10' 10+ymx\ ) #set xtics 0,2,10 set xtics (\ '0' 0+xsft0,'2' 2+xsft0,'4' 4+xsft0\ ,'6' 6+xsft0,'8' 8+xsft0,'10' 10+xsft0\ ,'0' 0+xsft1,'2' 2+xsft1,'4' 4+xsft1\ ,'6' 6+xsft1,'8' 8+xsft1,'10' 10+xsft1\ ,'0' 0+xsft2,'2' 2+xsft2,'4' 4+xsft2\ ,'6' 6+xsft2,'8' 8+xsft2,'10' 10+xsft2\ ) set parametric plot \ 'e100.cnv' u ($3+xsft0):($4):8 w err 1 13\ ,'e100.cnv-10' u ($3+xsft0):($4) w l 3 1\ ,'e100q.n' u ($1+xsft0):($9) w linesp 1 6\ ,'e100c.n' u ($1+xsft0):($9) w linesp 1 8\ ,'e100.cnv' u ($3+xsft0):($4+ymx):8 w err 1 13\ ,'e100.cnv-10' u ($3+xsft0):($4+ymx) w l 3 1\ ,'e100q.n' u ($1+xsft0):($5+ymx) w linesp 1 6\ ,'e100c.n' u ($1+xsft0):($5+ymx) w linesp 1 8\ ,'tics.dat' u ($1+xsft0):(xrat*$2+ymx) w l 1\ ,'e250.cnv' u ($3+xsft1):($4):8 w err 1 13\ ,'e250.cnv-10' u ($3+xsft1):($4) w l 3 1\ ,'e250q.n' u ($1+xsft1):($9) w linesp 1 6\ ,'e250c.n' u ($1+xsft1):($9) w linesp 1 8\ ,'e250.cnv' u ($3+xsft1):($4+ymx):8 w err 1 13\ ,'e250.cnv-10' u ($3+xsft1):($4+ymx) w l 3 1\ ,'e250q.n' u ($1+xsft1):($5+ymx) w linesp 1 6\ ,'e250c.n' u ($1+xsft1):($5+ymx) w linesp 1 8\ ,'tics.dat' u ($1+xsft1):(xrat*$2+ymx) w l 1\ ,'e400.cnv' u ($3+xsft2):($4):8 w err 1 13\ ,'e400.cnv-10' u ($3+xsft2):($4) w l 3 1\ ,'e400q.n' u ($1+xsft2):($9) w linesp 1 6\ ,'e400c.n' u ($1+xsft2):($9) w linesp 1 8\ ,'e400.cnv' u ($3+xsft2):($4+ymx):8 w err 1 13\ ,'e400.cnv-10' u ($3+xsft2):($4+ymx) w l 3 1\ ,'e400q.n' u ($1+xsft2):($5+ymx) w linesp 1 6\ ,'e400c.n' u ($1+xsft2):($5+ymx) w linesp 1 8\ ,'xtics.dat' u ($1+xsft2):(xrat*$2+ymx) w l 1\ ,'ytics.dat' u ($2*yrat+xsft1):($1+ymx) w l 1\ ,'ytics.dat' u ($2*yrat+xsft1):($1) w l 1\ ,'ytics.dat' u ($2*yrat+xsft2):($1+ymx) w l 1\ ,'ytics.dat' u ($2*yrat+xsft2):($1) w l 1\ ,t*xsft*3,ymx w l 1\ ,xsft1,t*ymx*2 w l 1\ ,xsft2,t*ymx*2 w l 1\