Problem 2445

Summary: G4NeutronCaptureXS isotopic cross-section problem
Product: Geant4 Reporter: Zoltan <zepofi>
Component: processes/hadronic/cross_sectionsAssignee: Vladimir.Ivantchenko
Status: RESOLVED FIXED    
Severity: minor CC: Alberto.Ribon
Priority: P4    
Version: 10.7   
Hardware: All   
OS: All   

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