gmuwork54

columns:
0  Iz(mA)
1  sigma_Iz(mA)
2  Bfield (Gauss)
3  sigma_B(Gauss)
4  f_low (kHz)
5 sig_f_low (kHz)
6 gmu_85(MeV/G)
7 sigma_gmu_85 (MeV)
8 Rel error
9 1/sig^2
10 xi/sig^2
11 Results

!io=24.72124161;sigio=0.0694;R=0.283;z=0.175

!set col 2 of gmuwork to 1.256637e-2*1e-3*180*abs(col 0 - 24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 3 of gmuwork to c*sqrt((R*col 1)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)


!set col 6 of gmuwork to 4.1375e-21*col 4*1e3/col 2


!h=4.1375e-21;
set col 7 of gmuwork=h*1e3/col 2*sqrt((col 5)**2+(col 4*col 3/col 2)**2)

!set col 8 of gmuwork=col 7/col 6*100
!set col 9 of gmuwork to 1/((col 7)**2)
!set col 10 of gmuwork to col 6*col 9

!do n = 1 to 15;
set col 11 row 1 of gmuwork=col 11 row 1 of gmuwork + col 9 row n;
end
!set col 11 row 6 of gmuwork=sqrt(1/col 11 row 1)

!do n = 1 to 15;
set col 11 row 2 of gmuwork=col 11 row 2 of gmuwork + col 10 row n;
end
!set col 11 row 7 of gmuwork=col 11 row 2 / col 11 row 1


!do n = 1 to 15;
set col 11 row 3 of gmuwork=col 11 row 3 of gmuwork + col 9 row n;
end
!set col 11 row 10 of gmuwork=sqrt(1/col 11 row 3)

!do n = 1 to 15;
set col 11 row 4 of gmuwork=col 11 row 4 of gmuwork + col 10 row n;
end
!set col 11 row 11 of gmuwork=col 11 row 4 / col 11 row 3


!set col 14 of gmuwork = col 7 - 2.945865e-15
!do n=1 to 15;
set col 16 row 1 of gmuwork=col 16 row 1+(col 14 row n)**2;
end

*** gmuwork is saved to gmuwork3
before following

---------------------------------------------------------------
columns:
0  Iz(mA)
1  sigma_Iz(mA)
2  Bfield (Gauss)
3  sigma_B(Gauss)
4  f_low (kHz)
5 sig_f_low (kHz)
6 gmu_85(MeV/G)
7 sigma_gmu_85 (MeV)
8 Rel error
9 1/sig^2
10 xi/sig^2
11 Results

!io=24.235;sigio=0.02;R=0.283;z=0.175

!set col 2 of gmuwork to 1.256637e-2*1e-3*180*abs(col 0 - 24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.02;io=-24.24;Rs=0.01;zs=0.005;
set col 3 of gmuwork to c*sqrt((R*0)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)


!set col 6 of gmuwork to 4.1375e-21*col 4*1e3/col 2


!h=4.1375e-21;
set col 7 of gmuwork=h*1e3/col 2*sqrt((col 4*col 3/col 2)**2)

!set col 8 of gmuwork=col 7/col 6*100
!set col 9 of gmuwork to 1/((col 7)**2)
!set col 10 of gmuwork to col 6*col 9

!do n = 1 to 15;
set col 11 row 1 of gmuwork=col 11 row 1 of gmuwork + col 9 row n;
end
!set col 11 row 6 of gmuwork=sqrt(1/col 11 row 1)

!do n = 1 to 15;
set col 11 row 2 of gmuwork=col 11 row 2 of gmuwork + col 10 row n;
end
!set col 11 row 7 of gmuwork=col 11 row 2 / col 11 row 1


!do n = 1 to 15;
set col 11 row 3 of gmuwork=col 11 row 3 of gmuwork + col 9 row n;
end
!set col 11 row 10 of gmuwork=sqrt(1/col 11 row 3)

!do n = 1 to 15;
set col 11 row 4 of gmuwork=col 11 row 4 of gmuwork + col 10 row n;
end
!set col 11 row 11 of gmuwork=col 11 row 4 / col 11 row 3


!set col 14 of gmuwork = col 7 - 2.945865e-15
!do n=1 to 15;
set col 16 row 1 of gmuwork=col 16 row 1+(col 14 row n)**2;
end


!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 3 of gmuwork to c*sqrt((R*col 1)**2+(R*ios)**2+(abs(io+col
0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 17 of gmuwork to c*sqrt((abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2)

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 18 of gmuwork to c*sqrt((abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 19 of gmuwork to c*sqrt((R*ios)**2)

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 20 of gmuwork to c*sqrt((R*col 1)**2)



4/30/90

!set col 1 of gmufin to col 6 of gmuwork2
!set col 2 of gmufin to col 7 of gmuwork2

!set col 3 of gmufin = col 1 + col 2
!set col 4 of gmufin = col 1 - col 2

!set col 5 of gmufin to col 6 of gmuwork
!set col 6 of gmufin to col 7 of gmuwork

!set col 7 of gmufin = col 5 + col 6
!set col 8 of gmufin = col 5 - col 6

plot graph ggmu2c noheader box notitle on file 'gmu.tek'


 # plot graph foo bottomkey on file 'foo.tek'


5/4/90

g85:(from gmuwork2)#########################################################
!io=24.235;sigio=0.02;R=0.283;z=0.175

!set col 2 of g85 to 1.256637e-2*1e-3*180*abs(col 0 - 24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5


!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.02;io=-24.24;Rs=0.0068;zs=0;
set col 3 of g85 to c*sqrt((R*0)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!set col 6 of g85 to 4.1375e-21*col 4*1e3/col 2


!h=4.1375e-21;
set col 7 of g85=h*1e3/col 2*sqrt((col 4*col 3/col 2)**2)

!set col 8 of g85=col 6-col 7
!set col 9 of g85=col 6+col 7

!set col 10 of g85=col 7/col 6*100

g87:(from gmuwork3)########################################################
!io=24.235;sigio=0.02;R=0.283;z=0.175

!set col 2 of g87 to 1.256637e-2*1e-3*180*abs(col 0 - 24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5


!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.02;io=-24.24;Rs=0.0068;zs=0;
set col 3 of g87 to c*sqrt((R*0)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!set col 6 of g87 to 4.1375e-21*col 4*1e3/col 2


!h=4.1375e-21;
set col 7 of g87=h*1e3/col 2*sqrt((col 4*col 3/col 2)**2)


!set col 8 of g87=col 6-col 7
!set col 9 of g87=col 6+col 7

!set col 10 of g87=col 7/col 6*100

gmufin:####################################################################
!set col 1 of gmufin to col 6 of g85
!set col 2 of gmufin to col 7 of g85

!set col 3 of gmufin = col 1 + col 2
!set col 4 of gmufin = col 1 - col 2

!set col 5 of gmufin to col 6 of g87
!set col 6 of gmufin to col 7 of g87

!set col 7 of gmufin = col 5 + col 6
!set col 8 of gmufin = col 5 - col 6


 # plot graph foo bottomkey on file 'foo.tek'
g85s:********************(from g85) Statistical errors
!set col 5 of g85s to col 2 of gmuratio

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.0;io=-24.24;Rs=0.00;zs=0.000;
set col 3 of g85s to c*sqrt((R*col 1)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!h=4.1375e-21;
set col 7 of g85s=h*1e3/col 2*sqrt((col 5)**2+(col 4*col 3/col 2)**2)

!set col 8 of g85s=col 6-col 7
!set col 9 of g85s=col 6+col 7
!set col 10 of g85s=col 7/col 6*100

;; Compute 1/sigma^2
!set col 11 row 1 of g85s = 0;
do n=1 to 14;
set col 11 row 1 of g85s = col 11 row 1 + 1/(col 7 row n)**2;
end

;; Compute x/sigma^2
!set col 11 row 2 of g85s = 0;
do n=1 to 14;
set col 11 row 2 of g85s = col 11 row 2 + (col 6 row n)/(col 7 row n)**2;
end

;; Compute weighted mean and deviation
!set col 11 row 6 of g85s = sqrt(1/col 11 row 1)
!set col 11 row 7 of g85s = col 11 row 2/col 11 row 1


;;Compute average
!set col 11 row 14 of g85s = 0;
do n=1 to 14;
set col 11 row 14 of g85s = col 11 row 14 + col 6 row n;
end
!set col 11 row 13 of g85s = col 11 row 14/14

g87s:********************(from g87) Statistical errors
!set col 5 of g87s to col 4 of gmuratio

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.0;io=-24.24;Rs=0.00;zs=0.000;
set col 3 of g87s to c*sqrt((R*col 1)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)

!h=4.1375e-21;
set col 7 of g87s=h*1e3/col 2*sqrt((col 5)**2+(col 4*col 3/col 2)**2)

!set col 8 of g87s=col 6-col 7
!set col 9 of g87s=col 6+col 7
!set col 10 of g87s=col 7/col 6*100

!set col 11 row 1 of g87s = 0;
do n=1 to 14;
set col 11 row 1 of g87s = col 11 row 1 + 1/(col 7 row n)**2;
end

!set col 11 row 2 of g87s = 0;
do n=1 to 14;
set col 11 row 2 of g87s = col 11 row 2 + (col 6 row n)/(col 7 row n)**2;
end

;; Compute weighted mean and deviation
!set col 11 row 6 of g87s = sqrt(1/col 11 row 1)
!set col 11 row 7 of g87s = col 11 row 2/col 11 row 1


plot graph ggmu2d noheader box bottomkey notitle on file '/usr/tmp/gmu.tek'




FINAL!!!

!set col 2 of gmuwork to 1.256637e-2*1e-3*180*abs(col 0 - 24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.01;zs=0.005;
set col 3 of gmuwork to c*sqrt((R*col 1)**2+(R*ios)**2+(abs(io+col 0)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+col 0)*R*z*zs/(z**2+R**2))**2)


Z field
!set col 0 row 1 of final to 1.256637e-2*1e-3*180*abs(24.24)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=0.2;io=-24.24;Rs=0.0068;zs=0.000;
set col 1 row 1 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+0)*R*z*zs/(z**2+R**2))**2)


!set col 0 row 1 of temp to 1.256637e-2*1e-3*180*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!set col 0 row 2 of temp to 1.256637e-2*1e-3*75*(0.226)**2/((0.143)**2+(0.226)**2)**1.5

!set col 0 row 3 of temp to 1.256637e-2*1e-3*50*(0.173)**2/((0.089)**2+(0.173)**2)**1.5


Y field
!set col 0 row 2 of final to 1.256637e-2*1e-3*75*abs(17.6)*(0.226)**2/((0.143)**2+(0.226)**2)**1.5

!R=0.226;z=0.143;mu=1.256637e-2;c=mu*1e-3*75*R/(z**2+R**2)**1.5;
ios=0.2;io=-17.6;Rs=0.005;zs=0.000;
set col 1 row 2 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io)*R*z*zs/(z**2+R**2))**2)

X field
!set col 0 row 3 of final to 1.256637e-2*1e-3*50*abs(124.0)*(0.173)**2/((0.089)**2+(0.173)**2)**1.5

!R=0.173;z=0.089;mu=1.256637e-2;c=mu*1e-3*50*R/(z**2+R**2)**1.5;
ios=0.9;io=-124.0;Rs=0.005;zs=0.000;
set col 1 row 3 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+0)*R*z*zs/(z**2+R**2))**2)



2nd Part
Z field
!set col 0 row 5 of final to 1.256637e-2*1e-3*180*abs(24.0)*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!R=0.283;z=0.175;mu=1.256637e-2;c=mu*1e-3*180*R/(z**2+R**2)**1.5;
ios=2.0;io=-24.0;Rs=0.0068;zs=0.000;
set col 1 row 5 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+0)*R*z*zs/(z**2+R**2))**2)

Y field
!set col 0 row 6 of final to 1.256637e-2*1e-3*75*abs(9.5)*(0.226)**2/((0.143)**2+(0.226)**2)**1.5

!R=0.226;z=0.143;mu=1.256637e-2;c=mu*1e-3*75*R/(z**2+R**2)**1.5;
ios=1.0;io=-9.5;Rs=0.005;zs=0.000;
set col 1 row 6 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io)*R*z*zs/(z**2+R**2))**2)

