<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=us-ascii"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="#0563C1" vlink="#954F72"> <div class="WordSection1"> Dear GATE users,<o:p></o:p> <o:p> </o:p> I still have this problem with GATE 7.1, although it was reported (previous message copied below) two years ago. <o:p></o:p> Yannick proposed to attach the phantom’s subvolumes to phantomSD as an alternative (see message below). Can anyone tell me how the subvolume-name is obtained? Such as the “AT_AG_SI1_31” in “/gate/patientAT_AG_SI1_31/attachPhantomSD”?<o:p></o:p> Thanks so much for your help!<o:p></o:p> Rutao<o:p></o:p> Rutao Yao, PhD<o:p></o:p> Professor<o:p></o:p> Dept. of Nuclear Medicine<o:p></o:p> Jacobs School of Medicine and Biomedical Sciences<o:p></o:p> University at Buffalo, State University of New York<o:p></o:p> <o:p> </o:p> =============================================<o:p></o:p> From: Berker, Yannick <<a href="http://gmane.org/get-address.php?address=Yannick.Berker%2d%2f%2ftT4aWmXyZBDgjK7y7TUQ%40public.gmane.org" target="_top">Yannick.Berker@...</a>> Subject: <a href="http://news.gmane.org/find-root.php?message_id=8B967663AF960949A7CDDBD35753F8123676C5D4%40AMSPRD9003MB065.MGDPHG.emi.philips.com" target="_top">Re: No comptonPhantom or RayleighPhantom entries for voxelized phantom</a> Newsgroups: <a href="http://news.gmane.org/gmane.comp.science.opengate.user" target="_top">gmane.comp.science.opengate.user</a> Date: 2014-03-06 13:02:34 GMT (2 years, 7 weeks, 5 days, 15 hours and 13 minutes ago)<o:p></o:p> Hi Tom,<o:p></o:p> which version of GATE are you using? I had my problems with early versions of GATE 6 (like 6.0 or 6.1), and I have managed to get seemingly correct values with GATE 6.2.<o:p></o:p> <o:p></o:p> First, I recall that for Rayleigh (not Compton) interactions, there was a bug in GATE 6.0 and 6.1 that was fixed in 6.2. In 6.0 and 6.1, an interaction from a “PenelopeRayleighScattering” (6.0/6.1) or “RayleighScattering” (6.1) process would not be counted as a Rayleigh interaction in the ROOT output since the code checks for a process name with “LowEnRayleigh” in it – this is due to these lines here:<o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.0.0/source/digits_hits/src/GateAnalysis.cc#L278">http://www.opengatecollaboration.org/lxr/source/6.0.0/source/digits_hits/src/GateAnalysis.cc#L278</a><o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.1.0/source/digits_hits/src/GateAnalysis.cc#L281">http://www.opengatecollaboration.org/lxr/source/6.1.0/source/digits_hits/src/GateAnalysis.cc#L281</a><o:p></o:p> <o:p></o:p> Here, is it fixed:<o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.2.0/source/digits_hits/src/GateAnalysis.cc#L292">http://www.opengatecollaboration.org/lxr/source/6.2.0/source/digits_hits/src/GateAnalysis.cc#L292</a><o:p></o:p> <o:p></o:p> Here are some physics processes which could have been affected by that bug:<o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.0.0/source/physics/src/GateRayleighPenelopePB.cc">http://www.opengatecollaboration.org/lxr/source/6.0.0/source/physics/src/GateRayleighPenelopePB.cc</a><o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.1.0/source/physics/src/GateRayleighPenelopePB.cc">http://www.opengatecollaboration.org/lxr/source/6.1.0/source/physics/src/GateRayleighPenelopePB.cc</a><o:p></o:p> <a href="http://www.opengatecollaboration.org/lxr/source/6.1.0/source/physics/src/GateRayleighPB.cc">http://www.opengatecollaboration.org/lxr/source/6.1.0/source/physics/src/GateRayleighPB.cc</a><o:p></o:p> <o:p></o:p> If you are already using GATE 6.2, the problem (or at least one solution) might be that you can attach the PhantomSD to the subvolumes which are created from your voxelized phantom. My current code contains this:<o:p></o:p> <o:p></o:p> /gate/patientAT_AG_SI1_31/attachPhantomSD<o:p></o:p> […, 4 more]<o:p></o:p> /gate/patientAT_AG_SI5_36/attachPhantomSD<o:p></o:p> /gate/patientAir_0/attachPhantomSD<o:p></o:p> /gate/patientAir_1/attachPhantomSD<o:p></o:p> /gate/patientConnectiveTissue_40/attachPhantomSD<o:p></o:p> /gate/patientLung_10/attachPhantomSD<o:p></o:p> […, 27 more]<o:p></o:p> /gate/patientLung_9/attachPhantomSD<o:p></o:p> /gate/patientMarrow_Bone01_41/attachPhantomSD<o:p></o:p> […, 13 more]<o:p></o:p> /gate/patientMarrow_Bone08_55/attachPhantomSD<o:p></o:p> /gate/patientSoftTissus_37/attachPhantomSD<o:p></o:p> /gate/patientSoftTissus_38/attachPhantomSD<o:p></o:p> /gate/patientSoftTissus_39/attachPhantomSD<o:p></o:p> /gate/patientworldDefaultAir/attachPhantomSD<o:p></o:p> <o:p></o:p> The complete list depends on your phantom – you should include all subvolumes, but you may only include existing ones. I think the list of volumes generated is part of the GATE output, so you can grab this from your current simulation, prepend “/gate/patient” (or whatever your generating volume name is), append “/attachPhantomSD” and save this as the phantomSD.mac specific to that voxelized input and the volume-generation parameters. Should any of these change without yout re-generating the mac file, in the best case you lose one of the volumes you had in your the mac file before and your script returns an error message. In the worst case, you will get inaccurate results without noticing, for example, when only additional volumes are generated and you miss scatter interactions in those without getting any error message.<o:p></o:p> <o:p></o:p> Best regards<o:p></o:p> <o:p></o:p> Yannick<o:p></o:p> <o:p></o:p> From: Tom Buitenhuis [mailto:h.j.t.buitenhuis-39IHFo8E5E0@public.gmane.org] Sent: Thursday, March 06, 2014 12:27 PM To: Berker, Yannick; gate-users-ZwoEplunGu2uWGBtOU8ma0fXCdDtWGEhSoIsB4E12gc@public.gmane.org Subject: No comptonPhantom or RayleighPhantom entries for voxelized phantom<o:p></o:p> <o:p></o:p> Dear GATE users,<o:p></o:p> <o:p></o:p> I am having the same problem as Yannick Berker described on 2011-10-21. I use the ImageNestedParametrisedVolume phantom for a PET simulation. The sources are described using interfiles. Everything seems to be working correctly, i.e. the sources and phantom are in the right place and coincidences are detected. The effect of phantom attenuation is also clearly visible.<o:p></o:p> <o:p></o:p> The problem that I’m having, is that in the ROOT coincidences tree there is no comptonPhantom and RayleighPhantom data. I cannot add the VoxelPhantomSD to the ImageNestedParametrisedVolume, because I get an error message saying the command is unavailable. I can add the ordinary PhantomSD, but this does not seem to add the comptonPhantom or RayleighPhantom information to the coincidence tree. The comptonCrystal and RayleighCrystal entries are correct, so I think the physics processes are registered and working correctly. <o:p></o:p> <o:p></o:p> I need the comptonPhantom and RayleighPhantom information to get scatter information on coincidences.<o:p></o:p> <o:p></o:p> Did anyone else come across this problem, and if so, how did you solve it? Or can anyone point me in the direction of a solution?<o:p></o:p> <o:p></o:p> Regards,<o:p></o:p> <o:p></o:p> Tom Buitenhuis<o:p></o:p> <o:p></o:p> PhD student<o:p></o:p> KVI - Center for Advanced Radiation Technology<o:p></o:p> University of Groningen<o:p></o:p> The Netherlands<o:p></o:p> <o:p> </o:p> Dear GATE users,<o:p></o:p> <o:p></o:p> I have some issues using the HounsfieldMaterialGenerator to simulate a PET scan of a patient. The<o:p></o:p> simulation itself is running fine, but I am lacking two features that I am used to from other simulations,<o:p></o:p> which are scatter information and color visualization. I mainly used the example from<o:p></o:p> <a href="http://www.creatis.insa-lyon.fr/gate/Example2/" target="_top">http://www.creatis.insa-lyon.fr/gate/Example2/</a> to generate the following phantom:<o:p></o:p> <o:p></o:p> /gate/HounsfieldMaterialGenerator/SetMaterialTable CT/Schneider2000MaterialsTable.txt<o:p></o:p> /gate/HounsfieldMaterialGenerator/SetDensityTable CT/Schneider2000DensitiesTable.txt<o:p></o:p> /gate/HounsfieldMaterialGenerator/SetDensityTolerance 0.03 g/cm3<o:p></o:p> /gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename CT/p971-HUmaterials.db<o:p></o:p> /gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename CT/p971-HU2mat.txt<o:p></o:p> /gate/HounsfieldMaterialGenerator/Generate<o:p></o:p> /gate/world/daughters/name patient<o:p></o:p> /gate/world/daughters/insert ImageRegionalizedVolume<o:p></o:p> /gate/geometry/setMaterialDatabase CT/p971-HUmaterials.db<o:p></o:p> /gate/patient/geometry/SetHUToMaterialFile CT/p971-HU2mat.txt<o:p></o:p> /gate/patient/geometry/SetImage CT/p971s3.hdr<o:p></o:p> ...<o:p></o:p> <o:p></o:p> The ROOT output shows me that the volume is placed correctly, and the distribution of detected<o:p></o:p> coincidences also shows the effect of photon attenuation by the phantom. However, I cannot attach a<o:p></o:p> VoxelPhantomSD to the patient - I receive an error message telling me that this command is unavailable:<o:p></o:p> [G4-cerr] ***** COMMAND NOT FOUND </gate/patient/attachVoxelPhantomSD> *****<o:p></o:p> Attaching an ordinary PhantomSD seems to work, but only generates 0 entries in ComptonPhantom and<o:p></o:p> RayleighPhantom. I also tried using a ImageNestedParametrisedVolume, with no effect. ComptonCrystal<o:p></o:p> entries, however, seem to be correct (>= 0).<o:p></o:p> <o:p></o:p> Apart from that, how can I visualize the voxelized phantom using colors? I was not able to add color<o:p></o:p> attributes (like "true 0.5 0.5 0.5 1.0") to p971-HU2mat.txt manually, since this resulted in the<o:p></o:p> following error message:<o:p></o:p> src/<a href="http://gatevimagevolume.cc/">GateVImageVolume.cc</a> (l.332): No Hounsfield material defined in the file<o:p></o:p> CT/p971s3_ToSim_OneBP_HU2mat_Colors.txt. Abort<o:p></o:p> <o:p></o:p> Does anyone have knowledge in this area, or can point to some helpful documentation on this topic? The GATE<o:p></o:p> User's Guide has a few, but only incomplete examples for the case of voxelized phantoms and sources.<o:p></o:p> <o:p></o:p> Best regards,<o:p></o:p> <o:p></o:p> Yannick Berker<o:p></o:p> <o:p> </o:p> </div> </body> </html>