[Gate-users] How to get crystalSD hits from ASCII output?

[Zhao, Wei]

Dear Jin,

You need a script to retrieve information from the ASCII output or root 
output. Answers regarded to your questions are as follows:

Q1: No. Phantoms do not mean PhantomSD. PhantomSD is used to record 
information that
Q2: No. This is crystalSD.
Q3: To enable hits output, you have to use this command 
(/gate/output/root/setRootHitFlag 1 or 
/gate/output/ascii/setOutFileHitsFlag 1).
Q4: I guess some of the content should be the same, while others are 
different. I never attach these two sensitive detector to the same object.
Q5: The first one may be used to restore interaction information and the 
second one may be used to record deposited energy.
Q6: To obatin counts of each pixel, you have to accumulate the photons 
on each pixel according to the corresponding pixelID.

For more details about how to retrieve information and the difference 
between the phantom SD and crystalSD, the following reference may be 
helpful.
http://iopscience.iop.org/article/10.1088/0031-9155/60/3/1339#citations

Best
Wei

On 10/18/2015 6:00 PM, 
gate-users-request at lists.opengatecollaboration.org wrote:
> Send Gate-users mailing list submissions to
> 	gate-users at lists.opengatecollaboration.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> 	http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
> or, via email, send a message with subject or body 'help' to
> 	gate-users-request at lists.opengatecollaboration.org
>
> You can reach the person managing the list at
> 	gate-users-owner at lists.opengatecollaboration.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Gate-users digest..."
>
>
> Today's Topics:
>
>     1. How to get crystalSD hits from ASCII output? (JinZhongmin)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sun, 18 Oct 2015 03:28:30 -0500
> From: JinZhongmin <jzm87 at live.com>
> To: "gate-users at lists.opengatecollaboration.org"
> 	<gate-users at lists.opengatecollaboration.org>
> Subject: [Gate-users] How to get crystalSD hits from ASCII output?
> Message-ID: <BLU181-W624EFD2D2AD5E9BA73FA90D13B0 at phx.gbl>
> Content-Type: text/plain; charset="gb2312"
>
> Dear Gate user and developers,
>
> I'm learning to use vGate 3.0 and simulate the SPECT system, but I cannot figure out which value is the counts detected by a real detector (crystalSD hits?) from the ASCII output. My simulation is based on the SPECT benchmark example, and both crystal sensitive detector (crystalSD) and phantom sensitive detector (phantomSD) are set as follows.
>
> # Crystal SD
> /gate/crystal/attachCrystalSD
>
> # Phantom SD
> /gate/Phantom/attachPhantomSD
> /gate/movsource/attachPhantomSD
> /gate/table/attachPhantomSD
> /gate/compartment/attachPhantomSD
> /gate/shielding/attachPhantomSD
> /gate/SPECThead/attachPhantomSD
> /gate/collimator/attachPhantomSD
>
> According to the user guide V7.0, hits output file (gateHITs.dat) includes number of Compton and Rayleigh interactions only in phantoms ( Q1: does it mean phantom SD?) before reaching the detector( Q2: does it mean phantom SD, too?), which are shown in columns N+15 and N+16.  I think a Hit stands for a trajectory step within any sensitive detector, why are the hits (interactions in crystals) of crystalSD not given? (Q3)
>
> For singles output file (gatesingles.dat). As I use both crystalSD and phantomSD, and the readout module of digitization is chosen by default, does it mean that I'll get two kinds of singles corresponding to the two separate sensitive detectors after this digitization?(Q4)  According to the user guide V7.0, singles output file includes "Number of Compton interactions in detectors before reaching the detector",  here what's the difference between those two 'detectors'?(Q5) From this output, which value means counts detected by the crystalSD?(Q6)
>
> Also, based on the interfile output of projection set, does the .sin file show us the image contains the real detector counts only?(Q7)
>
> Sorry for so many stupid questions, but I do want to learn from you guys!
>
> Thanks,
> Jin
>   		 	   		
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20151018/0a6ea950/attachment-0001.html>
>
> ------------------------------
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
> End of Gate-users Digest, Vol 113, Issue 16
> *******************************************