Message: BUT instruction of message 183 doesn't give the right SCDA range Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None BUT instruction of message 183 doesn't give the right SCDA range 

Forum: Physics List
Re: None Message no. 197 continued (Ma yunzhi)
Re: None Re: Message no. 197 continued (Huaimin Liu)
Date: 12 Jan, 2006
From: Ma yunzhi <Ma yunzhi>

Dear Huaiming and others,

     I followed the instruction of Message NO. 183. It helped to fix the 
bug of "eDep > KineE", but it is likely not to fix the bug of 
CSDA range.

     In the output of my new simulation, the TrackLength of each event
doesn't equal to each other, and they do not equal the right csda range
143 mum:

           1.55 mm, 0.537 mm, and 3.71 mm

     What does this mean? Under conditions contained 
in my macro file, the TrackLength of the event SHOULD be the very 
scda range. Is this right?

     Thank you for your prompt reply.


INFO:

My new low energy PhisList class looks like,
////////////////////////////////////////////////////////////////
    } else if (particleName == "e-") {
      //electron
      pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
      pmanager->AddProcess(new G4LowEnergyIonisation,        -1, 2,2);
      pmanager->AddProcess(new G4LowEnergyBremsstrahlung,-1,-1,3); //NEW
      //pmanager->AddProcess(new G4LowEnergyBremsstrahlung,-1,3,3); //OLD

    } else if (particleName == "e+") {
////////////////////////////////////////////////////////////////

and my micro-file contains,
////////////////////////////////////////////////////////////////
#
# Eloss only; no fluctuation; e- 1 GeV
#
/control/verbose 2
#
/testem/phys/addPhysics  lowen
#
/testem/phys/setGCut  95 cm
/testem/phys/setECut  95 cm
#
/run/initialize
#
/process/inactivate phot
/process/inactivate compt
/process/inactivate conv

/process/eLoss/integral true
/process/eLoss/fluct 0
#
/run/verbose 2
#
/gun/particle e-
/tracking/verbose 2
/gun/energy 100 keV
/run/beamOn 3
////////////////////////////////////////////////////////////////

BUT the CSDA ranges given by 3 events don't equal to each other.
look,
///////////////////////////////////////////////////////////////
---> 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    63.7 mum  -123 mum  8.45 keV  91.6 keV   537 mum   537 mum       Water  LowEnergyIoni
    2     -5 m    97.7 mum  -146 mum  58.3 eV   8.39 keV   506 mum  1.04 mm        Water  LowEnergyIoni
    3     -5 m    97.7 mum  -146 mum     0 eV   58.3 eV    505 mum  1.55 mm        Water  LowEnergyIoni
                                                                    OOOOOOOO

*********************************************************************************************************
* 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   537 mum   537 mum       Water  LowEnergyIoni
                                                                     OOOOOOOO

*********************************************************************************************************
* 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    -120 mum  69.7 mum  18.1 keV  81.9 keV   537 mum   537 mum       Water  LowEnergyIoni
    2     -5 m    -111 mum -12.5 mum  2.01 keV  16.1 keV   507 mum  1.04 mm        Water  LowEnergyIoni
    3     -5 m    -115 mum -10.9 mum  82.1 eV   1.93 keV   505 mum  1.55 mm        Water  LowEnergyIoni
    4     -5 m    -114 mum -11.6 mum  47.5 eV   34.6 eV    505 mum  2.05 mm        Water  LowEnergyIoni
    5     -5 m    -114 mum -10.4 mum  32.3 eV   15.1 eV    505 mum  2.56 mm        Water  LowEnergyIoni
    6     -5 m    -113 mum -10.1 mum  31.3 eV   1.03 eV    505 mum  3.06 mm        Water  LowEnergyIoni
    7     -5 m    -113 mum -11.1 mum 0.382 eV   30.9 eV    505 mum  3.57 mm        Water  LowEnergyIoni
    8     -5 m    -113 mum -11.1 mum     0 eV  0.382 eV    139 mum  3.71 mm        Water  LowEnergyIoni
                                                                    OOOOOOO
Run terminated.
Run Summary
  Number of events processed : 3
///////////////////////////////////////////////////////////////

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

1 None: Re: BUT instruction of message 183 doesn't give the right SCDA range   (Vladimir IVANTCHENKO - 15 Jan, 2006)
 Add Message Add Message
to: "BUT instruction of message 183 doesn't give the right SCDA 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 ]