Problem 1321 - Apparent Energy violation in ElectoNuclear scattering
Summary: Apparent Energy violation in ElectoNuclear scattering
Status: RESOLVED WONTFIX
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/chiral_inv_phase_space (show other problems)
Version: 9.5
Hardware: PC Linux
: P4 normal
Assignee: Alberto.Ribon
URL:
Depends on:
Blocks:
 
Reported: 2012-06-01 19:35 CEST by John Apostolakis
Modified: 2012-10-23 11:25 CEST (History)
2 users (show)

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description John Apostolakis 2012-06-01 19:35:29 CEST
When checking for energy violations in hadronic processes is activated, unit tests reveal that electro-nuclear scattering does not conserve energy.

Approximately 1 in 100 interactions for 100 GeV electrons in air has a problem.  ( 13 out of 1000 in one sample.)

The output shows apparent imbalances from 17 MeV to 17 GeV:

 Process: ElectroNuclear , Model: CHIPSElectroNuclear
 Primary: e- (11), E= 100000, target nucleus (7,14)
 relative limit 0.01 relative value = 0.00282248 pass
 absolute limit (MeV) 1 absolute value (MeV) = 282.247 fail

 Process: ElectroNuclear , Model: CHIPSElectroNuclear
 Primary: e- (11), E= 100000, target nucleus (8,16)
 relative limit 0.01 relative value = 0.250539 fail
 absolute limit (MeV) 1 absolute value (MeV) = 25053.9 fail

 Process: ElectroNuclear , Model: CHIPSElectroNuclear
 Primary: e- (11), E= 100000, target nucleus (7,14)
 relative limit 0.01 relative value = 0.000283411 pass
 absolute limit (MeV) 1 absolute value (MeV) = 28.341 fail

 Process: ElectroNuclear , Model: CHIPSElectroNuclear
 Primary: e- (11), E= 100000, target nucleus (7,14)
 relative limit 0.01 relative value = 0.000176907 pass
 absolute limit (MeV) 1 absolute value (MeV) = 17.6907 fail

To repeat this the following environment variables are needed:
   setenv G4Hadronic_epReportLevel 3
   setenv G4Hadronic_epCheckRelativeLevel 0.01
   setenv G4Hadronic_epCheckAbsoluteLevel 1.00
Comment 1 Alberto.Ribon 2012-10-23 11:25:31 CEST
The CHIPS package is no longer supported and will be dropped in the Geant4 version 10 .