Problems modifying particle decay table

Dominik_Werthmueller · 2019-09-11T11:53:07.822Z

Hi

I used to use the following commands to modify the particle decay table to my needs:

/particle/select kaon0S
/particle/property/decay/select 0
/particle/property/decay/br 0
/particle/property/decay/select 1
/particle/property/decay/br 1
/particle/property/decay/dump

/particle/select pi0
/particle/property/decay/select 0
/particle/property/decay/br 1
/particle/property/decay/select 1
/particle/property/decay/br 0
/particle/property/decay/dump

/particle/select lambda
/particle/property/decay/select 0
/particle/property/decay/br 0
/particle/property/decay/select 1
/particle/property/decay/br 1
/particle/property/decay/dump

This was working with Geant 9.6 (if I remember correctly) but fails now on 10.5.

It basically boils down to the problem that after selecting one particle and dumping its decay table the same does not work anymore for a second particle:

Idle> /particle/select lambda
/particle/select lambda
Idle> /particle/property/decay/dump
/particle/property/decay/dump
G4DecayTable: lambda
0: BR: 0.639 [Phase Space] : proton pi-
1: BR: 0.358 [Phase Space] : neutron pi0

Idle> /particle/select pi0
/particle/select pi0
Idle> /particle/property/decay/dump
/particle/property/decay/dump
The particle has no decay table !! Command ignored.
Idle>

Is there anything that I’m missing here or is this a bug?

Many thanks,
Dominik

civanch · 2019-09-26T09:51:14.953Z

Hello,

Can you clarify is PhysicsList is instatiated already?

VI

Dominik_Werthmueller · 2019-09-26T13:24:02.056Z

Hi Vladimir

I think I’m in the idle state when I’m doing this, the physics list has already been loaded.!

I attached a comparison of G4DecayTableMessenger.cc (9.6 vs 10.5) where I think a bug has been introduced by removing the indicated curly bracket. The currentDecayTable pointer will not be set correctly when the currently selected particle changed.

Cheers,
Dominik

decay_table_bug|690x378