X field
!set col 0 row 7 of final to 1.256637e-2*1e-3*50*abs(109.4)*(0.173)**2/((0.089)**2+(0.173)**2)**1.5

!R=0.173;z=0.089;mu=1.256637e-2;c=mu*1e-3*50*R/(z**2+R**2)**1.5;
ios=0.6;io=-109.4;Rs=0.005;zs=0.000;
set col 1 row 7 of final to c*sqrt((R*0)**2+(R*ios)**2+(abs(io)*(2*z**2+R**2)*Rs/(z**2+R**2))**2+(3*abs(io+0)*R*z*zs/(z**2+R**2))**2)


!set col 2 of final = col 1/col 0*100

columns:
0  Iz(mA)
1  sigma_Iz(mA)
2  Bfield (Gauss)
3  sigma_B(Gauss)
4  f_low (kHz)
5 sig_f_low (kHz)
6 gmu_85(MeV/G)
7 sigma_gmu_85 (MeV)
8 Rel error
9 1/sig^2
10 xi/sig^2
11 Results


fit function "c0 + c1/sqrt((x-x0)**2 + c2)" to curve 1 of gop4x

                           PARAMETER TABLE

0    1 PARAM.  2 FINAL      3 STD       4 T VALUE  5 SIG.    6 DEP.
                 VALUES       ERRORS                 LEV.      VALUES
