Discrepancy between totalenergy and detector energy deposition

[Graeme Scott]

I've been trying to use the detector output file to calculate the total neutron energy deposited in that detector, and I've attached a fairly minimal working example here.

My problem is that when I sum the energy deposited from the detector output file (5.4E3 MeV), I get quite a discrepancy with what the total energy file tells me is deposited in that volume (1.307e+04 MeV).

I wonder if I'm missing something obvious that would explain this discrepancy?

Thanks very much,

Graeme

[Graeme Scott]

So this appears to be a bug to me, and I'm providing the details below in case it might help the developers track down what's going on.

I seem to be getting the behaviour that the total energy deposited as reported by the detector ntuple is dependent on whether the 'totalenergy' command is present or not. 

I found this to be the case for running on both Ubuntu and Windows via the GUI. There also appears to be a weak dependence on the physics list that I use.

In the attached example with:

physics QGSP_BERT_HP disable=Decay

totalenergy OFF: Sum(energies in detector ntuple) = 1.438e+06 MeV

totalenergy ON:   Sum(energies in detector ntuple) = 5.754e+05 MeV, while

                               totalenergy file reports 1.438e+06 MeV deposited in detector

physics QGSP_BIC_AllHP disable=Decay

totalenergy OFF: Sum(energies in detector ntuple) = 1.438e+06 MeV

totalenergy ON:   Sum(energies in detector ntuple) = 7.447e+05 MeV, while

                               totalenergy file reports 1.438e+06 MeV deposited in detector

I also tested removing the command disable=Decay, which had minimal effect.

If there are any other details that might be helpful please just let me know,

All the best,

Graeme

[Lucien Cremaldi]

Maybe energy is escaping the detector volume?  Add a surrounding volume to check.

 Lucien

On Thu, 11 Jan 2024 07:17:27 -0800 (PST), Graeme Scott wrote

> --
> You received this message because you are subscribed to the Google Groups "G4Beamline" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to g4beamline+...@muonsinc.com.
> To view this discussion on the web visit https://groups.google.com/a/muonsinc.com/d/msgid/g4beamline/551f8886-e9ff-4909-a69b-cc24b1455ba6n%40muonsinc.com.


--
Open WebMail Project (http://openwebmail.org)
Debian Project (http://www.debian.org)