I have been working to simulate a scintillation detector with a heavy-ion beam and am struggling to show conservation of energy. In the attached flair file I am scoring the energy deposited in the scintillation volume and the OPTIPHOT’s that leave the scintillation volume. Combined, these 2 sources seem to account for maybe 1/4 the total energy that the beam should be depositing.
I am unsure if the error is in the quantities I am scoring or my interpretation of the results and hoping to get some assistance.
My understanding is the energy deposited + (the photon flux * the energy of the photons) should account for the bulk of the beams energy.
Energy deposited is determined from summing across a 1-D plot of USRBIN.
Photon flux is from USRBDX’s _sum.lis file’s (Part/pr) figure.
Actually, if this is done properly, one gets from your EneDep USRBIN scoring on unit 50 an energy deposition value that is almost 99% of the beam energy, in agreement with the energy conservation table that you can find at the the end of the standard output file.
In fact, the 1D projection plot gives you the average energy density [GeV/cm^3] profile, which - after being summed as you said - has to be multiplied by the transverse scoring area (34 cm x 34 cm) - over which it has been averaged - and the longitudinal bin size (0.07 cm), so getting to the exhaustive percentage I mentioned above.
The optical photon energy is already accounted in the deposited energy amount.
It shall be noted that energy is perfectly conserved, and the missing 1% is taken by particles reaching the blackhole material (typically neutrons) or is spent in nuclear reactions transforming kinetic energy into mass.