G4KM_OpticalEqRhs Class Reference

#include <G4KM_OpticalEqRhs.hh>

Inheritance diagram for G4KM_OpticalEqRhs:

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_OpticalEqRhs (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	SetFactor (G4double mass, G4double opticalParameter)
void	SetFieldObj (G4Field *pField)
	~G4KM_OpticalEqRhs ()

Private Attributes
G4double	fCof_val = 0.0
G4Field *	itsField = nullptr
G4double	theFactor
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_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;
}

References theFactor, and theMass.

~G4KM_OpticalEqRhs()

G4KM_OpticalEqRhs::~G4KM_OpticalEqRhs ( )

inline

Definition at line 69 of file G4KM_OpticalEqRhs.hh.

{ }

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

References source.hepunit::c_light, G4VNuclearDensity::GetDeriv(), G4V3DNucleus::GetNuclearDensity(), pos, theFactor, theMass, and theNucleus.

EvaluateRhsReturnB()

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

inlineinherited

Referenced by G4RKG3_Stepper::StepNoErr().

FCof()

G4double G4Mag_EqRhs::FCof ( ) const

inlineinherited

Definition at line 62 of file G4Mag_EqRhs.hh.

{ return fCof_val; }

References G4Mag_EqRhs::fCof_val.

Referenced by G4MagHelicalStepper::AdvanceHelix(), G4BFieldIntegrationDriver::CurvatureRadius(), G4ErrorMag_UsualEqRhs::EvaluateRhsGivenB(), G4Mag_SpinEqRhs::EvaluateRhsGivenB(), G4Mag_UsualEqRhs::EvaluateRhsGivenB(), G4ConstRK4::RightHandSideConst(), and G4TMagFieldEquation< T_Field >::TEvaluateRhsGivenB().

GetFieldObj() [1/2]

G4Field * G4EquationOfMotion::GetFieldObj ( )

inlineinherited

GetFieldObj() [2/2]

const G4Field * G4EquationOfMotion::GetFieldObj ( ) const

inlineinherited

Referenced by G4NystromRK4::GetField().

GetFieldValue()

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

inlineinherited

Referenced by G4BFieldIntegrationDriver::GetFieldValue(), and G4MultiLevelLocator::ReportFieldValue().

RightHandSide()

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

inlineinherited

Referenced by G4ModifiedMidpoint::DoStep(), and G4VFSALIntegrationStepper::RightHandSide().

SetChargeMomentumMass()

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

virtual

Reimplemented from G4Mag_EqRhs.

Definition at line 94 of file G4KM_OpticalEqRhs.cc.

{ 
}

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

References A, G4InuclSpecialFunctions::bindingEnergy(), G4lrint(), G4NucleiProperties::GetBindingEnergy(), G4V3DNucleus::GetCharge(), G4V3DNucleus::GetMassNumber(), source.hepunit::hbarc, source.hepunit::neutron_mass_c2, pi, source.hepunit::proton_mass_c2, theFactor, theMass, theNucleus, and Z.

Referenced by G4RKPropagation::Init().

SetFieldObj()

void G4EquationOfMotion::SetFieldObj ( G4Field *  pField )

inlineinherited

Referenced by G4FieldManager::SetDetectorField(), and G4NystromRK4::SetDistanceForConstantField().

Field Documentation

fCof_val

G4double G4Mag_EqRhs::fCof_val = 0.0

privateinherited

Definition at line 69 of file G4Mag_EqRhs.hh.

Referenced by G4Mag_EqRhs::FCof(), and G4Mag_EqRhs::SetChargeMomentumMass().

fUnitConstant

const G4double G4Mag_EqRhs::fUnitConstant = 0.299792458 * (GeV/(tesla*m))

staticprivateinherited

Definition at line 71 of file G4Mag_EqRhs.hh.

itsField

G4Field* G4EquationOfMotion::itsField = nullptr

privateinherited

Definition at line 90 of file G4EquationOfMotion.hh.

theFactor

G4double G4KM_OpticalEqRhs::theFactor

private

Definition at line 65 of file G4KM_OpticalEqRhs.hh.

Referenced by EvaluateRhsGivenB(), G4KM_OpticalEqRhs(), and SetFactor().

theMass

G4double G4KM_OpticalEqRhs::theMass

private

Definition at line 66 of file G4KM_OpticalEqRhs.hh.

Referenced by EvaluateRhsGivenB(), G4KM_OpticalEqRhs(), and SetFactor().

theNucleus

G4V3DNucleus* G4KM_OpticalEqRhs::theNucleus

private

Definition at line 64 of file G4KM_OpticalEqRhs.hh.

Referenced by EvaluateRhsGivenB(), and SetFactor().

---

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

• source/processes/hadronic/models/binary_cascade/include/G4KM_OpticalEqRhs.hh
• source/processes/hadronic/models/binary_cascade/src/G4KM_OpticalEqRhs.cc