#include <G4KM_OpticalEqRhs.hh>

Inheritance diagram for G4KM_OpticalEqRhs:

Public Member Functions
	G4KM_OpticalEqRhs (G4KM_DummyField field, G4V3DNucleus nucleus)

	~G4KM_OpticalEqRhs ()

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

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

void	SetFactor (G4double mass, G4double opticalParameter)

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

	~G4Mag_EqRhs () override

G4double	FCof () const

void	SetChargeMomentumMass (G4ChargeState particleCharge, G4double MomentumXc, G4double mass) override

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

G4Field *	GetFieldObj ()

void	SetFieldObj (G4Field *pField)

Detailed Description

Definition at line 49 of file G4KM_OpticalEqRhs.hh.

Constructor & Destructor Documentation

◆ G4KM_OpticalEqRhs()

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

Definition at line 45 of file G4KM_OpticalEqRhs.cc.

                                             :
  G4Mag_EqRhs(field), theNucleus(nucleus)
{
  theFactor = 0;
  theMass = 0;
}

◆ ~G4KM_OpticalEqRhs()

G4KM_OpticalEqRhs::~G4KM_OpticalEqRhs ( )

inline

Definition at line 69 of file G4KM_OpticalEqRhs.hh.

70{ }

Member Function Documentation

◆ EvaluateRhsGivenB()

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

virtual

Implements G4Mag_EqRhs.

Definition at line 71 of file G4KM_OpticalEqRhs.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]);
  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;   // 
 
// V=K*rho(r) ==> dydx[3] = -dV/dr*dr/dx = -K*d(rho)/dr*dr/dx.
// Idem for dydx[4] and dydx[5]
 
  const G4VNuclearDensity * nuclearDensity=theNucleus->GetNuclearDensity();
 
  G4ThreeVector pos(y[0],y[1],y[2]);
  G4double deriv = theFactor*nuclearDensity->GetDeriv(pos);
 
  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;
}

◆ SetChargeMomentumMass()

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

virtual

Implements G4EquationOfMotion.

Definition at line 94 of file G4KM_OpticalEqRhs.cc.

95{

96}

◆ SetFactor()

void G4KM_OpticalEqRhs::SetFactor	(	G4double	mass,
		G4double	opticalParameter )

Definition at line 54 of file G4KM_OpticalEqRhs.cc.

{
  G4double A = theNucleus->GetMassNumber();
  G4double Z = theNucleus->GetCharge();
  G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(G4lrint(A), G4lrint(Z));
  G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
  G4double reducedMass = mass*nucleusMass/(mass+nucleusMass);
 
  G4double nucleonMass = (proton_mass_c2+neutron_mass_c2)/2;
 
// _factor in (MeV*fermi)*fermi/MeV = fermi*fermi  -- need to have A as density normalized to 1
  theFactor = 2*pi*hbarc*hbarc*(1+mass/nucleonMass)* opticalParameter/reducedMass * A;
 
  theMass = mass;
}

Referenced by G4RKPropagation::Init().

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

Public Member Functions