Problem 2032

Summary:	Biasing mode in RadioactiveDecay
Product:	Geant4	Reporter:	michel.maire
Component:	processes/hadronic/models	Assignee:	dennis.herbert.wright
Status:	RESOLVED FIXED
Severity:	minor
Priority:	P4
Version:	10.4
Hardware:	All
OS:	All

Description michel.maire 2018-02-05 18:58:29 CET

Following macro is the decay Ca47 --> Sc47 --> Ti47
In biasing mode, Sc47 are never created and only ~30% of Ti47 are present.
What happens of the remaining 70 %

 /control/verbose 2
 /run/verbose 1
 #
 /gun/particle ion
 /gun/ion 20 47
 #
 ###/grdm/decayBiasProfile Ca47.time.data.txt
 ###/grdm/BRbias 0
 #
 /analysis/setFileName Ca47
 /analysis/h1/set 8  100  0. 30. d	#time of life
 #
 /run/printProgress 10000
 /run/beamOn 100000

Here the file Ca47.time.data.txt :
 0.	  1
 864000.  0

Comment 1 dennis.herbert.wright 2018-03-13 01:22:01 CET

The missing Ti47 nuclei were due to a bug in G4RadioactiveDecay::AddDeexcitationSpectrumForBiasMode: stable nuclei at the end of decay chains were not being added as secondaries.  This is now fixed in tag radioactive_decay-V10-04-01.

The missing Sc nuclei is not a bug.  Biased RDM decays all chains down to the stable nuclei and collects all the secondaries at the end of the tracking step.   Since the ground state of Sc47 is not stable, it does not appear.