Problem 1825

Summary:	Particle definition of isotopes of high Z
Product:	Geant4	Reporter:	Gloria Luzon <luzon>
Component:	particles	Assignee:	kurasige
Status:	RESOLVED INVALID
Severity:	normal	CC:	kurasige
Priority:	P5
Version:	10.2
Hardware:	All
OS:	All

Description Gloria Luzon 2016-02-08 13:54:08 CET

Using G4RadioactiveDecayPhysics in the PhysicsList, we tried to launch some isotopes and study their decays.

In case we try any isotope with Z>94 (i.e, above Pu) the code gives us a "G4Exception" message issued by G4ParticleGun::SetParticleDefinition() since null pointer was given. No problem with any other isotope we have tried (as Rn220, Th232, Co60 ..), it works fair.

We have checked the data base we are using: RadioactiveDecay4.2 and those isotopes are included.

Comment 1 kurasige 2016-02-14 06:15:48 CET

Hi,

 Thank you for problem report.
 
 I checked with examples/extended/hadronic/Hadr06 (which uses G4RadioactiveDecayPhysics), but I can not reproduce the problem.
 Cm242 is created successfully with commands below 
 and decays into Pu238 + alpha. (see below)
 
 Could you explain the detail of your application and environment?

Thanks and Cheers,
 

---- commands ---------------
/run/initialize
/gun/particle ion
/gun/ion 96 242 4
/tracking/verbose 1
/run/beamOn 1

--------------------
*********************************************************************************************************
* G4Track Information:   Particle = Cm242,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm      0 fm     14 MeV     0 meV     0 fm      0 fm   Water_ts   initStep
    1  -3.57 um  -2.18 um  -3.99 um      0 meV    14 MeV  5.78 um   5.78 um   Water_ts      ionIoni
    2  -3.57 um  -2.18 um  -3.99 um      0 meV     0 meV     0 fm   5.78 um   Water_ts   RadioactiveDecay

Comment 2 Gloria Luzon 2016-02-15 09:51:00 CET

Thans,

it was just our problem. We had put a limit in Z.

Sorry for the inconveniences

Comment 3 Gloria Luzon 2016-02-15 09:52:16 CET

(In reply to comment #1)
> Hi,
> 
>  Thank you for problem report.
> 
>  I checked with examples/extended/hadronic/Hadr06 (which uses
> G4RadioactiveDecayPhysics), but I can not reproduce the problem.
>  Cm242 is created successfully with commands below 
>  and decays into Pu238 + alpha. (see below)
> 
>  Could you explain the detail of your application and environment?
> 
> Thanks and Cheers,
> 
> 
> ---- commands ---------------
> /run/initialize
> /gun/particle ion
> /gun/ion 96 242 4
> /tracking/verbose 1
> /run/beamOn 1
> 
> --------------------
> *********************************************************************************************************
> * G4Track Information:   Particle = Cm242,   Track ID = 1,   Parent ID = 0
> *********************************************************************************************************
> 
> Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng   
> Volume     Process
>     0      0 fm      0 fm      0 fm     14 MeV     0 meV     0 fm      0 fm  
> Water_ts   initStep
>     1  -3.57 um  -2.18 um  -3.99 um      0 meV    14 MeV  5.78 um   5.78 um  
> Water_ts      ionIoni
>     2  -3.57 um  -2.18 um  -3.99 um      0 meV     0 meV     0 fm   5.78 um  
> Water_ts   RadioactiveDecay

Solved, just our problem.
Thanks.