Problem with restManager REST_ViewEvents

RestG4
1

System : Linux
REST version : v2.2.10
Hi:
My simulation is about the neutrinoless double-beta decay of Xe136 to excited
states of Ba136. So each event will emit two electrons and two gamma. The gamma
's energy are 760keV and 819keV and the total kinetic energy of two electrons are 879keV.
image
Here I used REST_ViewEvents to check the track information of my simulation.

image
image860×22 19.7 KB

But in the red frame l marked, l am confused with the e-(563.30keV)'s compositon.(eventID=1048)
image
image842×524 117 KB

In other words,why there are two gamma whose energy are 29.76keV and 4.13keV.
Does the energy conservation work here(563.3≠29.76+4.13)?
Or there is some bug with the code?

2

Hi, I believe is more a problem with the representation on the REST_EventViewer.

I believe, the 563.30keV given in the viewer is the initial kinetic energy of the electron, not deposited, at position (-220,302.4,-122.0). What it is confusing here is that the pink it is just a marker of the position and initial kinetic energy. Try to remove the pink lines mentally xD.

So, the 2 gammas are secondaries of the electron track, but the electron did NOT deposit 563.30keV. You must check the total energy deposited by the event, or that particular track. Again, what is given in the viewer is the initial kinetic energy of the particle. The electron then experienced an interaction where it loosed energy and gave it to the gamma.

Could you print out and attach the output of the PrintEvent method? This may help to understand the process.

You can achieve that, for example, by using this in a restRoot shell

TRestRun *run = new TRestRun();
run->OpenInputFile("G4file.root");
TRestG4Event *ev = new TRestG4Event();
run->SetInputEvent(ev);
run->GetEntry(0); // or any entry you want to retrieve
ev->PrintEvent();
3

G4event1048.txt (58.2 KB)
Yeah,l have used REST_PrintEvents.C to output the process

4

I cannot find REST_PrintEvents.C in the repository.

TRestRun *run = new TRestRun();
run->OpenInputFile("G4file.root");
TRestG4Event *ev = new TRestG4Event();
run->SetInputEvent(ev);
run->GetEntry(0); // I assume is your first entry
ev->GetTrackByID(13)->GetEnergy()

The energy you get should not be 563.3, but lower than 563.3 - 29.76 - 4.13, so that energy conservation is not violated.

5

REST_PrintEvents.txt (2.6 KB)
That’s what l used. l just replace the suffix .C with .txt to upload this one.

6

So, that’s your own script, or it is that in the repository?

7

I have updated the settings of the site so that it is allowed to upload macros with .C extension.

8

Yes,l added it by myself.

9

May be it would be good if you or @nkx upload it to the repository?