The user manual on RAD-BIOL is incomplete. I have gone through different threads in the forum ( even I have also asked different queries earlier in the forum). From the threads, I got some idea on the format of .dat files need to be provided. Still I have two queries for which I did not get any satisfactory solutions:
What will be the unit of LET if the alpha and beta is provided as a function of LET?
What will be the unit of Kinetic energy per nucleon if these are provided as a function of KE/u.
In the .dat file, should I mention the unit of LET or KE/u somewhere?
Kindly reply me. I am waiting for the reply. Thanking you in advance.
Input files
Please upload all relevant files. (FLUKA input file, Flair project file, user routines, and data files)
indeed the manual needs some improvement on the explanation of the RAD-BIOL card.
Using LET or energy doesn’t change the units of the alpha (Gy^{-1}) & beta (Gy^{-2}) coefficients. Internally the code will convert the kinetic energy per nucleon to LET to interpolate the coefficients.
MeV/n
bonus answer, since you mentioned you have only two queries The default is ENERGY for the tabulations (MeV/n). However if you want to specify it in LET (keV/\mu m) you need to add as last line literally the string LET
else it will consider that the last line will contain (if it is not present) ENERGY
It follows an example of a file with random/nonsense values (which should not be used). The comments as well, should not be there
0.1 0.05 α(Gy^-1) β(Gy^-2) for X-rays
2 number of ion species Nion=2 (Z=1 and 2)
5 number of energy Ne=5 (or let) intervals for Z=1
0.1 1.5 0 energy(MeV/n) α1(Gy^-1) β1(Gy^-2)
0.5 1.0 0
1.0 0.7 0
2.0 0.5 0.004
10 0.2 0.04 energy(MeV/n) α5(Gy^-1) β5(Gy^-2)
5 number of energy Ne=5 (or let) intervals for Z=2
0.1 1.9 0
0.5 1.1 0
1.0 1.0 0
2.0 0.7 0.0004
10 0.5 0.004
ENERGY to be replaced with LET if the first column is given in LET (keV/um)
The default is ENERGY for the tabulations (