<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Has anyone else noticed that the disk chopper component in 2.2a doesn’t block the beam?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The following doesn’t work on RHEL 6.4<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">DEFINE INSTRUMENT test(Par1=1)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">DECLARE<o:p></o:p></p>
<p class="MsoNormal">%{<o:p></o:p></p>
<p class="MsoNormal">%}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">INITIALIZE<o:p></o:p></p>
<p class="MsoNormal">%{<o:p></o:p></p>
<p class="MsoNormal">%}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">TRACE<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">COMPONENT Origin = Progress_bar()<o:p></o:p></p>
<p class="MsoNormal"> AT (0,0,0) ABSOLUTE<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">COMPONENT isis_mod = ISIS_moderator(<o:p></o:p></p>
<p class="MsoNormal"> Face = "crisp", dist = 6.3, focus_xw = 0.1, focus_yh = 0.1,<o:p></o:p></p>
<p class="MsoNormal"> xwidth = -1.0, yheight = -1.0, SAC = 1, Lmin = 0.25,<o:p></o:p></p>
<p class="MsoNormal"> Lmax = 20.0, verbose = 1)<o:p></o:p></p>
<p class="MsoNormal"> AT (0.0, 0.0, 0.0) RELATIVE Origin<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">COMPONENT chopper1 = DiskChopper(<o:p></o:p></p>
<p class="MsoNormal"> theta_0 = 77.0, radius = 0.4, yheight = 0.11, nu = 50.0,<o:p></o:p></p>
<p class="MsoNormal"> nslit = 1, jitter = 0.0, delay = 10000.0e-6, verbose = 1)<o:p></o:p></p>
<p class="MsoNormal"> AT (0.0, 0.0, 6.3) RELATIVE Origin<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">COMPONENT TOFmon1 = TOF_monitor(<o:p></o:p></p>
<p class="MsoNormal"> nt = 1000, filename = "TOFmon1.dat", xwidth = 0.15,<o:p></o:p></p>
<p class="MsoNormal"> yheight = 0.15, tmin = 0.0, tmax = 60000.0,<o:p></o:p></p>
<p class="MsoNormal"> restore_neutron = 1)<o:p></o:p></p>
<p class="MsoNormal"> AT (0.0, 0.0, 6.4) RELATIVE Origin<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">FINALLY<o:p></o:p></p>
<p class="MsoNormal">%{<o:p></o:p></p>
<p class="MsoNormal">%}<o:p></o:p></p>
<p class="MsoNormal">END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The error seems to be in the jitter part of the code.<o:p></o:p></p>
<p class="MsoNormal">This mod works at line 157 of diskchopper.comp for me.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> jitter=jitter*randnorm();<o:p></o:p></p>
<p class="MsoNormal"> //toff=fabs(t-atan2(x,yprime)/omega - delay - jitter ? jitter*randnorm():0);<o:p></o:p></p>
<p class="MsoNormal"> toff=fabs(t-atan2(x,yprime)/omega - delay - jitter );<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">===================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">Dr. Robert Dalgliesh<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">ISIS<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">STFC Rutherford Appleton Laboratory<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">R3 1.29<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">Harwell Oxford<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">Didcot<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">OX11 0QX<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">Office Tel: +44 (0) 1235 445687 or OffSpec Instrument: +44 (0) 1235 567045<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">Internal Mobile Extension: 1176<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black;mso-fareast-language:EN-GB">e-mail: <a href="mailto:robert.dalgliesh@stfc.ac.uk"><span style="color:blue">robert.dalgliesh@stfc.ac.uk</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>