Premiers exemples
25 minutes
Les exemples de Geant4
l
Des exemples sont dj disponibles dans Geant4
Situs dans $G4SRC
3 catgories
Novice : fonctionnalits de base
Extended : fonctionnalits spcifiques
Processus Physiques, biasing , champs lectromagntiques
Advanced : simulations compltes de cas dutilisation
Physique mdicale, spatiale, calorimtrie
Exemple novice N01
l
Gomtrie
Volume de gaz Ar contenant un cylindre
de Al et un bloc de Pb avec des
tranches de Al
Particules incidentes: geantino
pas dinteractions physiques
Seul le transport est actif
Gestion batch et verbosit
Essayons lexemple N01
l
Copier N01 vers votre home
cd
cp R $G4SRC/examples/novice/N01 .
ls
mkdir N01build
cd N01build
Lire le fichier README : exampleN01 est lexemple le plus simple parmi les exemples
novices. Il dcrit un tube de tracking et un calorimtre sandwich constitu de botes.
Il tire 3 geantinos (pseudo-particules qui nont pas dinteraction) par vnement
Permet de montrer le fonctionnement de base de Geant4
nedit ../N01/README &
Compiler & lier lexemple
cmake -DGeant4_DIR=$G4DIR ../N01
make j2
Une fois compil et li, lancer lexcutable avec
./exampleN01
Attention : il faut toujours que les variables denvironnement soient dfinies avant de
4
compiler et de lancer votre excutable
Exemple novice N03
l
Calorimtre chantillonnage avec des couches de Pb
(absorbeur) et des gaps de lAr (dtection, replicas)
Liste exhaustive de matriaux
Interface de commandes
Randomization du faisceau incident
Tous les processus EM + dcroissance avec coupures pour
, e+, e- ( utiliser pour ltude des gerbes)
Rponse du dtecteur: E dpose, longueur des traces
dans labsorbeur et les gaps
Visualisation
Gestion de la graine alatoire
Essayons lexemple N03
l
Copier N03 vers votre home
cd
cp R $G4SRC/examples/novice/N03 .
mkdir N03build
cd N03build
Lire le fichier README : simulation dun simple calorimtre
Compiler & lier lexample
cmake -DGeant4_DIR=$G4DIR ../N03
make j2
Excuter votre code SANS macros
l
Une fois compil et li, vous pouvez
excuter votre code avec
./exampleN03
Pas darguments aprs le nom
Utiliser la touche tabulation pour une
reconnaissance rapide des rpertoires
existants
l
Alors vous obtenez laffichage suivant
Excuter votre code SANS macros
Que sest-il pass ?
1.
Le run est initialis
1.
2.
3.
4.
5.
dfinition des matriaux,
construction de la gomtrie
dclaration des processus physiques
dclaration des coupures
...
2.
Un fichier macro [Link] est automatiquement lu pour activer les
pilotes de visualisation et le dispositif simul apparat dans une
fentre graphique (interface Qt)
3.
On peut saisir des commandes la main dans la fentre
Session
Par ex., changer la gomtrie, choisir quelles particules tirer,
quelle nergie, excuter un autre fichier macro, ...
9
Gomtrie par dfaut
Obtenu avec le
driver VRML
10 couches : 10 mm Pb + 5 mm L. Ar
pas de champ magntique
10
Excuter votre code sans macros
Essayons de tirer une particule
Dans la fentre Session: /run/beamOn 1
Par dfaut, un lectron de 50 MeV est tir en incidence normale
sur le calorimtre
phot:
for
gamma
SubType= 12
Initialisation des
tables de
physique
LambdaPrime table from 200 keV to 10 TeV in 54 bins
===== EM models for the G4Region
PhotoElectric :
Emin=
0 eV
DefaultRegionForTheWorld ======
Emax=
10 TeV
AngularGenSauterGavrila
[...]
Index : 1
used in the geometry : Yes
recalculation needed : No
Material : Lead
Range cuts
gamma
1 mm
e-
1 mm
Energy thresholds :
gamma
101.843 keV
e-
Region(s) which use this couple :
DefaultRegionForTheWorld
e+
1 mm
1.36749 MeV
Calcul des
coupures
proton 1 mm
e+
1.27862 MeV proton 100 keV
11
Excuter votre code sans macros
---> Begin of event: 0
---> End of event: 0
Absorber: total energy: 40.0567 MeV
total track length: 2.9093 cm
Gap: total energy: 5.52904 MeV
total track length: 2.63017 cm
Run terminated.
Run Summary
Number of events processed : 1
User=0.01s Real=0.12s Sys=0.02s
--------------------End of Run-----------------------------mean Energy in Absorber : 40.0567 MeV +- 0 eV
mean Energy in Gap
: 5.52904 MeV +- 0 eV
mean trackLength in Absorber : 2.9093 cm +- 0 fm
mean trackLength in Gap
: 2.63017 cm +- 0 fm
------------------------------------------------------------
Rsum de
levent
Rsum du
run
12
Excuter votre code sans macros
On visualise aussi lvnement tir
(lectron de 50 MeV)
gamma-rays
Code couleur
par dfaut
dans Geant4:
rouge = ngative
bleu = positive
vert = neutre
primaire ePilote DAWNFILE
13
Excuter votre code AVEC macros
Largument qui suit le nom de lexcutable peut tre un
fichier macro, par ex. [Link]
./exampleN03 [Link]
Les macros Geant4 sont des fichiers ASCII contentant
une squence de commandes Geant4:
#
/run/verbose 2
/event/verbose 0
/tracking/verbose 1
#
/gun/particle mu+
/gun/energy 300 MeV
/run/beamOn 3
Tire 3 + dnergie 300 MeV
14
Excuter votre code AVEC macros
---> End of event: 2
Absorber: total energy: 123.976 MeV
Gap: total energy: 10.7747 MeV
Run terminated.
Run Summary
Number of events processed : 3
User=0.01s Real=0.02s Sys=0.01s
total track length:
10.44 cm
total track length: 5.16631 cm
--------------------End of Run-----------------------------mean Energy in Absorber : 125.191 MeV +- 2.67303 MeV
mean Energy in Gap
: 9.94337 MeV +- 596.048 keV
Rsum du
deuxime
vnement
mean trackLength in Absorber : 10.5231 cm +- 1.14618 mm
mean trackLength in Gap
: 5.09232 cm +- 524.76 um
------------------------------------------------------------
Rsum du run
entier
15
Excuter votre code AVEC macros
Noter que
./exampleN03 [Link]
quivaut
./exampleN03
et dans la fentre Session:
/control/execute [Link]
Commande pour excuter les commandes contenues
dans le fichier macro [Link]
16
Excuter votre code AVEC macros
Capture dcran des trois vnements
primaire +
gamma-ray
delta-ray
Qt
17
Changer la gomtrie
/control/execute [Link]
/control/execute [Link]
1) La premire macro change la gomtrie
- seulement une tranche dabsorbeur (40 cm
deau), pas despace (thickness = 0 cm)
un bloc solide deau
- change les dimensions transverses,
/N03/det/setNbOfLayers 1
/N03/det/setAbsMat Water
/N03/det/setAbsThick 40 cm
/N03/det/setGapMat Air
/N03/det/setGapThick 0 cm
/N03/det/setSizeYZ
40 cm
/N03/det/setField 3 tesla
/N03/det/update
- ajoute un champ magntique de 3T
2) La seconde macro tire 3 + de 300 MeV, comme avant
18
Changer la gomtrie
y
primaire
Gometrie,
matriaux, champ
magntique et
particlues primaires
peuvent tre
slectionns dans
des macros ASCII,
sans recompiler le
code !
Pilote DAWNFILE
19
Autres exemples
20
Microfaisceau pour
lirradiation cellulaire
l
Simuler la ligne dirradiation cellulaire
par microfaisceau du CENBG, incluant
les aimants de focalisation
Rsolution faisceau lair ~1 m
Fantme cellulaire obtenu partir de
microscopie confocale et danalyse par
faisceau dions (PIXE, STIM) pour la
composition chimique
Dosimtrie lchelle sub-cellulaire
(noyau, cytoplasme)
Utilise les constructeurs de physique
CENBG/AIFIRA
HaCat
21
Dnaphysics
l
Explique comment utiliser
lextension Geant4-DNA
Ensemble de modles de
Physique jusquau sub-eV pour
la simulation des traces dans
leau liquide lchelle du
nanomtre
Dvelopp dans le cadre du
projet Geant4-DNA pour la
simulation des effets
biologiques des radiations
22
Nanobeam
l
OM50
Simulation de la ligne nanofaisceau
du CENBG
Rsolution du faisceau sous vide de
qlqs dizaines de nm
Nano-technologies
Physique EM standard
Description fine des quadruples
Champs analytiques avec/sans franges
Carte de champ interpole (OPERA3D)
Geant4 peut calculer
Images de spot faisceau
Images de grilles
Coefficients daberration intrinsques
23
En rsum
l Nombreux
exemples disponibles
dans Geant4 pour lapprentissage et
lillustration dutilisation dans divers
domaines
Novices, tendus, avancs
l Localiss
dans $G4SRC!
24