Back to home page

EIC code displayed by LXR

 
 

    


Warning, /geant4/examples/extended/medical/dna/chem5/plot.gp is written in an unsupported language. File is not indexed.

0001 set term postscript eps enhanced color "Helvetica" 46
0002 set output "plot.eps"
0003 set size 3.5,3.4
0004 
0005 set style line 1 dt 1 pt 5 lc 0 lw 10.0 ps 3.0 #1.5
0006 set style line 2 dt 2 pt 5 lc 0 lw 10.0 ps 3.0 #1.5
0007 set style line 3 dt 3 pt 5 lc 0 lw 10.0 ps 3.0 #1.5
0008 set style line 4 dt 4 pt 5 lc 0 lw 10.0 ps 3.0 #1.5
0009 set style line 5 dt 4 pt 5 lc 3 lw 10.0 ps 3.0 #1.5
0010 
0011 fileName1="Species.txt"
0012 
0013 set logscale x
0014 set format x "10^{%L}"
0015 set xrange[1.01:0.999e6]
0016 
0017 set xlabel "Time (ps)" font "Helvetica, 46" tc rgb "black" #white
0018 set ylabel "G(Species/100 eV)" tc lt 0 font "Helvetica, 46" offset -0, -8
0019 
0020 set multiplot
0021 ############### OH ##################
0022 set size 1.4,1.7
0023 set origin 0.0, 1.7
0024 set rmargin 0
0025 set lmargin 10
0026 set bmargin 0
0027 set tmargin 3
0028 
0029 set format x ""
0030 set xlabel ""
0031 set yrange[2.51:6]
0032 #^{\267}OH
0033 plot fileName1 u 1:(stringcolumn(4) eq "OH^0" ? $2 : 1/0) w l ls 1 title "Default",\
0034      "data/OH.txt" index 0 u 1:2:3 with errorbars pt 4 ps 4 lc 0 notitle,\
0035      "data/OH.txt" index 1 u 1:2:3 with errorbars pt 5 ps 4 lc 0 notitle,\
0036      "data/OH.txt" index 2 u 1:2:3 with errorbars pt 6 ps 4 lc 0 notitle
0037 
0038 ############### e_aq #################
0039 set size 1.4,1.7
0040 set origin 1.4, 1.7
0041 set lmargin 0
0042 set rmargin 10
0043 set bmargin 0
0044 set tmargin 3
0045 
0046 set format y ""
0047 set ylabel ""
0048 
0049 plot fileName1 u 1:(stringcolumn(4) eq "e_aq^-1" ? $2 : 1/0) w l ls 1 title "e-_{aq}",\
0050      "data/e_aq.txt" index 0 u 1:2:3 w errorbars pt 4 ps 4 lc 0 notitle,\
0051      "data/e_aq.txt" index 1 u 1:2:3 w errorbars pt 5 ps 4 lc 0 notitle,\
0052      "data/e_aq.txt" index 2 u 1:2:3 w errorbars pt 6 ps 4 lc 0 notitle,\
0053      "data/e_aq.txt" index 3 u 1:2:3 w errorbars pt 7 ps 4 lc 0 notitle
0054      
0055 
0056 ############### H3O #################
0057 
0058 set size 1.4,1.7
0059 set origin 2.415, 1.7
0060 set lmargin 0
0061 set rmargin 10
0062 set bmargin 0
0063 set tmargin 3
0064 
0065 plot fileName1 u 1:(stringcolumn(4) eq "H3O^1" ? $2 : 1/0) w l ls 1 title "H_{3}O^{+}",\
0066 
0067 ############### H2O2 #################
0068 
0069 set size 1.4,1.7
0070 set origin 0.0, 0.0
0071 set tmargin 0
0072 set lmargin 10
0073 set rmargin 0
0074 unset bmargin
0075 
0076 unset format y
0077 set format x "10^{%L}"
0078 set xlabel "Time (ps)" font "Helvetica, 46" tc rgb "black" #white
0079 
0080 set yrange[0:0.99]
0081 set ytics 0, 0.2, 1.
0082 
0083 plot fileName1 u 1:(stringcolumn(4) eq "H2O2^0" ? $2 : 1/0) w l ls 1 title "H_{2}O_{2}",\
0084 
0085 ############### H2 #################
0086 
0087 set size 1.4,1.7
0088 set origin 1.4, 0.0
0089 set tmargin 0
0090 set lmargin 0
0091 set rmargin 10
0092 unset bmargin
0093 
0094 set format y ""
0095 set ylabel ""
0096 
0097 plot fileName1 u 1:(stringcolumn(4) eq "H_2^0" ? $2 : 1/0) w l ls 1 title "H_{2}",\
0098 
0099 ############### H2 #################
0100 
0101 set size 1.4,1.7
0102 set origin 2.415, 0.0
0103 set tmargin 0
0104 set lmargin 0
0105 set rmargin 10
0106 unset bmargin
0107 
0108 plot fileName1 u 1:(stringcolumn(4) eq "H^0" ? $2 : 1/0) w l ls 1 title "H^{\267}",\
0109 
0110 set nomultiplot