<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div><br>
</div>
Dear Debora Lamia,
<div><br>
</div>
<div>I understood, that you are referring to the PET_Analyse.C in directory</div>
<div>example_ROOT_analyse!? The commented lines for Singles Hits Gate are</div>
<div>only useful, when using this module for analysing a Gate run, where you</div>
<div>have filled in these branches. </div>
<div><br>
</div>
<div>Hence, if you did NOT add lines like</div>
<div><br>
</div>
<div>
<div>/gate/digitizer/Singles/insert adder</div>
<div>/gate/digitizer/Singles/insert readout</div>
<div>/gate/digitizer/Singles/readout/setDepth 1</div>
<div><br>
</div>
<div>then you cannot expect to find singles in a singles branch of the ROOT-file.</div>
<div><br>
</div>
<div>Also you need lines like </div>
<div><br>
</div>
<div>
<div>#<span class="Apple-tab-span" style="white-space:pre"> </span>O U T P U T </div>
<div><br>
</div>
<div>/gate/output/root/enable</div>
<div>/gate/output/root/setFileName benchmarkPET</div>
<div>/gate/output/root/setRootHitFlag 0</div>
<div>/gate/output/root/setRootSinglesFlag 1</div>
<div>/gate/output/root/setRootCoincidencesFlag 1</div>
<div>/gate/output/root/setRootdelayFlag 1</div>
</div>
<div><br>
</div>
<div>to direct GATE to write the respective information into the ROOT-file</div>
<div><br>
</div>
<div>It may be helpful for you to also look into the benchmarks to get</div>
<div>an idea how the instructions to GATE match the analysis with</div>
<div>the root-macro.</div>
<div><br>
</div>
<div>Concerning the randoms you have to increase the activity considerably,</div>
<div>but in order to avoid long simulation times you may at the same time reduce the</div>
<div>total simulation duration like:</div>
<div><br>
</div>
<div>
<div>#<span class="Apple-tab-span" style="white-space:pre"> </span>S T A R T </div>
<div>/gate/application/setTimeSlice 0.1 s</div>
<div><br>
</div>
<div>#or even a lower value depending on the activity, like</div>
<div>/gate/application/setTimeSlice 0.01 s</div>
<div><br>
</div>
<div>/gate/application/setTimeStart 0. s</div>
<div>/gate/application/setTimeStop 1. s</div>
<div><br>
</div>
<div>/gate/application/startDAQ</div>
</div>
<div><br>
</div>
<div>in order to see the progress of the simulation!</div>
<div><br>
</div>
<div>Hope this helps!</div>
<div><br>
</div>
<div>Kind regards,</div>
<div><br>
</div>
<div>Uwe Pietrzyk</div>
<div><br>
</div>
<div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div><span style="orphans: 2; widows: 2;">/----------------------------------------------------------------------</span></div>
<div style="orphans: 2; widows: 2;">Prof. Dr. Uwe Pietrzyk<br>
Diplom-Physiker<br>
Institut für Neurowissenschaften und Medizin / INM-4<br>
Forschungszentrum Jülich GmbH<br>
52425 Jülich<br>
Tel: +49-2461-61 4763<br>
Fax: +49-2461-61 2820<br>
Mobil: +49-162-103-4871<br>
eMail: <a href="mailto:U.Pietrzyk@fz-juelich.de">U.Pietrzyk@fz-juelich.de</a><br>
<a href="http://www.fz-juelich.de/inm/inm-4/DE/Home/home_node.html">http://www.fz-juelich.de/inm/inm-4/DE/Home/home_node.html</a><br>
&<br>
Bergische Universität Wuppertal<br>
Fakultät für Mathematik und Naturwissenschaften</div>
<div style="orphans: 2; widows: 2;">Büro: D.08.01<br>
Gaußstr. 20<br>
42097 Wuppertal<br>
Tel: +49-(0)202-439-3523<br>
Fax: +49-(0)202-439-2811<br>
eMail: <a href="mailto:Uwe.Pietrzyk@physik.uni-wuppertal.de">Uwe.Pietrzyk@uni-wuppertal.de</a><br>
<a href="http://www.medizinphysik.uni-wuppertal.de/">http://www.medizinphysik.uni-wuppertal.de</a><br>
----------------------------------------------------------------------/</div>
</div>
</div>
<br>
<div>
<div>On 13 Jun 2016, at 12:41, Debora Lamia <<a href="mailto:debora.lamia@ibfm.cnr.it">debora.lamia@ibfm.cnr.it</a>> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">Dear Prof. Pietrzyk,<br>
<br>
I am an Italian research fellow at the institute of molecular bioimaging and physiology (IBFM) of the National Research Council.<br>
I am a physicist and I'm dealing with medical physics research, in particular dosimetry and Geant4 Monte Carlo simulations for radiotherapy applications.<br>
Lately I've started to use GATE to simulate the PET system.<br>
I would extract the NEC curves, so I would obtain the trues, scattered and random events from the simulations.<br>
Please, I would ask you some information about the PET simulation analysis.<br>
I used your file PET_Analyse.C to analyze my root file obtained carrying out example_PET GATE simulations.<br>
I set activity until 3 kBq. I obtained prompt and trues counts, but the random counts are results zero. I detected only one random count with activity 3 kBq.<br>
I know that the random counts increase with the square of the activity.<br>
My question is about the analysis of the Singles, Hits and Gate (24, 25 and 26 lines of the PET_Analyse.C).<br>
In your file these lines were commented, but if I want detect the random counts, do I also have to analyze singles events?<br>
<br>
Do I have to write a similar analysis file, adding the variables of the singles, hits and gate?<br>
<br>
<br>
Thank you very much for your help<br>
<br>
Best regards,<br>
<br>
Debora Lamia<br>
<br>
<br>
-- <br>
Debora Lamia<br>
Research Fellow @ IBFM - CNR UOS Cefalù<br>
c/da Pietrapollastra Pisciotto, 90015 Cefalù (PA)<br>
tel. +39 0921 920707<br>
<br>
</blockquote>
</div>
<br>
</div>
<br>
<font face="Arial" color="Black" size="1"><br>
------------------------------------------------------------------------------------------------<br>
------------------------------------------------------------------------------------------------<br>
Forschungszentrum Juelich GmbH<br>
52425 Juelich<br>
Sitz der Gesellschaft: Juelich<br>
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498<br>
Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher<br>
Geschaeftsfuehrung: Prof. Dr.-Ing. Wolfgang Marquardt (Vorsitzender),<br>
Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,<br>
Prof. Dr. Sebastian M. Schmidt<br>
------------------------------------------------------------------------------------------------<br>
------------------------------------------------------------------------------------------------<br>
<br>
</font>
</body>
</html>