00001 /*! 00002 00003 $Id$ 00004 00005 Last update to: 00006 Id: Readme-Questions 155 2008-09-10 11:08:48Z maire 00007 00008 00009 \file "jour3e/.Readme-Questions" 00010 \brief L'exercise de jour3e 00011 00012 \mainpage jour3e 00013 00014 \section s1 GEOMETRIE 00015 00016 Une navette spatiale est modelisee par un cylindre d'Aluminium. 00017 A l'interieur de la cabine se trouve un plancher, egalement en Aluminium. 00018 La navette est placée dans l'espace ('World') : une boite aussi petite que 00019 possible. 00020 00021 Un (mini)astronaute est placé sur le plancher. 00022 00023 00024 \section s2 COMMENT DEMARRER ? 00025 00026 - lire le code : jour3e.cc (main) et les 6 classes DetectorConstruction, 00027 PhysicsList, PrimaryGeneratorAction, RunAction, EventAction, SteppingAction 00028 00029 - compiler et generer un executable : 00030 \verbatim 00031 % gmake 00032 \endverbatim 00033 00034 - executer jour3e en mode interactif avec visualisation : 00035 \verbatim 00036 % $G4WORKDIR/bin/$G4SYSTEM/jour3e 00037 .... 00038 Idle> /run/beamOn 1 00039 .... 00040 Idle> exit 00041 \endverbatim 00042 00043 - executer jour3e en mode batch avec une macro : 00044 \verbatim 00045 % $G4WORKDIR/bin/$G4SYSTEM/jour3e run1.mac 00046 \endverbatim 00047 00048 00049 \section s3 EXERCICES 00050 00051 Il s'agit de creer et de remplir des histogrammes via l'interface AIDA. 00052 Il faut : 00053 - tache 1 : creer ou ouvrir un fichier d'histogrammes 00054 - tache 2 : definir les histogrammes 00055 - tache 3 : remplir ces histogrammes avec les variables que l'on veut etudier 00056 - tache 4 : fermer le fichier et l'ecrire sur disque 00057 00058 Les taches 1-2 sont a effectuer une seule fois, en debut d'execution du job 00059 ou en debut de run. 00060 La tache 4 est aussi a effectuer une seule fois, en fin de run ou en fin 00061 d'execution du job. 00062 La tache 3 se fait autant de fois que necessaire en cours de run, evenement, 00063 trace, step ... 00064 00065 \subsection sub1 Exercice A 00066 00067 00068 On veut histogrammer l'energie totale recue par l'astronaute, par particule 00069 incidente. Le code est deja implementé dans jour3e et sera expliqué en 00070 sceance. Remarquer le flag de compilation G4ANALYSIS_USE. 00071 00072 - taches 1-2 : voir RunAction::BookHistoFile(). 00073 Noter: 00074 - choix du format d'ecriture (root ou hbook) 00075 - definition de l'histogramme en histo[0] 00076 00077 - tache 4 : voir RunAction::SaveHistoFile(). 00078 00079 \subsubsection subsub1 Question 1 : 00080 - quand ces fonctions sont-elles appelées ? 00081 - tache 3 : l'energie recue est calculée dans EventAction::EndOfEventAction(). 00082 Voir comment on remplit un histogramme avec cette information. 00083 Quelle est la signification des 2 arguments de la fonction fill() ? 00084 00085 \subsubsection subsub2 Question 2 : 00086 - activer le flag G4ANALYSIS_USE dans GNUmakefile 00087 - compiler 00088 - executer en mode batch avec la macro run1.mac 00089 le programme produit un fichier histograms.root (ou .hbook) 00090 - lancer ROOT (ou PAW) et executer le script plotHisto.C (ou .kumac) 00091 00092 \subsection sub2 Exercice B 00093 00094 00095 On veut histogrammer le profile longitudinal de l'energie totale recue par 00096 l'astronaute au cours d'un run : Erun = f(zlocal) 00097 (information similaire a celle contenue dans le tableau depthDose.ascii) 00098 00099 - quel type d'histogramme faut-il definir ? 00100 - ajouter ce second histogramme dans BookHistoFile(), en histo[1] 00101 - remplir l'histogramme dans SteppingAction::UserSteppingAction() 00102 - compiler et executer jour3e 00103 - completer le script plotHisto.C pour analyser le fichier histograms.root 00104 00105 \section s4 Corrige : voir code jour3f 00106 00107 */ 00108