00001 /*! 00002 00003 $Id$ 00004 00005 Last update to: 00006 Id: Readme-Questions 172 2008-09-22 10:23:14Z maire 00007 00008 \file "jour4d/.Readme-Questions" 00009 \brief L'exercise de jour4d 00010 00011 \mainpage jour4d 00012 00013 \section s1 GEOMETRIE 00014 00015 Une navette spatiale est modelisee par un cylindre d'Aluminium. 00016 A l'interieur de la cabine se trouve un plancher, egalement en Aluminium. 00017 La navette est placée dans l'espace ('World') : une boite aussi petite que 00018 possible. 00019 00020 Un (mini)astronaute est placé sur le plancher. 00021 00022 00023 \section s2 COMMENT DEMARRER ? 00024 00025 - lire le code : jour4d.cc (main) et les 6 classes DetectorConstruction, 00026 PhysicsList, PrimaryGeneratorAction, RunAction, EventAction, SteppingAction 00027 00028 - compiler et generer un executable : 00029 \verbatim 00030 % gmake 00031 \endverbatim 00032 00033 - executer jour4d en mode interactif avec visualisation : 00034 \verbatim 00035 % $G4WORKDIR/bin/$G4SYSTEM/jour4d 00036 00037 Idle> /tutorial/gun/rndm off 00038 .... 00039 Idle> /run/beamOn 1 00040 .... 00041 Idle> exit 00042 \endverbatim 00043 00044 - executer jour4d en mode batch avec une macro : 00045 \verbatim 00046 % $G4WORKDIR/bin/$G4SYSTEM/jour4d run1.mac 00047 \endverbatim 00048 00049 00050 \section s3 EXERCICES 00051 00052 (Re-)Lire la classe PhysicsList, en particulier les fonctions 00053 ConstructParticles() et ConstructEM() 00054 00055 Y-a-t-il des processus pour les interactions hadroniques des protons avec 00056 les noyaux des atomes ? 00057 00058 \subsection sub1 Exercice 1 00059 00060 On se propose d'ecrire dans la classe PhysicsList une fonction ConstructHad() 00061 On s'inspirera du cours et du fichier HadronPhysics de ce directoire (jour4d) 00062 00063 a) Ajouter un processus pour la diffusion elastique des protons. 00064 Utiliser le modele G4LElastic. 00065 Compiler et executer en interactif 00066 00067 b) Ajouter un processus inelastique pour les protons. Utiliser une cascade 00068 (ex. G4BinaryCascade) comme modele. 00069 Compiler et executer en interactif 00070 00071 c) Utiliser une meilleure description des sections efficaces en utilisant la 00072 classe G4NucleonNuclearCrossSection. 00073 Compiler et executer en interactif 00074 00075 Note-1 : il faut aussi completer la liste des particules par les hadrons et 00076 ions. 00077 00078 Note-2 : Le cascade s'applique aux protons au dessous de 10 GeV maximum. 00079 Pour aller au dela, il faut ajouter un autre modele, par exemple FTF ou QGS. 00080 00081 d) Comme b) pour les neutrons. En plus, pour les neutrons, il faut ajouter 00082 - un processus de "capture at rest" : utiliser le modele G4LCapture 00083 - un processus de fission: utiliser le modele G4LFission. 00084 00085 Compiler et executer en interactif 00086 00087 \subsection sub2 Exercice 2 proton incident 00088 00089 a) visualiser qqe dizaines d'evenments; l'energie des protons doit etre 00090 celle utilisee dans jour4c, exercice 1b. 00091 Comparer les evenements avec ceux obtenus a jour4c/1b; commentaires ? 00092 La dose absorbee par l'astronaute a-t-elle change? Pourquoi? 00093 00094 b) On etudie la dose absorbee en fonction de l'energie de proton. 00095 Quelle est cette dose a 500 MeV ? 1 GeV ? 2 GeV ? 00096 Note: on peut utiliser run1.mac 00097 00098 c) L'importance des interactions hadronique. 00099 desactiver les interactions hadronique: 00100 \verbatim 00101 /process/list 00102 /process/inactivate HadronElastic 00103 /process/inactivate ProtonInelastic 00104 /process/inactivate NeutronInelastic 00105 /process/inactivate nCapture 00106 /process/inactivate nFission 00107 \endverbatim 00108 et refaire l'exercise b). On peut utilier run2.mac 00109 00110 d) Comparer les resultats de b) et c). Explications? 00111 00112 Si le temps le permet, on peut aussi evaluer l'importance des interactions 00113 des neutrons, ou de certains processus en deactivant les processus des 00114 neutrons, ou des processus specifiques. 00115 00116 \subsection sub3 Exercise 3 PhysicsLists pre-definies 00117 00118 a) Dans la PhysicsList creee, certains processus manquent, lesquels ? 00119 00120 Note-1: La cascade est limitee en energie. La solution (jour4e) inclut en 00121 commentaire une possibilite d'inclure des modeles a plus haute energie. 00122 00123 Note-2: Dans un interaction hadronique, des hadrons autre que proton 00124 et neutron peuvent etre crees. 00125 00126 b) Une alternative a la creation d'une PhysicsList par l'utilisateur est 00127 d'utiliser une PhysicsList distribuee avec Geant4: 00128 voir $G4INSTALL/geant4/source/physics_lists/lists 00129 00130 La PhysicsList de cet exercice ressemble a la liste de reference FTF_BIC, 00131 sauf que cette derniere couvre toutes les particules jusqu'aux hautes energies. 00132 00133 Modifier jour4d.cc pour utiliser FTF_BIC ( ou un autre ). 00134 (solution: voir jour4e-FTF_BIC.cc) 00135 00136 Compiler et executer avec cette PhysicsList. Verifier la liste des processus. 00137 00138 */ 00139 00140 00141
1.5.7.1