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