G4KM_OpticalEqRhs Class Reference

#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)
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_OpticalEqRhs.hh.

Constructor & Destructor Documentation

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

inline

Definition at line 69 of file G4KM_OpticalEqRhs.hh.

{ }

Member Function Documentation

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.

References python.hepunit::c_light, 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]);
  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;
}

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

virtual

Reimplemented from G4Mag_EqRhs.

Definition at line 94 of file G4KM_OpticalEqRhs.cc.

{ 
}

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

Definition at line 54 of file G4KM_OpticalEqRhs.cc.

References G4InuclSpecialFunctions::bindingEnergy(), G4lrint(), G4NucleiProperties::GetBindingEnergy(), G4V3DNucleus::GetCharge(), G4V3DNucleus::GetMassNumber(), python.hepunit::hbarc, python.hepunit::neutron_mass_c2, python.hepunit::pi, and python.hepunit::proton_mass_c2.

Referenced by G4RKPropagation::Init().

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

---

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

• G4KM_OpticalEqRhs.hh
• G4KM_OpticalEqRhs.cc