G4KM_NucleonEqRhs Class Reference

#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 (G4ChargeState particleCharge, G4double MomentumXc, G4double MassXc2)
void	SetMass (G4double aMass)
Public Member Functions inherited from G4Mag_EqRhs
	G4Mag_EqRhs (G4MagneticField *magField)
virtual	~G4Mag_EqRhs ()
G4double	FCof () const
Public Member Functions inherited from G4EquationOfMotion
	G4EquationOfMotion (G4Field *Field)
virtual	~G4EquationOfMotion ()
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_DummyField *	field,
		G4V3DNucleus *	nucleus
	)

Definition at line 44 of file G4KM_NucleonEqRhs.cc.

References G4V3DNucleus::GetMassNumber(), python.hepunit::hbarc, and python.hepunit::pi2.

                                             :
  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 ( )

inline

Definition at line 81 of file G4KM_NucleonEqRhs.hh.

{ }

Member Function Documentation

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.

References python.hepunit::c_light, density, G4VNuclearDensity::GetDensity(), G4VNuclearDensity::GetDeriv(), and G4V3DNucleus::GetNuclearDensity().

{
  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;
 */
}

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

virtual

Reimplemented from G4Mag_EqRhs.

Definition at line 111 of file G4KM_NucleonEqRhs.cc.

{ 
}

void G4KM_NucleonEqRhs::SetMass ( G4double  aMass )

inline

Definition at line 74 of file G4KM_NucleonEqRhs.hh.

Referenced by G4RKPropagation::Init().

{
  theMass = aMass;
}

---

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

• G4KM_NucleonEqRhs.hh
• G4KM_NucleonEqRhs.cc