Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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 85013 2014-10-23 09:45:07Z gcosmo $
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 "G4EmParameters.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
101  G4double emin = std::min(param->MinKinEnergy(),ion->LowEnergyLimit());
103  G4int bin = G4lrint(param->NumberOfBinsPerDecade()*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 
120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
virtual void InitialiseEnergyLossProcess(const G4ParticleDefinition *, const G4ParticleDefinition *) override
G4int NumberOfBinsPerDecade() const
const XML_Char * name
Definition: expat.h:151
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:640
virtual void PrintInfo() override
G4double MaxKinEnergy() const
static constexpr double mm
Definition: G4SIunits.hh:115
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:633
const char * p
Definition: xmltok.h:285
void SetParticle(const G4ParticleDefinition *p)
int G4int
Definition: G4Types.hh:78
void SetStepFunction(G4double v1, G4double v2, G4bool lock=true)
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:724
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=nullptr)
virtual G4bool IsApplicable(const G4ParticleDefinition &p) override
bool G4bool
Definition: G4Types.hh:79
static constexpr double eplus
Definition: G4SIunits.hh:199
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
void SetMaxKinEnergy(G4double e)
G4double MinKinEnergy() const
void SetSecondaryParticle(const G4ParticleDefinition *p)
static const G4double emax
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
static G4EmParameters * Instance()
virtual ~G4mplIonisation()
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static G4Electron * Electron()
Definition: G4Electron.cc:94
double G4double
Definition: G4Types.hh:76
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:731
static constexpr double fine_structure_const
void SetBaseParticle(const G4ParticleDefinition *p)
void SetVerboseLevel(G4int value)
Definition: G4VProcess.hh:437
void SetDEDXBinning(G4int nbins)
void SetMinKinEnergy(G4double e)