52Cr is a stable isotope of Cr (Z=24 and A=52). If neutrons undergo inelastic interaction with 52Cr, then it will be exited, and as a result, it will be deexcited to the ground state (Same situation for 12C).
Please look at the screenshot:
I know it is possible to simulate neutrons and a 52Cr target with a high energy resolution gamma detector to detect those gamma-rays, but as the probability of neutron inelastic interaction with 52Cr is not so high, therefore the efficiency of the simulation will be too low and needs a huge number of particles to achieve a suitable precision. I decided to simulate the decay scheme of 52Cr (in the screenshot) directly. I mean, I am interested in simulating those gamma-rays (with the exact cascades and energies) without using neutrons. In other words, those gamma-rays will be the primary particles. I tried to simulate this thing using the gps source definition like this:
/gps/ion 24 52 0 0
But I think it is not correct, as I could not see any gamma-rays in the defined radiation detector. I think the problem can be due to the fact that 52Cr is a stable isotope, and my gps source definition (above) is related to the natural spontaneous isotopes. Any help please?
If I am not mistaken, this lines implement the de-excitation of a 10MeV excited 52Cr nucleus per event
/gps/ion 24 52 0 10000
/gps/pos/radius 0.5 mm
/gps/pos/centre 0 0 0 mm
However, biasing of the physics may be the way to go,
Here, an example of macro for examples/extended/radioactivedecay/rdecay01
52Cr.mac.txt (355 Bytes)
52Cr.out.txt (8.5 KB)
Thank you for your reply,
Is it possible to use /rdecay01/fullChain true in gps source definition? I tried but I got this error message:
**** COMMAND NOT FOUND </rdecay01/fullChain true> ****
**** Batch is interrupted!! ****
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : UIMAN0123
issued by : G4UImanager::ApplyCommand
Command aborted (100)
Error code : 100
** This is just a warning message. **
-------- WWWW -------- G4Exception-END --------- WWWW -------
This command is specific to rdecay01 example.
You do not need it: full decay chain is the default in G4RadioactiveDecay. Make sur that this process is registered in your physics list.
I would recommend to run few events with verbosity :
make sur that you have in your main :
G4int precision = 4;