Geant4  10.00.p02
G4DNAExcitation.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: G4DNAExcitation.cc 70171 2013-05-24 13:34:18Z gcosmo $
27 
28 #include "G4DNAExcitation.hh"
29 #include "G4SystemOfUnits.hh"
30 
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 
33 using namespace std;
34 
36  G4ProcessType type):G4VEmProcess (processName, type),
37  isInitialised(false)
38 {
40 }
41 
42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 
45 {}
46 
47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
48 
50 {
51 
54 
55  return
56  (
57  &p == G4Electron::Electron()
59  || &p == instance->GetIon("hydrogen")
60  || &p == instance->GetIon("alpha++")
61  || &p == instance->GetIon("alpha+")
62  || &p == instance->GetIon("helium")
63  );
64 }
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
67 
69 {
70  if(!isInitialised)
71  {
72  isInitialised = true;
73  SetBuildTableFlag(false);
74 
76 
77  if(name == "e-")
78  {
79 
80  // Emfietzoglou model
81 /*
82  if(!EmModel()) SetEmModel(new G4DNAEmfietzoglouExcitationModel);
83  EmModel()->SetLowEnergyLimit(8.23*eV);
84  EmModel()->SetHighEnergyLimit(10*MeV);
85 */
86  // Born model
87 
91 
92  AddEmModel(1, EmModel());
93  }
94 
95  if(name == "proton")
96  {
100 
102  EmModel(2)->SetLowEnergyLimit(500*keV);
103  EmModel(2)->SetHighEnergyLimit(100*MeV);
104 
105  AddEmModel(1, EmModel(1));
106  AddEmModel(2, EmModel(2));
107  }
108 
109  if(name == "hydrogen")
110  {
114 
115  AddEmModel(1, EmModel());
116  }
117 
118 
119  if( name == "alpha" || name == "alpha+" || name == "helium" )
120  {
124 
125  AddEmModel(1, EmModel());
126  }
127 
128  }
129 }
130 
131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132 
134 {
135  if (EmModel(2))
136  {
137  G4cout
138  << " Total cross sections computed from "
139  << EmModel(1)->GetName()
140  << " and "
141  << EmModel(2)->GetName()
142  << " models"
143  << G4endl;
144  }
145  else
146  {
147  G4cout
148  << " Total cross sections computed from "
149  << EmModel()->GetName()
150  << G4endl;
151  }
152 }
153 
154 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static const double MeV
Definition: G4SIunits.hh:193
void SetBuildTableFlag(G4bool val)
virtual void InitialiseProcess(const G4ParticleDefinition *)
G4DNAExcitation(const G4String &processName="DNAExcitation", G4ProcessType type=fElectromagnetic)
G4VEmModel * EmModel(G4int index=1) const
G4String name
Definition: TRTMaterials.hh:40
virtual void PrintInfo()
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:683
virtual ~G4DNAExcitation()
void SetEmModel(G4VEmModel *, G4int index=1)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
static G4DNAGenericIonsManager * Instance(void)
static const double eV
Definition: G4SIunits.hh:194
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
const G4String & GetName() const
Definition: G4VEmModel.hh:753
static MCTruthManager * instance
static const double keV
Definition: G4SIunits.hh:195
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
G4ParticleDefinition * GetIon(const G4String &name)
G4ProcessType
virtual G4bool IsApplicable(const G4ParticleDefinition &)