Difference between revisions of "Simulations of Particle Interactions with Matter"
| Line 78: | Line 78: | ||
==Building GEANT4.10== | ==Building GEANT4.10== | ||
| + | |||
| + | ===4.11.2=== | ||
| + | [[TF_GEANT4.11]] | ||
===4.10.02=== | ===4.10.02=== | ||
Revision as of 22:30, 15 January 2024
Class Admin
Homework Problems
HomeWork_Simulations_of_Particle_Interactions_with_Matter
Introduction
Energy Loss
Ann. Phys. vol. 5, 325, (1930)
Interactions of Electrons and Photons with Matter
Hadronic Interactions
Final Project
A final project will be submitted that will be graded with the following metrics:
1.) The document must be less than 15 pages.
2.) The document must contain references in a bibliography (5 points) .
3.) A comparison must be made between GEANT4's prediction and either the prediction of someone else or an experimental result(30 points).
4.) The graphs must be of publication quality with font sizes similar or larger than the 12 point font (10 points).
5.) The document must be grammatically correct (5 points).
6.) The document format must contain the following sections: An abstract of 5 sentences (5 points) , an Introduction(10 points), a Theory section (20 points) , if applicable a section describing the experiment that was simulated, a section delineating the comparisons that were made, and a conclusion( 15 points).
Resources
Saving/restoring Random number seed
You save the current state of the random number generator with the command
/random/setSavingFlag 1
/run/beamOn 100
/random/saveThisRun
A file is created called
currentEvent.rndm
/control/shell mv currentEvent.rndm currentEvent10.rndm
You can restore the random number generator and begin generating random number from the last save time
/random/resetEngineFrom currentEvent.rndm
Building GEANT4.10
4.11.2
4.10.02
4.10.01
Building GEANT4.9.6
Building GEANT4.9.5
An old version of Installation notes for versions prior to 9.5
Visualization Libraries:
Compiling G4 with ROOT
These instruction describe how you can create a tree within ExN02SteppingVerbose to store tracking info in an array (max number of steps in a track is set to 100 for the desired particle)
Using SLURM
http://slurm.schedmd.com/quickstart.html
https://rc.fas.harvard.edu/resources/documentation/convenient-slurm-commands/
simple batch script for one process job
create the file submit.sbatch below
#!/bin/sh #SBATCH --time=1 cd src/PI ./PI_MC 100000000000000
the execute
- sbatch submit.sbatch
check if its running with
- squeue
to kill a batch job
- scancel JOBID
On minerve
Sample script to submit 10 batch jobs.
the filename is minervesubmit and you run like
source minervesubmit
cd /home/foretony/src/GEANT4/geant4.9.5/Simulations/N02wROOT/batch qsub submit10mil qsub submit20mil qsub submit30mil qsub submit40mil qsub submit50mil qsub submit60mil qsub submit70mil qsub submit80mil qsub submit90mil qsub submit100mil
The file submit10mil looks like this
#!/bin/sh #PBS -l nodes=1 #PBS -A FIAC #PBS -M foretony@isu.edu #PBS -m abe # source /home/foretony/src/GEANT4/geant4.9.5/geant4.9.6-install/bin/geant4.sh cd /home/foretony/src/GEANT4/geant4.9.5/Simulations/N02wROOT/batch/10mil ../../exampleN02 run1.mac > /dev/null
use
qstat
to check that the process is still running
use
qdel jobID#
if you want to kill the batch job, the jobID number shows up when you do stat.
for example
[foretony@minerve HW10]$ qstat Job id Name User Time Use S Queue ------------------------- ---------------- --------------- -------- - ----- 27033.minerve submit foretony 00:41:55 R default [foretony@minerve HW10]$ qdel 27033 [foretony@minerve HW10]$ qstat
Definitions of Materials
File:MCNP Compendium of Material Composition.pdf
Minerve2 GEANT 4.10.1 Xterm error
On OS X El Capitan V 10.11.4 using XQuartz
~/src/GEANT4/geant4.10.1/Simulations/B2/B2a/exsmpleB2a
# Use this open statement to create an OpenGL view: /vis/open OGL 600x600-0+0 /vis/sceneHandler/create OGL /vis/viewer/create ! ! 600x600-0+0 libGL error: failed to load driver: swrast X Error of failed request: BadValue (integer parameter out of range for operation) Major opcode of failed request: 150 (GLX) Minor opcode of failed request: 3 (X_GLXCreateContext) Value in failed request: 0x0 Serial number of failed request: 25 Current serial number in output stream: 26