From farhi at ill.fr  Tue Jan  9 23:08:37 2007
From: farhi at ill.fr (Emmanuel FARHI)
Date: Tue, 09 Jan 2007 23:08:37 +0100 (CET)
Subject: [neutron-mc] ILL position: MCNPX permanent position
Message-ID: <58455.193.49.43.124.1168380517.squirrel@195.83.126.3>

Dear MC simulators,

The ILL has opened a permanent position for a physicist-engineer in order
to work with MCNPX. Deadline is Jan 21st 2007.
Details are to be found at
<http://www9.profilsearch.com/recrute/voir_offre.asp?diffusion_id=165&mission_id=47&Lang=GB>
or <http://www.ill.fr> in Jobs/Join our team/Job vacancies.

You are welcome to apply !

E. Farhi.

-- 
FARHI Emmanuel <farhi at ill.fr>
Groupe DS/CS, ILL4/156, Tel 04 76 20 71 35
ILL, Grenoble



From percival at physics.queensu.ca  Wed Jan 10 16:27:44 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Wed, 10 Jan 2007 10:27:44 -0500
Subject: [neutron-mc] Simple Question
Message-ID: <006801c734cb$e2f941a0$a21e0f82@amgap>

Good Day,
 
I have a simple question about exporting McStas data from Matlab to another
platform.  I'll give an example of what I want to do:
 
I have a type Monitor_nD detector with the following options set, options =
"banana, theta, limits=[60,150], bins=100".
 
I'd like to take the data from this detector and export it to, say Origin
for further analysis.  I've found in the detector file the array for the
data and the errors, however, I cannot find the array for the x-values (in
this case theta values).  However, I did find something called stepx:
0.9091.  Are the x-values generated by starting at (in this case) 60, then
increasing each increment by 0.9091, 100 times (since there are 100 bins)?
 
Regards,
 
Aaron M. Percival
 
*************************************************
Aaron M. Percival
M.Sc. Candidate
Dept. of Physics, Engineering Physics & Astronomy
Queen's University
Kingston, Ontario, Canada, K7L 3N6
Office: 613-533-6000 ext. 74789
Fax: 613-533-6463
*************************************************
 
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070110/6aa4a9d7/attachment.html>

From farhi at ill.fr  Wed Jan 10 16:48:15 2007
From: farhi at ill.fr (Emmanuel FARHI)
Date: Wed, 10 Jan 2007 16:48:15 +0100
Subject: [neutron-mc] Simple Question
In-Reply-To: <006801c734cb$e2f941a0$a21e0f82@amgap>
References: <006801c734cb$e2f941a0$a21e0f82@amgap>
Message-ID: <45A50ABF.7090703@ill.fr>

An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070110/5fd05be3/attachment.html>

From percival at physics.queensu.ca  Wed Jan 10 19:40:16 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Wed, 10 Jan 2007 13:40:16 -0500
Subject: [neutron-mc] Simple Question
In-Reply-To: <45A50ABF.7090703@ill.fr>
Message-ID: <009901c734e6$c8b50f20$a21e0f82@amgap>

Emmanuel,
 
Your suggestion worked.  I was able to convert the data file from Matlab
format to the McStas format (which seems to be simple text) using
mcformat.exe.  I then used excel to import only the data columns (in this
case there was theta, counts, error and intensity).  I took this data from
excel into Origin, and used Origin to perform analysis.  I do notice in
doing this that the Gaussian fit found in Origin gives slightly different
values for X0 and dX (the variable McStas used in Matlab) for the centre and
width of the curve respectively.  The actual numbers are X0 = 99.852 dX =
6.024 (from Matlab) and X0 = 97.3, dX = 7.7, which as you can see is quite
significant.
 
Another interesting point to note is that the curve produced in Origin is
slightly different than the one produced in Matlab by mcplot, but they
should be produced from the same data.
 
Aaron.
 
  _____  

From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On
Behalf Of Emmanuel FARHI
Sent: Wednesday, January 10, 2007 10:48 AM
To: McStas users list at neutron.risoe.dk
Subject: Re: [neutron-mc] Simple Question
 
Hi AAron,

The idea is that axes are always regular, so that the min, max and number of
steps is sufficient to generate the axis vector. The min and max is given in
the 'xylimits' field of Matlab data structure, or on '# xylimits:' line of
PGPLOT data format file. The number of steps is written in the 'type' field
and '# type:' line as well. Then that's really easy to generate the axes
with 'linspace' function in Matlab.

the best for you is to generate the data using the PGPLOT format. If you do
not have PGPLOT you will not be able to plot it. Anyway you may use the
'mcformat' tool to convert it to e.g. Matlab format. You may also convert -
still with mcformat - your Matlab data into PGPLOT if the simulation is too
long to re-run.

E.

Aaron M. Percival wrote: 
Good Day,
 
I have a simple question about exporting McStas data from Matlab to another
platform.  I'll give an example of what I want to do:
 
I have a type Monitor_nD detector with the following options set, options =
"banana, theta, limits=[60,150], bins=100".
 
I'd like to take the data from this detector and export it to, say Origin
for further analysis.  I've found in the detector file the array for the
data and the errors, however, I cannot find the array for the x-values (in
this case theta values).  However, I did find something called stepx:
0.9091.  Are the x-values generated by starting at (in this case) 60, then
increasing each increment by 0.9091, 100 times (since there are 100 bins)?
 
Regards,
 
Aaron M. Percival
 
*************************************************
Aaron M. Percival
M.Sc. Candidate
Dept. of Physics, Engineering Physics & Astronomy
Queen's University
Kingston, Ontario, Canada, K7L 3N6
Office: 613-533-6000 ext. 74789
Fax: 613-533-6463
*************************************************
 
 
 





  _____  



 
_______________________________________________
neutron-mc mailing list
neutron-mc at risoe.dk
http://mailman.risoe.dk/mailman/listinfo/neutron-mc
  



-- 
What's up Doc ?
--------------------------------------------
Emmanuel FARHI, http://www.ill.fr/Computing/people/Farhi  \|/ ____ \|/
CS-Group BP2/22,   Institut Laue-Langevin (ILL) Grenoble  ~@-/ oO \-@~
6 rue J. Horowitz, BP 156, 38042 Grenoble Cedex 9,France  /_( \__/ )_\
Work :Tel (33/0) 4 76 20 71 35. Fax (33/0) 4 76 48 39 06     \__U_/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070110/fce2c2e8/attachment.html>

From farhi at ill.fr  Wed Jan 10 21:54:28 2007
From: farhi at ill.fr (Emmanuel FARHI)
Date: Wed, 10 Jan 2007 21:54:28 +0100 (CET)
Subject: [neutron-mc] Simple Question
In-Reply-To: <009901c734e6$c8b50f20$a21e0f82@amgap>
References: <45A50ABF.7090703@ill.fr> <009901c734e6$c8b50f20$a21e0f82@amgap>
Message-ID: <45775.193.49.43.124.1168462468.squirrel@195.83.126.3>

Hi again Aaron,

the difference comes from the way the axis is generated. Consider that the
data is stored into bins which have a given extent dx. The first left bin
value is min, the right last bin value is max. This means that the axis
values should better be the center of the bins for data analysis and
plotting. That's the way McStas plot the data. So you should generate the
bins between  min+dx/2 and max-dx/2, with dx=(max-min)/bins.

E.



