<div dir="ltr">Hello again Tom,<br><br>In taking a closer look at your instrument, I would personally be focusing to slightly in front of the M2Exit slit.  My preference would be the following:<br><br>xw=0.08m , yh= 0.13m, dist=3.9m<br>
<br>This will assure that the focusing parameters do not influence the results of the simulation.  Looking at your project 3, these parameters are probably fine.  However, in project 4 you are probably under illuminated the collimator and definitely under illuminating the guide.<br>
<br>Now, to continue on with your simulation, I would suggest using the virtual_input after collimator1 along with the same assortment of detectors.  Run a simulation such that you have a good distribution of divergence and histories on your detectors.  You can then continue the simulation form the point after the collimator1 using the virtual_output.  You will find that the repeat parameter will come in handy.  Remember to use the SPLIT 10 keyword on the powderN component.<br>
<br>Hope this helps, feel free to ask further questions.<br><br>Aaron<br><br><div class="gmail_quote">On Wed, Jul 30, 2008 at 4:14 AM, Peter Willendrup <span dir="ltr"><<a href="mailto:peter.willendrup@risoe.dk">peter.willendrup@risoe.dk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div style="">Hello Tom,<div><br><div><div class="Ih2E3d"><div>On Jul 30, 2008, at 6:17 AM, wokaoyan1981 wrote:</div>
<blockquote type="cite"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<div><div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">Thank you for your wonderful suggestions. I delete components behind collimator 1 and set a series of<span> </span><span> </span>monitors instead . Parameters dist, xw and yh are set as the following:</span></div>
<div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">1: xw=0.06m</span><span>¡¢</span><span lang="EN-US">yh=0.11m</span><span>¡¢</span><span lang="EN-US">dist=2.000m</span><span>£¨</span><span lang="EN-US">prototype</span><span>£©</span><span lang="EN-US">;</span><span></span></div>
<div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">2: xw=0.08m</span><span>¡¢</span><span lang="EN-US">yh=0.13m</span><span>¡¢</span><span lang="EN-US">dist=2.000m;</span></div>
<div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">3: xw=0.08m</span><span>¡¢</span><span lang="EN-US">yh=0.13m</span><span>¡¢</span><span lang="EN-US">dist=4.049m<span> </span>(focusing to the M2Exit slit)</span><span>£»</span></div>
<div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">4: xw=0.08m</span><span>¡¢</span><span lang="EN-US">yh=0.13m</span><span>¡¢</span><span lang="EN-US">dist=10.649m<span> </span>(focusing to the end of the collimator1)</span><span>£»</span></div>
<div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span><br></span></div><div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;">
<span lang="EN-US">The simulation suggests: N increases remarkably while ERR decreases.<span> </span></span></div></div></div></span></blockquote><div><br></div></div><div>Yes, this is indeed the case. Consider a given McStas detector bin. N neutron rays (not actual neutrons but neutron "rays",   cf. chapter 4 in the manual) hit this bin.</div>
<div>The N neutron rays each have a given weight p_j, from which we determine the intensity (in units of neutrons/second) by Sum(p_j). ERR is nothing but the RMS error in units of neutrons/second, given that N rays of individual weight p_j hit the bin. So yes, ERR will decrease with increased N.</div>
<div><br></div><div>In summary, I and ERR is your intensity with errorbars while N describes your statistics (and hence defines the errorbars).</div><div class="Ih2E3d"><div><br></div><blockquote type="cite"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<div><div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><br></div><div style="margin: 0cm 0cm 0.0001pt; text-indent: 15.75pt; text-align: justify; font-size: 10.5pt;"><span lang="EN-US">I don't know which project to choose since the results differ in neutron intenstiy, divergence and distribution. See I, y_div, x_I in enclosure.</span></div>
</div></div></span></blockquote></div></div></div><div><br></div><div>I did not look at your instrument file in detail, so exuse me if I get something wrong. The table above suggests that you are experimenting with focusing from the source? Source focusing is only meant to reduce the solid angle of the produced neutrons, ie. it does not make sense to emit neutrons in full 4PI (as a real world neutron source) when we are only interested in neutrons that hit our guide entrance.</div>
<div><br></div><div>Beware however; If you reduce the solid angle so much that eg. the first guide element is under illuminated you WILL get faulty results, in terms of underestimated divergence and wrong intensity.</div>
<div><br></div><div>Will have a real look at your instrument file as soon as time allows.</div><div><br></div><div><br></div><div>Hope this helps,</div><div><br></div><font color="#888888"><div>Peter Willendrup</div><div>
<br></div></font></div><br>_______________________________________________<br>
mcstas-users mailing list<br>
<a href="mailto:mcstas-users@mcstas.org">mcstas-users@mcstas.org</a><br>
<a href="http://mailman.mcstas.org/cgi-bin/mailman/listinfo/mcstas-users" target="_blank">http://mailman.mcstas.org/cgi-bin/mailman/listinfo/mcstas-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>*************************************************<br>Aaron M. Percival<br>M.Sc. Candidate<br>Dept. of Physics, Engineering Physics & Astronomy<br>Queen's University<br>
Kingston, Ontario, Canada, K7L 3N6<br>Office: 613-533-6000 ext. 74789<br>Fax: 613-533-6463<br>*************************************************<br>
</div>