[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