Geant4  10.03
GammaRayTelIonPhysics.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 //
27 // $Id: GammaRayTelIonPhysics.cc 93678 2015-10-28 10:03:43Z gcosmo $
28 //
29 //
30 
31 #include <iomanip>
32 
33 #include "GammaRayTelIonPhysics.hh"
34 
35 #include "globals.hh"
36 #include "G4ios.hh"
37 #include "G4SystemOfUnits.hh"
38 
39 
41  : G4VPhysicsConstructor(name)
42 {;}
43 
45 {;}
46 
48 {;}
49 
50 
51 #include "G4ProcessManager.hh"
52 
53 
55 {
56  G4ProcessManager * pManager = 0;
57 
58  const G4double theBERTMin = 0.0*GeV;
59  const G4double theBERTMax = 5.0*GeV;
60  const G4double theFTFMin = 4.0*GeV;
61  const G4double theFTFMax = 100.0*TeV;
62 
63  G4FTFModel* theStringModel = new G4FTFModel;
65  theStringModel->SetFragmentationModel( theStringDecay );
66  G4PreCompoundModel* thePreEquilib = new G4PreCompoundModel( new G4ExcitationHandler );
67  G4GeneratorPrecompoundInterface* theCascade = new G4GeneratorPrecompoundInterface( thePreEquilib );
68 
69  G4TheoFSGenerator* theModel = new G4TheoFSGenerator( "FTFP" );
70  theModel->SetHighEnergyGenerator( theStringModel );
71  theModel->SetTransport( theCascade );
72  theModel->SetMinEnergy( theFTFMin );
73  theModel->SetMaxEnergy( theFTFMax );
74 
75  G4CascadeInterface * theBERTModel = new G4CascadeInterface;
76  theBERTModel->SetMinEnergy( theBERTMin );
77  theBERTModel->SetMaxEnergy( theBERTMax );
78 
79  // Elastic Process
82 
83  // Generic Ion
85  // add process
87 
88  pManager->AddProcess(&fIonIonisation, ordInActive, 2, 2);
89 
93 
94  G4ComponentGGNuclNuclXsc * ggNuclNuclXsec = new G4ComponentGGNuclNuclXsc();
95  G4VCrossSectionDataSet * theGGNuclNuclData = new G4CrossSectionInelastic(ggNuclNuclXsec);
96 
97  // Deuteron
99  // add process
101  fDeuteronProcess.AddDataSet(theGGNuclNuclData);
102  fDeuteronProcess.RegisterMe(theBERTModel);
103  fDeuteronProcess.RegisterMe(theModel);
105 
106  pManager->AddProcess(&fDeuteronIonisation, ordInActive, 2, 2);
107 
111 
112  // Triton
113  pManager = G4Triton::Triton()->GetProcessManager();
114  // add process
116  fTritonProcess.AddDataSet(theGGNuclNuclData);
117  fTritonProcess.RegisterMe(theBERTModel);
118  fTritonProcess.RegisterMe(theModel);
120 
121  pManager->AddProcess(&fTritonIonisation, ordInActive, 2, 2);
122 
126 
127  // Alpha
128  pManager = G4Alpha::Alpha()->GetProcessManager();
129  // add process
131 
132  fAlphaProcess.AddDataSet(theGGNuclNuclData);
133  fAlphaProcess.RegisterMe(theBERTModel);
134  fAlphaProcess.RegisterMe(theModel);
135  pManager->AddDiscreteProcess(&fAlphaProcess);
136 
137  pManager->AddProcess(&fAlphaIonisation, ordInActive, 2, 2);
138 
142 
143  // He3
144  pManager = G4He3::He3()->GetProcessManager();
145  // add process
147 
148  pManager->AddProcess(&fHe3Ionisation, ordInActive, 2, 2);
149 
153 
154 }
155 
156 
157 
G4hMultipleScattering fDeuteronMultipleScattering
void SetFragmentationModel(G4VStringFragmentation *aModel)
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4TritonInelasticProcess fTritonProcess
const char * name(G4int ptype)
void SetHighEnergyGenerator(G4VHighEnergyGenerator *const value)
void RegisterMe(G4HadronicInteraction *a)
static constexpr double TeV
Definition: G4SIunits.hh:218
G4HadronElastic * theElasticModel
void SetMinEnergy(G4double anEnergy)
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
G4DeuteronInelasticProcess fDeuteronProcess
GammaRayTelIonPhysics(const G4String &name="ion")
G4HadronElasticProcess theElasticProcess
G4hMultipleScattering fTritonMultipleScattering
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
static G4Triton * Triton()
Definition: G4Triton.cc:95
void SetProcessOrdering(G4VProcess *aProcess, G4ProcessVectorDoItIndex idDoIt, G4int ordDoIt=ordDefault)
G4AlphaInelasticProcess fAlphaProcess
G4hMultipleScattering fHe3MultipleScattering
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
static G4GenericIon * GenericIon()
Definition: G4GenericIon.cc:93
G4ProcessManager * GetProcessManager() const
static constexpr double GeV
Definition: G4SIunits.hh:217
void SetMaxEnergy(const G4double anEnergy)
void SetTransport(G4VIntraNuclearTransportModel *const value)
G4hMultipleScattering fAlphaMultipleScattering
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89
double G4double
Definition: G4Types.hh:76
G4hMultipleScattering fIonMultipleScattering
static G4He3 * He3()
Definition: G4He3.cc:94