[Gate-users] [Gate-devel] problem with parametrised collimator setup ...
Alex Vergara Gil
alexvergaragil at gmail.com
Tue Jul 7 15:09:28 CEST 2015
Dear Uwe
Can you send me the full macro so I can check what is happening here.
I see no visible errors in the macros you have shown, and the error
might be in the code.
Regards
Alex
2015-07-07 8:44 GMT-04:00, Pietrzyk, Uwe <u.pietrzyk at fz-juelich.de>:
>
> 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
> ------------------------------------------------------------------------------------------------
> ------------------------------------------------------------------------------------------------
>
>
More information about the Gate-users
mailing list