source: trunk/examples/advanced/microbeam/README @ 1319

-------------------------------------------------------------------
$Id: README,v 1.10 2007/02/27 12:02:09 sincerti Exp $
-------------------------------------------------------------------

     =========================================================
                  Geant4 - Microbeam example
     =========================================================

                                README file
                          ----------------------

                           CORRESPONDING AUTHOR 

S. Incerti (a, *) et al.
a. Centre d'Etudes Nucleaires de Bordeaux-Gradignan 
(CENBG), IN2P3 / CNRS / Bordeaux 1 University, 33175 Gradignan, France
* e-mail:incerti@cenbg.in2p3.fr

Last modified by S. Incerti, 27/02/2007

---->0. INTRODUCTION.                                                    
                                                                       
The microbeam example simulates the cellular irradiation beam line 
installed on the AIFIRA electrostatic accelerator facility located at 
CENBG, Bordeaux-Gradignan, France. For more information on this facility, 
please visit :
http://www.cenbg.in2p3.fr/

An overall description of this example is also available in this directory:
to access it, simply open the microbeam.htm file with your internet browser.

---->1. GEOMETRY SET-UP.
 
The elements simulated are:

1. A switching dipole magnet with fringing field, to deflect the 3 MeV alpha 
beam generated by the electrostatic accelerator into the microbeam line, 
oriented at 10 degrees from the main beam direction;

2. A circular collimator object, defining the incident beam size at the 
microbeam line entrance;

3. A quadrupole based magnetic symmetric focusing system allowing equal 
transverse demagnifications of 10. Fringe fields are calculated from Enge's 
model.

4. A dedicated cellular irradiation chamber setup;

5. A set of horizontal and vertical electrostatic deflecting plates which can 
be turned on or off to deflect the beam on target; 

6. A realistic human keratinocyte voxellized cell observed from confocal 
microscopy and taking into account realistic nucleus and cytoplasm chemical 
compositions


---->2. EXPERIMENTAL SET-UP.      
                                 
The beam is defined at the microbeam line entrance through a collimator 
5 micrometer in diameter. The beam is then focused onto target using a 
quadruplet of quadrupoles in the so-called Dymnikov magnetic configuration. 
The beam is sent to the irradiation chamber where it travels through a 
isobutane gas detector for counting purpose before reaching the polypropylene 
culture foil of the target cell which is immersed in the growing medium and 
enclosed within a dish.  

A cell is placed on the polypropylene foil and is irradiated using the 
microbeam. The cell is represented through a 3D phantom (G4PVParameterization) 
obtained from confocal microscopy. In the provided example, the voxels sizes 
are : 359 nm (X) x 359 nm (Y) x 163 nm (Z)

The primary particle beam parameters are generated from experimental 
measurements performed on the AIFIRA facility. Incident particle used for 
cellular irradiation are 3 MeV alpha particles.

