Problem 937 - QGSP_BIC_HP physics list not using neutron HPFission, HPInelastic, or HPCapture
Summary: QGSP_BIC_HP physics list not using neutron HPFission, HPInelastic, or HPCapture
Status: RESOLVED FIXED
Alias: None
Product: Examples/Extended
Classification: Unclassified
Component: hadronic/Hadr00,01,02 (show other problems)
Version: 8.2
Hardware: All All
: P3 normal
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2007-04-17 23:19 CEST by Lauren Hsu
Modified: 2007-04-18 08:46 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 Lauren Hsu 2007-04-17 23:19:56 CEST 
This looks like a typo to me.  In the /Hadr01/src/PhysicsList.cc we have:

229 } else if (name == "QGSP_BIC_HP") {
230 
231     SetStandardList(true);
232     hadronPhys.push_back( new HadronPhysicsQGSP_BIC());
233     dump = true;

It seems to me that line 232 should be adding HadronPhysicsQGSP_BIC_HP instead of HadronPhysicsQGSP, which does not use neutron HP classes for fission, inelastic or capture processes.   

Meanwhile SetStandardList(true) ensures that HP is used for elastic scattering.  The way the example currently works, the user only gets HP elastic scattering, and all other neutron processes default to non-HP models.  This seems a little inconsistent to me (and not what I want for my application).    Is this what the author intended?
Comment 1 Vladimir.Ivantchenko 2007-04-18 08:46:41 CEST 
Thanks for pointing to this bug. It is fixed as you suggest and will be available in the next version of Geant4. For now you can fix it in your private version.

VI