source: MML/trunk/machine/SOLEIL/doc_html/machine/Soleil/StorageRing/modelemit.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: 4.6 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 modelemit</title>
6  <meta name="keywords" content="modelemit">
7  <meta name="description" content="MODELEMIT - Returns emittance using Ohmi's enveloppe formalism">
8  <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
9  <meta name="generator" content="m2html &copy; 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> &gt;  <a href="#">machine</a> &gt; <a href="#">Soleil</a> &gt; <a href="index.html">StorageRing</a> &gt; modelemit.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>modelemit
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>MODELEMIT - Returns emittance using Ohmi's enveloppe formalism</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 [emit sigmas] = modelemit(TwissData) </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">MODELEMIT - Returns emittance using Ohmi's enveloppe formalism
31
32  INPUTS
33  1. TwissData
34
35  OUPUTS
36  1. emit - emittance vector in nm.rad
37  2. sigmas - Sigma dimension around the ring</pre></div>
38
39<!-- crossreference -->
40<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
41This function calls:
42<ul style="list-style-image:url(../../../matlabicon.gif)">
43</ul>
44This function is called by:
45<ul style="list-style-image:url(../../../matlabicon.gif)">
46</ul>
47<!-- crossreference -->
48
49
50<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
51<div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [emit sigmas] = modelemit(TwissData)</a>
520002 <span class="comment">%MODELEMIT - Returns emittance using Ohmi's enveloppe formalism</span>
530003 <span class="comment">%</span>
540004 <span class="comment">%  INPUTS</span>
550005 <span class="comment">%  1. TwissData</span>
560006 <span class="comment">%</span>
570007 <span class="comment">%  OUPUTS</span>
580008 <span class="comment">%  1. emit - emittance vector in nm.rad</span>
590009 <span class="comment">%  2. sigmas - Sigma dimension around the ring</span>
600010
610011 <span class="comment">% TODO complete the documentation</span>
620012
630013 <span class="comment">% Written by M. Munoz</span>
640014 <span class="comment">% Modified by Laurent S. Nadolski</span>
650015
660016 <span class="keyword">global</span>  THERING;
670017 <span class="keyword">if</span> (nargin==0)
680018     [TwissData, tune]  = twissring(THERING,0,1:length(THERING)+1,<span class="string">'chrom'</span>);
690019 <span class="keyword">end</span>
700020
710021
720022 ati = atindex(THERING);
730023 ibend = sort([ati.BEND]);
740024 <span class="keyword">for</span> i = ibend
750025     THERING{i}.PassMethod = <span class="string">'BndMPoleSymplectic4RadPass'</span>;
760026     <span class="comment">%THERING{i}.Energy = getenergymodel*1e9;</span>
770027 <span class="keyword">end</span>
780028
790029 RADELEMINDEX = sort([ibend]);
800030
810031 [ENV, DP, DL] = ohmienvelope(THERING,RADELEMINDEX, 1:length(THERING)+1);
820032 sigmas = cat(2,ENV.Sigma);
830033
840034
850035 [TwissData, tune]  = twissring(THERING,0,1:length(THERING)+1,<span class="string">'chrom'</span>);
860036 beta = cat(1,TwissData.beta);
870037 eta  = cat(2,TwissData.Dispersion);
880038
890039 epsx = (sigmas(1,:).^2-eta(1,:).^2*DP^2)./beta(:,1)';
900040 epsy = (sigmas(2,:).^2-eta(3,:).^2*DP^2)./beta(:,2)';
910041 emit(1) = median(epsx);
920042 emit(2) = median(epsy);
930043 <span class="keyword">for</span> i = ibend
940044     THERING{i}.PassMethod = <span class="string">'BndMPoleSymplectic4Pass'</span>;
950045     <span class="comment">%THERING{i}.Energy = getenergymodel*1e9;</span>
960046 <span class="keyword">end</span>
970047
980048 emit = emit*1e9;</pre></div>
99<hr><address>Generated on Mon 21-May-2007 15:35:27 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> &copy; 2003</address>
100</body>
101</html>
Note: See TracBrowser for help on using the repository browser.