Problem 753 - G4NeutronHPInelastic swaps D2A and T2A final states
Summary: G4NeutronHPInelastic swaps D2A and T2A final states
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models (show other problems)
Version: 7.0
Hardware: All All
: P2 normal
Assignee: Hans-Peter.Wellisch
URL:
Depends on:
Blocks:
 
Reported: 2005-06-12 22:49 CEST by naranjo
Modified: 2005-06-13 04:01 CEST (History)
1 user (show)

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description naranjo 2005-06-12 22:49:39 CEST 
In G4NeutronHPInelastic.cc,
	theInelastic[i].Register(&theD2AFS, "F32");
	theInelastic[i].Register(&theT2AFS, "F33");
should be:
	theInelastic[i].Register(&theT2AFS, "F32");
	theInelastic[i].Register(&theD2AFS, "F33");

Using G4NDL3.7, this problem only affects Boron-10, causing the production of a
D+2A+residual final state instead of the correct T+2A+residual final state.

Brian
Comment 1 Hans-Peter.Wellisch 2005-06-13 04:01:59 CEST 
Hi,

  thank you for notifying this. It turns out that this is not a problem in the
code, but an error in the data evaluation. The data formats for the two
channels are the same, so the code will run, even if the data have been located
incorrectly. If you are only interested in Boron, your fix will do the job very
well.

I'll take this into account in the new data evaluation I am preparing.

I also have to state that for reasons of copyright I have no longer the
possibility to distribute G4NDL publicly, and cannot be seen in support its
distribution as such...

Many greetings,

Hans-Peter.