Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

measchro.html @ 4

Last change on this file since 4 was 4, checked in by zhangj, 10 years ago
Initial import--MML version from SOLEIL@2013
File size: 50.3 KB

Line
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 measchro</title>
6	<meta name="keywords" content="measchro">
7	<meta name="description" content="MEASCHRO - measures the chromaticity function emperically">
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> > measchro.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>measchro
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>MEASCHRO - measures the chromaticity function emperically</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 [Chromaticity, FileName] = measchro(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">MEASCHRO - measures the chromaticity function emperically
31	Chrom = measchro(DeltaRF, WaitFlag);
32	ChromHardware = measchro(DeltaRF, WaitFlag, 'Hardware');
33	ChromPhysics = measchro(DeltaRF, WaitFlag, 'Physics');
34	ChromStruct = measchro(DeltaRF, WaitFlag, 'Struct');
35
36	INPUTS
37	1. DeltaRF - Vector of master oscillator values to scan over
38	{Default: [-.4% -.2% 0 .2% .4%] energy change}
39	2. WaitFlag >= 0, WaitFlag seconds before measuring the tune (sec)
40	= -1, wait until the magnets are done ramping
41	= -3, wait until the magnets are done ramping + a delay of 2.2*getfamilydata('TuneDelay') {default}
42	= -4, wait until keyboard input
43	= -5, input the tune measurement manually by keyboard input
44	4. 'Hardware' - Returns chromaticity in hardware units (typically, Tune/MHz or Tune/MHz)
45	'Physics' - Returns chromaticity in physics units (Tune/(dp/p)) {Default}
46	5. 'Struct' - Will return a two element dispersion data structure array {Default, unless Mode='Model'}
47	'Numeric' - Will return vector outputs
48	6. Optional override of the mode:
49	'Online' - Set/Get data online
50	'Simulator' - Set/Get data on the simulated accelerator (ie, same commands as 'Online')
51	'Model' - Get the model chromaticity directly from the model (uses modelchro, DeltaRF is ignored)
52	'Manual' - Set/Get data manually
53	7. 'Archive' - Save a chromaticity data structure to \<Directory.ChroData>\Chromaticity\
54	with filename <ChroArchiveFile><Date><Time>.mat {Default, unless Mode='Model'}
55	To change the filename, included the filename after the 'Archive', '' to browse
56	'NoArchive' - No file archive {Default}
57	8. 'Display' - Prints status information to the command window {Default, unless Mode='Model'}
58	'NoDisplay' - Nothing is printed to the command window
59
60
61	OUTPUT
62	\| Horizontal Chromaticity \|
63	ChromHardware = \| \| [Delta Tune / Delta Frequency]
64	\| Vertical Chromaticity \| (Hardware Units)
65
66
67	\| Horizontal Chromaticity \|
68	ChromPhysics = \| \| [Delta Tune / Delta Energy]
69	\| Vertical Chromaticity \| (Physics Units)
70
71	When computing physics units the momentum compaction factor is required. The default MCF is
72	found using getmcf. To override the default enter the new value after the 'Physics' input.
73	For example, ChromPhysics = measchro(DeltaRF, WaitFlag, 'Physics', .0011);
74
75	Tune vs RF frequency or momentum are plotted to the screen
76
77	Fields for structure outputs:
78	Data: [2x1] Chromaticity vector
79	FamilyName: 'Chromaticity'
80	Monitor: Tune structure
81	Actuator: RF frequency structure
82	DeltaActuator: Vector of frequency shifts in Hz
83	TimeStamp: Timestamp
84	DataDescriptor: 'Chromacity'
85	CreatedBy: 'measchro'
86	MCF: Momentum compaction factor/linear
87	RF: Vector of frequency settings in Hz
88	X: Reference orbit
89	Y: Reference orbit
90	Tune0: Initial tune
91	Tune: Tune change with RF frequency, 2 row vectors
92	dp: Vector of normalized momentum shifts
93	PolyFit: Polynomial fit of chromaticity in terms of rf shift or momentum
94
95	NOTE
96	1. 'Hardware', 'Physics', 'Eta', 'Archive', 'Numeric', and 'Struct' are not case sensitive
97	2. 'Zeta' can be used instead of 'Physics'
98	3. All inputs are optional
99	4. One reason FamilyName is added to the output structure so that getdata can be
100	used to locate archived dispersion measurements.
101	5. Units for DeltaRF depend on the 'Physics' or 'Hardware' flags
102	6. Beware of what units you are working in. The default units for chromaticity
103	are physics units. This is an exception to the normal middle layer convention.
104	Hardware units for "chromaticity" is in tune change per change in RF frequency.
105	Since this is an unusual unit to work with, the default units for chromaticity
106	is physics units. Note that goal chromaticity is also stored in physics units.
107	plotchro can switch between 'Hardware' and 'Physics' after the measurement is taken.
108	As an example of the difference between the units, at Spear3 1 unit of chromaticity
109	in physics units corresponds to roughly -1.8 units in hardware units.
110
111	See also <a href="plotchro.html" class="code" title="function [c, FileName] = plotchro(varargin)">plotchro</a>, <a href="measdisp.html" class="code" title="function [Dx, Dy, FileName] = measdisp(varargin);">measdisp</a>
112
113	Written by Greg Portmann and Jeff Corbett</pre></div>
114
115	<!-- crossreference -->
116	<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
117	This function calls:
118	<ul style="list-style-image:url(../matlabicon.gif)">
119	<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="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a> GETAM - Gets monitor channels</li><li><a href="getenergy.html" class="code" title="function [Energy, HCMEnergy] = getenergy(varargin)">getenergy</a> GETENERGY - Returns the beam energy base on the bend magnet</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="gethbpmfamily.html" class="code" title="function Family = gethbpmfamily">gethbpmfamily</a> GETHBPMFAMILY - Return the default horizontal BPM family</li><li><a href="getmcf.html" class="code" title="function Alpha = getmcf(ModelString)">getmcf</a> GETMCF - Returns the momentum compaction factor (MCF) stored in the AD or the model</li><li><a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a> GETRF - Gets the RF frequency</li><li><a href="gettune.html" class="code" title="function [Tune, tout, DataTime, ErrorFlag] = gettune(varargin)">gettune</a> GETTUNE - Returns the betatron tunes</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="getx.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getx(varargin)">getx</a> GETX - Returns the horizontal orbit</li><li><a href="gety.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = gety(varargin)">gety</a> GETY - Returns the vertical orbit</li><li><a href="hw2physics.html" class="code" title="function S = hw2physics(Family, Field, value, DeviceList, Energy)">hw2physics</a> HW2PHYSICS - Converts from 'Hardware' units to 'Physics' units</li><li><a href="isfamily.html" class="code" title="function [FamilyFlag, AO] = isfamily(Family, Field)">isfamily</a> ISFAMILY - True for family names</li><li><a href="plotchro.html" class="code" title="function [c, FileName] = plotchro(varargin)">plotchro</a> PLOTCHRO - Plot the chromaticity function</li><li><a href="setrf.html" class="code" title="function setrf(RF, varargin)">setrf</a> SETRF - Sets the RF frequency</li><li><a href="setsp.html" class="code" title="function ErrorFlag = setsp(Family, varargin)">setsp</a> SETSP - Makes an absolute setpoint change to the 'Setpoint' field</li></ul>
120	This function is called by:
121	<ul style="list-style-image:url(../matlabicon.gif)">
122	<li><a href="getchro.html" class="code" title="function [Data, FileName] = getchro(varargin)">getchro</a> GETCHRO - Return the chromaticity function (from file)</li><li><a href="measchroresp.html" class="code" title="function [Rmat, OutputFileName] = measchroresp(varargin)">measchroresp</a> MEASCHRORESP - measures the response from sextupoles to chromaticity</li><li><a href="setchro.html" class="code" title="function [DelSext, ActuatorFamily] = setchro(varargin)">setchro</a> SETCHRO - Measures then sets the chromaticity</li></ul>
123	<!-- crossreference -->
124
125
126	<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
127	<div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [Chromaticity, FileName] = measchro(varargin)</a>
128	0002 <span class="comment">%MEASCHRO - measures the chromaticity function emperically</span>
129	0003 <span class="comment">% Chrom = measchro(DeltaRF, WaitFlag);</span>
130	0004 <span class="comment">% ChromHardware = measchro(DeltaRF, WaitFlag, 'Hardware');</span>
131	0005 <span class="comment">% ChromPhysics = measchro(DeltaRF, WaitFlag, 'Physics');</span>
132	0006 <span class="comment">% ChromStruct = measchro(DeltaRF, WaitFlag, 'Struct');</span>
133	0007 <span class="comment">%</span>
134	0008 <span class="comment">% INPUTS</span>
135	0009 <span class="comment">% 1. DeltaRF - Vector of master oscillator values to scan over</span>
136	0010 <span class="comment">% {Default: [-.4% -.2% 0 .2% .4%] energy change}</span>
137	0011 <span class="comment">% 2. WaitFlag >= 0, WaitFlag seconds before measuring the tune (sec)</span>
138	0012 <span class="comment">% = -1, wait until the magnets are done ramping</span>
139	0013 <span class="comment">% = -3, wait until the magnets are done ramping + a delay of 2.2*getfamilydata('TuneDelay') {default}</span>
140	0014 <span class="comment">% = -4, wait until keyboard input</span>
141	0015 <span class="comment">% = -5, input the tune measurement manually by keyboard input</span>
142	0016 <span class="comment">% 4. 'Hardware' - Returns chromaticity in hardware units (typically, Tune/MHz or Tune/MHz)</span>
143	0017 <span class="comment">% 'Physics' - Returns chromaticity in physics units (Tune/(dp/p)) {Default}</span>
144	0018 <span class="comment">% 5. 'Struct' - Will return a two element dispersion data structure array {Default, unless Mode='Model'}</span>
145	0019 <span class="comment">% 'Numeric' - Will return vector outputs</span>
146	0020 <span class="comment">% 6. Optional override of the mode:</span>
147	0021 <span class="comment">% 'Online' - Set/Get data online</span>
148	0022 <span class="comment">% 'Simulator' - Set/Get data on the simulated accelerator (ie, same commands as 'Online')</span>
149	0023 <span class="comment">% 'Model' - Get the model chromaticity directly from the model (uses modelchro, DeltaRF is ignored)</span>
150	0024 <span class="comment">% 'Manual' - Set/Get data manually</span>
151	0025 <span class="comment">% 7. 'Archive' - Save a chromaticity data structure to \<Directory.ChroData>\Chromaticity\</span>
152	0026 <span class="comment">% with filename <ChroArchiveFile><Date><Time>.mat {Default, unless Mode='Model'}</span>
153	0027 <span class="comment">% To change the filename, included the filename after the 'Archive', '' to browse</span>
154	0028 <span class="comment">% 'NoArchive' - No file archive {Default}</span>
155	0029 <span class="comment">% 8. 'Display' - Prints status information to the command window {Default, unless Mode='Model'}</span>
156	0030 <span class="comment">% 'NoDisplay' - Nothing is printed to the command window</span>
157	0031 <span class="comment">%</span>
158	0032 <span class="comment">%</span>
159	0033 <span class="comment">% OUTPUT</span>
160	0034 <span class="comment">% \| Horizontal Chromaticity \|</span>
161	0035 <span class="comment">% ChromHardware = \| \| [Delta Tune / Delta Frequency]</span>
162	0036 <span class="comment">% \| Vertical Chromaticity \| (Hardware Units)</span>
163	0037 <span class="comment">%</span>
164	0038 <span class="comment">%</span>
165	0039 <span class="comment">% \| Horizontal Chromaticity \|</span>
166	0040 <span class="comment">% ChromPhysics = \| \| [Delta Tune / Delta Energy]</span>
167	0041 <span class="comment">% \| Vertical Chromaticity \| (Physics Units)</span>
168	0042 <span class="comment">%</span>
169	0043 <span class="comment">% When computing physics units the momentum compaction factor is required. The default MCF is</span>
170	0044 <span class="comment">% found using getmcf. To override the default enter the new value after the 'Physics' input.</span>
171	0045 <span class="comment">% For example, ChromPhysics = measchro(DeltaRF, WaitFlag, 'Physics', .0011);</span>
172	0046 <span class="comment">%</span>
173	0047 <span class="comment">% Tune vs RF frequency or momentum are plotted to the screen</span>
174	0048 <span class="comment">%</span>
175	0049 <span class="comment">% Fields for structure outputs:</span>
176	0050 <span class="comment">% Data: [2x1] Chromaticity vector</span>
177	0051 <span class="comment">% FamilyName: 'Chromaticity'</span>
178	0052 <span class="comment">% Monitor: Tune structure</span>
179	0053 <span class="comment">% Actuator: RF frequency structure</span>
180	0054 <span class="comment">% DeltaActuator: Vector of frequency shifts in Hz</span>
181	0055 <span class="comment">% TimeStamp: Timestamp</span>
182	0056 <span class="comment">% DataDescriptor: 'Chromacity'</span>
183	0057 <span class="comment">% CreatedBy: 'measchro'</span>
184	0058 <span class="comment">% MCF: Momentum compaction factor/linear</span>
185	0059 <span class="comment">% RF: Vector of frequency settings in Hz</span>
186	0060 <span class="comment">% X: Reference orbit</span>
187	0061 <span class="comment">% Y: Reference orbit</span>
188	0062 <span class="comment">% Tune0: Initial tune</span>
189	0063 <span class="comment">% Tune: Tune change with RF frequency, 2 row vectors</span>
190	0064 <span class="comment">% dp: Vector of normalized momentum shifts</span>
191	0065 <span class="comment">% PolyFit: Polynomial fit of chromaticity in terms of rf shift or momentum</span>
192	0066 <span class="comment">%</span>
193	0067 <span class="comment">% NOTE</span>
194	0068 <span class="comment">% 1. 'Hardware', 'Physics', 'Eta', 'Archive', 'Numeric', and 'Struct' are not case sensitive</span>
195	0069 <span class="comment">% 2. 'Zeta' can be used instead of 'Physics'</span>
196	0070 <span class="comment">% 3. All inputs are optional</span>
197	0071 <span class="comment">% 4. One reason FamilyName is added to the output structure so that getdata can be</span>
198	0072 <span class="comment">% used to locate archived dispersion measurements.</span>
199	0073 <span class="comment">% 5. Units for DeltaRF depend on the 'Physics' or 'Hardware' flags</span>
200	0074 <span class="comment">% 6. Beware of what units you are working in. The default units for chromaticity</span>
201	0075 <span class="comment">% are physics units. This is an exception to the normal middle layer convention.</span>
202	0076 <span class="comment">% Hardware units for "chromaticity" is in tune change per change in RF frequency.</span>
203	0077 <span class="comment">% Since this is an unusual unit to work with, the default units for chromaticity</span>
204	0078 <span class="comment">% is physics units. Note that goal chromaticity is also stored in physics units.</span>
205	0079 <span class="comment">% plotchro can switch between 'Hardware' and 'Physics' after the measurement is taken.</span>
206	0080 <span class="comment">% As an example of the difference between the units, at Spear3 1 unit of chromaticity</span>
207	0081 <span class="comment">% in physics units corresponds to roughly -1.8 units in hardware units.</span>
208	0082 <span class="comment">%</span>
209	0083 <span class="comment">% See also plotchro, measdisp</span>
210	0084 <span class="comment">%</span>
211	0085 <span class="comment">% Written by Greg Portmann and Jeff Corbett</span>
212	0086
213	0087 <span class="comment">%</span>
214	0088 <span class="comment">% Written by Gregory J. Portmann</span>
215	0089 <span class="comment">% Modified by Laurent S. Nadolski</span>
216	0090
217	0091 NRFSteps = 1;
218	0092 WaitFlag = 4;
219	0093 MCF = [];
220	0094 BPMxFamily = <a href="gethbpmfamily.html" class="code" title="function Family = gethbpmfamily">gethbpmfamily</a>; <span class="comment">% Just an extra monitor</span>
221	0095 BPMyFamily = <a href="getvbpmfamily.html" class="code" title="function Family = getvbpmfamily">getvbpmfamily</a>; <span class="comment">% Just an extra monitor</span>
222	0096 StructOutputFlag = 0;
223	0097 FileName = -1;
224	0098 ArchiveFlag = 0;
225	0099 DisplayFlag = -1;
226	0100 ModeFlag = <span class="string">''</span>; <span class="comment">% model, online, manual, or '' for default mode</span>
227	0101 UnitsFlag = <span class="string">'Physics'</span>; <span class="comment">% hardware, physics, or '' for default units</span>
228	0102
229	0103
230	0104 <span class="comment">% Look if 'struct' or 'numeric' in on the input line</span>
231	0105 <span class="keyword">for</span> i = length(varargin):-1:1
232	0106 <span class="keyword">if</span> strcmpi(varargin{i},<span class="string">'Struct'</span>)
233	0107 StructOutputFlag = 1;
234	0108 varargin(i) = [];
235	0109 <span class="keyword">elseif</span> iscell(varargin{i})
236	0110 <span class="comment">% Ignor cells</span>
237	0111 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Numeric'</span>)
238	0112 StructOutputFlag = 0;
239	0113 varargin(i) = [];
240	0114 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Archive'</span>)
241	0115 ArchiveFlag = 1;
242	0116 <span class="keyword">if</span> length(varargin) > i
243	0117 <span class="comment">% Look for a filename as the next input</span>
244	0118 <span class="keyword">if</span> ischar(varargin{i+1})
245	0119 FileName = varargin{i+1};
246	0120 varargin(i+1) = [];
247	0121 <span class="keyword">end</span>
248	0122 <span class="keyword">end</span>
249	0123 varargin(i) = [];
250	0124 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'NoArchive'</span>)
251	0125 ArchiveFlag = 0;
252	0126 varargin(i) = [];
253	0127 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Zeta'</span>) \|\| strcmpi(varargin{i},<span class="string">'Physics'</span>)
254	0128 UnitsFlag = <span class="string">'Physics'</span>;
255	0129 varargin(i) = [];
256	0130 <span class="keyword">if</span> length(varargin) >= i
257	0131 <span class="keyword">if</span> isnumeric(varargin{i})
258	0132 MCF = varargin{i};
259	0133 varargin(i) = [];
260	0134 <span class="keyword">if</span> any(size(MCF)>1)
261	0135 error(<span class="string">'Input MCF must be a scalar'</span>);
262	0136 <span class="keyword">end</span>
263	0137 <span class="keyword">end</span>
264	0138 <span class="keyword">end</span>
265	0139 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Hardware'</span>)
266	0140 UnitsFlag = <span class="string">'Hardware'</span>;
267	0141 varargin(i) = [];
268	0142 <span class="keyword">if</span> length(varargin) >= i
269	0143 <span class="keyword">if</span> isnumeric(varargin{i})
270	0144 MCF = varargin{i};
271	0145 varargin(i) = [];
272	0146 <span class="keyword">if</span> any(size(MCF)>1)
273	0147 error(<span class="string">'Input MCF must be a scalar'</span>);
274	0148 <span class="keyword">end</span>
275	0149 <span class="keyword">end</span>
276	0150 <span class="keyword">end</span>
277	0151 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Simulator'</span>) \|\| strcmpi(varargin{i},<span class="string">'Model'</span>) \|\| strcmpi(varargin{i},<span class="string">'Online'</span>) \|\| strcmpi(varargin{i},<span class="string">'Manual'</span>)
278	0152 ModeFlag = varargin{i};
279	0153 varargin(i) = [];
280	0154 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'NoDisplay'</span>)
281	0155 DisplayFlag = 0;
282	0156 varargin(i) = [];
283	0157 <span class="keyword">elseif</span> strcmpi(varargin{i},<span class="string">'Display'</span>)
284	0158 DisplayFlag = 1;
285	0159 varargin(i) = [];
286	0160 <span class="keyword">end</span>
287	0161 <span class="keyword">end</span>
288	0162
289	0163 <span class="comment">% Default for Model is no display, no archive</span>
290	0164 <span class="keyword">if</span> DisplayFlag == -1 && strcmpi(ModeFlag,<span class="string">'Model'</span>)
291	0165 DisplayFlag = 0;
292	0166 <span class="keyword">end</span>
293	0167 <span class="keyword">if</span> ArchiveFlag == -1 && strcmpi(ModeFlag,<span class="string">'Model'</span>)
294	0168 ArchiveFlag = 0;
295	0169 <span class="keyword">end</span>
296	0170
297	0171
298	0172 <span class="comment">% DeltaRF input</span>
299	0173 <span class="keyword">if</span> length(varargin) >= 1
300	0174 <span class="keyword">if</span> isnumeric(varargin{1})
301	0175 DeltaRF = varargin{1};
302	0176 <span class="keyword">else</span>
303	0177 DeltaRF = [];
304	0178 <span class="keyword">end</span>
305	0179 <span class="keyword">else</span>
306	0180 DeltaRF = [];
307	0181 <span class="keyword">end</span>
308	0182
309	0183 <span class="comment">% WaitFlag input</span>
310	0184 <span class="keyword">if</span> length(varargin) >= 2
311	0185 WaitFlag = varargin{2};
312	0186 <span class="keyword">end</span>
313	0187 <span class="keyword">if</span> isempty(WaitFlag) \|\| WaitFlag == -3
314	0188 WaitFlag = 2.2 * <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'TuneDelay'</span>);
315	0189 <span class="keyword">end</span>
316	0190 <span class="keyword">if</span> isempty(WaitFlag)
317	0191 WaitFlag = input(<span class="string">' Delay for Tune Measurement (Seconds, Keyboard Pause = -4, or Manual Tune Input = -5) = '</span>);
318	0192 <span class="keyword">end</span>
319	0193
320	0194
321	0195 <span class="comment">% Archive data structure</span>
322	0196 <span class="keyword">if</span> ArchiveFlag
323	0197 <span class="keyword">if</span> isempty(FileName)
324	0198 FileName = appendtimestamp(<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Default'</span>, <span class="string">'ChroArchiveFile'</span>));
325	0199 DirectoryName = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Directory'</span>,<span class="string">'ChroData'</span>);
326	0200 <span class="keyword">if</span> isempty(DirectoryName)
327	0201 DirectoryName = [<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Directory'</span>,<span class="string">'DataRoot'</span>) <span class="string">'Chromaticity'</span>, filesep];
328	0202 <span class="keyword">else</span>
329	0203 <span class="comment">% Make sure default directory exists</span>
330	0204 DirStart = pwd;
331	0205 [DirectoryName, ErrorFlag] = gotodirectory(DirectoryName);
332	0206 cd(DirStart);
333	0207 <span class="keyword">end</span>
334	0208 [FileName, DirectoryName] = uiputfile(<span class="string">'*.mat'</span>, <span class="string">'Select Chromaticity File'</span>, [DirectoryName FileName]);
335	0209 <span class="keyword">if</span> FileName == 0
336	0210 ArchiveFlag = 0;
337	0211 disp(<span class="string">' Chromaticity measurement canceled.'</span>);
338	0212 Chromaticity=[]; FileName=<span class="string">''</span>;
339	0213 <span class="keyword">return</span>
340	0214 <span class="keyword">end</span>
341	0215 FileName = [DirectoryName, FileName];
342	0216 <span class="keyword">elseif</span> FileName == -1
343	0217 FileName = appendtimestamp(<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Default'</span>, <span class="string">'ChroArchiveFile'</span>));
344	0218 DirectoryName = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Directory'</span>,<span class="string">'ChroData'</span>);
345	0219 <span class="keyword">if</span> isempty(DirectoryName)
346	0220 DirectoryName = [<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'Directory'</span>,<span class="string">'DataRoot'</span>) <span class="string">'Chromaticity'</span>, filesep];
347	0221 <span class="keyword">end</span>
348	0222 FileName = [DirectoryName, FileName];
349	0223 <span class="keyword">end</span>
350	0224 <span class="keyword">end</span>
351	0225
352	0226
353	0227 <span class="comment">% Get units from the RF frequency</span>
354	0228 <span class="keyword">if</span> isempty(UnitsFlag)
355	0229 UnitsFlag = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'Units'</span>);
356	0230 <span class="keyword">end</span>
357	0231
358	0232
359	0233 <span class="comment">% if isempty(ModeFlag)</span>
360	0234 <span class="comment">% if strcmpi(getfamilydata('RF','Setpoint','Mode'), getfamilydata(BPMxFamily,'Monitor','Mode'))</span>
361	0235 <span class="comment">% ModeFlag = getfamilydata(BPMxFamily,'Monitor','Mode');</span>
362	0236 <span class="comment">% else</span>
363	0237 <span class="comment">% error('Mix Mode for RF and orbits');</span>
364	0238 <span class="comment">% end</span>
365	0239 <span class="comment">% end</span>
366	0240
367	0241
368	0242
369	0243 <span class="keyword">if</span> strcmpi(UnitsFlag,<span class="string">'Hardware'</span>)
370	0244 RFUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'HWUnits'</span>);
371	0245 <span class="keyword">elseif</span> strcmpi(UnitsFlag,<span class="string">'Physics'</span>)
372	0246 RFUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'PhysicsUnits'</span>);
373	0247 <span class="keyword">else</span>
374	0248 error(<span class="string">'RF units unknown. Inputs DeltaRF directly.'</span>);
375	0249 <span class="keyword">end</span>
376	0250
377	0251 <span class="comment">% DeltaRF default</span>
378	0252 <span class="keyword">if</span> isempty(DeltaRF)
379	0253 <span class="comment">% Get the default from the AD is in Hardware units</span>
380	0254 DeltaRF = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'DeltaRFChro'</span>);
381	0255
382	0256 <span class="comment">% If the default is not in the AD</span>
383	0257 <span class="keyword">if</span> isempty(DeltaRF)
384	0258 DeltaRF = <a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a>(<span class="string">'Hardware'</span>) * <a href="getmcf.html" class="code" title="function Alpha = getmcf(ModelString)">getmcf</a> * [-.004 -.002 0 .002 .004] ; <span class="comment">% .2% energy change per step</span>
385	0259 <span class="comment">% DeltaRF = [-100 -50 0 50 150]; % Hz</span>
386	0260 <span class="comment">%if strcmpi(RFUnitsString, 'Hz')</span>
387	0261 <span class="comment">% % Default units OK</span>
388	0262 <span class="comment">%elseif strcmpi(RFUnitsString, 'kHz')</span>
389	0263 <span class="comment">% % Change to kHz</span>
390	0264 <span class="comment">% DeltaRF = DeltaRF / 1e3;</span>
391	0265 <span class="comment">%elseif strcmpi(RFUnitsString, 'MHz')</span>
392	0266 <span class="comment">% % Change to MHz</span>
393	0267 <span class="comment">% DeltaRF = DeltaRF / 1e6;</span>
394	0268 <span class="comment">%else</span>
395	0269 <span class="comment">% error('RF units unknown. Input DeltaRF directly or put the default in AD.DeltaRFChro.');</span>
396	0270 <span class="comment">%end</span>
397	0271 <span class="keyword">else</span>
398	0272 <span class="keyword">if</span> strcmpi(UnitsFlag,<span class="string">'Physics'</span>)
399	0273 <span class="comment">% Since the default from the AO must be in hardware units, change to physics units</span>
400	0274 DeltaRF = <a href="hw2physics.html" class="code" title="function S = hw2physics(Family, Field, value, DeviceList, Energy)">hw2physics</a>(<span class="string">'RF'</span>, <span class="string">'Setpoint'</span>, DeltaRF, [1 1], ModeFlag);
401	0275 <span class="keyword">end</span>
402	0276 <span class="keyword">end</span>
403	0277 <span class="keyword">end</span>
404	0278
405	0279
406	0280 <span class="comment">% Check DeltaRF for resonable values</span>
407	0281 <span class="keyword">if</span> strcmpi(RFUnitsString, <span class="string">'MHz'</span>)
408	0282 <span class="keyword">if</span> abs(max(DeltaRF)-min(DeltaRF)) > .020; <span class="comment">% .020 MHz</span>
409	0283 tmp = questdlg(sprintf(<span class="string">'%f MHz is a large RF change. Do you want to continue?'</span>, abs(max(DeltaRF)-min(DeltaRF))),<span class="string">'Dispersion Measurement'</span>,<span class="string">'YES'</span>,<span class="string">'NO'</span>,<span class="string">'YES'</span>);
410	0284 <span class="keyword">if</span> strcmp(tmp,<span class="string">'NO'</span>)
411	0285 Chromaticity=[];
412	0286 <span class="keyword">return</span>
413	0287 <span class="keyword">end</span>
414	0288 <span class="keyword">end</span>
415	0289 <span class="keyword">elseif</span> strcmpi(RFUnitsString, <span class="string">'kHz'</span>)
416	0290 <span class="keyword">if</span> abs(max(DeltaRF)-min(DeltaRF)) > 20; <span class="comment">% kHz</span>
417	0291 tmp = questdlg(sprintf(<span class="string">'%f kHz is a large RF change. Do you want to continue?'</span>, abs(max(DeltaRF)-min(DeltaRF))),<span class="string">'Dispersion Measurement'</span>,<span class="string">'YES'</span>,<span class="string">'NO'</span>,<span class="string">'YES'</span>);
418	0292 <span class="keyword">if</span> strcmp(tmp,<span class="string">'NO'</span>)
419	0293 Chromaticity=[];
420	0294 <span class="keyword">return</span>
421	0295 <span class="keyword">end</span>
422	0296 <span class="keyword">end</span>
423	0297 <span class="keyword">elseif</span> strcmpi(RFUnitsString, <span class="string">'Hz'</span>)
424	0298 <span class="keyword">if</span> abs(max(DeltaRF)-min(DeltaRF)) > 20000; <span class="comment">% Hz</span>
425	0299 tmp = questdlg(sprintf(<span class="string">'%f Hz is a large RF change. Do you want to continue?'</span>, abs(max(DeltaRF)-min(DeltaRF))),<span class="string">'Dispersion Measurement'</span>,<span class="string">'YES'</span>,<span class="string">'NO'</span>,<span class="string">'YES'</span>);
426	0300 <span class="keyword">if</span> strcmp(tmp,<span class="string">'NO'</span>)
427	0301 Chromaticity=[];
428	0302 <span class="keyword">return</span>
429	0303 <span class="keyword">end</span>
430	0304 <span class="keyword">end</span>
431	0305 <span class="keyword">else</span>
432	0306 <span class="comment">% Don't who how to check, hence no check made</span>
433	0307 <span class="keyword">end</span>
434	0308
435	0309 <span class="comment">% DeltaRF must be in "RFUnitsString" units at this point</span>
436	0310
437	0311
438	0312 RFsp = <a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a>(<span class="string">'Struct'</span>, UnitsFlag, ModeFlag);
439	0313
440	0314 <span class="keyword">if</span> isempty(MCF)
441	0315 MCF = <a href="getmcf.html" class="code" title="function Alpha = getmcf(ModelString)">getmcf</a>(ModeFlag);
442	0316 <span class="keyword">end</span>
443	0317
444	0318
445	0319 <span class="comment">% Fill the chromaticity structure (response matrix structure + some fields)</span>
446	0320 c.Data = [];
447	0321 c.FamilyName = <span class="string">'Chromaticity'</span>;
448	0322 <span class="keyword">if</span> <a href="isfamily.html" class="code" title="function [FamilyFlag, AO] = isfamily(Family, Field)">isfamily</a>(<span class="string">'TUNE'</span>)
449	0323 c.Monitor = <a href="family2datastruct.html" class="code" title="function [DataStruct, ErrorFlag] = family2datastruct(varargin)">family2datastruct</a>(<span class="string">'TUNE'</span>,<span class="string">'Monitor'</span>,[1 1;1 2]);
450	0324 <span class="keyword">else</span>
451	0325 c.Monitor = <a href="gettune.html" class="code" title="function [Tune, tout, DataTime, ErrorFlag] = gettune(varargin)">gettune</a>(<span class="string">'Struct'</span>, <span class="string">'Model'</span>); <span class="comment">% Just to fill the structure</span>
452	0326 c.Monitor.Data = NaN * c.Monitor.Data;
453	0327 <span class="keyword">end</span>
454	0328 c.Monitor = <a href="gettune.html" class="code" title="function [Tune, tout, DataTime, ErrorFlag] = gettune(varargin)">gettune</a>(<span class="string">'Struct'</span>, ModeFlag); <span class="comment">% Should WaitFlag before to be sure to get a good measurement ?????</span>
455	0329 c.Actuator = RFsp;
456	0330 c.ActuatorDelta = DeltaRF;
457	0331 c.GeV = <a href="getenergy.html" class="code" title="function [Energy, HCMEnergy] = getenergy(varargin)">getenergy</a>(ModeFlag);
458	0332 c.DCCT = <a href="getam.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getam(varargin)">getam</a>(<span class="string">'DCCT'</span>, ModeFlag);
459	0333 c.ModulationMethod = <span class="string">'Unipolar'</span>;
460	0334 c.WaitFlag = WaitFlag;
461	0335 c.TimeStamp = clock;
462	0336 c.Mode = ModeFlag;
463	0337 c.Units = UnitsFlag;
464	0338 c.UnitsString = [];
465	0339 c.DataDescriptor = <span class="string">'Chromaticity'</span>;
466	0340 c.CreatedBy = <span class="string">'measchro'</span>;
467	0341 c.OperationalMode = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'OperationalMode'</span>);
468	0342
469	0343 <span class="comment">% Nonstandard response matrix fields</span>
470	0344 <span class="keyword">if</span> strcmpi(ModeFlag,<span class="string">'Manual'</span>)
471	0345 c.X = NaN;
472	0346 c.Y = NaN;
473	0347 <span class="keyword">else</span>
474	0348 c.X = <a href="getx.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = getx(varargin)">getx</a>(<span class="string">'Struct'</span>, UnitsFlag, ModeFlag);
475	0349 c.Y = <a href="gety.html" class="code" title="function [AM, tout, DataTime, ErrorFlag] = gety(varargin)">gety</a>(<span class="string">'Struct'</span>, UnitsFlag, ModeFlag);
476	0350 <span class="keyword">end</span>
477	0351 c.MCF = MCF;
478	0352 RF0 = RFsp.Data(1);
479	0353 c.dp = -DeltaRF / (RF0*MCF);
480	0354
481	0355
482	0356 <span class="keyword">if</span> strcmpi(ModeFlag,<span class="string">'Model'</span>) \|\| strcmpi(ModeFlag,<span class="string">'Simulator'</span>)
483	0357 <span class="comment">% No need for delays with the model</span>
484	0358 WaitFlag = 0;
485	0359 ExtraDelay = 0;
486	0360 <span class="keyword">end</span>
487	0361
488	0362
489	0363 <span class="keyword">if</span> strcmpi(ModeFlag,<span class="string">'Model'</span>)
490	0364 c.Data = modelchro(<span class="string">'Physics'</span>);
491	0365
492	0366 <span class="keyword">if</span> strcmpi(UnitsFlag,<span class="string">'Physics'</span>)
493	0367 TuneUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'TUNE'</span>,<span class="string">'Monitor'</span>,<span class="string">'PhysicsUnits'</span>);
494	0368 <span class="keyword">if</span> isempty(TuneUnitsString)
495	0369 TuneUnitsString = <span class="string">'Fractional Tune'</span>;
496	0370 <span class="keyword">end</span>
497	0371 c.UnitsString = [TuneUnitsString, <span class="string">'/(dp/p)'</span>];
498	0372 <span class="keyword">else</span>
499	0373 <span class="comment">% Tune Shift vs. RF Frequency</span>
500	0374 c.Data(1,1) = c.Data(1,1) / (-RF0 * MCF);
501	0375 c.Data(2,1) = c.Data(2,1) / (-RF0 * MCF);
502	0376 TuneUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'TUNE'</span>,<span class="string">'Monitor'</span>,<span class="string">'HWUnits'</span>);
503	0377 <span class="keyword">if</span> isempty(TuneUnitsString)
504	0378 TuneUnitsString = <span class="string">'Fractional Tune'</span>;
505	0379 <span class="keyword">end</span>
506	0380 c.UnitsString = [TuneUnitsString, <span class="string">'/'</span>,<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'HWUnits'</span>)];
507	0381 <span class="keyword">end</span>
508	0382
509	0383 <span class="keyword">else</span>
510	0384 <span class="comment">% Online or Simulator</span>
511	0385 <span class="comment">% Start measurement</span>
512	0386 <span class="keyword">if</span> DisplayFlag
513	0387 fprintf(<span class="string">' Begin chromaticity measurement\n'</span>);
514	0388 <span class="keyword">end</span>
515	0389 <span class="keyword">for</span> i = 1:length(DeltaRF)
516	0390 <span class="comment">%setrf(RF0 + DeltaRF(i), UnitsFlag, ModeFlag);</span>
517	0391 <span class="keyword">if</span> (isempty(ModeFlag) && strcmpi(<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'Mode'</span>),<span class="string">'Manual'</span>)) \|\| strcmpi(ModeFlag,<span class="string">'Manual'</span>)
518	0392 <span class="comment">% One shot setting of RF</span>
519	0393 <a href="setrf.html" class="code" title="function setrf(RF, varargin)">setrf</a>(RF0 + DeltaRF(i), UnitsFlag, ModeFlag);
520	0394 <span class="keyword">else</span>
521	0395 <span class="comment">% Slow setting of RF</span>
522	0396 rf = <a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a>(UnitsFlag, ModeFlag);
523	0397 <span class="keyword">for</span> k = 1:NRFSteps
524	0398 <a href="setsp.html" class="code" title="function ErrorFlag = setsp(Family, varargin)">setsp</a>(<span class="string">'RF'</span>, rf + k/NRFSteps * (RF0+DeltaRF(i)-rf), [], -1, UnitsFlag, ModeFlag);
525	0399 pause(0.1);
526	0400 <span class="keyword">end</span>
527	0401 <span class="keyword">end</span>
528	0402
529	0403 RF(:,i) = <a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a>(UnitsFlag, ModeFlag);
530	0404 <span class="keyword">if</span> DisplayFlag
531	0405 fprintf(<span class="string">' %d. RF frequency is %.5f\n'</span>, i, RF(:,i));
532	0406 <span class="keyword">end</span>
533	0407
534	0408 <span class="comment">% Wait for tune monitor to have fresh data</span>
535	0409 <span class="keyword">if</span> WaitFlag >= 0
536	0410 <span class="keyword">if</span> DisplayFlag && ~strcmpi(ModeFlag,<span class="string">'Manual'</span>)
537	0411 fprintf(<span class="string">' Pausing %f seconds for the tune measurement\n'</span>, WaitFlag);
538	0412 pause(0);
539	0413 <span class="keyword">end</span>
540	0414 sleep(WaitFlag);
541	0415 Tune(:,i) = <a href="gettune.html" class="code" title="function [Tune, tout, DataTime, ErrorFlag] = gettune(varargin)">gettune</a>(ModeFlag);
542	0416 <span class="keyword">elseif</span> WaitFlag == -4
543	0417 tmp = input(<span class="string">' Hit return when the tune measurement is ready. '</span>);
544	0418 Tune(:,i) = <a href="gettune.html" class="code" title="function [Tune, tout, DataTime, ErrorFlag] = gettune(varargin)">gettune</a>(ModeFlag);
545	0419 <span class="keyword">elseif</span> WaitFlag == -5
546	0420 Tune(1,i) = input(<span class="string">' Input the horizontal tune = '</span>);
547	0421 Tune(2,i) = input(<span class="string">' Input the vertical tune = '</span>);
548	0422 <span class="keyword">else</span>
549	0423 error(<span class="string">'Tune delay method unknown'</span>);
550	0424 <span class="keyword">end</span>
551	0425
552	0426 <span class="comment">%if any(isnan(Tune))</span>
553	0427 <span class="comment">% fprintf(' Chromaticity measurement failed. RF frequency reset.\n');</span>
554	0428 <span class="comment">% setrf(RF0, UnitsFlag, ModeFlag);</span>
555	0429 <span class="comment">% Chromaticity = [NaN; NaN];</span>
556	0430 <span class="comment">% return;</span>
557	0431 <span class="comment">%end</span>
558	0432 <span class="keyword">end</span>
559	0433
560	0434
561	0435 <span class="comment">% Reset RF</span>
562	0436 <span class="comment">%setrf(RF0, UnitsFlag, ModeFlag);</span>
563	0437 <span class="keyword">if</span> isempty(ModeFlag) && strcmpi(<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'Mode'</span>),<span class="string">'Manual'</span>)
564	0438 <span class="comment">% One shot setting of RF</span>
565	0439 <a href="setrf.html" class="code" title="function setrf(RF, varargin)">setrf</a>(RF0, UnitsFlag, ModeFlag);
566	0440 <span class="keyword">else</span>
567	0441 <span class="comment">% Slow setting of RF</span>
568	0442 rf = <a href="getrf.html" class="code" title="function [RFsp, RFam, DataTime, ErrorFlag] = getrf(varargin)">getrf</a>(UnitsFlag, ModeFlag);
569	0443 <span class="keyword">for</span> k = 1:NRFSteps
570	0444 <a href="setsp.html" class="code" title="function ErrorFlag = setsp(Family, varargin)">setsp</a>(<span class="string">'RF'</span>, rf + k/NRFSteps * (RF0-rf), [], -1, UnitsFlag, ModeFlag);
571	0445 pause(0.1);
572	0446 <span class="keyword">end</span>
573	0447 <span class="keyword">end</span>
574	0448
575	0449
576	0450 <span class="comment">% Load Tune measurements into the chromaticy structure</span>
577	0451 c.Tune = Tune;
578	0452
579	0453 <span class="keyword">if</span> strcmpi(UnitsFlag,<span class="string">'Physics'</span>)
580	0454 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
581	0455 <span class="comment">% Tune Shift vs. Momentum %</span>
582	0456 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
583	0457
584	0458 <span class="comment">% Horizontal tune vs. momentum</span>
585	0459 p = polyfit(c.dp, Tune(1,:), 2); <span class="comment">%2nd order polynomial fit to data</span>
586	0460 c.PolyFit(1,:) = p;
587	0461 c.Data(1,1) = p(2);
588	0462
589	0463 <span class="comment">% Vertical tune vs. rf frequency</span>
590	0464 p = polyfit(c.dp, Tune(2,:), 2);
591	0465 c.PolyFit(2,:) = p;
592	0466 c.Data(2,1) = p(2);
593	0467
594	0468 TuneUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'TUNE'</span>,<span class="string">'Monitor'</span>,<span class="string">'PhysicsUnits'</span>);
595	0469 <span class="keyword">if</span> isempty(TuneUnitsString)
596	0470 c.UnitsString = [<span class="string">'Fractional Tune/(dp/p)'</span>];
597	0471 <span class="keyword">else</span>
598	0472 c.UnitsString = [TuneUnitsString,<span class="string">'/(dp/p)'</span>];
599	0473 <span class="keyword">end</span>
600	0474
601	0475 <span class="comment">%fprintf('\n Horizontal Chromaticity (Un-normalized) = %f \n', c.Data(1));</span>
602	0476 <span class="comment">%fprintf(' Vertical Chromaticity (Un-normalized) = %f \n' , c.Data(2));</span>
603	0477 <span class="keyword">else</span>
604	0478 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
605	0479 <span class="comment">% Tune Shift vs. RF Frequency %</span>
606	0480 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
607	0481
608	0482 <span class="comment">% Horizontal tune vs. rf frequency</span>
609	0483 p = polyfit(DeltaRF, Tune(1,:), 2); <span class="comment">% 2nd order polynomial fit to data</span>
610	0484 c.PolyFit(1,:) = p;
611	0485 c.Data(1,1) = p(2);
612	0486
613	0487 <span class="comment">% Vertical tune vs. rf frequency</span>
614	0488 p = polyfit(DeltaRF, Tune(2,:), 2);
615	0489 c.PolyFit(2,:) = p;
616	0490 c.Data(2,1) = p(2);
617	0491
618	0492 TuneUnitsString = <a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'TUNE'</span>,<span class="string">'Monitor'</span>,<span class="string">'HWUnits'</span>);
619	0493 <span class="keyword">if</span> isempty(TuneUnitsString)
620	0494 c.UnitsString = [<span class="string">'Fractional Tune/'</span>,<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'HWUnits'</span>)];
621	0495 <span class="keyword">else</span>
622	0496 c.UnitsString = [TuneUnitsString,<span class="string">'/'</span>,<a href="getfamilydata.html" class="code" title="function [Data, ErrorFlag] = getfamilydata(Family, Field1, Field2, DeviceList)">getfamilydata</a>(<span class="string">'RF'</span>,<span class="string">'Setpoint'</span>,<span class="string">'HWUnits'</span>)];
623	0497 <span class="keyword">end</span>
624	0498 <span class="keyword">end</span>
625	0499 <span class="keyword">end</span>
626	0500
627	0501 <span class="keyword">if</span> DisplayFlag
628	0502 fprintf(<span class="string">' Chromaticity = %f [%s]\n'</span>, c.Data(1), c.UnitsString);
629	0503 fprintf(<span class="string">' Chromaticity = %f [%s]\n'</span>, c.Data(2), c.UnitsString);
630	0504 <span class="keyword">end</span>
631	0505
632	0506 <span class="keyword">if</span> DisplayFlag && ~strcmpi(ModeFlag,<span class="string">'Model'</span>)
633	0507 <span class="comment">%figure;</span>
634	0508 <a href="plotchro.html" class="code" title="function [c, FileName] = plotchro(varargin)">plotchro</a>(c);
635	0509 <span class="keyword">end</span>
636	0510
637	0511
638	0512 <span class="comment">% Archive data structure</span>
639	0513 <span class="keyword">if</span> ArchiveFlag
640	0514 <span class="comment">% If the filename contains a directory then make sure it exists</span>
641	0515 [DirectoryName, FileName, Ext] = fileparts(FileName);
642	0516 DirStart = pwd;
643	0517 [DirectoryName, ErrorFlag] = gotodirectory(DirectoryName);
644	0518 Chromaticity = c;
645	0519 save(FileName, <span class="string">'Chromaticity'</span>);
646	0520 <span class="keyword">if</span> DisplayFlag
647	0521 fprintf(<span class="string">' Chromaticity data saved to %s.mat\n'</span>, [DirectoryName FileName]);
648	0522 <span class="keyword">if</span> ErrorFlag
649	0523 fprintf(<span class="string">' Warning: %s was not the desired directory\n'</span>, DirectoryName);
650	0524 <span class="keyword">end</span>
651	0525 <span class="keyword">end</span>
652	0526 cd(DirStart);
653	0527 FileName = [DirectoryName, FileName, <span class="string">'.mat'</span>];
654	0528 <span class="keyword">end</span>
655	0529 <span class="keyword">if</span> FileName == -1
656	0530 FileName = <span class="string">''</span>;
657	0531 <span class="keyword">end</span>
658	0532
659	0533
660	0534 <span class="comment">% Load output data</span>
661	0535 <span class="keyword">if</span> StructOutputFlag
662	0536 Chromaticity = c;
663	0537 <span class="keyword">else</span>
664	0538 Chromaticity = c.Data;
665	0539 <span class="keyword">end</span>
666	0540
667	0541
668	0542 <span class="keyword">if</span> DisplayFlag
669	0543 fprintf(<span class="string">' Chromaticity measurement is complete.\n'</span>);
670	0544 <span class="keyword">end</span>
671	0545
672	0546</pre></div>
673	<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>
674	</body>
675	</html>

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: MML/trunk/mml/doc_html/mml/measchro.html @ 4

Download in other formats: