Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

'Forbid' does not work with quadrant repeater #556

Open
jcaravaca opened this issue Dec 8, 2022 · 0 comments
Open

'Forbid' does not work with quadrant repeater #556

jcaravaca opened this issue Dec 8, 2022 · 0 comments

Comments

@jcaravaca
Copy link

jcaravaca commented Dec 8, 2022

Describe the bug
When trying to 'Forbid' particle generation in a quadrant repeater volume, only particles in the first volume are 'Forbidden'.
I am trying to model a cold-rod Derenzo phantom by only one of the rods per each set is 'cold' (see picture).
image

Desktop (please complete the following information):

  • OS: Ubuntu 18.4
  • Gate version or commit hash: v9.0
  • Geant4 version: 10.6
  • Root version: 6-16-00
  • compiler version: gcc 7.5.0

Minimal example

P H A N T O M

/gate/world/daughters/name IPhantom
/gate/world/daughters/insert cylinder
/gate/IPhantom/geometry/setRmax 108 mm # 216 mm
/gate/IPhantom/geometry/setRmin 0. mm
/gate/IPhantom/geometry/setHeight 88 mm
/gate/IPhantom/placement/setTranslation 0. 0. 0. mm
/gate/IPhantom/setMaterial Polyethylene
/gate/IPhantom/vis/setColor blue
/gate/IPhantom/vis/forceWireframe
/gate/IPhantom/attachPhantomSD

R O D part

/gate/IPhantom/daughters/name rods1
/gate/IPhantom/daughters/insert cylinder
/gate/rods1/geometry/setRmax 9.55 mm # 19.1 mm
/gate/rods1/geometry/setRmin 0. mm
/gate/rods1/geometry/setHeight 88 mm
/gate/rods1/placement/setTranslation 0. 29. 0. mm
/gate/rods1/placement/setPhiOfTranslation 150. deg
/gate/rods1/setMaterial Water
/gate/rods1/vis/setColor red
/gate/rods1/repeaters/insert quadrant
/gate/rods1/quadrant/setCopySpacing 31 mm
/gate/rods1/quadrant/setOrientation 150. deg
/gate/rods1/quadrant/setLineNumber 3
/gate/rods1/quadrant/setMaxRange 90. mm
/gate/rods1/attachPhantomSD

D E F I N E T H E S O U R C E

/gate/source/addSource background
/gate/source/background/gps/type Volume
/gate/source/background/gps/shape Cylinder
/gate/source/background/gps/radius 108 mm
/gate/source/background/gps/halfz 44. mm
/gate/source/background/gps/centre 0. 0. 0. mm
/gate/source/background/gps/particle gamma
/gate/source/background/gps/energy 140. keV
/gate/source/background/setActivity 350000 Bq #350 000 000
/gate/source/background/gps/angtype iso
/gate/source/background/gps/Forbid rods1

Expected behavior
According to the documentation, 'Forbid' should prevent from primary particles being generated in the chosen volumes, but in practice, it only applies to the first volume in the repeater (rods1 in example). The result should be something like this (achieved by including each rod manually without using the repeater):
image

Additional context
The similar command 'confine' works as expected. This is implemented by Geant4 while it seems that 'Forbid' is exclusive from GATE.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant