Geant4-11
Public Member Functions | Data Fields | Private Types | Private Attributes | Static Private Attributes
G4TMagFieldEquation< T_Field > Class Template Reference

#include <G4TMagFieldEquation.hh>

Inheritance diagram for G4TMagFieldEquation< T_Field >:
G4Mag_UsualEqRhs G4Mag_EqRhs G4EquationOfMotion

Public Member Functions

 __attribute__ ((always_inline)) void RightHandSide(const G4double y[]
 
void EvaluateRhsGivenB (const G4double y[], const G4double B[3], G4double dydx[]) const
 
void EvaluateRhsReturnB (const G4double y[], G4double dydx[], G4double Field[]) const
 
G4double FCof () const
 
 G4TMagFieldEquation (T_Field *f)
 
G4FieldGetFieldObj ()
 
const G4FieldGetFieldObj () const
 
void GetFieldValue (const G4double Point[4], G4double Field[]) const
 
 GetFieldValue (PositionAndTime, Field)
 
void RightHandSide (const G4double y[], G4double dydx[]) const
 
virtual void SetChargeMomentumMass (G4ChargeState particleCharge, G4double MomentumXc, G4double mass)
 
void SetFieldObj (G4Field *pField)
 
void TEvaluateRhsGivenB (const G4double y[], const G4double B[3], G4double dydx[]) const
 
 TEvaluateRhsGivenB (y, Field, dydx)
 
virtual ~G4TMagFieldEquation ()
 

Data Fields

G4double dydx []
 
 PositionAndTime [0] = y[0]
 
G4double PositionAndTime [4]
 

Private Types

enum  { G4maximum_number_of_field_components = 24 }
 

Private Attributes

G4double fCof_val = 0.0
 
T_Field * itsField
 

Static Private Attributes

static const G4double fUnitConstant = 0.299792458 * (GeV/(tesla*m))
 

Detailed Description

template<class T_Field>
class G4TMagFieldEquation< T_Field >

Definition at line 44 of file G4TMagFieldEquation.hh.

Member Enumeration Documentation

◆ anonymous enum

template<class T_Field >
anonymous enum
private
Enumerator
G4maximum_number_of_field_components 

Definition at line 96 of file G4TMagFieldEquation.hh.

Constructor & Destructor Documentation

◆ G4TMagFieldEquation()

template<class T_Field >
G4TMagFieldEquation< T_Field >::G4TMagFieldEquation ( T_Field *  f)
inline

Definition at line 48 of file G4TMagFieldEquation.hh.

50 {
51 itsField = f;
52 }
G4Mag_UsualEqRhs(G4MagneticField *MagField)

References G4TMagFieldEquation< T_Field >::itsField.

◆ ~G4TMagFieldEquation()

template<class T_Field >
virtual G4TMagFieldEquation< T_Field >::~G4TMagFieldEquation ( )
inlinevirtual

Definition at line 54 of file G4TMagFieldEquation.hh.

54{;}

Member Function Documentation

◆ __attribute__()

template<class T_Field >
G4TMagFieldEquation< T_Field >::__attribute__ ( (always_inline)  ) const

◆ EvaluateRhsGivenB()

void G4Mag_UsualEqRhs::EvaluateRhsGivenB ( const G4double  y[],
const G4double  B[3],
G4double  dydx[] 
) const
virtualinherited

Implements G4Mag_EqRhs.

Definition at line 46 of file G4Mag_UsualEqRhs.cc.

49{
50 G4double momentum_mag_square = y[3]*y[3] + y[4]*y[4] + y[5]*y[5];
51 G4double inv_momentum_magnitude = 1.0 / std::sqrt( momentum_mag_square );
52
53 G4double cof = FCof()*inv_momentum_magnitude;
54
55 dydx[0] = y[3]*inv_momentum_magnitude; // (d/ds)x = Vx/V
56 dydx[1] = y[4]*inv_momentum_magnitude; // (d/ds)y = Vy/V
57 dydx[2] = y[5]*inv_momentum_magnitude; // (d/ds)z = Vz/V
58
59 dydx[3] = cof*(y[4]*B[2] - y[5]*B[1]) ; // Ax = a*(Vy*Bz - Vz*By)
60 dydx[4] = cof*(y[5]*B[0] - y[3]*B[2]) ; // Ay = a*(Vz*Bx - Vx*Bz)
61 dydx[5] = cof*(y[3]*B[1] - y[4]*B[0]) ; // Az = a*(Vx*By - Vy*Bx)
62
63 return;
64}
G4double B(G4double temperature)
double G4double
Definition: G4Types.hh:83
G4double FCof() const
Definition: G4Mag_EqRhs.hh:62

References B(), and G4Mag_EqRhs::FCof().

Referenced by G4ErrorMag_UsualEqRhs::EvaluateRhsGivenB().

◆ EvaluateRhsReturnB()

void G4EquationOfMotion::EvaluateRhsReturnB ( const G4double  y[],
G4double  dydx[],
G4double  Field[] 
) const
inlineinherited

◆ FCof()

G4double G4Mag_EqRhs::FCof ( ) const
inlineinherited

◆ GetFieldObj() [1/2]

G4Field * G4EquationOfMotion::GetFieldObj ( )
inlineinherited

◆ GetFieldObj() [2/2]

const G4Field * G4EquationOfMotion::GetFieldObj ( ) const
inlineinherited

Referenced by G4NystromRK4::GetField().

◆ GetFieldValue() [1/2]

template<class T_Field >
void G4TMagFieldEquation< T_Field >::GetFieldValue ( const G4double  Point[4],
G4double  Field[] 
) const
inline

Definition at line 56 of file G4TMagFieldEquation.hh.

58 {
59 itsField->T_Field::GetFieldValue(Point, Field);
60 }

References G4TMagFieldEquation< T_Field >::itsField.

◆ GetFieldValue() [2/2]

template<class T_Field >
G4TMagFieldEquation< T_Field >::GetFieldValue ( PositionAndTime  ,
Field   
)

◆ RightHandSide()

void G4EquationOfMotion::RightHandSide ( const G4double  y[],
G4double  dydx[] 
) const
inlineinherited

◆ SetChargeMomentumMass()

void G4Mag_UsualEqRhs::SetChargeMomentumMass ( G4ChargeState  particleCharge,
G4double  MomentumXc,
G4double  mass 
)
virtualinherited

Reimplemented from G4Mag_EqRhs.

Definition at line 67 of file G4Mag_UsualEqRhs.cc.

71{
72 G4Mag_EqRhs::SetChargeMomentumMass( particleCharge, MomentumXc, mass);
73}
virtual void SetChargeMomentumMass(G4ChargeState particleCharge, G4double MomentumXc, G4double mass)
Definition: G4Mag_EqRhs.cc:49

References G4Mag_EqRhs::SetChargeMomentumMass().

◆ SetFieldObj()

void G4EquationOfMotion::SetFieldObj ( G4Field pField)
inlineinherited

◆ TEvaluateRhsGivenB() [1/2]

template<class T_Field >
void G4TMagFieldEquation< T_Field >::TEvaluateRhsGivenB ( const G4double  y[],
const G4double  B[3],
G4double  dydx[] 
) const
inline

Definition at line 62 of file G4TMagFieldEquation.hh.

65 {
66 G4double momentum_mag_square = y[3]*y[3] + y[4]*y[4] + y[5]*y[5];
67 G4double inv_momentum_magnitude = 1.0 / std::sqrt( momentum_mag_square );
68 G4double cof = FCof()*inv_momentum_magnitude;
69
70 dydx[0] = y[3]*inv_momentum_magnitude; // (d/ds)x = Vx/V
71 dydx[1] = y[4]*inv_momentum_magnitude; // (d/ds)y = Vy/V
72 dydx[2] = y[5]*inv_momentum_magnitude; // (d/ds)z = Vz/V
73
74 dydx[3] = cof*(y[4]*B[2] - y[5]*B[1]) ; // Ax = a*(Vy*Bz - Vz*By)
75 dydx[4] = cof*(y[5]*B[0] - y[3]*B[2]) ; // Ay = a*(Vz*Bx - Vx*Bz)
76 dydx[5] = cof*(y[3]*B[1] - y[4]*B[0]) ; // Az = a*(Vx*By - Vy*Bx)
77
78 return ;
79 }

References B(), G4TMagFieldEquation< T_Field >::dydx, and G4Mag_EqRhs::FCof().

◆ TEvaluateRhsGivenB() [2/2]

template<class T_Field >
G4TMagFieldEquation< T_Field >::TEvaluateRhsGivenB ( ,
Field  ,
dydx   
)

Field Documentation

◆ dydx

template<class T_Field >
G4double G4TMagFieldEquation< T_Field >::dydx[]

◆ fCof_val

G4double G4Mag_EqRhs::fCof_val = 0.0
privateinherited

Definition at line 69 of file G4Mag_EqRhs.hh.

Referenced by G4Mag_EqRhs::FCof(), and G4Mag_EqRhs::SetChargeMomentumMass().

◆ fUnitConstant

const G4double G4Mag_EqRhs::fUnitConstant = 0.299792458 * (GeV/(tesla*m))
staticprivateinherited

Definition at line 71 of file G4Mag_EqRhs.hh.

◆ itsField

template<class T_Field >
T_Field* G4TMagFieldEquation< T_Field >::itsField
private

◆ PositionAndTime [1/2]

template<class T_Field >
G4TMagFieldEquation< T_Field >::PositionAndTime[3] = y[0]

Definition at line 87 of file G4TMagFieldEquation.hh.

◆ PositionAndTime [2/2]

template<class T_Field >
G4double G4TMagFieldEquation< T_Field >::PositionAndTime[4]

Definition at line 86 of file G4TMagFieldEquation.hh.


The documentation for this class was generated from the following file: