lcsim-cal-calib
diff -u -r1.1 -r1.2
--- sf.plt 20 May 2008 21:34:38 -0000 1.1
+++ sf.plt 27 May 2008 18:24:36 -0000 1.2
@@ -5,29 +5,39 @@
#plot "sf.txt" using 1:2 title 'Theta vs EM1 SF'
# break fitting region into barrel and endcap
-sf(x)=sf0+a*sin(2*pi*(x-90)/360-pi/2)
-sf_2(x)=sf0_2+a_2*sin(2*pi*(x-90)/360)
+
# EM1 Sampling Fractions
+em1_barrel(x)= em1b_0 + em1b_1*sin(2*pi*(x-90)/360-pi/2)
+em1_endcap(x)= em1e_0 + em1e_1*sin(2*pi*(x-90)/360)
+
-fit [80.:145.] sf(x) "sf.txt" using 1:2 via sf0, a
-fit [145.:180.] sf_2(x) "sf.txt" using 1:2 via sf0_2, a_2
+fit [80.:145.] em1_barrel(x) "sf.txt" using 1:2 via em1b_0, em1b_1
+fit [145.:180.] em1_endcap(x) "sf.txt" using 1:2 via em1e_0, em1e_1
#only plot function within range of validity
set ylabel "EM1 Sampling Fraction"
-plot "sf.txt" using 1:2 title 'Theta vs EM1 SF', (x>145. ? 1/0 : sf(x)), (x<145. ? 1/0 : sf_2(x))
+plot "sf.txt" using 1:2 title 'Theta vs EM1 SF', (x>145. ? 1/0 : em1_barrel(x)), (x<145. ? 1/0 : em1_endcap(x))
+
+
# EM2 Sampling Fractions
+em2_barrel(x)= em2b_0 + em2b_1*sin(2*pi*(x-90)/360-pi/2)
+em2_endcap(x)= em2e_0 + em2e_1*sin(2*pi*(x-90)/360)
-fit [80.:145.] sf(x) "sf.txt" using 1:3 via sf0, a
-fit [145.:180.] sf_2(x) "sf.txt" using 1:3 via sf0_2, a_2
+fit [80.:145.] em2_barrel(x) "sf.txt" using 1:3 via em2b_0, em2b_1
+fit [145.:180.] em2_endcap(x) "sf.txt" using 1:3 via em2e_0, em2e_1
set ylabel "EM2 Sampling Fraction"
-plot "sf.txt" using 1:3 title 'Theta vs EM2 SF', (x>145. ? 1/0 : sf(x)), (x<145. ? 1/0 : sf_2(x))
+plot "sf.txt" using 1:3 title 'Theta vs EM2 SF', (x>145. ? 1/0 : em2_barrel(x)), (x<145. ? 1/0 : em2_endcap(x))
+
+
# HCal Sampling Fractions
+had_barrel(x)= hadb_0 + hadb_1*sin(2*pi*(x-90)/360-pi/2)
+had_endcap(x)= hade_0 + hade_1*sin(2*pi*(x-90)/360)
-fit [80.:145.] sf(x) "sf.txt" using 1:4 via sf0, a
-fit [145.:180.] sf_2(x) "sf.txt" using 1:4 via sf0_2, a_2
+fit [80.:145.] had_barrel(x) "sf.txt" using 1:4 via hadb_0, hadb_1
+fit [145.:180.] had_endcap(x) "sf.txt" using 1:4 via hade_0, hade_1
set ylabel "HCal Sampling Fraction"
-plot "sf.txt" using 1:4 title 'Theta vs HCal SF', (x>145. ? 1/0 : sf(x)), (x<145. ? 1/0 : sf_2(x))
+plot "sf.txt" using 1:4 title 'Theta vs HCal SF', (x>145. ? 1/0 : had_barrel(x)), (x<145. ? 1/0 : had_endcap(x))