PageRenderTime 28ms CodeModel.GetById 11ms RepoModel.GetById 0ms app.codeStats 0ms

/src/utils/get_DR_quads.m

http://slac-lucretia.googlecode.com/
MATLAB | 326 lines | 183 code | 35 blank | 108 comment | 20 complexity | 42f359c5119e701467ff82d369d9cda6 MD5 | raw file
Possible License(s): BSD-2-Clause 
    

  1. function varargout=get_DR_quads(opCode,energy,varargin)

    2. 

  2. %

    3. 

  3. % [Iquad,KL]=get_DR_quads(1,energy,Imain,Itrim,fflag);

    4. 

  4. %

    5. 

  5. % Compute "effective" main coil currents and normalized integrated strengths

    6. 

  6. % (GL/brho) for individual DR quadrupoles from main coil currents, trim coil

    7. 

  7. % currents, and beam energy

    8. 

  8. %

    9. 

  9. % Iquad=get_DR_quads(2,energy,KL,fflag);

    10. 

  10. %

    11. 

  11. % Compute main coil currents for DR quadrupole families from normalized

    12. 

  12. % integrated strengths (GL/brho) and beam energy

    13. 

  13. %

    14. 

  14. % INPUTs (opCode=1):

    15. 

  15. %

    16. 

  16. %   opCode = 1 (I to KL conversion)

    17. 

  17. %   energy = DR energy (GeV) [scalar]

    18. 

  18. %   Imain  = DR quadrupole family main coil currents (amps) [26 element array]

    19. 

  19. %   Itrim  = DR quadrupole trim coil currents (amps) [100 element array]

    20. 

  20. %   fflag  = (optional) if present and nonzero, fudge factors will be used

    21. 

  21. %

    22. 

  22. % INPUTs (opCode=2):

    23. 

  23. %

    24. 

  24. %   opCode = 2 (KL to I conversion)

    25. 

  25. %   energy = DR energy (GeV) [scalar]

    26. 

  26. %   KL     = normalized integrated strengths (1/m; positive means "QF")

    27. 

  27. %             [26 element array]

    28. 

  28. %   fflag  = (optional) if present and nonzero, fudge factors will be used

    29. 

  29. %

    30. 

  30. % The order of quadrupole families in the Imain or KL arrays must be:

    31. 

  31. %

    32. 

  32. %   QF1R,QF2R,

    33. 

  33. %   QM1R,QM2R,QM3R,QM4R,QM5R,

    34. 

  34. %   QM6R,QM7R,QM8R,QM9R,QM10R,

    35. 

  35. %   QM11R,QM12R,QM13R,QM14R,QM15R,

    36. 

  36. %   QM16R,QM17R,QM18R,QM19R,QM20R,

    37. 

  37. %   QM21R,QM22R,QM23R,QM7AR

    38. 

  38. %

    39. 

  39. % The order of quadrupole trims in the Itrim array must be:

    40. 

  40. %

    41. 

  41. %   QF1R.1:28,QF2R.1:26,

    42. 

  42. %   QM1R.1:2,QM2R.1:2,QM3R.1:2,QM4R.1:2,QM5R.1:2,

    43. 

  43. %   QM6R.1:2,QM7R.1:2,QM8R.1:2,QM9R.1:2,QM10R.1:2,

    44. 

  44. %   QM11R.1:2,QM12R.1:2,QM13R.1:2,QM14R.1:2,QM15R.1:2,

    45. 

  45. %   QM16R.1:2,QM17R.1:2,QM18R.1:2,QM19R.1:2,QM20R.1:2,

    46. 

  46. %   QM21R.1:2,QM22R.1:2,QM23R.1:2

    47. 

  47. %

    48. 

  48. % OUTPUTs (opCode=1):

    49. 

  49. %

    50. 

  50. %   Iquad = effective DR quadrupole currents (amps) [100 element array,

    51. 

  51. %            same order as Itrim]

    52. 

  52. %   KL    = normalized integrated strengths (1/m; positive means "QF")

    53. 

  53. %            [26 element array]

    54. 

  54. %

    55. 

  55. % OUTPUT (opCode=2):

    56. 

  56. %

    57. 

  57. %   Iquad = DR quadrupole family currents (amps) [26 element array]

    58. 

  58. 

    59. 

  59. % ==============================================================================

    60. 

  60. % 26-FEB-2009, M. Woodley

    61. 

  61. %    Correct Nm(=49) and Nt(=20) for QEA quads per IHEP drawing

    62. 

  62. % 22-FEB-2009, M. Woodley

    63. 

  63. %    Add support for QM7AR and QEAs; do either I-to-KL and KL-to-I conversion

    64. 

  64. % 22-MAY-2007, M. Woodley

    65. 

  65. %    Data for IHEP "QEA" quads (QM12R, QM13R, QM14R) from control system file

    66. 

  66. % 06-FEB-2007, M. Woodley

    67. 

  67. %    Tokin 3581 quads replaced with IHEP "QEA" quads (QM12R, QM13R, QM14R)

    68. 

  68. % ==============================================================================

    69. 

  69. 

    70. 

  70. Nmain=26;

    71. 

  71. Ntrim=100;

    72. 

  72. 

    73. 

  73. % check input args

    74. 

  74. 

    75. 

  75. if (~exist('opCode','var'))

    76. 

  76.   error('opCode argument not supplied')

    77. 

  77. end

    78. 

  78. if (~exist('energy','var'))

    79. 

  79.   error('energy argument not supplied')

    80. 

  80. else

    81. 

  81.   if (energy<=0)

    82. 

  82.     error('Bad energy value')

    83. 

  83.   end

    84. 

  84. end

    85. 

  85. if (opCode==1)

    86. 

  86.   if (nargin<4)

    87. 

  87.     error('At least 4 input arguments required for opCode=1')

    88. 

  88.   end

    89. 

  89.   Imain=varargin{1};

    90. 

  90.   if (length(Imain)~=Nmain)

    91. 

  91.     error('Incorrect Imain length')

    92. 

  92.   end

    93. 

  93.   Itrim=varargin{2};

    94. 

  94.   if (length(Itrim)~=Ntrim)

    95. 

  95.     error('Incorrect Itrim length')

    96. 

  96.   end

    97. 

  97.   if (nargin>4)

    98. 

  98.     fflag=varargin{3};

    99. 

  99.   else

    100. 

  100.     fflag=0;

    101. 

  101.   end

    102. 

  102. elseif (opCode==2)

    103. 

  103.   if (nargin<3)

    104. 

  104.     error('At least 3 input arguments required for opCode=2')

    105. 

  105.   end

    106. 

  106.   KL=varargin{1};

    107. 

  107.   if (length(KL)~=Nmain)

    108. 

  108.     error('Incorrect KL length')

    109. 

  109.   end

    110. 

  110.   if (nargin>3)

    111. 

  111.     fflag=varargin{2};

    112. 

  112.   else

    113. 

  113.     fflag=0;

    114. 

  114.   end

    115. 

  115. else

    116. 

  116.   error('Invalid opCode')

    117. 

  117. end

    118. 

  118. fudge=(fflag~=0);

    119. 

  119. 

    120. 

  120. % ==============================================================================

    121. 

  121. % DR QUAD MAGNET DATA

    122. 

  122. % ==============================================================================

    123. 

  123. 

    124. 

  124. % there are 11 DR quadrupole types:

    125. 

  125. 

    126. 

  126. %    1 = Hitachi type 1

    127. 

  127. %    2 = Hitachi type 2

    128. 

  128. %    3 = Hitachi type 3

    129. 

  129. %    4 = Hitachi type 4

    130. 

  130. %    5 = Tokin type 3325 ("18 cm")

    131. 

  131. %    6 = Tokin type 3582 ("18 cm (42p)")

    132. 

  132. %    7 = Tokin type 3393 ("6 cm")

    133. 

  133. %    8 = Tokin type 3581 ("6 cm")

    134. 

  134. %    9 = QEA-13 (IHEP type D32L180 # 13) ... in series with QEA-14

    135. 

  135. %   10 = QEA-19 (IHEP type D32L180 # 19) ... in series with QEA-17

    136. 

  136. %   11 = QEA-21 (IHEP type D32L180 # 21) ... in series with QEA-18

    137. 

  137. 

    138. 

  138. % magnet data  (I is amp, G is T/m, leff is m)

    139. 

  139. % (from: ATF$MAG:MAG_KI_Q_HITACHI_1.FOR

    140. 

  140. %        ATF$MAG:MAG_KI_Q_HITACHI_2.FOR

    141. 

  141. %        ATF$MAG:MAG_KI_Q_HITACHI_3.FOR

    142. 

  142. %        ATF$MAG:MAG_KI_Q_HITACHI_4.FOR

    143. 

  143. %        ATF$MAG:MAG_KI_Q_TOKIN_3325.FOR

    144. 

  144. %        ATF$MAG:MAG_KI_Q_TOKIN_3582.FOR

    145. 

  145. %        ATF$MAG:MAG_KI_Q_TOKIN_3393.FOR

    146. 

  146. %        ATF$MAG:MAG_KI_Q_TOKIN_3581.FOR

    147. 

  147. %        mainamp50AmaxAbs.txt (M.Masuzawa))

    148. 

  148. 

    149. 

  149. Nmms=[11,11,11,11,11,11,8,11,16,16,16];

    150. 

  150. 

    151. 

  151. qI=[ ...

    152. 

  152.     0.0,  0.0,  0.0,  0.0,  0.00,  0.0,  0.0,  0.0,0.0000000,0.0000000,0.0000000; ...

    153. 

  153.    14.2, 10.2, 20.4, 20.2,100.00, 50.0, 20.0, 50.0,3.3333001,3.3333001,3.3333001; ...

    154. 

  154.    28.6, 20.2, 40.0, 40.2,200.00,100.0, 40.0, 80.0,6.6666999,6.6666999,6.6666999; ...

    155. 

  155.    42.2, 30.2, 60.0, 60.4,300.00,150.0, 60.0,100.0,10.000000,10.000000,10.000000; ...

    156. 

  156.    56.0, 40.0, 80.2, 80.4,400.00,170.0, 80.0,130.0,13.333000,13.333000,13.333000; ...

    157. 

  157.    70.2, 50.0,100.8,100.0,417.24,190.0,100.0,150.0,16.667000,16.667000,16.667000; ...

    158. 

  158.    84.2, 60.2,120.2,120.2,434.48,210.0,120.0,170.0,20.000000,20.000000,20.000000; ...

    159. 

  159.    98.2, 70.2,140.0,140.2,451.72,220.0,139.0,190.0,23.333000,23.333000,23.333000; ...

    160. 

  160.   112.2, 80.2,160.2,160.2,468.97,230.0,  0  ,210.0,26.667000,26.667000,26.667000; ...

    161. 

  161.   126.0, 90.2,180.2,180.0,503.45,240.0,  0  ,230.0,30.000000,30.000000,30.000000; ...

    162. 

  162.   140.2,100.2,200.2,200.4,512.07,245.0,  0  ,245.0,33.333000,33.333000,33.333000; ...

    163. 

  163.     0  ,  0  ,  0  ,  0  ,  0   ,  0  ,  0  ,  0  ,36.667000,36.667000,36.667000; ...

    164. 

  164.     0  ,  0  ,  0  ,  0  ,  0   ,  0  ,  0  ,  0  ,40.000000,40.000000,40.000000; ...

    165. 

  165.     0  ,  0  ,  0  ,  0  ,  0   ,  0  ,  0  ,  0  ,43.333000,43.333000,43.333000; ...

    166. 

  166.     0  ,  0  ,  0  ,  0  ,  0   ,  0  ,  0  ,  0  ,46.667000,46.667000,46.667000; ...

    167. 

  167.     0  ,  0  ,  0  ,  0  ,  0   ,  0  ,  0  ,  0  ,50.000000,50.000000,50.000000; ...

    168. 

  168. ];                                                 

    169. 

  169. 

    170. 

  170. qG=[ ...

    171. 

  171.    0.000, 0.000, 0.000, 0.000, 0.00, 0.0, 0.0, 0.0,0         ,0         ,0         ; ...

    172. 

  172.    2.350, 3.832, 5.628, 4.612,10.80, 7.6, 3.4, 7.5,0.32538533,0.33133936,0.32236442; ...

    173. 

  173.    4.712, 7.646,11.109, 9.201,21.45,14.9, 6.8,11.8,0.63208896,0.63420987,0.62648714; ...

    174. 

  174.    6.981,11.498,16.656,13.858,32.05,22.2,10.0,14.8,0.94148004,0.93937200,0.93356729; ...

    175. 

  175.    9.256,15.202,22.243,18.440,42.65,25.0,13.2,19.1,1.2529521 ,1.2468803 ,1.2427585 ; ...

    176. 

  176.   11.604,18.927,27.866,22.884,44.70,28.0,16.6,21.8,1.5639043 ,1.5545483 ,1.5514039 ; ...

    177. 

  177.   13.892,22.657,33.122,27.464,46.65,30.8,19.8,24.4,1.8741443 ,1.8627108 ,1.8593891 ; ...

    178. 

  178.   16.167,26.285,38.487,31.925,48.10,32.1,23.0,26.4,2.1841798 ,2.1716371 ,2.1672122 ; ...

    179. 

  179.   18.410,29.896,43.828,36.386,49.70,33.3, 0  ,28.0,2.5015478 ,2.4888871 ,2.4827976 ; ...

    180. 

  180.   20.629,33.438,48.575,40.726,52.35,34.2, 0  ,29.3,2.8104367 ,2.7979240 ,2.7898581 ; ...

    181. 

  181.   22.898,36.758,51.971,45.143,52.80,34.7, 0  ,30.1,3.1178350 ,3.1058538 ,3.0956228 ; ...

    182. 

  182.    0    , 0    , 0    , 0    , 0   , 0  , 0  , 0  ,3.4247386 ,3.4134352 ,3.4008279 ; ...

    183. 

  183.    0    , 0    , 0    , 0    , 0   , 0  , 0  , 0  ,3.7313569 ,3.7210751 ,3.7056735 ; ...

    184. 

  184.    0    , 0    , 0    , 0    , 0   , 0  , 0  , 0  ,4.0364571 ,4.0274525 ,4.0088239 ; ...

    185. 

  185.    0    , 0    , 0    , 0    , 0   , 0  , 0  , 0  ,4.3408098 ,4.3334475 ,4.3110003 ; ...

    186. 

  186.    0    , 0    , 0    , 0    , 0   , 0  , 0  , 0  ,4.6463265 ,4.6411600 ,4.6146355 ; ...

    187. 

  187. ];

    188. 

  188. 

    189. 

  189. qleff=[0.078765,0.07867,0.19847,0.198745,0.19886,0.202628,0.07890677,0.084339,0.19849,0.19849,0.19849];

    190. 

  190. 

    191. 

  191. Nm=[17,39,29,24,11,26,17,26,49,49,49]; % main coil turns

    192. 

  192. Nt=[20,20,20,20,20,20,20,20,20,20,20]; % trim coil turns

    193. 

  193. 

    194. 

  194. % NOTE: data for types 9-11 (IHEP type D32L180 ("QEA")) are integrated strength

    195. 

  195. %       in Tesla ... divide by leff for standard handling

    196. 

  196. 

    197. 

  197. for n=9:11

    198. 

  198.   qG(:,n)=qG(:,n)/qleff(n);

    199. 

  199. end

    200. 

  200. 

    201. 

  201. % quadrupole family names

    202. 

  202. 

    203. 

  203. qname=['QF1R ';'QF2R ';'QM1R ';'QM2R ';'QM3R ';'QM4R ';'QM5R ';'QM6R '; ...

    204. 

  204.        'QM7R ';'QM8R ';'QM9R ';'QM10R';'QM11R';'QM12R';'QM13R';'QM14R'; ...

    205. 

  205.        'QM15R';'QM16R';'QM17R';'QM18R';'QM19R';'QM20R';'QM21R';'QM22R'; ...

    206. 

  206.        'QM23R';'QM7AR'];

    207. 

  207. 

    208. 

  208. % quadrupole family types

    209. 

  209. % (from: ATF$MAG:MAG_KI_MAIN.FOR)

    210. 

  210. 

    211. 

  211. qtype=[1;5;3;2;2;3;4;4;7;7;4;4;6;9;10;11;6;4;4;4;4;3;2;2;3;8];

    212. 

  212. 

    213. 

  213. % polarities

    214. 

  214. 

    215. 

  215. qsgn=[1;1;1;-1;-1;1;1;-1;1;1;-1;1;-1;1;-1;1;-1;1;-1;-1;1;1;-1;-1;1;1];

    216. 

  216. 

    217. 

  217. % "Kubo" fudge factors

    218. 

  218. % (NOTE: set fudge factors to zero ... who knows what they are now)

    219. 

  219. 

    220. 

  220. qfudge=[ ...

    221. 

  221.   0.0; ... % QF1R

    222. 

  222.   0.0; ... % QF2R

    223. 

  223.   0.0; ... % QM1R

    224. 

  224.   0.0; ... % QM2R

    225. 

  225.   0.0; ... % QM3R

    226. 

  226.   0.0; ... % QM4R

    227. 

  227.   0.0; ... % QM5R

    228. 

  228.   0.0; ... % QM6R

    229. 

  229.   0.0; ... % QM7R

    230. 

  230.   0.0; ... % QM8R

    231. 

  231.   0.0; ... % QM9R

    232. 

  232.   0.0; ... % QM10R

    233. 

  233.   0.0; ... % QM11R

    234. 

  234.   0.0; ... % QM12R

    235. 

  235.   0.0; ... % QM13R

    236. 

  236.   0.0; ... % QM14R

    237. 

  237.   0.0; ... % QM15R

    238. 

  238.   0.0; ... % QM16R

    239. 

  239.   0.0; ... % QM17R

    240. 

  240.   0.0; ... % QM18R

    241. 

  241.   0.0; ... % QM19R

    242. 

  242.   0.0; ... % QM20R

    243. 

  243.   0.0; ... % QM21R

    244. 

  244.   0.0; ... % QM22R

    245. 

  245.   0.0; ... % QM23R

    246. 

  246.   0.0; ... % QM7AR

    247. 

  247. ];

    248. 

  248. 

    249. 

  249. % quadrupole family pointers

    250. 

  250. 

    251. 

  251. idq=[ ...

    252. 

  252.    1*ones(28,1); ... % QF1R

    253. 

  253.    2*ones(26,1); ... % QF2R

    254. 

  254.    3*ones( 2,1); ... % QM1R

    255. 

  255.    4*ones( 2,1); ... % QM2R

    256. 

  256.    5*ones( 2,1); ... % QM3R

    257. 

  257.    6*ones( 2,1); ... % QM4R

    258. 

  258.    7*ones( 2,1); ... % QM5R

    259. 

  259.    8*ones( 2,1); ... % QM6R

    260. 

  260.          [26;9]; ... % QM7AR,QM7R

    261. 

  261.   10*ones( 2,1); ... % QM8R

    262. 

  262.   11*ones( 2,1); ... % QM9R

    263. 

  263.   12*ones( 2,1); ... % QM10R

    264. 

  264.   13*ones( 2,1); ... % QM11R

    265. 

  265.   14*ones( 2,1); ... % QM12R

    266. 

  266.   15*ones( 2,1); ... % QM13R

    267. 

  267.   16*ones( 2,1); ... % QM14R

    268. 

  268.   17*ones( 2,1); ... % QM15R

    269. 

  269.   18*ones( 2,1); ... % QM16R

    270. 

  270.   19*ones( 2,1); ... % QM17R

    271. 

  271.   20*ones( 2,1); ... % QM18R

    272. 

  272.   21*ones( 2,1); ... % QM19R

    273. 

  273.   22*ones( 2,1); ... % QM20R

    274. 

  274.   23*ones( 2,1); ... % QM21R

    275. 

  275.   24*ones( 2,1); ... % QM22R

    276. 

  276.   25*ones( 2,1); ... % QM23R

    277. 

  277. ];

    278. 

  278. 

    279. 

  279. % ==============================================================================

    280. 

  280. 

    281. 

  281. % compute rigidity

    282. 

  282. 

    283. 

  283. Cb=1e10/2.99792458e8; % kG-m/GeV

    284. 

  284. brho=Cb*energy;

    285. 

  285. 

    286. 

  286. % perform requested conversion

    287. 

  287. % (NOTE: polynomial evaluation is not presently used to convert current to

    288. 

  288. %        gradient ... linear interpolation of MMS data is used)

    289. 

  289. 

    290. 

  290. if (opCode==1)

    291. 

  291.   

    292. 

  292. % convert current to KL

    293. 

  293.   

    294. 

  294.   Iquad=zeros(Ntrim,1);

    295. 

  295.   KL=zeros(Ntrim,1);

    296. 

  296.   for n=1:Ntrim

    297. 

  297.     m=idq(n); % quad family (1-26)

    298. 

  298.     t=qtype(m); % quad type (1-11)

    299. 

  299.     Iquad(n)=Imain(m)+(Nt(t)/Nm(t))*Itrim(n);

    300. 

  300.     if (Iquad(n)<=0)

    301. 

  301.       G=0;

    302. 

  302.     else

    303. 

  303.       G=interp1(qI(1:Nmms(t),t),qG(1:Nmms(t),t),Iquad(n),'linear');

    304. 

  304.     end

    305. 

  305.     KL(n)=qsgn(m)*(10*G)*qleff(t)/brho;

    306. 

  306.     if (fudge),KL(n)=KL(n)/(1+qfudge(m));end

    307. 

  307.   end

    308. 

  308.   varargout{1}=Iquad;

    309. 

  309.   varargout{2}=KL;

    310. 

  310. else

    311. 

  311.   

    312. 

  312. % convert KL to current

    313. 

  313.   

    314. 

  314.   Iquad=zeros(Nmain,1);

    315. 

  315.   for n=1:Nmain

    316. 

  316.     if (KL(n)~=0)

    317. 

  317.       t=qtype(n); % quad type (1-11)

    318. 

  318.       G=0.1*brho*abs(KL(n))/qleff(t); % T/m

    319. 

  319.       if (fudge),G=G*(1+qfudge(t));end

    320. 

  320.       Iquad(n)=interp1(qG(1:Nmms(t),t),qI(1:Nmms(t),t),G,'linear');

    321. 

  321.     end

    322. 

  322.   end

    323. 

  323.   varargout{1}=Iquad;

    324. 

  324. end

    325. 

  325. 

    326. 

  326. end

    327. 

  327.