unset multiplot set term pdf fontscale 0.6 size 4.95in,3.3in norm = 1 all = 1 if(norm == 0){outfile = 'compiled_data-Ir-Zr-Ru-Db-Sn-Nb-Ge-Cm-Pr_G.pdf'} if(norm == 1){outfile = 'compiled_data-Ir-Zr-Ru-Db-Sn-Nb-Ge-Cm-Pr_norm.pdf'} if(norm == 1 && all == 1){outfile = 'compiled_data-Ir-Db_mod.pdf'} set output outfile column1 = 11 # kicker 1 set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6' set label 'plotting_scripts/compile_data.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance nturns = 10000 set yrange [0:3] set ylabel "captured [%]" set xlabel 'kick [G]' if(norm == 1){set xlabel 'Blumlein voltage [arb]'; set xrange [0.6:2]; set xtics 0.4,0.2} set label at graph 0.1,0.95 'Inflector tilt = 0, Injection angle = 0' dir1 = 'Ge' N=2 dir1= 'Cm' # Cm = 228.59 f1=228.59 dir2 = 'Db' Db=201.95 f2=201.95 dir3 = 'Ge' Ge=216.94 f3=216.94 dir4 = 'Ir' Ir=194.62 dir5 = 'Nb' Nb=226.11 dir6='Pr' Pr=228.53 dir7='Ru' Ru=194.63 dir8='Sn' Sn=197.99 dir9='Zr' Zr=204.92 dir10='Db_mod' Db_mod=201.95 f1='Cm:CapMatch+2sections' f2='Db:+2sections' f3 = 'Ge:PerfectMatch' f4 = 'Ir:AsBuilt' f5 = 'Nb:CapMatch' f6='Pr:CapMatch+3sections' f7='Ru:StaggerCabLengths' f8='Sn:FourKickPlates' f9='Zr:+4sections' f10='Db_mod' B0=292 B0=1. if(norm == 0){Cm=1;Db=1;Ge=1;Ir=1;Nb=1;Pr=1;Ru=1;Sn=1;Zr=1;Db_mod=1.} filename1 = dir1.'/compilation_norm.dat' filename2 = dir2.'/compilation_norm.dat' filename3 = dir3.'/compilation_norm.dat' filename4 = dir4.'/compilation_norm.dat' filename5 = dir5.'/compilation_norm.dat' filename6 = dir6.'/compilation_norm.dat' filename7 = dir7.'/compilation_norm.dat' filename8 = dir8.'/compilation_norm.dat' filename8 = dir8.'/20180517_093852/VparamDependence.dat' filename8 = dir8.'/20180517_123525/VparamDependence.dat' filename8 = dir8.'/20180518_064219/VparamDependence.dat' filename9 = dir9.'/compilation_norm.dat' filename10=dir10.'/20180519_111400/VparamDependence.dat' set key left Left bottom font ',8' noenhance set yrange [0:3] if(norm == 0){set xrange [120:360]} #set xrange [0.5:1.9] #plot for [i=1:N] u ($11*10000):($2 == nturns ? ($18*100/$2):1/0) w p ps 0.5 pt 5 t dir(i) plot \ filename4 u ($11*10000/Ir*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 t f4 ,\ filename10 u ($19/Sn*B0):($16*100/$14) w lp ps 0.5 pt 9 t f10 #,\ #filename9 u ($11*10000/Zr*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 6 t f9 ,\ #filename7 u ($11*10000/Ru*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 t f7 ,\ #filename2 u ($11*10000/Db*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 t f2 ,\ #filename8 u ($19/Sn*B0):($16*100/$14) w lp ps 0.5 pt 7 lc 6 t f8 ,\ #filename5 u ($11*10000/Nb*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 lc 5 t f5 ,\ #filename3 u ($11*10000/Ge*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 lc 7 t f3 ,\ #filename1 u ($11*10000/Cm*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 5 lc 8 t f1 ,\ #filename6 u ($11*10000/Pr*B0):($2 == nturns ? ($18*100/$2):1/0) w lp ps 0.5 pt 8 lc 9 t f6 reset set xlabel 'kick [G]' if(norm == 0){set xrange [120:360]} if(norm == 1){set xlabel 'Blumlein voltage [arb]'; set xrange [0.6:2]; set xtics 0.4,0.2} set ylabel "{/Symbol s}_p [%]" set label at graph 0.1,0.05 'Inflector tilt = 0, Injection angle = 0' set key Left left font ',8' noenhance set yrange [0.1:0.18] #set xrange [0.5:1.9] plot \ filename4 u ($11*10000/Ir*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 t f4 ,\ filename10 u ($19/Sn*B0):($13*100) w lp ps 0.5 pt 9 t f10 #,\ #filename9 u ($11*10000/Zr*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 6 t f9 ,\ #filename7 u ($11*10000/Ru*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 t f7 ,\ #filename2 u ($11*10000/Db*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 t f2 ,\ #filename8 u ($19/Sn*B0):($13*100) w lp ps 0.5 pt 7 lc 6 t f8 ,\ #filename5 u ($11*10000/Nb*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 lc 5 t f5 ,\ #filename3 u ($11*10000/Ge*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 t f3 ,\ #filename1 u ($11*10000/Cm*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 5 t f1 ,\ #filename6 u ($11*10000/Pr*B0):($2 == nturns ? ($24*100):1/0) w lp ps 0.5 pt 8 t f6 reset set xlabel 'kick [G]' if(norm == 0){set xrange [120:360]} if(norm == 1){set xlabel 'Blumlein voltage [arb]'; set xrange [0.6:2]; set xtics 0.4,0.2} set ylabel "<{/Symbol D}p/p> [%]" set label at graph 0.1,0.05 'Inflector tilt = 0, Injection angle = 0' set key right Left font ',8' noenhance set yrange [-0.06:0.14] #set xrange [0.5:1.9] plot \ filename4 u ($11*10000/Ir*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 t f4 ,\ filename10 u ($19/Sn*B0):($35*100) w lp ps 0.5 pt 9 t f10 #,\ #filename9 u ($11*10000/Zr*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 6 t f9 ,\ #filename7 u ($11*10000/Ru*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 t f7 ,\ #filename2 u ($11*10000/Db*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 t f2 ,\ #filename8 u ($19/Sn*B0):($35*100) w lp ps 0.5 pt 7 lc 6 t f8 ,\ #filename5 u ($11*10000/Nb*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 lc 5 t f5 ,\ #filename3 u ($11*10000/Ge*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 t f3 ,\ #filename1 u ($11*10000/Cm*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 5 t f1 ,\ #filename6 u ($11*10000/Pr*B0):($2 == nturns ?($25*100):1/0) w lp ps 0.5 pt 8 t f6 reset set xlabel 'kick [G]' set ylabel "CBO p-p [mm]" set key left Left bottom font ',8' noenhance if(norm == 0){set xrange [120:360]} if(norm == 1){set xlabel 'Blumlein voltage [arb]'; set xrange [0.6:2]; set xtics 0.4,0.2; set key right Left bottom font ',8'} set yrange [0:30] #set xrange [0.5:1.9] set label at graph 0.2,0.95 'Inflector tilt = 0, Injection angle = 0' plot \ filename4 u ($11*10000/Ir*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 t f4 ,\ filename10 u ($19/Sn*B0):($4) w lp ps 0.5 pt 9 t f10 #,\ #filename9 u ($11*10000/Zr*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 6 t f9 ,\ #filename7 u ($11*10000/Ru*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 t f7 ,\ #filename2 u ($11*10000/Db*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 t f2 ,\ #filename8 u ($19/Sn*B0):($4) w lp ps 0.5 pt 7 lc 6 t f8 ,\ #filename5 u ($11*10000/Nb*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 lc 5 t f5 ,\ #filename3 u ($11*10000/Ge*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 t f3 ,\ #filename1 u ($11*10000/Cm*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 5 t f1 ,\ #filename6 u ($11*10000/Pr*B0):($2 == nturns ?($26*1000):1/0) w lp ps 0.5 pt 8 t f6 print outfile