Back to home page

EIC code displayed by LXR

 
 

    


Warning, /geant4/examples/extended/medical/fanoCavity2/README is written in an unsupported language. File is not indexed.

0001 -------------------------------------------------------------------
0002 
0003      =========================================================
0004      Geant4 - an Object-Oriented Toolkit for Simulation in HEP
0005      =========================================================
0006 
0007                             fanoCavity2
0008                             -----------
0009 
0010     This program computes the dose deposited in an ionization chamber by an
0011     extended (one dimensional) monoenergetic electron source.
0012     The geometry of the chamber satisfies the conditions of charged particle
0013     equilibrium. Hence, under idealized conditions, the ratio of the dose 
0014     deposited over the beam energy fluence must be equal to 1.
0015     This variante of the Fano cavity test make use of an reciprocity theorem.
0016     
0017     J.Sempau and P.Andreo, Phys. Med. Biol. 51 (2006) 3533        
0018  
0019  1- GEOMETRY
0020  
0021     The chamber is modelized as a cylinder with a cavity in it.
0022     
0023     5 parameters define the geometry :
0024       - the radius of the chamber (must be big)
0025       - the material of the wall
0026       - the thickness of the wall
0027       - the material of the cavity
0028       - the thickness of the cavity
0029 
0030     Wall and cavity must be made of the same material, but with different
0031     density.
0032     Radius must be bigger than range of electrons in cavity.            
0033 
0034     All above parameters can be redifined via the UI commands built in 
0035     DetectorMessenger class.
0036 
0037                         _________________
0038      radius (infinite)  |     |   |     |
0039                         |     |   |     |
0040                         |     |   |     |
0041                         |     |   |     |
0042                         |     | <-+-----+--- cavity
0043                         |     |   |     |
0044                         |     |   |     |
0045                  ---------------------------- cylinder axis = e- source
0046                         |     |   |     |
0047                         |     |   |     |
0048                         |     |   |     |
0049                         |wall |   |wall |
0050                         |     |   |     |
0051                         |     |   |     |
0052                         |     |   |     |
0053                         -----------------
0054 
0055  2- BEAM
0056   
0057     Monoenergetic (E0) incident electron source is uniformly distribued along
0058     cylinder axis, within wall and cavity, with constant lineic density
0059     per mass: I.
0060     An effective wall thickness is defined from the range of e- at energy E0.
0061      
0062     Beam_energy_fluence is E0*I
0063     
0064  3- PURPOSE OF THE PROGRAM
0065     
0066     The program computes the dose deposited in the cavity and the ratio
0067     Dose/Beam_energy_fluence. This ratio must be 1.
0068  
0069     The program needs high statistic to reach precision on the computed dose.
0070     The UI command /run/printProgress allows to survey the convergence of
0071     the dose calculation.
0072     
0073     The simplest way to study the effect of e- tracking parameters on dose 
0074     deposition is to use the command /testem/stepMax.
0075     
0076  4- PHYSICS
0077  
0078     The physics list contains the standard electromagnetic processes, with few 
0079     modifications listed here.
0080     
0081     - Bremsstrahlung : Fano conditions imply no energy transfer via
0082     bremsstrahlung radiation. Therefore this process is not registered in the
0083     physics list. However, it is always possible to include it. 
0084     See PhysListEm classes.
0085     
0086     - Ionization : In order to have same stopping power in wall and cavity, one
0087     must cancel the density correction term in the dedx formula. This is done in
0088     a specific MollerBhabha model (MyMollerBhabhaModel) which inherites from 
0089     G4MollerBhabhaModel.
0090     
0091     To prevent explicit generation of delta-rays, the default production
0092     threshold (i.e. cut) is set to 10 km (CSDA condition).
0093     
0094     The finalRange of the step function is set to 10 um, which more on less
0095     correspond to a tracking cut in water of about 20 keV. See emOptions.
0096     Once again, the above parameters can be controled via UI commands.
0097     
0098     - Multiple scattering : is switched in single Coulomb scattering mode near
0099     boundaries. This is selected via EM options in PhysicsList, and can be
0100     controled with UI commands.
0101     
0102     - All PhysicsTables are built with 100 bins per decade.  
0103         
0104  5- HISTOGRAMS
0105  
0106    fanoCavity2 has several predefined 1D histograms : 
0107   
0108       1 : emission point of e+-
0109       2 : energy spectrum of e+-
0110       3 : theta distribution of e+-
0111       4 : emission point of e+- hitting cavity
0112       5 : energy spectrum of e+- when entering in cavity
0113       6 : theta distribution of e+- before enter in cavity
0114       7 : theta distribution of e+- at first step in cavity      
0115       8 : track segment of e+- in cavity
0116       9 : step size of e+- in wall
0117      10 : step size of e+- in cavity
0118      11 : energy deposit in cavity per track
0119    
0120    The histograms are managed by G4AnalysisManager class and its Messenger. 
0121    The histos can be individually activated with the command :
0122    /analysis/h1/set id nbBins  valMin valMax unit 
0123    where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
0124    
0125    One can control the name of the histograms file with the command:
0126    /analysis/setFileName  name  (default fanocavity2)
0127    
0128    It is possible to choose the format of the histogram file : root (default),
0129    hdf5, xml, csv, by changing the default file type in HistoManager.cc
0130    
0131    It is also possible to print selected histograms on an ascii file:
0132    /analysis/h1/setAscii id
0133    All selected histos will be written on a file name.ascii 
0134    (default fanocavity2) 
0135  
0136  6- HOW TO START ?
0137  
0138     - execute fanoCavity2 in 'batch' mode from macro files
0139         % fanoCavity2   run01.mac
0140  
0141     - execute fanoCavity2 in 'interactive mode' with visualization
0142         % fanoCavity2
0143         ....
0144         Idle> type your commands
0145         ....
0146         Idle> exit
0147 
0148    Alternative macro files:
0149    basic.mac - disabled  multiple scattering and fluctuations of energy loss
0150    essai.mac - run WVI EM physics configuration 
0151    stepfunction.mac - the step function optimisation using histogram