[mcstas-users] Ncount doubled? McStas 3.3

Tue Apr 18 13:19:08 CEST 2023

I can confirm this also holds for McXtrace, but only for the 'legacy' code.

McXtrace 3.x is OK.

Emmanuel.

Le 18/04/2023 à 13:07, Richard Wagner a écrit :
>
> Hi Peter,
>
> as an example I took the BTsimple.instr from the template section.
>
> The command to run McStas is: mcrun BTsimple.instr -d 
> BTsimple_20230418_115509 *-n 1000* lambda=10 dlambda=9.9 maxhd=3 
> maxvd=3 gw=0.1 gh=0.1 gL=50 gm=6 delta1=0 delta2=0
>
> and the reads (just as an example)
>
> *BT_in.dat *
>
> # Format: McCode with text headers # URL: http://www.mccode.org # 
> Creator: # Instrument: BTsimple.instr *# Ncount:**2000* # Trace: no # 
> Gravitation: no # Seed: 1681811710726221 # Directory: 
> /Users/wagnerrichard/BTsimple_20230418_115509 # Param: lambda=10 # 
> Param: dlambda=9.9 # Param: maxhd=3 # Param: maxvd=3 # Param: gw=0.1 # 
> Param: gh=0.1 # Param: gL=50 # Param: gm=6 # Param: delta1=0 # Param: 
> delta2=0 # Date: Tue Apr 18 11:55:10 2023 (1681811710) # type: 
> array_1d(101) # Source: BTsimple (BTsimple.instr) # component: BT_in # 
> position: 0 0 2 # title: Wavelength monitor # Ncount: 2000 # filename: 
> BT_in.dat # statistics: X0=2.95806; dX=1.79349; # signal: Min=0; 
> Max=1.42052e+11; Mean=7.26371e+09; # values: 7.33634e+11 2.44642e+11 
> 34 # xvar: L # yvar: (I,I_err) # xlabel: Wavelength [AA] # ylabel: 
> Intensity # xlimits: 0.1 19.9 # variables: L I I_err N 0.198019802 
> 0.06290653129 0.06290653129 1 0.3940594059 0 0 0 0.5900990099 0 0 0 
> 0.7861386139 0 0 0 0.9821782178 8.634588395e+10 8.634588395e+10 1 
> 1.178217822 1.129377287e+11 1.129377287e+11 1 1.374257426 0 0 0 
> 1.57029703 0 0 0 1.766336634 0 0 0 1.962376238 0 0 0 2.158415842 0 0 0 
> 2.354455446 0 0 0 2.55049505 1.045974734e+11 1.045974734e+11 1 
> 2.746534653 0 0 0 2.942574257 8.997598489e+10 8.997598489e+10 1 
> 3.138613861 7.835386438e+10 7.835386438e+10 1 3.334653465 
> 1.420521026e+11 1.004534555e+11 2 3.530693069 0 0 0 3.726732673 0 0 0 
> 3.922772277 4.930911027e+10 4.930911027e+10 1 4.118811881 0 0 0 
> 4.314851485 0 0 0 4.510891089 0 0 0 4.706930693 0 0 0 4.902970297 0 0 
> 0 5.099009901 0 0 0 5.295049505 0 0 0 5.491089109 0 0 0 5.687128713 0 
> 0 0 5.883168317 2.834814784e+10 2.005320461e+10 2 6.079207921 0 0 0 
> 6.275247525 0 0 0 6.471287129 0 0 0 6.667326733 1.878195926e+10 
> 1.328092168e+10 2 6.863366337 0 0 0 7.059405941 0 0 0 7.255445545 0 0 
> 0 7.451485149 0 0 0 7.647524752 5150986989 5150986989 1 7.843564356 0 
> 0 0 8.03960396 4013292994 4013292994 1 8.235643564 3595007505 
> 3595007505 1 8.431683168 0 0 0 8.627722772 0 0 0 8.823762376 0 0 0 
> 9.01980198 0 0 0 9.215841584 2189836065 2189836065 1 9.411881188 
> 1988533204 1988533204 1 9.607920792 0 0 0 9.803960396 0 0 0 10 
> 1492967095 1492967095 1 10.1960396 0 0 0 10.39207921 0 0 0 10.58811881 
> 0 0 0 10.78415842 0 0 0 10.98019802 0 0 0 11.17623762 0 0 0 
> 11.37227723 848524385.4 848524385.4 1 11.56831683 776540913.4 
> 776540913.4 1 11.76435644 704975245.4 704975245.4 1 11.96039604 0 0 0 
> 12.15643564 0 0 0 12.35247525 0 0 0 12.54851485 0 0 0 12.74455446 0 0 
> 0 12.94059406 0 0 0 13.13663366 0 0 0 13.33267327 0 0 0 13.52871287 0 
> 0 0 13.72475248 361109386.4 361109386.4 1 13.92079208 0 0 0 
> 14.11683168 0 0 0 14.31287129 0 0 0 14.50891089 272324491.7 
> 272324491.7 1 14.7049505 0 0 0 14.9009901 241800376.7 241800376.7 1 
> 15.0970297 224371882.4 224371882.4 1 15.29306931 203464827.5 
> 203464827.5 1 15.48910891 0 0 0 15.68514851 182210050 182210050 1 
> 15.88118812 173421371.5 173421371.5 1 16.07722772 0 0 0 16.27326733 0 
> 0 0 16.46930693 0 0 0 16.66534653 140822062.6 140822062.6 1 
> 16.86138614 0 0 0 17.05742574 124243190.4 124243190.4 1 17.25346535 0 
> 0 0 17.44950495 0 0 0 17.64554455 0 0 0 17.84158416 0 0 0 18.03762376 
> 91358644.93 91358644.93 1 18.23366337 90222429.45 90222429.45 1 
> 18.42970297 0 0 0 18.62574257 0 0 0 18.82178218 0 0 0 19.01782178 0 0 
> 0 19.21386139 0 0 0 19.40990099 65966180.97 65966180.97 1 19.60594059 
> 0 0 0 19.8019802 0 0 0
>
> Holds for the rest of the dat-files and for simulations with higer n 
> numbers
> I just took n = 1000 for demonstration purposes.
>
> Best
> Richard
>
> On 18/04/2023 12:07, Peter Kjær Willendrup wrote:
>> Hi Richard,
>>
>>
>> Hmm, no this sounds strange?
>>
>> Could you help by providing a little bit more context? I.e.
>>
>> * What are you running?
>> * How?
>> * Where do you read off the resulting Ncount?
>>
>> Best,
>> Peter
>>
>>> On 18 Apr 2023, at 11.58, Richard Wagner <wagnerrichard at ill.fr> wrote:
>>>
>>> Dear McStas Experts and Community
>>>
>>> Is there a reason that the in Mcstas 3.3 Ncount is twice as high 
>>> then the "Ray count" selected?
>>>
>>> eg.:/-n 1000/ -> /# Ncount: 2000 /
>>>
>>> Best
>>>
>>> Richard
>>>
>>> -- 
>>> *Richard Wagner*
>>> Postdoctoral Research Engineer
>>> Nuclear and Particle Physics Group
>>> Institut Laue-Langevin - ILL
>>> 71, avenue des Martyrs
>>> CS 20156
>>> 38042 Grenoble Cedex 9
>>> France
>>>
>>> www.ill.eu <x-msg://1/www.ill.eu>
>>>
>>> _______________________________________________
>>> mcstas-users mailing list
>>> mcstas-users at mcstas.org
>>> https://mailman2.mcstas.org/mailman/listinfo/mcstas-users
>>
> -- 
> *Richard Wagner*
> Postdoctoral Research Engineer
> Nuclear and Particle Physics Group
> Institut Laue-Langevin - ILL
> 71, avenue des Martyrs
> CS 20156
> 38042 Grenoble Cedex 9
> France
>
> www.ill.eu <www.ill.eu>
>
>
> _______________________________________________
> mcstas-users mailing list
> mcstas-users at mcstas.org
> https://mailman2.mcstas.org/mailman/listinfo/mcstas-users
-- 
signature
/ ___|__/\_| |   | ____|_ _| |      FARHI Emmanuel
\___ \\ | |   |  _|  | || |      Div Exp/Data Reduction and Analysis Team
  ___) /_  _| |___| |___ | || |___   Tel : +33 (1) 69 35 96 04 
<callto:+33%20%281%29%2069%2035%2096%2004>
|____/  \/ |_____|_____|___|_____|  Saint-Aubin BP 48 - 91192 GIF/YVETTE 
CEDEX
SYNCHROTRON http://www.synchrotron-soleil.fr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20230418/7c862dd7/attachment-0001.html>