Geant4  10.00.p03
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 // This is the *BASIC* version of Hadrontherapy, a Geant4-based application
27 // See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy
28 //
29 // Visit the Hadrontherapy web site (http://www.lns.infn.it/link/Hadrontherapy) to request
30 // the *COMPLETE* version of this program, together with its documentation;
31 // Hadrontherapy (both basic and full version) are supported by the Italian INFN
32 // Institute in the framework of the MC-INFN Group
33 //
34 //
35 // In this class the models for ion-ion interactions at intermediate energies (0 - 1 GeV per nucleon)
36 // can be activate. This class can be used alternatively to the "binary_ion" physics list
37 //
38 // The usefullness of this class is that you can explicitally see the total inelastic sections
39 // activated and the models called. Moreover you can choose to activate for ions (from deuteron
40 // to heavier nucleus) three different and exclusive models: the Binary Light Ion cascade, the QMD
41 // and The Wilson.
42 
43 // For hadrotherapy pouposes, where distributions of produced fragments is importante we strongly
44 // suggest to use Binary or QMD. The Binary model is the default and at moment, you can swith beetween models decommenting
45 // the line of code and recompiling
46 
48 #include "G4SystemOfUnits.hh"
49 #include "G4ParticleDefinition.hh"
50 #include "G4ProcessManager.hh"
51 
52 // Total cross section for inelastic processes
56 
61 #include "G4QMDReaction.hh"
62 #include "G4WilsonAbrasionModel.hh"
63 #include "G4IonInelasticProcess.hh"
65 
69 {
70  G4cout << G4endl
71  << "A local inelastic model is activated for all ions"
72  << G4endl;
73 }
74 
77 {}
78 
81 {
82  G4ParticleDefinition* particle = 0;
83  G4ProcessManager* processManager = 0;
84 
85  // ****************************************************************
86  // *** Ion-Ion models definition ***
87  // ****************************************************************
88  G4QMDReaction* JQMDmodel = new G4QMDReaction();
89  JQMDmodel -> SetMinEnergy(0*MeV);
90  JQMDmodel -> SetMaxEnergy(10*GeV);
91 
93  ligthBinary -> SetMinEnergy(0*MeV);
94  ligthBinary -> SetMaxEnergy(10*GeV);
95 
96  G4WilsonAbrasionModel* WilsonModel = new G4WilsonAbrasionModel();
97  WilsonModel -> SetUseAblation(true);
98  WilsonModel -> SetMinEnergy(0*MeV);
99  WilsonModel -> SetMaxEnergy(10 *GeV);
100 
101  G4TripathiCrossSection* TripatiCrossSections = new G4TripathiCrossSection;
102  G4TripathiLightCrossSection* TripatiLightCrossSections = new G4TripathiLightCrossSection;
103  G4IonsShenCrossSection* ShenCrossSections = new G4IonsShenCrossSection;
104 
105  // ****************
106  // *** Deuteron ***
107  // ****************
108  G4DeuteronInelasticProcess* deuteronInelasticProcess = new G4DeuteronInelasticProcess;
109 
110  deuteronInelasticProcess -> AddDataSet(ShenCrossSections);
111  deuteronInelasticProcess -> AddDataSet(TripatiCrossSections);
112  deuteronInelasticProcess -> AddDataSet(TripatiLightCrossSections);
113 
114  deuteronInelasticProcess -> RegisterMe(ligthBinary);
115  //deuteronInelasticProcess -> RegisterMe(JQMDmodel);
116  //deuteronInelasticProcess -> RegisterMe(WilsonModel);
117 
118  particle = G4Deuteron::Deuteron();
119  processManager = particle -> GetProcessManager();
120  processManager -> AddDiscreteProcess(deuteronInelasticProcess);
121 
122  // **************
123  // *** Triton ***
124  // **************
125  G4TritonInelasticProcess* tritonInelasticProcess = new G4TritonInelasticProcess;
126 
127  tritonInelasticProcess -> AddDataSet(ShenCrossSections);
128  tritonInelasticProcess -> AddDataSet(TripatiCrossSections);
129  tritonInelasticProcess -> AddDataSet(TripatiLightCrossSections);
130 
131  tritonInelasticProcess -> RegisterMe(ligthBinary);
132  //tritonInelasticProcess -> RegisterMe(JQMDmodel);
133  //tritonInelasticProcess -> RegisterMe(WilsonModel);
134 
135  particle = G4Triton::Triton();
136  processManager = particle -> GetProcessManager();
137  processManager -> AddDiscreteProcess(tritonInelasticProcess);
138 
139  // *************
140  // *** Alpha ***
141  // *************
142  G4AlphaInelasticProcess* alphaInelasticProcess = new G4AlphaInelasticProcess;
143 
144  alphaInelasticProcess -> AddDataSet(ShenCrossSections);
145  alphaInelasticProcess -> AddDataSet(TripatiCrossSections);
146  alphaInelasticProcess -> AddDataSet(TripatiLightCrossSections);
147 
148  alphaInelasticProcess -> RegisterMe(ligthBinary);
149  //alphaInelasticProcess -> RegisterMe(JQMDmodel);
150  //alphaIonInelasticProcess -> RegisterMe(WilsonModel);
151 
152  particle = G4Alpha::Alpha();
153  processManager = particle -> GetProcessManager();
154  processManager -> AddDiscreteProcess(alphaInelasticProcess);
155 
156  // *******************
157  // *** Generic Ion ***
158  // *******************
159  G4IonInelasticProcess* genericIonInelasticProcess = new G4IonInelasticProcess();
160 
161  genericIonInelasticProcess -> AddDataSet(ShenCrossSections);
162  genericIonInelasticProcess -> AddDataSet(TripatiCrossSections);
163  genericIonInelasticProcess -> AddDataSet(TripatiLightCrossSections);
164 
165  genericIonInelasticProcess -> RegisterMe(ligthBinary);
166  //genericIonInelasticProcess -> RegisterMe(JQMDmodel);
167  //genericIonInelasticProcess -> RegisterMe(WilsonModel);
168 
169  particle = G4GenericIon::GenericIon();
170  processManager = particle -> GetProcessManager();
171  processManager -> AddDiscreteProcess(genericIonInelasticProcess);
172 }
173 
174 
175 
static const double MeV
Definition: G4SIunits.hh:193
LocalIonIonInelasticPhysic(const G4String &name="local_ion_ion_inelastic")
G4String name
Definition: TRTMaterials.hh:40
G4GLOB_DLL std::ostream G4cout
static G4Triton * Triton()
Definition: G4Triton.cc:95
static const double GeV
Definition: G4SIunits.hh:196
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
static G4GenericIon * GenericIon()
Definition: G4GenericIon.cc:93
#define G4endl
Definition: G4ios.hh:61
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89