Dear @riya,
Thanks for your questions.
-
No, You will obtain just Bq. And as you said, if you add the mass, you can obtain Bq/kg. Considering that you want the total activity, it will be given to you when you do not add any normalization. (ref. * page 26)
-
What you obtain from RESNUCLEi is the activity. If you are interested in the mass, you could obtain first the number of particles using the following formula:
A: Total activity
λ : Decay constant
N: Number of particles
and from there you an obtain the mass.
If this is not what you are looking for, please let me know what do you mean with the ‘mass’.
Look at this reference for more details: Radioactive decay - Wikipedia (rates seccion)
- Yes, it is possible. You have to use negative times. There you can obtain the information during those irradiation stages (ref * page 21)
4 & 5.
As you say, you could add a USRTRACK or USRBDX in the region of uranium. But, of course, your results will be ‘affected’ by the interaction of the gammas with the material.
Finally, regarding question 5 (and 4 as well), I’m still trying to find a better apporach, but I have two options for you, making use of userroutines:
Option 1:
For this case, you have to use the routines called USDRAW in mgdraw.f and USRRNC in usrrnc.f.
- Subroutine USRRNC is called every time a residual nucleus is produced.
- USDRAW is called after each particle interaction, but in particular, with the ICODE = 110, you can obtain just the products from the radiactive decay.
Then, you can dump the information of the residual nuclei when they are produced in USRRNC, and in USDRAW you can dump the products from the radiactive decay. You can even filter in USDRAW when there is a gamma generated using the variable JTRACK = 7.
Dumping this information, and a bit of processing of the data, you should be able to obtain the contribution of gammas from each radionuclides.
PROs: You are dumping all the necessary information to do the study
CONTRAs: The post-processing could be not so simple
Option 2:
Using the routine mentioned above (USRRNC) you can do the following. There is a parameter related to the weight of the particle ‘WEE’. What you can do is to modify the value setting WEE = 0 to all the isotopes except one (filtering by Z and A). In that case, your simulation will continue, but taking account only for the cases where WEE!=0. Then, all the gammas generated (if there is any) will be from that radionucleus in particular.
PROs: Easy to apply. Simple to do the scoring
CONTRAs: You will have to run a simulation for each case. This is particularly demanding considering that the material of your target is Uranium…
Regarding the thresholds, I would recommend you to use the option PROMPT in the EMFCUT card in order to just get rid of prompt ones, but not affecting the generated from the delay part.
If I find a better approach I let you know asap.
In case you need more clarification do not hesitate to ask.
Best,
André