> Emmanuel,
>
> Your suggestion worked.  I was able to convert the data file from Matlab
> format to the McStas format (which seems to be simple text) using
> mcformat.exe.  I then used excel to import only the data columns (in this
> case there was theta, counts, error and intensity).  I took this data from
> excel into Origin, and used Origin to perform analysis.  I do notice in
> doing this that the Gaussian fit found in Origin gives slightly different
> values for X0 and dX (the variable McStas used in Matlab) for the centre
> and
> width of the curve respectively.  The actual numbers are X0 = 99.852 dX =
> 6.024 (from Matlab) and X0 = 97.3, dX = 7.7, which as you can see is quite
> significant.
>
> Another interesting point to note is that the curve produced in Origin is
> slightly different than the one produced in Matlab by mcplot, but they
> should be produced from the same data.
>
> Aaron.
>
>   _____
>
> From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On
> Behalf Of Emmanuel FARHI
> Sent: Wednesday, January 10, 2007 10:48 AM
> To: McStas users list at neutron.risoe.dk
> Subject: Re: [neutron-mc] Simple Question
>
> Hi AAron,
>
> The idea is that axes are always regular, so that the min, max and number
> of
> steps is sufficient to generate the axis vector. The min and max is given
> in
> the 'xylimits' field of Matlab data structure, or on '# xylimits:' line of
> PGPLOT data format file. The number of steps is written in the 'type'
> field
> and '# type:' line as well. Then that's really easy to generate the axes
> with 'linspace' function in Matlab.
>
> the best for you is to generate the data using the PGPLOT format. If you
> do
> not have PGPLOT you will not be able to plot it. Anyway you may use the
> 'mcformat' tool to convert it to e.g. Matlab format. You may also convert
> -
> still with mcformat - your Matlab data into PGPLOT if the simulation is
> too
> long to re-run.
>
> E.
>
> Aaron M. Percival wrote:
> Good Day,
>
> I have a simple question about exporting McStas data from Matlab to
> another
> platform.  I'll give an example of what I want to do:
>
> I have a type Monitor_nD detector with the following options set, options
> =
> "banana, theta, limits=[60,150], bins=100".
>
> I'd like to take the data from this detector and export it to, say Origin
> for further analysis.  I've found in the detector file the array for the
> data and the errors, however, I cannot find the array for the x-values (in
> this case theta values).  However, I did find something called stepx:
> 0.9091.  Are the x-values generated by starting at (in this case) 60, then
> increasing each increment by 0.9091, 100 times (since there are 100 bins)?
>
> Regards,
>
> Aaron M. Percival
>
> *************************************************
> Aaron M. Percival
> M.Sc. Candidate
> Dept. of Physics, Engineering Physics & Astronomy
> Queen's University
> Kingston, Ontario, Canada, K7L 3N6
> Office: 613-533-6000 ext. 74789
> Fax: 613-533-6463
> *************************************************
>
>
>
>
>
>
>
>
>   _____
>
>
>
>
> _______________________________________________
> neutron-mc mailing list
> neutron-mc at risoe.dk
> http://mailman.risoe.dk/mailman/listinfo/neutron-mc
>
>
>
>
> --
> What's up Doc ?
> --------------------------------------------
> Emmanuel FARHI, http://www.ill.fr/Computing/people/Farhi  \|/ ____ \|/
> CS-Group BP2/22,   Institut Laue-Langevin (ILL) Grenoble  ~@-/ oO \-@~
> 6 rue J. Horowitz, BP 156, 38042 Grenoble Cedex 9,France  /_( \__/ )_\
> Work :Tel (33/0) 4 76 20 71 35. Fax (33/0) 4 76 48 39 06     \__U_/
> _______________________________________________
> neutron-mc mailing list
> neutron-mc at risoe.dk
> http://mailman.risoe.dk/mailman/listinfo/neutron-mc
>


-- 
FARHI Emmanuel <farhi at ill.fr>
Groupe DS/CS, ILL4/156, Tel 04 76 20 71 35
ILL, Grenoble



From percival at physics.queensu.ca  Wed Jan 10 22:19:27 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Wed, 10 Jan 2007 16:19:27 -0500
Subject: [neutron-mc] Simple Question
In-Reply-To: <45775.193.49.43.124.1168462468.squirrel@195.83.126.3>
Message-ID: <00ad01c734fd$067a5840$a21e0f82@amgap>

Emmanuel,

Your explanation is clear and is certainly the logical way to set up the
bins.  I compared the formula that you gave min+dx/2 and max-dx/2, with
dx=(max-min)/bins, to the numbers that the data file had in it (the one that
I converted from Matlab format to text using mcformat) and found the x-value
matched what they should be (i.e. the centre of each bin).  My concern is
that the graphical data from Matlab seems a little off.  As I mentioned
before the Gaussian centre and HWHM do not agree with the values that I
found in Origin using the same data.  Also, and maybe more importantly, the
plots generated by Matlab seem 'skewed' from the ones generated by Origin.
It's difficult to explain what I mean by skewed, but the plot generated in
Matlab look more Maxwellian than Gaussian (i.e. sharper on one side), the
ones generated from Origin look almost Gaussian.

On the first point, my guess would be that the fitting functions used are
different, in that case, Origin's fit give the centre of the Gaussian where
I expect it to be, Matlab is about 2.5 degrees off.

As for the second point, I do not know why the plots themselves would have
different shapes.  This seems very strange to me.

Aaron

-----Original Message-----
From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On
Behalf Of Emmanuel FARHI
Sent: Wednesday, January 10, 2007 3:54 PM
To: McStas users list at neutron.risoe.dk
Subject: RE: [neutron-mc] Simple Question

Hi again Aaron,

the difference comes from the way the axis is generated. Consider that the
data is stored into bins which have a given extent dx. The first left bin
value is min, the right last bin value is max. This means that the axis
values should better be the center of the bins for data analysis and
plotting. That's the way McStas plot the data. So you should generate the
bins between  min+dx/2 and max-dx/2, with dx=(max-min)/bins.

E.



> Emmanuel,
>
> Your suggestion worked.  I was able to convert the data file from Matlab
> format to the McStas format (which seems to be simple text) using
> mcformat.exe.  I then used excel to import only the data columns (in this
> case there was theta, counts, error and intensity).  I took this data from
> excel into Origin, and used Origin to perform analysis.  I do notice in
> doing this that the Gaussian fit found in Origin gives slightly different
> values for X0 and dX (the variable McStas used in Matlab) for the centre
> and
> width of the curve respectively.  The actual numbers are X0 = 99.852 dX =
> 6.024 (from Matlab) and X0 = 97.3, dX = 7.7, which as you can see is quite
> significant.
>
> Another interesting point to note is that the curve produced in Origin is
> slightly different than the one produced in Matlab by mcplot, but they
> should be produced from the same data.
>
> Aaron.
>
>   _____
>
> From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On
> Behalf Of Emmanuel FARHI
> Sent: Wednesday, January 10, 2007 10:48 AM
> To: McStas users list at neutron.risoe.dk
> Subject: Re: [neutron-mc] Simple Question
>
> Hi AAron,
>
> The idea is that axes are always regular, so that the min, max and number
> of
> steps is sufficient to generate the axis vector. The min and max is given
> in
> the 'xylimits' field of Matlab data structure, or on '# xylimits:' line of
> PGPLOT data format file. The number of steps is written in the 'type'
> field
> and '# type:' line as well. Then that's really easy to generate the axes
> with 'linspace' function in Matlab.
>
> the best for you is to generate the data using the PGPLOT format. If you
> do
> not have PGPLOT you will not be able to plot it. Anyway you may use the
> 'mcformat' tool to convert it to e.g. Matlab format. You may also convert
> -
> still with mcformat - your Matlab data into PGPLOT if the simulation is
> too
> long to re-run.
>
> E.
>
> Aaron M. Percival wrote:
> Good Day,
>
> I have a simple question about exporting McStas data from Matlab to
> another
> platform.  I'll give an example of what I want to do:
>
> I have a type Monitor_nD detector with the following options set, options
> =
> "banana, theta, limits=[60,150], bins=100".
>
> I'd like to take the data from this detector and export it to, say Origin
> for further analysis.  I've found in the detector file the array for the
> data and the errors, however, I cannot find the array for the x-values (in
> this case theta values).  However, I did find something called stepx:
> 0.9091.  Are the x-values generated by starting at (in this case) 60, then
> increasing each increment by 0.9091, 100 times (since there are 100 bins)?
>
> Regards,
>
> Aaron M. Percival
>
> *************************************************
> Aaron M. Percival
> M.Sc. Candidate
> Dept. of Physics, Engineering Physics & Astronomy
> Queen's University
> Kingston, Ontario, Canada, K7L 3N6
> Office: 613-533-6000 ext. 74789
> Fax: 613-533-6463
> *************************************************
>
>
>
>
>
>
>
>
>   _____
>
>
>
>
> _______________________________________________
> neutron-mc mailing list
> neutron-mc at risoe.dk
> http://mailman.risoe.dk/mailman/listinfo/neutron-mc
>
>
>
>
> --
> What's up Doc ?
> --------------------------------------------
> Emmanuel FARHI, http://www.ill.fr/Computing/people/Farhi  \|/ ____ \|/
> CS-Group BP2/22,   Institut Laue-Langevin (ILL) Grenoble  ~@-/ oO \-@~
> 6 rue J. Horowitz, BP 156, 38042 Grenoble Cedex 9,France  /_( \__/ )_\
> Work :Tel (33/0) 4 76 20 71 35. Fax (33/0) 4 76 48 39 06     \__U_/
> _______________________________________________
> neutron-mc mailing list
> neutron-mc at risoe.dk
> http://mailman.risoe.dk/mailman/listinfo/neutron-mc
>


