Skip to content

Commit

Permalink
Merge pull request #90 from OpenGATE/t30_isotope
Browse files Browse the repository at this point in the history 
Add t30 tests for isotopes
  • Loading branch information
@tbaudier
tbaudier authored Feb 7, 2024
2 parents 3d7dae6 + f284252 commit 87b7219
Show file tree
Hide file tree
Showing 153 changed files with 1,760 additions and 1 deletion.
3 changes: 2 additions & 1 deletion .github/workflows/main.yml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,8 @@ jobs:
t28_Singles_2SD,
t29_optical_digi,
t30_dna,
t31_vpgTLE-tt]
t31_vpgTLE-tt,
t32_isotopes]

steps:
- name: Checkout github repo
Expand Down
18 changes: 18 additions & 0 deletions t32_isotopes/data/GateMaterials.db
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,18 @@
[Isotopes]
C14: Z=6 ; N=14 ; A= 14.0032419884 g/mole

[Elements]
Hydrogen: S= H ; Z= 1. ; A= 1.01 g/mole
Carbon: S= C ; Z= 6. ; A= 12.01 g/mole
C14: n=1 ; S= C14
+iso: name=auto ; f=1

[Materials]
Vacuum: d=0.000001 mg/cm3 ; n=1
+el: name=Hydrogen ; n=1

C: d=2.1 g/cm3; n=1; state=solid
+el: name=Carbon; n=1

C14: d=2.1 g/cm3; n=1; state=solid
+el: name=auto; n=1
121 changes: 121 additions & 0 deletions t32_isotopes/mac/main.mac
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
#/control/alias mainMaterial "C"
#/control/alias npart "1000000"
#/control/alias suffix "1M"
#/control/alias physicList "QGSP_BIC_HP"

/control/execute mac/verbose.mac

#=====================================================
# GEOMETRY
#=====================================================

/gate/geometry/setMaterialDatabase data/GateMaterials.db

# World
/gate/world/geometry/setXLength 2 m
/gate/world/geometry/setYLength 2 m
/gate/world/geometry/setZLength 4 m
/gate/world/setMaterial Vacuum


# Sphere (virtual detector)
/gate/world/daughters/name shellcyl1
/gate/world/daughters/insert sphere
/gate/shellcyl1/geometry/setRmin 75 cm
/gate/shellcyl1/geometry/setRmax 75.1 cm
/gate/shellcyl1/placement/setTranslation 0.0 0.0 25 cm
/gate/shellcyl1/setMaterial Vacuum

# Global Box
/gate/world/daughters/name mainbox1
/gate/world/daughters/insert box
/gate/mainbox1/geometry/setXLength 20 cm
/gate/mainbox1/geometry/setYLength 20 cm
/gate/mainbox1/geometry/setZLength 100 cm
/gate/mainbox1/placement/setTranslation 0.0 0.0 25 cm
/gate/mainbox1/setMaterial {mainMaterial}
/gate/mainbox1/vis/setVisible 1
/gate/mainbox1/vis/setColor yellow


#=====================================================
# PHYSICS
#=====================================================

/gate/physics/addPhysicsList {physicList}

/gate/physics/Gamma/SetCutInRegion world 10 mm
/gate/physics/Electron/SetCutInRegion world 10 mm
/gate/physics/Positron/SetCutInRegion world 10 mm

/gate/physics/Proton/SetCutInRegion mainbox1 0.5 mm

/gate/physics/SetMaxStepSizeInRegion world 10 mm
/gate/physics/SetMaxStepSizeInRegion mainbox1 1 mm

/gate/physics/ActivateStepLimiter proton

/gate/physics/displayCuts


#=====================================================
# OUTPUT
#=====================================================

/gate/actor/addActor KillActor MyActor
/gate/actor/MyActor/save output/part-{mainMaterial}-{physicList}-{suffix}.txt
/gate/actor/MyActor/attachTo world
/gate/actor/MyActor/addFilter particleFilter
/gate/actor/MyActor/particleFilter/addParticle neutron

