In a heavy ion interaction simulation (with Li7 beam on Ta target), I have scored DOSE-EQ, DOSE and ENERGY using USRBIN and R-phi-z settings. A DCYSCORE with cooling time = 0 was added for all 3 of them. The dimension of this region is R = 5 cm and H = 10 cm and the medium inside = air.
The dose-eq output at a cooling time of 0 = 4.30E+03 pSv/s
The dose output at the same condition = 1.09E-02 GeV/g per s, which after multiplication by 1.6022E-07, leads to 1.74E+03 pGy/s.
The energy output at the same condition = 2.0860E-14 GeV/cc per s, which will lead to a very small value as compared to above two.
Can you please explain the reason of these discrepancies ? I expected the difference between DOSE-EQ and DOSE since the method of estimating these are different and also DOSE-EQ is by definition, equialent dose i.e. absorbed dose x radiation weighting factor and DOSE is by definition, absorbed dose. However, I expected that DOSE and ENERGY scoring will give similar result after proper conversion. It will be helpful if you please correct me.
From the .out file, Number of secondaries generated in inelastic interactions per beam particle: 49% is photon and 50% is neutron. Hence, can we infer that the dose observed at different cooling period is due to these particles? But the neutron life time is small. So at different cooling period, the dose is mainly because of the activation of air and hence the gamma emitted from the active air in the room where this R-phi-z geometry is located. In this case, I expected that DOSE-EQ and DOSE would be similar since radiation weighting factor for gamma = 1.