magnetif fields examples

maureeungaro · maureeungaro · commit e926cb171109 · 2026-06-24T15:01:09.000-04:00
diff --git a/examples/basic/b1/annotations.yaml b/examples/basic/b1/annotations.yaml
@@ -1,3 +1,6 @@
+show_field_lines: 0
+show_auxiliary_edges: false
+
 g4decoration:
   scale: true
   axes: true
diff --git a/examples/basic/b1/b1.yaml b/examples/basic/b1/b1.yaml
@@ -1,6 +1,8 @@
 experiment: examples
 runno: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: b1
diff --git a/examples/basic/b2/b2.yaml b/examples/basic/b2/b2.yaml
@@ -1,6 +1,8 @@
 experiment: examples
 runno: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: b2
@@ -27,4 +29,3 @@ gmultipoles:
     strength: 2
     rotaxis: z
     rotation_angle: "30*deg"
-
diff --git a/examples/basic/b3/b3.yaml b/examples/basic/b3/b3.yaml
@@ -2,6 +2,8 @@ experiment: examples
 runno: 1
 n: 1
 nthreads: 4
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gparticle:
   - name: e-
diff --git a/examples/basic/material/annotations.yaml b/examples/basic/material/annotations.yaml
@@ -1,3 +1,6 @@
+show_field_lines: 0
+show_auxiliary_edges: false
+
 g4decoration:
   scale: true
   scaleLength: 10
diff --git a/examples/basic/material/materials.yaml b/examples/basic/material/materials.yaml
@@ -1,6 +1,8 @@
 experiment: examples
 runno: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: material
diff --git a/examples/basic/pyvista/pyvista.yaml b/examples/basic/pyvista/pyvista.yaml
@@ -1,6 +1,8 @@
 experiment: examples
 runno: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: pyvista
diff --git a/examples/basic/scintillator_barrel/scintillator_barrel.yaml b/examples/basic/scintillator_barrel/scintillator_barrel.yaml
@@ -2,6 +2,8 @@ experiment: examples
 runno: 1
 n: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: scintillator_barrel
@@ -25,4 +27,4 @@ g4camera:
 
 gstreamer:
   - format: csv
-    filename: barrel
+    filename: barrel
diff --git a/examples/basic/simple_flux/simple_flux.yaml b/examples/basic/simple_flux/simple_flux.yaml
@@ -1,6 +1,8 @@
 experiment: examples
 runno: 1
 nthreads: 1
+show_field_lines: 0
+show_auxiliary_edges: false
 
 gsystem:
   - name: simple_flux
@@ -22,4 +24,4 @@ gparticle:
     vz: -3*cm
     delta_vx: 0.1*cm
     delta_vy: 0.1*cm
