Message: Re: PiMinusAbsorptionBertini multiple calls Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: PiMinusAbsorptionBertini multiple calls 

Forum: Hadronic Processes
Re: Question PiMinusAbsorptionBertini multiple calls (Sam)
Re: None Re: PiMinusAbsorptionBertini multiple calls (Michael H. Kelsey)
Date: 11 Jun, 2012
From: Sam <Sam>

I am using 9.5.p01. Here is the output for /tracking/verbose 2.

---> Begin of Event: 0

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0    873 um  -4.11 mm   -1.3 cm      1 MeV     0 eV      0 fm      0 fm     D2target    initStep
    1    996 um  -4.22 mm  -1.32 cm    668 keV   332 keV   240 um    240 um     D2target       hIoni
    2   1.08 mm  -4.29 mm  -1.33 cm    281 keV   387 keV   172 um    412 um     D2target       hIoni
    3   1.11 mm  -4.31 mm  -1.33 cm      0 eV    281 keV  51.8 um    463 um     D2target       hIoni
    4   1.11 mm  -4.31 mm  -1.33 cm      0 eV      0 eV      0 fm    463 um     D2target  PiMinusAbsorptionBertini
    :----- List of 2ndaries - #SpawnInStep=  3(Rest= 3,Along= 0,Post= 0), #SpawnTotal=  3 ---------------
    :   1.11 mm  -4.31 mm  -1.33 cm   29.8 MeVanti_nu_mu
    :   1.11 mm  -4.31 mm  -1.33 cm   4.12 MeV       mu-
    :   1.11 mm  -4.31 mm  -1.33 cm      1 keV       pi-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = pi-,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0   1.11 mm  -4.31 mm  -1.33 cm      1 keV     0 eV      0 fm      0 fm     D2target    initStep
    1   1.11 mm  -4.31 mm  -1.33 cm      0 eV      1 keV   714 nm    714 nm     D2target       hIoni
    2   1.11 mm  -4.31 mm  -1.33 cm      0 eV      0 eV      0 fm    714 nm     D2target  PiMinusAbsorptionBertini
    :----- List of 2ndaries - #SpawnInStep=  3(Rest= 3,Along= 0,Post= 0), #SpawnTotal=  3 ---------------
    :   1.11 mm  -4.31 mm  -1.33 cm   29.8 MeVanti_nu_mu
    :   1.11 mm  -4.31 mm  -1.33 cm   4.12 MeV       mu-
    :   1.11 mm  -4.31 mm  -1.33 cm      1 keV       pi-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

I am only seeing the PiMinusAbsorptionBertini process used once the pi- comes to rest. The reactions I expect to happen are pi- + d -> 2n and pi- + p -> n g.

The processes loaded onto pi- are

G4ProcessManager:  particle[pi-]
[0]=== process[Transportation :Transportation] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:       2:       0:       7:       0:
  parameter             -1:      -1:       0:       0:       0:       0:
[1]=== process[msc :Electromagnetic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:       1:       1:       6:       1:
  parameter             -1:      -1:       1:       1:       1:       1:
[2]=== process[hIoni :Electromagnetic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:       0:       2:       5:       2:
  parameter             -1:      -1:       2:       2:       2:       2:
[3]=== process[hBrems :Electromagnetic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:      -1:      -1:       4:       3:
  parameter             -1:      -1:      -1:      -1:       3:       3:
[4]=== process[hPairProd :Electromagnetic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:      -1:      -1:       3:       4:
  parameter             -1:      -1:      -1:      -1:       4:       4:
[5]=== process[Decay :Decay] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                  1:       0:      -1:      -1:       2:       5:
  parameter           1000:    1000:      -1:      -1:    1000:    1000:
[6]=== process[hadElastic :Hadronic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:      -1:      -1:       1:       6:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:
[7]=== process[PionMinusInelastic :Hadronic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                 -1:      -1:      -1:      -1:       0:       7:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:
[8]=== process[PiMinusAbsorptionBertini :Hadronic] Active 
  Ordering::             AtRest             AlongStep          PostStep   
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt 
  Ordering::      
  index                  0:       1:      -1:      -1:      -1:      -1:
  parameter           1000:    1000:      -1:      -1:      -1:      -1:

Thank you,

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

1 None: Re: PiMinusAbsorptionBertini multiple calls   (Michael H. Kelsey - 11 Jun, 2012)
1 None: Re: PiMinusAbsorptionBertini multiple calls   (Sam - 11 Jun, 2012)
(_ None: Re: PiMinusAbsorptionBertini multiple calls   (Mike Kelsey - 11 Jun, 2012)
2 None: Re: PiMinusAbsorptionBertini multiple calls -- bug in 9.5-patch-01 version   (Michael H. Kelsey - 11 Jun, 2012)
1 Ok: Re: PiMinusAbsorptionBertini multiple calls -- bug in 9.5-patch-01 version   (Michael H. Kelsey - 12 Jun, 2012)
 Add Message Add Message
to: "Re: PiMinusAbsorptionBertini multiple calls"

 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 ]