Message: Re: How is Base Material Optimized? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: How is Base Material Optimized? 

Forum: Geometry
Re: None How is Base Material Optimized? (Youming Yang)
Re: None RE: How is Base Material Optimized? (Vladimir Ivanchenko)
Re: None Re: How is Base Material Optimized? (Youming Yang)
Date: 21 Jun, 2012
From: Youming Yang <Youming Yang>

Hello John and VI,
I believe I might be getting mixed conclusions from the two responses (which I am very grateful for!) that I have received:

2a.  For a parameterized geometry where the Material scan needs to create all 
the tables, how should I give it materials?  Should I give it my 18 base 
materials, or should I give it my 4000+?

VI said:
>all 4000 materials should be known to geometry. Even if you do not use 
>base materials approach you need this.

And John:
>You should need to inform the G4 system only of the base materials which you are using in the parameterisation.
>The EM physics tables will include entries only to the base materials.
>In case you find evidence otherwise, we would be very interested to see it.

Feeding GetMaterial all 4000+ materials took prohibitively long to initialize and made it impossible to debug, so I attempted the second suggestion. My code is implemented such that:

-I pass two G4Materials** vectors into my parameterisation: 28 base materials (fbMaterials), and 4072 materials created with BuildMaterialWithNewDensity (fMaterials)

The relevant pieces of code within my parameterisation include: -GetNumberofMaterials returns the number of base materials (28) -GetMaterial returns the base material fbMaterial[i] -ComputeMaterial can return one of the 4072 materials fMaterial[i]

My output (with verbosity 4) is appended. It initializes and loads very quickly, but seems to crash when it receives material1076 (which is built off of a base material)
John, was I incorrect in my interpretation of your response?

I have attached the text output of the crash, as well as a screenshot of where it crashes within the kernel. As with one of my previous posts, the material couple that was input is the null pointer 0x00000000

Thanks, Ming

 [ MIME part of type text/html without a name stripped ]

*************************************************************
 Geant4 version Name: geant4-09-05-ref-00    (2-December-2011)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

Verbosity in physics: 0
Visualization Manager instantiating with verbosity "warnings (3)"...
Visualization Manager initialising...
Registering graphics systems...

You have successfully registered the following graphics systems.
Current available graphics systems are:
  ASCIITree (ATree)
  DAWNFILE (DAWNFILE)
  G4HepRep (HepRepXML)
  G4HepRepFile (HepRepFile)
  RayTracer (RayTracer)
  VRML1FILE (VRML1FILE)
  VRML2FILE (VRML2FILE)
  gMocrenFile (gMocrenFile)

Registering model factories...

You have successfully registered the following model factories.
Registered model factories:
  generic
  drawByCharge
  drawByParticleID
  drawByOriginVolume
  drawByAttribute

Registered filter factories:
  chargeFilter
  particleFilter
  originVolumeFilter
  attributeFilter

Some /vis commands (optionally) take a string to specify colour.
Available colours:
  black, blue, cyan, gray, green, grey, magenta, red, white, yellow

/run/verbose 1
/event/verbose 1
/tracking/verbose 1
/geometry/navigator/verbose 1
/physics/addPhysics LowE_Livermore
/run/initialize
Number of VoxX: 512
Number of VoxY: 512
Number of VoxZ: 121
world size: (522,522,131)
ctx[cm]: 51.2
cty[cm]: 51.2
ctz[cm]: 12.1
mat: 0mat: 1mat: 2mat: 3mat: 4mat: 5mat: 6mat: 7mat: 8mat: 9mat: 10mat: 11mat: 12mat: 13mat: 14mat: 15mat: 16mat: 17mat: 18mat: 19mat: 20mat: 21mat: 2
2mat: 23mat: 24mat: 25mat: 26mat: 27mat: 28### G4EmConfigurator::AddModels n= 0
/gps/particle e-
/gps/pos/type Plane
/gps/pos/shape Square
/gps/pos/centre 0. 0. 0. cm
/gps/pos/halfx 0. cm
/gps/pos/halfz 0. cm
/gps/direction 0. 0. 1.
/gps/ene/mono 20 MeV
/analysis/filename testout.out
/analysis/numvox 512 512 121
/run/beamOn 1
mat: 0mat: 1mat: 2mat: 3mat: 4mat: 5mat: 6mat: 7mat: 8mat: 9mat: 10mat: 11mat: 12mat: 13mat: 14mat: 15mat: 16mat: 17mat: 18mat: 19mat: 20mat: 21mat: 2
2mat: 23mat: 24mat: 25mat: 26mat: 27mat: 28
phot:   for  gamma    SubType= 12
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
 LivermorePhElectric :  Emin=        0 eV    Emax=        1 GeV  FluoActive
       PhotoElectric :  Emin=        1 GeV   Emax=       10 TeV  FluoActive

