MonopoleEquation.h

Go to the documentation of this file.

//
// =======================================================================
//
// class G4MonopoleEquation 
//
// Created:  30 April 2010, S. Burdin, B. Bozsogi 
//                       G4MonopoleEquation class for
//                       Geant4 extended example "monopole"
//
// Adopted for CMSSW by V.Ivanchenko 30 April 2018
// from Geant4 global tag geant4-10-04-ref-03                   
//
// =======================================================================
//
// Class description:
//
// This is the right-hand side of equation of motion in a 
// magnetic field only for magnetic monopoles.
//
// -------------------------------------------------------------------

#ifndef SimG4Core_MagneticField_MonopoleEquation_h
#define SimG4Core_MagneticField_MonopoleEquation_h 1

#include "G4EquationOfMotion.hh"
#include "G4MagneticField.hh"

class MonopoleEquation : public G4EquationOfMotion
{
public:  // with description

  MonopoleEquation(G4MagneticField *emField );

  ~MonopoleEquation() override;

  void  SetChargeMomentumMass( G4ChargeState particleChargeState,
                               G4double      momentum, 
                               G4double      mass) override;
  // magnetic charge in e+ units
                                 
  void EvaluateRhsGivenB(const G4double y[],
                         const G4double Field[],
                         G4double dydx[] ) const override;
  // Given the value of the electromagnetic field, this function 
  // calculates the value of the derivative dydx.

private:

  G4double  fMagCharge ;
  G4double  fElCharge;
  G4double  fMassCof;
};

#endif