Message: G4 example gives STRANGE energy loss and csda range Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None G4 example gives STRANGE energy loss and csda range 

Forum: Physics List
Date: 11 Jan, 2006
From: Ma yunzhi <Ma yunzhi>

Experts,
  A simple simulation adapted from "geant4.7.0/examples/extended/
electromagnetic/TestEm1" shows STRANGE energy loss and csda range
with low-energy physcics package and penelope package. The simulation 
with standard physics package goes rightly. Relevant aspects are 
listed below. 

Could anyone give some suggest or instruction?
Is there something wrong with the G4LowEnergyIonization class and 
G4PenelopeIonization class ?

Thanks in advance.

os: redhat 9.0
G4: 4.7.1.p01

adaption: When TestEm1 has been compiled, thr micro file run01.mac
was used to calculate the csda range.

DETAILS below



Simulation with STANDARD package ===>>> things OK

//////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////
#run01.mac

/testem/phys/setGCut  95 cm
/testem/phys/setECut  95 cm

/run/initialize

 The Box is 10 m   of Water
PhysicsList::SetCuts:CutLength : 1 mm 

/process/inactivate phot
/process/inactivate compt
/process/inactivate conv
/process/inactivate msc

#/process/eLoss/StepFunction 0.03 0.1 micrometer
/process/eLoss/integral true    <<<= added by me
/process/eLoss/fluct 0

/run/verbose 2

/gun/particle e-
/tracking/verbose 2
/gun/energy 100 keV
/run/beamOn 3

///////////////////////////////////////////////////////////////////////
---> Begin of Event: 0

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    100 keV   143 mum   143 mum       Water       eIoni

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    100 keV   143 mum   143 mum       Water       eIoni

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    100 keV   143 mum   143 mum       Water       eIoni
Run terminated.
Run Summary
  Number of events processed : 3
  User=0s Real=0.01s Sys=0s

 ======================== run summary ======================

 The run was: 3 e- of 100.00 keV through 10.00 m   of Water (density: 1.00 g/cm3 )

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

 total energy deposit: 100.000 keV

////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////


Simulation with LOW ENERGY package ===>>> things STRANGE
////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////

/testem/phys/addPhysics  lowen

/testem/phys/setGCut  95 cm
/testem/phys/setECut  95 cm

/run/initialize

 The Box is 10 m   of Water
PhysicsList::SetCuts:CutLength : 1 mm 

/process/inactivate phot
/process/inactivate compt
/process/inactivate conv
/process/inactivate msc

/process/eLoss/integral true
/process/eLoss/fluct 0

/run/verbose 2

/gun/particle e-
/tracking/verbose 2
/gun/energy 100 keV
/run/beamOn 3
////////////////////////////////////////////////////////////////
---> Begin of Event: 0

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    178 keV   537 mum   537 mum       Water   LowEnBrem
                                        OOOOOOOOOOOOOOOOOOO
                                           dEStep > KineE
                                        OOOOOOOOOOOOOOOOOOO

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    188 keV   537 mum   537 mum       Water   LowEnBrem

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 m       0 fm      0 fm    100 keV     0 eV      0 fm      0 fm        Water    initStep
    1     -5 m       0 fm      0 fm      0 eV    187 keV   537 mum   537 mum       Water   LowEnBrem
Run terminated.
Run Summary
  Number of events processed : 3
  User=0s Real=0s Sys=0s

 ======================== run summary ======================

 The run was: 3 e- of 100.00 keV through 10.00 m   of Water (density: 1.00 g/cm3 )

 ============================================================
////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////

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

 Add Message Add Message
to: "G4 example gives STRANGE energy loss and csda range"

 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 ]