<div dir="ltr">Dear Gate users,<div><br></div><div>I am simulating an OpticalSystem in Gate. I need to set a threshold and an uphold to the digitizer. However, I found that no matter what value I set as optical photon energy threshold, this threshold does not work (all photon are detected even I set the threshold to an impossibly high value). Moreover, no matter what value I set as uphold, it will result in no detection of photon at all in the final projection output. I cannot understand this. Can anyone help me? Thank you in advance.</div>
<div><br></div><div>I use Gate6.2 in vGate. Below is my macro, and the (I think) relevant code is in bold.</div><div><br></div><div>Best,</div><div>Beien<br><div><br></div><div><div># V I S U A L I Z A T I O N</div><div>
/vis/open OGLSX</div><div>/vis/ogl/set/displayListLimit 8000000</div><div>/vis/viewer/reset</div><div>/vis/viewer/set/viewpointThetaPhi 90 0</div><div>/vis/viewer/zoom 1</div><div>/vis/viewer/set/style surface</div><div>/vis/drawVolume</div>
<div>/vis/scene/add/trajectories</div><div>/vis/scene/add/hits</div><div>/tracking/storeTrajectory 1</div><div>/vis/scene/endOfEventAction accumulate 100</div><div>/vis/viewer/update</div><div><br></div><div>/gate/geometry/setMaterialDatabase GateMaterials.db</div>
<div><br></div><div># W O R L D </div><div>/gate/world/geometry/setXLength 20. cm </div><div>/gate/world/geometry/setYLength 20. cm </div><div>/gate/world/geometry/setZLength 20. cm</div><div><br></div><div>#=========================</div>
<div># S Y S T E M</div><div>/gate/world/daughters/name OpticalSystem</div><div>/gate/world/daughters/insert box</div><div>/gate/OpticalSystem/setMaterial Air</div><div>/gate/OpticalSystem/geometry/setXLength 10. cm</div>
<div>/gate/OpticalSystem/geometry/setYLength 10. cm</div><div>/gate/OpticalSystem/geometry/setZLength 3. cm</div><div>/gate/OpticalSystem/placement/setTranslation 0. 0. -1.5 cm </div><div>/gate/OpticalSystem/vis/setLineStyle dotted</div>
<div>/gate/OpticalSystem/vis/setColor red</div><div>/gate/OpticalSystem/vis/forceWireframe</div><div><br></div><div># S C I N T I L L A T O R but no scintillation in this case</div><div>/gate/OpticalSystem/daughters/name scintillator</div>
<div>/gate/OpticalSystem/daughters/insert box</div><div>/gate/scintillator/setMaterial Air</div><div>/gate/scintillator/geometry/setXLength 10. cm</div><div>/gate/scintillator/geometry/setYLength 10. cm</div><div>/gate/scintillator/geometry/setZLength 2. cm</div>
<div>/gate/scintillator/placement/setTranslation 0. 0. 0.5 cm</div><div>/gate/scintillator/vis/setColor red</div><div>/gate/scintillator/vis/forceSolid</div><div><b>/gate/scintillator/attachCrystalSD</b></div><div><br></div>
<div># E L E C T R O N I C S </div><div>/gate/OpticalSystem/daughters/name electronics</div><div>/gate/OpticalSystem/daughters/insert box</div><div>/gate/electronics/setMaterial Silicon</div><div>/gate/electronics/geometry/setXLength 10. cm</div>
<div>/gate/electronics/geometry/setYLength 10. cm</div><div>/gate/electronics/geometry/setZLength 1 cm</div><div>/gate/electronics/placement/setTranslation 0. 0. -1 cm</div><div>/gate/electronics/vis/setColor yellow</div>
<div>/gate/electronics/vis/forceSolid</div><div><br></div><div><b>/gate/systems/OpticalSystem/crystal/attach scintillator</b></div><div>#===================================</div><div><b># P H Y S I C S</b></div><div><b>/gate/physics/addProcess OpticalAbsorption</b></div>
<div><b>/gate/physics/addProcess OpticalRayleigh</b></div><div><b>/gate/physics/addProcess OpticalBoundary</b></div><div><b><br></b></div><div><b>/gate/physics/processList Enabled</b></div><div><b>/gate/physics/processList Initialized</b></div>
<div><br></div><div>#=====================================================</div><div># S U R F A C E</div><div><br></div><div>/gate/electronics/surfaces/name scintillator_electronics</div><div>/gate/electronics/surfaces/insert scintillator</div>
<div><b>/gate/electronics/surfaces/scintillator_electronics/setSurface perfect_apd</b></div><div><b><br></b></div><div>/gate/scintillator/surfaces/name electronics_scintillator</div><div>/gate/scintillator/surfaces/insert electronics</div>
<div>/gate/scintillator/surfaces/electronics_scintillator/setSurface perfect_apd</div><div>#=====================================================</div><div><b># D I G I T I Z E R</b></div><div><b>/gate/digitizer/Singles/insert opticaladder</b></div>
<div><b>/gate/digitizer/Singles/insert readout</b></div><div><b>/gate/digitizer/Singles/readout/setDepth 1</b></div><div><b>/gate/digitizer/Singles/insert thresholder</b></div><div><b>/gate/digitizer/Singles/thresholder/setThreshold 1. eV</b></div>
<div><b>/gate/digitizer/Singles/insert upholder</b></div><div><b>/gate/digitizer/Singles/upholder/setUphold 10. eV</b></div><div><br></div><div><br></div><div>#======================================================</div><div>
<br></div><div># I N I T I A L I Z E </div><div>/gate/run/initialize</div><div><br></div><div>#===========================================</div><div># D E F I N E T H E S O U R C E</div><div>#####</div><div>/gate/source/addSource opticalsource</div>
<div>/gate/source/opticalsource/gps/type Point</div><div>/gate/source/opticalsource/gps/centre 0. 0. -1. cm</div><div>/gate/source/opticalsource/gps/particle opticalphoton</div><div>/gate/source/opticalsource/gps/polarization 1 0 0</div>
<div><b>/gate/source/opticalsource/gps/monoenergy 3.45 eV</b></div><div>/gate/source/opticalsource/setActivity 3000. Bq</div><div>/gate/source/opticalsource/gps/angtype iso</div><div>/gate/source/opticalsource/gps/mintheta 0. deg </div>
<div>/gate/source/opticalsource/gps/maxtheta 90. deg </div><div>/gate/source/opticalsource/gps/minphi 0. deg </div><div>/gate/source/opticalsource/gps/maxphi 360. deg </div><div>/gate/source/opticalsource/visualize 1 yellow 1</div>
<div><br></div><div>#===================================================</div><div># O U T P U T</div><div>/gate/output/analysis/disable</div><div>/gate/output/fastanalysis/enable</div><div><br></div><div>/gate/output/projection/enable</div>
<div>/gate/output/projection/setFileName opticalprojection</div><div>/gate/output/projection/pixelSizeX <span class="" style="white-space:pre"> </span> 0.02 cm</div><div>/gate/output/projection/pixelSizeY <span class="" style="white-space:pre"> </span> 0.02 cm</div>
<div>/gate/output/projection/pixelNumberX <span class="" style="white-space:pre"> </span> 500</div><div>/gate/output/projection/pixelNumberY <span class="" style="white-space:pre"> </span> 500</div><div>/gate/output/projection/projectionPlane XY</div>
<div><br></div><div>/gate/output/root/enable</div><div>/gate/output/root/setFileName OpticalSimulation</div><div>/gate/output/root/setRootOpticalFlag 1</div><div>#====================================================</div>
<div><br></div><div># S T A R T the A C Q U I S I T I O N</div><div>/gate/application/setTimeSlice 0.1 s</div><div>/gate/application/setTimeStart 0. s</div><div>/gate/application/setTimeStop 0.1 s</div><div>
<br></div><div>/gate/application/startDAQ</div></div>-- <br><div dir="ltr"><font face="courier new, monospace">Beien Wang</font><div><font face="courier new, monospace">Master's student at KTH, Royal Institute of Technology</font></div>
<div><br></div></div>
</div></div>