-- 
FARHI Emmanuel <farhi at ill.fr>
Groupe DS/CS, ILL4/156, Tel 04 76 20 71 35
ILL, Grenoble

_______________________________________________
neutron-mc mailing list
neutron-mc at risoe.dk
http://mailman.risoe.dk/mailman/listinfo/neutron-mc





From farhi at ill.fr  Thu Jan 11 11:18:44 2007
From: farhi at ill.fr (Emmanuel FARHI)
Date: Thu, 11 Jan 2007 11:18:44 +0100
Subject: [neutron-mc] Simple Question
In-Reply-To: <00ad01c734fd$067a5840$a21e0f82@amgap>
References: <00ad01c734fd$067a5840$a21e0f82@amgap>
Message-ID: <45A60F04.8060302@ill.fr>

An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070111/1fc4f739/attachment.html>

From percival at physics.queensu.ca  Tue Jan 23 17:23:26 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Tue, 23 Jan 2007 11:23:26 -0500
Subject: [neutron-mc] Trouble with Single_Crystal
Message-ID: <000001c73f0a$d34f7050$a21e0f82@amgap>

Hello all,
 
I'm having a little bit of trouble with the single_crystal component.  I've
created a .lau file in Crystallographica to simulated an Si analyzer
crystal.  I simply want to reflect a neutron beam of wavelength 2.866 A
using the 311 plane (d = 1.637) of Si.  So I've set up the crystal to
satisfy the Bragg condition, giving theta = 61.823, and thus 2theta =
123.646.  I put a banana shaped detector around this 2theta, and I get no
reflections.
 
Here is the crystal:
 
COMPONENT analyzer = Single_crystal(
    reflections = "Si.lau", xwidth = 0.1, yheight = 0.1,
    zthick = 0.01, delta_d_d = 0, mosaic = 18, order = 0,
    p_transmit = 0, absorption = 0, incoherent = 0, ax = 5.430530,
    by = 5.430530, cz = 5.430530, aa = 90, bb = 90, cc = 90)
  AT (0, 0, 0) RELATIVE a1
  ROTATED ( 0, theta311, 0) RELATIVE sample_arm
 
Can any one offer any suggestions as to why this might not be working?
I've also attached the Si.lau file.
 
Thanks, 
 
Aaron
 
*************************************************
Aaron M. Percival
M.Sc. Candidate
Dept. of Physics, Engineering Physics & Astronomy
Queen's University
Kingston, Ontario, Canada, K7L 3N6
Office: 613-533-6000 ext. 74789
Fax: 613-533-6463
*************************************************
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070123/f4b8328d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Si.lau
Type: application/octet-stream
Size: 30426 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070123/f4b8328d/attachment.obj>

From percival at physics.queensu.ca  Tue Jan 23 19:16:22 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Tue, 23 Jan 2007 13:16:22 -0500
Subject: [neutron-mc] Trouble with Single_Crystal
Message-ID: <000b01c73f1a$993a1a90$a21e0f82@amgap>

Hello all,
 
I'm having a little bit of trouble with the single_crystal component.  I've
created a .lau file in Crystallographica to simulate an Si analyzer crystal.
I simply want to reflect a neutron beam of wavelength 2.866 A using the 311
plane (d = 1.637) of Si.  So I've set up the crystal to satisfy the Bragg
condition, giving theta = 61.823, and thus 2theta = 123.646.  I put a banana
shaped detector around this 2theta, and I get no reflections.
 
Here is the crystal:
 
COMPONENT analyzer = Single_crystal(
    reflections = "Si.lau", xwidth = 0.1, yheight = 0.1,
    zthick = 0.01, delta_d_d = 0, mosaic = 18, order = 0,
    p_transmit = 0, absorption = 0, incoherent = 0, ax = 5.430530,
    by = 5.430530, cz = 5.430530, aa = 90, bb = 90, cc = 90)
  AT (0, 0, 0) RELATIVE a1
  ROTATED ( 0, theta311, 0) RELATIVE sample_arm
 
Where theta311 = 61.823 relative to the incident beam.
 
Can any one offer any suggestions as to why this might not be working?
I've also attached the Si.lau file.
 
Thanks, 
 
Aaron
 
 
*************************************************
Aaron M. Percival
M.Sc. Candidate
Dept. of Physics, Engineering Physics & Astronomy
Queen's University
Kingston, Ontario, Canada, K7L 3N6
Office: 613-533-6000 ext. 74789
Fax: 613-533-6463
*************************************************
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070123/60254726/attachment.html>

From fliang at ciae.ac.cn  Thu Jan 25 08:30:41 2007
From: fliang at ciae.ac.cn (Feng Liang)
Date: Thu, 25 Jan 2007 15:30:41 +0800
Subject: [neutron-mc] Could you please give me a hand?
Message-ID: <006b01c74052$b7760710$39410a0a@6311fc60030c48d>

Hallo! Thank you for reading my message. I'm a student form Neutron Scattering Laboratory in China Institute of Atomic Energy, CHN. I have konwn the great M-C program MCSTAS for a long time. 

Now, we will design and install some neutron scattering instruments at our new reactor CARR, but I still have some problems about runing the MCSTAS under the windows system. And I have no experience about using the program. 

So could you please send me a example of a simple SANS instrument simulation just including the cold neutron source, neutron guide, collimation guide, sample station and the dector? I am so appreciate about that help!

Thank you very much! Looking forward your help.

Best regards

Feng Liang
Neutron Scattering Laboratory
China Institute of Atomic Energy
   
  Beijing P.O. Box 275(30)
  102413 Beijing
  China
Tel:+86-10-69357727
Fax:+86-10-69357787
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070125/f991954f/attachment.html>

From peter.willendrup at risoe.dk  Thu Jan 25 10:46:36 2007
From: peter.willendrup at risoe.dk (peter.willendrup at risoe.dk)
Date: Thu, 25 Jan 2007 10:46:36 +0100
Subject: [neutron-mc] Trouble with Single_Crystal
In-Reply-To: <000b01c73f1a$993a1a90$a21e0f82@amgap>
Message-ID: <6463F256A85DC14CAB07F087A2479972282203@EXCHG-VS1.risoe.dk>

Hi Aaron,

 

 

I have not had time to fully investigate your problem. From an earlier attempt to use Single_crystal it seems to me that we might indeed have a problem concerning definition of the internal coordinate system.

 

Will get back to you (and this mailinglist) when I have looked more carefully.

 

Cheers,

 

Peter Willendrup

 

-- 

-----------------------------------------------

Peter Willendrup, MSc

Materials Research Department

Building 227

Ris? National Laboratory

Technical University of Denmark

Frederiksborgvej 399

DK-4000 Roskilde

Denmark

 

Tel:  (+45) 4677 5862

Fax:  (+45) 4677 5758

Email: Peter.Willendrup at risoe.dk

 

-----------------------------------------------

 

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for

Food and Veterinary Research, the Danish Institute for Fisheries

Research, the Danish National Space Center and the Danish Transport

Research Institute have been merged with the Technical University of

Denmark (DTU) with DTU as the continuing unit. 

 

________________________________

From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On Behalf Of Aaron M. Percival
Sent: Tuesday, January 23, 2007 7:16 PM
To: 'McStas users list at neutron.risoe.dk'
Subject: [neutron-mc] Trouble with Single_Crystal

 

Hello all,

 