---------------------------------------------------------------------
1    C0          -2.280897    0.643314     -3.546  0.012130   0.99480
2    C1         370.821981   52.581705      7.052  0.000407   0.99820
3    X0         109.420010    0.671055    163.100  0.000100   0.05426
4    C2        1458.887974  239.853053      6.082  0.000898   0.99120

fit function "c0 + c1/sqrt((x-x0)**2 + c2**2)" to curve 1 of gop4y
Result of FIT is
    4.795285+26.30486/SQRT((X-9.546274)**2+85.704523)

FIT converged.  This does NOT necessarily imply a good fit.
Final SUM OF SQUARES is 0.069879.
2 iterations used.
                          PARAMETER TABLE

0    1 PARAM.  2 FINAL    3 STD      4 T VALUE  5 SIG.    6 DEP.
                 VALUES     ERRORS                LEV.      VALUES
------------------------------------------------------------------
1    C0         4.795285   1.169291      4.101  0.026230    0.9976
2    C1        26.304860  23.094678      1.139  0.337400    0.9993
3    X0         9.546274   0.968659      9.855  0.002221    0.4540
4    C2        85.704523  84.801840      1.011  0.386500    0.9969

plot graph gop4x on file '/usr/tmp/a'
plot graph gop4y on file '/usr/tmp/b'


