source: MML/trunk/machine/SOLEIL/common/magnetcoefficients.m

function [C, Leff, MagnetType, A] = magnetcoefficients(MagnetCoreType, Amps, InputType)
%MAGNETCOEFFICIENTS - Retrieves coefficient for conversion between Physics and Hardware units
%[C, Leff, MagnetType, A] = magnetcoefficients(MagnetCoreType)
%
% INPUTS
% 1. MagnetCoreType - Family name or type of magnet
%
% OUTPUTS
% 1. C vector coefficients for the polynomial expansion of the magnet field
%    based on magnet measurements
% 2. Leff - Effective length ie, which is used in AT
% 3. MagnetType
% 4. A - vector coefficients for the polynomial expansion of the curviline 
%        integral of the magnet field based on magnet measurements
%
% C and A are vector coefficients for the polynomial expansion of the magnet field
% based on magnet measurements.
%
% The amp2k and k2amp functions convert between the two types of units.
%   amp2k returns BLeff, B'Leff, or B"Leff scaled by Brho if A-coefficients are used.
%   amp2k returns B    , B'    , or B"     scaled by Brho if C-coefficients are used.
%
% The A coefficients are direct from magnet measurements with a DC term:
%   a8*I^8+a7*I^7+a6*I^6+a5*I^5+a4*I^4+a3*I^3+a2*I^2+a1*I+a0 = B*Leff or B'*Leff or B"*Leff
%   A = [a8 a7 a6 a5 a4 a3 a2 a1 a0]
%
% C coefficients have been scaled to field (AT units, except correctors) and includes a DC term:
%   c8 * I^8+ c7 * I^7+ c6 * I^6 + c5 * I^5 + c4 * I^4 + c3 * I^3 + c2 * I^2 + c1*I + c0 = B or B' or B"
%   C = A/Leff
%
% For dipole:      k = B / Brho      (for AT: KickAngle = BLeff / Brho)
% For quadrupole:  k = B'/ Brho
% For sextupole:   k = B"/ Brho / 2  (to be compatible with AT)
%                  (all coefficients all divided by 2 for sextupoles)
%
% MagnetCoreType is the magnet measurements name for the magnet core (string, string matrix, or cell)
%   For SOLEIL:   BEND
%                 Q1 - Q10 S1 - S10,
%                 QT, HCOR, VCOR, FHCOR, FVCOR
%
% Leff is the effective length of the magnet
%
% See Also amp2k, k2amp

%
% Written by M. Yoon 4/8/03
% Adapted By Laurent S. Nadolski354.09672
%
% Partie Anneau modifiï¿œe par P. Brunelle et A. Nadji le 31/03/06
%
% Add a switch on accelerator 

% NOTE: Make sure the sign on the 'C' coefficients is reversed where positive current generates negative K-values
% Or use Tango K value set to -1
%
% NOTE: Q11 and Q12 are not well configured yet, just quick fix


if nargin < 1
    error('MagnetCoreType input required');
end

if nargin < 2
    Amps = 230;  % not sure!!!
end

if nargin < 3
    InputType = 'Amps';
end



% For a string matrix
if iscell(MagnetCoreType)
    for i = 1:size(MagnetCoreType,1)
        for j = 1:size(MagnetCoreType,2)
            [C{i,j}, Leff{i,j}, MagnetType{i,j}, A{i,j}] = magnetcoefficients(MagnetCoreType{i});
        end
    end
    return
end

% For a string matrix
if size(MagnetCoreType,1) > 1
    C=[]; Leff=[]; MagnetType=[]; A=[];
    for i = 1:size(MagnetCoreType,1)
        [C1, Leff1, MagnetType1, A1] = magnetcoefficients(MagnetCoreType(i,:));
        C(i,:) = C1;
        Leff(i,:) = Leff1;
        MagnetType = strvcat(MagnetType, MagnetType1);
        A(i,:) = A1;
    end
    return
end

%% get accelerator name
AcceleratorName = getfamilydata('SubMachine');

