[Gate-users] optical system threshold/uphold

[Beien Wang]

Dear Gate users,

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.

I use Gate6.2 in vGate. Below is my macro, and the (I think) relevant code
is in bold.

Best,
Beien

#  V I S U A L I Z A T I O N
/vis/open OGLSX
/vis/ogl/set/displayListLimit 8000000
/vis/viewer/reset
/vis/viewer/set/viewpointThetaPhi 90 0
/vis/viewer/zoom 1
/vis/viewer/set/style surface
/vis/drawVolume
/vis/scene/add/trajectories
/vis/scene/add/hits
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate 100
/vis/viewer/update

/gate/geometry/setMaterialDatabase GateMaterials.db

# W O R L D
/gate/world/geometry/setXLength 20. cm
/gate/world/geometry/setYLength 20. cm
/gate/world/geometry/setZLength 20. cm

#=========================
# S Y S T E M
/gate/world/daughters/name OpticalSystem
/gate/world/daughters/insert box
/gate/OpticalSystem/setMaterial Air
/gate/OpticalSystem/geometry/setXLength 10. cm
/gate/OpticalSystem/geometry/setYLength 10. cm
/gate/OpticalSystem/geometry/setZLength 3. cm
/gate/OpticalSystem/placement/setTranslation 0. 0. -1.5 cm
/gate/OpticalSystem/vis/setLineStyle dotted
/gate/OpticalSystem/vis/setColor red
/gate/OpticalSystem/vis/forceWireframe

# S C I N T I L L A T O R  but no scintillation in this case
/gate/OpticalSystem/daughters/name scintillator
/gate/OpticalSystem/daughters/insert box
/gate/scintillator/setMaterial Air
/gate/scintillator/geometry/setXLength 10. cm
/gate/scintillator/geometry/setYLength 10. cm
/gate/scintillator/geometry/setZLength 2. cm
/gate/scintillator/placement/setTranslation 0. 0. 0.5 cm
/gate/scintillator/vis/setColor red
/gate/scintillator/vis/forceSolid
*/gate/scintillator/attachCrystalSD*

# E L E C T R O N I C S
/gate/OpticalSystem/daughters/name electronics
/gate/OpticalSystem/daughters/insert box
/gate/electronics/setMaterial Silicon
/gate/electronics/geometry/setXLength 10. cm
/gate/electronics/geometry/setYLength 10. cm
/gate/electronics/geometry/setZLength 1 cm
/gate/electronics/placement/setTranslation 0. 0. -1 cm
/gate/electronics/vis/setColor yellow
/gate/electronics/vis/forceSolid

*/gate/systems/OpticalSystem/crystal/attach scintillator*
#===================================
*#  P H Y S I C S*
*/gate/physics/addProcess OpticalAbsorption*
*/gate/physics/addProcess OpticalRayleigh*
*/gate/physics/addProcess OpticalBoundary*

*/gate/physics/processList Enabled*
*/gate/physics/processList Initialized*

#=====================================================
# S U R F A C E

/gate/electronics/surfaces/name scintillator_electronics
/gate/electronics/surfaces/insert scintillator
*/gate/electronics/surfaces/scintillator_electronics/setSurface perfect_apd*

/gate/scintillator/surfaces/name electronics_scintillator
/gate/scintillator/surfaces/insert electronics
/gate/scintillator/surfaces/electronics_scintillator/setSurface perfect_apd
#=====================================================
*# D I G I T I Z E R*
*/gate/digitizer/Singles/insert opticaladder*
*/gate/digitizer/Singles/insert readout*
*/gate/digitizer/Singles/readout/setDepth 1*
*/gate/digitizer/Singles/insert thresholder*
*/gate/digitizer/Singles/thresholder/setThreshold 1. eV*
*/gate/digitizer/Singles/insert upholder*
*/gate/digitizer/Singles/upholder/setUphold 10. eV*


#======================================================

# I N I T I A L I Z E
/gate/run/initialize

#===========================================
# D E F I N E  T H E  S O U R C E
#####
/gate/source/addSource opticalsource
/gate/source/opticalsource/gps/type Point
/gate/source/opticalsource/gps/centre 0. 0. -1. cm
/gate/source/opticalsource/gps/particle opticalphoton
/gate/source/opticalsource/gps/polarization   1 0 0
*/gate/source/opticalsource/gps/monoenergy 3.45 eV*
/gate/source/opticalsource/setActivity 3000. Bq
/gate/source/opticalsource/gps/angtype iso
/gate/source/opticalsource/gps/mintheta 0. deg
/gate/source/opticalsource/gps/maxtheta 90. deg
/gate/source/opticalsource/gps/minphi 0. deg
/gate/source/opticalsource/gps/maxphi 360. deg
/gate/source/opticalsource/visualize 1 yellow 1

#===================================================
# O U T P U T
/gate/output/analysis/disable
/gate/output/fastanalysis/enable

/gate/output/projection/enable
/gate/output/projection/setFileName opticalprojection
/gate/output/projection/pixelSizeX       0.02 cm
/gate/output/projection/pixelSizeY       0.02 cm
/gate/output/projection/pixelNumberX     500
/gate/output/projection/pixelNumberY     500
/gate/output/projection/projectionPlane XY

/gate/output/root/enable
/gate/output/root/setFileName     OpticalSimulation
/gate/output/root/setRootOpticalFlag      1
#====================================================

# S T A R T  the A C Q U I S I T I O N
/gate/application/setTimeSlice     0.1  s
/gate/application/setTimeStart     0.   s
/gate/application/setTimeStop      0.1  s

/gate/application/startDAQ
-- 
Beien Wang
Master's student at KTH, Royal Institute of Technology
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20140512/fdd48595/attachment.html>