!set col 2 of optical1 = col 1 - 2
!set col 3 of optical1 = col 1 + 2


!set col 5 of optical1 = col 4 - 3
!set col 6 of optical1 = col 4 +3

plot graph gop1f1 on file '/usr/tmp/x.tek'
plot graph gmuplot on file '/usr/tmp/z.tek'
plot graph ggmu on file '/usr/tmp/zabs.tek'


!set col 0 row 1 of temp to 1.256637e-2*1e-3*180*(0.283)**2/((0.175)**2+(0.283)**2)**1.5

!set col 0 row 2 of temp to 1.256637e-2*1e-3*75*(0.226)**2/((0.143)**2+(0.226)**2)**1.5

!set col 0 row 3 of temp to 1.256637e-2*1e-3*50*(0.173)**2/((0.089)**2+(0.173)**2)**1.5


gmuratio:

!set col 0 of gmuratio to col 0 of g85
!set col 1 of gmuratio to col 4 of g85
!set col 2 of gmuratio to col 5 of g85
!set col 3 of gmuratio to col 4 of g87
!set col 4 of gmuratio to col 5 of g87

!set col 5 of gmuratio = col 3/col 1
!set col 6 of gmuratio = sqrt((col 2/col 3)**2+((col 4*col 1/(col 3)**2))**2)
!set col 7 of gmuratio = col 6/col 5*100
!set col 8 of gmuratio = col 5-col 6
!set col 9 of gmuratio = col 5+col 6



col 8 of gmuratio
col 9 of gmuratio
