Radioactive decay Am-241 and Pb-210

Physics Processes, Models and Cross Sections
3

Hello @svanlang , yes I use a macro to define the radionuclides:

/gps/particle ion
/gps/ion 82 210 0 0
/gps/energy 0. eV

I found this old topic Incorrect Am241 decay spectrum and I read @mkelsey reply in wich he suggested to limit the radionuclide by the /grdm/nucleusLimits command then I added in my code

for Pb210
/grdm/nucleusLimits 210 210 82 82

and for Am241
/grdm/nucleusLimits 241 241 95 95

and now I get

for Pb210

Pb210
Pb2101276×996 6.06 KB

and Am241

Am241
Am2411276×996 7.12 KB

I’m almost sure about Pb210 because on Table of Isotopes decay data I read only a peak at 45.5keV…
Regarding the Am241 I’ve some doubts because on Table of Isotopes decay data it gives so many peaks. Nevertheless, the Igamma% is at energy 59.6keV and I see a peak for this energy value.

I also tried to follow the radioactivedecay example, as suggested by @maire , then:

  1. I wrote trackingaction.cc , trackingaction.hh, trackingmessenger.cc, trackingactionmessenger.cc

TrackingAction.cc (3.9 KB)
TrackingMessenger.cc (2.9 KB)

  1. I added the code in actionvisualization
  TrackingAction* trackingAction = new TrackingAction(eventAction);
  SetUserAction(trackingAction);

B1ActionInitialization.cc (3.1 KB)

  1. I added in the macro

/rdecay01/fullChain false

Using this method, I don’t get errors during compiling and method, but I don’t get the deposited energy.

Thank you

1 Like