Message: No 2.22 MeV gamma after n capture on proton Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question No 2.22 MeV gamma after n capture on proton 

Forum: Hadronic Processes
Date: 15 May, 2012
From: Yuri Shitov <Yuri Shitov>

Hi to all,

I simulate low neutrons in large plastic scintillator. I've discovered, that with new geant4 version:

*************************************************************
 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
*************************************************************

<<< Geant4 Physics List simulation engine: QGSP_BIC_HP 1.3

No 2.22 MeV gammas are generated after capture in proton, deuteron is here.

---------------------------------------------------

# original track #1

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0      446    -26.5      452      0.02        0        0         0 scinstripcore.phys initStep
    1      451    -31.7      443   0.00307        0     11.1      11.1 scinstripcore.phys hadElastic
...
  222      404     -175      386         0        0     4.97       842 scinstripcore.phys nCapture

* G4Track Information: Particle = deuteron, Track ID = 5, Parent ID = 1

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0      404     -175      386    0.0013        0        0         0 scinstripcore.phys initStep
    1      404     -175      386         0   0.0013 0.000209  0.000209 scinstripcore.phys hIoni
...
# secondaries tracks connected to primary track #0
* G4Track Information:   Particle = deuteron,   Track ID = 4,   Parent ID = 1
...
* G4Track Information:   Particle = proton,   Track ID = 3,   Parent ID = 1
...
* G4Track Information:   Particle = proton,   Track ID = 2,   Parent ID = 1
...
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
...
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0

# next primary

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

---------------------------------------------------

The SAME code compiled with previously used version of geant works well:

---------------------------------------------------

*************************************************************
 Geant4 version Name: geant4-09-02-patch-02    (28-August-2009)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

<<< Geant4 Physics List engine packaging library: PACK 5.4
<<< Geant4 Physics List simulation engine: QGSP_BIC_HP

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
  312     68.1     -128     -191  1.08e-08        0        0  1.78e+03 scinstripcore.phys initStep
  313     75.4     -127     -194         0        0     7.84  1.79e+03 scinstripcore.phys nCapture

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 288,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     75.4     -127     -194      2.22        0        0         0 scinstripcore.phys initStep
    1     64.7     -100     -167      2.22        0     39.4      39.4 scinstriprefl.phys Transportation

---------------------------------------------------

Any ideas how one can restore normal behavior?

Cheers, Yuri

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

1 None: Re: No 2.22 MeV gamma after n capture on proton   (Yuri Shitov - 15 May, 2012)
(_ None: Re: No 2.22 MeV gamma after n capture on proton   (Koi, Tatsumi - 15 May, 2012)
(_ Ok: Re: No 2.22 MeV gamma after n capture on proton   (Yuri Shitov - 16 May, 2012)
 Add Message Add Message
to: "No 2.22 MeV gamma after n capture on proton"

 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 ]