Problem 2445 - G4NeutronCaptureXS isotopic cross-section problem
Summary: G4NeutronCaptureXS isotopic cross-section problem
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/cross_sections (show other problems)
Version: 10.7
Hardware: All All
: P4 minor
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2021-11-17 13:48 CET by Zoltan
Modified: 2021-12-15 17:19 CET (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 Zoltan 2021-11-17 13:48:46 CET 
Hello,

I guess I found a bug in G4NeutronCaptureXS file.

In Initialise(), it is written:

    for(G4int A=amin[Z]; A<=amax[Z]; ++A) {
      std::ostringstream ost1;
      ost1 << gDataDirectory << Z << "_" << A;
      v = RetrieveVector(ost1, false);
      data->AddComponent(Z, A, v); 
    }

So, isotopic cross-sections are added with ID number equal to the mass.

However, in IsoCrossSection(), ID number = A-amin[Z] is used

    G4PhysicsVector* pviso = data->GetComponentDataByID(Z, A - amin[Z]);



Regards,
Zoltan
Comment 1 Vladimir.Ivantchenko 2021-12-15 17:19:10 CET 
Hello,

thank you for this report. The bug is already fixed in Geant4 10.7patch03. Please, switch to this patch.

VI