exe hallc#setup
 
 fortran/file 66 sd2000_both_sys.eps
 meta 66 -113
 
 label 1 7 ' ' ' ' ' ' ' ' ' ' ' ' ' '
 set ndvx 6.11
 set ygti 1.0

*  Order is: Mott,A Moller, B Moller, C Moller, A Compton
 ve/create pol(5) r 72.10 76.76 69.84 73.30 73.33
 ve/create epol(5) r 0.28 1.04 0.64 0.62 1.42
 
 ve/create polb(5) r 72.27 75.40 70.01 73.37 73.85
 ve/create epolb(5) r 0.23 0.33 1.80 0.42 2.15
 
 
*  these are relative - need to be multiplied be abs(pol)
 ve/create esys(5) r 1.06 2.98 1.61 0.94 1.0
 sigma esys = esys/100.0*pol
 
  ve/create esysb(5) r 1.06 2.98 1.61 0.94 1.0
  sigma esysb = esysb/100.0*polb
 
 
 ve/cr hall(5) r 1 2 3 4 5 
 ve/cr hallb(5) r 1.1 2.1 3.1 4.1 5.1 
 
 norm=pol(1)
 
 sigma pol=pol/[norm]
 sigma epol=epol/[norm]
 sigma esys=esys/[norm]
 
 
 sigma polb=polb/[norm]
 sigma epolb=epolb/[norm]
 sigma esysb=esysb/[norm]
 
 sigma etot = sqrt(epol**2+esys**2)
 sigma etotb = sqrt(epolb**2+esysb**2)
 
 ve/create ex(5) r
 
 null 0 6. 0.9 1.1
 
 do i=1,5
   
   if [i]=5 then
     set hcol 6
     set pmci 6
   else
     set hcol [i] 
     set pmci [i]
   endif
   hplot/errors hall([i]:[i]) pol([i]:[i]) ex([i]:[i]) epol([i]:[i]) 1 20 0.2 e1
   hplot/errors hallb([i]:[i]) polb([i]:[i]) ex([i]:[i]) epolb([i]:[i]) 1 24 0.2 e1
   
   hplot/errors hall([i]:[i]) pol([i]:[i]) ex([i]:[i]) etot([i]:[i]) 1 20 0.2
   hplot/errors hallb([i]:[i]) polb([i]:[i]) ex([i]:[i]) etotb([i]:[i]) 1 24 0.2
   
 enddo

 
 atit ' ' 'P?polarimeter!/P?Mott!'
 set ltyp 2
 line 0 1 6.5 1
 
  set tang 20.0
 
 ypos =0.91
 
 
 xpos = 0.6
 set txci 1
 itx [xpos] [ypos] 'Mott'
 
 
 xpos = [xpos]+1
 set txci 2
 itx [xpos] [ypos] 'Moller A'
 
 xpos = [xpos]+1
 set txci 3
 itx [xpos] [ypos] 'Moller B'
 
 xpos = [xpos]+1
 set txci 4
 itx [xpos] [ypos] 'Moller C'
 
 xpos = [xpos]+1
 set txci 6
 itx [xpos] [ypos] 'Compton A'
 
 close 66