I'm having a little bit of trouble with the single_crystal component.  I've created a .lau file in Crystallographica to simulate an Si analyzer crystal.  I simply want to reflect a neutron beam of wavelength 2.866 A using the 311 plane (d = 1.637) of Si.  So I've set up the crystal to satisfy the Bragg condition, giving theta = 61.823, and thus 2theta = 123.646.  I put a banana shaped detector around this 2theta, and I get no reflections.

 

Here is the crystal:

 

COMPONENT analyzer = Single_crystal(

    reflections = "Si.lau", xwidth = 0.1, yheight = 0.1,

    zthick = 0.01, delta_d_d = 0, mosaic = 18, order = 0,

    p_transmit = 0, absorption = 0, incoherent = 0, ax = 5.430530,

    by = 5.430530, cz = 5.430530, aa = 90, bb = 90, cc = 90)

  AT (0, 0, 0) RELATIVE a1

  ROTATED ( 0, theta311, 0) RELATIVE sample_arm

 

Where theta311 = 61.823 relative to the incident beam.

 

Can any one offer any suggestions as to why this might not be working?

I've also attached the Si.lau file.

 

Thanks, 

 

Aaron

 

 

*************************************************

Aaron M. Percival

M.Sc. Candidate

Dept. of Physics, Engineering Physics & Astronomy

Queen's University

Kingston, Ontario, Canada, K7L 3N6

Office: 613-533-6000 ext. 74789

Fax: 613-533-6463

*************************************************

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070125/0faabd80/attachment.html>

From peter.willendrup at risoe.dk  Thu Jan 25 10:56:08 2007
From: peter.willendrup at risoe.dk (peter.willendrup at risoe.dk)
Date: Thu, 25 Jan 2007 10:56:08 +0100
Subject: [neutron-mc] Could you please give me a hand?
In-Reply-To: <006b01c74052$b7760710$39410a0a@6311fc60030c48d>
Message-ID: <6463F256A85DC14CAB07F087A2479972282204@EXCHG-VS1.risoe.dk>

Hello Feng,

 

 

The McStas installation includes a simple SANS instrument file as the one you are looking for. Once installed,

it is available in c:/McStas/lib/examples/SANS.instr.

 

Regarding installation and use of McStas, refer to the instructions given in the installer and the manual. The simplest way to install McStas on Windows is to download and execute

 

http://www.mcstas.org/downloads/McStas-1.10-i686-Win32.exe

 

(Say 'yes' everywhere and click 'next' all the way).

 

Hope this helps,

 

Peter Willendrup

 

-- 

-----------------------------------------------

Peter Willendrup, MSc

Materials Research Department

Building 227

Ris? National Laboratory

Technical University of Denmark

Frederiksborgvej 399

DK-4000 Roskilde

Denmark

 

Tel:  (+45) 4677 5862

Fax:  (+45) 4677 5758

Email: Peter.Willendrup at risoe.dk

 

-----------------------------------------------

 

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for

Food and Veterinary Research, the Danish Institute for Fisheries

Research, the Danish National Space Center and the Danish Transport

Research Institute have been merged with the Technical University of

Denmark (DTU) with DTU as the continuing unit. 

 

________________________________

From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On Behalf Of Feng Liang
Sent: Thursday, January 25, 2007 8:31 AM
To: neutron-mc at risoe.dk
Subject: [neutron-mc] Could you please give me a hand?

 

Hallo! Thank you for reading my message. I'm a student form Neutron Scattering Laboratory in China Institute of Atomic Energy, CHN. I have konwn the great M-C program MCSTAS for a long time. 

 

Now, we will design and install some neutron scattering instruments at our new reactor CARR, but I still have some problems about runing the MCSTAS under the windows system. And I have no experience about using the program. 

 

So could you please send me a example of a simple SANS instrument simulation just including the cold neutron source, neutron guide, collimation guide, sample station and the dector? I am so appreciate about that help!

 

Thank you very much! Looking forward your help.

 

Best regards

 

Feng Liang
Neutron Scattering Laboratory
China Institute of Atomic Energy
   
  Beijing P.O. Box 275(30)
  102413 Beijing
  China
Tel:+86-10-69357727
Fax:+86-10-69357787

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070125/7288cb1d/attachment.html>

From fliang at ciae.ac.cn  Fri Jan 26 14:47:30 2007
From: fliang at ciae.ac.cn (Feng Liang)
Date: Fri, 26 Jan 2007 21:47:30 +0800
Subject: [neutron-mc] Re: neutron-mc Digest, Vol 25, Issue 6
References: <200701250956.l0P9uRNt012462@ftp.risoe.dk>
Message-ID: <003601c74150$87593cb0$39410a0a@6311fc60030c48d>

Thank you very much!

Feng Liang
Neutron Scattering Laboratory
China Institute of Atomic Energy
   
  Beijing P.O. Box 275(30)
  102413 Beijing
  China
Tel:+86-10-69357727
Fax:+86-10-69357787
----- Original Message ----- 
From: <neutron-mc-request at risoe.dk>
To: <neutron-mc at risoe.dk>
Sent: Thursday, January 25, 2007 5:56 PM
Subject: neutron-mc Digest, Vol 25, Issue 6


Send neutron-mc mailing list submissions to
neutron-mc at risoe.dk

To subscribe or unsubscribe via the World Wide Web, visit
http://mailman.risoe.dk/mailman/listinfo/neutron-mc
or, via email, send a message with subject or body 'help' to
neutron-mc-request at risoe.dk

You can reach the person managing the list at
neutron-mc-owner at risoe.dk

When replying, please edit your Subject line so it is more specific
than "Re: Contents of neutron-mc digest..."


Today's Topics:

   1. Could you please give me a hand? (Feng Liang)
   2. RE: Trouble with Single_Crystal (peter.willendrup at risoe.dk)
   3. RE: Could you please give me a hand? (peter.willendrup at risoe.dk)


----------------------------------------------------------------------

Message: 1
Date: Thu, 25 Jan 2007 15:30:41 +0800
From: Feng Liang <fliang at ciae.ac.cn>
Subject: [neutron-mc] Could you please give me a hand?
To: neutron-mc at risoe.dk
Message-ID: <006b01c74052$b7760710$39410a0a at 6311fc60030c48d>
Content-Type: text/plain; charset="gb2312"

Hallo! Thank you for reading my message. I'm a student form Neutron Scattering Laboratory in China Institute of Atomic Energy, CHN. I have konwn the great M-C program MCSTAS for a long time. 

Now, we will design and install some neutron scattering instruments at our new reactor CARR, but I still have some problems about runing the MCSTAS under the windows system. And I have no experience about using the program. 

So could you please send me a example of a simple SANS instrument simulation just including the cold neutron source, neutron guide, collimation guide, sample station and the dector? I am so appreciate about that help!

Thank you very much! Looking forward your help.

Best regards

Feng Liang
Neutron Scattering Laboratory
China Institute of Atomic Energy
   
  Beijing P.O. Box 275(30)
  102413 Beijing
  China
Tel:+86-10-69357727
Fax:+86-10-69357787
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.risoe.dk/pipermail/neutron-mc/attachments/20070125/f991954f/attachment-0001.htm

------------------------------

Message: 2
Date: Thu, 25 Jan 2007 10:46:36 +0100
From: peter.willendrup at risoe.dk
Subject: RE: [neutron-mc] Trouble with Single_Crystal
To: neutron-mc at risoe.dk
Message-ID:
<6463F256A85DC14CAB07F087A2479972282203 at EXCHG-VS1.risoe.dk>
Content-Type: text/plain; charset="iso-8859-1"

Hi Aaron,

 

 

I have not had time to fully investigate your problem. From an earlier attempt to use Single_crystal it seems to me that we might indeed have a problem concerning definition of the internal coordinate system.

 

Will get back to you (and this mailinglist) when I have looked more carefully.

 

Cheers,

 

Peter Willendrup

 

-- 

-----------------------------------------------

Peter Willendrup, MSc

Materials Research Department

Building 227

Ris? National Laboratory

Technical University of Denmark

Frederiksborgvej 399

DK-4000 Roskilde

Denmark

 

Tel:  (+45) 4677 5862

Fax:  (+45) 4677 5758

Email: Peter.Willendrup at risoe.dk

 

-----------------------------------------------

 

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for

Food and Veterinary Research, the Danish Institute for Fisheries

Research, the Danish National Space Center and the Danish Transport