compt:   for  gamma    SubType= 13
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
    LivermoreCompton :  Emin=        0 eV    Emax=        1 GeV  FluoActive
       Klein-Nishina :  Emin=        1 GeV   Emax=       10 TeV

conv:   for  gamma    SubType= 14
      Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
 LivermoreConversion :  Emin=        0 eV    Emax=        1 GeV
        BetheHeitler :  Emin=        1 GeV   Emax=       10 TeV

Rayl:   for  gamma    SubType= 11
      Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh :  Emin=        0 eV    Emax=        1 GeV
   LivermoreRayleigh :  Emin=        1 GeV   Emax=      100 GeV

msc:   for e-    SubType= 10
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          UrbanMsc95 :  Emin=        0 eV    Emax=       10 TeV

### ===  Deexcitation model UAtomDeexcitation is activated for regions:
          DefaultRegionForTheWorld

eIoni:   for  e-    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
       LowEnergyIoni :  Emin=        0 eV    Emax=      100 keV
        MollerBhabha :  Emin=      100 keV   Emax=       10 TeV

eBrem:   for  e-    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      LPM flag: 1 for E > 1 GeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           LowEnBrem :  Emin=        0 eV    Emax=       25 MeV   AngularGen2BS
             eBremSB :  Emin=       25 MeV   Emax=        1 GeV   AngularGenUrban
            eBremLPM :  Emin=        1 GeV   Emax=       10 TeV   AngularGenUrban

eIoni:   for  e+    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        MollerBhabha :  Emin=        0 eV    Emax=       10 TeV

eBrem:   for  e+    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      LPM flag: 1 for E > 1 GeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB :  Emin=        0 eV    Emax=        1 GeV   AngularGenUrban
            eBremLPM :  Emin=        1 GeV   Emax=       10 TeV   AngularGenUrban

annihil:   for  e+    SubType= 5
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            eplus2gg :  Emin=        0 eV    Emax=       10 TeV

msc:   for proton    SubType= 10
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          UrbanMsc90 :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  proton    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg :  Emin=        0 eV    Emax=        2 MeV
          BetheBloch :  Emin=        2 MeV   Emax=       10 TeV

hBrems:   for  proton    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem :  Emin=        0 eV    Emax=       10 TeV

hPairProd:   for  proton    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd :  Emin=        0 eV    Emax=       10 TeV

msc:   for GenericIon    SubType= 10
      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          UrbanMsc90 :  Emin=        0 eV    Emax=       10 TeV

ionIoni:   for  GenericIon    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
      Stopping Power data for 17 ion/material pairs
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
         ParamICRU73 :  Emin=        0 eV    Emax=       10 TeV

nuclearStopping:   for  GenericIon    SubType= 8
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   ICRU49NucStopping :  Emin=        0 eV    Emax=       10 TeV

nuclearStopping:   for  alpha    SubType= 8
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   ICRU49NucStopping :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  anti_proton    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO :  Emin=        0 eV    Emax=        2 MeV
          BetheBloch :  Emin=        2 MeV   Emax=       10 TeV

msc:   for kaon+    SubType= 10
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          UrbanMsc90 :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  kaon+    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg :  Emin=        0 eV    Emax=  1.05231 MeV
          BetheBloch :  Emin=  1.05231 MeV   Emax=       10 TeV

hBrems:   for  kaon+    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem :  Emin=        0 eV    Emax=       10 TeV

hPairProd:   for  kaon+    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  kaon-    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO :  Emin=        0 eV    Emax=  1.05231 MeV
          BetheBloch :  Emin=  1.05231 MeV   Emax=       10 TeV

hBrems:   for  kaon-    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem :  Emin=        0 eV    Emax=       10 TeV

hPairProd:   for  kaon-    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd :  Emin=        0 eV    Emax=       10 TeV

muMsc:   for mu+    SubType= 10
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni :  Emin=        0 eV    Emax=       10 TeV

muIoni:   for  mu+    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg :  Emin=        0 eV    Emax=      200 keV
          BetheBloch :  Emin=      200 keV   Emax=        1 GeV
        MuBetheBloch :  Emin=        1 GeV   Emax=       10 TeV

muBrems:   for  mu+    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem :  Emin=        0 eV    Emax=       10 TeV

muPairProd:   for  mu+    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd :  Emin=        0 eV    Emax=       10 TeV

