Problem 1862 - Absurdly large doses generated in example B3a
Summary: Absurdly large doses generated in example B3a
Status: RESOLVED FIXED
Alias: None
Product: Examples/Basic and Novice
Classification: Unclassified
Component: B3 (show other problems)
Version: 10.2
Hardware: PC Linux
: P5 normal
Assignee: Ivana Hrivnacova
URL:
Depends on:
Blocks:
 
Reported: 2016-05-18 23:30 CEST by James Matta
Modified: 2016-05-20 12:40 CEST (History)
2 users (show)

See Also:


Attachments
This is the output from running run1.mac, run2.mac and the manual stuff described above, I separate them with "OUTPUT FROM ..." where ... is run1.mac, run2.mac, or MANUAL (59.47 KB, text/plain)
2016-05-18 23:30 CEST, James Matta
Details

Note You need to log in before you can comment on or make changes to this problem.
Description James Matta 2016-05-18 23:30:09 CEST
Created attachment 400 [details]
This is the output from running run1.mac, run2.mac and the manual stuff described above, I separate them with "OUTPUT FROM ..." where ... is run1.mac, run2.mac, or MANUAL

I started example B3a in interactive mode. Then I ran the following commands:
  /vis/disable
  /run/initialize
  /run/beamOn 1000000

The output dose was insane:

G4WT1 > --------------------End of Local Run------------------------
G4WT1 >   The run was 499444 F18; Nb of 'good' e+ annihilations: 64456
G4WT1 >  Total dose in patient : 14.5659 nanoGy
G4WT1 > ------------------------------------------------------------
G4WT1 >
G4WT0 > Thread-local run terminated.
G4WT0 > Run Summary
G4WT0 >   Number of events processed : 500556
G4WT0 >   User=104.78s Real=52.45s Sys=0.01s
G4WT0 > Go to merge parameters
G4WT0 >
G4WT0 > --------------------End of Local Run------------------------
G4WT0 >   The run was 500556 F18; Nb of 'good' e+ annihilations: 64614
G4WT0 >  Total dose in patient : 1.45725e+17 Gy
G4WT0 > ------------------------------------------------------------
G4WT0 >
Run terminated.
Run Summary
Number of events processed : 1000000
User=104.8s Real=52.47s Sys=0.03s
--------------------End of Global Run-----------------------
The run was 1000000 events ; Nb of 'good' e+ annihilations: 129070
Total dose in patient : 1.45725e+17 Gy
------------------------------------------------------------

As you can see, one worker thread has a sane does, the other is insane. 145PetaGreys would vaporize (with prejudice) any material that received such a dose. There is no way that a mere 1,000,000 positrons and 2,000,000 511 keV photons could give that dose.

This effect scales, running "/run/beamOn 2000000" yields:

G4WT0 > --------------------End of Local Run------------------------
G4WT0 >   The run was 1001000 F18; Nb of 'good' e+ annihilations: 129484
G4WT0 >  Total dose in patient : 2.91575e+17 Gy
G4WT0 > ------------------------------------------------------------
G4WT0 >
G4WT1 > Thread-local run terminated.
G4WT1 > Run Summary
G4WT1 >   Number of events processed : 999000
G4WT1 >   User=215.26s Real=107.73s Sys=0.05s
G4WT1 > Go to merge parameters
G4WT1 >
G4WT1 > --------------------End of Local Run------------------------
G4WT1 >   The run was 999000 F18; Nb of 'good' e+ annihilations: 128768
G4WT1 >  Total dose in patient : 29.1472 nanoGy
G4WT1 > ------------------------------------------------------------
G4WT1 >
Run terminated.
Run Summary
Number of events processed : 2000000
User=215.27s Real=107.73s Sys=0.05s
--------------------End of Global Run-----------------------
The run was 2000000 events ; Nb of 'good' e+ annihilations: 258252
Total dose in patient : 2.91575e+17 Gy
------------------------------------------------------------

When run in batch mode on run1.mac and run2.mac, it yields good results for run1.mac but does this insane behavior in run2.mac
Comment 1 Andrea Dotti 2016-05-19 00:39:00 CEST
Hi Ivana,
it looks like the problem is the new G4Parameter stuff. Thread 0 seems to have a correct initialization value, but not thread 1...
Andrea
P.S. The reports come from our tutorial
Comment 2 Ivana Hrivnacova 2016-05-20 12:40:19 CEST
The problem was in the initialization of the hits collections IDs in B3aEventAction.
It is now fixed in the development version and the fix will be included in the next patch to the 10.2 release.