Problem 1126 - Unphysical secondary particles in 2 GeV gamma-nuclear interaction in lead
Summary: Unphysical secondary particles in 2 GeV gamma-nuclear interaction in lead
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/chiral_inv_phase_space (show other problems)
Version: 9.2
Hardware: All All
: P5 trivial
Assignee: dennis.herbert.wright
URL: https://savannah.cern.ch/bugs/?68785
Depends on:
Blocks:
 
Reported: 2010-07-08 17:15 CEST by John Apostolakis
Modified: 2010-09-16 11:02 CEST (History)
0 users

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description John Apostolakis 2010-07-08 17:15:45 CEST
ATLAS has reported anomalous event which has high pT jets not matching any MC truth particle. 

2GeV gamma produced 90GeV pi0 and 95GeV lambda !!

gamma e=1944.0247 -> 
Process: PhotonInelastic at Eta=-1.0890 Phi=0.9310 R=1504.9776 

Resulting particles:
       * pi0        e=90219.733   px=-2821.7829, py=-66930.4036, pz=60632.6144 
       * lambda e=95020.0396  px=5196.8159, py=70747.6641, pz=-64873.6737 
       *  ... 

This information is extracted from the ATLAS Savannah report (URL above).

The material of in the volume of the PhotonInelastic interaction is 

material: Thinabs 
mass fraction of elements are 
Pb 0.7493 
Fe 0.1837 
C 0.0251 
H 0.0034 
O 0.0268 
Si 0.0117
Comment 1 John Apostolakis 2010-07-08 17:23:56 CEST
(From Andrea Dotti's report, 16 June 2010):

The problem has been investigated, reproduced in a simple setup.

The bug is confirmed to occur in the reaction: 
Gamma + A -> X 
(Energy of the gamma of 2 GeV), 
for which at least one of the secondaries has too much kinetic energy (up to 90 GeV).   Typically 

The bug has been confirmed to affect A=Pb (but other materials may be effected). 
We have traced down the problem to the relevant parts of the code and we are investigating it. 
The nature of the bug makes it difficult to measure the frequency of the problematic events, it can be estimated to occur once every 10^6 interactions. 

A candidate fix has been identified and undergone testing.
Comment 2 Gabriele Cosmo 2010-09-16 11:02:55 CEST
Fixed. The fix, which relates only 9.2 series release, has been
provided to ATLAS and will be included in release 9.2.p04.