[Gate-users] Jaszczak phantom, hot & cold phantom example

라이노 wnsghd4491 at gmail.com
Wed Jun 3 07:46:36 CEST 2020


Now, i am simulating two phantom, Jaszczak & hot and cold phantom.
In the case of  Jaszczak phantom, it was  made as the code below.

/gate/world/daughters/name jaszczak
/gate/world/daughters/insert cylinder
/gate/jaszczak/geometry/setRmax 30 mm
/gate/jaszczak/geometry/setHeight 52.96 mm
/gate/jaszczak/setMaterial Air
/gate/jaszczak/vis/forceWireframe

/vis/viewer/set/viewpointThetaPhi 0 0 deg

/gate/jaszczak/daughters/name rods1
/gate/jaszczak/daughters/insert cylinder
/gate/rods1/geometry/setRmax 0.5 mm
/gate/rods1/geometry/setHeight 52.96 mm
/gate/rods1/setMaterial Air
/gate/rods1/placement/setTranslation 7.5 0. 0. mm
/gate/rods1/placement/setPhiOfTranslation 30. deg
/gate/rods1/repeaters/insert quadrant
/gate/rods1/quadrant/setCopySpacing 3. mm
/gate/rods1/quadrant/setOrientation 30. deg
/gate/rods1/quadrant/setLineNumber 8
/gate/rods1/quadrant/setMaxRange 26.5 mm

/gate/jaszczak/daughters/name rods2
/gate/jaszczak/daughters/insert cylinder
/gate/rods2/geometry/setRmax 0.55 mm
/gate/rods2/geometry/setHeight 52.96 mm
/gate/rods2/setMaterial Air
/gate/rods2/placement/setTranslation 7.5 0. 0. mm
/gate/rods2/placement/setPhiOfTranslation 90. deg
/gate/rods2/repeaters/insert quadrant
/gate/rods2/quadrant/setCopySpacing 4. mm
/gate/rods2/quadrant/setOrientation 90. deg
/gate/rods2/quadrant/setLineNumber 6
/gate/rods2/quadrant/setMaxRange 26.5 mm

/gate/jaszczak/daughters/name rods3
/gate/jaszczak/daughters/insert cylinder
/gate/rods3/geometry/setRmax 0.6 mm
/gate/rods3/geometry/setHeight 52.96 mm
/gate/rods3/setMaterial Air
/gate/rods3/placement/setTranslation 7.5 0. 0. mm
/gate/rods3/placement/setPhiOfTranslation 150. deg
/gate/rods3/repeaters/insert quadrant
/gate/rods3/quadrant/setCopySpacing 5. mm
/gate/rods3/quadrant/setOrientation 150. deg
/gate/rods3/quadrant/setLineNumber 5
/gate/rods3/quadrant/setMaxRange 26.5 mm

/gate/jaszczak/daughters/name rods4
/gate/jaszczak/daughters/insert cylinder
/gate/rods4/geometry/setRmax 0.65 mm
/gate/rods4/geometry/setHeight 52.96 mm
/gate/rods4/setMaterial Air
/gate/rods4/placement/setTranslation 7.5 0. 0. mm
/gate/rods4/placement/setPhiOfTranslation 210. deg
/gate/rods4/repeaters/insert quadrant
/gate/rods4/quadrant/setCopySpacing 6. mm
/gate/rods4/quadrant/setOrientation 210. deg
/gate/rods4/quadrant/setLineNumber 4
/gate/rods4/quadrant/setMaxRange 26.5 mm

/gate/jaszczak/daughters/name rods5
/gate/jaszczak/daughters/insert cylinder
/gate/rods5/geometry/setRmax 0.7 mm
/gate/rods5/geometry/setHeight 52.96 mm
/gate/rods5/setMaterial Air
/gate/rods5/placement/setTranslation 7.5 0. 0. mm
/gate/rods5/placement/setPhiOfTranslation 270. deg
/gate/rods5/repeaters/insert quadrant
/gate/rods5/quadrant/setCopySpacing 7. mm
/gate/rods5/quadrant/setOrientation 270. deg
/gate/rods5/quadrant/setLineNumber 4
/gate/rods5/quadrant/setMaxRange 26.5 mm

/gate/jaszczak/daughters/name rods6
/gate/jaszczak/daughters/insert cylinder
/gate/rods6/geometry/setRmax 0.75 mm
/gate/rods6/geometry/setHeight 52.96 mm
/gate/rods6/setMaterial Air
/gate/rods6/placement/setTranslation 7.5 0. 0. mm
/gate/rods6/placement/setPhiOfTranslation 330. deg
/gate/rods6/repeaters/insert quadrant
/gate/rods6/quadrant/setCopySpacing 8. mm
/gate/rods6/quadrant/setOrientation 330. deg
/gate/rods6/quadrant/setLineNumber 3
/gate/rods6/quadrant/setMaxRange 23. mm

