[4] | 1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" |
---|
| 2 | "http://www.w3.org/TR/REC-html40/loose.dtd"> |
---|
| 3 | <html> |
---|
| 4 | <head> |
---|
| 5 | <title>Description of setorbitbump</title> |
---|
| 6 | <meta name="keywords" content="setorbitbump"> |
---|
| 7 | <meta name="description" content="SETORBITBUMP - Local bump program (uses setorbit)"> |
---|
| 8 | <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> |
---|
| 9 | <meta name="generator" content="m2html © 2003 Guillaume Flandin"> |
---|
| 10 | <meta name="robots" content="index, follow"> |
---|
| 11 | <link type="text/css" rel="stylesheet" href="../m2html.css"> |
---|
| 12 | </head> |
---|
| 13 | <body> |
---|
| 14 | <a name="_top"></a> |
---|
| 15 | <div><a href="../index.html">Home</a> > <a href="index.html">mml</a> > setorbitbump.m</div> |
---|
| 16 | |
---|
| 17 | <!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td> |
---|
| 18 | <td align="right"><a href="index.html">Index for mml <img alt=">" border="0" src="../right.png"></a></td></tr></table>--> |
---|
| 19 | |
---|
| 20 | <h1>setorbitbump |
---|
| 21 | </h1> |
---|
| 22 | |
---|
| 23 | <h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 24 | <div class="box"><strong>SETORBITBUMP - Local bump program (uses setorbit)</strong></div> |
---|
| 25 | |
---|
| 26 | <h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 27 | <div class="box"><strong>function [OCS, OCS0, V, S, ErrorFlag] = setorbitbump(varargin) </strong></div> |
---|
| 28 | |
---|
| 29 | <h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 30 | <div class="fragment"><pre class="comment">SETORBITBUMP - Local bump program (uses setorbit) |
---|
| 31 | [OCS, OCS0, V, S] = setorbitbump(BPMFamily, BPMDeviceList, GoalOrbit, CMFamily, CMIncrementList, NIter, SVDIndex, BPMWeight) |
---|
| 32 | [OCS, OCS0, V, S] = setorbitbump(OCS, GoalOrbit, CMIncrementList, NIter, SVDIndex, BPMWeight) |
---|
| 33 | |
---|
| 34 | INPUTS |
---|
| 35 | 1. BPMFamily - BPM family |
---|
| 36 | 2. BPMDeviceList - BPM device list |
---|
| 37 | 3. GoalOrbit - Goal orbit position (Referenced to an absolute or incremental orbit |
---|
| 38 | depending on 'Incremental' or 'Absolute' flag) |
---|
| 39 | 4. CMFamily - Corrector magnet family |
---|
| 40 | 5. CMIncrementList - Magnet list upstream and downstream of the BPMs |
---|
| 41 | For instance, [-2 1 3] uses the 2nd correctors upstream and |
---|
| 42 | 1st and 3rd corrector downstream from the BPMs in the bump. |
---|
| 43 | or |
---|
| 44 | CMIncrementList can be an ordinary device list used for specifying correctors. |
---|
| 45 | 6. NIter - Number of iterations {Default: 3} |
---|
| 46 | (set NIter to 0 if no BPMs are functioning, see setorbit) |
---|
| 47 | 7. SVDIndex - vector, maximum number, 'All', or |
---|
| 48 | base on a threshold of min/max singular value {Default: see setorbit} |
---|
| 49 | 8. BPMWeight - Weight applied to the BPMs inside the bump (GoalOrbit). BPM weighting |
---|
| 50 | should only be required when the corrector are in a region of dispersion |
---|
| 51 | and the RF frequency is not used in the orbit correction. |
---|
| 52 | {Default BPMWeight is to weight the leakage twice as much as the GoalOrbit} |
---|
| 53 | 9. Optional flags (any flags used in setorbit can also be used here) |
---|
| 54 | 'Incremental' {Default} or 'Absolute' - Type of orbit change |
---|
| 55 | 'Display' - Display bump characteristics (no display is the default) |
---|
| 56 | 'SetPV' or 'SetSP' - Make the setpoint change {Default} |
---|
| 57 | 'NoSetPV' or 'NoSetSP' - Don't set the magnets, just return the OCS |
---|
| 58 | 'FitRF' or 'SetRF' - Flag to include the RF frequency as part of orbit correction. |
---|
| 59 | 'LeakageCorrection' - Option to try to clean up the leakage after the bump is applied. |
---|
| 60 | Note: 1. NIter for leakage correction is the same for the bump |
---|
| 61 | 2. It's just like setting the weight to zero for a second set of iterations. |
---|
| 62 | 3. For the horizontal plane, usually one should include the RF frequency |
---|
| 63 | along with leakage correction. Otherwise, the goal orbit may not be attained. |
---|
| 64 | 'Tolerance', Tol - Quit correction if the std(error) < Tol {Tol: 0} |
---|
| 65 | 'RampSteps', NSteps - Number of steps to ramp in the correction {NSteps: 1} |
---|
| 66 | This is used to avoid a large transient between setpoint changes. |
---|
| 67 | |
---|
| 68 | OUTPUTS (same as setorbit) |
---|
| 69 | 1. OCS - Orbit correction structure (OCS) usable by setorbit |
---|
| 70 | 2. OCS0 - Starting OCS |
---|
| 71 | 3. V - Corrector magnet singular vectors |
---|
| 72 | 4. Svalue - Singular values (vector) |
---|
| 73 | |
---|
| 74 | NOTES |
---|
| 75 | 1. setorbitbump creates an OCS structure and uses setorbit to actually change the orbit. |
---|
| 76 | |
---|
| 77 | EXAMPLES |
---|
| 78 | 1. 4 magnet horizontal incremental bump at BPM(6,6) and BPM(7,1) including the RF frequency |
---|
| 79 | OCS = setorbitbump('BPMx',[6 6;7 1],[1;-1],'HCM',[-2 -1 1 2],'FitRF', 'Display'); |
---|
| 80 | On 12-16-2003, [-2 -1 1 2] corresponded to HCM(6,3), HCM(6,4), HCM(7,1), HCM(7,3) |
---|
| 81 | Display plots information before applying the bump. |
---|
| 82 | |
---|
| 83 | 2. If the response matrix is not perfect, a clean (minimal rms orbit error) bump can be created |
---|
| 84 | by increasing the iteration (note: in regions of dispersion, all bumps will have leakage unless |
---|
| 85 | the RF frequency is included). |
---|
| 86 | [OCS, OCS0] = setorbitbump('BPMx',[6 6;7 1],[1;-1],'HCM',[-2 -1 1 2], 4, 'Display'); |
---|
| 87 | DeltaHCM = OCS.CM.Data - OCS0.CM.Data; |
---|
| 88 | BPMxDeviceList = OCS.BPM.DeviceList; |
---|
| 89 | HCMDeviceList = OCS.CM.DeviceList; |
---|
| 90 | |
---|
| 91 | |
---|
| 92 | See also <a href="setorbit.html" class="code" title="function [OCS, OCS0, V, S, ErrorFlag] = setorbit(varargin)">setorbit</a>, <a href="orbitcorrectionmethods.html" class="code" title="function [OCS, SmatNoWeights, S, U, V] = orbitcorrectionmethods(OCS, Smat, S, U, V)">orbitcorrectionmethods</a>, setorbitgui, <a href="rmdisp.html" class="code" title="function [DeltaRF, BPM, c, DispOrbit] = rmdisp(varargin)">rmdisp</a>, <a href="plotcm.html" class="code" title="function [DeltaRF, HCMEnergyChangeTotal, DeltaL] = plotcm(varargin)">plotcm</a>, setorbitbumpgui |
---|
| 93 | |
---|
| 94 | Written by Greg Portmann</pre></div> |
---|
| 95 | |
---|
| 96 | <!-- crossreference --> |
---|
| 97 | <h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 98 | This function calls: |
---|
| 99 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 100 | <li><a href="family2datastruct.html" class="code" title="function [DataStruct, ErrorFlag] = family2datastruct(varargin)">family2datastruct</a> FAMILY2DATASTRUCTURE - Returns a datastructure corresponding to a Family</li><li><a href="family2dev.html" class="code" title="function DeviceList = family2dev(Family, varargin);">family2dev</a> FAMILY2DEV - Return the device list for a family</li><li><a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a> GETAM - Gets monitor channels</li><li><a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a> GETFAMILYDATA - Gets data associated with the accelerator control</li><li><a href="getgolden.html" class="code" title="function Data = getgolden(varargin)">getgolden</a> GETGOLDEN - Returns the golden values for a family</li><li><a href="gethbpmfamily.html" class="code" title="function Family = gethbpmfamily">gethbpmfamily</a> GETHBPMFAMILY - Return the default horizontal BPM family</li><li><a href="gethcmfamily.html" class="code" title="function Family = gethcmfamily">gethcmfamily</a> GETHCMFAMILY - Returns the default horizontal corrector family</li><li><a href="getmode.html" class="code" title="function Mode = getmode(Family, Field)">getmode</a> GETMODE - Returns the present family mode ('Online', 'Simulator', 'Manual', 'Special', etc)</li><li><a href="getoffset.html" class="code" title="function Data = getoffset(varargin)">getoffset</a> GETOFFSET - Returns the offset values for a family</li><li><a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a> GETSPOS - Returns the longitudinal position in meters</li><li><a href="getunits.html" class="code" title="function [Units, UnitsString] = getunits(Family, Field)">getunits</a> GETUNITS - Return the present family units and units string</li><li><a href="getvbpmfamily.html" class="code" title="function Family = getvbpmfamily">getvbpmfamily</a> GETVBPMFAMILY - Return the default vertical BPM family</li><li><a href="getvcmfamily.html" class="code" title="function Family = getvcmfamily">getvcmfamily</a> GETVCMFAMILY - Returns the default vertical corrector family</li><li><a href="setorbit.html" class="code" title="function [OCS, OCS0, V, S, ErrorFlag] = setorbit(varargin)">setorbit</a> SETORBIT - Orbit correction function</li></ul> |
---|
| 101 | This function is called by: |
---|
| 102 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 103 | <li><a href="meascmhysteresis.html" class="code" title="function meascmhysteresis(BPMList, CMFamily, CMList, MaxChange, NSteps)">meascmhysteresis</a> MEASCMHYSTERESIS - Measure corrector magnet hysteresis</li><li><a href="quadcenterfit.html" class="code" title="function [XOffset, YOffset] = quadcenterfit(varargin)">quadcenterfit</a> BBASEARCH - Model search method to find a quarupole center</li><li><a href="quadcenterinit.html" class="code" title="function QMS = quadcenterinit(QuadFamily, QuadDev, QuadPlane)">quadcenterinit</a> QMS = quadcenterinit(Family, Device, QuadPlane)</li></ul> |
---|
| 104 | <!-- crossreference --> |
---|
| 105 | |
---|
| 106 | |
---|
| 107 | <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 108 | <div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [OCS, OCS0, V, S, ErrorFlag] = setorbitbump(varargin)</a> |
---|
| 109 | 0002 <span class="comment">%SETORBITBUMP - Local bump program (uses setorbit)</span> |
---|
| 110 | 0003 <span class="comment">% [OCS, OCS0, V, S] = setorbitbump(BPMFamily, BPMDeviceList, GoalOrbit, CMFamily, CMIncrementList, NIter, SVDIndex, BPMWeight)</span> |
---|
| 111 | 0004 <span class="comment">% [OCS, OCS0, V, S] = setorbitbump(OCS, GoalOrbit, CMIncrementList, NIter, SVDIndex, BPMWeight)</span> |
---|
| 112 | 0005 <span class="comment">%</span> |
---|
| 113 | 0006 <span class="comment">% INPUTS</span> |
---|
| 114 | 0007 <span class="comment">% 1. BPMFamily - BPM family</span> |
---|
| 115 | 0008 <span class="comment">% 2. BPMDeviceList - BPM device list</span> |
---|
| 116 | 0009 <span class="comment">% 3. GoalOrbit - Goal orbit position (Referenced to an absolute or incremental orbit</span> |
---|
| 117 | 0010 <span class="comment">% depending on 'Incremental' or 'Absolute' flag)</span> |
---|
| 118 | 0011 <span class="comment">% 4. CMFamily - Corrector magnet family</span> |
---|
| 119 | 0012 <span class="comment">% 5. CMIncrementList - Magnet list upstream and downstream of the BPMs</span> |
---|
| 120 | 0013 <span class="comment">% For instance, [-2 1 3] uses the 2nd correctors upstream and</span> |
---|
| 121 | 0014 <span class="comment">% 1st and 3rd corrector downstream from the BPMs in the bump.</span> |
---|
| 122 | 0015 <span class="comment">% or</span> |
---|
| 123 | 0016 <span class="comment">% CMIncrementList can be an ordinary device list used for specifying correctors.</span> |
---|
| 124 | 0017 <span class="comment">% 6. NIter - Number of iterations {Default: 3}</span> |
---|
| 125 | 0018 <span class="comment">% (set NIter to 0 if no BPMs are functioning, see setorbit)</span> |
---|
| 126 | 0019 <span class="comment">% 7. SVDIndex - vector, maximum number, 'All', or</span> |
---|
| 127 | 0020 <span class="comment">% base on a threshold of min/max singular value {Default: see setorbit}</span> |
---|
| 128 | 0021 <span class="comment">% 8. BPMWeight - Weight applied to the BPMs inside the bump (GoalOrbit). BPM weighting</span> |
---|
| 129 | 0022 <span class="comment">% should only be required when the corrector are in a region of dispersion</span> |
---|
| 130 | 0023 <span class="comment">% and the RF frequency is not used in the orbit correction.</span> |
---|
| 131 | 0024 <span class="comment">% {Default BPMWeight is to weight the leakage twice as much as the GoalOrbit}</span> |
---|
| 132 | 0025 <span class="comment">% 9. Optional flags (any flags used in setorbit can also be used here)</span> |
---|
| 133 | 0026 <span class="comment">% 'Incremental' {Default} or 'Absolute' - Type of orbit change</span> |
---|
| 134 | 0027 <span class="comment">% 'Display' - Display bump characteristics (no display is the default)</span> |
---|
| 135 | 0028 <span class="comment">% 'SetPV' or 'SetSP' - Make the setpoint change {Default}</span> |
---|
| 136 | 0029 <span class="comment">% 'NoSetPV' or 'NoSetSP' - Don't set the magnets, just return the OCS</span> |
---|
| 137 | 0030 <span class="comment">% 'FitRF' or 'SetRF' - Flag to include the RF frequency as part of orbit correction.</span> |
---|
| 138 | 0031 <span class="comment">% 'LeakageCorrection' - Option to try to clean up the leakage after the bump is applied.</span> |
---|
| 139 | 0032 <span class="comment">% Note: 1. NIter for leakage correction is the same for the bump</span> |
---|
| 140 | 0033 <span class="comment">% 2. It's just like setting the weight to zero for a second set of iterations.</span> |
---|
| 141 | 0034 <span class="comment">% 3. For the horizontal plane, usually one should include the RF frequency</span> |
---|
| 142 | 0035 <span class="comment">% along with leakage correction. Otherwise, the goal orbit may not be attained.</span> |
---|
| 143 | 0036 <span class="comment">% 'Tolerance', Tol - Quit correction if the std(error) < Tol {Tol: 0}</span> |
---|
| 144 | 0037 <span class="comment">% 'RampSteps', NSteps - Number of steps to ramp in the correction {NSteps: 1}</span> |
---|
| 145 | 0038 <span class="comment">% This is used to avoid a large transient between setpoint changes.</span> |
---|
| 146 | 0039 <span class="comment">%</span> |
---|
| 147 | 0040 <span class="comment">% OUTPUTS (same as setorbit)</span> |
---|
| 148 | 0041 <span class="comment">% 1. OCS - Orbit correction structure (OCS) usable by setorbit</span> |
---|
| 149 | 0042 <span class="comment">% 2. OCS0 - Starting OCS</span> |
---|
| 150 | 0043 <span class="comment">% 3. V - Corrector magnet singular vectors</span> |
---|
| 151 | 0044 <span class="comment">% 4. Svalue - Singular values (vector)</span> |
---|
| 152 | 0045 <span class="comment">%</span> |
---|
| 153 | 0046 <span class="comment">% NOTES</span> |
---|
| 154 | 0047 <span class="comment">% 1. setorbitbump creates an OCS structure and uses setorbit to actually change the orbit.</span> |
---|
| 155 | 0048 <span class="comment">%</span> |
---|
| 156 | 0049 <span class="comment">% EXAMPLES</span> |
---|
| 157 | 0050 <span class="comment">% 1. 4 magnet horizontal incremental bump at BPM(6,6) and BPM(7,1) including the RF frequency</span> |
---|
| 158 | 0051 <span class="comment">% OCS = setorbitbump('BPMx',[6 6;7 1],[1;-1],'HCM',[-2 -1 1 2],'FitRF', 'Display');</span> |
---|
| 159 | 0052 <span class="comment">% On 12-16-2003, [-2 -1 1 2] corresponded to HCM(6,3), HCM(6,4), HCM(7,1), HCM(7,3)</span> |
---|
| 160 | 0053 <span class="comment">% Display plots information before applying the bump.</span> |
---|
| 161 | 0054 <span class="comment">%</span> |
---|
| 162 | 0055 <span class="comment">% 2. If the response matrix is not perfect, a clean (minimal rms orbit error) bump can be created</span> |
---|
| 163 | 0056 <span class="comment">% by increasing the iteration (note: in regions of dispersion, all bumps will have leakage unless</span> |
---|
| 164 | 0057 <span class="comment">% the RF frequency is included).</span> |
---|
| 165 | 0058 <span class="comment">% [OCS, OCS0] = setorbitbump('BPMx',[6 6;7 1],[1;-1],'HCM',[-2 -1 1 2], 4, 'Display');</span> |
---|
| 166 | 0059 <span class="comment">% DeltaHCM = OCS.CM.Data - OCS0.CM.Data;</span> |
---|
| 167 | 0060 <span class="comment">% BPMxDeviceList = OCS.BPM.DeviceList;</span> |
---|
| 168 | 0061 <span class="comment">% HCMDeviceList = OCS.CM.DeviceList;</span> |
---|
| 169 | 0062 <span class="comment">%</span> |
---|
| 170 | 0063 <span class="comment">%</span> |
---|
| 171 | 0064 <span class="comment">% See also setorbit, orbitcorrectionmethods, setorbitgui, rmdisp, plotcm, setorbitbumpgui</span> |
---|
| 172 | 0065 <span class="comment">%</span> |
---|
| 173 | 0066 <span class="comment">% Written by Greg Portmann</span> |
---|
| 174 | 0067 |
---|
| 175 | 0068 |
---|
| 176 | 0069 <span class="comment">%%%%%%%%%%%%</span> |
---|
| 177 | 0070 <span class="comment">% Defaults %</span> |
---|
| 178 | 0071 <span class="comment">%%%%%%%%%%%%</span> |
---|
| 179 | 0072 OCS.BPM.FamilyName = []; |
---|
| 180 | 0073 OCS.BPM.DeviceList = []; |
---|
| 181 | 0074 GoalOrbit = []; |
---|
| 182 | 0075 OCS.CM.FamilyName = []; |
---|
| 183 | 0076 OCS.CM.DeviceList = []; |
---|
| 184 | 0077 CMIncrementList = []; |
---|
| 185 | 0078 OCS.NIter = 1; |
---|
| 186 | 0079 BPMWeight = []; |
---|
| 187 | 0080 OCS.Incremental = 1; |
---|
| 188 | 0081 DisplayFlag = 0; |
---|
| 189 | 0082 SetSPFlag = 1; |
---|
| 190 | 0083 NoLeakageFlag = 0; |
---|
| 191 | 0084 UnitsFlag = <span class="string">''</span>; |
---|
| 192 | 0085 ModeFlag = <span class="string">''</span>; |
---|
| 193 | 0086 |
---|
| 194 | 0087 L = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Circumference'</span>); |
---|
| 195 | 0088 |
---|
| 196 | 0089 |
---|
| 197 | 0090 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 198 | 0091 <span class="comment">% Input parsing and checking %</span> |
---|
| 199 | 0092 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 200 | 0093 Flags = {}; |
---|
| 201 | 0094 <span class="keyword">for</span> i = length(varargin):-1:1 |
---|
| 202 | 0095 <span class="keyword">if</span> isstruct(varargin{i}) |
---|
| 203 | 0096 <span class="comment">% Ignor structures</span> |
---|
| 204 | 0097 <span class="keyword">elseif</span> iscell(varargin{i}) |
---|
| 205 | 0098 <span class="comment">% Ignor cells</span> |
---|
| 206 | 0099 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Golden'</span>) |
---|
| 207 | 0100 <span class="comment">% Use the golden orbit</span> |
---|
| 208 | 0101 <span class="comment">%GoalOrbit = 'Golden';</span> |
---|
| 209 | 0102 <span class="comment">%varargin(i) = [];</span> |
---|
| 210 | 0103 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Offset'</span>) |
---|
| 211 | 0104 <span class="comment">% Use the offset orbit</span> |
---|
| 212 | 0105 <span class="comment">%GoalOrbit = 'Offset';</span> |
---|
| 213 | 0106 <span class="comment">%varargin(i) = [];</span> |
---|
| 214 | 0107 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Display'</span>) |
---|
| 215 | 0108 Flags = [Flags varargin(i)]; |
---|
| 216 | 0109 DisplayFlag = 1; |
---|
| 217 | 0110 varargin(i) = []; |
---|
| 218 | 0111 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'NoDisplay'</span>) |
---|
| 219 | 0112 Flags = [Flags varargin(i)]; |
---|
| 220 | 0113 DisplayFlag = 0; |
---|
| 221 | 0114 varargin(i) = []; |
---|
| 222 | 0115 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Inc'</span>) || strcmpi(varargin{i},<span class="string">'Incremental'</span>) |
---|
| 223 | 0116 OCS.Incremental = 1; |
---|
| 224 | 0117 varargin(i) = []; |
---|
| 225 | 0118 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Abs'</span>) || strcmpi(varargin{i},<span class="string">'Absolute'</span>) |
---|
| 226 | 0119 OCS.Incremental = 0; |
---|
| 227 | 0120 varargin(i) = []; |
---|
| 228 | 0121 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'simulator'</span>) || strcmpi(varargin{i},<span class="string">'model'</span>) |
---|
| 229 | 0122 ModeFlag = <span class="string">'SIMULATOR'</span>; |
---|
| 230 | 0123 varargin(i) = []; |
---|
| 231 | 0124 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Online'</span>) |
---|
| 232 | 0125 ModeFlag = <span class="string">'Online'</span>; |
---|
| 233 | 0126 varargin(i) = []; |
---|
| 234 | 0127 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Manual'</span>) |
---|
| 235 | 0128 ModeFlag = <span class="string">'Manual'</span>; |
---|
| 236 | 0129 varargin(i) = []; |
---|
| 237 | 0130 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'physics'</span>) |
---|
| 238 | 0131 UnitsFlag = <span class="string">'Physics'</span>; |
---|
| 239 | 0132 varargin(i) = []; |
---|
| 240 | 0133 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'hardware'</span>) |
---|
| 241 | 0134 UnitsFlag = <span class="string">'Hardware'</span>; |
---|
| 242 | 0135 varargin(i) = []; |
---|
| 243 | 0136 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'SetPV'</span>,<span class="string">'SetSP'</span>})) |
---|
| 244 | 0137 SetSPFlag = 1; |
---|
| 245 | 0138 varargin(i) = []; |
---|
| 246 | 0139 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'NoSetPV'</span>,<span class="string">'NoSetSP'</span>})) |
---|
| 247 | 0140 SetSPFlag = 0; |
---|
| 248 | 0141 Flags = [Flags varargin(i)]; |
---|
| 249 | 0142 varargin(i) = []; |
---|
| 250 | 0143 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'LeakageCorrection'</span>) |
---|
| 251 | 0144 NoLeakageFlag = 1; |
---|
| 252 | 0145 varargin(i) = []; |
---|
| 253 | 0146 |
---|
| 254 | 0147 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'archive'</span>) || strcmpi(varargin{i},<span class="string">'noarchive'</span>) |
---|
| 255 | 0148 <span class="comment">% Just remove</span> |
---|
| 256 | 0149 varargin(i) = []; |
---|
| 257 | 0150 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'struct'</span>) || strcmpi(varargin{i},<span class="string">'numeric'</span>) |
---|
| 258 | 0151 <span class="comment">% Just remove</span> |
---|
| 259 | 0152 varargin(i) = []; |
---|
| 260 | 0153 |
---|
| 261 | 0154 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'CorrectorGain'</span>) |
---|
| 262 | 0155 Flags = [Flags varargin(i) varargin(i+1)]; |
---|
| 263 | 0156 varargin(i+1) = []; |
---|
| 264 | 0157 varargin(i) = []; |
---|
| 265 | 0158 |
---|
| 266 | 0159 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Tolerance'</span>) |
---|
| 267 | 0160 Flags = [Flags varargin(i) varargin(i+1)]; |
---|
| 268 | 0161 varargin(i+1) = []; |
---|
| 269 | 0162 varargin(i) = []; |
---|
| 270 | 0163 |
---|
| 271 | 0164 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'RampSteps'</span>) |
---|
| 272 | 0165 Flags = [Flags varargin(i) varargin(i+1)]; |
---|
| 273 | 0166 varargin(i+1) = []; |
---|
| 274 | 0167 varargin(i) = []; |
---|
| 275 | 0168 |
---|
| 276 | 0169 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'ModelResp'</span>) |
---|
| 277 | 0170 Flags = [Flags varargin(i)]; |
---|
| 278 | 0171 varargin(i) = []; |
---|
| 279 | 0172 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'ModelDisp'</span>) |
---|
| 280 | 0173 Flags = [Flags varargin(i)]; |
---|
| 281 | 0174 varargin(i) = []; |
---|
| 282 | 0175 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'MeasDisp'</span>) |
---|
| 283 | 0176 Flags = [Flags varargin(i)]; |
---|
| 284 | 0177 varargin(i) = []; |
---|
| 285 | 0178 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'GoldenDisp'</span>) |
---|
| 286 | 0179 Flags = [Flags varargin(i)]; |
---|
| 287 | 0180 varargin(i) = []; |
---|
| 288 | 0181 |
---|
| 289 | 0182 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'FitRF'</span>,<span class="string">'SetRF'</span>,<span class="string">'RFCorrector'</span>})) |
---|
| 290 | 0183 Flags = [Flags varargin(i)]; |
---|
| 291 | 0184 varargin(i) = []; |
---|
| 292 | 0185 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'FitRFHCM0'</span>,<span class="string">'FitRFEnergy'</span>,<span class="string">'SetRFHCM0'</span>,<span class="string">'SetRFEnergy'</span>})) |
---|
| 293 | 0186 Flags = [Flags varargin(i)]; |
---|
| 294 | 0187 varargin(i) = []; |
---|
| 295 | 0188 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'FitRFDeltaHCM0'</span>,<span class="string">'SetRFDeltaHCM0'</span>})) |
---|
| 296 | 0189 Flags = [Flags varargin(i)]; |
---|
| 297 | 0190 varargin(i) = []; |
---|
| 298 | 0191 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'FitRFDispersionOrbit'</span>,<span class="string">'SetRFDispersionOrbit'</span>})) |
---|
| 299 | 0192 Flags = [Flags varargin(i)]; |
---|
| 300 | 0193 varargin(i) = []; |
---|
| 301 | 0194 <span class="keyword">elseif</span> any(strcmpi(varargin{i},{<span class="string">'FitRFDispersionHCM'</span>,<span class="string">'SetRFDispersionHCM'</span>})) |
---|
| 302 | 0195 Flags = [Flags varargin(i)]; |
---|
| 303 | 0196 varargin(i) = []; |
---|
| 304 | 0197 |
---|
| 305 | 0198 <span class="keyword">end</span> |
---|
| 306 | 0199 <span class="keyword">end</span> |
---|
| 307 | 0200 |
---|
| 308 | 0201 |
---|
| 309 | 0202 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 310 | 0203 <span class="keyword">if</span> isstruct(varargin{1}) |
---|
| 311 | 0204 OCS = varargin{1}; |
---|
| 312 | 0205 varargin(1) = []; |
---|
| 313 | 0206 <span class="keyword">if</span> length(varargin) < 2 |
---|
| 314 | 0207 error(<span class="string">'An OCS plus at least 2 inputs are required'</span>); |
---|
| 315 | 0208 <span class="keyword">end</span> |
---|
| 316 | 0209 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 317 | 0210 GoalOrbit = varargin{1}; |
---|
| 318 | 0211 varargin(1) = []; |
---|
| 319 | 0212 <span class="keyword">end</span> |
---|
| 320 | 0213 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 321 | 0214 CMIncrementList = varargin{1}; |
---|
| 322 | 0215 varargin(1) = []; |
---|
| 323 | 0216 <span class="keyword">end</span> |
---|
| 324 | 0217 <span class="keyword">else</span> |
---|
| 325 | 0218 OCS.BPM.FamilyName = varargin{1}; |
---|
| 326 | 0219 varargin(1) = []; |
---|
| 327 | 0220 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 328 | 0221 OCS.BPM.DeviceList = varargin{1}; |
---|
| 329 | 0222 varargin(1) = []; |
---|
| 330 | 0223 <span class="keyword">end</span> |
---|
| 331 | 0224 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 332 | 0225 GoalOrbit = varargin{1}; |
---|
| 333 | 0226 varargin(1) = []; |
---|
| 334 | 0227 <span class="keyword">end</span> |
---|
| 335 | 0228 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 336 | 0229 OCS.CM.FamilyName = varargin{1}; |
---|
| 337 | 0230 varargin(1) = []; |
---|
| 338 | 0231 <span class="keyword">end</span> |
---|
| 339 | 0232 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 340 | 0233 CMIncrementList = varargin{1}; |
---|
| 341 | 0234 varargin(1) = []; |
---|
| 342 | 0235 <span class="keyword">end</span> |
---|
| 343 | 0236 <span class="keyword">end</span> |
---|
| 344 | 0237 |
---|
| 345 | 0238 <span class="comment">% Get NIter</span> |
---|
| 346 | 0239 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 347 | 0240 <span class="keyword">if</span> isnumeric(varargin{1}) |
---|
| 348 | 0241 OCS.NIter = varargin{1}; |
---|
| 349 | 0242 varargin(1) = []; |
---|
| 350 | 0243 <span class="keyword">end</span> |
---|
| 351 | 0244 <span class="keyword">end</span> |
---|
| 352 | 0245 |
---|
| 353 | 0246 <span class="comment">% Get SVDIndex (can be a fractional number, vector, or 'All')</span> |
---|
| 354 | 0247 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 355 | 0248 <span class="comment">%if isnumeric(varargin{1})</span> |
---|
| 356 | 0249 OCS.SVDIndex = varargin{1}; |
---|
| 357 | 0250 varargin(1) = []; |
---|
| 358 | 0251 <span class="comment">%end</span> |
---|
| 359 | 0252 <span class="keyword">end</span> |
---|
| 360 | 0253 |
---|
| 361 | 0254 <span class="comment">% Get BPMWeight</span> |
---|
| 362 | 0255 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 363 | 0256 <span class="keyword">if</span> isnumeric(varargin{1}) |
---|
| 364 | 0257 BPMWeight = varargin{1}; |
---|
| 365 | 0258 varargin(1) = []; |
---|
| 366 | 0259 <span class="keyword">end</span> |
---|
| 367 | 0260 <span class="keyword">end</span> |
---|
| 368 | 0261 <span class="keyword">end</span> |
---|
| 369 | 0262 |
---|
| 370 | 0263 |
---|
| 371 | 0264 <span class="comment">% Pass the extra inputs on</span> |
---|
| 372 | 0265 <span class="keyword">if</span> length(varargin) >= 1 |
---|
| 373 | 0266 Flags = [Flags varargin]; |
---|
| 374 | 0267 <span class="keyword">end</span> |
---|
| 375 | 0268 |
---|
| 376 | 0269 |
---|
| 377 | 0270 <span class="comment">% Fill the empty fields</span> |
---|
| 378 | 0271 <span class="keyword">if</span> isempty(OCS.BPM.FamilyName) |
---|
| 379 | 0272 i = menu(<span class="string">'Select a plane'</span>,<span class="string">'Horizontal'</span>,<span class="string">'Vertical'</span>); |
---|
| 380 | 0273 <span class="keyword">if</span> i == 1 |
---|
| 381 | 0274 OCS.BPM.FamilyName = <a href="gethbpmfamily.html" class="code" title="function Family = gethbpmfamily">gethbpmfamily</a>; |
---|
| 382 | 0275 OCS.CM.FamilyName = <a href="gethcmfamily.html" class="code" title="function Family = gethcmfamily">gethcmfamily</a>; |
---|
| 383 | 0276 <span class="keyword">elseif</span> i == 2 |
---|
| 384 | 0277 OCS.BPM.FamilyName = <a href="getvbpmfamily.html" class="code" title="function Family = getvbpmfamily">getvbpmfamily</a>; |
---|
| 385 | 0278 OCS.CM.FamilyName = <a href="getvcmfamily.html" class="code" title="function Family = getvcmfamily">getvcmfamily</a>; |
---|
| 386 | 0279 <span class="keyword">end</span> |
---|
| 387 | 0280 <span class="keyword">end</span> |
---|
| 388 | 0281 <span class="keyword">if</span> isempty(OCS.BPM.DeviceList) |
---|
| 389 | 0282 OCS.BPM.DeviceList = editlist(<a href="family2dev.html" class="code" title="function DeviceList = family2dev(Family, varargin);">family2dev</a>(OCS.BPM.FamilyName), OCS.BPM.FamilyName, 0); |
---|
| 390 | 0283 <span class="comment">% Getting a list can messes up the order for bumps going from the last sector to the first</span> |
---|
| 391 | 0284 <span class="keyword">for</span> i = 2:length(OCS.BPM.DeviceList) |
---|
| 392 | 0285 <span class="keyword">if</span> <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.BPM.FamilyName,OCS.BPM.DeviceList(i,:))-<a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.BPM.FamilyName,OCS.BPM.DeviceList(i-1,:)) > L/2 |
---|
| 393 | 0286 <span class="keyword">break</span>; |
---|
| 394 | 0287 <span class="keyword">end</span> |
---|
| 395 | 0288 <span class="keyword">end</span> |
---|
| 396 | 0289 <span class="keyword">if</span> i < length(OCS.BPM.DeviceList) |
---|
| 397 | 0290 OCS.BPM.DeviceList = [OCS.BPM.DeviceList(i:<span class="keyword">end</span>,:); OCS.BPM.DeviceList(1:i-1,:)]; |
---|
| 398 | 0291 <span class="keyword">end</span> |
---|
| 399 | 0292 <span class="keyword">end</span> |
---|
| 400 | 0293 <span class="keyword">if</span> isempty(OCS.CM.FamilyName) |
---|
| 401 | 0294 <span class="keyword">if</span> strcmpi(OCS.BPM.FamilyName, <a href="gethbpmfamily.html" class="code" title="function Family = gethbpmfamily">gethbpmfamily</a>) |
---|
| 402 | 0295 OCS.CM.FamilyName = <a href="gethcmfamily.html" class="code" title="function Family = gethcmfamily">gethcmfamily</a>; |
---|
| 403 | 0296 <span class="keyword">else</span> |
---|
| 404 | 0297 OCS.CM.FamilyName = <a href="getvcmfamily.html" class="code" title="function Family = getvcmfamily">getvcmfamily</a>; |
---|
| 405 | 0298 <span class="keyword">end</span> |
---|
| 406 | 0299 <span class="keyword">end</span> |
---|
| 407 | 0300 |
---|
| 408 | 0301 <span class="keyword">if</span> isempty(UnitsFlag) |
---|
| 409 | 0302 UnitsFlag = <a href="getunits.html" class="code" title="function [Units, UnitsString] = getunits(Family, Field)">getunits</a>(OCS.BPM.FamilyName); |
---|
| 410 | 0303 <span class="keyword">end</span> |
---|
| 411 | 0304 |
---|
| 412 | 0305 <span class="keyword">if</span> isempty(ModeFlag) |
---|
| 413 | 0306 ModeFlag = <a href="getmode.html" class="code" title="function Mode = getmode(Family, Field)">getmode</a>(OCS.BPM.FamilyName); |
---|
| 414 | 0307 <span class="keyword">end</span> |
---|
| 415 | 0308 |
---|
| 416 | 0309 <span class="keyword">if</span> NoLeakageFlag && OCS.NIter==0 |
---|
| 417 | 0310 error(<span class="string">'Leakage cleanup flag cannot be used with zero iterations.'</span>); |
---|
| 418 | 0311 <span class="keyword">end</span> |
---|
| 419 | 0312 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 420 | 0313 <span class="comment">% End of input parsing %</span> |
---|
| 421 | 0314 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 422 | 0315 |
---|
| 423 | 0316 |
---|
| 424 | 0317 <span class="comment">% Possibly convert the GoalOrbit</span> |
---|
| 425 | 0318 <span class="keyword">if</span> any(isnan(GoalOrbit)) |
---|
| 426 | 0319 error(<span class="string">'Goal orbit has a NaN'</span>); |
---|
| 427 | 0320 <span class="keyword">end</span> |
---|
| 428 | 0321 <span class="keyword">if</span> isempty(GoalOrbit) |
---|
| 429 | 0322 <span class="comment">% Input the goal orbit</span> |
---|
| 430 | 0323 <span class="keyword">for</span> i = 1:size(OCS.BPM.DeviceList,1) |
---|
| 431 | 0324 labels{i} = sprintf(<span class="string">'Input the goal orbit at %s(%d,%d)'</span>,OCS.BPM.FamilyName, OCS.BPM.DeviceList(i,1), OCS.BPM.DeviceList(i,2)); |
---|
| 432 | 0325 StartingPoint{i} = <span class="string">'0'</span>; |
---|
| 433 | 0326 <span class="keyword">end</span> |
---|
| 434 | 0327 answer = inputdlg(labels,<span class="string">'Local Bump'</span>,1,StartingPoint); |
---|
| 435 | 0328 <span class="keyword">if</span> isempty(answer) |
---|
| 436 | 0329 fprintf(<span class="string">' Local bump cancelled'</span>); |
---|
| 437 | 0330 <span class="keyword">return</span> |
---|
| 438 | 0331 <span class="keyword">end</span> |
---|
| 439 | 0332 <span class="keyword">for</span> i = 1:size(OCS.BPM.DeviceList,1) |
---|
| 440 | 0333 GoalOrbit(i,1) = str2num(answer{i}); |
---|
| 441 | 0334 <span class="keyword">end</span> |
---|
| 442 | 0335 <span class="keyword">end</span> |
---|
| 443 | 0336 <span class="keyword">if</span> ischar(GoalOrbit) |
---|
| 444 | 0337 <span class="keyword">if</span> strcmpi(GoalOrbit, <span class="string">'Golden'</span>) |
---|
| 445 | 0338 GoalOrbit = <a href="getgolden.html" class="code" title="function Data = getgolden(varargin)">getgolden</a>(OCS.BPM.FamilyName, OCS.BPM.DeviceList); |
---|
| 446 | 0339 <span class="keyword">elseif</span> strcmpi(GoalOrbit, <span class="string">'Offset'</span>) |
---|
| 447 | 0340 GoalOrbit = <a href="getoffset.html" class="code" title="function Data = getoffset(varargin)">getoffset</a>(OCS.BPM.FamilyName, OCS.BPM.DeviceList); |
---|
| 448 | 0341 <span class="keyword">else</span> |
---|
| 449 | 0342 error(<span class="string">'Goal unknown'</span>); |
---|
| 450 | 0343 <span class="keyword">end</span> |
---|
| 451 | 0344 <span class="keyword">end</span> |
---|
| 452 | 0345 GoalOrbit = GoalOrbit(:); |
---|
| 453 | 0346 |
---|
| 454 | 0347 <span class="comment">% Check the length</span> |
---|
| 455 | 0348 <span class="keyword">if</span> length(GoalOrbit) ~= size(OCS.BPM.DeviceList,1) |
---|
| 456 | 0349 error(<span class="string">'Length of the GoalOrbit must equal the number of devices in BPMList'</span>); |
---|
| 457 | 0350 <span class="keyword">end</span> |
---|
| 458 | 0351 |
---|
| 459 | 0352 |
---|
| 460 | 0353 <span class="comment">% CMIncrementList must be a vector, no zeros, no repeats</span> |
---|
| 461 | 0354 CheckCM = 0; |
---|
| 462 | 0355 <span class="keyword">if</span> isempty(CMIncrementList) |
---|
| 463 | 0356 CheckCM = 1; |
---|
| 464 | 0357 Ncm = size(OCS.BPM.DeviceList,1)+2; |
---|
| 465 | 0358 <span class="keyword">if</span> rem(Ncm,2) == 0 |
---|
| 466 | 0359 NcmDiv2 = Ncm / 2; |
---|
| 467 | 0360 CMIncrementList = [-(Ncm/2):-1 1:(Ncm/2)]; |
---|
| 468 | 0361 <span class="keyword">else</span> |
---|
| 469 | 0362 NcmDiv2 = Ncm / 2; |
---|
| 470 | 0363 CMIncrementList = [-floor(Ncm/2):-1 1:(floor(Ncm/2)+1)]; |
---|
| 471 | 0364 <span class="keyword">end</span> |
---|
| 472 | 0365 <span class="keyword">end</span> |
---|
| 473 | 0366 |
---|
| 474 | 0367 <span class="comment">% Check for a device list input</span> |
---|
| 475 | 0368 <span class="keyword">if</span> ~(size(CMIncrementList,2) == 2 && size(CMIncrementList,1) > 1) |
---|
| 476 | 0369 CMIncrementList = CMIncrementList(:); |
---|
| 477 | 0370 CMIncrementList = sort(CMIncrementList); |
---|
| 478 | 0371 CMIncrementList(find(CMIncrementList==0)) = []; |
---|
| 479 | 0372 CMIncrementList(find(diff(CMIncrementList)==0)) = []; |
---|
| 480 | 0373 <span class="keyword">end</span> |
---|
| 481 | 0374 |
---|
| 482 | 0375 <span class="comment">% Get BPM positions</span> |
---|
| 483 | 0376 BPMListTotal = <a href="family2dev.html" class="code" title="function DeviceList = family2dev(Family, varargin);">family2dev</a>(OCS.BPM.FamilyName, 1); |
---|
| 484 | 0377 BPMsposTotal = <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.BPM.FamilyName, BPMListTotal); |
---|
| 485 | 0378 BPMspos = <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.BPM.FamilyName, OCS.BPM.DeviceList); |
---|
| 486 | 0379 |
---|
| 487 | 0380 |
---|
| 488 | 0381 <span class="comment">% Stack 3 rings so you you don't have to worry about the L to 0 transition</span> |
---|
| 489 | 0382 CMListTotal = <a href="family2dev.html" class="code" title="function DeviceList = family2dev(Family, varargin);">family2dev</a>(OCS.CM.FamilyName, 1); |
---|
| 490 | 0383 CMsposTotal = <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.CM.FamilyName, CMListTotal); |
---|
| 491 | 0384 CMListTotal = [CMListTotal; CMListTotal; CMListTotal]; |
---|
| 492 | 0385 CMsposTotal = [CMsposTotal-L; CMsposTotal; CMsposTotal+L]; |
---|
| 493 | 0386 |
---|
| 494 | 0387 |
---|
| 495 | 0388 <span class="comment">% Find the correctors</span> |
---|
| 496 | 0389 <span class="comment">% Check for a device list input</span> |
---|
| 497 | 0390 <span class="keyword">if</span> size(CMIncrementList,2) == 2 && size(CMIncrementList,1) > 1 |
---|
| 498 | 0391 <span class="comment">% DeviceList</span> |
---|
| 499 | 0392 OCS.CM.DeviceList = CMIncrementList; |
---|
| 500 | 0393 <span class="keyword">else</span> |
---|
| 501 | 0394 <span class="keyword">for</span> i = 1:length(CMIncrementList) |
---|
| 502 | 0395 <span class="keyword">if</span> CMIncrementList(i) <= 0 |
---|
| 503 | 0396 j = find(CMsposTotal <= BPMspos(1)); |
---|
| 504 | 0397 OCS.CM.DeviceList(i,:) = CMListTotal(j(end)+CMIncrementList(i)+1,:); |
---|
| 505 | 0398 <span class="keyword">else</span> |
---|
| 506 | 0399 j = find(CMsposTotal >= BPMspos(end)); |
---|
| 507 | 0400 OCS.CM.DeviceList(i,:) = CMListTotal(j(1)+CMIncrementList(i)-1,:); |
---|
| 508 | 0401 <span class="keyword">end</span> |
---|
| 509 | 0402 <span class="keyword">end</span> |
---|
| 510 | 0403 <span class="keyword">end</span> |
---|
| 511 | 0404 |
---|
| 512 | 0405 <span class="keyword">if</span> CheckCM |
---|
| 513 | 0406 DevTotal = <a href="family2dev.html" class="code" title="function DeviceList = family2dev(Family, varargin);">family2dev</a>(OCS.CM.FamilyName); |
---|
| 514 | 0407 checkvec = zeros(size(DevTotal,1),1); |
---|
| 515 | 0408 checkvec(findrowindex(OCS.CM.DeviceList, DevTotal)) = 1; |
---|
| 516 | 0409 OCS.CM.DeviceList = editlist(DevTotal, OCS.CM.FamilyName, checkvec); |
---|
| 517 | 0410 |
---|
| 518 | 0411 <span class="comment">% Getting a list can messes up the order for bumps going from the last sector to the first</span> |
---|
| 519 | 0412 <span class="keyword">for</span> i = 2:length(OCS.CM.DeviceList) |
---|
| 520 | 0413 <span class="keyword">if</span> <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.CM.FamilyName,OCS.CM.DeviceList(i,:)) - <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.CM.FamilyName,OCS.CM.DeviceList(i-1,:)) > L/2 |
---|
| 521 | 0414 <span class="keyword">break</span>; |
---|
| 522 | 0415 <span class="keyword">end</span> |
---|
| 523 | 0416 <span class="keyword">end</span> |
---|
| 524 | 0417 <span class="keyword">if</span> i < length(OCS.CM.DeviceList) |
---|
| 525 | 0418 OCS.CM.DeviceList = [OCS.CM.DeviceList(i:<span class="keyword">end</span>,:); OCS.CM.DeviceList(1:i-1,:)]; |
---|
| 526 | 0419 <span class="keyword">end</span> |
---|
| 527 | 0420 <span class="keyword">end</span> |
---|
| 528 | 0421 |
---|
| 529 | 0422 |
---|
| 530 | 0423 <span class="comment">% Find all BPMs outside the bump (leakage control BPMs)</span> |
---|
| 531 | 0424 CMspos = <a href="getspos.html" class="code" title="function S = getspos(Family, DeviceList)">getspos</a>(OCS.CM.FamilyName, OCS.CM.DeviceList); |
---|
| 532 | 0425 <span class="keyword">if</span> CMspos(1) > CMspos(end) |
---|
| 533 | 0426 j1 = intersect(find(BPMsposTotal < CMspos(1)), find(BPMsposTotal > CMspos(end))); |
---|
| 534 | 0427 j2 = []; |
---|
| 535 | 0428 <span class="keyword">else</span> |
---|
| 536 | 0429 j1 = find(BPMsposTotal < CMspos(1)); |
---|
| 537 | 0430 j2 = find(BPMsposTotal > CMspos(end)); |
---|
| 538 | 0431 <span class="keyword">end</span> |
---|
| 539 | 0432 <span class="keyword">if</span> isempty(j1) && isempty(j2) |
---|
| 540 | 0433 error(<span class="string">'Cound not find any leakage control BPMs'</span>); |
---|
| 541 | 0434 <span class="keyword">end</span> |
---|
| 542 | 0435 OCS.BPM.DeviceList = [BPMListTotal(j1,:); OCS.BPM.DeviceList; BPMListTotal(j2,:);]; |
---|
| 543 | 0436 BPMDeviceListLeakage = [BPMListTotal(j1,:); BPMListTotal(j2,:);]; |
---|
| 544 | 0437 |
---|
| 545 | 0438 |
---|
| 546 | 0439 <span class="keyword">if</span> OCS.Incremental |
---|
| 547 | 0440 <span class="comment">% Incremental</span> |
---|
| 548 | 0441 OCS.GoalOrbit = [zeros(length(j1),1); GoalOrbit; zeros(length(j2),1)]; |
---|
| 549 | 0442 |
---|
| 550 | 0443 <span class="keyword">if</span> OCS.NIter == 0 |
---|
| 551 | 0444 LeakageGoalOrbit = zeros(size(BPMDeviceListLeakage,1),1); |
---|
| 552 | 0445 <span class="keyword">else</span> |
---|
| 553 | 0446 <span class="keyword">if</span> isempty(j1) |
---|
| 554 | 0447 LeakageOrbit1 = []; |
---|
| 555 | 0448 <span class="keyword">else</span> |
---|
| 556 | 0449 LeakageOrbit1 = <a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a>(OCS.BPM.FamilyName,BPMListTotal(j1,:), UnitsFlag, ModeFlag); |
---|
| 557 | 0450 <span class="keyword">end</span> |
---|
| 558 | 0451 <span class="keyword">if</span> isempty(j2) |
---|
| 559 | 0452 LeakageOrbit2 = []; |
---|
| 560 | 0453 <span class="keyword">else</span> |
---|
| 561 | 0454 LeakageOrbit2 = <a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a>(OCS.BPM.FamilyName,BPMListTotal(j2,:), UnitsFlag, ModeFlag); |
---|
| 562 | 0455 <span class="keyword">end</span> |
---|
| 563 | 0456 LeakageGoalOrbit = [LeakageOrbit1; LeakageOrbit2]; |
---|
| 564 | 0457 <span class="keyword">end</span> |
---|
| 565 | 0458 <span class="keyword">elseif</span> ~OCS.Incremental |
---|
| 566 | 0459 <span class="comment">% Absolute</span> |
---|
| 567 | 0460 <span class="keyword">if</span> isempty(j1) |
---|
| 568 | 0461 LeakageOrbit1 = []; |
---|
| 569 | 0462 <span class="keyword">else</span> |
---|
| 570 | 0463 LeakageOrbit1 = <a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a>(OCS.BPM.FamilyName,BPMListTotal(j1,:), UnitsFlag, ModeFlag); |
---|
| 571 | 0464 <span class="keyword">end</span> |
---|
| 572 | 0465 <span class="keyword">if</span> isempty(j2) |
---|
| 573 | 0466 LeakageOrbit2 = []; |
---|
| 574 | 0467 <span class="keyword">else</span> |
---|
| 575 | 0468 LeakageOrbit2 = <a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a>(OCS.BPM.FamilyName,BPMListTotal(j2,:), UnitsFlag, ModeFlag); |
---|
| 576 | 0469 <span class="keyword">end</span> |
---|
| 577 | 0470 OCS.GoalOrbit = [LeakageOrbit1; GoalOrbit; LeakageOrbit2]; |
---|
| 578 | 0471 LeakageGoalOrbit = [LeakageOrbit1; LeakageOrbit2]; |
---|
| 579 | 0472 <span class="keyword">else</span> |
---|
| 580 | 0473 error(<span class="string">'Absolute or incremental input unknown.'</span>); |
---|
| 581 | 0474 <span class="keyword">end</span> |
---|
| 582 | 0475 |
---|
| 583 | 0476 |
---|
| 584 | 0477 <span class="comment">% Add a weight</span> |
---|
| 585 | 0478 <span class="keyword">if</span> isempty(BPMWeight) |
---|
| 586 | 0479 <span class="comment">% Default BPMWeight is to weight the leakage twice as much as the bump goal</span> |
---|
| 587 | 0480 BPMWeight = .5 * (length(j1) + length(j2)) / length(GoalOrbit); |
---|
| 588 | 0481 OCS.BPMWeight = [ones(length(j1),1); BPMWeight .* ones(length(GoalOrbit),1); ones(length(j2),1)]; |
---|
| 589 | 0482 <span class="keyword">else</span> |
---|
| 590 | 0483 OCS.BPMWeight = [ones(length(j1),1); BPMWeight .* ones(length(GoalOrbit),1); ones(length(j2),1)]; |
---|
| 591 | 0484 <span class="keyword">end</span> |
---|
| 592 | 0485 |
---|
| 593 | 0486 |
---|
| 594 | 0487 <span class="comment">% Just to fill the data structure fields</span> |
---|
| 595 | 0488 OCS.BPM = <a href="family2datastruct.html" class="code" title="function [DataStruct, ErrorFlag] = family2datastruct(varargin)">family2datastruct</a>(OCS.BPM.FamilyName, <span class="string">'Monitor'</span>, OCS.BPM.DeviceList, UnitsFlag, ModeFlag); |
---|
| 596 | 0489 OCS.CM = <a href="family2datastruct.html" class="code" title="function [DataStruct, ErrorFlag] = family2datastruct(varargin)">family2datastruct</a>(OCS.CM.FamilyName, <span class="string">'Setpoint'</span>, OCS.CM.DeviceList, UnitsFlag, ModeFlag); |
---|
| 597 | 0490 |
---|
| 598 | 0491 |
---|
| 599 | 0492 <span class="comment">% Bumps stats</span> |
---|
| 600 | 0493 <span class="keyword">if</span> DisplayFlag |
---|
| 601 | 0494 <span class="comment">% Corrector strength (meters/radian) and (mm/amp)</span> |
---|
| 602 | 0495 <span class="comment">% Warn on when to add RF frequency (generated dispersion is greater</span> |
---|
| 603 | 0496 <span class="comment">% than mm per amp or % of mm/amp bump)</span> |
---|
| 604 | 0497 <span class="keyword">end</span> |
---|
| 605 | 0498 |
---|
| 606 | 0499 |
---|
| 607 | 0500 <span class="comment">% Set the bump</span> |
---|
| 608 | 0501 [OCS, OCS0, V, S, ErrorFlag] = <a href="setorbit.html" class="code" title="function [OCS, OCS0, V, S, ErrorFlag] = setorbit(varargin)">setorbit</a>(OCS, UnitsFlag, ModeFlag, Flags{:}); |
---|
| 609 | 0502 |
---|
| 610 | 0503 |
---|
| 611 | 0504 <span class="comment">% Clean up the leakage</span> |
---|
| 612 | 0505 <span class="keyword">if</span> NoLeakageFlag |
---|
| 613 | 0506 <span class="keyword">if</span> DisplayFlag |
---|
| 614 | 0507 fprintf(<span class="string">' Correcting the leakage\n'</span>); |
---|
| 615 | 0508 <span class="keyword">end</span> |
---|
| 616 | 0509 OCS1 = OCS; |
---|
| 617 | 0510 OCS1.GoalOrbit = LeakageGoalOrbit; |
---|
| 618 | 0511 OCS1.BPM = <a href="family2datastruct.html" class="code" title="function [DataStruct, ErrorFlag] = family2datastruct(varargin)">family2datastruct</a>(OCS.BPM.FamilyName, <span class="string">'Monitor'</span>, BPMDeviceListLeakage); |
---|
| 619 | 0512 OCS1.BPMWeight = []; |
---|
| 620 | 0513 OCS1.Incremental = 0; |
---|
| 621 | 0514 OCS1 = <a href="setorbit.html" class="code" title="function [OCS, OCS0, V, S, ErrorFlag] = setorbit(varargin)">setorbit</a>(OCS1, UnitsFlag, ModeFlag); |
---|
| 622 | 0515 <span class="keyword">end</span> |
---|
| 623 | 0516</pre></div> |
---|
| 624 | <hr><address>Generated on Mon 21-May-2007 15:29:18 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address> |
---|
| 625 | </body> |
---|
| 626 | </html> |
---|