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

Add t30 tests for isotopes #90

Merged
merged 1 commit into from
Feb 7, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension


Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion .github/workflows/main.yml
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Loading