Dear all
I have a question about how to get inelastic cross sections for antihelium-4. Taking Hadr09 example as a basis, I use anti-alpha particles as a projectile and protons as a target for [1.E-3, 1.E4] (GeV) energy range. In /src/HadronicGenerator.cc i ask for the inelastic cross sections as follows:
//Get cross section
static bool una_vez = true;
if(una_vez){
G4cout << "####################################" << G4endl;
//G4cout << "4Hebar + p" << G4endl;
G4cout << "4Hebar + Fe" << G4endl;
G4cout << "####################################" << G4endl;
G4cout << "projectileEnergy [GeV]" << " " << "Inelastic cross section [mbarn]" <<G4endl; //LFGC022025
una_vez = false;
}
G4DynamicParticle * dynParticle = new G4DynamicParticle( projectileDefinition, projectileDirection, projectileEnergy );
const G4Element* element = targetMaterial->GetElement(0);
//G4cout << "####################" << G4endl;
//G4cout << "Inelastic cross section = " <<theProcess->GetElementCrossSection(dynParticle, element, targetMaterial)*1e25 << " [mbarn]" << G4endl;
//G4cout << "####################" << G4endl;
G4cout << projectileEnergy/CLHEP::GeV << " " << theProcess->GetElementCrossSection(dynParticle, element, targetMaterial)*1e25 << G4endl; //LFGC 022025
return aChange;
}
Is this correct?, it works well for antideuteron and antihelium-3 cases, but not for 4He-bar because the values are much lower (see attached figure below). I hope you can help me
Thank you very much
Specifications:
OS: Ubuntu 24.04.2 LTS
G4: 11.2.0
cmake: 3.28.3
gcc: 13.3.0