switch AcceleratorName
    case 'LT1'
        %%%%
        switch upper(deblank(MagnetCoreType))

            case 'BEND'    
                Leff = 0.30; % 300 mm
                % B = 1e-4 * (0.0004 Iￜ + 16.334 I + 1.7202)
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  4.8861e-4;
                a0 =  1.19e-4;

                A = [a8 a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'BEND';

            case {'QP'}   % 150 mm quadrupole
                % Find the current from the given polynomial for B'Leff
                Leff=0.150; % 162 mm;
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
%                 a4 =  1.49e-6;
%                 a3 =  2.59e-5;
%                 a2 =  1.93e-4;
%                 a1 =  4.98e-2;
%                 a0 =  0.0;
                a4 =  -1.49e-6;
                a3 =  2.59e-5;
                a2 =  -1.93e-4;
                a1 =  4.98e-2;
                a0 =  8.13e-4;              
                
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'QUAD';

            case {'CH','CV'}    % 16 cm horizontal corrector
                % Magnet Spec: Theta = 0.8e-3 radians @ 2.75 GeV and 10 amps
                % Theta = BLeff / Brho    [radians]
                % Therefore,
                %       Theta = ((BLeff/Amp)/ Brho) * I
                %       BLeff/Amp = 0.8e-3 * getbrho(2.75) / 10
                %       B*Leff = a0 * I   => a0 = 0.8e-3 * getbrho(2.75) / 10
                %
                % The C coefficients are w.r.t B
                %       B = c0 + c1*I = (0 + a0*I)/Leff
                % However, AT uses Theta in radians so the A coefficients
                % must be used for correctors with the middle layer with
                % the addition of the DC term

                % Find the current from the given polynomial for BLeff and B
                % NOTE: AT used BLeff (A) for correctors
                MagnetType = 'COR';
                
                Leff = 1e-6; % 0.1577 m
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  4.49e-4;
                a0 =  0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
            otherwise
                error(sprintf('MagnetCoreType %s is not unknown', MagnetCoreType));
                %k = 0;
                %MagnetType = '';
                %return
        end

        % compute B-field = int(Bdl)/Leff
        C = A/ Leff;

        MagnetType = upper(MagnetType);
    
    case 'StorageRing'
        %%%%
        switch upper(deblank(MagnetCoreType))

            case 'BEND'   
                % Moyenne des longueurs magnï¿œtiques mesurï¿œes = 1055.548mm
                % Décalage en champ entre le dipￜle de référence et les
                % dipï¿œles de l'Anneau = DB/B= +1.8e-03.
                % On part de l'ￜtalonnage B(I) effectuￜ sur le dipￜle de
                % rï¿œfï¿œrence dans la zone de courant 516 - 558 A 
                % les coefficients du fit doivent ï¿œtre affectï¿œs du facteur
                % (1-1.8e-3) pour passer du dipￜle de rￜfￜrence ￜ l'Anneau
                % et du facteur Leff pour passer ￜ l'intￜgrale de champ.
                % 
                
                % B=1.7063474 T correspond ￜ 2.75 GeV 
                % ?  longueur magnétique du model : Leff = 1.052433;
%                 Leff=1.055548;
%                 a7= 0.0;
%                 a6=-0.0;
%                 a5= 0.0;
%                 a4=-0.0;
%                 a3= 0.0;
%                 a2=-9.7816E-6*(1-1.8e-3)*Leff;
%                 a1= 1.26066E-02*(1-1.8E-3)*Leff;
%                 a0= -2.24944*(1-1.8E-3)*Leff;
%                 A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                Leff=1.052433;
                a7= 0.0;
                a6=-0.0;
                a5= 0.0;
                a4=-0.0;
                a3= 0.0;
                a2=-9.7816E-6*(1-1.8e-3)*Leff*(1.055548/1.052433);
                a1= 1.26066E-02*(1-1.8E-3)*Leff*(1.055548/1.052433);
                a0= -2.24944*(1-1.8E-3)*Leff*(1.055548/1.052433);
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                

                MagnetType = 'BEND';

            case {'Q1','Q6', 'Q11'}  
                % Familles Q1 et Q6 l= 320 mm 
                % Etalonnage GL(I) sur 90 - 150 A quadrupï¿œle court 
                % le courant remontￜ est nￜgatif car Q1 et Q6 dￜfocalisants
                % il faut donc un k < 0. Les coefficients du fit a0, a2,
                % a4,...sont multipliï¿œs par -1.
                
                % Correction des coefficients des QC de + 3 10-3 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);
