Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
HadronPhysicsQGSC_BERT.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: HadronPhysicsQGSC_BERT
31 //
32 // Author: 2007 G.Folger
33 // created from HadronPhysicsQGSC, created by J.P. Wellisch
34 //
35 // Modified:
36 //
37 //----------------------------------------------------------------------------
38 //
39 #include <iomanip>
40 
41 #include "HadronPhysicsQGSC_BERT.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"
51 #include "G4ShortLivedConstructor.hh"
52 
53 // factory
54 #include "G4PhysicsConstructorFactory.hh"
55 //
56 G4_DECLARE_PHYSCONSTR_FACTORY(HadronPhysicsQGSC_BERT);
57 
58 HadronPhysicsQGSC_BERT::HadronPhysicsQGSC_BERT(G4int)
59  : G4VPhysicsConstructor("hInelastic QGSC_BERT")
60  , theNeutrons(0)
61  , theQGSCNeutron(0)
62  , theBertiniNeutron(0)
63  , theLEPNeutron(0)
64  , thePiK(0)
65  , theQGSCPiK(0)
66  , theBertiniPiK(0)
67  , thePro(0)
68  , theQGSCPro(0)
69  , theBertiniPro(0)
70  , theMiscLHEP(0)
71  , QuasiElastic(true)
72 {}
73 
74 HadronPhysicsQGSC_BERT::HadronPhysicsQGSC_BERT(const G4String& name, G4bool quasiElastic)
75  : G4VPhysicsConstructor(name)
76  , theNeutrons(0)
77  , theQGSCNeutron(0)
78  , theBertiniNeutron(0)
79  , theLEPNeutron(0)
80  , thePiK(0)
81  , theQGSCPiK(0)
82  , theBertiniPiK(0)
83  , thePro(0)
84  , theQGSCPro(0)
85  , theBertiniPro(0)
86  , theMiscLHEP(0)
87  , QuasiElastic(quasiElastic)
88 {}
89 
90 void HadronPhysicsQGSC_BERT::CreateModels()
91 {
92  theNeutrons=new G4NeutronBuilder;
93  theNeutrons->RegisterMe(theQGSCNeutron=new G4QGSCNeutronBuilder(QuasiElastic));
94  theNeutrons->RegisterMe(theBertiniNeutron=new G4BertiniNeutronBuilder);
95 
96  theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
97  theLEPNeutron->SetMinInelasticEnergy(0.0*eV); // no inelastic from LEP
98  theLEPNeutron->SetMaxInelasticEnergy(0.0*eV);
99 
100  theQGSCNeutron->SetMinEnergy(6*GeV);
101  theBertiniNeutron->SetMinEnergy(0.0*GeV);
102  theBertiniNeutron->SetMaxEnergy(9.0*GeV);
103 
104  thePro=new G4ProtonBuilder;
105  thePro->RegisterMe(theQGSCPro=new G4QGSCProtonBuilder(QuasiElastic));
106  thePro->RegisterMe(theBertiniPro=new G4BertiniProtonBuilder);
107 
108  theQGSCPro->SetMinEnergy(6*GeV);
109  theBertiniPro->SetMaxEnergy(9.0*GeV);
110 
111  thePiK=new G4PiKBuilder;
112  thePiK->RegisterMe(theQGSCPiK=new G4QGSCPiKBuilder(QuasiElastic));
113  thePiK->RegisterMe(theBertiniPiK=new G4BertiniPiKBuilder);
114 
115  theQGSCPiK->SetMinEnergy(6*GeV);
116  theBertiniPiK->SetMaxEnergy(9.0*GeV);
117 
118  theMiscLHEP=new G4MiscLHEPBuilder;
119 }
120 
121 HadronPhysicsQGSC_BERT::~HadronPhysicsQGSC_BERT()
122 {
123  delete theQGSCNeutron;
124  delete theBertiniNeutron;
125  delete theLEPNeutron;
126  delete theNeutrons;
127  delete theQGSCPro;
128  delete theBertiniPro;
129  delete thePro;
130  delete theQGSCPiK;
131  delete theBertiniPiK;
132  delete thePiK;
133  delete theMiscLHEP;
134 }
135 
136 void HadronPhysicsQGSC_BERT::ConstructParticle()
137 {
138  G4MesonConstructor pMesonConstructor;
139  pMesonConstructor.ConstructParticle();
140 
141  G4BaryonConstructor pBaryonConstructor;
142  pBaryonConstructor.ConstructParticle();
143 
144  G4ShortLivedConstructor pShortLivedConstructor;
145  pShortLivedConstructor.ConstructParticle();
146 }
147 
148 #include "G4ProcessManager.hh"
149 void HadronPhysicsQGSC_BERT::ConstructProcess()
150 {
151  CreateModels();
152  theNeutrons->Build();
153  thePro->Build();
154  thePiK->Build();
155  theMiscLHEP->Build();
156 }
157 
Generated on Sat May 25 2013 14:33:28 for Geant4 by   doxygen 1.8.4