ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
G4RKG3_Stepper.hh
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4RKG3_Stepper.hh
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
// G4RKG3_Stepper
27
//
28
// Class description:
29
//
30
// Runga-Kutta integrator stepper from Geant3.
31
32
// Created: J.Apostolakis, V.Grichine - 30.01.1997
33
// -------------------------------------------------------------------
34
#ifndef G4RKG3_STEPPER_HH
35
#define G4RKG3_STEPPER_HH
36
37
#include "
G4Types.hh
"
38
#include "
G4MagIntegratorStepper.hh
"
39
#include "
G4ThreeVector.hh
"
40
41
class
G4Mag_EqRhs
;
42
43
class
G4RKG3_Stepper
:
public
G4MagIntegratorStepper
44
{
45
public
:
// with description
46
47
G4RKG3_Stepper
(
G4Mag_EqRhs
* EqRhs);
48
// Integrate over 6 variables only: position & velocity.
49
// Not implemented yet !
50
51
~G4RKG3_Stepper
();
52
53
void
Stepper
(
const
G4double
yIn[],
54
const
G4double
dydx[],
55
G4double
h
,
56
G4double
yOut[],
57
G4double
yErr[] );
58
// The method which must be provided, even if less efficient.
59
60
G4double
DistChord
()
const
;
61
62
void
StepNoErr
(
const
G4double
tIn[8],
63
const
G4double
dydx[6],
64
G4double
Step,
65
G4double
tOut[8],
66
G4double
B
[3] );
67
// Integrator RK Stepper from G3 with only two field evaluation per
68
// Step. It is used in propagation initial Step by small substeps
69
// after solution error and delta geometry considerations.
70
// B[3] is magnetic field which is passed from substep to substep.
71
72
void
StepWithEst
(
const
G4double
tIn[8],
73
const
G4double
dydx[6],
74
G4double
Step,
75
G4double
tOut[8],
76
G4double
& alpha2,
77
G4double
& beta2,
78
const
G4double
B1[3],
79
G4double
B2[3] );
80
// Integrator for RK from G3 with evaluation of error in solution and
81
// delta geometry based on naive similarity with the case of uniform
82
// magnetic field.
83
// B1[3] is input and is the first magnetic field values
84
// B2[3] is output and is the final magnetic field values.
85
86
public
:
// without description
87
88
inline
G4int
IntegratorOrder
()
const
{
return
4; }
89
90
private
:
91
92
G4ThreeVector
fyInitial
,
93
fyMidPoint
,
94
fyFinal
;
95
G4ThreeVector
fpInitial
;
96
G4ThreeVector
BfldIn
;
97
G4double
hStep
= 0.0;
98
};
99
100
#endif
geant4
tree
geant4-10.6-release
source
geometry
magneticfield
include
G4RKG3_Stepper.hh
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:16
using
1.8.2 with
ECCE GitHub integration