Research Institute have been merged with the Technical University of

Denmark (DTU) with DTU as the continuing unit. 

 

________________________________

From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On Behalf Of Aaron M. Percival
Sent: Tuesday, January 23, 2007 7:16 PM
To: 'McStas users list at neutron.risoe.dk'
Subject: [neutron-mc] Trouble with Single_Crystal

 

Hello all,

 

I'm having a little bit of trouble with the single_crystal component.  I've created a .lau file in Crystallographica to simulate an Si analyzer crystal.  I simply want to reflect a neutron beam of wavelength 2.866 A using the 311 plane (d = 1.637) of Si.  So I've set up the crystal to satisfy the Bragg condition, giving theta = 61.823, and thus 2theta = 123.646.  I put a banana shaped detector around this 2theta, and I get no reflections.

 

Here is the crystal:

 

COMPONENT analyzer = Single_crystal(

    reflections = "Si.lau", xwidth = 0.1, yheight = 0.1,

    zthick = 0.01, delta_d_d = 0, mosaic = 18, order = 0,

    p_transmit = 0, absorption = 0, incoherent = 0, ax = 5.430530,

    by = 5.430530, cz = 5.430530, aa = 90, bb = 90, cc = 90)

  AT (0, 0, 0) RELATIVE a1

  ROTATED ( 0, theta311, 0) RELATIVE sample_arm

 

Where theta311 = 61.823 relative to the incident beam.

 

Can any one offer any suggestions as to why this might not be working?

I've also attached the Si.lau file.

 

Thanks, 

 

Aaron

 

 

*************************************************

Aaron M. Percival

M.Sc. Candidate

Dept. of Physics, Engineering Physics & Astronomy

Queen's University

Kingston, Ontario, Canada, K7L 3N6

Office: 613-533-6000 ext. 74789

Fax: 613-533-6463

*************************************************

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.risoe.dk/pipermail/neutron-mc/attachments/20070125/0faabd80/attachment-0001.htm

------------------------------

Message: 3
Date: Thu, 25 Jan 2007 10:56:08 +0100
From: peter.willendrup at risoe.dk
Subject: RE: [neutron-mc] Could you please give me a hand?
To: neutron-mc at risoe.dk
Message-ID:
<6463F256A85DC14CAB07F087A2479972282204 at EXCHG-VS1.risoe.dk>
Content-Type: text/plain; charset="iso-8859-1"

Hello Feng,

 

 

The McStas installation includes a simple SANS instrument file as the one you are looking for. Once installed,

it is available in c:/McStas/lib/examples/SANS.instr.

 

Regarding installation and use of McStas, refer to the instructions given in the installer and the manual. The simplest way to install McStas on Windows is to download and execute

 

http://www.mcstas.org/downloads/McStas-1.10-i686-Win32.exe

 

(Say 'yes' everywhere and click 'next' all the way).

 

Hope this helps,

 

Peter Willendrup

 

-- 

-----------------------------------------------

Peter Willendrup, MSc

Materials Research Department

Building 227

Ris? National Laboratory

Technical University of Denmark

Frederiksborgvej 399

DK-4000 Roskilde

Denmark

 

Tel:  (+45) 4677 5862

Fax:  (+45) 4677 5758

Email: Peter.Willendrup at risoe.dk

 

-----------------------------------------------

 

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for

Food and Veterinary Research, the Danish Institute for Fisheries

Research, the Danish National Space Center and the Danish Transport

Research Institute have been merged with the Technical University of

Denmark (DTU) with DTU as the continuing unit. 

 

________________________________

From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On Behalf Of Feng Liang
Sent: Thursday, January 25, 2007 8:31 AM
To: neutron-mc at risoe.dk
Subject: [neutron-mc] Could you please give me a hand?

 

Hallo! Thank you for reading my message. I'm a student form Neutron Scattering Laboratory in China Institute of Atomic Energy, CHN. I have konwn the great M-C program MCSTAS for a long time. 

 

Now, we will design and install some neutron scattering instruments at our new reactor CARR, but I still have some problems about runing the MCSTAS under the windows system. And I have no experience about using the program. 

 

So could you please send me a example of a simple SANS instrument simulation just including the cold neutron source, neutron guide, collimation guide, sample station and the dector? I am so appreciate about that help!

 

Thank you very much! Looking forward your help.

 

Best regards

 

Feng Liang
Neutron Scattering Laboratory
China Institute of Atomic Energy
   
  Beijing P.O. Box 275(30)
  102413 Beijing
  China
Tel:+86-10-69357727
Fax:+86-10-69357787

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.risoe.dk/pipermail/neutron-mc/attachments/20070125/7288cb1d/attachment.htm

------------------------------

_______________________________________________
neutron-mc mailing list
neutron-mc at risoe.dk
http://mailman.risoe.dk/mailman/listinfo/neutron-mc


End of neutron-mc Digest, Vol 25, Issue 6
*****************************************



From peter.willendrup at risoe.dk  Thu Feb  1 09:54:17 2007
From: peter.willendrup at risoe.dk (Willendrup, Peter Kj=?ISO-8859-1?B?5g==?=r)
Date: Thu, 01 Feb 2007 09:54:17 +0100
Subject: [neutron-mc] McStas info mail, 2007 february
Message-ID: <C1E76949.DBB%PKWI@exchange.Risoe.dk>

Hello McStas users,


Here is a little information from the McStas team:

---------------------------------------------------------------------------
* Bugs/annoyances in current (mcstas-1.10) release:
---------------------------------------------------------------------------
- DiskChopper.comp: A bug in the shipped component makes it impossible
  to have multiple instances of this component, having different radius.
$ Solution: Get the bug-fixed component from
http://cvs.mcstas.org/cgi-bin/viewcvs.cgi/mcstas/lib/optics/DiskChopper.comp

- Threading: We recommend NOT to use the threading mechanism shipped with
mcstas-1.10 for these reasons: A) Performance is not as good as expected. B)
Threading will not work when writing virtual source files using
Virtual_Output. C) Threading sometimes produces corrupted output when using
mcdisplay.
$ Solution: Install MPI (preferably mpich) on your multicore/multiprocessor
system (Please still take remarks on Virtual_input/output and MPI from the
manual into account).

