#include <G4KM_NucleonEqRhs.hh>

Inheritance diagram for G4KM_NucleonEqRhs:

Public Member Functions
virtual 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

	G4KM_NucleonEqRhs (G4KM_DummyField field, G4V3DNucleus nucleus)

G4Field *	GetFieldObj ()

const G4Field *	GetFieldObj () const

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

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

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

void	SetFieldObj (G4Field *pField)

void	SetMass (G4double aMass)

	~G4KM_NucleonEqRhs ()

Private Attributes
G4int	A

G4double	factor

G4double	fCof_val = 0.0

G4Field *	itsField = nullptr

G4double	theMass

G4V3DNucleus *	theNucleus

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

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 45 of file G4KM_NucleonEqRhs.cc.

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

References A, G4Pow::A23(), factor, G4Pow::GetInstance(), G4V3DNucleus::GetMassNumber(), source.hepunit::hbarc, pi2, theMass, and theNucleus.

◆ ~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 55 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)/
            G4Pow::GetInstance()->A13(density)*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;
 */
}