.Readme-Questions

Go to the documentation of this file.

/*!

$Id$

Last update to:
Id: Readme-Questions 122 2008-06-30 15:25:12Z maire 


\file "jour3b/.Readme-Questions"
\brief L'exercise de jour3b

\mainpage jour3b
        
\section s1 GEOMETRIE

   Une navette spatiale est modelisee par un cylindre d'Aluminium.
   A l'interieur de la cabine se trouve un plancher, egalement en Aluminium.
   La navette est placée dans l'espace ('World') : une boite aussi petite que
   possible.
   
   Un astronaute est placé sur le plancher.
        
                
\section s2 COMMENT DEMARRER ?
 
  - lire le code : jour3b.cc (main) et les 5 classes DetectorConstruction,
    DetectorMessenger, PhysicsList, PrimaryGeneratorAction, 
    PrimaryGeneratorMessenger.
 
  - compiler et generer un executable : 
\verbatim  
        % gmake
\endverbatim  
                
  - executer jour3b en mode interactif avec visualisation :
\verbatim  
        % $G4WORKDIR/bin/$G4SYSTEM/jour3b
        ....
        Idle> /run/beamOn 1
        ....
        Idle> exit
\endverbatim  
  
\section s3 CORRIGE et COMMENTAIRES de jour3a
 
   - RunAction.hh : remarquer le constructeur et la fonction in-line AddEdep()
   - RunAction.cc : pour le calcul de la dose, on a besoin de la masse de 
                    l'astronaute. On peut obtenir cette information via
                    logicalVolume de l'astronaute. D'ou la necessité d'un
                    pointeur sur DetectorConstruction.
                    Remarquer aussi les differents modes d'impression.
   - SteppingAction : on a modifié le constructeur pour avoir un pointeur sur
                      RunAction.
   - jour3b.cc : la classe RunAction est enregistrée aupres de RunManager
                
\section s4 EXERCICES
   
   On veut calculer la precision sur la dose cumulée par l'astronaute. 
   Pour cela, on a besoin de la variable aleatoire "energie recue par evenement",
   pour calculer sa valeur moyenne et son ecart-type.
   
   (1 evenement = 1 proton initial qui 'touche' l'astronaute.)
   
   Il faut donc sommer l'energie recue evenement par evenement avant 
   de la sommer sur tout le run. 

   On prendra modele sur examples/extended/electromagnetic/TestEm11
     
\subsection sub1  PARTIE 1 : sommer edep dans EventAction
     
   a) Ajouter une classe EventAction (vide pour l'instant). 
   Dans les data membres privées, prevoir un pointeur sur RunAction 
   
   b) Prevoir la variable EdepPerEvent et la fonction inline AddEdep()
   
   c) Instancier la classe EventAction dans le main programme: jour3b.cc
   
   d) Modifier le constructeur de SteppingAction pour passer un pointeur sur
   EventAction, en plus de celui sur RunAction
   
   e) Modifier la fonction SteppingAction::UserSteppingAction() pour effectuer 
   la somme sur edep dans EventAction
   
   f) Imprimer EdepPerEvent depuis EventAction::EndOfEventAction()

\subsection sub2  PARTIE 2 : sommer EdepPerEvent dans RunAction
  
   Pour calculer la valeur moyenne et la variance de EdepPerEvent, on a besoin
   de la somme sur EdepPerEvent et de la somme des carrés.
   
   a) dans RunAction.hh, prevoir 3 données membres privées :
   sumEdep, sumEdep2 (pour les carrés), nbEntry (pour le nombre d'evenements
   avec EdepPerEvent non nul)
   
   b) prevoir aussi une fonction inline SumEvents()
   (voir dans TestEm11 la fonction inline RunAction::AddEdep(e)).
   Y ajouter l'incrementation du compteur nbEntry.
     
   c) Dans EventAction::EndOfEventAction(), si EdepPerEvent est non nul, la
   passer a RunAction::SumEvents()
   
   d) Dans RunAction::EndOfRunAction(), calculer la moyenne et l'ecart type de
   EdepPerEvent (voir TestEm11).
   
   L'erreur relative sur la dose cumulée est (1/sqrt(nbEntry))*(rms/mean)
   Imprimer la dose cumulée et son erreur.
   
   e) supprimer l'impression dans EventAction et faire une macro pour executer
   un job en mode 'batch' avec 50000 evenements.   

\section s5 SOLUTIONS : voir le code jour3c
  
*/

---