-    multiplicity: 10
+    multiplicity: 10
diff --git a/examples/fields/constant/constant.py b/examples/fields/constant/constant.py
@@ -0,0 +1,22 @@
+#!/usr/bin/env python3
+from pygemc import autogeometry
+from pygemc import GVolume
+
+cfg = autogeometry("examples", "constant")
+
+world = GVolume("root")
+world.description = "World for constant field display"
+world.make_box(3.0, 3.0, 3.0, "m")
+world.material = "G4_AIR"
+world.visible = 0
+world.publish(cfg)
+
+field_volume = GVolume("constant_field_box")
+field_volume.description = "Box enclosing the uniform constant field"
+field_volume.make_box(1.0, 1.0, 1.0, "m")
+field_volume.material = "G4_AIR"
+field_volume.mfield = "constant"
+field_volume.color = "orchid"
+field_volume.style = 0
+field_volume.opacity = 0.35
+field_volume.publish(cfg)
diff --git a/examples/fields/constant/constant.yaml b/examples/fields/constant/constant.yaml
@@ -0,0 +1,21 @@
+experiment: examples
+runno: 1
+n: 0
+nthreads: 1
+show_field_lines: 30
+show_auxiliary_edges: true
+
+gsystem:
+  - name: constant
+    factory: sqlite
+
+gmultipoles:
+  - name: constant
+    pole_number: 2
+    strength: 1
+    rotaxis: z
+    longitudinal: true
+
+g4camera:
+  theta: 70*deg
+  phi: 45*deg
diff --git a/examples/fields/dipole/dipole.py b/examples/fields/dipole/dipole.py
@@ -0,0 +1,22 @@
+#!/usr/bin/env python3
+from pygemc import autogeometry
+from pygemc import GVolume
+
+cfg = autogeometry("examples", "dipole")
+
+world = GVolume("root")
+world.description = "World for dipole field display"
+world.make_box(1.8, 1.2, 1.2, "m")
+world.material = "G4_AIR"
+world.visible = 0
+world.publish(cfg)
+
+field_volume = GVolume("dipole_field_box")
+field_volume.description = "Box enclosing the ideal dipole field region"
+field_volume.make_box(75, 35, 35, "cm")
+field_volume.material = "G4_AIR"
+field_volume.mfield = "dipole"
+field_volume.color = "cornflowerblue"
+field_volume.style = 0
+field_volume.opacity = 0.35
+field_volume.publish(cfg)
diff --git a/examples/fields/dipole/dipole.yaml b/examples/fields/dipole/dipole.yaml
@@ -0,0 +1,21 @@
+experiment: examples
+runno: 1
+n: 0
+nthreads: 1
+show_field_lines: 80
+show_auxiliary_edges: true
+
+gsystem:
+  - name: dipole
+    factory: sqlite
+
+gmultipoles:
+  - name: dipole
+    pole_number: 2
+    strength: 2
+    rotaxis: z
+    rotation_angle: "30*deg"
+
+g4camera:
+  theta: 75*deg
+  phi: 35*deg
diff --git a/examples/fields/solenoid/solenoid.py b/examples/fields/solenoid/solenoid.py
@@ -0,0 +1,22 @@
+#!/usr/bin/env python3
+from pygemc import autogeometry
+from pygemc import GVolume
+
+cfg = autogeometry("examples", "solenoid")
+
+world = GVolume("root")
+world.description = "World for solenoid field display"
+world.make_box(7.2, 7.2, 7.2, "m")
+world.material = "G4_AIR"
+world.visible = 0
+world.publish(cfg)
+
+field_volume = GVolume("solenoid_field_tube")
+field_volume.description = "Cylindrical volume enclosing the solenoid ASCII field map"
+field_volume.make_tube(0, 3.0, 3.0, 0, 360, "m", "deg")
+field_volume.material = "G4_AIR"
+field_volume.mfield = "solenoid"
+field_volume.color = "seagreen"
+field_volume.style = 0
+field_volume.opacity = 0.3
+field_volume.publish(cfg)
diff --git a/examples/fields/solenoid/solenoid.yaml b/examples/fields/solenoid/solenoid.yaml
@@ -0,0 +1,30 @@
+experiment: examples
+runno: 1
+n: 0
+nthreads: 1
+show_field_lines: 20
+show_auxiliary_edges: true
+
+gsystem:
+  - name: solenoid
+    factory: sqlite
+
+gfields:
+  - name: solenoid
+    type: asciimap
+    symmetry: cylindrical-z
+    map: solenoid_map.txt
+    field_unit: T
+    interpolation: linear
+    coordinate1: "transverse,   4, 0*m, 3*m"
+    coordinate2: "longitudinal, 7, -3*m, 3*m"
+    vx: 0*cm
+    vy: 0*cm
+    vz: 0*cm
+    rx: 0*deg
+    ry: 0*deg
+    rz: 0*deg
+
+g4camera:
+  theta: 80*deg
+  phi: 25*deg
diff --git a/examples/fields/solenoid/solenoid_map.txt b/examples/fields/solenoid/solenoid_map.txt
@@ -0,0 +1,31 @@
+# Data-only cylindrical-z solenoid map for the asciimap plugin (see solenoid.yaml).
+# Complete 4 x 7 grid (made-up but smooth values for a runnable example).
+# Columns: transverse[m]  longitudinal[m]  Bt[T]  Bl[T]
+ 0.000   -3.000    0.000000    0.526996
+ 0.000   -2.000    0.000000    1.839397
+ 0.000   -1.000    0.000000    3.894004
+ 0.000    0.000   -0.000000    5.000000
+ 0.000    1.000   -0.000000    3.894004
+ 0.000    2.000   -0.000000    1.839397
+ 0.000    3.000   -0.000000    0.526996
+ 1.000   -3.000    0.066667    0.500646
+ 1.000   -2.000    0.044444    1.747427
+ 1.000   -1.000    0.022222    3.699304
+ 1.000    0.000   -0.000000    4.750000
+ 1.000    1.000   -0.022222    3.699304
+ 1.000    2.000   -0.044444    1.747427
+ 1.000    3.000   -0.066667    0.500646
+ 2.000   -3.000    0.133333    0.474297
+ 2.000   -2.000    0.088889    1.655457
+ 2.000   -1.000    0.044444    3.504604
+ 2.000    0.000   -0.000000    4.500000
+ 2.000    1.000   -0.044444    3.504604
+ 2.000    2.000   -0.088889    1.655457
+ 2.000    3.000   -0.133333    0.474297
+ 3.000   -3.000    0.200000    0.447947
+ 3.000   -2.000    0.133333    1.563488
+ 3.000   -1.000    0.066667    3.309903
+ 3.000    0.000   -0.000000    4.250000
+ 3.000    1.000   -0.066667    3.309903
+ 3.000    2.000   -0.133333    1.563488
+ 3.000    3.000   -0.200000    0.447947
diff --git a/examples/fields/torus/torus.py b/examples/fields/torus/torus.py
@@ -0,0 +1,24 @@
+#!/usr/bin/env python3
+from pygemc import autogeometry
+from pygemc import GVolume
+
+cfg = autogeometry("examples", "torus")
+
+world = GVolume("root")
+world.description = "World for torus field display"
+world.make_box(7.0, 7.0, 7.0, "m")
+world.material = "G4_AIR"
+world.visible = 0
+world.publish(cfg)
+
+field_volume = GVolume("torus_field_volume")
+field_volume.description = "Toroidal volume enclosing the torus ASCII field map"
+field_volume.solid = "G4Torus"
+field_volume.parameters = "0*cm, 260*cm, 275*cm, 0*deg, 360*deg"
+field_volume.material = "G4_AIR"
+field_volume.mfield = "torus"
+field_volume.set_position(0, 0, 350, "cm")
+field_volume.color = "gold"
+field_volume.style = 0
+field_volume.opacity = 0.35
+field_volume.publish(cfg)
diff --git a/examples/fields/torus/torus.yaml b/examples/fields/torus/torus.yaml
@@ -0,0 +1,27 @@
+experiment: examples
+runno: 1
+n: 0
+nthreads: 1
+show_field_lines: 20
+show_auxiliary_edges: true
+
+gsystem:
+  - name: torus
+    factory: sqlite
+
+gfields:
+  - name: torus
+    type: asciimap
+    symmetry: phi-segmented
+    map: torus_map.txt
+    field_unit: kilogauss
+    interpolation: none
+    integration_stepper: G4ClassicalRK4
+    minimum_step: 1*mm
+    coordinate1: "azimuthal,    4, 0*deg,  30*deg"
+    coordinate2: "transverse,   6, 0*cm,   500*cm"
+    coordinate3: "longitudinal, 6, 100*cm, 600*cm"
+
+g4camera:
+  theta: 65*deg
+  phi: 35*deg
diff --git a/examples/fields/torus/torus_map.txt b/examples/fields/torus/torus_map.txt