---------------------------------------------------------------------------
* US visit:
---------------------------------------------------------------------------
- The McStas team recently visited the SNS in Oak Ridge, TN and met with SNS
staff and Jiao Lin from the DANSE
(http://wiki.cacr.caltech.edu/danse/index.php/Main_Page) project. Our time
there was very productive and gave results that the McStas project and user
base will benefit from:
$ The next release will include the possibility to generate NeXus files if
the  NeXus libraries are available at install time. SNS is part of the final
tuning of the NeXus implementation in McStas.
$ The collaboration with DANSE has been strengthened: Standard McStas
releases can now be utilized in DANSE without modification of the core
sourcecode. 

---------------------------------------------------------------------------
* Next release:
---------------------------------------------------------------------------
- The plan is to release mcstas-1.11 late february/early march this year. We
here list a few of the release highlights:
$ Bugfixes, see above.
$ NeXus output if NeXus libs are available at install time
$ On Unix systems (including MacOS X), possibility to 'detach' simulations
from mcgui (start simulations in the background and be ready to continue
work in mcgui)
$ .deb (and possibly .rpm) distribution packages for Linux on Intel
processors will be provided
$ Improvements of support on MacOS X (possibly a .dmg style installer with
all needed support tools)


Best regards,

The McStas team:

Peter Willendrup (Risoe, DTU)
Emmanuel Farhi (ILL)
Kim Lefmann (Risoe, DTU)


--  
-----------------------------------------------
Peter Willendrup, MSc
Materials Research Department
Building 227
Ris? National Laboratory
Technical University of Denmark
Frederiksborgvej 399
DK-4000 Roskilde
Denmark

Tel:  (+45) 4677 5862
Fax:  (+45) 4677 5758
Email: Peter.Willendrup at risoe.dk

-----------------------------------------------

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for
Food and Veterinary Research, the Danish Institute for Fisheries
Research, the Danish National Space Center and the Danish Transport
Research Institute have been merged with the Technical University of
Denmark (DTU) with DTU as the continuing unit.



Sent using the Microsoft Entourage 2004 for Mac Test Drive.




From peter.willendrup at risoe.dk  Fri Feb  9 09:10:47 2007
From: peter.willendrup at risoe.dk (Willendrup, Peter Kj=?ISO-8859-1?B?5g==?=r)
Date: Fri, 09 Feb 2007 09:10:47 +0100
Subject: [neutron-mc] McStas Mac OS X beta installer
Message-ID: <C1F1EB17.EA2%PKWI@exchange.Risoe.dk>

Hello everyone,


Good news for Mac OS X users:

We now have a beta quality installer available for download at:

http://osx.mcstas.org/

(File size ~ 60 Mb)

(Currently tested on OS X 10.4 Intel, should work directly on PPC, may
need tuning to run on earlier releases of Mac OS X.)

The installer is more 'Unix' like than 'Mac' like, but to my opinion
has the required functionality. - See the included README.TXT for details.


I will be very happy to receive your feedback directly at

peter.willendrup at risoe.dk



Regards,

Peter Willendrup


--  
-----------------------------------------------
Peter Willendrup, MSc
Materials Research Department
Building 227
Ris? National Laboratory
Technical University of Denmark
Frederiksborgvej 399
DK-4000 Roskilde
Denmark

Tel:  (+45) 4677 5862
Fax:  (+45) 4677 5758
Email: Peter.Willendrup at risoe.dk

-----------------------------------------------

>From 1 January 2007, Ris? National Laboratory, the Danish Institute for
Food and Veterinary Research, the Danish Institute for Fisheries
Research, the Danish National Space Center and the Danish Transport
Research Institute have been merged with the Technical University of
Denmark (DTU) with DTU as the continuing unit. 




From percival at physics.queensu.ca  Mon Feb 12 17:23:07 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Mon, 12 Feb 2007 11:23:07 -0500
Subject: [neutron-mc] Source_Adapt problem
Message-ID: <000001c74ec2$18b22ae0$4a1680a0$@queensu.ca>

Hello all,

 

I have a question concerning the Source_Adapt component and its energy (or
wavelength) distribution.  My simplified understanding is that this
component initially emits a uniform energy (or wavelength) distribution,
then throughout a simulation adapts this distribution to increase counts in
the adapt_check component.  However, I am having a problem getting an
initial uniform distribution in wavelength, in fact, it seems as though the
distribution models some sort of exponential decay between E0 -dE and E0 +
dE.  In order to test this, I did two things.  I first set up a simulation
where no adaptation in energy was required by creating a simple instrument
with a Source_Adapt, two slits, a collimator and an adapt_check, I've
attached a file showing the wavelength distribution from an L_monitor placed
right after the source (pos1_plot.jpg).   I then added a powder sample to
the previous instrument and moved the adapt_check behind it, but set beta =
1 on the source_adapt (to force the component to distribute all neutrons
uniformly among all bins, i.e. no adaptation).  The wavelength distribution
seen from the source is in pos2_plot.jpg and it the same as the previous.

 

Here is the component definition for the Souce_Adapt in my instrument:

 

COMPONENT source = Source_Adapt(

    N_E = 20, N_xpos = 20, N_xdiv = 20, xmin = -0.005,

    xmax = 0.005, ymin = -0.005, ymax = 0.005, dist = 6.134,

    xw = 0.01, yh = 0.01, E0 =180, dE = 172, alpha = 0.25,

    beta = 0.25)

  AT (0, 0, 0) RELATIVE origin

 

E0 and dE chosen to distribute wavelengths between approximately 0.5AA and
3.2AA.

 

What would like is for the Source_adapt to initially give a uniform
distribution, much like the Source_Simple component.   Is this possible with
my current parameters for E0 and dE?

 

Regards,

 

Aaron

 

*************************************************

Aaron M. Percival

M.Sc. Candidate

Dept. of Physics, Engineering Physics & Astronomy

Queen's University

Kingston, Ontario, Canada, K7L 3N6

Office: 613-533-6000 ext. 74789

Fax: 613-533-6463

*************************************************

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070212/6a965c7d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pos2_plot.jpg
Type: image/jpeg
Size: 29412 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070212/6a965c7d/attachment.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pos1_plot.jpg
Type: image/jpeg
Size: 29585 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070212/6a965c7d/attachment-0001.jpg>

From peter.willendrup at risoe.dk  Wed Feb 14 15:13:08 2007
From: peter.willendrup at risoe.dk (Willendrup, Peter Kj=?ISO-8859-1?B?5g==?=r)
Date: Wed, 14 Feb 2007 15:13:08 +0100
Subject: [neutron-mc] Source_Adapt problem
In-Reply-To: <000001c74ec2$18b22ae0$4a1680a0$@queensu.ca>
Message-ID: <C1F8D78B.EFC%PKWI@exchange.Risoe.dk>

Hello Aaron,


On 2/12/07 5:23 PM, "Aaron M. Percival" <percival at physics.queensu.ca> wrote:

> I have a question concerning the Source_Adapt component and its energy (or
> wavelength) distribution.  My simplified understanding is that this
> component initially emits a uniform energy (or wavelength) distribution,
> then throughout a simulation adapts this distribution to increase counts in
> the adapt_check component.  However, I am having a problem getting an
> initial uniform distribution in wavelength, in fact, it seems as though the
> distribution models some sort of exponential decay between E0 -dE and E0 +
> dE.  

The Source_adapt component allows you to specify a flat 'start out'
distribution in energy, defined by E0 and dE, it does not have a set of
similar handles to specify a flat wavelength distribution. As the
relationship between energy and the (de-Broglie) wavelength of the neutron
is not linear, you will never get a flat distribution in wavelength by
specifying one in energy.

> In order to test this, I did two things.  I first set up a simulation
> where no adaptation in energy was required by creating a simple instrument
> with a Source_Adapt, two slits, a collimator and an adapt_check, I've
> attached a file showing the wavelength distribution from an L_monitor placed
> right after the source (pos1_plot.jpg).   I then added a powder sample to
> the previous instrument and moved the adapt_check behind it, but set beta =
> 1 on the source_adapt (to force the component to distribute all neutrons
> uniformly among all bins, i.e. no adaptation).  The wavelength distribution
> seen from the source is in pos2_plot.jpg and it the same as the previous.

As I understand the written meaning of the 'beta' parameter, it defines the
fraction of neutrons to emit uniformly, hence _not_ to adapt using
Adapt_check. This in effect means that setting beta=1 is equivalent to not
having an Adapt_check component at all. I might be missing something here,
but it does occur to me that your to test cases are in fact the same
situation?

> What would like is for the Source_adapt to initially give a uniform
> distribution, much like the Source_Simple component.   Is this possible with
> my current parameters for E0 and dE?

I guess that the easiest way to achieve that would be to write your own
component inspired from the existing Source_adapt, sampling uniformly in
wavelength instead of energy.

To summarize the behaviour of Source_adapt etc, I have written a small
example instrument (test.instr) which I attach that might explain things
slightly better. It has two input parameters, Beta is a handle to control
the 'beta' value of the Source_adapt, and CHECK is used in a WHEN statement
to control if the Adapt_check component is active. There is a powder sample
in there as in your case, which should let neutron distributions 'converge'
to a single value of energy/wavelength.

I have also attached graphical output from 4 simulation runs with

A) Beta = 1, CHECK = 0 (No Adapt_check and only 'uniform neutrons')
B) Beta = 1, CHECK = 1 (Adapt_check and only 'uniform neutrons')
C) Beta = 0, CHECK = 0 (No Adapt_check and no 'uniform neutrons')
D) Beta = 0, CHECK = 1 (Adapt_check and no 'uniform neutrons')

The graphics from cases A->C are identical, only when using a low Beta (zero
is ofcourse quite extreme) the incoming energy and wavelength distributions
are affected.


Hope this helps,

Peter

-------------- next part --------------
A non-text attachment was scrubbed...
Name: beta_0_check_0.png
Type: application/octet-stream
Size: 14047 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070214/3ba5ae61/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: beta_0_check_1.png
Type: application/octet-stream
Size: 13858 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070214/3ba5ae61/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: beta_1_check_0.png
Type: application/octet-stream
Size: 14437 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070214/3ba5ae61/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.instr
Type: application/octet-stream
Size: 2968 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070214/3ba5ae61/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: beta_1_check_1.png
Type: application/octet-stream
Size: 14257 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070214/3ba5ae61/attachment-0004.obj>

From percival at physics.queensu.ca  Wed Feb 14 15:28:34 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Wed, 14 Feb 2007 09:28:34 -0500
Subject: [neutron-mc] Source_Adapt problem
In-Reply-To: <C1F8D78B.EFC%PKWI@exchange.Risoe.dk>
References: <000001c74ec2$18b22ae0$4a1680a0$@queensu.ca>
	<C1F8D78B.EFC%PKWI@exchange.Risoe.dk>
Message-ID: <001801c75044$6bb3c450$431b4cf0$@queensu.ca>

Thanks for your reply Peter.

Taking your suggestion, I will attempt to modify the Source_Apapt component
to emit a flat 'start-out' distribution in wavelength rather than in energy.

For my purposes, this is what I need.

Regards,

Aaron

-----Original Message-----
From: neutron-mc-bounces at risoe.dk [mailto:neutron-mc-bounces at risoe.dk] On
Behalf Of Willendrup, Peter Kj?r
Sent: Wednesday, February 14, 2007 9:13 AM
To: McStas users list at neutron.risoe.dk
Subject: Re: [neutron-mc] Source_Adapt problem

Hello Aaron,


On 2/12/07 5:23 PM, "Aaron M. Percival" <percival at physics.queensu.ca> wrote:

> I have a question concerning the Source_Adapt component and its energy 
> (or
> wavelength) distribution.  My simplified understanding is that this 
> component initially emits a uniform energy (or wavelength) 
> distribution, then throughout a simulation adapts this distribution to 
> increase counts in the adapt_check component.  However, I am having a 
> problem getting an initial uniform distribution in wavelength, in 
> fact, it seems as though the distribution models some sort of 
> exponential decay between E0 -dE and E0 + dE.

The Source_adapt component allows you to specify a flat 'start out'
distribution in energy, defined by E0 and dE, it does not have a set of
similar handles to specify a flat wavelength distribution. As the
relationship between energy and the (de-Broglie) wavelength of the neutron
is not linear, you will never get a flat distribution in wavelength by
specifying one in energy.

> In order to test this, I did two things.  I first set up a simulation 
> where no adaptation in energy was required by creating a simple 
> instrument with a Source_Adapt, two slits, a collimator and an 
> adapt_check, I've attached a file showing the wavelength distribution from
an L_monitor placed
> right after the source (pos1_plot.jpg).   I then added a powder sample to
> the previous instrument and moved the adapt_check behind it, but set 
> beta =
> 1 on the source_adapt (to force the component to distribute all 
> neutrons uniformly among all bins, i.e. no adaptation).  The 
> wavelength distribution seen from the source is in pos2_plot.jpg and it
the same as the previous.

As I understand the written meaning of the 'beta' parameter, it defines the
fraction of neutrons to emit uniformly, hence _not_ to adapt using
Adapt_check. This in effect means that setting beta=1 is equivalent to not
having an Adapt_check component at all. I might be missing something here,
but it does occur to me that your to test cases are in fact the same
situation?

> What would like is for the Source_adapt to initially give a uniform
> distribution, much like the Source_Simple component.   Is this possible
with
> my current parameters for E0 and dE?

I guess that the easiest way to achieve that would be to write your own
component inspired from the existing Source_adapt, sampling uniformly in
wavelength instead of energy.

To summarize the behaviour of Source_adapt etc, I have written a small
example instrument (test.instr) which I attach that might explain things
slightly better. It has two input parameters, Beta is a handle to control
the 'beta' value of the Source_adapt, and CHECK is used in a WHEN statement
to control if the Adapt_check component is active. There is a powder sample
in there as in your case, which should let neutron distributions 'converge'
to a single value of energy/wavelength.

I have also attached graphical output from 4 simulation runs with

A) Beta = 1, CHECK = 0 (No Adapt_check and only 'uniform neutrons')
B) Beta = 1, CHECK = 1 (Adapt_check and only 'uniform neutrons')
C) Beta = 0, CHECK = 0 (No Adapt_check and no 'uniform neutrons')
D) Beta = 0, CHECK = 1 (Adapt_check and no 'uniform neutrons')

The graphics from cases A->C are identical, only when using a low Beta (zero
is ofcourse quite extreme) the incoming energy and wavelength distributions
are affected.


Hope this helps,

Peter







From ajain at barc.gov.in  Sun Feb 18 11:34:02 2007
From: ajain at barc.gov.in (ajain at barc.gov.in)
Date: Sun, 18 Feb 2007 16:04:02 +0530 (IST)
Subject: [neutron-mc] Elastically bent perfect crystal as monochromator
Message-ID: <1171794842.45d82b9aebac1@bts.barc.ernet.in>

Hello all,
               I am using current version (McStas-1.10-i686-Win32.exe )of 
McStas for simulating an instrument for measuring residual stress .I wnat to 
know whether it is possible to use Elastically bent perfect crystal(silicon) 
as monochromator in the simulation and what is the process. 

Look forward to hear you.
With kind regards,
Anil Jain

-------------------------------------------------




From peter.willendrup at risoe.dk  Wed Feb 21 10:52:13 2007
From: peter.willendrup at risoe.dk (Willendrup, Peter Kj=?ISO-8859-1?B?5g==?=r)
Date: Wed, 21 Feb 2007 10:52:13 +0100
Subject: [neutron-mc] Elastically bent perfect crystal as monochromator
In-Reply-To: <1171794842.45d82b9aebac1@bts.barc.ernet.in>
Message-ID: <C201D4DD.1294%PKWI@exchange.Risoe.dk>

Hello Anil,


On 2/18/07 11:34 AM, "ajain at barc.gov.in" <ajain at barc.gov.in> wrote:

>                I am using current version (McStas-1.10-i686-Win32.exe )of
> McStas for simulating an instrument for measuring residual stress .I wnat to
> know whether it is possible to use Elastically bent perfect crystal(silicon)
> as monochromator in the simulation and what is the process.

What you are describing here would certainly be very useful for the
community. However:

* We have no model of elastically bent crystals what so ever.
* The single crystal component (which includes secondary extinction but does
_not_ handle primary extinction if you need that?) would in principle be a
starting point. But, as recently reported by Aaron Percival
( http://mailman.risoe.dk/pipermail/neutron-mc/2007q1/000689.html ), there
are currently problems when trying to use it as a monochromator. (We are
currently debugging to investigate further).

If you decide to try to write your own component for your project, we would
certainly welcome it as a contribution to the project.


Best regards,

Peter Willendrup


--  
-----------------------------------------------
Peter Willendrup, MSc
Materials Research Department
Building 227
Ris? National Laboratory
Technical University of Denmark
Frederiksborgvej 399
DK-4000 Roskilde
Denmark

Tel:  (+45) 4677 5862
Fax:  (+45) 4677 5758
Email: Peter.Willendrup at risoe.dk

-----------------------------------------------
>From 1 January 2007, Ris? National Laboratory, the Danish Institute for
Food and Veterinary Research, the Danish Institute for Fisheries
Research, the Danish National Space Center and the Danish Transport
Research Institute have been merged with the Technical University of
Denmark (DTU) with DTU as the continuing unit.





From percival at physics.queensu.ca  Mon Feb 26 16:58:21 2007
From: percival at physics.queensu.ca (Aaron M. Percival)
Date: Mon, 26 Feb 2007 10:58:21 -0500
Subject: [neutron-mc] Modified Source_Adapt
Message-ID: <000701c759be$f5d2f510$e178df30$@queensu.ca>

Hello all,

 

As a follow up to a previous posting, I have modified the Source_Adapt
component to allow for an initial distribution that is uniform in wavelength
to be used and have added it to the existing Source_Adapt component.  The
component now allows for either mean energy and energy spread or mean
wavelength and wavelength spread to be inputted but the user, then will give
a 'start out' neutron distribution that is either uniform in energy or in
wavelength (as was defined in the input).

 

I have attached the modified component, and revised the documentation to
reflect the changes.

 

Aaron

 

*************************************************

Aaron M. Percival

M.Sc. Candidate

Dept. of Physics, Engineering Physics & Astronomy

Queen's University

Kingston, Ontario, Canada, K7L 3N6

Office: 613-533-6000 ext. 74789

Fax: 613-533-6463

*************************************************

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070226/bf7eaddb/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Source_adapt.comp
Type: application/octet-stream
Size: 9172 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070226/bf7eaddb/attachment.obj>

From peter.willendrup at risoe.dk  Wed Feb 28 09:19:59 2007
From: peter.willendrup at risoe.dk (peter.willendrup at risoe.dk)
Date: Wed, 28 Feb 2007 09:19:59 +0100
Subject: [neutron-mc] Modified Source_Adapt
References: <000701c759be$f5d2f510$e178df30$@queensu.ca>
Message-ID: <6463F256A85DC14CAB07F087A247997228A422@EXCHG-VS1.risoe.dk>

Hi Aaron and neutron-mc list,


Thanks for the contribution. We will let it undergo the usual test procedure and
merge the edits into the official Source_adapt of the next release.


Thanks,

Peter

-----Original Message-----
From: neutron-mc-bounces at risoe.dk on behalf of Aaron M. Percival
Sent: Mon 26/02/2007 16:58
To: 'McStas users list at neutron.risoe.dk'
Subject: [neutron-mc] Modified Source_Adapt
 
Hello all,

 

As a follow up to a previous posting, I have modified the Source_Adapt
component to allow for an initial distribution that is uniform in wavelength
to be used and have added it to the existing Source_Adapt component.  The
component now allows for either mean energy and energy spread or mean
wavelength and wavelength spread to be inputted but the user, then will give
a 'start out' neutron distribution that is either uniform in energy or in
wavelength (as was defined in the input).

 

I have attached the modified component, and revised the documentation to
reflect the changes.

 

Aaron

 

*************************************************

Aaron M. Percival

M.Sc. Candidate

Dept. of Physics, Engineering Physics & Astronomy

Queen's University

Kingston, Ontario, Canada, K7L 3N6

Office: 613-533-6000 ext. 74789

Fax: 613-533-6463

*************************************************

 


-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 3167 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20070228/a304fa2e/attachment.bin>

From r.m.dalgliesh at rl.ac.uk  Thu Mar 15 15:36:13 2007
From: r.m.dalgliesh at rl.ac.uk (Robert Dalgliesh)
Date: Thu, 15 Mar 2007 14:36:13 +0000
Subject: [neutron-mc] scilab 4.1
Message-ID: <1173969373.24794.15.camel@ndv140.isis.cclrc.ac.uk>

Hi Peter, Emmanuel,
I've been trying to use scilab 4.1 with a trace using mcdisplay from
mcgui (mcstas-1.10) on a fedora core 6 system.
It seems to be broken as no output window appear at all and the
simulation simply completes as if you were not running a trace. 

I've not been able to trace the cause as no error messages appear.

having reread the instructions I've now installed scilab 4.0 and it is
perfectly happy as are Matlab and PGPLOT.

As I can't find any mention of a specific problem with scilab 4.1 in the
bugs or the installation instructions I thought I had better let you
know that there may be a problem.

All the best

Rob

---------------------------------------------
Dr. Robert Dalgliesh
ISIS
Rutherford Appleton Laboratory
Chilton
Didcot
Oxfordshire
OX11 0QX
Tel: +44(0)1235 445687
Fax: +44(0)1235 445720
e-mail: r.m.dalgliesh at rl.ac.uk 




From peter.willendrup at risoe.dk  Thu Mar 15 15:38:31 2007
From: peter.willendrup at risoe.dk (Willendrup, Peter Kj=?ISO-8859-1?B?5g==?=r)
Date: Thu, 15 Mar 2007 15:38:31 +0100
Subject: [neutron-mc] scilab 4.1
In-Reply-To: <1173969373.24794.15.camel@ndv140.isis.cclrc.ac.uk>
Message-ID: <C21F18F7.1548%PKWI@exchange.Risoe.dk>

Hi Rob and neutron-mc users,

On 3/15/07 3:36 PM, "Robert Dalgliesh" <r.m.dalgliesh at rl.ac.uk> wrote:

> I've been trying to use scilab 4.1 with a trace using mcdisplay from
> mcgui (mcstas-1.10) on a fedora core 6 system.
> It seems to be broken as no output window appear at all and the
> simulation simply completes as if you were not running a trace.
> 
> I've not been able to trace the cause as no error messages appear.
> 
> having reread the instructions I've now installed scilab 4.0 and it is
> perfectly happy as are Matlab and PGPLOT.
> 
> As I can't find any mention of a specific problem with scilab 4.1 in the
> bugs or the installation instructions I thought I had better let you
> know that there may be a problem.

Thanks for reporting! You are absolutely right that there is indeed a
problem with the new Scilab version in regard to McStas and also right in
the fact that this bug has not been documented anywhere - sorry about that.

For now, the solution will be to stick with Scilab 4.0 as you propose above.
A fix for McStas to Scilab 4.1 is _not_ likely to be in the next release as
other functioning plotters exist.

Note also that the next McStas release will bring the classical PGPLOT
McStas plotters to Windows!


Best regards,

Peter

--  
-----------------------------------------------
Peter Willendrup, MSc
Materials Research Department
Building 227
Ris? National Laboratory
Technical University of Denmark
Frederiksborgvej 399
DK-4000 Roskilde
Denmark

Tel:  (+45) 4677 5862
Fax:  (+45) 4677 5758
Email: Peter.Willendrup at risoe.dk

-----------------------------------------------
>From 1 January 2007, Ris? National Laboratory, the Danish Institute for
Food and Veterinary Research, the Danish Institute for Fisheries
Research, the Danish National Space Center and the Danish Transport
Research Institute have been merged with the Technical University of
Denmark (DTU) with DTU as the continuing unit.





From farhi at ill.fr  Thu Mar 15 15:39:51 2007
From: farhi at ill.fr (Emmanuel FARHI)
Date: Thu, 15 Mar 2007 15:39:51 +0100
Subject: [neutron-mc] scilab 4.1
In-Reply-To: <1173969373.24794.15.camel@ndv140.isis.cclrc.ac.uk>
References: <1173969373.24794.15.camel@ndv140.isis.cclrc.ac.uk>
Message-ID: <45F95AB7.1050707@ill.fr>

Thanks Robert,

It seems Scilab has decided to make totally obsolete the 'old graphics' 
mode that we have been using since Scilab 2.7. We'll have to work on 
this (and it's in our dev plan, but not for the next release).
Thanks for reporting.

Emmanuel.

Robert Dalgliesh wrote:
> Hi Peter, Emmanuel,
> I've been trying to use scilab 4.1 with a trace using mcdisplay from
> mcgui (mcstas-1.10) on a fedora core 6 system.
> It seems to be broken as no output window appear at all and the
> simulation simply completes as if you were not running a trace. 
>
> I've not been able to trace the cause as no error messages appear.
>
> having reread the instructions I've now installed scilab 4.0 and it is
> perfectly happy as are Matlab and PGPLOT.
>
> As I can't find any mention of a specific problem with scilab 4.1 in the
> bugs or the installation instructions I thought I had better let you
> know that there may be a problem.
>
> All the best
>
> Rob
>
> ---------------------------------------------
> Dr. Robert Dalgliesh
> ISIS
> Rutherford Appleton Laboratory
> Chilton
> Didcot
> Oxfordshire
> OX11 0QX
> Tel: +44(0)1235 445687
> Fax: +44(0)1235 445720
> e-mail: r.m.dalgliesh at rl.ac.uk 
>
>
> _______________________________________________
> neutron-mc mailing list
> neutron-mc at risoe.dk
> http://mailman.risoe.dk/mailman/listinfo/neutron-mc
>
>   

-- 
What's up Doc ?
--------------------------------------------
Emmanuel FARHI, http://www.ill.fr/Computing/people/Farhi  \|/ ____ \|/
CS-Group BP2/22,   Institut Laue-Langevin (ILL) Grenoble  ~@-/ oO \-@~
6 rue J. Horowitz, BP 156, 38042 Grenoble Cedex 9,France  /_( \__/ )_\
Work :Tel (33/0) 4 76 20 71 35. Fax (33/0) 4 76 48 39 06     \__U_/