# 2.  A T T A C H    C R Y S T A L  SD
/gate/jaszczak/attachPhantomSD
/gate/rods1/attachPhantomSD
/gate/rods2/attachPhantomSD
/gate/rods3/attachPhantomSD
/gate/rods4/attachPhantomSD
/gate/rods5/attachPhantomSD
/gate/rods6/attachPhantomSD

#================================================
# CUTS
#================================================

/gate/physics/Gamma/SetCutInRegion   topcrystal 1.0 cm
/gate/physics/Electron/SetCutInRegion  topcrystal 1.0 cm
/gate/physics/Positron/SetCutInRegion  topcrystal 1.0 cm

/gate/physics/Gamma/SetCutInRegion   bottomcrystal 1.0 cm
/gate/physics/Electron/SetCutInRegion  bottomcrystal 1.0 cm
/gate/physics/Positron/SetCutInRegion  bottomcrystal 1.0 cm

/gate/physics/Gamma/SetCutInRegion    jaszczak 0.1 mm
/gate/physics/Electron/SetCutInRegion   jaszczak 0.1 mm
/gate/physics/Positron/SetCutInRegion   jaszczak 0.1 mm

/gate/physics/SetMaxStepSizeInRegion    jaszczak 0.01 mm

#================================================
# INITIALISATION
#================================================

/gate/run/initialize

# ============================================
#  D E F I N I T I O N   O F   Y O U R   A C Q U I S I T I O N    #
# ============================================

#   A D D E R
/gate/digitizer/Singles/insert  adder

#       R E A D O U T
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1

#      E N E R G Y   B L U R R I N G
/gate/digitizer/Singles/insert  blurring
/gate/digitizer/Singles/blurring/setResolution 0.146
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV


#       E N E R G Y   C U T
/gate/digitizer/Singles/insert  thresholder
/gate/digitizer/Singles/thresholder/setThreshold 350. keV
/gate/digitizer/Singles/insert  upholder
/gate/digitizer/Singles/upholder/setUphold 625. keV

#       D E A D    T I M E
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 300. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume rsector

#   C O I N C I    S O R T E R
/gate/digitizer/Coincidences/setWindow 3.42 ns
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/describe
/gate/digitizer/name delay
/gate/digitizer/insert coincidenceSorter
/gate/digitizer/delay/setWindow 3.42 ns
/gate/digitizer/delay/setOffset 100. ns
/gate/digitizer/delay/describe
/gate/digitizer/name finalCoinc
/gate/digitizer/insert coincidenceChain
/gate/digitizer/finalCoinc/addInputName delay
/gate/digitizer/finalCoinc/addInputName Coincidences
/gate/digitizer/finalCoinc/usePriority true
/gate/digitizer/finalCoinc/insert deadtime
/gate/digitizer/finalCoinc/deadtime/setDeadTime 500. ns
/gate/digitizer/finalCoinc/deadtime/setMode paralysable
/gate/digitizer/finalCoinc/describe

# jaszczak Phantom


/gate/source/addSource rods1
/gate/source/rods1/setActivity     37000 kBq
/gate/source/rods1/gps/particle    gamma
/gate/source/rods1/gps/energy      140 keV
/gate/source/rods1/gps/angtype     iso
/gate/source/rods1/gps/type        Volume
/gate/source/rods1/gps/shape       Cylinder
/gate/source/rods1/gps/radius      0.5 mm
/gate/source/rods1/gps/halfz         26.48 mm
/gate/source/rods1/gps/centre      0. 0. 0. mm
/gate/source/rods1/gps/confine     rods1_P
/gate/source/rods1/gps/mintheta     79. deg
/gate/source/rods1/gps/maxtheta    101. deg
/gate/source/rods1/gps/minphi        0  deg
/gate/source/rods1/gps/maxphi      360. deg
/gate/source/list

/gate/source/addSource             rods2
/gate/source/rods2/setActivity     37000 kBq
/gate/source/rods2/gps/particle    gamma
/gate/source/rods2/gps/energy      140 keV
/gate/source/rods2/gps/angtype     iso
/gate/source/rods2/gps/type        Volume
/gate/source/rods2/gps/shape       Cylinder
/gate/source/rods2/gps/radius      0.55 mm
/gate/source/rods2/gps/halfz           26.48 mm
/gate/source/rods2/gps/centre      0. 0. 0. mm
/gate/source/rods2/gps/confine     rods2_P
/gate/source/rods2/gps/mintheta     79. deg
/gate/source/rods2/gps/maxtheta    101. deg
/gate/source/rods2/gps/minphi        0  deg
/gate/source/rods2/gps/maxphi      360. deg
/gate/source/list

