Message: Re: internal conversion process with example of radioactivedecay(exrdm) Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: internal conversion process with example of radioactivedecay(exrdm) 

Forum: Hadronic Processes
Re: None internal conversion process with example of radioactivedecay(exrdm) (Yeongduk Kim)
Re: None Re: internal conversion process with example of radioactivedecay(exrdm) (Fan Lei)
Date: 23 Sep, 2007
From: Yeongduk Kim <Yeongduk Kim>

Dear Fan,

Thank you for your comments.

According to your suggestion, I have commented out line 279 of G4VGammaDeexcitation.cc. But the result was the same.

I have attached my macro file, and output of an event for which pb210 decays to bi210[46.5] state. As shown below in detail, there is a line as ; G4VGammaDeexcitation::GenerateGamma - Gamma energy 0.04253199999489397 ** New excitation 0.004006999999999997.

I expected an internal conversion electron with energy 42.5 keV and an x-ray with energy 4.0 keV. But there was only an electron of 42.5 keV and bi210[0.0] from bi210[46.5] state.

Would you comment further on this ?

Thank you. Yeongduk Kim

My macro file :

################################################################################
#
/exrdm/phys/SelectPhysics LowEnergy_EM /run/initialize /run/verbose 0 /event/verbose 0 /tracking/verbose 2

/grdm/BRbias 0 /grdm/fBeta 0 /grdm/verbose 5 /grdm/analogueMC 1 # /histo/fileName one-iso.aida # /gps/position 0 0 0 /gps/energy 2. eV /gps/ion 82 210

/run/beamOn 100 #

Output of the above macro file with exRDM. only one event decaying to bi210[46.5] is shown.

*************************************************************
 Geant4 version Name: geant4-08-00-patch-01    (10-February-2006)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

---> warning from G4Element::GetElement(). The element:  H does not exist in the table. Return NULL pointer.
 Material:Vacuum with formula: H added! 
     Nb of Material = 1
     Nb of Isotope =  0
     Nb of Element =  1

---> warning from G4Element::GetElement(). The element: N does not exist in the table. Return NULL pointer.

---> warning from G4Element::GetElement(). The element:  O does not exist in the table. Return NULL pointer.
 Material:Air with formula: N0.78-O0.22 added! 
     Nb of Material = 2
     Nb of Isotope =  0
     Nb of Element =  3

---> warning from G4Element::GetElement(). The element: Al does not exist in the table. Return NULL pointer.
 Material:Aluminium with formula: Al added! 
     Nb of Material = 3
     Nb of Isotope =  0
     Nb of Element =  4

---> warning from G4Element::GetElement(). The element: Si does not exist in the table. Return NULL pointer.
 Material:Silicon with formula: Si added! 
     Nb of Material = 4
     Nb of Isotope =  0
     Nb of Element =  5

---> warning from G4Element::GetElement(). The element: Pb does not exist in the table. Return NULL pointer.
 Material:Lead with formula: Pb added! 
     Nb of Material = 5
     Nb of Isotope =  0
     Nb of Element =  6

---> warning from G4Element::GetElement(). The element: Ge does not exist in the table. Return NULL pointer.
 Material:Germanium with formula: Ge added! 
     Nb of Material = 6
     Nb of Isotope =  0
     Nb of Element =  7

---> warning from G4Element::GetElement(). The element: Cs does not exist in the table. Return NULL pointer.

---> warning from G4Element::GetElement(). The element:  I does not exist in the table. Return NULL pointer.
 Material:CsI with formula: Cs-I added! 
     Nb of Material = 7
     Nb of Isotope =  0
     Nb of Element =  9
G4VUserPhysicsList::SetVerboseLevel  : Verbose level is set to 2
Visualization Manager instantiating...
Visualization Manager initialising...
Registering graphics systems...