More details on the experimental setup and its simulation with Geant4 can 
be found in the following papers, which may be found on the SLAC-SPIRES
online database (http://www.slac.stanford.edu/spires/) :

- MONTE CARLO MICRODOSIMETRY FOR TARGETED IRRADIATION OF INDIVIDUAL CELLS USING 
A MICROBEAM FACILITY 
By S. Incerti, T. Pouthier, H. Seznec, Ph. Moretto, O. Boissonnade, 
T. M. H. Ha, F. Andersson, Ph. Barberet, C. Habchi and D. T. Nguyen 
In preparation (2007)

- MONTE CARLO SIMULATION OF THE CENBG MICROBEAM AND NANOBEAM LINES WITH THE
GEANT4 TOOLKIT
By S. Incerti, Q. Zhang, F. Andersson, Ph. Moretto, G.W. Grime, 
M.J. Merchant, D.T. Nguyen, C. Habchi, T. Pouthier and H. Seznec
In press in Nucl.Instrum.Meth.B, 2007

- A COMPARISON OF CELLULAR IRRADIATION TECHNIQUES WITH ALPHA PARTICLES USING 
THE GEANT4 MONTE CARLO SIMULATION TOOLKIT
By S. Incerti, N. Gault, C. Habchi, J.L.. Lefaix, Ph. Moretto, J.L.. Poncy, 
T. Pouthier, H. Seznec. Dec 2006. 3pp.
Published in Rad.Prot.Dos.,1-3,2006 (Micros 2005 special issue).

- GEANT4 SIMULATION OF THE NEW CENBG MICRO AND NANO PROBES FACILITY
By S. Incerti, C. Habchi, Ph. Moretto, J. Olivier and H. Seznec. May 2006. 5pp.
Published in Nucl.Instrum.Meth.B249:738-742, 2006

- A COMPARISON OF RAY-TRACING SOFTWARE FOR THE DESIGN OF QUADRUPOLE MICROBEAM 
SYSTEMS
By S. Incerti et al., 
Published in Nucl.Instrum.Meth.B231:76-85, 2005

- DEVELOPMENT OF A FOCUSED CHARGED PARTICLE MICROBEAM FOR THE IRRADIATION OF 
INDIVIDUAL CELLS.
By Ph. Barberet, A. Balana, S. Incerti, C. Michelet-Habchi, Ph. Moretto, 
Th. Pouthier. Dec 2004. 6pp. 
Published in Rev.Sci.Instrum.76:015101, 2005

- SIMULATION OF CELLULAR IRRADIATION WITH THE CENBG MICROBEAM LINE USING 
GEANT4.
By S. Incerti, Ph. Barberet, R. Villeneuve, P. Aguer, E. Gontier, 
C. Michelet-Habchi, Ph. Moretto, D.T. Nguyen, T. Pouthier, R.W. Smith. Oct 2003. 6pp. 
Published in IEEE Trans.Nucl.Sci.51:1395-1401, 2004

- SIMULATION OF ION PROPAGATION IN THE MICROBEAM LINE OF CENBG USING 
GEANT4.
By S. Incerti, Ph. Barberet, B. Courtois, C. Michelet-Habchi, 
Ph. Moretto. Sep 2003. 
Published in Nucl.Instrum.Meth.B210:92-97, 2003


---->3. SET-UP 
                                                                        
- a standard Geant4 example GNUmakefile is provided                     

setup with:                                                             
compiler = gcc-3.2.3
G4SYSTEM = linux-g++                                                    

The following section gives the necessary environment variables.                     

------->>3.1  ENVIRONMENT VARIABLES

All variables are defined with their default value.

 - G4SYSTEM = Linux-g++

 - G4INSTALL              points to the installation directory of GEANT4;

 - G4LIB                  point to the compiled libraries of GEANT4;

 - G4WORKDIR              points to the work directory;

 - CLHEP_BASE_DIR         points to the installation directory of CHLEP; 

 - G4LEDATA               points to the low energy electromagnetic libraries;

 - LD_LIBRARY_PATH = $CLHEP_BASE_DIR/lib

 - G4LEVELGAMMADATA       points to the photoevaporation library;

 - NeutronHPCrossSections points to the neutron data files;

 - G4RADIOACTIVEDATA      points to the libraries for radio-active decay 
                          hadronic processes;
 
However, the $G4LEVELGAMMADATA, $NeutronHPCrossSections and $G4RADIOACTIVEDATA
variables do not need to be defined for this example.

Once these variables have been set, simply type gmake to compile the Microbeam
example. 

------->>3.2  VISUALIZATION

The user can visualize the targeted cell with OpenGL, DAWN and vrml, 
as chosen in the microbeam.mac file. OpenGL is the default viewer. The 
cytoplasm in shown in red and the nucleus in green.


---->4. HOW TO RUN THE EXAMPLE                                         

In interactive mode, run:

> $G4WORDIR/bin/Linux-g++/Microbeam

The macro microbeam.mac is executed by default. To get vizualisation, make
sure to uncomment the /vis/... lines in the microbeam.mac macro.
The Microbeam code reads the phantom.dat file containing all the necessary 
information describing the cell phantom. 10 alphas particles are generated.


---->5. PHYSICS

Low energy electromagnetic processes (for alphas, electrons, photons) and 
hadronic elastic and inelastic scattering for alphas are activated by default. 
Low energy electromagnetic electronic and nuclear stopping power are computed 
from ICRU tables.
        

---->6. SIMULATION OUTPUT AND RESULT ANALYZIS                                    

This example does not need any external analysis package. 
The output results consists in several .txt files:

* dose.txt : gives the total deposited dose in the cell nucleus and in the cell 
cytoplasm by each incident alpha particle;

* 3DDose.txt : gives the average on the whole run of the dose deposited per 
Voxel per incident alpha particle;

* range.txt : indicates the final stopping (x,y,z) position of the incident 
alpha particle within the irradiated medium (cell or culture medium)

* stoppingPower.txt : gives the actual stopping power dE/dx of the incident 
alpha particle just before penetrating into the targeted cell;

* beamPosition.txt : gives the beam transverse position distribution(X and Y) 
just before penetrating into the targeted cell;

These files can be easily analyzed using for example the provided ROOT macro 
file plot.C; to do so :
* be sure to have ROOT installed on your machine
* be sure to be in the microbeam directory
* launch ROOT by typing root
* under your ROOT session, type in : .X plot.C to execute the macro file

A graphical output obtained with this macro for 40000 incident alpha particles 
is shown in the file microbeam.gif

The simulation predicts that 95% of the incident alpha particles detected by the
gas detector are located within a circle of 10 um in diameter on the target, in 
nice agreement with experimental measurements performed on the CENBG setup.

---------------------------------------------------------------------------

Should you have any enquiry, please do not hesitate to contact: 
incerti@cenbg.in2p3.fr