Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HadronPhysicsQGS_BIC.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: HadronPhysicsQGS_BIC
31 //
32 // Author: 2007 Gunter Folger
33 // created from HadronPhysicsQGSP_BIC by H.P.Wellisch
34 //
35 // Modified:
36 //
37 //----------------------------------------------------------------------------
38 //
39 #include <iomanip>
40 
41 #include "HadronPhysicsQGS_BIC.hh"
42 
43 #include "globals.hh"
44 #include "G4ios.hh"
45 #include "G4SystemOfUnits.hh"
46 #include "G4ParticleDefinition.hh"
47 #include "G4ParticleTable.hh"
48 
49 #include "G4MesonConstructor.hh"
50 #include "G4BaryonConstructor.hh"
52 
53 // factory
55 //
57 
59  : G4VPhysicsConstructor("hInelastic QGS_BIC")
60  , theNeutrons(0)
61  , theLEPNeutron(0)
62  , theQGSBinaryNeutron(0)
63  , theBinaryNeutron(0)
64  , thePiK(0)
65  , theBICPiK(0)
66  , theLEPPiK(0)
67  , theQGSBinaryPiK(0)
68  , thePro(0)
69  , theLEPPro(0)
70  , theQGSBinaryPro(0)
71  , theBinaryPro(0)
72  , theMisc(0)
73  , QuasiElastic(true)
74 {}
75 
77  : G4VPhysicsConstructor(name)
78  , theNeutrons(0)
79  , theLEPNeutron(0)
80  , theQGSBinaryNeutron(0)
81  , theBinaryNeutron(0)
82  , thePiK(0)
83  , theBICPiK(0)
84  , theLEPPiK(0)
85  , theQGSBinaryPiK(0)
86  , thePro(0)
87  , theLEPPro(0)
88  , theQGSBinaryPro(0)
89  , theBinaryPro(0)
90  , theMisc(0)
91 , QuasiElastic(quasiElastic)
92 {}
93 
94 void HadronPhysicsQGS_BIC::CreateModels()
95 {
96  theNeutrons=new G4NeutronBuilder;
97 
98  theNeutrons->RegisterMe(theQGSBinaryNeutron=new G4QGSBinaryNeutronBuilder(QuasiElastic));
99  theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
100  theLEPNeutron->SetMinInelasticEnergy(9.5*GeV);
101  theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
102 
103  theNeutrons->RegisterMe(theBinaryNeutron=new G4BinaryNeutronBuilder);
104  theBinaryNeutron->SetMaxEnergy(9.9*GeV);
105 
106  thePro=new G4ProtonBuilder;
107  thePro->RegisterMe(theQGSBinaryPro=new G4QGSBinaryProtonBuilder(QuasiElastic));
108  thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
109  theLEPPro->SetMinEnergy(9.5*GeV);
110  theLEPPro->SetMaxEnergy(25*GeV);
111 
112  thePro->RegisterMe(theBinaryPro=new G4BinaryProtonBuilder);
113  theBinaryPro->SetMaxEnergy(9.9*GeV);
114 
115  thePiK=new G4PiKBuilder;
116  thePiK->RegisterMe(theQGSBinaryPiK=new G4QGSBinaryPiKBuilder(QuasiElastic));
117  thePiK->RegisterMe(theBICPiK = new G4BinaryPiKBuilder);
118  thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
119  theLEPPiK->SetMinPionEnergy(1.2*GeV);
120  theLEPPiK->SetMaxEnergy(25*GeV);
121 
122  theMisc=new G4MiscBuilder;
123 }
124 
126 {
127  delete theMisc;
128  delete theQGSBinaryNeutron;
129  delete theLEPNeutron;
130  delete theBinaryNeutron;
131  delete theQGSBinaryPro;
132  delete theLEPPro;
133  delete thePro;
134  delete theBinaryPro;
135  delete theQGSBinaryPiK;
136  delete theBICPiK;
137  delete theLEPPiK;
138  delete thePiK;
139 }
140 
142 {
143  G4MesonConstructor pMesonConstructor;
144  pMesonConstructor.ConstructParticle();
145 
146  G4BaryonConstructor pBaryonConstructor;
147  pBaryonConstructor.ConstructParticle();
148 
149  G4ShortLivedConstructor pShortLivedConstructor;
150  pShortLivedConstructor.ConstructParticle();
151 }
152 
153 #include "G4ProcessManager.hh"
155 {
156  CreateModels();
157  theNeutrons->Build();
158  thePro->Build();
159  thePiK->Build();
160  theMisc->Build();
161 }
162