/gate/actor/addActor SimulationStatisticActor stat
/gate/actor/stat/save output/stat-{mainMaterial}-{physicList}-{suffix}.txt

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

/gate/run/initialize
# Enable the following lines to display available and enabled processes
# /gate/physics/processList Available
# /gate/physics/processList Enabled


#=====================================================
# BEAMS
#=====================================================

/gate/source/addSource mybeam gps

/gate/source/mybeam/gps/particle proton
/gate/source/mybeam/gps/pos/type Beam
/gate/source/mybeam/gps/pos/shape Cylinder
/gate/source/mybeam/gps/pos/centre 0 0 -75 cm
/gate/source/mybeam/gps/pos/radius 1 mm
/gate/source/mybeam/gps/pos/halfz 1 cm
/gate/source/mybeam/gps/ene/mono 100 MeV
/gate/source/mybeam/gps/ene/type Gauss
/gate/source/mybeam/gps/ene/sigma 2.0 MeV
/gate/source/mybeam/gps/direction 0 0 1

#=====================================================
# VISUALISATION
#=====================================================

#/control/execute mac/visu.mac

#=====================================================
# START BEAMS
#=====================================================

#MersenneTwister
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto

/material/g4/printMaterial {mainMaterial}

/gate/application/setTotalNumberOfPrimaries {npart}
/gate/application/start

exit

28 changes: 28 additions & 0 deletions t32_isotopes/mac/verbose.mac
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@

# G4 verbose
/control/verbose 0
/run/verbose 0
/event/verbose 0
/tracking/verbose 0

# gate verbose

/gate/application/verbose 0
/gate/generator/verbose 0
/gate/source/verbose 0
/gate/random/verbose 0
/gate/output/verbose 0
/gate/verbose Physic 0
/gate/verbose Cuts 1
/gate/verbose SD 0
/gate/verbose Actions 0
/gate/verbose Actor 0
/gate/verbose Step 0
/gate/verbose Error 0
/gate/verbose Warning 0
/gate/verbose Output 0
/gate/verbose Beam 0
/gate/verbose Volume 2
/gate/verbose Image 0
/gate/verbose Geometry 1
/gate/verbose Core 1
17 changes: 17 additions & 0 deletions t32_isotopes/mac/visu.mac
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,17 @@
/vis/open OGLIQt
/vis/drawVolume
/vis/viewer/flush
/tracking/storeTrajectory 1
/vis/scene/add/trajectories
/vis/scene/endOfEventAction accumulate

/vis/scene/add/axes 0 0 0 500 mm
/vis/scene/add/text 10 0 0 cm 20 0 0 X
/vis/scene/add/text 0 10 0 cm 20 0 0 Y
/vis/scene/add/text 0 0 10 cm 20 0 0 Z

#/vis/viewer/set/viewpointThetaPhi 90 -90
#/vis/viewer/panTo 0 0
#/vis/viewer/zoom 50

