| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //---------------------------------------------------------------------------
28
29 #include <iomanip>
30 #include "G4HadronPhysicsQGSP_BIC_AllHP.hh"
31 #include "globals.hh"
32 #include "G4ios.hh"
33 #include "G4SystemOfUnits.hh"
34 #include "G4ParticleDefinition.hh"
35 #include "G4ParticleTable.hh"
36 #include "G4ProtonBuilder.hh"
37 #include "G4FTFPProtonBuilder.hh"
38 #include "G4QGSPProtonBuilder.hh"
39 #include "G4BinaryProtonBuilder.hh"
40 #include "G4ProtonPHPBuilder.hh"
41 #include "G4ProcessVector.hh"
42 #include "G4ProcessManager.hh"
43 #include "G4CrossSectionDataSetRegistry.hh"
44 #include "G4PhysListUtil.hh"
45 #include "G4HadronicParameters.hh"
46 // factory
47 #include "G4PhysicsConstructorFactory.hh"
48 //
49 G4_DECLARE_PHYSCONSTR_FACTORY( G4HadronPhysicsQGSP_BIC_AllHP );
50
51
52 G4HadronPhysicsQGSP_BIC_AllHP::G4HadronPhysicsQGSP_BIC_AllHP(G4int verb)
53 : G4HadronPhysicsQGSP_BIC_AllHP( "hInelastic QGSP_BIC_AllHP" )
54 {
55 G4HadronicParameters::Instance()->SetVerboseLevel(verb);
56 }
57
58 G4HadronPhysicsQGSP_BIC_AllHP::G4HadronPhysicsQGSP_BIC_AllHP( const G4String& name, G4bool quasiElastic )
59 : G4HadronPhysicsQGSP_BIC_HP( name, quasiElastic )
60 {
61 minBIC_proton = 190.0*CLHEP::MeV;
62 maxHP_proton = 200.0*CLHEP::MeV;
63 }
64
65
66 void G4HadronPhysicsQGSP_BIC_AllHP::Proton() {
67 G4HadronicParameters* param = G4HadronicParameters::Instance();
68 G4bool useFactorXS = param->ApplyFactorXS();
69
70 auto pro = new G4ProtonBuilder;
71 AddBuilder( pro );
72 auto qgs = new G4QGSPProtonBuilder( QuasiElasticQGS );
73 AddBuilder( qgs );
74 qgs->SetMinEnergy( minQGSP_proton );
75 pro->RegisterMe( qgs );
76 auto ftf = new G4FTFPProtonBuilder( QuasiElasticFTF );
77 AddBuilder( ftf );
78 ftf->SetMinEnergy( minFTFP_proton );
79 ftf->SetMaxEnergy( maxFTFP_proton );
80 pro->RegisterMe( ftf );
81 auto bic = new G4BinaryProtonBuilder;
82 AddBuilder( bic );
83 bic->SetMinEnergy( minBIC_proton );
84 bic->SetMaxEnergy( maxBIC_proton );
85 pro->RegisterMe( bic );
86 auto hp = new G4ProtonPHPBuilder;
87 AddBuilder( hp );
88 hp->SetMaxEnergy( maxHP_proton );
89 pro->RegisterMe( hp );
90 pro->Build();
91
92 const G4ParticleDefinition* proton = G4Proton::Proton();
93 G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(proton);
94 if(nullptr != inel) {
95 if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
96 }
97 }
98