You have successfully chosen to use the following graphics systems.
Current available graphics systems are:
  ASCIITree (ATree)
  DAWNFILE (DAWNFILE)
  GAGTree (GAGTree)
  G4HepRep (HepRepXML)
  G4HepRepFile (HepRepFile)
  RayTracer (RayTracer)
  VRML1FILE (VRML1FILE)
  VRML2FILE (VRML2FILE)
  FukuiRenderer (DAWN)
  OpenGLImmediateX (OGLIX)
  OpenGLStoredX (OGLSX)
  OpenGLImmediateXm (OGLIXm)
  OpenGLStoredXm (OGLSXm)
  RayTracerX (RayTracerX)
  VRML1 (VRML1)
  VRML2 (VRML2)

Registering model factories... exrdmPhysicsList::SelectPhysicsList: <LowEnergy_EM> G4VUserPhysicsList::Construct() Construct processes world cuts are set Target cuts are set Detector cuts are set

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
G4AugerData for Element no. 7 are loaded
G4AugerData for Element no. 8 are loaded
G4AugerData for Element no. 13 are loaded
G4AugerData for Element no. 14 are loaded
G4AugerData for Element no. 32 are loaded
G4AugerData for Element no. 53 are loaded
G4AugerData for Element no. 55 are loaded
G4AugerData for Element no. 82 are loaded
G4AugerData for Element no. 83 are loaded
AugerTransitionTable complete

LowEnergyIoni:  Total cross sections from EEDL database.
      Gamma energy sampled from a parametrised formula.
      Implementation of the continuous dE/dx part.
      At present it can be used for electrons in the energy range [250eV,100GeV].
      The process must work with G4LowEnergyBremsstrahlung.

LowEnBrem:  Total cross sections from EEDL database.
      Gamma energy sampled from a parameterised formula.
      Implementation of the continuous dE/dx part.
      At present it can be used for electrons in the energy range [250eV,100GeV].
      The process must work with G4LowEnergyIonisation.

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.

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

hLowEIoni:    Knock-on electron cross sections . 
        Good description above the mean excitation energy.
        Delta ray energy sampled from  differential Xsection.
        PhysicsTables from 10 eV  to 0.1 TeV  in 360 bins.
        Electronic stopping power model is  ICRU_R49p
        from 1 keV  to 2 MeV .

        Parametrisation model for antiprotons is  ICRU_R49p
        from 25 keV  to 2 MeV .
        Parametrization of the Barkas effect is switched on.
        Nuclear stopping power model is ICRU_R49

msc:  Model variant of multiple scattering for Pb210[0.0]

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 

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 : Air
 Range cuts        :  gamma 1 mum    e- 1 mum    e+ 1 mum
 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 : CsI
 Range cuts        :  gamma 1 mum    e- 1 mum    e+ 1 mum
 Energy thresholds :  gamma 990 eV     e- 1.32069 keV    e+ 1.30453 keV
 Region(s) which use this couple : 
    Target

Index : 2     used in the geometry : Yes     recalculation needed : No 
 Material : Germanium
 Range cuts        :  gamma 1 mum    e- 1 mum    e+ 1 mum
 Energy thresholds :  gamma 990 eV     e- 1.91672 keV    e+ 1.89327 keV
 Region(s) which use this couple : 
    Detector

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

### Run : 0

[comment] 1st event deleted.....

*********************************************************************************************************
* G4Track Information:   Particle = Pb210[0.0],   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0     2e-06        0        0         0      Target initStep
G4RadioactiveDecay::GetMeanFreePath() 
KineticEnergy:1.999999994950486e-09[GeV]Mass:195.5576999980841[GeV]c*Tau:3.043567276416114e+17[m]
mean free path: 1376497543757.474[m]
    1        0        0 1.37e-10         0    2e-06 3.07e-09  3.07e-09      Target hLowEIoni
G4RadioactiveDecay::GetMeanLifeTime() 
KineticEnergy:0[GeV]Mass:195.5576999980841[GeV]Life time: 1.015224764732445e+18[ns]
mean life time: 1.015224764732445e+18[ns]
DecayIt:  Analogue MC version 
Begin of DoDecay...
Selecte a channel...
G4VDecayChannel::FillDaughters()
0:e-:0xab09f38
1:Bi210[46.5]:0xb3ae4d8
2:anti_nu_e:0xab10b40
G4NuclearDecayChannel::DecayIt the decay mass = 195557.7153799711
G4Decay::BetaDecayIt()
     daughter 0:0.0001033450413194893[GeV/c]
     daughter 1:0.0001051242598857712[GeV/c]
     daughter 2:6.215388644023064e-06[GeV/c]
