17 Geant4 ExamplesDeBase

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 24

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 vis.mac 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. run1.mac
./exampleN03 run1.mac
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 run1.mac
quivaut
./exampleN03
et dans la fentre Session:
/control/execute run1.mac
Commande pour excuter les commandes contenues
dans le fichier macro run1.mac
16

Excuter votre code AVEC macros


Capture dcran des trois vnements
primaire +
gamma-ray

delta-ray
Qt

17

Changer la gomtrie
/control/execute newgeom.mac
/control/execute run1.mac
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

Vous aimerez peut-être aussi