[Gate-users] problem with parametrised collimator setup ...
Pietrzyk, Uwe
u.pietrzyk at fz-juelich.de
Tue Jul 7 14:44:48 CEST 2015
Hello Gate-Users,
at least since the GATE version 6.2 there is the option of a parametrised collimator setup,
which allows to define the collimator fan beam and parallel beam geometries.
I wanted to compare both collimations to eventually include them as a new example
within the EduGATE-project.
However, I experience a strange problem:
fan beam:
GateVVolume.cc<http://GateVVolume.cc> (l.189): The material of the volume hole is not defined.
parallel beam:
GateVVolume.cc<http://GateVVolume.cc> (l.189): The material of the volume side hole 1 is not defined.
Obviously, the collimator is considered as a bloc of lead with hole “drilled” into and it seems,
that GATE does not know which material the holes a remade of. And I simply don’t know how
to specify the material of the holes.
Defining the collimator in the “classic” way was to specify lead for the collimator volume and
air for the holes being defined as daughters of the collimator volume.
There are also subtle differences in defining the material:
fan beam wants: /gate/fanbeam/setMaterial Lead
parallel beam wants: /gate/colli/setMaterialName Lead
Did anyone ever succeed to work with the with parametrised collimator setup for either type
or is there any example available??
Many thanks for any help,
Uwe Pietrzyk
——————————————————
[G4]
[G4] *************************************************************
[G4] Geant4 version Name: geant4-10-01 (5-December-2014)
[G4] Copyright : Geant4 Collaboration
[G4] Reference : NIM A 506 (2003), 250-303
[G4] WWW : http://cern.ch/geant4
[G4] *************************************************************
[G4]
[Core-0] Initialization of geometry
[Core-0] Initialization of physics
[Core-0] Initialization of actors
[Core-0]
[Core-0] **********************************************************************
[Core-0] GATE version name: gate_v7.1
[Core-0] Copyright : OpenGATE Collaboration
[Core-0] Reference : Phys. Med. Biol. 49 (2004) 4543-4561
[Core-0] Reference : Phys. Med. Biol. 56 (2011) 881-901
[Core-0] WWW : http://www.opengatecollaboration.org
[Core-0] **********************************************************************
[Core-0]
[Core-0] You are using Geant4 version 10.1.0
[Core-0] Starting macro Gamma_Camera.mac
GateDetectorConstruction::insertARFSD entered
created a ARF Sensivitive Detector
created command /gate/systems/SPECThead/arf/setARFStage
Constructor GateParameterisedCollimator - begin fan beam <——————————— start fan beam
m_InnerRadius = 0.5
ok
Nx 271
Ny 146
Dy 1.5
Dx 1.44338
OffsetX -194.856
OffsetY -109.25
Constructor GateParameterisedCollimator - end
ok
Nx 251
Ny 146
Dy 1.5
Dx 1.44338
OffsetX -180.422
OffsetY -109.25
ok
Nx 251
Ny 346
Dy 1.5
Dx 1.44338
OffsetX -180.422
OffsetY -259.25
ok
Nx 251
Ny 346
Dy 1.5
Dx 1.44338
OffsetX -180.422
OffsetY -259.25
ok
Nx 249
Ny 346
Dy 1.5
Dx 1.44338
OffsetX -178.979
OffsetY -259.25
ok
Nx 244
Ny 346
Dy 1.5
Dx 1.44338
OffsetX -175.37
OffsetY -259.25
ok
Nx 64
Ny 61
Dy 8.5
Dx 5.48483
OffsetX -172.772
OffsetY -259
ok
Nx 10
Ny 20
Dy 25
Dx 34.0637
OffsetX -153.286
OffsetY -241.5
GATE object: 'systems/SPECThead'
Components:
GATE object: 'systems/SPECThead/base'
Attached to volume: SPECThead
Nb of children: 1
GATE object: 'systems/SPECThead/crystal'
Attached to volume: crystal
Nb of children: 1
GATE object: 'systems/SPECThead/pixel'
Attached to volume: ---
[Physic-0] The following Geant4's physic-list is enabled :emstandard
[Core-0] Initialization of geometry
GateVVolume.cc<http://GateVVolume.cc> (l.189): The material of the volume hole is not defined. <—————— ????????
Nb of children: 0
root [0]
Processing Gamma_Camera.C...
root [1] medmac23 39 -->
This is the relevant macro: -------------------------
#
# C A M E R A
#
# Type: SPECThead
/gate/world/daughters/name SPECThead
/gate/world/daughters/insert box
/gate/SPECThead/geometry/setXLength 15. cm
/gate/SPECThead/geometry/setYLength 40. cm
/gate/SPECThead/geometry/setZLength 54. cm
/gate/SPECThead/placement/setTranslation {x_placement} 0. 0. cm
/gate/SPECThead/setMaterial Air
/gate/SPECThead/vis/setColor white
/gate/SPECThead/vis/forceWireframe
# -------------
# FanBeam COLLIMATOR
# Create the collimator volume
/gate/SPECThead/daughters/name fanbeam
/gate/SPECThead/daughters/insert collimator
#set the material for the collimator
/gate/fanbeam/setMaterial Lead
#/gate/fanbeam/setMaterial Air
#/gate/fanbeam/setMaterialName Lead
#/gate/fanbeam/setMaterialName Air
#define the X and Y size of the collimator
/gate/fanbeam/geometry/setDimensionX 38 cm
/gate/fanbeam/geometry/setDimensionY 52 cm
#specify the focal length
/gate/fanbeam/geometry/setFocalDistanceY 0.0 cm
/gate/fanbeam/geometry/setFocalDistanceX 35.0 cm
#specify the thickness of the collimator
/gate/fanbeam/geometry/setHeight 5.8 cm
#set the septal thickness to the required distance between the holes
/gate/fanbeam/geometry/setSeptalThickness 0.8 cm
#specify the hole radius
/gate/fanbeam/geometry/setInnerRadius 1.70 cm
/gate/fanbeam/placement/setRotationAxis 0 0 1
/gate/fanbeam/placement/setRotationAngle -90 deg
/gate/fanbeam/vis/setVisible 1
/gate/fanbeam/vis/setColor red
#/gate/fanbeam/vis/forceWireframe
/gate/fanbeam/vis/forceSolid
/gate/fanbeam/attachPhantomSD
# -------------
# CRYSTAL
# Create the crystal volume
/gate/SPECThead/daughters/name crystal
/gate/SPECThead/daughters/insert box
/gate/crystal/geometry/setXLength 1. cm
/gate/crystal/geometry/setYLength 19. cm
/gate/crystal/geometry/setZLength 28. cm
/gate/crystal/placement/setTranslation 0. 0. 0. cm
/gate/crystal/setMaterial {CrystalMaterial}
/gate/crystal/vis/setColor yellow
#/gate/crystal/vis/forceWireframe
/gate/crystal/vis/forceSolid
/gate/systems/SPECThead/crystal/attach crystal
/gate/crystal/attachCrystalSD
/gate/SPECThead/attachPhantomSD
/gate/systems/SPECThead/describe
---------------------#### the following part is for parallel hole collimator!!!!!!!! <----------------------------------------
[G4]
[G4] *************************************************************
[G4] Geant4 version Name: geant4-10-01 (5-December-2014)
[G4] Copyright : Geant4 Collaboration
[G4] Reference : NIM A 506 (2003), 250-303
[G4] WWW : http://cern.ch/geant4
[G4] *************************************************************
[G4]
[Core-0] Initialization of geometry
[Core-0] Initialization of physics
[Core-0] Initialization of actors
[Core-0]
[Core-0] **********************************************************************
[Core-0] GATE version name: gate_v7.1
[Core-0] Copyright : OpenGATE Collaboration
[Core-0] Reference : Phys. Med. Biol. 49 (2004) 4543-4561
[Core-0] Reference : Phys. Med. Biol. 56 (2011) 881-901
[Core-0] WWW : http://www.opengatecollaboration.org
[Core-0] **********************************************************************
[Core-0]
[Core-0] You are using Geant4 version 10.1.0
[Core-0] Starting macro Gamma_Camera.mac
/tracking/storeTrajectory 1
Attributes available for modeling and filtering with
"/vis/modeling/trajectories/create/drawByAttribute" and
"/vis/filtering/trajectories/create/attributeFilter" commands:
G4TrajectoriesModel:
Event ID (EventID): G4int
Run ID (RunID): G4int
G4Trajectory:
Charge (Ch): unit: e+ (G4double)
Track ID (ID): G4int
Initial kinetic energy (IKE): G4BestUnit (G4double)
Initial momentum magnitude (IMag): G4BestUnit (G4double)
Initial momentum (IMom): G4BestUnit (G4ThreeVector)
No. of points (NTP): G4int
PDG Encoding (PDG): G4int
Parent ID (PID): G4int
Particle Name (PN): G4String
G4TrajectoryPoint:
Position (Pos): G4BestUnit (G4ThreeVector)
WARNING: Trajectory storing has been requested. This action may be
reversed with "/tracking/storeTrajectory 0".
WARNING: The vis manager will keep up to 20 events.
This may use a lot of memory.
It may be changed with, e.g., "/vis/scene/endOfEventAction accumulate 10".
GateDetectorConstruction::insertARFSD entered
created a ARF Sensivitive Detector
created command /gate/systems/SPECThead/arf/setARFStage
GetDirectoryName()+ setHeight = /gate/central hole/geometry/setHeight
GATE object: 'systems/SPECThead'
Components:
GATE object: 'systems/SPECThead/base'
Attached to volume: SPECThead
Nb of children: 1
GATE object: 'systems/SPECThead/crystal'
Attached to volume: crystal
Nb of children: 1
GATE object: 'systems/SPECThead/pixel'
Attached to volume: ---
[Physic-0] The following Geant4's physic-list is enabled :emstandard
[Core-0] Initialization of geometry
GateVVolume.cc<http://GateVVolume.cc> (l.189): The material of the volume side hole 1 is not defined. <---------------------
Nb of children: 0
root [0]
Processing Gamma_Camera.C...
root [1] medmac23 40 -->
This is the relevant macro: -------------------------
#
# C A M E R A
#
# Type: SPECThead
/gate/world/daughters/name SPECThead
/gate/world/daughters/insert box
/gate/SPECThead/geometry/setXLength 15. cm
/gate/SPECThead/geometry/setYLength 40. cm
/gate/SPECThead/geometry/setZLength 54. cm
/gate/SPECThead/placement/setTranslation {x_placement} 0. 0. cm
/gate/SPECThead/setMaterial Air
/gate/SPECThead/vis/setColor white
/gate/SPECThead/vis/forceWireframe
#
# -------------
# ParallelBeam COLLIMATOR
# Create the collimator volume
/gate/SPECThead/daughters/name colli
#specify that the parallel beam collimator setup must be used
/gate/SPECThead/daughters/insert parallelbeam
#set the collimator material
/gate/colli/setMaterialName Lead
#set the collimator dimensions
/gate/colli/geometry/setDimensionX 38 cm
/gate/colli/geometry/setDimensionY 52 cm
#set the thickness of the collimator
/gate/colli/geometry/setHeight 3 cm
#specify the hole radius
/gate/colli/geometry/setInnerRadius 0.5 cm
#set the septal thickness to the required
#distance between the holes
/gate/colli/geometry/setSeptalThickness 0.2 cm
/gate/colli/placement/alignToX
/gate/colli/placement/setRotationAxis 0 0 1
/gate/colli/placement/setRotationAngle -90 deg
/gate/colli/placement/setTranslation -2. 0. 0. cm
/gate/colli/vis/setVisible 1
/gate/colli/vis/setColor red
#/gate/colli/vis/forceWireframe
/gate/colli/vis/forceSolid
/gate/colli/attachPhantomSD
# -------------
# CRYSTAL
# Create the crystal volume
/gate/SPECThead/daughters/name crystal
/gate/SPECThead/daughters/insert box
/gate/crystal/geometry/setXLength 1. cm
/gate/crystal/geometry/setYLength 19. cm
/gate/crystal/geometry/setZLength 28. cm
/gate/crystal/placement/setTranslation 0. 0. 0. cm
/gate/crystal/setMaterial {CrystalMaterial}
/gate/crystal/vis/setColor yellow
#/gate/crystal/vis/forceWireframe
/gate/crystal/vis/forceSolid
/gate/systems/SPECThead/crystal/attach crystal
/gate/crystal/attachCrystalSD
/gate/SPECThead/attachPhantomSD
/gate/systems/SPECThead/describe
/----------------------------------------------------------------------
Prof. Dr. Uwe Pietrzyk
Diplom-Physiker
Institut für Neurowissenschaften und Medizin / INM-4
Forschungszentrum Jülich GmbH
52425 Jülich
Tel: +49-2461-61 4763
Fax: +49-2461-61 2820
Mobil: +49-162-103-4871
eMail: U.Pietrzyk at fz-juelich.de<mailto:U.Pietrzyk at fz-juelich.de>
http://www.fz-juelich.de/inm/inm-4/DE/Home/home_node.html
&
Bergische Universität Wuppertal
Fachbereich C - Physik / D.08.01
Gaußstr. 20
42097 Wuppertal
Tel: +49-(0)202-439-3523
Fax: +49-(0)202-439-2811
eMail: Uwe.Pietrzyk at physik.uni-wuppertal.de<mailto:Uwe.Pietrzyk at physik.uni-wuppertal.de>
http://www.medizinphysik.uni-wuppertal.de<http://www.medizinphysik.uni-wuppertal.de/>
----------------------------------------------------------------------/
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich
Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher
Geschaeftsfuehrung: Prof. Dr.-Ing. Wolfgang Marquardt (Vorsitzender),
Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,
Prof. Dr. Sebastian M. Schmidt
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150707/406204c2/attachment-0001.html>
More information about the Gate-users
mailing list