[mcstas-users] elliptic mirror problem

Mon Sep 17 13:15:08 CEST 2018

Dear Ekatarina,

I have had a quick look at your instrument, which describes the attached geometry (top view) when visualised with mcdisplay.

Is the large horizontal offset between the two guides intentional?

Could you supply some more information about what you are trying to simulate?

Best,

Peter Willendrup

[cid:6BF1FF73-E488-4562-91BE-32509504C243 at win.dtu.dk]
On 16 Sep 2018, at 21.56, Екатерина Рувинская <ruvinskaia.e.e at mail.ru<mailto:ruvinskaia.e.e at mail.ru>> wrote:

Dear all,

I have a problem. I need to neutron beam pass through elliptical mirrors system. I tried approach through component "Guide_tapering" and component "Mirror_Elliptic". But in both cases neutron beam does not move after the mirror. Could you say what the source of the problem and how I can fix this. Code with "Guide_tapering" for example are attached.

Thanking in advance.
Yours faithfully, Ekaterina.

DEFINE INSTRUMENT NSE_res(double slit_w = 0.1)

DECLARE
%{

double slit_h;
int iter;
double arm_lenght = 5;
%}

INITIALIZE
%{

%}

TRACE

COMPONENT Origin = Progress_bar()
  AT (0,0,0) ABSOLUTE

COMPONENT source =  Source_gen(
    yheight=0.00001,xwidth=slit_w, lambda0=9, dlambda=0.000000000000001,I1=1e13,
    verbose=1,focus_xw=slit_w,focus_yh=0.0001, dist = arm_lenght/4)
    AT (0, 0, 0)  RELATIVE Origin
  EXTEND
  %{
     t = rand0max(1e-20); // set time from 0 to 1 ms for TOF instruments.
  %}

COMPONENT mon_TOF_source = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits=[0.000000001 0.3], bins=180, lambda, limits=[8.9 9.1], bins=180")
    AT (0, 0, 0.00001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_source = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "lambda, limits = [8.99, 9.01], bins = 180")
    AT (0, 0, 0.00001) RELATIVE PREVIOUS

COMPONENT first_arm = Arm()
    AT (0, 0, 0.000001) RELATIVE PREVIOUS

    COMPONENT slit_first = Slit(xwidth = slit_w, yheight = 0.00001)
        AT (0, 0, 1E-4) RELATIVE first_arm

    COMPONENT elliptical_1 = Guide_tapering(
        option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
        linw = 0.5, loutw = 0.5,  linh = 0.5, louth = 0.5,
         segno = 100, R0 = 0.995,
        Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
        alphay = 6.07, W = 0.003, mx = 3, my = 3)
        AT (-0.5+slit_w/2, 0,1E-3) RELATIVE first_arm
        EXTEND
        %{
           if(!SCATTERED) ABSORB;
        %}

COMPONENT second_arm = Arm()
    AT (0, 0, arm_lenght) RELATIVE first_arm

    COMPONENT elliptical_2 = Guide_tapering(
        option = "elliptical", w1 = 1, h1 = 0.5, l = arm_lenght,
        linw = 0.5, loutw = 0.5,  linh = 0.5, louth = 0.5,
         segno = 10, R0 = 0.995,
        Qcx = 0.0217, Qcy = 0.0217, alphax = 6.07,
        alphay = 6.07, W = 0.003, mx = 3, my = 3)
        AT (0.5-slit_w/2, 0,1E-3) RELATIVE second_arm
        EXTEND
        %{
           if(!SCATTERED) ABSORB;
        %}

    COMPONENT slit_f = Slit(xwidth = slit_w, yheight = 0.00001)
        AT (0, 0, arm_lenght) RELATIVE second_arm

COMPONENT mon_TOF_end = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time bins=180, lambda, limits=[29 31], bins=180")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_lambda_end = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "lambda, limits = [8.99, 9.01], bins = 180")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_time_end_9 = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits = [0.022784, 0.022792] bins = 200")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

COMPONENT monitor_time_end_30 = Monitor_nD(
    xwidth = 0.4, yheight = 0.4, restore_neutron = 1,
    options = "time, limits = [0.0455305, 0.0455307] bins = 200")
    AT (0, 0, 0.001) RELATIVE PREVIOUS

FINALLY
%{
%}

END

_______________________________________________
mcstas-users mailing list
mcstas-users at mcstas.org<mailto:mcstas-users at mcstas.org>
https://mailman2.mcstas.org/mailman/listinfo/mcstas-users

Peter Kjær Willendrup
Forskningsingeniør, Speciakonsulent
Næstformand for DTU Fysik LSU

DTU Physics

[cid:b6339554-7f28-441d-9f31-5944f811397c at win.dtu.dk]

Technical University of Denmark

[cid:e108e480-fcab-46e2-9531-b38165079572 at win.dtu.dk]

Department of Physics
Fysikvej
Building 307
DK-2800 Kongens Lyngby
Direct +45 2125 4612
Mobil +45 2125 4612
Fax +45 4593 2399
pkwi at fysik.dtu.dk<mailto:pkwi at fysik.dtu.dk>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2018-09-17 at 13.11.26.png
Type: image/png
Size: 66203 bytes
Desc: Screen Shot 2018-09-17 at 13.11.26.png
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.gif
Type: image/gif
Size: 58 bytes
Desc: image001.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0002.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.gif
Type: image/gif
Size: 1055 bytes
Desc: image002.gif
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20180917/73c7bbf4/attachment-0003.gif>