FEDRA VMC ROOT On Screen

From AstroNuWiki
Jump to: navigation, search

naumov@naumov-laptop(17:35)> root

 *******************************************
 *                                         *
 *        W E L C O M E  to  R O O T       *
 *                                         *
 *   Version   5.13/01      11 July 2006   *
 *                                         *
 *  You are welcome to visit our Web site  *
 *          http://root.cern.ch          *
 *                                         *
 *******************************************
FreeType Engine v2.1.9 used to render TrueType fonts. 
Compiled on 20 July 2006 for linux with thread support.
CINT/ROOT C/C++ Interpreter version 5.16.13, June 8, 2006
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
Loading root libraries:
Note: File "/cern/root/lib/libGui.so" already loaded
        LibGeom, LibVMC,LibGui,LibGed
Loading FEDRA libraries:
libEmath        loaded...
libEmath, vt++ version /1-35 loaded (Fri Jan 18 14:33:54 2002).
libvt, libEdb, libEdr
Loading EVMC libraries:
       LibVMCGeometry, LibVMCInput, LibVMCOutput, LibVMCKernel, LibVMCEvent, LibVMCViewer, LibVMCTools,    
LibRecoKernel
Set macros path: ./macros
root [0] .x run_g4.C
Loading Geant4 granular libraries                    ... done
Loading VGM libraries                                ... done
Loading Geant4 VMC libraries                         ... done
*************************************************************
 Geant4 version Name: geant4-08-01    (30-June-2006)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************
Visualization Manager instantiating...
Visualization Manager initialising...
Registering graphics systems...
You have successfully registered the following graphics systems.
Current available graphics systems are:
 ASCIITree (ATree)
 DAWNFILE (DAWNFILE)
 GAGTree (GAGTree)
 G4HepRepFile (HepRepFile)
 G4HepRep (HepRepXML)
 RayTracer (RayTracer)
 VRML1FILE (VRML1FILE)
 VRML2FILE (VRML2FILE)
 OpenGLImmediateX (OGLIX)
 OpenGLStoredX (OGLSX)
 OpenGLImmediateXm (OGLIXm)
 OpenGLStoredXm (OGLSXm)
Registering model factories...
You have successfully registered the following model factories.
Registered model factories:
 drawByCharge
 drawByParticleID
Registered models:
 None
Registered filter factories:
 None
Registered filters:
 None
Geant4 has been created.
Info in <EMCApplication::ConstructGeometry>: Geometry will be defined via TGeo
Info in <TGeoManager::TGeoManager>: Geometry DetectorGeometry, DetectorGeometry created
Info in <EdbBrickGeometry::EdbBrickGeometry()>: EdbBrickGeometry is enabled with the following settings
--->   Half of lenght along X is              5.1
--->   Half of lenght along Y is              6.35
--->   Width of first emultion layer in film  0.0045
--->   Width of base layer in film            0.02
--->   Width of second emultion layer in film 0.0045
--->   Width of lead layer in brick           0.1
--->   Width of last plastic layer in film    0.2
--->   Number of emultion films in the brick  58
--->   Shift between plates                   0
Info in <EdbBrickGeometry::MakeMaterials()>: Sum of weights for plastic is 1.000000
Info in <EdbBrickGeometry::MakeMaterials()>: Sum of weights for emulsion is 1.000000
Info in <TGeoManager::SetTopVolume>: Top volume is TOP. Master volume is TOP
Info in <TGeoManager::CheckGeometry>: Fixing runtime shapes...
Info in <TGeoManager::CheckGeometry>: ...Nothing to fix
Info in <TGeoManager::CloseGeometry>: Counting nodes...
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CloseGeometry>: Building cache...
Info in <TGeoManager::BuildCache>: --- Maximum geometry depth is 3
Info in <TGeoManager::CloseGeometry>: 285 nodes/ 285 volume UID's in DetectorGeometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
### General physics constructed.
### Electromagnetic physics constructed.
### Muon physics constructed.
### EM physics for ions constructed.
### Hadron physics for ions constructed.
###
You are using the Geant4 VMC hadronics physics
equivalent to ExN04HadronPhysics from Geant4 (version 6.0)
Note that this hadronic physics list is not optimized for any particular usage.
If you wish to have a starting point tailored for a particular area of work,
please use one of the available physics lists by use-case.
More information can also be found from the Geant4 HyperNews.
###
### EM physics for hadrons constructed.
### Hadron physics for hadrons constructed.
### Step limiter physics constructed.
**********************************************
    GEANT4 Geometry statistics:
            285 logical volumes
            284 physical volumes
              6 materials
              5 user limits
            285 sensitive detectors