%                   bob=1-8.0e-3-2.3e-3+3e-3;
                % Find the current from the given polynomial for B'Leff
                Leff=0.320;

        % G. Portmann
%                 % Just a way to get the correct polynomial for different k-values
%                 if strcmpi(InputType, 'K')
%                     if Amps < 0.18928717429288
%                         Amps = 10;
%                     elseif Amps < 0.75086308092911
%                         Amps = 50;
%                     elseif Amps < 1.31229873800730
%                         Amps = 100;
%                     elseif Amps < 1.68337408687106
%                         Amps = 150;
%                     elseif Amps < 2.04021595285003
%                         Amps = 200;
%                     else
%                         Amps = 230;
%                     end
%                 end
        % G. Portmann
                
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  0.0;
                a2=  -8.6E-7*(-1)*(1.003)*bob;
                a1=  2.7664E-2*(1.003)*bob;
                a0=  -3.3E-3*(-1)*(1.003)*bob;
%                 a7=  0.0;
%                 a6=  0.0;
%                 a5=  0.0;
%                 a4=  0.0;
%                 a3=  0.0;
%                 a2=  -8.6E-7*(-1)*bob;
%                 a1=  2.7664E-2*bob;
%                 a0=  -3.3E-3*(-1)*bob;
       % G. Portmann                
