#!/usr/bin/gnuplot -persist

load "~/Script/gnuplot/ps_landscape.gpi"



#########################################################################################
# Shift versus position

set key bottom left box

#__________________________________

set title "X shift (Measured position - hole position) versus X measured, B2 = 18000 DAC"
set xlabel "X measured [cm]"
set ylabel "X shift [cm]"

set output "ShiftVsPos_Temp_u34.ps"

plot "ShiftVsPos_Temp.dat" index 4 using 3:4 ps 2 pt 3 lw 2 t "Upstream collimator", "ShiftVsPos_Temp.dat" index 5 using 3:4 ps 2 pt 8 lw 2 t "Downstream collimator"


set key screen 0.55, screen 0.5

set title "X shift (Measured position - hole position) versus Y measured, B2 = 18000 DAC"
set xlabel "Y measured [cm]"
set ylabel "X shift [cm]"

set output "ShiftVsPos_Temp_u64.ps"

# Upstream collimator rotation around Z
f(x) = 0.16076
g(x) = -0.24

set label "y = 0.16076" at screen 0.4, screen 0.77
set label "y = -0.24" at screen 0.4, screen 0.25

plot "ShiftVsPos_Temp.dat" index 4 using 6:4 ps 2 pt 3 lw 2 t "Upstream collimator", f(x) t "Upstream collimator rotation" lw 2, "ShiftVsPos_Temp.dat" index 5 using 6:4 ps 2 pt 8 lw 2 t "Downstream collimator", g(x) t "Upstream collimator rotation" lw 2

unset label
set key bottom left

set title "Y shift (Measured position - hole position) versus Y measured, B2 = 18000 DAC"
set xlabel "Y measured [cm]"
set ylabel "Y shift [cm]"

set output "ShiftVsPos_Temp_u67.ps"

plot "ShiftVsPos_Temp.dat" index 4 using 6:7 ps 2 pt 3 lw 2 t "Upstream collimator", "ShiftVsPos_Temp.dat" index 5 using 6:7 ps 2 pt 8 lw 2 t "Downstream collimator"


set key screen 0.57, screen 0.45

set title "Y shift (Measured position - hole position) versus X measured, B2 = 18000 DAC"
set xlabel "X measured [cm]"
set ylabel "Y shift [cm]"

set output "ShiftVsPos_Temp_u37.ps"

# Upstream collimator rotation around Z
f(x) = 0.10857
g(x) = 0.015

set label "y = 0.10857" at screen 0.8, screen 0.73
set label "y = 0.015" at screen 0.8, screen 0.54

plot "ShiftVsPos_Temp.dat" index 4 using 3:7 ps 2 pt 3 lw 2 t "Upstream collimator", f(x) t "Upstream collimator rotation" lw 2, "ShiftVsPos_Temp.dat" index 5 using 3:7 ps 2 pt 8 lw 2 t "Downstream collimator", g(x) t "Downstream collimator rotation" lw 2

unset label

#########################################################################################
# Shift versus B2

set key screen 0.92, screen 0.889
set rmargin 16

#__________________________________

set output "ShiftVsB2_Temp_idx_0_2.ps"

set title "Upstream collimator, X shift (Measured position - hole position) versus B2"
set xlabel "B2 [DAC]"
set ylabel "X shift [cm]"

plot "ShiftVsB2_Temp.dat" index 0 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 0 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 2 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints

set output "ShiftVsB2_Temp_idx_1_3.ps"

set title "Upstream collimator, Y shift (Measured position - hole position) versus B2"
set xlabel "B2 [DAC]"
set ylabel "Y shift [cm]"

plot "ShiftVsB2_Temp.dat" index 1 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 1 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 3 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_4_6.ps"

set title "Downstream collimator, X shift (Measured position - hole position) versus B2"
set xlabel "B2 [DAC]"
set ylabel "X shift [cm]"

plot [:][-0.8:0.6] "ShiftVsB2_Temp.dat" index 4 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 4 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 6 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints

set output "ShiftVsB2_Temp_idx_5_7.ps"

set title "Downstream collimator, Y shift (Measured position - hole position) versus B2"
set xlabel "B2 [DAC]"
set ylabel "Y shift [cm]"

plot "ShiftVsB2_Temp.dat" index 5 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 5 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 7 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_8_9.ps"

set title "Upstream collimator, ring heights versus B2"
set xlabel "B2 [DAC]"
set ylabel "Height [counts]"

plot "ShiftVsB2_Temp.dat" index 8 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 8 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 9 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_10_11.ps"

set title "Downstream collimator, ring heights versus B2"
set xlabel "B2 [DAC]"
set ylabel "Height [counts]"

plot "ShiftVsB2_Temp.dat" index 10 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 10 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 11 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_12_13.ps"

set title "Upstream collimator, peak heights versus B2"
set xlabel "B2 [DAC]"
set ylabel "Height [counts]"

plot "ShiftVsB2_Temp.dat" index 12 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 12 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 13 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_14_15.ps"

set title "Downstream collimator, peak heights versus B2"
set xlabel "B2 [DAC]"
set ylabel "Height [counts]"

plot "ShiftVsB2_Temp.dat" index 14 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 14 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 15 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_16_17.ps"

set title "Upstream collimator, ring chisquare versus B2"
set xlabel "B2 [DAC]"
set ylabel "Chisquare"

plot "ShiftVsB2_Temp.dat" index 16 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 16 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 17 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "ShiftVsB2_Temp_idx_18_19.ps"

set title "Downstream collimator, ring chisquare versus B2"
set xlabel "B2 [DAC]"
set ylabel "Chisquare"

plot "ShiftVsB2_Temp.dat" index 18 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 18 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "ShiftVsB2_Temp.dat" index 19 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints






#########################################################################################
# Position and angle wrt the yoke


set key screen 0.92, screen 0.889
set rmargin 16

#__________________________________

set output "CentAngVsB2_Temp_idx_0_2.ps"

set title "X position of the center of the detector wrt the yoke versus B2"
set xlabel "B2 [DAC]"
set ylabel "X position [cm]"

plot "CentAngVsB2_Temp.dat" index 0 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 0 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 2 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "CentAngVsB2_Temp_idx_1_3.ps"

set title "Y position of the center of the detector wrt the yoke versus B2"
set xlabel "B2 [DAC]"
set ylabel "Y position [cm]"

plot "CentAngVsB2_Temp.dat" index 1 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 1 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 3 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints





set output "CentAngVsB2_Temp_idx_4_6.ps"

set title "X angle of the detector wrt the yoke versus B2"
set xlabel "B2 [DAC]"
set ylabel "X angle [mrad]"

plot "CentAngVsB2_Temp.dat" index 4 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 4 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 6 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints



set output "CentAngVsB2_Temp_idx_5_7.ps"

set title "Y angle of the detector wrt the yoke versus B2"
set xlabel "B2 [DAC]"
set ylabel "Y angle [mrad]"

plot "CentAngVsB2_Temp.dat" index 5 using 1:2 t "Hole -3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:3 t "Hole -2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:4 t "Hole -1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:5 t "Hole 0,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:6 t "Hole 1,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:7 t "Hole 2,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 5 using 1:8 t "Hole 3,0" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:2 t "Hole 0,-3" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:3 t "Hole 0,-2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:4 t "Hole 0,-1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:6 t "Hole 0,1" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:7 t "Hole 0,2" ps 2 lw 2 w linespoints, "CentAngVsB2_Temp.dat" index 7 using 1:8 t "Hole 0,3" ps 2 lw 2 w linespoints