Geant4-11
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4VRangeToEnergyConverter.hh
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// G4VRangeToEnergyConverter
27//
28// Class description:
29//
30// Base class for Range to Energy Converters.
31// Cut in energy corresponding to given cut value in range
32// is calculated for a material by using Convert() method.
33
34// Author: H.Kurashige, 05 October 2002 - First implementation
35// --------------------------------------------------------------------
36#ifndef G4VRangeToEnergyConverter_hh
37#define G4VRangeToEnergyConverter_hh 1
38
39#include <vector>
40
41#include "globals.hh"
43#include "G4Material.hh"
44#include "G4Threading.hh"
45
47{
48public:
49
51
53
54 // operators are not used
57 (const G4VRangeToEnergyConverter &r) = delete;
60
61 // Calculate energy cut from given range cut for the material
62 virtual G4double Convert(const G4double rangeCut, const G4Material* material);
63
64 // Set energy range for all particle type
65 // if highedge > 10 GeV, highedge value is not changed
66 static void SetEnergyRange(const G4double lowedge, const G4double highedge);
67
68 // Get energy range for all particle type
71
72 // Get/set max cut energy for all particle type
73 // No check on the value
75 static void SetMaxEnergyCut(const G4double value);
76
77 // Return pointer to the particle type which this converter takes care of
78 inline const G4ParticleDefinition* GetParticleType() const;
79
80 inline void SetVerboseLevel(G4int value);
81 inline G4int GetVerboseLevel() const;
82 // control flag for output message
83 // 0: Silent
84 // 1: Warning message
85 // 2: More
86
87protected:
88
89 virtual G4double ComputeValue(const G4int Z, const G4double kinEnergy) = 0;
90
91private:
92
93 static void FillEnergyVector(const G4double emin, const G4double emax);
94
95 G4double ConvertForGamma(const G4double rangeCut, const G4Material* material);
96
97 G4double ConvertForElectron(const G4double rangeCut,
98 const G4Material* material);
99
100 inline G4double LiniearInterpolation(const G4double e1, const G4double e2,
101 const G4double r1, const G4double r2,
102 const G4double r);
103
104protected:
105
106#ifdef G4MULTITHREADED
107 static G4Mutex theMutex;
108#endif
109
111 static G4double Emax;
112 static std::vector<G4double>* Energy;
114 static G4int Nbin;
115
117
121};
122
123// ------------------
124// Inline methods
125// ------------------
126
127inline
129{
130 verboseLevel = value;
131}
132
133inline
135{
136 return verboseLevel;
137}
138
139inline
141{
142 return theParticle;
143}
144
146 const G4double e1, const G4double e2,
147 const G4double r1, const G4double r2, const G4double r)
148{
149 return (r1 == r2) ? e1 : e1 + (e2 - e1)*(r - r1)/(r2 - r1);
150}
151
152#endif
static const G4double e1[44]
static const G4double e2[44]
static const G4double emax
std::mutex G4Mutex
Definition: G4Threading.hh:81
double G4double
Definition: G4Types.hh:83
bool G4bool
Definition: G4Types.hh:86
int G4int
Definition: G4Types.hh:85
const G4int Z[17]
static void SetMaxEnergyCut(const G4double value)
virtual G4double ComputeValue(const G4int Z, const G4double kinEnergy)=0
G4double LiniearInterpolation(const G4double e1, const G4double e2, const G4double r1, const G4double r2, const G4double r)
G4double ConvertForGamma(const G4double rangeCut, const G4Material *material)
static std::vector< G4double > * Energy
const G4ParticleDefinition * GetParticleType() const
G4double ConvertForElectron(const G4double rangeCut, const G4Material *material)
static void FillEnergyVector(const G4double emin, const G4double emax)
virtual G4double Convert(const G4double rangeCut, const G4Material *material)
G4bool operator==(const G4VRangeToEnergyConverter &r) const =delete
static void SetEnergyRange(const G4double lowedge, const G4double highedge)
G4VRangeToEnergyConverter(const G4VRangeToEnergyConverter &r)=delete
const G4ParticleDefinition * theParticle
G4bool operator!=(const G4VRangeToEnergyConverter &r) const =delete
string material
Definition: eplot.py:19