<div dir="ltr">Hello, <div><br></div><div>Disclaimer: I am not used to tracking of optical photons, probably more experienced users can bring more advice. </div><div><br></div><div>However, conventional PET simulation with GATE does NOT track optical photons. In the scintillator crystal, the deposited energy is considered as a surrogate of the number of optical photons (that are not produced, not tracked), eventually weighted by some efficiency factor. If you plan a more deep simulation that really tracks optical photons in such a large system than a PET, be prepare to face enormous computation time. Simulation with optical photons is, to my knowledge, limited to the characterisation of one single detector.</div><div><br></div><div>HTH, </div><div>David</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Nov 1, 2019 at 3:25 PM Christodoulou, Theodoros (2017) <<a href="mailto:Theodoros.Christodoulou.2017@live.rhul.ac.uk">Theodoros.Christodoulou.2017@live.rhul.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Dear D.Sarrut,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I have a question related to GATE. I want to simulate a PET scanner but instead of using the standard scintillation mediums (LSO e.t.c) I want to use liquid Argon. Therefore I have to attach the scintillation properties of <span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">liquid
argon to G4_lAr material. </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Moreover, the vessel which will contain liquid argon will be stainless_steel and I want to add a reflective coating on it so the scintillation photons will
bounce inside the vessel and then detected by the sensitive detector. </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Can Gate do this thing? I have managed to do that in a Geant4 simulation I wrote for a single cell but I want to use Gate for a complete detector.</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Best</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Theo </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Example of G4 code I wrote and tested for the two question I have above. </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"></span></div>
<div style="color:rgb(0,0,0);background-color:rgb(255,255,254);font-family:SFMono-Medium,"SF Mono","Segoe UI Mono","Roboto Mono","Ubuntu Mono",Menlo,monospace;font-weight:normal;font-size:13px;line-height:20px">
<div><span style="color:rgb(160,160,160);font-style:italic">/Materials</span></div>
<div><span> </span><span style="color:rgb(160,160,160);font-style:italic">//Liquid Argon</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">G4NistManager</span><span> *</span><span style="color:rgb(32,32,32)">mat</span><span> =
</span><span style="color:rgb(32,32,32)">G4NistManager</span><span>::</span><span style="color:rgb(32,32,32)">Instance</span><span>();</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">G4Material</span><span> *</span><span style="color:rgb(32,32,32)">LAr</span><span> =
</span><span style="color:rgb(32,32,32)">mat</span><span>-></span><span style="color:rgb(32,32,32)">FindOrBuildMaterial</span><span>(</span><span style="color:rgb(191,38,0)">"G4_lAr"</span><span>) ;</span></div>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<div style="color:rgb(0,0,0);background-color:rgb(255,255,254);font-family:SFMono-Medium,"SF Mono","Segoe UI Mono","Roboto Mono","Ubuntu Mono",Menlo,monospace;font-weight:normal;font-size:13px;line-height:20px">
<div><span style="color:rgb(160,160,160);font-style:italic">//LAr properties</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">G4double</span><span> </span><span style="color:rgb(32,32,32)">lar_Energy</span><span>[]= {
</span><span style="color:rgb(101,84,192)">20.60</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,
</span><span style="color:rgb(101,84,192)">15.5</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,
</span><span style="color:rgb(101,84,192)">11.92</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">11.2</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">10.87</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">10.69</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">10.5</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">10.3</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">10.13</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span></div>
<div><span> </span><span style="color:rgb(101,84,192)">9.97</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.81</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.73</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.66</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.58</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.51</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.44</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.36</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.22</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">9.09</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">8.96</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span></div>
<div><span> </span><span style="color:rgb(101,84,192)">8.83</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">8.70</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">8.52</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">8.35</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">8.13</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">7.02</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>,</span><span style="color:rgb(101,84,192)">1.55</span><span>*</span><span style="color:rgb(32,32,32)">eV</span><span>};</span></div>
</div>
<br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<div style="color:rgb(0,0,0);background-color:rgb(255,255,254);font-family:SFMono-Medium,"SF Mono","Segoe UI Mono","Roboto Mono","Ubuntu Mono",Menlo,monospace;font-weight:normal;font-size:13px;line-height:20px">
<div><span style="color:rgb(32,32,32)">G4MaterialPropertiesTable</span><span> *</span><span style="color:rgb(32,32,32)">LAr_mt</span><span> =
</span><span style="color:rgb(9,30,66);font-weight:bold">new</span><span> </span><span style="color:rgb(32,32,32)">G4MaterialPropertiesTable</span><span>();</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddProperty</span><span>(</span><span style="color:rgb(191,38,0)">"FASTCOMPONENT"</span><span>,
</span><span style="color:rgb(32,32,32)">lar_Energy</span><span>, </span><span style="color:rgb(32,32,32)">lar_fSCINT</span><span>,
</span><span style="color:rgb(32,32,32)">larnum</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddProperty</span><span>(</span><span style="color:rgb(191,38,0)">"SLOWCOMPONENT"</span><span>,
</span><span style="color:rgb(32,32,32)">lar_Energy</span><span>, </span><span style="color:rgb(32,32,32)">lar_sSCINT</span><span>,
</span><span style="color:rgb(32,32,32)">larnum</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddProperty</span><span>(</span><span style="color:rgb(191,38,0)">"RINDEX"</span><span>, </span><span style="color:rgb(32,32,32)">lar_Energy</span><span>,
</span><span style="color:rgb(32,32,32)">lar_RIND</span><span>, </span><span style="color:rgb(32,32,32)">larnum</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddProperty</span><span>(</span><span style="color:rgb(191,38,0)">"ABSLENGTH"</span><span>,
</span><span style="color:rgb(32,32,32)">lar_Energy</span><span>, </span><span style="color:rgb(32,32,32)">lar_ABSL</span><span>, </span><span style="color:rgb(32,32,32)">larnum</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddConstProperty</span><span>(</span><span style="color:rgb(191,38,0)">"SCINTILLATIONYIELD"</span><span>,</span><span style="color:rgb(101,84,192)">40</span><span>./</span><span style="color:rgb(32,32,32)">keV</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddConstProperty</span><span>(</span><span style="color:rgb(191,38,0)">"RESOLUTIONSCALE"</span><span>,</span><span style="color:rgb(101,84,192)">0.11</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddConstProperty</span><span>(</span><span style="color:rgb(191,38,0)">"FASTTIMECONSTANT"</span><span>,</span><span style="color:rgb(101,84,192)">6</span><span>.*</span><span style="color:rgb(32,32,32)">ns</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddConstProperty</span><span>(</span><span style="color:rgb(191,38,0)">"SLOWTIMECONSTANT"</span><span>,</span><span style="color:rgb(101,84,192)">1600</span><span>*</span><span style="color:rgb(32,32,32)">ns</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr_mt</span><span>-></span><span style="color:rgb(32,32,32)">AddConstProperty</span><span>(</span><span style="color:rgb(191,38,0)">"YIELDRATIO"</span><span>,</span><span style="color:rgb(101,84,192)">0.67</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LAr</span><span>-></span><span style="color:rgb(32,32,32)">SetMaterialPropertiesTable</span><span>(</span><span style="color:rgb(32,32,32)">LAr_mt</span><span>);</span></div>
</div>
<br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">For stainless steel, I created a reflective surface between steel and LAr</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"></span></div>
<div style="color:rgb(0,0,0);background-color:rgb(255,255,254);font-family:SFMono-Medium,"SF Mono","Segoe UI Mono","Roboto Mono","Ubuntu Mono",Menlo,monospace;font-weight:normal;font-size:13px;line-height:20px">
<div><span style="color:rgb(32,32,32)">G4OpticalSurface</span><span>* </span><span style="color:rgb(32,32,32)">OpLarSurface</span><span> =
</span><span style="color:rgb(9,30,66);font-weight:bold">new</span><span> </span><span style="color:rgb(32,32,32)">G4OpticalSurface</span><span>(</span><span style="color:rgb(191,38,0)">"LarSurface"</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">G4LogicalBorderSurface</span><span>*
</span><span style="color:rgb(32,32,32)">LarSurface</span><span> = </span><span style="color:rgb(9,30,66);font-weight:bold">new</span><span> </span><span style="color:rgb(32,32,32)">G4LogicalBorderSurface</span><span> (</span><span style="color:rgb(191,38,0)">"LarSurface"</span><span>,</span><span style="color:rgb(32,32,32)">PVLAr</span><span>,</span><span style="color:rgb(32,32,32)">PVsteel</span><span>,</span><span style="color:rgb(32,32,32)">OpLarSurface</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">OpLarSurface</span><span>-></span><span style="color:rgb(32,32,32)">SetType</span><span>(</span><span style="color:rgb(32,32,32)">dielectric_metal</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">OpLarSurface</span><span>-></span><span style="color:rgb(32,32,32)">SetFinish</span><span>(</span><span style="color:rgb(32,32,32)">polished</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">OpLarSurface</span><span>-></span><span style="color:rgb(32,32,32)">SetModel</span><span>(</span><span style="color:rgb(32,32,32)">glisur</span><span>);</span></div>
<div><span> </span><span style="color:rgb(32,32,32)">LarSurface</span><span>-></span><span style="color:rgb(32,32,32)">DumpInfo</span><span>();</span></div>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"> </span></div>
</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">David Sarrut, Phd<br>Directeur de recherche CNRS<br>CREATIS, UMR CNRS 5220, Inserm U1206<div>Centre de lutte contre le cancer Léon Bérard<br>28 rue Laënnec, 69373 Lyon cedex 08<br>Tel : 04 78 78 51 51 / 06 74 72 05 42<br><a href="http://www.creatis.insa-lyon.fr/~dsarrut" target="_blank">http://www.creatis.insa-lyon.fr/~dsarrut</a><br>_________________________________</div><div> "2 + 2 = 5, for extremely large values of 2"<br>_________________________________</div></div></div></div></div></div></div>