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

     =========================================================
 		  Geant4 - Nanobeam 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

---->1. INTRODUCTION.                                                    
                                                                       
The nanobeam example simulates the beam optics of the nanobeam 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/

The code can be used to calculate : 
1) intrinsic aberration coefficients of the nanobeam line
2) beam image from a relasitic primary emittance distribution
3) grid shadow images

Three quadrupole field models can be used :
- a simple square field model
- a 3D mesh field model computed from OPERA3D
- an analytical model based on Enge's model

---->2. GEOMETRY SET-UP.
 
The full magnetic configuration of the nanobeam line is simulated.
This configuration is made of a combination of a doublet and triplet of 
5 Oxford Microbeams Ltd. OM50 quadrupoles. 

More details on the experimental setup and its simulation with Geant4 can 
be found in the following papers:

- A DETAILED RAY-TRACING SIMULATION OF THE HIGH RESOLUTION MICROBEAM AT THE
AIFIRA FACILITY
By F. Andersson, Ph. Barberet, S. Incerti, Ph. Moretto
Published in Nucl.Instrum.Meth.B266:1653-1658, 2008

- 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
Published in Nucl.Instrum.Meth.B260:20-27, 2007

- GEANT4 SIMULATION OF THE NEW CENBG MICRO AND NANO PROBES FACILITY
By S. Incerti, C. Habchi, Ph. Moretto, J. Olivier and H. Seznec
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

---->3  VISUALIZATION

Visualization has not been implemented. 

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

1) You must have compiled your Geant4 installation with the FULL version of the
CLHEP library which can handle matrix operations.

2) The code should be compiled cmake and run with :

./nanobeam 

The macro file default.mac is read by default.

Several macro files are provided:

1) for the computation of intrinsic aberration coefficients :
coef-square.mac : using square magnetic field model 
coef-map.mac : using 3D map magnetic field model
coef-enge.mac : using Enge's analytical field model

2) for the simulation of the beam image with a realistic emittance :
image-square.mac : using square magnetic field model (=default.mac)
image-map.mac : using 3D map magnetic field model
image-enge.mac : using Enge's analytical field model

3) for the simulation of grid shadow images
grid-square.mac : using square magnetic field model
grid-map.mac : using 3D map magnetic field model
grid-enge.mac : using Enge's analytical field model

These macros files are stored in the ./macros directory.

To run macros which include *map* in their name, copy the file OM50.grid
into the directory in which you run ./nanobeam.

The code can be run in MT mode, for high statistics image simulation.
Do not use MT for aberration coefficients calculation (32 rays only are shot).
The switch to MT can be made in nanobeam.cc.

---->5. PHYSICS

The example runs with protons with fluctuating energies around 3 MeV.
Standard electromagnetic processes are activated by default (corresponding to the
Physics builder G4EmStandardPhysics), including the G4StepLimiter process. 

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

All results are stored in the nanobeam.root file and can be displayed with the provided
ROOT macro file plot.C: 
* be sure to have ROOT installed on your machine
* be sure to be in the directory where ROOT output files are generated
* copy plot.C into this directory
* launch ROOT by typing root, then under your ROOT session, type in : .X plot.C 
  to execute the macro file
* or type directly: root plot.X

This macro file shows :
- the beam profile along the nanobeam line (only for the computation of intrinsic 
  coefficients)
- the beam image (Y vs X) on target
- the beam emittance (THETA vs X) and (PHY vs Y) on target
- the grid shadow image (option)

The output ntuples can be written as xml or csv files, by changing the G4AnalysisManager default file type in RunAction::BeginOfRunAction().

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

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