[Gate-users] Problem with detection of optical photons
Nuray Yavuzkanat
ny535 at york.ac.uk
Tue Dec 2 13:33:27 CET 2014
Dear Gate users,
I am simulating an OpticalSystem in Gate6.2 for obtaining number of
detected photon on scintillator volume. I changed many parameters such as;
scintillation material which has already defined in Materials.xml file,
source etc. but nothing change in Root output file. I can't see any photon
produced by plastic scintillator, no detection of photon. Also I simulated
OpticalSystem examples and still I can't see any photon detected. Can
anyone help me in advance ? Thank you .
Below is my macro.
Best Regards
Nuray
#=====================================================
# VISUALISATION
#=====================================================
/vis/open OGLSQt
/vis/viewer/set/viewpointThetaPhi 30 60
/vis/viewer/zoom 4
#/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 10 cm
/gate/world/geometry/setYLength 10 cm
/gate/world/geometry/setZLength 10 cm
/gate/world/setMaterial Air
#=========================================================
# G E O M E T R Y
#=========================================================
# === optical system ===
/gate/world/daughters/name OpticalSystem
/gate/world/daughters/insert box
/gate/OpticalSystem/geometry/setXLength 10. mm
/gate/OpticalSystem/geometry/setYLength 7.0 mm
/gate/OpticalSystem/geometry/setZLength 7.0 mm
/gate/OpticalSystem/placement/setTranslation 0 0 0 mm
/gate/OpticalSystem/setMaterial Air
/gate/OpticalSystem/vis/setColor white
/gate/OpticalSystem/vis/forceWireframe
/gate/OpticalSystem/vis/setVisible 1
# === crystal ===
/gate/OpticalSystem/daughters/name scintillator
/gate/OpticalSystem/daughters/insert box
/gate/scintillator/geometry/setXLength 3 mm
/gate/scintillator/geometry/setYLength 7 mm
/gate/scintillator/geometry/setZLength 7 mm
/gate/scintillator/placement/setTranslation -1.5 0 0 mm
/gate/scintillator/setMaterial Scinti-detector-C10H11
/gate/scintillator/vis/setColor gray
/gate/scintillator/vis/forceSolid
/gate/scintillator/attachCrystalSD
# === electronic board ===
/gate/OpticalSystem/daughters/name Electronics
/gate/OpticalSystem/daughters/insert box
/gate/Electronics/geometry/setXLength 2 mm
/gate/Electronics/geometry/setYLength 7 mm
/gate/Electronics/geometry/setZLength 7 mm
/gate/Electronics/setMaterial Air
/gate/Electronics/placement/setTranslation -4 0 0 mm
/gate/Electronics/vis/setColor blue
/gate/Electronics/vis/forceSolid
/gate/systems/OpticalSystem/crystal/attach scintillator
# ===surface=====
/gate/scintillator/surfaces/name Detection1
/gate/scintillator/surfaces/insert Electronics
/gate/scintillator/surfaces/Detection1/setSurface perfect_apd
#/gate/Electronics/surfaces/name Detection2
#/gate/Electronics/surfaces/insert scintillator
#/gate/Electronics/surfaces/Detection2/setSurface perfect_apd
#=========================================================
# P H Y S I C S
#=========================================================
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel PenelopeModel
/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel PenelopeModel
/gate/physics/addProcess RayleighScattering
/gate/physics/processes/RayleighScattering/setModel PenelopeModel
/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+
/gate/physics/addProcess OpticalAbsorption
/gate/physics/addProcess PositronAnnihilation
/gate/physics/addProcess Scintillation
/gate/physics/addProcess OpticalBoundary
/gate/physics/addProcess OpticalRayleigh
/gate/physics/addProcess RadioactiveDecay
/gate/physics/processList Enabled
/gate/physics/processList Initialized
# I N I T I A L I Z E
/gate/run/initialize
#=====================================================
# DIGITIZER
#=====================================================
/gate/digitizer/Singles/insert opticaladder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1
#=====================================================
# SOURCE
#=====================================================
/gate/source/addSource PETsource
/gate/source/PETsource/gps/particle e+
/gate/source/PETsource/gps/energytype Fluor18
/gate/source/PETsource/setActivity 100 becquerel
/gate/source/PETsource/gps/angtype iso
/gate/source/PETsource/gps/type Volume
/gate/source/PETsource/gps/shape Cylinder
/gate/source/PETsource/gps/radius 0.3 mm
/gate/source/PETsource/gps/halfz 7 mm
/gate/source/PETsource/setForcedHalfLife 6586 s
/gate/source/PETsource/gps/direction -10. 0. 0.
/gate/source/PETsource/gps/centre 4. 0. 0. mm
/gate/source/PETsource/gps/angtype iso
/gate/source/PETsource/visualize 2000 red 10
/gate/source/list
#=====================================================
# D A T A O U T P U T
#=====================================================
/gate/output/root/enable
/gate/output/root/setFileName OutputTest
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootOpticalFlag 1
/gate/output/root/setRootRecordFlag 0
/gate/output/root/setSaveRndmFlag 0
/gate/output/root/setRootNtupleFlag 1
/gate/output/analysis/disable
/gate/output/fastanalysis/enable
#=====================================================
# M E A S U R E M E N T S E T T I N G S
#=====================================================
/gate/application/setTimeSlice 4 s
/gate/application/setTimeStart 0 s
/gate/application/setTimeStop 4 s
#/gate/output/allowNoOutput
/gate/application/startDAQ
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20141202/83816b9d/attachment-0001.html>
More information about the Gate-users
mailing list