SLAC Magnetic Measurements
Date: 06-23-1995
Time: 08:34:16

Magnet Name (Serial #): 000
Bar Code Number: 001
Project: PEP II High Energy Ring Sextupoles
Test Stand: IR8, S1
Measurement Coil: # 1
Operator: pdr
Run Number: 1
Comment: Check sextupole measurement program

Standardization Currents (A):
   250.0    10.0   200.0    10.0   200.0    10.0   200.0
    10.0

Test Currents (A):
    25.0    50.0    75.0   100.0   125.0   150.0   175.0
   200.0   175.0   125.0    75.0    25.0



            INTEGRATED SEXTUPOLE STRENGTH VS CURRENT

Coil: Radius =  .044882  m, # Turns =  20 
Average: # Rotations/Measurement =  10 , # Measurements =  4 

  Imag    sigImag     SL     sigSL      SL/I   sigSL/I 
  (A)       (A)     (T/m)    (T/m)    (T/m/kA) (T/m/kA)
--------+--------  --------+--------  --------+--------
  25.622    0.001   5.23619  0.00186  204.3652   0.0732
  50.605    0.001  10.20650  0.00256  201.6892   0.0507
  75.612    0.001  15.23926  0.00137  201.5456   0.0184
 100.599    0.001  20.25043  0.00539  201.2993   0.0536
 125.577    0.003  25.26371  0.00466  201.1813   0.0374
 150.572    0.002  30.29137  0.00790  201.1758   0.0526
 175.557    0.003  35.26151  0.01010  200.8551   0.0576
 200.565    0.002  40.16702  0.01310  200.2693   0.0654
 175.559    0.002  35.46939  0.00740  202.0367   0.0422
 125.570    0.003  25.54362  0.00338  203.4218   0.0273
  75.603    0.004  15.47060  0.00379  204.6283   0.0511
  25.611    0.002   4.08111  2.23642  159.3472  87.3212


           SUMMARY OF THE CALCULATIONS AND CONVENTIONS USED


Field Expansion:
The expansion of the radial and azimuthal field in polar
coordinates is
 Br(r,th) = Sum Bn (r/rref)^(n-1) cos(n(th-thspole))
 Bth(r,th) = -Sum Bn (r/rref)^(n-1) sin(n(th-thspole))
Our convention is to set rref = Rcoil.
Thspole is the angle of the first magnetic south pole
with respect to the horizontal, measured ccw by a shaft
encoder in the system.

Coil Voltage:
The coil voltage from each field harmonic is
 Vn(th) = Nturns * velocity * Brn(Rcoil,th) * L
L is the magnet effective length, or BLn is the integrated
field strength of the n'th harmonic.
At the coil radius, the radial field as a function of angle is,
 Brn(Rcoil,th) * L = BLn * cos(n*(th - thspole))
The coil voltage is
 Vn(th) = Nturns * velocity * BLn * cos(n*(th - thspole))
        = Nturns * Rcoil * ang_freq * BLn * cos(n*(th - THspole))
An FFT of the coil voltage gives Vn and PhiVn according to the formula
 Vn(i) = Vn * cos(n*2pi*i/N + PhiVn)

Multipole Field Calculations:
To find the multipole field magnitudes and phases,
the measured voltage harmonics are related to their values
calculated from the field harmonics:
 Nturns * Rcoil * ang_freq * BLn = Vn
 -n * thspole = PhiVn
Or,
 BLn = Vn / (Nturns * Rcoil * ang_freq)
 thspole = -PhiVn / n

Harmonic Strength Ratios:
The main field, denoted by capital N, is the field harmonic
with the largest strength at the coil radius.
The field strength ratio is defined by
 Rn = BLn / BLN
It gives the ratio of each harmonic field strength to the
main field strength at the coil radius.

Calculation Of SL:
The sextupole strength S is defined by the vertical
field on the x-axis (median plane).
It is the quadratic term in the Taylor expansion.
 By(x) = 1/2 S x^2
On the x-axis, By(x) = Bth(r=x,th=0)
From the expression for Bth above,
 By(x) = - B3 (x/rref)^2 sin(-3 thspole)
Take rref = Rcoil, thspole = pi/6, then at x=Rcoil,
 By(Rcoil) = B3 = 1/2 S Rcoil^2
So,
 S = 2 B3 / Rcoil^2
The integrated sextupole strength is
 SL = 2 BL3 / Rcoil^2

Calculation Of The Sextupole Center:
In the sextupole's frame,
 Bx' = S * x' * y'
 By' = 1/2 * S * (x'^2 - y'^2).
In the coil's frame (unprimed frame) the magnetic center is at (x0, y0).
In the coil's frame,
 Bx = S * (x - x0) * (y - y0)
 By = 1/2 * S * [(x - x0)^2 - (y - y0)^2].
The magnetic center can be found in terms of the measured
quadrupole field.  Compare Br and Btheta on the x-axis to
 Bxquad = -S * (x0 * y + y0 * x)
 Byquad = -S * (x0 * x - y0 * y).
Evaluate at x = Rcoil, y = 0 and compare to
 Br(theta = 0) = Bx, Btheta(theta = 0) = By.
In terms of the measured integrated strengths, this gives
 Xcenter = x0 = - (1/(SL * Rcoil)) * BL2 * sin(2 * THspole2)
 Ycenter = y0 = - (1/(SL * Rcoil)) * BL2 * cos(2 * THspole2)