/vis/viewer/set/auxiliaryEdge true
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-FTFP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 40449
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-FTFP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 20131
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-FTFP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 201701
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 40472
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 20184
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 201175
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BIC_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 65827
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BIC_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 32702
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_BIC_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 329993
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_INCLXX_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 56790
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_INCLXX_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 28493
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-QGSP_INCLXX_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 285411
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-Shielding-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 40052
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-Shielding-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 20143
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-Shielding-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 202176
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-ShieldingM-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 40445
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-ShieldingM-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 20198
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C-ShieldingM-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 202023
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-FTFP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 192157
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-FTFP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 97032
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-FTFP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 964771
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 193587
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 96955
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 963899
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BIC_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 177851
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BIC_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 88970
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_BIC_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 890113
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_INCLXX_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 183609
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_INCLXX_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 91819
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-QGSP_INCLXX_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 920165
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-Shielding-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 191741
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-Shielding-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 96886
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-Shielding-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 960867
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-ShieldingM-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 192092
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-ShieldingM-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 96289
1 change: 1 addition & 0 deletions t32_isotopes/output-9.1/part-C14-ShieldingM-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
# NumberOfKillTracks = 961615
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-FTFP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 1000000
# NumberOfTracks = 2032192
# NumberOfSteps = 52263541
# NumberOfGeometricalSteps = 3794166
# NumberOfPhysicalSteps = 48469375
# ElapsedTime = 144.108
# ElapsedTimeWoInit = 142.591
# StartDate = Fri Apr 8 14:45:31 2022
# EndDate = Fri Apr 8 14:47:56 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 7013.09
# TPS (Track per sec) = 14251.9
# SPS (Step per sec) = 366529
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-FTFP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 500000
# NumberOfTracks = 1014288
# NumberOfSteps = 26126714
# NumberOfGeometricalSteps = 1894803
# NumberOfPhysicalSteps = 24231911
# ElapsedTime = 68.5382
# ElapsedTimeWoInit = 67.6891
# StartDate = Fri Apr 8 14:33:53 2022
# EndDate = Fri Apr 8 14:35:01 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 7386.71
# TPS (Track per sec) = 14984.5
# SPS (Step per sec) = 385981
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-FTFP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 5000000
# NumberOfTracks = 10151895
# NumberOfSteps = 261278146
# NumberOfGeometricalSteps = 18964827
# NumberOfPhysicalSteps = 242313319
# ElapsedTime = 722.52
# ElapsedTimeWoInit = 721.084
# StartDate = Fri Apr 8 14:08:13 2022
# EndDate = Fri Apr 8 14:20:15 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 6934
# TPS (Track per sec) = 14078.7
# SPS (Step per sec) = 362341
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-QGSP_BERT_HP-1M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 1000000
# NumberOfTracks = 2033401
# NumberOfSteps = 52266079
# NumberOfGeometricalSteps = 3797536
# NumberOfPhysicalSteps = 48468543
# ElapsedTime = 174.084
# ElapsedTimeWoInit = 173.206
# StartDate = Fri Apr 8 14:39:32 2022
# EndDate = Fri Apr 8 14:42:26 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 5773.47
# TPS (Track per sec) = 11739.8
# SPS (Step per sec) = 301757
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-QGSP_BERT_HP-500K.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 500000
# NumberOfTracks = 1013890
# NumberOfSteps = 26122767
# NumberOfGeometricalSteps = 1893699
# NumberOfPhysicalSteps = 24229068
# ElapsedTime = 88.3376
# ElapsedTimeWoInit = 87.1457
# StartDate = Fri Apr 8 14:30:50 2022
# EndDate = Fri Apr 8 14:32:18 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 5737.52
# TPS (Track per sec) = 11634.4
# SPS (Step per sec) = 299760
19 changes: 19 additions & 0 deletions t32_isotopes/output-9.1/stat-C-QGSP_BERT_HP-5M.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# NumberOfRun = 1
# NumberOfEvents = 5000000
# NumberOfTracks = 10149242
# NumberOfSteps = 261275107
# NumberOfGeometricalSteps = 18962252
# NumberOfPhysicalSteps = 242312855
# ElapsedTime = 752.39
# ElapsedTimeWoInit = 751.235
# StartDate = Fri Apr 8 13:40:04 2022
# EndDate = Fri Apr 8 13:52:36 2022
# StartSimulationTime = 0
# StopSimulationTime = 1
# CurrentSimulationTime = 1
# VirtualStartSimulationTime = 0
# VirtualStopSimulationTime = 1
# ElapsedSimulationTime = 1
# PPS (Primary per sec) = 6655.7
# TPS (Track per sec) = 13510.1
# SPS (Step per sec) = 347794
Loading

0 comments on commit 87b7219

Please sign in to comment.