#from the run/ directory
#make a histogram of run lengths for all exp, and by q2 setting
awk 'BEGIN{BIN_SIZE=120} {bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times.dat | sort -n > all

awk 'BEGIN{BIN_SIZE=120} $1<=41611{bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times.dat | sort -n > q05

awk 'BEGIN{BIN_SIZE=120} $1>=41794{bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times.dat | sort -n > q10

graph -m 1 -S 0 -m 2 q05 -m 4 q10 -m 1 -W 0.005 all -L "Gen01 Runs" -X "Length (sec)" -Tps > length_histogram.ps

#make a histogram of old run lengths for all exp, and by q2 setting
awk 'BEGIN{BIN_SIZE=120} {bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times_old.dat | sort -n > all

awk 'BEGIN{BIN_SIZE=120} $1<=41611{bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times_old.dat | sort -n > q05

awk 'BEGIN{BIN_SIZE=120} $1>=41794{bins[int(($3-$2)/BIN_SIZE)]++} END{for(i in bins) print i*BIN_SIZE+BIN_SIZE/2,bins[i]}' run_times_old.dat | sort -n > q10

graph -m 1 -S 0 -m 2 q05 -m 4 q10 -m 1 -W 0.005 all -L "Gen01 Runs" -X "Length (sec)" -Tps > length_histogram_old.ps

#from the data/ directory
#compare different versions of the polarization file
join target_pol_v{6,7}.dat -o 1.1,1.11,2.13 | awk '$2!=0{print $1,$3/$2-1}' | g41 -m 0 -L "Ver 6 to Ver 7" -Tps > ver6_ver7.ps

join target_pol_v{7,8}.dat -o 1.1,1.13,2.13 | awk '$2!=0{print $1-40000,$3/$2-1}' | g41 -m 0 -L "Ver 7 to Ver 8" -y -.2 .2 -Y "(new - old)/old" -X "Run Number - 40000" -r .2 -Tps >ver7_ver8.ps

#look at possilbe problem runs
join target_pol_v6.dat ../runpol/target_pol.dat -o 1.1,1.2,1.5,2.2,2.5,1.11,2.13 | awk '$6!=0{printf"%5d  %5d   | %4d %6.2f | %4d %6.2f |\n",$1,($5-$4)-($3-$2),($3-$2),$6,($5-$4),$7}' | sort -n -k 2 | awk 'BEGIN{print"               |   Ver 6     |   Ver 10    |";print"Run #  Len Chg |  Len   Pol  |  Len   Pol  |";print"---------------+-------------+-------------+"} $2>=4000 && $2<5000' | mpage -1 -o - > big_change_runs.ps