Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HadronPhysicsQGSP_BIC_HP.cc
Go to the documentation of this file.
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 // $Id$
27 //
28 //---------------------------------------------------------------------------
29 //
30 // ClassName: HadronPhysicsQGSP_BIC_HP
31 //
32 // Author: 2006 G.Folger
33 //
34 // Based on HadronPhysicsQGSP_BIC
35 //
36 // Modified:
37 // 25.04.2007 G.Folger: Add code for quasielastic
38 // 31.10.2012 A.Ribon: Use G4MiscBuilder
39 //
40 //----------------------------------------------------------------------------
41 //
42 #include <iomanip>
43 
44 #include "HadronPhysicsQGSP_BIC_HP.hh"
45 
46 #include "globals.hh"
47 #include "G4ios.hh"
48 #include "G4SystemOfUnits.hh"
49 #include "G4ParticleDefinition.hh"
50 #include "G4ParticleTable.hh"
51 
52 #include "G4MesonConstructor.hh"
53 #include "G4BaryonConstructor.hh"
54 #include "G4ShortLivedConstructor.hh"
55 
56 // factory
57 #include "G4PhysicsConstructorFactory.hh"
58 //
59 G4_DECLARE_PHYSCONSTR_FACTORY(HadronPhysicsQGSP_BIC_HP);
60 
61 HadronPhysicsQGSP_BIC_HP::HadronPhysicsQGSP_BIC_HP(G4int)
62  : G4VPhysicsConstructor("hInelastic QGSP_BIC_HP")
63  , theNeutrons(0)
64  , theLEPNeutron(0)
65  , theQGSPNeutron(0)
66  , theBinaryNeutron(0)
67  , theHPNeutron(0)
68  , thePiK(0)
69  , theLEPPiK(0)
70  , theQGSPPiK(0)
71  , thePro(0)
72  , theLEPPro(0)
73  , theQGSPPro(0)
74  , theBinaryPro(0)
75  , theMisc(0)
76  , QuasiElastic(true)
77 {}
78 
79 HadronPhysicsQGSP_BIC_HP::HadronPhysicsQGSP_BIC_HP(const G4String& name, G4bool quasiElastic)
80  : G4VPhysicsConstructor(name)
81  , theNeutrons(0)
82  , theLEPNeutron(0)
83  , theQGSPNeutron(0)
84  , theBinaryNeutron(0)
85  , theHPNeutron(0)
86  , thePiK(0)
87  , theLEPPiK(0)
88  , theQGSPPiK(0)
89  , thePro(0)
90  , theLEPPro(0)
91  , theQGSPPro(0)
92  , theBinaryPro(0)
93  , theMisc(0)
94  , QuasiElastic(quasiElastic)
95 {}
96 
97 void HadronPhysicsQGSP_BIC_HP::CreateModels()
98 {
99  theNeutrons=new G4NeutronBuilder;
100 
101  theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic));
102  theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
103  theLEPNeutron->SetMinEnergy(19.9*MeV);
104  theLEPNeutron->SetMinInelasticEnergy(9.5*GeV);
105  theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
106 
107  theNeutrons->RegisterMe(theBinaryNeutron=new G4BinaryNeutronBuilder);
108  theBinaryNeutron->SetMinEnergy(19.9*MeV);
109  theBinaryNeutron->SetMaxEnergy(9.9*GeV);
110 
111  theNeutrons->RegisterMe(theHPNeutron=new G4NeutronHPBuilder);
112 
113  thePro=new G4ProtonBuilder;
114  thePro->RegisterMe(theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic));
115  thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
116  theLEPPro->SetMinEnergy(9.5*GeV);
117  theLEPPro->SetMaxEnergy(25*GeV);
118 
119  thePro->RegisterMe(theBinaryPro=new G4BinaryProtonBuilder);
120  theBinaryPro->SetMaxEnergy(9.9*GeV);
121 
122  thePiK=new G4PiKBuilder;
123  thePiK->RegisterMe(theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic));
124  thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
125  theLEPPiK->SetMaxEnergy(25*GeV);
126 
127  theMisc=new G4MiscBuilder;
128 }
129 
130 HadronPhysicsQGSP_BIC_HP::~HadronPhysicsQGSP_BIC_HP()
131 {
132  delete theMisc;
133  delete theQGSPNeutron;
134  delete theLEPNeutron;
135  delete theBinaryNeutron;
136  delete theHPNeutron;
137  delete theQGSPPro;
138  delete theLEPPro;
139  delete thePro;
140  delete theBinaryPro;
141  delete theQGSPPiK;
142  delete theLEPPiK;
143  delete thePiK;
144 }
145 
146 void HadronPhysicsQGSP_BIC_HP::ConstructParticle()
147 {
148  G4MesonConstructor pMesonConstructor;
149  pMesonConstructor.ConstructParticle();
150 
151  G4BaryonConstructor pBaryonConstructor;
152  pBaryonConstructor.ConstructParticle();
153 
154  G4ShortLivedConstructor pShortLivedConstructor;
155  pShortLivedConstructor.ConstructParticle();
156 }
157 
158 #include "G4ProcessManager.hh"
159 void HadronPhysicsQGSP_BIC_HP::ConstructProcess()
160 {
161  CreateModels();
162  theNeutrons->Build();
163  thePro->Build();
164  thePiK->Build();
165  theMisc->Build();
166 }
167 
Generated on Sat May 25 2013 14:33:28 for Geant4 by   doxygen 1.8.4