CoulombScat:   for  mu+    SubType= 1
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      180 < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering :  Emin=        0 eV    Emax=       10 TeV

muIoni:   for  mu-    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO :  Emin=        0 eV    Emax=      200 keV
          BetheBloch :  Emin=      200 keV   Emax=        1 GeV
        MuBetheBloch :  Emin=        1 GeV   Emax=       10 TeV

muBrems:   for  mu-    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem :  Emin=        0 eV    Emax=       10 TeV

muPairProd:   for  mu-    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd :  Emin=        0 eV    Emax=       10 TeV

CoulombScat:   for  mu-    SubType= 1
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      180 < Theta(degree) < 180; pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  pi+    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg :  Emin=        0 eV    Emax=  297.505 keV
          BetheBloch :  Emin=  297.505 keV   Emax=       10 TeV

hBrems:   for  pi+    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem :  Emin=        0 eV    Emax=       10 TeV

hPairProd:   for  pi+    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd :  Emin=        0 eV    Emax=       10 TeV

msc:   for pi-    SubType= 10
      Lambda tables from 100 eV  to 10 TeV in 220 bins, spline: 1
      RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          UrbanMsc90 :  Emin=        0 eV    Emax=       10 TeV

hIoni:   for  pi-    SubType= 2
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO :  Emin=        0 eV    Emax=  297.505 keV
          BetheBloch :  Emin=  297.505 keV   Emax=       10 TeV

hBrems:   for  pi-    SubType= 3
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem :  Emin=        0 eV    Emax=       10 TeV

hPairProd:   for  pi-    SubType= 4
      dE/dx and range tables from 100 eV  to 10 TeV in 220 bins
      Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd :  Emin=        0 eV    Emax=       10 TeV

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No
 Material : Water
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.94056 keV    e-  351.877 keV    e+  342.545 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No
 Material : material0
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 2     used in the geometry : Yes     recalculation needed : No
 Material : material50
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  990 eV     e-  57.3324 keV    e+  56.567 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 3     used in the geometry : Yes     recalculation needed : No
 Material : material880
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.3478 keV    e-  324.616 keV    e+  316.006 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 4     used in the geometry : Yes     recalculation needed : No
 Material : material917
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.46088 keV    e-  333.46 keV    e+  324.616 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 5     used in the geometry : Yes     recalculation needed : No
 Material : material947
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.58809 keV    e-  337.972 keV    e+  331.226 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 6     used in the geometry : Yes     recalculation needed : No
 Material : material977
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.67652 keV    e-  344.855 keV    e+  337.972 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 7     used in the geometry : Yes     recalculation needed : No
 Material : material1007
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.78665 keV    e-  351.877 keV    e+  342.545 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 8     used in the geometry : Yes     recalculation needed : No
 Material : material1018
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.91107 keV    e-  356.639 keV    e+  344.855 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 9     used in the geometry : Yes     recalculation needed : No
 Material : material1080
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  2.92086 keV    e-  363.901 keV    e+  356.639 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 10     used in the geometry : Yes     recalculation needed : No
 Material : material1120
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.00044 keV    e-  371.312 keV    e+  356.639 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 11     used in the geometry : Yes     recalculation needed : No
 Material : material1200
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.18751 keV    e-  376.336 keV    e+  368.825 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 12     used in the geometry : Yes     recalculation needed : No
 Material : material1300
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.38624 keV    e-  391.82 keV    e+  376.336 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 13     used in the geometry : Yes     recalculation needed : No
 Material : material1400
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.60947 keV    e-  402.496 keV    e+  389.196 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 14     used in the geometry : Yes     recalculation needed : No
 Material : material1500
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.77063 keV    e-  413.461 keV    e+  397.122 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 15     used in the geometry : Yes     recalculation needed : No
 Material : material1600
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  3.97889 keV    e-  424.726 keV    e+  410.692 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 16     used in the geometry : Yes     recalculation needed : No
 Material : material1700
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  4.17053 keV    e-  436.298 keV    e+  419.056 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 17     used in the geometry : Yes     recalculation needed : No
 Material : material1800
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  4.34212 keV    e-  442.201 keV    e+  433.376 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 18     used in the geometry : Yes     recalculation needed : No
 Material : material1900
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  4.52077 keV    e-  457.312 keV    e+  442.201 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 19     used in the geometry : Yes     recalculation needed : No
 Material : material2000
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  4.64394 keV    e-  466.625 keV    e+  448.185 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 20     used in the geometry : Yes     recalculation needed : No
 Material : material2100
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  4.83501 keV    e-  476.127 keV    e+  460.395 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 21     used in the geometry : Yes     recalculation needed : No
 Material : material2200
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.00023 keV    e-  485.824 keV    e+  472.939 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 22     used in the geometry : Yes     recalculation needed : No
 Material : material2300
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.13646 keV    e-  499.06 keV    e+  482.57 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 23     used in the geometry : Yes     recalculation needed : No
 Material : material2400
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.2764 keV    e-  512.657 keV    e+  492.397 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 24     used in the geometry : Yes     recalculation needed : No
 Material : material2500
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.42016 keV    e-  519.593 keV    e+  499.06 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 25     used in the geometry : Yes     recalculation needed : No
 Material : material2600
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.4935 keV    e-  530.175 keV    e+  512.657 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 26     used in the geometry : Yes     recalculation needed : No
 Material : material3000
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  5.75809 keV    e-  570.85 keV    e+  548.291 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 27     used in the geometry : Yes     recalculation needed : No
 Material : Titanium
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  13.2486 keV    e-  815.083 keV    e+  782.873 keV proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 28     used in the geometry : Yes     recalculation needed : No
 Material : CT_Air
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 29     used in the geometry : Yes     recalculation needed : No
 Material : vacuum
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm
 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

