Problem 1988 - G4AtomicTransitionManager::Initialise() not called by RadioactiveDecay
Summary: G4AtomicTransitionManager::Initialise() not called by RadioactiveDecay
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/radioactive_decay (show other problems)
Version: 10.3
Hardware: PC Linux
: P4 major
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2017-06-28 22:54 CEST by Pedro.Arce
Modified: 2017-12-18 22:23 CET (History)
2 users (show)

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description Pedro.Arce 2017-06-28 22:54:40 CEST
I got a crash in geant4.10.03.p01 using QGSP_BIC_AllHP at

   issued by : G4AtomicTransitionManager::Shell()
No de-excitation for Z= 3  shellIndex= 0. AtomicShell not found - check if data are uploaded

I found out that it was due that G4AtomicTransitionManager::Initialise() is not called as it used to be by

G4RadioactiveDecay::BuildPhysicsTable(const G4ParticleDefinition&)

because this code is commented out:

    /*
      p = new G4UAtomicDeexcitation();
      p->SetFluo(true);
      p->SetAuger(true);
      p->InitialiseAtomicDeexcitation();
      theManager->SetAtomDeexcitation(p);
      */

and then there is no shellTable when a G4ECDecay happens.

Pedro
Comment 1 Vladimir.Ivantchenko 2017-07-01 16:59:43 CEST
Hello Pedro,

thank you for the bug report. Please, note, that G4RadioactiveDecay itself must not change EM physics parameters, it is why these code is commented out. EM parameters should be changed inside physics list classes. 

I will try reproduce the problem.

Laurent, you may switch bug report to me.

Vladimir
Comment 2 Vladimir.Ivantchenko 2017-12-12 15:28:08 CET
Hello Pedro,

In Geant4 10.3 and 10.4 the Physics List QGSP_BIC_AllHP does not have radioactive decay enabled. Independently on your application, we checked, that G4RadioactiveDecayPhysics may be enabled on top of any Physics List and atomic de-excitation parameters are defined. Can you, please, confirm 

   a) what physics list was used in reality?
   b) is there a problem in 10.4?

Vledimir
Comment 3 Vladimir.Ivantchenko 2017-12-18 22:23:01 CET
I do not think that the problem exist in 10.4. Let us close this bug report.

Vladimir