Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HadronPhysicsQGSP_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: HadronPhysicsQGSP_BIC
31 //
32 // Author: 2002 J.P. Wellisch
33 //
34 // Modified:
35 // 23.11.2005 G.Folger: migration to non static particles
36 // 08.06.2006 V.Ivanchenko: remove stopping
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.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);
60 
61 HadronPhysicsQGSP_BIC::HadronPhysicsQGSP_BIC(G4int)
62  : G4VPhysicsConstructor("hInelastic QGSP_BIC")
63  , theNeutrons(0)
64  , theLEPNeutron(0)
65  , theQGSPNeutron(0)
66  , theBinaryNeutron(0)
67  , thePiK(0)
68  , theLEPPiK(0)
69  , theQGSPPiK(0)
70  , thePro(0)
71  , theLEPPro(0)
72  , theQGSPPro(0)
73  , theBinaryPro(0)
74  , theMisc(0)
75  , QuasiElastic(true)
76 {}
77 
78 HadronPhysicsQGSP_BIC::HadronPhysicsQGSP_BIC(const G4String& name, G4bool quasiElastic)
79  : G4VPhysicsConstructor(name)
80  , theNeutrons(0)
81  , theLEPNeutron(0)
82  , theQGSPNeutron(0)
83  , theBinaryNeutron(0)
84  , thePiK(0)
85  , theLEPPiK(0)
86  , theQGSPPiK(0)
87  , thePro(0)
88  , theLEPPro(0)
89  , theQGSPPro(0)
90  , theBinaryPro(0)
91  , theMisc(0)
92  , QuasiElastic(quasiElastic)
93 {}
94 
95 void HadronPhysicsQGSP_BIC::CreateModels()
96 {
97  theNeutrons=new G4NeutronBuilder;
98 
99  theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic));
100  theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
101  theLEPNeutron->SetMinInelasticEnergy(9.5*GeV);
102  theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
103 
104  theNeutrons->RegisterMe(theBinaryNeutron=new G4BinaryNeutronBuilder);
105  theBinaryNeutron->SetMaxEnergy(9.9*GeV);
106 
107  thePro=new G4ProtonBuilder;
108  thePro->RegisterMe(theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic));
109  thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
110  theLEPPro->SetMinEnergy(9.5*GeV);
111  theLEPPro->SetMaxEnergy(25*GeV);
112 
113  thePro->RegisterMe(theBinaryPro=new G4BinaryProtonBuilder);
114  theBinaryPro->SetMaxEnergy(9.9*GeV);
115 
116  thePiK=new G4PiKBuilder;
117  thePiK->RegisterMe(theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic));
118  thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
119  theLEPPiK->SetMaxEnergy(25*GeV);
120 
121  theMisc=new G4MiscBuilder;
122 }
123 
124 HadronPhysicsQGSP_BIC::~HadronPhysicsQGSP_BIC()
125 {
126  delete theMisc;
127  delete theQGSPNeutron;
128  delete theLEPNeutron;
129  delete theBinaryNeutron;
130  delete theQGSPPro;
131  delete theLEPPro;
132  delete thePro;
133  delete theBinaryPro;
134  delete theQGSPPiK;
135  delete theLEPPiK;
136  delete thePiK;
137 }
138 
139 void HadronPhysicsQGSP_BIC::ConstructParticle()
140 {
141  G4MesonConstructor pMesonConstructor;
142  pMesonConstructor.ConstructParticle();
143 
144  G4BaryonConstructor pBaryonConstructor;
145  pBaryonConstructor.ConstructParticle();
146 
147  G4ShortLivedConstructor pShortLivedConstructor;
148  pShortLivedConstructor.ConstructParticle();
149 }
150 
151 #include "G4ProcessManager.hh"
152 void HadronPhysicsQGSP_BIC::ConstructProcess()
153 {
154  CreateModels();
155  theNeutrons->Build();
156  thePro->Build();
157  thePiK->Build();
158  theMisc->Build();
159 }
160 
Generated on Sat May 25 2013 14:33:28 for Geant4 by   doxygen 1.8.4