[Gate-users] optical system threshold/uphold

Tue May 13 05:44:12 CEST 2014

Dear Wang,

I cannot test your macro on vgate. Can you send to me your "Surfaces.xml"
file so that I can run it?

deba999

On Tue, May 13, 2014 at 2:55 AM, Beien Wang <benjaminwang91 at gmail.com>wrote:

> 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
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>

-- 
*Dr Debendranath Sahoo*
Honorary Professor
Physics Department
Sri Sathya Sai Institute of Higher Learning
Prashanthi Nilayam
Puttaparthy
A.P. 515 134
e-mail: deba999 at gmail.com
Mobile: *094945-43026*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20140513/818eba6b/attachment-0001.html>