/gate/source/addSource             rods3
/gate/source/rods3/setActivity     37000 kBq
/gate/source/rods3/gps/particle    gamma
/gate/source/rods3/gps/energy      140 keV
/gate/source/rods3/gps/angtype     iso
/gate/source/rods3/gps/type        Volume
/gate/source/rods3/gps/shape       Cylinder
/gate/source/rods3/gps/radius      0.6 mm
/gate/source/rods3/gps/halfz           26.48 mm
/gate/source/rods3/gps/centre      0. 0. 0. mm
/gate/source/rods3/gps/confine     rods3_P
/gate/source/rods3/gps/mintheta     79. deg
/gate/source/rods3/gps/maxtheta    101. deg
/gate/source/rods3/gps/minphi        0  deg
/gate/source/rods3/gps/maxphi      360. deg
/gate/source/list

/gate/source/addSource             rods4
/gate/source/rods4/setActivity     37000 kBq
/gate/source/rods4/gps/particle    gamma
/gate/source/rods4/gps/energy      140 keV
/gate/source/rods4/gps/angtype     iso
/gate/source/rods4/gps/type        Volume
/gate/source/rods4/gps/shape       Cylinder
/gate/source/rods4/gps/radius      0.65 mm
/gate/source/rods4/gps/halfz          26.48 mm
/gate/source/rods4/gps/centre      0. 0. 0. mm
/gate/source/rods4/gps/confine     rods4_P
/gate/source/rods4/gps/mintheta     79. deg
/gate/source/rods4/gps/maxtheta    101. deg
/gate/source/rods4/gps/minphi        0  deg
/gate/source/rods4/gps/maxphi      360. deg
/gate/source/list

/gate/source/addSource             rods5
/gate/source/rods5/setActivity     37000 kBq
/gate/source/rods5/gps/particle    gamma
/gate/source/rods5/gps/energy      140 keV
/gate/source/rods5/gps/angtype     iso
/gate/source/rods5/gps/type        Volume
/gate/source/rods5/gps/shape       Cylinder
/gate/source/rods5/gps/radius      0.7 mm
/gate/source/rods5/gps/halfz           26.48 mm
/gate/source/rods5/gps/centre      0. 0. 0. mm
/gate/source/rods5/gps/confine     rods5_P
/gate/source/rods5/gps/mintheta     79. deg
/gate/source/rods5/gps/maxtheta    101. deg
 /gate/source/rods5/gps/minphi        0  deg
/gate/source/rods5/gps/maxphi      360. deg
/gate/source/list

/gate/source/addSource             rods6
/gate/source/rods6/setActivity     37000 kBq
/gate/source/rods6/gps/particle    gamma
/gate/source/rods6/gps/energy      140 keV
/gate/source/rods6/gps/angtype     iso
/gate/source/rods6/gps/type        Volume
/gate/source/rods6/gps/shape       Cylinder
/gate/source/rods6/gps/radius      0.75 mm
/gate/source/rods6/gps/halfz           26.48 mm
/gate/source/rods6/gps/centre      0. 0. 0. mm
/gate/source/rods6/gps/confine     rods6_P
/gate/source/rods6/gps/mintheta     79. deg
/gate/source/rods6/gps/maxtheta    101. deg
/gate/source/rods6/gps/minphi        0  deg
/gate/source/rods6/gps/maxphi      360. deg
/gate/source/list
/vis/viewer/zoom 2

but strangely, it is noe being stored to .dat file.

In the case of the second, hot&cold phantom, I couldn't find an example of
how .mac file was made.

Question 1. which part of the Jaszczak phantom code is problematic?
Question 2. how to make hot and cold phantom using a GATE .mac file? if you
have example, please let me know.

[image: image.png]
답장전달
<https://drive.google.com/u/0/settings/storage?hl=ko&utm_medium=web&utm_source=gmail&utm_campaign=manage_storage>
<https://www.google.com/intl/ko/policies/terms/>
<https://www.google.com/intl/ko/policies/privacy/>
<https://www.google.com/gmail/about/policy/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200603/61fbb27f/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 212 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200603/61fbb27f/attachment-0001.png>


More information about the Gate-users mailing list