[Gate-users] How to extract scatter fraction from the ROOT output?

[Ville-Veikko Wettenhovi]

Hi,

I'm not sure if this is of any help since you use CT-data, but the I've made the following Matlab mex-file(s) for reading PET simulated ROOT data and it can also extract scatter data:
https://github.com/villekf/OMEGA/blob/master/source/GATE_root_matlab_C.cpp for the C-api (crashes Matlab unless used with -nojvm)
https://github.com/villekf/OMEGA/blob/master/source/GATE_root_matlab.cpp for the C++-api (2019a and up).

Regards,
Ville-Veikko Wettenhovi
________________________________
From: Gate-users <gate-users-bounces at lists.opengatecollaboration.org> on behalf of SL H <hsl1196997465 at gmail.com>
Sent: Tuesday, September 22, 2020 11:48
To: Xinjie Cao <xinjie.cao at stonybrook.edu>; gate-users at lists.opengatecollaboration.org <gate-users at lists.opengatecollaboration.org>
Subject: Re: [Gate-users] How to extract scatter fraction from the ROOT output?

Hi, Xinjie.
I have followed your suggestion and converted the ROOT file to a text (it is so large for MATLAB to read) using the code on this page:http://root.cern.ch/phpBB3/viewtopic.php?f=3&t=16590<https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Froot.cern.ch%2FphpBB3%2Fviewtopic.php%3Ff%3D3%26t%3D16590&data=02%7C01%7C%7Cb45d4cb772484c299e2e08d85ed45bb1%7C87879f2e73044bf2baf263e7f83f3c34%7C0%7C1%7C637363613465524116&sdata=Gv6ZODouU9zVzfkcEorF3PK5vkaLlMOJ9B33W4G2%2FpA%3D&reserved=0>.
After that I got something like the attached picture1. I used TBrowser to check my ROOT file and found there are two Singles and two Hits( attached picture2), which one should I convert to text for filtering scatter data?
And in the Gate user doc, I can not find the look up table which you said, there is not much information about the ROOT output, could you please give me some instructions or an example to filter the scatter data?
Thank you!



Xinjie Cao <xinjie.cao at stonybrook.edu<mailto:xinjie.cao at stonybrook.edu>> 于2020年9月22日周二 上午1:58写道：
My suggestion is to convert the ROOT output file into the txt file and filter all scatters you need by MATLAB (with the look up table for ROOT output file format in the GATE instruction file). Then if you still need the ROOT output file you can convert the filtered txt data into ROOT file reversely. We achieved this before, but I am wondering if anyone else has a better solution, like filtering the scatter in ROOT command line work. Good luck!

Best,
Xinjie

On Mon, Sep 21, 2020 at 6:35 AM SL H <hsl1196997465 at gmail.com<mailto:hsl1196997465 at gmail.com>> wrote:
Dear gate users.  I am learning to use GATE for CT simulation, I know that primray+scatter=total (.dat projection). But how can I get the Compton, Rayleigh scatter fraction  from the ROOT output file?
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users<https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.opengatecollaboration.org%2Fmailman%2Flistinfo%2Fgate-users&data=02%7C01%7C%7Cb45d4cb772484c299e2e08d85ed45bb1%7C87879f2e73044bf2baf263e7f83f3c34%7C0%7C1%7C637363613465534108&sdata=HnwZCdo3r4yrflU06BDVauYe4gL%2B42a4APYYvI55haE%3D&reserved=0>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200922/796bff8b/attachment.html>