G4NuclearDecayChannel::BetaDecayIt   create decay products in rest frame 
 ----- List of DecayProducts  -----
 ------ Parent Particle ----------
 Particle type - Pb210[0.0]
   mass:        195.5576999980841[GeV]
   charge:      82[e]
   Direction x: 0, y: 0, z: 0
   Total Momentum = 0[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: 0[GeV]
   Total Energy   = 195.5576999980841[GeV]
   Kinetic Energy = 0[GeV]
   ProperTime     = 0[ns]
 ------ Daughter Particles  ------
 ----------1 -------------
 Particle type - e-
   mass:        0.0005109990599999999[GeV]
   charge:      -1[e]
   Direction x: 0.2839294414374162, y: -0.4430019818554305, z: -0.8503724574309763
   Total Momentum = 0.0001033450413194893[GeV]
   Momentum: 2.934269985716929e-05[GeV], y: -4.578205811946511e-05[GeV], z: -8.788177675015989e-05[GeV]
   Total Energy   = 0.0005213446430972611[GeV]
   Kinetic Energy = 1.034558309726109e-05[GeV]
   ProperTime     = 0[ns]
 ----------2 -------------
 Particle type - anti_nu_e
   mass:        0[GeV]
   charge:      0[e]
   Direction x: 0.6238531394629775, y: 0.6497316505949959, z: -0.4343455336449182
   Total Momentum = 6.215388644023064e-06[GeV]
   Momentum: 3.877489718556327e-06[GeV], y: 4.038334722770499e-06[GeV], z: -2.699626297398762e-06[GeV]
   Total Energy   = 6.215388644023064e-06[GeV]
   Kinetic Energy = 6.215388644023064e-06[GeV]
   ProperTime     = 0[ns]
 ----------3 -------------
 Particle type - Bi210[46.5]
   mass:        195.5571878199111[GeV]
   charge:      83[e]
   Direction x: -0.3160087843835754, y: 0.3970893440016001, z: 0.8616603165243215
   Total Momentum = 0.0001051242584969309[GeV]
   Momentum: -3.322018913683987e-05[GeV], y: 4.174372284520091e-05[GeV], z: 9.058140185085006e-05[GeV]
   Total Energy   = 195.5571878199393[GeV]
   Kinetic Energy = 2.825544243023614e-11[GeV]
   ProperTime     = 0[ns]
 ----- End List of DecayProducts  -----

G4ContinuumGammaDeexcitation::CreateTransition - Created
G4ContinuumGammaDeexcitation::CanDoTransition -  Excitation 0.04653899999489397 below max discrete level 3.24431
G4VGammaDeexcitation::DoChain - Transition deleted, end of chain 
 = BreakItUp = 0 gammas from ContinuumDeexcitation 
G4DiscreteGammaDeexcitation::CreateTransition - (A,Z) is valid 
G4DiscreteGammaDeexcitation::CreateTransition - Created from level energy 0.046539, excitation is 0.04653899999489397
G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation 0.04653899999489397, Min-Max are 0.046539 3.24431
G4DiscreteGammaDeexcitation::CanDoTransition -  Halflife 3, Calling from RDM  False  Max-HL = 1000
G4DiscreteGammaDeexcitation::CanDoTransition - CanDo: True 
G4VGammaDeexcitation::GenerateGamma - Gamma energy 0.04253199999489397 ** New excitation 0.004006999999999997
G4VGammaDeexcitation::GenerateGamma -  Gamma fragment generated:   Electron
G4VGammaDeexcitation::Update - Transition deleted 
G4DiscreteGammaDeexcitation::CreateTransition - (A,Z) is valid 
G4DiscreteGammaDeexcitation::CreateTransition - Created from level energy 0.046539, excitation is 0
G4VGammaDeexcitation::Update modified -- ydkim 070923 
G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation <= 0
G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation 0, Min-Max are 0.046539 3.24431
G4DiscreteGammaDeexcitation::CanDoTransition - CanDo: False 
G4VGammaDeexcitation::DoChain - Transition deleted, end of chain 
 = BreakItUp = 1 gammas from DiscreteDeexcitation 
 = BreakItUp = 0 particles from G4AtomicDeexcitation 
 = BreakItUp = Nucleus added to products
*-*-* Photon evaporation: 2
G4RadioactiveDecay::DecayIt : Decay vertex : Time: 4.654116989262185e+16[ns] X:0[cm] Y:0[cm] Z:1.366039968016557e-11[cm]
G4Decay::DecayIt  : decay products in Lab. Frame
 ----- List of DecayProducts  -----
 ------ Parent Particle ----------
 Particle type - Pb210[0.0]
   mass:        195.5576999980841[GeV]
   charge:      82[e]
   Direction x: 0, y: 0, z: 0
   Total Momentum = 0[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: 0[GeV]
   Total Energy   = 195.5576999980841[GeV]
   Kinetic Energy = 0[GeV]
   ProperTime     = 0[ns]
 ------ Daughter Particles  ------
 ----------1 -------------
 Particle type - e-
   mass:        0.0005109990599999999[GeV]
   charge:      -1[e]
   Direction x: 0.2839294414374162, y: -0.4430019818554305, z: -0.8503724574309763
   Total Momentum = 0.0001033450413194893[GeV]
   Momentum: 2.934269985716929e-05[GeV], y: -4.578205811946511e-05[GeV], z: -8.788177675015989e-05[GeV]
   Total Energy   = 0.0005213446430972611[GeV]
   Kinetic Energy = 1.034558309726109e-05[GeV]
   ProperTime     = 0[ns]
 ----------2 -------------
 Particle type - anti_nu_e
   mass:        0[GeV]
   charge:      0[e]
   Direction x: 0.6238531394629775, y: 0.6497316505949959, z: -0.4343455336449182
   Total Momentum = 6.215388644023064e-06[GeV]
   Momentum: 3.877489718556327e-06[GeV], y: 4.038334722770499e-06[GeV], z: -2.699626297398762e-06[GeV]
   Total Energy   = 6.215388644023064e-06[GeV]
   Kinetic Energy = 6.215388644023064e-06[GeV]
   ProperTime     = 0[ns]
 ----------3 -------------
 Particle type - e-
   mass:        0.0005109990599999999[GeV]
   charge:      -1[e]
   Direction x: -0.4662656887707136, y: -0.7360201976994417, z: -0.4907856722171571
   Total Momentum = 0.0002127829764294768[GeV]
   Momentum: -9.921340106357254e-05[GeV], y: -0.0001566125683786992[GeV], z: -0.0001044308361233083[GeV]
   Total Energy   = 0.0005535310599948939[GeV]
   Kinetic Energy = 4.253199999489397e-05[GeV]
   ProperTime     = 3.030743039898416[ns]
 ----------4 -------------
 Particle type - Bi210[0.0]
   mass:        195.5571412809111[GeV]
   charge:      83[e]
   Direction x: -0.3160087843835754, y: 0.3970893440016001, z: 0.8616603165243215
   Total Momentum = 0.0001051242724238209[GeV]
   Momentum: -3.322019353785948e-05[GeV], y: 4.174372837542056e-05[GeV], z: 9.058141385109857e-05[GeV]
   Total Energy   = 195.5571412809393[GeV]
   Kinetic Energy = 2.825545664109086e-11[GeV]
   ProperTime     = 0[ns]
 ----- End List of DecayProducts  -----

    2        0        0 1.37e-10         0        0        0  3.07e-09      Target RadioactiveDecay
    :----- List of 2ndaries - #SpawnInStep=  4(Rest= 4,Along= 0,Post= 0), #SpawnTotal=  4 ---------------
    :         0         0  1.37e-10  2.83e-08         Bi210[0.0]
    :         0         0  1.37e-10    0.0425                 e-
    :         0         0  1.37e-10   0.00622          anti_nu_e
    :         0         0  1.37e-10    0.0103                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

[comment] further tracking follows...

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

 Add Message Add Message
to: "Re: internal conversion process with example of radioactivedecay(exrdm)"

 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 ]