| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file BiasingOperation.cc
27 /// \brief Implementation of the BiasingOperat
28 //
29 //
30
31 //....oooOO0OOooo........oooOO0OOooo........oo
32 //....oooOO0OOooo........oooOO0OOooo........oo
33
34 #include "BiasingOperation.hh"
35
36 #include "G4BGGNucleonInelasticXS.hh"
37 #include "G4BGGPionInelasticXS.hh"
38 #include "G4BiasingProcessInterface.hh"
39 #include "G4CascadeInterface.hh"
40 #include "G4CrossSectionDataStore.hh"
41 #include "G4ExcitedStringDecay.hh"
42 #include "G4FTFModel.hh"
43 #include "G4GeneratorPrecompoundInterface.hh"
44 #include "G4HadronInelasticProcess.hh"
45 #include "G4HadronicParameters.hh"
46 #include "G4INCLXXInterface.hh"
47 #include "G4LundStringFragmentation.hh"
48 #include "G4NeutronInelasticXS.hh"
49 #include "G4TheoFSGenerator.hh"
50 #include "G4VParticleChange.hh"
51
52 //....oooOO0OOooo........oooOO0OOooo........oo
53
54 BiasingOperation::BiasingOperation(G4String na
55 {
56 // Create the inelastic processes for p , n
57 fProtonInelasticProcess = new G4HadronInelas
58 fNeutronInelasticProcess =
59 new G4HadronInelasticProcess("neutronInela
60 fPionPlusInelasticProcess =
61 new G4HadronInelasticProcess("pi+Inelastic
62 fPionMinusInelasticProcess =
63 new G4HadronInelasticProcess("pi-Inelastic
64
65 // Set the energy ranges
66 const G4double maxBERT = 41.0 * CLHEP::MeV;
67 const G4double maxINCLXX = 12.0 * CLHEP::GeV
68 const G4double minINCLXX = 40.0 * CLHEP::MeV
69 const G4double minFTFP = 3.0 * CLHEP::GeV;
70 const G4double maxFTFP = G4HadronicParameter
71
72 // Create the hadronic models (to replace FT
73 // keeping the same energy ranges for the tr
74 // Notice that it is better to create the mo
75 // instead of reusing the existing ones, bec
76 // existing ones associated to the wrong par
77 // --- FTFP model ---
78 G4FTFModel* theStringModel = new G4FTFModel;
79 G4LundStringFragmentation* theLund = new G4L
80 G4ExcitedStringDecay* theStringDecay = new G
81 theStringModel->SetFragmentationModel(theStr
82 G4GeneratorPrecompoundInterface* thePrecoInt
83 G4TheoFSGenerator* theHighEnergyModel = new
84 theHighEnergyModel->SetHighEnergyGenerator(t
85 theHighEnergyModel->SetTransport(thePrecoInt
86 theHighEnergyModel->SetMinEnergy(minFTFP);
87 theHighEnergyModel->SetMaxEnergy(maxFTFP);
88 // Bertini : create a new model to be used b
89 G4CascadeInterface* theBertiniModel = new G4
90 theBertiniModel->SetMinEnergy(0.0);
91 theBertiniModel->SetMaxEnergy(maxBERT);
92 // --- Preco ---
93 // --- INCLXX model ---
94 // The instance for nucleons:
95 G4INCLXXInterface* theInclxxModel = new G4IN
96 theInclxxModel->SetMinEnergy(minINCLXX);
97 theInclxxModel->SetMaxEnergy(maxINCLXX); //
98
99 // Register the models
100 fProtonInelasticProcess->RegisterMe(theHighE
101 fProtonInelasticProcess->RegisterMe(theInclx
102 fProtonInelasticProcess->RegisterMe(theBerti
103 fNeutronInelasticProcess->RegisterMe(theHigh
104 fNeutronInelasticProcess->RegisterMe(theIncl
105 fNeutronInelasticProcess->RegisterMe(theBert
106 fPionPlusInelasticProcess->RegisterMe(theHig
107 fPionPlusInelasticProcess->RegisterMe(theInc
108 fPionPlusInelasticProcess->RegisterMe(theBer
109 fPionMinusInelasticProcess->RegisterMe(theHi
110 fPionMinusInelasticProcess->RegisterMe(theIn
111 fPionMinusInelasticProcess->RegisterMe(theBe
112
113 G4VCrossSectionDataSet* theProtonXSdata = ne
114 theProtonXSdata->BuildPhysicsTable(*(G4Proto
115 fProtonInelasticProcess->AddDataSet(theProto
116 G4VCrossSectionDataSet* theNeutronXSdata = n
117 theNeutronXSdata->BuildPhysicsTable(*(G4Neut
118 fNeutronInelasticProcess->AddDataSet(theNeut
119 G4VCrossSectionDataSet* thePionPlusXSdata =
120 thePionPlusXSdata->BuildPhysicsTable(*(G4Pio
121 fPionPlusInelasticProcess->AddDataSet(thePio
122 G4VCrossSectionDataSet* thePionMinusXSdata =
123 thePionMinusXSdata->BuildPhysicsTable(*(G4Pi
124 fPionMinusInelasticProcess->AddDataSet(thePi
125 }
126
127 //....oooOO0OOooo........oooOO0OOooo........oo
128
129 BiasingOperation::~BiasingOperation() {}
130
131 //....oooOO0OOooo........oooOO0OOooo........oo
132
133 G4VParticleChange* BiasingOperation::ApplyFina
134
135
136 {
137 if (track->GetParticleDefinition() == G4Prot
138 auto particle = track->GetDynamicParticle(
139 auto material = track->GetMaterial();
140 fProtonInelasticProcess->GetCrossSectionDa
141 return fProtonInelasticProcess->PostStepDo
142 }
143 else if (track->GetParticleDefinition() == G
144 auto particle = track->GetDynamicParticle(
145 auto material = track->GetMaterial();
146 fNeutronInelasticProcess->GetCrossSectionD
147 return fNeutronInelasticProcess->PostStepD
148 }
149 else if (track->GetParticleDefinition() == G
150 auto particle = track->GetDynamicParticle(
151 auto material = track->GetMaterial();
152 fPionPlusInelasticProcess->GetCrossSection
153 return fPionPlusInelasticProcess->PostStep
154 }
155 else if (track->GetParticleDefinition() == G
156 auto particle = track->GetDynamicParticle(
157 auto material = track->GetMaterial();
158 fPionMinusInelasticProcess->GetCrossSectio
159 return fPionMinusInelasticProcess->PostSte
160 }
161 else {
162 G4cerr << "ERROR in BiasingOperation::Appl
163 << track->GetParticleDefinition()->
164 return 0;
165 }
166 }
167
168 //....oooOO0OOooo........oooOO0OOooo........oo
169