tips

# 1-D
cat outFile.csv | awk -F, 'BEGIN{pi=4*atan2(1,1)}($1!="#"){phi1=$3;phi2=$5;dphi=phi2+phi1;if(dphi>pi)dphi-=2*pi;if(dphi<-pi)dphi+=2*pi;print(1,dphi*180/pi)}' | whist.pl -200 200 400 >x.hist
gnuplot> plot "x.hist" with steps

# 2-D
cat outFile.csv | awk -F, 'BEGIN{pi=4*atan2(1,1)}($1!="#"){phi1=$3;phi2=$5;print(1,phi1*180/pi,phi2*180/pi)}' | whist.pl -200 200 40 -200 200 40 >x.hist
gnuplot> splot "x.hist" with steps

cat outFile.csv | awk -F, 'BEGIN{pi=4*atan2(1,1)}($1!="#"){theta1=$2;phi1=$3;theta2=$4;phi2=$5;dphi=phi2+phi1;if(dphi>pi)dphi-=2*pi;if(dphi<-pi)dphi+=2*pi;if(theta1*180/pi>80&&theta1*180/pi<90&&theta2*180/pi>80&&theta2*180/pi<90)print(1,dphi*180/pi)}' | whist.pl -200 200 400 >x.hist