%                 if Amps < 20
%                     a2=  0.0;
%                     a1=  0.027473;
%                     a0=  0.0;
%                 elseif Amps < 80
%                     a2=  0.0;
%                     a1=  0.027473;
%                     a0=  0.006270;
%                  elseif Amps < 140
%                     a2= -5.6000e-7;
%                     a1=  2.7598e-2;
%                     a0=  2.1000e-4;
%                  elseif Amps < 180
%                     a3= -3.25400e-8;
%                     a2=  1.05300e-5;
%                     a1=  2.63758e-2;
%                     a0=  4.30800e-2;
%                  elseif Amps < 220
%                     a4= -1.14374e-8;
%                     a3=  8.69384e-6;
%                     a2= -2.49129e-3;
%                     a1=  3.45626e-1;
%                     a0= -1.52486e+1;
%                 else
%                     a5= -3.417777770e-8;
%                     a4=  3.992090910e-5;
%                     a3= -1.864441566e-2;
%                     a2=  4.351937350e+0;
%                     a1= -5.076535920e+2;
%                     a0=  2.367859440e+4;
%                 end
         % G. Portmann
         
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'quad';
                
             case {'Q8','Q9'}  
                % Familles Q8 et Q9 l= 320 mm 
                % Etalonnage GL(I) sur 160 - 200 A quadrupï¿œle court 
                % le courant remontￜ est nￜgatif car Q8 et Q9 dￜfocalisants
                % il faut donc un k < 0. Les coefficients du fit a0, a2,
                % a4,...sont multipliï¿œs par -1.
                
                % Correction des coefficients des QC de + 3 10-3 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);

                % Find the current from the given polynomial for B'Leff
                Leff=0.320;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  -8.843E-8*(1.003)*bob;
                a2=  3.6389E-5*(-1)*(1.003)*bob;
                a1=  2.2448E-2*(1.003)*bob;
                a0=  2.382E-1*(-1)*(1.003)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'quad';  
                
            case {'Q3'}  
                % Famille Q3 l= 320 mm 
                % Etalonnage GL(I) sur 50 - 100 A quadrupï¿œle court 
                % le courant remontￜ est nￜgatif car Q3 est dￜfocalisant
                % il faut donc un k < 0. Les coefficients du fit a0, a2,
                % a4,...sont multipliï¿œs par -1.
                
                %Correction des coefficients des QC de + 3 10-3 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);

                % Find the current from the given polynomial for B'Leff
                Leff=0.320;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  0.;
                a2=  1.4E-7*(-1)*(1.003)*bob;
                a1=  2.7471E-2*(1.003)*bob;
                a0=  5.83E-3*(-1)*(1.003)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'quad'; 
                
             case {'Q4'}  
                % Famille Q4 l= 320 mm 
                % Etalonnage GL(I) sur 120 - 170 A quadrupï¿œle court 
                % le courant remontￜ est nￜgatif car Q4 est dￜfocalisant
                % il faut donc un k < 0. Les coefficients du fit a0, a2,
                % a4,...sont multipliï¿œs par -1.
                
                %Correction des coefficients des QC de + 3 10-3 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);
                
                % Find the current from the given polynomial for B'Leff
                Leff=0.320;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  -5.2680E-8*(1.003)*bob;
                a2=  1.9620E-5*(-1)*(1.003)*bob;
                a1=  2.5016E-2*(1.003)*bob;
                a0=  1.1046E-1*(-1)*(1.003)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'quad'; 
                
            case {'Q5','Q10'}   % 320 mm quadrupole
                 % Familles Q5 et Q10 l= 320 mm 
                % Etalonnage GL(I) sur 180 - 230 A quadrupï¿œle court 
                % le courant remontￜ est nￜgatif car Q5 et Q10 sont
                % focalisants
                
                %Correction des coefficients des QC de + 3 10-3 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);
                
                % Find the current from the given polynomial for B'Leff
                Leff=0.320;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  -8.0497E-09*(1.003)*bob;
                a3= 6.01284E-06*(1.003)*bob;
                a2=  -1.696898E-03*(1.003)*bob;
                a1=  2.41175E-01*(1.003)*bob;
                a0=  -1.01064E+01*(1.003)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];

                MagnetType = 'quad';    
                    
            case {'Q2', 'Q12'}   % l= 460 mm 
                % quadrupï¿œle focalisant
                % Etalonnage GL(I) sur 140 - 190 A quadrupï¿œle long
                
                
                %Correction des coefficients des QL de + 1.55 10-2 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);

                % Find the current from the given polynomial for B'Leff
                Leff=0.460;
                a7=  0.0;
                a6=  0.0;
                a5= -0.0;
                a4=  0.0;
                a3= -2.7609E-7*(1.0155)*bob;
                a2=  1.17098E-4*(1.0155)*bob;
                a1=  2.7718E-2*(1.0155)*bob;
                a0=  8.2470E-1*(1.0155)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'quad';

            case {'Q7'}   % l= 460 mm 
                % quadrupï¿œle focalisant
                % Etalonnage GL(I) sur 190 - 230 A quadrupï¿œle long
                
                %Correction des coefficients des QL de + 1.55 10-2 (manque
                % capteur BMS)
                
                %correction offset capteur BMS -2.310-3 P. Brunelle 30/05/06 
                bob=0.9977*(1-8e-3);

                % Find the current from the given polynomial for B'Leff
                Leff=0.460;
                a7=  0.0;
                a6=  0.0;
                a5= 1.50427350E-9*(1.0155)*bob;
                a4=  -1.52722610E-6*(1.0155)*bob;
                a3= 6.16874120E-4*(1.0155)*bob;
                a2=  -1.24044936E-1*(1.0155)*bob;
                a1=  1.24707096E+01*(1.0155)*bob;
                a0=  -4.96304380E+02*(1.0155)*bob;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'quad';

                % Sextupï¿œles : on multiplie les coefficients par 2 car ils
                % sont exprimï¿œs en B"L et non B"L/2
                
            case {'S1','S10', 'S11', 'S12'}    
                % l= 160 mm focalisants
                % Etalonnage HL(I) sur 40 - 160 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  0.0;
                a2=  -3.773E-6;
                a1=  1.5476E-1;
                a0=  2.36991E-1;
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT';
                
                
            case {'S3','S9'}    
                % l= 160 mm dï¿œfocalisants
                % Etalonnage HL(I) sur 80 - 250 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  -2.6735E-8;
                a2=  5.8793E-6*(-1);
                a1=  1.5364E-1;
                a0=  2.7867E-1*(-1);
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT';
                
             case {'S6'}    
                % l= 160 mm focalisant
                % Etalonnage HL(I) sur 80 - 250 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  0.0;
                a3=  -2.6735E-8;
                a2=  5.8793E-6;
                a1=  1.5364E-1;
                a0=  2.7867E-1;
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT';

                
              case {'S4','S8'}    
                % l= 160 mm focalisants
                % Etalonnage HL(I) sur 170 - 300 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  -8.8836E-10;
                a3=  7.1089E-7;
                a2=  -2.2277E-4;
                a1=  1.8501E-1;
                a0=  -1.329;
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT';

             case {'S2','S5'}    
                % l= 160 mm dï¿œfocalisants
                % Etalonnage HL(I) sur 170 - 300 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  0.0;
                a4=  -8.8836E-10*(-1);
                a3=  7.1089E-7;
                a2=  -2.2277E-4*(-1);
                a1=  1.8501E-1;
                a0=  -1.329*(-1);
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT';    
        
              case {'S7'}    
                % l= 160 mm dï¿œfocalisant
                % Etalonnage HL(I) sur 250 - 350 A
                % Find the current from the given polynomial for B''Leff
                Leff=1e-8; % modeled as thin length;
                a7=  0.0;
                a6=  0.0;
                a5=  -2.613556E-10;
                a4=  3.730258E-7*(-1);
                a3=  -2.1301205E-4;
                a2=  6.077561E-2*(-1);
                a1=  -8.5069349;
                a0=  4.933E+2*(-1);
                A = [a7 a6 a5 a4 a3 a2 a1 a0]*2;
                MagnetType = 'SEXT'; 
                
            case 'QT'    % 160 mm dans sextupole
                % Etalonnage: moyenne sur les 32 sextupï¿œles incluant un QT.
                % Efficacite = 3 G.m/A @ R=32mm; soit 93.83 G/A
                % Le signe du courant est donnￜ par le DeviceServer (Tango)
                % Find the currAO.(ifam).Monitor.HW2PhysicsParams{1}(1,:) = magnetcoefficients(AO.(ifam).FamilyName );
                Leff = 1e-8;
                a7= 0.0;
                a6= 0.0;
                a5= 0.0;
                a4= 0.0;
                a3= 0.0;
                a2= 0.0;
                a1= 93.83E-4;
                a0= 0.0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];

                MagnetType = 'QT';
                
            case 'SQ'   % 160 mm dans sextupole
                % Etalonnage: moyenne sur les 32 sextupï¿œles incluant un QT.
                % Efficacitee = 3 G.m/A @ R=32mm; soit 93.83 G/A
                % Le signe du courant est donnee par le DeviceServer (Tango)
                % Find the currAO.(ifam).Monitor.HW2PhysicsParams{1}(1,:) = magnetcoefficients(AO.(ifam).FamilyName );
                Leff = 1e-8;
                a7= 0.0;
                a6= 0.0;
                a5= 0.0;
                a4= 0.0;
                a3= 0.0;
                a2= 0.0;
                a1= 93.83E-4;
                a0= 0.0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'QT';
                
            case {'HCOR'}    % 16 cm horizontal corrector
                % Etalonnage: moyenne sur les 56 sextupï¿œles incluant un CORH.
                % Efficacitￜ = 8.143 G.m/A
                % Le signe du courant est donnￜ par le DeviceServer (Tango)
                % Find the currAO.(ifam).Monitor.HW2PhysicsParams{1}(1,:) = magnetcoefficients(AO.(ifam).FamilyName );
                Leff = 0.16;
                a7= 0.0;
                a6= 0.0;
                a5= 0.0;
                a4= 0.0;
                a3= 0.0;
                a2= 0.0;
                a1= 8.143E-4;
                a0= 0.0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'COR';
                
                
            case {'FHCOR'}    % 10 cm horizontal corrector
                % Magnet Spec: Theta = 280e-6 radians @ 2.75 GeV and 10 amps
                % Theta = BLeff / Brho    [radians]
                % Therefore,
                %       Theta = ((BLeff/Amp)/ Brho) * I
                %       BLeff/Amp = 280e-6 * getbrho(2.75) / 10
                %       B*Leff = a0 * I   => a0 = 0.8e-3 * getbrho(2.75) / 10
                %
                % The C coefficients are w.r.t B
                %       B = c0 + c1*I = (0 + a0*I)/Leff
                % However, AT uses Theta in radians so the A coefficients
                % must be used for correctors with the middle layer with
                % the addition of the DC term
                
                % Find the current from the given polynomial for BLeff and B
                % NOTE: AT used BLeff (A) for correctors
                Leff = .10;
                imax = 10;
                cormax = 28e-6 ; % 28 urad for imax = 10 A
                MagnetType = 'COR';
                A = [0 cormax*getbrho(2.75)/imax 0];
                
            case {'VCOR'}    % 16 cm vertical corrector
                % Etalonnage: moyenne sur les 56 sextupï¿œles incluant un CORV.
                % Efficacitￜ = 4.642 G.m/A
                % Le signe du courant est donnￜ par le DeviceServer (Tango)
                % Find the currAO.(ifam).Monitor.HW2PhysicsParams{1}(1,:) = magnetcoefficients(AO.(ifam).FamilyName );
                Leff = 0.16;
                a7= 0.0;
                a6= 0.0;
                a5= 0.0;
                a4= 0.0;
                a3= 0.0;
                a2= 0.0;
                a1= 4.642E-4;
                a0= 0.0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
                MagnetType = 'COR';
                
            case {'FVCOR'}    % 10 cm vertical corrector
                % Find the current from the given polynomial for BLeff and B
                Leff = .10;
                imax = 10;
                cormax = 23e-6 ; % 23 urad for imax = 10 A
                MagnetType = 'COR';
                A = [0 cormax*getbrho(2.75)/imax 0];

            case {'K_INJ'}
                % Kicker d'injection
                % étalonnage provisoire
                % attention l'element n'etant pas dans le modele,definition
                % de A ambigue
                Leff = .6;
                vmax = 8000;
                alphamax = 8e-3 ; % 8 mrad pour 8000 V
                MagnetType = 'K_INJ';
                A = [0 alphamax*getbrho(2.75)/vmax 0]*Leff;
                
             case {'K_INJ1'}
                % Kickers d'injection 1 et 4
                Leff = .6;
                vmax = 7500; % tension de mesure
                SBDL = 75.230e-3 ; % somme de Bdl mesurée
                MagnetType = 'K_INJ1';
                A = [0 -SBDL/vmax 0]*Leff; 
                
             case {'K_INJ2'}
                % Kickers d'injection 2 et 3
                Leff = .6;
                vmax = 7500;% tension de mesure
                SBDL = 74.800e-3 ; % somme de Bdl mesurée
                MagnetType = 'K_INJ2';  
                A = [0 SBDL/vmax 0]*Leff; 
                
            case {'SEP_P'}
                % Septum passif d'injection
                Leff = .6;
                vmax = 547; % tension de mesure V
                SBDL = 263e-3; % somme de Bdl mesurée
                MagnetType = 'SEP_P';
                A = [0 SBDL/vmax 0]*Leff; 
                
             case {'SEP_A'}
                % Septum actif d'injection
                Leff = 1.;
                vmax = 111;
                MagnetType = 'SEP_A';
                SBDL = 1147.8e-3 ; % Somme de Bdl mesurée à 111 V
                A = [0 SBDL/vmax 0]*Leff; 

             case {'KEMH'}
                % Horizontal pinger (KEM)
                Leff = 1.;
                vmax = 1000;
                MagnetType = 'KEMH';
                SBDL = 0.13e-3 ; 
                A = [0 SBDL/vmax 0]*Leff; 

             case {'KEMV'}
                % vertical pinger (KEM)
                Leff = 1.;
                vmax = 1000;
                MagnetType = 'KEMV';
                SBDL = 0.055e-3 ; 
                A = [0 SBDL/vmax 0]*Leff; 

            otherwise
                error(sprintf('MagnetCoreType %s is not unknown', MagnetCoreType));
                k = 0;
                MagnetType = '';
                return
        end

        % compute B-field = int(Bdl)/Leff
        C = A / Leff;

        MagnetType = upper(MagnetType);


        % Power Series Denominator (Factoral) be AT compatible
        if strcmpi(MagnetType,'SEXT')
            C = C / 2;
        end
        if strcmpi(MagnetType,'OCTO')
            C = C / 6;
        end
        return;
    case 'Booster'
        %%%%
        switch upper(deblank(MagnetCoreType))

            case 'BEND'    
                % B[T] = 0.00020 + 0.0013516 I[A]
                % B[T] = 0.00020 + (0.0013051 + 0.00005/540 I) I[A] Alex
                Leff = 2.160; % 2160 mm                
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  9.2e-8*Leff;
                a1 =  0.0013051*Leff;
                a0 =  2.0e-3*Leff;

                A = [a8 a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'BEND';

            case {'QF'}   % 400 mm quadrupole
                % Find the current from the given polynomial for B'Leff
                % G[T/m] = 0.0465 + 0.0516 I[A] Alex
                Leff=0.400; 
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  0.0516*Leff;
                a0 =  0.0465*Leff;
                
                A = [a7 a6 a5 a4 a3 a2 a1 a0]; %*getbrho(0.1);
                MagnetType = 'QUAD';

            case {'QD'}   % 400 mm quadrupole
                % Find the current from the given polynomial for B'Leff
                % G[T/m] = 0.0485 + 0.0518 I[A] Alex
                Leff=0.400; 
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  -0.0518*Leff;
                a0 =  -0.0485*Leff;
                
                A = [a7 a6 a5 a4 a3 a2 a1 a0]; %*getbrho(0.1);
                MagnetType = 'QUAD';

            case {'SF', 'SD'}   % 150 mm sextupole
                % Find the current from the given polynomial for B'Leff
                % HL [T/m] = 0.2 I [A] (deja intￜgrￜ)
                Leff=1.e-8; % thin lens;
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  0.2*2;
                a0 =  0.0;
                
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                MagnetType = 'SEXT';
                
            case {'HCOR','VCOR'}    % ?? cm horizontal corrector
                % Magnet Spec: Theta = 0.8e-3 radians @ 2.75 GeV and 10 amps
                % Theta = BLeff / Brho    [radians]
                % Therefore,
                %       Theta = ((BLeff/Amp)/ Brho) * I
                %       BLeff/Amp = 0.8e-3 * getbrho(2.75) / 10
                %       B*Leff = a0 * I   => a0 = 0.8e-3 * getbrho(2.75) / 10
                %
                % The C coefficients are w.r.t B
                %       B = c0 + c1*I = (0 + a0*I)/Leff
                % However, AT uses Theta in radians so the A coefficients
                % must be used for correctors with the middle layer with
                % the addition of the DC term

                % Find the current from the given polynomial for BLeff and B
                % NOTE: AT used BLeff (A) for correctors
                MagnetType = 'COR';
                % theta [mrad] = 1.34 I[A] @ 0.1 GeV
                Leff = 1e-6;
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  1.34e-3*getbrho(0.1);
                a0 =  0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
            otherwise
                error(sprintf('MagnetCoreType %s is not unknown', MagnetCoreType));
                %k = 0;
                %MagnetType = '';
                %return
        end

        % compute B-field = int(Bdl)/Leff
        C = A/ Leff;
 
        % Power Series Denominator (Factoral) be AT compatible
        if strcmpi(MagnetType,'SEXT')
            C = C / 2;
        end
 
        MagnetType = upper(MagnetType);

    case 'LT2'
        %%%%
        switch upper(deblank(MagnetCoreType))

            case 'BEND'    
                % les coefficients et longueur magnétique sont recopiés de l'anneau
                Leff=1.052433;
                a7= 0.0;
                a6=-0.0;
                a5= 0.0;
                a4=-0.0;
                a3= 0.0;
                a2=-9.7816E-6*(1-1.8e-3)*Leff*(1.055548/1.052433);
                a1= 1.26066E-02*(1-1.8E-3)*Leff*(1.055548/1.052433);
                a0= -2.24944*(1-1.8E-3)*Leff*(1.055548/1.052433);
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                

                MagnetType = 'BEND';

            case {'QP'}   % 400 mm quadrupole
                % Find the current from the given polynomial for B'Leff
                
                % G[T/m] = 0.1175 + 0.0517 I[A]
                % le rémanent est + fort que pour les quad Booster car les
                % courants max sont + eleves
                Leff=0.400; 
%                 a8 =  0.0;
%                 a7 =  0.0;
%                 a6 =  0.0;
%                 a5 =  0.0;
%                 a4 =  0.0;
%                 a3 =  0.0;
%                 a2 =  0.0;
%                 a1 =  0.0517*Leff;
%                 a0 =  0.1175*Leff;
                
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  -1.3345e-10;
                a3 =  8.1746e-8;
                a2 =  -1.6548e-5;
                a1 =  2.197e-2;
                a0 =  2.73e-2;
                A = [a7 a6 a5 a4 a3 a2 a1 a0]; 
                MagnetType = 'QUAD';

            case {'CH','CV'}    % 16 cm horizontal corrector
                

                
                % Magnet Spec: Theta = environ 1 mradians @ 2.75 GeV and 10 amps
                % Theta = BLeff / Brho    [radians]
                % Therefore,
                %       Theta = ((BLeff/Amp)/ Brho) * I
                %       BLeff/Amp = 1.e-3 * getbrho(2.75) / 10
                %       B*Leff = a1 * I   => a1 = 1.e-3 * getbrho(2.75) / 10
                %
                % The C coefficients are w.r.t B
                %       B = c0 + c1*I = (0 + a0*I)/Leff
                % However, AT uses Theta in radians so the A coefficients
                % must be used for correctors with the middle layer with
                % the addition of the DC term

                % Find the current from the given polynomial for BLeff and B
                % NOTE: AT used BLeff (A) for correctors
                
                % environ 32 cm  corrector
                % Efficacitￜ = 11.06 G.m/A 
                % Le signe du courant est donnￜ par le DeviceServer (Tango)
                % Find the currAO.(ifam).Monitor.HW2PhysicsParams{1}(1,:) =
                % magnetcoefficien
                
                MagnetType = 'COR';
                
                Leff = 1e-6; % 0.1577 m
                a8 =  0.0;
                a7 =  0.0;
                a6 =  0.0;
                a5 =  0.0;
                a4 =  0.0;
                a3 =  0.0;
                a2 =  0.0;
                a1 =  110.6e-4/10;
                a0 =  0;
                A = [a7 a6 a5 a4 a3 a2 a1 a0];
                
            otherwise
                error(sprintf('MagnetCoreType %s is not unknown', MagnetCoreType));
                %k = 0;
                %MagnetType = '';
                %return
        end

        % compute B-field = int(Bdl)/Leff
        C = A/ Leff;

        MagnetType = upper(MagnetType);

    otherwise
        error('Unknown accelerator name %s', AcceleratorName);
end