Geant4_10
G4mplIonisation.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: G4mplIonisation.cc 66241 2012-12-13 18:34:42Z gunter $
27 //
28 // -------------------------------------------------------------------
29 //
30 // GEANT4 Class file
31 //
32 //
33 // File name: G4mplIonisation
34 //
35 // Author: Vladimir Ivanchenko
36 //
37 // Creation date: 25.08.2005
38 //
39 // Modifications:
40 //
41 //
42 // -------------------------------------------------------------------
43 //
44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46 
47 #include "G4mplIonisation.hh"
48 #include "G4PhysicalConstants.hh"
49 #include "G4SystemOfUnits.hh"
50 #include "G4Electron.hh"
51 #include "G4mplIonisationModel.hh"
53 #include "G4LossTableManager.hh"
54 
55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
56 
57 using namespace std;
58 
60  : G4VEnergyLossProcess(name),
61  magneticCharge(mCharge),
62  isInitialised(false)
63 {
64  // By default classical magnetic charge is used
65  if(magneticCharge == 0.0) { magneticCharge = eplus*0.5/fine_structure_const; }
66 
67  SetVerboseLevel(0);
69  SetStepFunction(0.2, 1*mm);
71 }
72 
73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
74 
76 {}
77 
78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
79 
81 {
82  return true;
83 }
84 
85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
86 
88  const G4ParticleDefinition*)
89 {
90  if(isInitialised) { return; }
91 
92  SetBaseParticle(0);
93 
94  // monopole model is responsible both for energy loss and fluctuations
96  new G4mplIonisationWithDeltaModel(magneticCharge,"PAI");
97  ion->SetParticle(p);
98 
99  // define size of dedx and range tables
100  G4double emin = std::min(MinKinEnergy(),ion->LowEnergyLimit());
102  G4int bin = G4lrint(G4LossTableManager::Instance()->GetNumberOfBinsPerDecade()
103  *std::log10(emax/emin));
104  ion->SetLowEnergyLimit(emin);
105  ion->SetHighEnergyLimit(emax);
106  SetMinKinEnergy(emin);
107  SetMaxKinEnergy(emax);
108  SetDEDXBinning(bin);
109 
110  AddEmModel(1,ion,ion);
111 
112  isInitialised = true;
113 }
114 
115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
116 
118 {
119  G4cout << " No delta-electron production, only dE/dx"
120  << G4endl;
121 }
122 
123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:599
static G4LossTableManager * Instance()
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:592
virtual void PrintInfo()
const char * p
Definition: xmltok.h:285
void SetStepFunction(G4double v1, G4double v2)
void SetParticle(const G4ParticleDefinition *p)
float bin[41]
Definition: plottest35.C:14
const XML_Char * name
Definition: expat.h:151
int G4int
Definition: G4Types.hh:78
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:683
virtual void InitialiseEnergyLossProcess(const G4ParticleDefinition *, const G4ParticleDefinition *)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
void SetMaxKinEnergy(G4double e)
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=0)
void SetSecondaryParticle(const G4ParticleDefinition *p)
G4mplIonisation(G4double mCharge=0.0, const G4String &name="mplIoni")
int G4lrint(double ad)
Definition: templates.hh:163
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual ~G4mplIonisation()
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double MinKinEnergy() const
G4double MaxKinEnergy() const
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
virtual G4bool IsApplicable(const G4ParticleDefinition &p)
void SetBaseParticle(const G4ParticleDefinition *p)
void SetVerboseLevel(G4int value)
Definition: G4VProcess.hh:437
void SetDEDXBinning(G4int nbins)
void SetMinKinEnergy(G4double e)