Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
LocalIonIonInelasticPhysic.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 // Hadrontherapy advanced example for Geant4
27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
28 //
29 //
30 // In this class the models for ion-ion interactions at intermediate energies (0 - 1 GeV per nucleon)
31 // can be activate. This class can be used alternatively to the "binary_ion" physics list
32 //
33 // The usefullness of this class is that you can explicitally see the total inelastic sections
34 // activated and the models called. Moreover you can choose to activate for ions (from deuteron
35 // to heavier nucleus) three different and exclusive models: the Binary Light Ion cascade, the QMD
36 // and The Wilson.
37 
38 // For hadrotherapy pouposes, where distributions of produced fragments is importante we strongly
39 // suggest to use Binary or QMD. The Binary model is the default and at moment, you can swith beetween models decommenting
40 // the line of code and recompiling
41 
43 #include "G4SystemOfUnits.hh"
44 #include "G4ParticleDefinition.hh"
45 #include "G4ProcessManager.hh"
46 
47 // Total cross section for inelastic processes
51 
56 #include "G4QMDReaction.hh"
57 #include "G4WilsonAbrasionModel.hh"
58 #include "G4IonInelasticProcess.hh"
60 
64 {
65  G4cout << G4endl
66  << "A local inelastic model is activated for all ions"
67  << G4endl;
68 }
69 
72 {}
73 
76 {
77  G4ParticleDefinition* particle = 0;
78  G4ProcessManager* processManager = 0;
79 
80  // ****************************************************************
81  // *** Ion-Ion models definition ***
82  // ****************************************************************
83  G4QMDReaction* JQMDmodel = new G4QMDReaction();
84  JQMDmodel -> SetMinEnergy(0*MeV);
85  JQMDmodel -> SetMaxEnergy(10*GeV);
86 
88  ligthBinary -> SetMinEnergy(0*MeV);
89  ligthBinary -> SetMaxEnergy(10*GeV);
90 
91  G4WilsonAbrasionModel* WilsonModel = new G4WilsonAbrasionModel();
92  WilsonModel -> SetUseAblation(true);
93  WilsonModel -> SetMinEnergy(0*MeV);
94  WilsonModel -> SetMaxEnergy(10 *GeV);
95 
96  G4TripathiCrossSection* TripatiCrossSections = new G4TripathiCrossSection;
97  G4TripathiLightCrossSection* TripatiLightCrossSections = new G4TripathiLightCrossSection;
98  G4IonsShenCrossSection* ShenCrossSections = new G4IonsShenCrossSection;
99 
100  // ****************
101  // *** Deuteron ***
102  // ****************
103  G4DeuteronInelasticProcess* deuteronInelasticProcess = new G4DeuteronInelasticProcess;
104 
105  deuteronInelasticProcess -> AddDataSet(ShenCrossSections);
106  deuteronInelasticProcess -> AddDataSet(TripatiCrossSections);
107  deuteronInelasticProcess -> AddDataSet(TripatiLightCrossSections);
108 
109  deuteronInelasticProcess -> RegisterMe(ligthBinary);
110  //deuteronInelasticProcess -> RegisterMe(JQMDmodel);
111  //deuteronInelasticProcess -> RegisterMe(WilsonModel);
112 
113  particle = G4Deuteron::Deuteron();
114  processManager = particle -> GetProcessManager();
115  processManager -> AddDiscreteProcess(deuteronInelasticProcess);
116 
117  // **************
118  // *** Triton ***
119  // **************
120  G4TritonInelasticProcess* tritonInelasticProcess = new G4TritonInelasticProcess;
121 
122  tritonInelasticProcess -> AddDataSet(ShenCrossSections);
123  tritonInelasticProcess -> AddDataSet(TripatiCrossSections);
124  tritonInelasticProcess -> AddDataSet(TripatiLightCrossSections);
125 
126  tritonInelasticProcess -> RegisterMe(ligthBinary);
127  //tritonInelasticProcess -> RegisterMe(JQMDmodel);
128  //tritonInelasticProcess -> RegisterMe(WilsonModel);
129 
130  particle = G4Triton::Triton();
131  processManager = particle -> GetProcessManager();
132  processManager -> AddDiscreteProcess(tritonInelasticProcess);
133 
134  // *************
135  // *** Alpha ***
136  // *************
137  G4AlphaInelasticProcess* alphaInelasticProcess = new G4AlphaInelasticProcess;
138 
139  alphaInelasticProcess -> AddDataSet(ShenCrossSections);
140  alphaInelasticProcess -> AddDataSet(TripatiCrossSections);
141  alphaInelasticProcess -> AddDataSet(TripatiLightCrossSections);
142 
143  alphaInelasticProcess -> RegisterMe(ligthBinary);
144  //alphaInelasticProcess -> RegisterMe(JQMDmodel);
145  //alphaIonInelasticProcess -> RegisterMe(WilsonModel);
146 
147  particle = G4Alpha::Alpha();
148  processManager = particle -> GetProcessManager();
149  processManager -> AddDiscreteProcess(alphaInelasticProcess);
150 
151  // *******************
152  // *** Generic Ion ***
153  // *******************
154  G4IonInelasticProcess* genericIonInelasticProcess = new G4IonInelasticProcess();
155 
156  genericIonInelasticProcess -> AddDataSet(ShenCrossSections);
157  genericIonInelasticProcess -> AddDataSet(TripatiCrossSections);
158  genericIonInelasticProcess -> AddDataSet(TripatiLightCrossSections);
159 
160  genericIonInelasticProcess -> RegisterMe(ligthBinary);
161  //genericIonInelasticProcess -> RegisterMe(JQMDmodel);
162  //genericIonInelasticProcess -> RegisterMe(WilsonModel);
163 
164  particle = G4GenericIon::GenericIon();
165  processManager = particle -> GetProcessManager();
166  processManager -> AddDiscreteProcess(genericIonInelasticProcess);
167 }
168 
169 
170 
const XML_Char * name
Definition: expat.h:151
LocalIonIonInelasticPhysic(const G4String &name="local_ion_ion_inelastic")
G4GLOB_DLL std::ostream G4cout
static G4Triton * Triton()
Definition: G4Triton.cc:95
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
static G4GenericIon * GenericIon()
Definition: G4GenericIon.cc:93
static constexpr double GeV
Definition: G4SIunits.hh:217
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89