I’m getting some confusing results. As a test I’ve tried creating a physics list using
GetReferencePhysList(“QGSP_BERT”);
and then replacing physics with
PhysicsList->ReplacePhysics(new G4HadronPhysicsQGSP_BERT())
which I exected to be identical, however I get different results when I analyse what is produced.
Is this expected?
In the print out for hadronic proceses I get duplicated sections like this:
Process: pi+Inelastic
Model: QGSP: 12 GeV —> 100 TeV
Model: FTFP: 9.5 GeV —> 25 GeV
Model: BertiniCascade: 0 eV —> 9.9 GeV
Cr_sctns: G4CrossSectionPairGG: 0 eV —> 100 TeV
G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections
below 91 GeV, Glauber-Gribov above
Cr_sctns: G4CrossSectionPairGG: 0 eV —> 100 TeV
G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections
below 91 GeV, Glauber-Gribov above
Cr_sctns: GheishaInelastic: 0 eV —> 100 TeV
Process: pi+Inelastic
Model: QGSP: 12 GeV —> 100 TeV
Model: FTFP: 9.5 GeV —> 25 GeV
Model: BertiniCascade: 0 eV —> 9.9 GeV
Cr_sctns: G4CrossSectionPairGG: 0 eV —> 100 TeV
G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections
below 91 GeV, Glauber-Gribov above
Cr_sctns: G4CrossSectionPairGG: 0 eV —> 100 TeV
G4CrossSectionPairGG: G4PiNuclearCrossSection cross sections
below 91 GeV, Glauber-Gribov above
Cr_sctns: GheishaInelastic: 0 eV —> 100 TeV