#include <G4KM_NucleonEqRhs.hh>

Inheritance diagram for G4KM_NucleonEqRhs:

Public Member Functions
	G4KM_NucleonEqRhs (G4KM_DummyField field, G4V3DNucleus nucleus)

	~G4KM_NucleonEqRhs ()

virtual void	EvaluateRhsGivenB (const G4double y[], const G4double B[3], G4double dydx[]) const

virtual void	SetChargeMomentumMass (G4double particleCharge, G4double MomentumXc, G4double MassXc2)

void	SetMass (G4double aMass)

Public Member Functions inherited from G4Mag_EqRhs
	G4Mag_EqRhs (G4MagneticField *magField)

virtual	~G4Mag_EqRhs ()

virtual void	EvaluateRhsGivenB (const G4double y[], const G4double B[3], G4double dydx[]) const =0

G4double	FCof () const

virtual void	SetChargeMomentumMass (G4double particleCharge, G4double MomentumXc, G4double mass)

Public Member Functions inherited from G4EquationOfMotion
	G4EquationOfMotion (G4Field *Field)

virtual	~G4EquationOfMotion ()

virtual void	EvaluateRhsGivenB (const G4double y[], const G4double B[3], G4double dydx[]) const =0

virtual void	SetChargeMomentumMass (G4double particleCharge, G4double MomentumXc, G4double MassXc2)=0

void	RightHandSide (const G4double y[], G4double dydx[]) const

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

void	GetFieldValue (const G4double Point[4], G4double Field[]) const

const G4Field *	GetFieldObj () const

void	SetFieldObj (G4Field *pField)

Detailed Description

Definition at line 49 of file G4KM_NucleonEqRhs.hh.

Constructor & Destructor Documentation

◆ G4KM_NucleonEqRhs()

G4KM_NucleonEqRhs::G4KM_NucleonEqRhs	(	G4KM_DummyField *	field,
		G4V3DNucleus *	nucleus
	)

Definition at line 44 of file G4KM_NucleonEqRhs.cc.

                                             :
  G4Mag_EqRhs(field), theNucleus(nucleus)
{
  theMass = 0.;
  A = theNucleus->GetMassNumber();
  factor = hbarc*hbarc*std::pow(3.*pi2*A,2./3.)/3.;
}

◆ ~G4KM_NucleonEqRhs()

G4KM_NucleonEqRhs::~G4KM_NucleonEqRhs ( )

inline

Definition at line 81 of file G4KM_NucleonEqRhs.hh.

82{ }

Member Function Documentation

◆ EvaluateRhsGivenB()

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

virtual

Implements G4Mag_EqRhs.

Definition at line 54 of file G4KM_NucleonEqRhs.cc.

{
  G4double yMod = std::sqrt(y[0]*y[0]+y[1]*y[1]+y[2]*y[2]);
  G4double e = std::sqrt(theMass*theMass+y[3]*y[3]+y[4]*y[4]+y[5]*y[5]);
 
// y[0..2] is position
// y[3..5] is momentum (and not mom.direction)
    
  dydx[0] = c_light*y[3]/e;   //
  dydx[1] = c_light*y[4]/e;   //  dq/dt=dH/dp = c*p/e
  dydx[2] = c_light*y[5]/e;   // 
 
/*
 * // debug
 *  G4cout << " Nucleon RHS : 0..2(dpos/dt) " << 
 *       dydx[0] << " " <<
 *       dydx[1] << " " <<
 *       dydx[2] << " " << G4endl;
 */
 
      
// V=K*rho(r) ==> dydx[3] = -dV/dr*dr/dx = -K*d(rho)/dr*dr/dx.
// GF should be V=K*rho(r) ==> dydx[3] = -dV/dr*dr/dx = -K*d(rho)/dr*dr/dt
// GF  and dV/dt = dE/dt ==> dp/dt = dE/dt * dp/dE = dE/dt *e/p
// Idem for dydx[4] and dydx[5]
 
  G4ThreeVector pos(y[0],y[1],y[2]);
 
  const G4VNuclearDensity * nuclearDensity=theNucleus->GetNuclearDensity();
 
// do not check for theMass != 0 : it is an error and core dump will signal it
 
  G4double density=  nuclearDensity->GetDensity(pos);
  G4double deriv(0);
  if (density > 0 ) deriv = (factor/theMass)*
            std::pow(density, -1./3.)*nuclearDensity->GetDeriv(pos);
 
//  dydx[3] = yMod == 0 ? 0 : -deriv*y[0]/yMod;
//  dydx[4] = yMod == 0 ? 0 : -deriv*y[1]/yMod;
//  dydx[5] = yMod == 0 ? 0 : -deriv*y[2]/yMod;
  dydx[3] = yMod == 0 ? 0 : deriv*y[0]/yMod*c_light;
  dydx[4] = yMod == 0 ? 0 : deriv*y[1]/yMod*c_light;
  dydx[5] = yMod == 0 ? 0 : deriv*y[2]/yMod*c_light;
 
 
/*
 * // debug
 * G4cout << " Nucleon RHS : 3..5(dE/dt) " << 
 *       dydx[3] << " " <<
 *       dydx[4] << " " <<
 *       dydx[5] << " " << G4endl;
 */
}

◆ SetChargeMomentumMass()

void G4KM_NucleonEqRhs::SetChargeMomentumMass	(	G4double	particleCharge,
		G4double	MomentumXc,
		G4double	MassXc2
	)

inlinevirtual

Reimplemented from G4Mag_EqRhs.

Definition at line 87 of file G4KM_NucleonEqRhs.hh.

88{ }

◆ SetMass()

void G4KM_NucleonEqRhs::SetMass ( G4double aMass )

inline

Definition at line 74 of file G4KM_NucleonEqRhs.hh.

{
  theMass = aMass;
}

Referenced by G4RKPropagation::Init().

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

Public Member Functions