====================================================================

### Run 0 start.
/vis/scene/notifyHandlers
WARNING: Scene "refresh" not found.
  /vis/scene/list to see scenes.
Starting allocate memory for dose output
Start Run processing.
HU is: 71. Material is: material1071
=====================================
  G4EventManager::ProcessOneEvent()
=====================================

---> Begin of event: 0
G4PrimaryTransformer::PrimaryVertex (0(mm),0(mm),0(mm),0(nsec))
1 primaries are passed from G4EventTransformer.
!!!!!!! Now start processing an event !!!!!!!

*********************************************************************************************************
* G4Track Information:   Particle = e-,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng   Density      Volume     Process
    0      0 fm      0 fm      0 fm     20 MeV     0 eV      0 fm      0 fm         RepX    initStep
*** G4Navigator::ComputeStep: ***
    Volume = RepX - Proposed step length = 113.673
*** G4Navigator::ComputeStep: ***
    Volume = RepX - Proposed step length = 1.78122
    1      0 fm      0 fm   1.78 mm   19.7 MeV   246 keV  1.78 mm   1.78 mm         RepX     1 g/cc         eBrem
Copy Number is: 256
History Depth is: 3

*** G4Navigator::ComputeStep: ***
    Volume = RepX - Proposed step length = 112.227
*** G4Navigator::ComputeStep: ***
    Volume = RepX - Proposed step length = 0.00150415
HU is: 80. Material is: material1080
    2  -35.9 fm  -1.42 fm   1.78 mm   19.7 MeV0.000174 eV 1e+003 fm   1.78 mm         RepX     1 g/cc    Transportation
Copy Number is: 256
History Depth is: 3

*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 106.064
*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 0.152636
    3  -5.34 um  -59.4 nm   1.93 mm   19.7 MeV    24 keV   153 um   1.94 mm     ParamVox  1.08 g/cc           msc
Copy Number is: 62
History Depth is: 4

*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 105.933
*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 0.293533
    4  -14.5 um   -220 nm   2.23 mm   19.6 MeV  87.7 keV   294 um   2.23 mm     ParamVox  1.08 g/cc           msc
Copy Number is: 62
History Depth is: 4

*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 105.458
*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 0.271428
    5    -21 um  -10.6 nm    2.5 mm   19.6 MeV  53.1 keV   271 um    2.5 mm     ParamVox  1.08 g/cc           msc
Copy Number is: 62
History Depth is: 4

*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 105.17
*** G4Navigator::ComputeStep: ***
    Volume = ParamVox - Proposed step length = 0.00144106
HU is: 76. Material is: material1076
    6    -21 um 1.27e-005 fm    2.5 mm   19.6 MeV  95.5 eV    790 nm    2.5 mm     ParamVox  1.07 g/cc    Transportation
Copy Number is: 62
History Depth is: 4


   Attachment:
     http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/06/21/10.50-2707-crash.png

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: How is Base Material Optimized?   (John Apostolakis - 21 Jun, 2012)
1 None: Re: How is Base Material Optimized?   (Ming - 21 Jun, 2012)
2 None: Re: How is Base Material Optimized?   (Youming Yang - 22 Jun, 2012)
1 None: Re: How is Base Material Optimized?   (John Apostolakis - 25 Jun, 2012)
3 None: Re: How is Base Material Optimized?   (Ming - 22 Jun, 2012)
 Add Message Add Message
to: "Re: How is Base Material Optimized?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]