source: MML/trunk/machine/SOLEIL/doc_html/machine/Soleil/StorageRing/stepalpha2.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: 7.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 stepalpha2</title>
6  <meta name="keywords" content="stepalpha2">
7  <meta name="description" content="FUNCTION STEPALPHA2 - step chromaticities and alpha2 using sextupole families">
8  <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
9  <meta name="generator" content="m2html v1.5 &copy; 2003-2005 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> &gt;  <a href="#">machine</a> &gt; <a href="#">Soleil</a> &gt; <a href="index.html">StorageRing</a> &gt; stepalpha2.m</div>
16
17<!--<table width="100%"><tr><td align="left"><a href="../../../index.html"><img alt="<" border="0" src="../../../left.png">&nbsp;Master index</a></td>
18<td align="right"><a href="index.html">Index for machine/Soleil/StorageRing&nbsp;<img alt=">" border="0" src="../../../right.png"></a></td></tr></table>-->
19
20<h1>stepalpha2
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>FUNCTION STEPALPHA2 - step chromaticities and alpha2 using sextupole families</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 stepalpha2(dxix, dxiz, dalpha2) </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"> FUNCTION STEPALPHA2 - step chromaticities and alpha2 using sextupole families
31
32  INPUTS
33  1. dxix - horizontal tune variation
34  2. dxiz - vertical tune variation
35  3. dalpha2 - second order momentum compaction variation
36
37  Notes
38  1. Matrix generation
39  [DSx DSz DchroVal]   = modelchrosensitivity;
40  [SextuRmcf1 SextuRmcf2]=modelmcfsensitivity('Sextu');
41
42  See Also <a href="stepalpha1.html" class="code" title="function stepalpha1(dnux, dnuz, dalpha1, dalpha2)">stepalpha1</a>, <a href="physics_mcf.html" class="code" title="function [alpha1 alpha2 alpha3] = physics_mcf(varargin)">physics_mcf</a></pre></div>
43
44<!-- crossreference -->
45<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
46This function calls:
47<ul style="list-style-image:url(../../../matlabicon.gif)">
48</ul>
49This function is called by:
50<ul style="list-style-image:url(../../../matlabicon.gif)">
51</ul>
52<!-- crossreference -->
53
54
55
56<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
57<div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function stepalpha2(dxix, dxiz, dalpha2)</a>
580002 <span class="comment">% FUNCTION STEPALPHA2 - step chromaticities and alpha2 using sextupole families</span>
590003 <span class="comment">%</span>
600004 <span class="comment">%  INPUTS</span>
610005 <span class="comment">%  1. dxix - horizontal tune variation</span>
620006 <span class="comment">%  2. dxiz - vertical tune variation</span>
630007 <span class="comment">%  3. dalpha2 - second order momentum compaction variation</span>
640008 <span class="comment">%</span>
650009 <span class="comment">%  Notes</span>
660010 <span class="comment">%  1. Matrix generation</span>
670011 <span class="comment">%  [DSx DSz DchroVal]   = modelchrosensitivity;</span>
680012 <span class="comment">%  [SextuRmcf1 SextuRmcf2]=modelmcfsensitivity('Sextu');</span>
690013 <span class="comment">%</span>
700014 <span class="comment">%  See Also stepalpha1, physics_mcf</span>
710015
720016
730017   
740018
750019 <span class="comment">%</span>
760020 <span class="comment">%% Written by Laurent S. Nadolski</span>
770021
780022 InteractiveFlag = 1;
790023
800024 <span class="keyword">if</span> nargin &lt; 3
810025     error(<span class="string">'Number of input parameters is at least than 3!'</span>)
820026 <span class="keyword">end</span>
830027
840028 <span class="comment">% Default families to used</span>
850029 familyList = [1 4 9];
860030
870031 <span class="comment">% Load linear response matrices</span>
880032 <span class="comment">% These are the sensibility matrices expressed in amperes</span>
890033 <span class="comment">% DchroVal</span>
900034 <span class="comment">%</span>
910035
920036 dirName = getfamilydata(<span class="string">'Directory'</span>,<span class="string">'OpsData'</span>);
930037 <span class="comment">%load([dirName '..' filesep 'lowalpha_dec08' filesep 'Rmatrix_alphaby20_nouveau_modele_dec08_opt_nonlin']);</span>
940038 load([dirName <span class="string">'..'</span> filesep <span class="string">'lowalpha_dec08'</span> filesep <span class="string">'Rmatrix_alphaby40_nouveau_modele_janv10_opt_nonlin'</span>]);
950039
960040 SextuFamNumber = length(familyList);
970041 <span class="keyword">if</span> nargin == 3
980042     DI= ones(SextuFamNumber,1)*NaN;
990043     DI = pinv([DchroVal(:,familyList); SextuRmcf2(familyList)])*[dxix dxiz dalpha2]';
1000044 <span class="keyword">else</span>
1010045     error(<span class="string">'Wrong number of parameters'</span>)
1020046 <span class="keyword">end</span>
1030047
1040048 <span class="comment">%Display variation and ask confirmation</span>
1050049 fprintf(<span class="string">'\n\n Sextupole current values to be applied\n'</span>);
1060050 <span class="keyword">for</span> iSextu = 1:SextuFamNumber,
1070051     fprintf(sprintf(<span class="string">'Sextupole step for S%d is %f A\n'</span>, familyList(iSextu), DI(iSextu)));
1080052 <span class="keyword">end</span>
1090053
1100054 <span class="keyword">if</span> InteractiveFlag
1110055
1120056     <span class="comment">% Ask user for applying correction</span>
1130057     reply = input(<span class="string">'Apply values (Y/n) \n'</span>, <span class="string">'s'</span>);
1140058     <span class="keyword">if</span> isempty(reply)
1150059         reply = <span class="string">'Y'</span>;
1160060     <span class="keyword">end</span>
1170061
1180062     <span class="keyword">switch</span> reply
1190063         <span class="keyword">case</span> {<span class="string">'y'</span>, <span class="string">'Y'</span>}
1200064             <span class="keyword">for</span> iSextu = 1:SextuFamNumber,
1210065                 sprintf(<span class="string">'S%d'</span>,iSextu);
1220066                 stepsp(sprintf(<span class="string">'S%d'</span>,familyList(iSextu)), DI(iSextu));
1230067             <span class="keyword">end</span>
1240068             fprintf(<span class="string">'Correction applied\n'</span>);
1250069         <span class="keyword">otherwise</span>
1260070             fprintf(<span class="string">'Aborting ...\n'</span>);
1270071     <span class="keyword">end</span>
1280072 <span class="keyword">else</span>
1290073     <span class="keyword">for</span> iSextu = 1:SextuFamNumber,
1300074         sprintf(<span class="string">'Q%d'</span>,iSextu);
1310075         stepsp(sprintf(<span class="string">'Q%d'</span>,familyList(iSextu)), DI(iSextu));
1320076     <span class="keyword">end</span>
1330077     fprintf(<span class="string">'Correction applied\n'</span>);
1340078 <span class="keyword">end</span>
1350079</pre></div>
136<hr><address>Generated on Fri 23-Jul-2010 00:42:07 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/" title="Matlab Documentation in HTML">m2html</a></strong> &copy; 2005</address>
137</body>
138</html>
Note: See TracBrowser for help on using the repository browser.