Problem 1836 - G4RadioactivityTable not filled in MT mode
Summary: G4RadioactivityTable not filled in MT mode
Status: RESOLVED WONTFIX
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models (show other problems)
Version: 10.2
Hardware: All All
: P5 trivial
Assignee: dennis.herbert.wright
URL:
Depends on:
Blocks:
 
Reported: 2016-02-29 11:33 CET by michel.maire
Modified: 2016-03-18 23:18 CET (History)
0 users

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description michel.maire 2016-02-29 11:33:11 CET

    
Comment 1 michel.maire 2016-02-29 11:37:56 CET
It seems that G4RadioactivityTable is not filled in MT mode (it is filled in sequential mode)
The problem can be seen in example rdecay02, Run::WriteActivity(), with proton-beam.mac
Comment 2 dennis.herbert.wright 2016-03-18 23:18:19 CET
Th immediate problem is the line

  if (isMaster) fRun->EndOfRun();

in RunAction::EndOfRunAction(const G4Run*)

The master RDM process carries no data, while each of the worker RDM threads do.
So no data gets written to the tables in MT mode.  A temporary fix would be to remove the test for master and just do 

fRun->EndOfRun();

Then the tables will be filled, but there will be a table for each thread at the end of the run so that some post-processing of the output will need to be done to collect everything.   

The proper fix, however, is deeper and requires a re-design of the RDM biased 
mode.  Basically radioactive decay collects run data in a table in the process, which in general, processes should not do.

So I will propose this redesign of the radioactive decay biased mode as a project and close this report as "WONTFIX"