Problem 1825 - Particle definition of isotopes of high Z
Summary: Particle definition of isotopes of high Z
Status: RESOLVED INVALID
Alias: None
Product: Geant4
Classification: Unclassified
Component: particles (show other problems)
Version: 10.2
Hardware: All All
: P5 normal
Assignee: kurasige
URL:
Depends on:
Blocks:
 
Reported: 2016-02-08 13:54 CET by Gloria Luzon
Modified: 2016-02-15 09:57 CET (History)
1 user (show)

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
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.