**********************************************
Info in <OperaHbookGen::OperaHbookGen>: Added file 1) /home/naumov/work/fedra/data/numucc_dis/1/neg.root
Info in <OperaHbookGen::OperaHbookGen>: Added file 2) /home/naumov/work/fedra/data/numucc_dis/2/neg.root
--- Drawing      227 nodes with 5 visible levels
<TCanvas::MakeDefCanvas>: created default TCanvas with name c1
phot:  Total cross sections from Sandia parametrisation.
compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
      Sampling according Klein-Nishina model
     tables are built for  gamma
     Lambda tables from 100 eV  to 100 GeV in 90 bins.
conv:  Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
     sampling secondary e+e- according to the Bethe-Heitler model
     tables are built for  gamma
     Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
msc:  Model variant of multiple scattering for e-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1
eIoni:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
eBrem:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Total cross sections from a parametrisation based on the EEDL data library.
      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV.
eIoni:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
eBrem:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Total cross sections from a parametrisation based on the EEDL data library.
      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV.
annihil:       Heilter model of formula of annihilation into 2 photons
      tables are built for  e+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
 hIoni:   tables are built for  proton
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used to proton dE/dx and range
      Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc:  Model variant of multiple scattering for proton
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1 
msc:  Model variant of multiple scattering for GenericIon
     Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1
hIoni:   tables are built for  anti_proton
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used to proton dE/dx and range
      Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc:  Model variant of multiple scattering for mu+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1
muIoni:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
muBrems:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model
muPairProd:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model
muIoni:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
 muBrems:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model
 muPairProd:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model
 hIoni:   tables are built for  pi+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used to proton dE/dx and range
      Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
 hIoni:   tables are built for  pi-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used to proton dE/dx and range
      Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49 parametrisation for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc:  Model variant of multiple scattering for pi-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1
========= Table of registered couples ==============================
Index : 0     used in the geometry : Yes     recalculation needed : No
 Material : TOP
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
 Region(s) which use this couple :
   DefaultRegionForTheWorld
Index : 1     used in the geometry : Yes     recalculation needed : No
 Material : Vacuum
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
 Region(s) which use this couple :
   DefaultRegionForTheWorld
Index : 2     used in the geometry : Yes     recalculation needed : No
 Material : Emulsion
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 15.1169 keV    e- 514.737 keV    e+ 496.074 keV
 Region(s) which use this couple :
   DefaultRegionForTheWorld
Index : 3     used in the geometry : Yes     recalculation needed : No
 Material : Plastic
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 2.36895 keV    e- 351.438 keV    e+ 342.891 keV
 Region(s) which use this couple :
   DefaultRegionForTheWorld
Index : 4     used in the geometry : Yes     recalculation needed : No
 Material : Pb
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 100.511 keV    e- 1.37814 MeV    e+ 1.28002 MeV
 Region(s) which use this couple :
   DefaultRegionForTheWorld
====================================================================
### Run 0 start.
---> Begin of event: 1---> Generated Energy 17.531763 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 2---> Generated Energy 24.619991 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 3---> Generated Energy 20.860819 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 4---> Generated Energy 189.815521 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 5---> Generated Energy 22.343647 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 6---> Generated Energy 12.526907 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 7---> Generated Energy 17.070768 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 8---> Generated Energy 60.537205 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 9---> Generated Energy 20.140907 GeV
--- Drawing      227 nodes with 5 visible levels
---> Begin of event: 10---> Generated Energy 13.580837 GeV 
--- Drawing      227 nodes with 5 visible levels
Time of this run:   User=25.42s Real=29.35s Sys=0.38s
Number of events processed: 10