[Gate-users] Splitting ROOT files according to time sampling times

[Ida Häggström]

Dear Gaters,

I've started looking at running dynamic phantom/sources, i.e. assigning a
TAC to each voxel. I haven't really figured out how to handle the ROOT
output though. Here I'm working with TACs with 401 posts, one per second
from 0 to 400 sec.

I have a time acquisition of  the form:

/gate/application/setTimeSlice     1.  s
/gate/application/setTimeStart     0.  s
/gate/application/setTimeStop      400. s
/gate/application/startDAQ

and for the source:

/gate/source/voxelSource/interfileReader/SetTimeSampling 1. s

and finally for the phantom:

/gate/RTVPhantom/SetNumberOfFrames 401
/gate/RTVPhantom/SetTimePerFrame 1. s

The TACs are tables of the form:

401
0.0 74.100
1. 935.300
2. 507.500
3. 249.300
4. 140.700
5. 98.900
6. 83.200
7. 77.400
...

What I wonder is if it is possible to split the single ROOT output file I
get from the run into several ones according to the time sampling? I.e. one
file containing all coincidences  between 0-1 s, next file coincidences
between 1-2 s and so on. Actually, can you here set any sampling times you
want (that are larger than the timeSampling you defined in your macro and
TACs)? This would be similar to LMF acquisition I guess, but in ROOT format.
Here I could for example split the large ROOT file to one with coincidences
0-10 s, 10-20 s and so on.

I don't really know how to sort the data and match an event to the right
time and so on.
Thank you in advance!
Ida