Geant4
v4-10.4-release
메인 페이지
관련된 페이지
모듈
네임스페이스
클래스
파일들
파일 목록
파일 멤버
모두
클래스
네임스페이스들
파일들
함수
변수
타입정의
열거형 타입
열거형 멤버
Friends
매크로
그룹들
페이지들
examples
extended
medical
dna
wholeNuclearDNA
src
extended/medical/dna/wholeNuclearDNA/src/G4ElectronCapture.cc
이 파일의 문서화 페이지로 가기
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
//---------------------------------------------------------------------------
27
//
28
// ClassName: G4ElectronCapture
29
//
30
// Description: The process to kill particles to save CPU
31
//
32
// Author: V.Ivanchenko 31 August 2010
33
//
34
//----------------------------------------------------------------------------
35
//
36
// This example is provided by the Geant4-DNA collaboration
37
// Any report or published results obtained using the Geant4-DNA software
38
// and the DNA geometry given in the Geom_DNA example
39
// shall cite the following Geant4-DNA collaboration publications:
40
// [1] NIM B 298 (2013) 47-54
41
// [2] Med. Phys. 37 (2010) 4692-4708
42
// The Geant4-DNA web site is available at http://geant4-dna.org
43
//
44
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45
//
46
// $ID$
49
50
#include "G4ElectronCapture.hh"
51
#include "
G4SystemOfUnits.hh
"
52
#include "
G4ParticleDefinition.hh
"
53
#include "
G4Step.hh
"
54
#include "
G4Track.hh
"
55
#include "
G4Region.hh
"
56
#include "
G4RegionStore.hh
"
57
#include "
G4Electron.hh
"
58
59
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60
61
G4ElectronCapture::G4ElectronCapture
(
const
G4String
& regName,
G4double
ekinlim)
62
:
G4VDiscreteProcess
(
"eCapture"
,
fElectromagnetic
),
63
fKinEnergyThreshold(ekinlim),
64
fRegionName(regName),
65
region(0)
66
{
67
if
(regName ==
""
|| regName ==
"world"
)
68
{
69
fRegionName
=
"DefaultRegionForTheWorld"
;
70
}
71
pParticleChange
= &
fParticleChange
;
72
}
73
74
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
75
76
G4ElectronCapture::~G4ElectronCapture
()
77
{
78
}
79
80
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
81
82
void
83
G4ElectronCapture::SetKinEnergyLimit
(
G4double
val)
84
{
85
fKinEnergyThreshold
= val;
86
if
(
verboseLevel
> 0)
87
{
88
G4cout
<<
"### G4ElectronCapture: Tracking cut E(MeV) = "
89
<<
fKinEnergyThreshold
/
MeV
<<
G4endl
;
90
}
91
}
92
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
93
94
void
95
G4ElectronCapture::BuildPhysicsTable
(
const
G4ParticleDefinition
&)
96
{
97
region
= (
G4RegionStore::GetInstance
())->GetRegion(
fRegionName
);
98
if
(
region
&&
verboseLevel
> 0)
99
{
100
G4cout
<<
"### G4ElectronCapture: Tracking cut E(MeV) = "
101
<<
fKinEnergyThreshold
/
MeV
<<
" is assigned to "
<<
fRegionName
102
<<
G4endl
;
103
}
104
}
105
106
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107
108
G4bool
109
G4ElectronCapture::IsApplicable
(
const
G4ParticleDefinition
&)
110
{
111
return
true
;
112
}
113
114
G4double
115
G4ElectronCapture::PostStepGetPhysicalInteractionLength
(
const
G4Track
& aTrack,
116
G4double
,
117
G4ForceCondition
*
condition
)
118
{
119
// condition is set to "Not Forced"
120
*condition =
NotForced
;
121
122
G4double
limit =
DBL_MAX
;
123
if
(
region
)
124
{
125
if
(aTrack.
GetVolume
()->
GetLogicalVolume
()->
GetRegion
() ==
region
126
&& aTrack.
GetKineticEnergy
() <
fKinEnergyThreshold
)
127
{
128
limit = 0.0;
129
}
130
}
131
return
limit;
132
}
133
134
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
135
136
G4VParticleChange
*
137
G4ElectronCapture::PostStepDoIt
(
const
G4Track
& aTrack,
const
G4Step
&)
138
{
139
pParticleChange
->
Initialize
(aTrack);
140
pParticleChange
->
ProposeTrackStatus
(
fStopAndKill
);
141
pParticleChange
->
ProposeLocalEnergyDeposit
(aTrack.
GetKineticEnergy
());
142
fParticleChange
.
SetProposedKineticEnergy
(0.0);
143
return
pParticleChange
;
144
}
145
146
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
147
148
G4double
149
G4ElectronCapture::GetMeanFreePath
(
const
G4Track
&,
G4double
,
G4ForceCondition
*)
150
{
151
return
DBL_MAX
;
152
}
153
154
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
155
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4Electron.hh
G4VPhysicalVolume::GetLogicalVolume
G4LogicalVolume * GetLogicalVolume() const
MeV
static constexpr double MeV
Definition:
G4SIunits.hh:214
G4ElectronCapture::PostStepGetPhysicalInteractionLength
virtual G4double PostStepGetPhysicalInteractionLength(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:98
G4VDiscreteProcess
Definition:
G4VDiscreteProcess.hh:58
G4endl
#define G4endl
Definition:
G4ios.hh:61
G4ElectronCapture::BuildPhysicsTable
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:80
G4ElectronCapture::G4ElectronCapture
G4ElectronCapture(const G4String ®Name, G4double ekinlimit)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:53
fStopAndKill
Definition:
G4TrackStatus.hh:56
G4String
Definition:
examples/extended/parallel/TopC/ParN02/AnnotatedFiles/G4String.hh:45
G4ParticleChangeForGamma::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition:
G4ParticleChangeForGamma.hh:129
G4ElectronCapture::GetMeanFreePath
virtual G4double GetMeanFreePath(const G4Track &, G4double, G4ForceCondition *)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:127
fElectromagnetic
Definition:
G4ProcessType.hh:47
G4VParticleChange
Definition:
G4VParticleChange.hh:94
condition
G4double condition(const G4ErrorSymMatrix &m)
G4ElectronCapture::PostStepDoIt
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:115
G4ElectronCapture::region
G4Region * region
Definition:
advanced/microelectronics/include/G4ElectronCapture.hh:97
G4VParticleChange::Initialize
virtual void Initialize(const G4Track &)
G4ElectronCapture::fKinEnergyThreshold
G4double fKinEnergyThreshold
Definition:
extended/medical/dna/microdosimetry/include/G4ElectronCapture.hh:98
G4double
double G4double
Definition:
G4Types.hh:76
G4bool
bool G4bool
Definition:
G4Types.hh:79
G4VProcess::pParticleChange
G4VParticleChange * pParticleChange
Definition:
G4VProcess.hh:283
G4Track
Definition:
G4Track.hh:76
G4ElectronCapture::~G4ElectronCapture
virtual ~G4ElectronCapture()
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:65
NotForced
Definition:
G4ForceCondition.hh:56
G4SystemOfUnits.hh
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:73
G4Step
Definition:
G4Step.hh:76
G4ElectronCapture::fRegionName
G4String fRegionName
Definition:
extended/medical/dna/microdosimetry/include/G4ElectronCapture.hh:99
G4Track.hh
G4RegionStore::GetInstance
static G4RegionStore * GetInstance()
Definition:
G4RegionStore.cc:163
G4VProcess::verboseLevel
G4int verboseLevel
Definition:
G4VProcess.hh:371
G4ForceCondition
G4ForceCondition
Definition:
G4ForceCondition.hh:49
G4cout
G4GLOB_DLL std::ostream G4cout
G4ElectronCapture::IsApplicable
virtual G4bool IsApplicable(const G4ParticleDefinition &)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:92
G4RegionStore.hh
G4Step.hh
G4ParticleDefinition.hh
G4ElectronCapture::fParticleChange
G4ParticleChangeForGamma fParticleChange
Definition:
advanced/microelectronics/include/G4ElectronCapture.hh:98
G4Region.hh
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
DBL_MAX
#define DBL_MAX
Definition:
templates.hh:83
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4LogicalVolume::GetRegion
G4Region * GetRegion() const
G4ElectronCapture::SetKinEnergyLimit
void SetKinEnergyLimit(G4double)
Definition:
advanced/microelectronics/src/G4ElectronCapture.cc:70
G4Track::GetVolume
G4VPhysicalVolume * GetVolume() const
다음에 의해 생성됨 :
1.8.5