Geant4-11
source
processes
electromagnetic
adjoint
include
G4AdjointAlongStepWeightCorrection.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
//
27
// Class: G4AdjointAlongStepWeightCorrection
28
// Author: L. Desorgher
29
// Organisation: SpaceIT GmbH
31
//
32
// Documentation:
33
//
34
// Continuous processes act on adjoint particles to continuously correct their
35
// weight during the adjoint reverse tracking. This process is needed when
36
// the adjoint cross sections are not scaled such that the total adjoint cross
37
// section matches the total forward cross section. By default the mode where
38
// the total adjoint cross section is equal to the total forward cross section
39
// is used and therefore this along step weightcorrection factor is 1. However
40
// in some cases (some energy ranges) the total forward cross section or the
41
// total adjoint cross section can be zero. In this case the along step weight
42
// correction is needed and is given by exp(-(Sigma_tot_adj-Sigma_tot_fwd).dx)
43
//
44
//-------------------------------------------------------------
45
46
#ifndef G4AdjointAlongStepWeightCorrection_h
47
#define G4AdjointAlongStepWeightCorrection_h 1
48
49
#include "
globals.hh
"
50
#include "
G4VContinuousProcess.hh
"
51
52
class
G4AdjointCSManager
;
53
class
G4MaterialCutsCouple
;
54
class
G4ParticleDefinition
;
55
class
G4ParticleChange
;
56
class
G4Step
;
57
class
G4Track
;
58
59
class
G4AdjointAlongStepWeightCorrection
:
public
G4VContinuousProcess
60
{
61
public
:
62
explicit
G4AdjointAlongStepWeightCorrection
(
63
const
G4String
&
name
=
"ContinuousWeightCorrection"
,
64
G4ProcessType
type =
fElectromagnetic
);
65
66
~G4AdjointAlongStepWeightCorrection
()
override
;
67
68
G4VParticleChange
*
AlongStepDoIt
(
const
G4Track
&,
const
G4Step
&)
override
;
69
70
void
ProcessDescription
(std::ostream&)
const override
;
71
void
DumpInfo
()
const override
{
ProcessDescription
(
G4cout
); };
72
73
G4AdjointAlongStepWeightCorrection
(
G4AdjointAlongStepWeightCorrection
&) =
74
delete
;
75
G4AdjointAlongStepWeightCorrection
&
operator=
(
76
const
G4AdjointAlongStepWeightCorrection
& right) =
delete
;
77
78
protected
:
79
G4double
GetContinuousStepLimit
(
const
G4Track
& track,
80
G4double
previousStepSize,
81
G4double
currentMinimumStep,
82
G4double
& currentSafety)
override
;
83
84
private
:
85
void
DefineMaterial
(
const
G4MaterialCutsCouple
* couple);
86
87
const
G4MaterialCutsCouple
*
fCurrentCouple
=
nullptr
;
88
G4AdjointCSManager
*
fCSManager
=
nullptr
;
89
G4ParticleChange
*
fParticleChange
;
90
91
G4double
fPreStepKinEnergy
= 1.;
92
};
93
94
inline
void
G4AdjointAlongStepWeightCorrection::DefineMaterial
(
95
const
G4MaterialCutsCouple
* couple)
96
{
97
if
(couple !=
fCurrentCouple
)
98
{
99
fCurrentCouple
= couple;
100
}
101
}
102
103
#endif
G4ProcessType
G4ProcessType
Definition:
G4ProcessType.hh:38
fElectromagnetic
@ fElectromagnetic
Definition:
G4ProcessType.hh:41
G4double
double G4double
Definition:
G4Types.hh:83
G4VContinuousProcess.hh
G4cout
G4GLOB_DLL std::ostream G4cout
G4AdjointAlongStepWeightCorrection
Definition:
G4AdjointAlongStepWeightCorrection.hh:60
G4AdjointAlongStepWeightCorrection::fParticleChange
G4ParticleChange * fParticleChange
Definition:
G4AdjointAlongStepWeightCorrection.hh:89
G4AdjointAlongStepWeightCorrection::G4AdjointAlongStepWeightCorrection
G4AdjointAlongStepWeightCorrection(G4AdjointAlongStepWeightCorrection &)=delete
G4AdjointAlongStepWeightCorrection::GetContinuousStepLimit
G4double GetContinuousStepLimit(const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double ¤tSafety) override
Definition:
G4AdjointAlongStepWeightCorrection.cc:106
G4AdjointAlongStepWeightCorrection::ProcessDescription
void ProcessDescription(std::ostream &) const override
Definition:
G4AdjointAlongStepWeightCorrection.cc:51
G4AdjointAlongStepWeightCorrection::DefineMaterial
void DefineMaterial(const G4MaterialCutsCouple *couple)
Definition:
G4AdjointAlongStepWeightCorrection.hh:94
G4AdjointAlongStepWeightCorrection::fPreStepKinEnergy
G4double fPreStepKinEnergy
Definition:
G4AdjointAlongStepWeightCorrection.hh:91
G4AdjointAlongStepWeightCorrection::~G4AdjointAlongStepWeightCorrection
~G4AdjointAlongStepWeightCorrection() override
Definition:
G4AdjointAlongStepWeightCorrection.cc:45
G4AdjointAlongStepWeightCorrection::DumpInfo
void DumpInfo() const override
Definition:
G4AdjointAlongStepWeightCorrection.hh:71
G4AdjointAlongStepWeightCorrection::AlongStepDoIt
G4VParticleChange * AlongStepDoIt(const G4Track &, const G4Step &) override
Definition:
G4AdjointAlongStepWeightCorrection.cc:68
G4AdjointAlongStepWeightCorrection::G4AdjointAlongStepWeightCorrection
G4AdjointAlongStepWeightCorrection(const G4String &name="ContinuousWeightCorrection", G4ProcessType type=fElectromagnetic)
Definition:
G4AdjointAlongStepWeightCorrection.cc:36
G4AdjointAlongStepWeightCorrection::operator=
G4AdjointAlongStepWeightCorrection & operator=(const G4AdjointAlongStepWeightCorrection &right)=delete
G4AdjointAlongStepWeightCorrection::fCurrentCouple
const G4MaterialCutsCouple * fCurrentCouple
Definition:
G4AdjointAlongStepWeightCorrection.hh:87
G4AdjointAlongStepWeightCorrection::fCSManager
G4AdjointCSManager * fCSManager
Definition:
G4AdjointAlongStepWeightCorrection.hh:88
G4AdjointCSManager
Definition:
G4AdjointCSManager.hh:59
G4MaterialCutsCouple
Definition:
G4MaterialCutsCouple.hh:53
G4ParticleChange
Definition:
G4ParticleChange.hh:57
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:61
G4Step
Definition:
G4Step.hh:62
G4String
Definition:
G4String.hh:62
G4Track
Definition:
G4Track.hh:67
G4VContinuousProcess
Definition:
G4VContinuousProcess.hh:46
G4VParticleChange
Definition:
G4VParticleChange.hh:71
globals.hh
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition:
G4InuclParticleNames.hh:76
Generated by
1.9.3