In this Run there was generated 10 cycles with 10 primary photons in each. Bismuth target was irradiated. I want to calculate photoyield of isotopes. Here are obtained results from _sum.lis file:
A_min: 1 - A_max: 223
A: 209 5.5761835E-06 +/- 3.2775269E+01 %
A: 208 1.5990599E-04 +/- 3.7025284E+01 %
A: 207 1.8750001E-05 +/- 5.0917508E+01 %
A: 206 9.9875002E-05 +/- 2.6683871E+01 %
A: 1 6.2500003E-06 +/- 9.9000000E+01 %
Tell please is this result mean that there are induced for instance 9.9875002E-05 +/- 2.6683871E+01% isotopes with A=206 per one impinging photon in one cubic centimetre of bismuth target? Here is .inp file: 20200206_isotope_production_gold(1).inp (1.9 KB)
Correct (since you provided the target volume in the RESNUCLEi card). Note that in the produced output you can also find single Z,A isotopes (and not just the total over all given A isobars). Moreover, as you can see from the large statistical errors, 100 primaries is a by far too low statistics.
Ok, in one of Runs there was generated 10 cycles with 10E+07 primary photons in each. Here are results:
A_min: 1 - A_max: 223
A: 210 1.1649408E-08 +/- 4.3943631E+00 %
A: 209 1.9306620E-06 +/- 3.7562329E-01 %
A: 208 1.0068768E-05 +/- 4.8764563E-01 %
A: 207 3.8707355E-05 +/- 3.2845620E-01 %
A: 206 2.8991645E-04 +/- 9.3203536E-02 %
A: 205 2.9599643E-08 +/- 9.6229443E+00 %
A: 204 5.7292608E-08 +/- 7.4525954E+00 %
A: 203 6.6835541E-08 +/- 5.2215063E+00 %
A: 4 1.1617040E-07 +/- 5.7975322E+00 %
A: 3 6.3661991E-10 +/- 9.9000000E+01 %
A: 2 1.3591517E-07 +/- 6.3137558E+00 %
A: 1 7.6199301E-06 +/- 5.3171963E-01 %
I want obtain the total number of produced isotopes with for instance A=205. For this purpose should I multiply the volume of bismuth target (0,314159265 cubic centimetres) by total number of photons (10E+08) and by 2.9599643E-08 ? If yes it’s obtained that total number of isotopes A=205 produced in total volume equals to 0,929900209. How could this number be not integer ?
Here is .inp file:
20200206_isotope_production_gold(1).inp (2.1 KB)
The number of primary photons you have run has absolutely nothing to do with the proper normalization of the result. You have to multiply by the number of photons you have in reality, not by the number of photons you simulated (that is usually much lower). The latter simply affects the quality of the simulation result, via its statistical error.
Then, what’s the point of dividing by the target volume (putting it in the RESNUCLEi card) and later multiplying back? Just leave empty the respective RESNUCLEi parameter (implying a default value of 1 cm^3, for normalization purposes only) and be aware that the result is then not per cm^3 but referring to the whole target volume, as you aim to.
Last, the final number is not an integer because it’s the average of a statistical variable. Of course, for each primary photon one gets an integer number (often zero, sometimes one, hardly more).
OK, thank you, I’ll leave volume field empty to obtain result in whole volume. Also can you tell please what does mean numbers in WHAT(2) of RESNCULEi card ? For example what’s the difference between -21 and -54 ?
That’s just the unit number Fortran uses to write the file for RESNUCLEi output. The only difference you’ll see is that the corresponding output file name ends with 21 or 54.