[mcstas-users] Ncount doubled? McStas 3.3

Tue Apr 18 15:13:35 CEST 2023

Hi both,

OK, thanks for testing.

Looking on my machine here with McStas 3.3, it seems the the first entry in mccode.sim (i.e. the “simulation level”) is OK, whereas the entries for any monitor after has been doubled:

McStas simulation description file for BNL_H8.

Date:    Tue Apr 18 12:02:56 2023

Program:

begin instrument: BNL_H8

  File: BNL_H8_20230418_120254/mccode

  Source: BNL_H8.instr

  Parameters:  lambda(double)

  Trace_enabled: yes

  Default_main: yes

  Embedded_runtime: yes

end instrument

begin simulation: BNL_H8_20230418_120254

  Format: McCode with text headers

  URL: http://www.mccode.org

  Creator:

  Instrument: BNL_H8.instr

  Ncount: 1000000

  Trace: no

  Gravitation: no

  Seed: 1681812175442876

  Directory: BNL_H8_20230418_120254

  Param: lambda=2.36

end simulation

begin data

  Date: Tue Apr 18 12:02:57 2023 (1681812177)

  type: array_2d(20, 20)

  Source: BNL_H8 (BNL_H8.instr)

  component: D0_Source

  position: 0 0 0.0001

  title: PSD monitor

  Ncount: 2000000

  filename: D0_Source.psd

  statistics: X0=-0.00187139; dX=0.864864; Y0=0.0044862; dY=1.56205;

  signal: Min=6.0781e-06; Max=1.04372e-05; Mean=7.97501e-06;

  values: 0.00319 1.39945e-05 51960

  xvar: X

  yvar: Y

  xlabel: X position [cm]

  ylabel: Y position [cm]

  zvar: I

  zlabel: Signal per bin

  xylimits: -1.5 1.5 -2.7 2.7

  variables: I I_err N

end data

...

I would think this “only" a header / summary data - issue, but of course annoying. I don’t believe the actual number of simulated events is actually doubled.

Will investigate further during the next days and report back.

Best,
Peter

On 18 Apr 2023, at 13.19, Emmanuel FARHI <emmanuel.farhi at synchrotron-soleil.fr> wrote:

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<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><mailto: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<mailto: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<x-msg://3/www.ill.eu>

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

--
/ ___|__/\_| |   | ____|_ _| |      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<http://www.synchrotron-soleil.fr/>
_______________________________________________
mcstas-users mailing list
mcstas-users at mcstas.org<mailto:mcstas-users at mcstas.org>
https://mailman2.mcstas.org/mailman/listinfo/mcstas-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20230418/3d097f29/attachment-0001.html>