BDSIntegratorMag.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2023.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSINTEGRATORMAG_H
#define BDSINTEGRATORMAG_H
 
#include "BDSAuxiliaryNavigator.hh"
#include "BDSIntegratorDrift.hh"
#include "BDSMagnetStrength.hh"
 
#include "globals.hh" // geant4 types / globals
#include "G4MagIntegratorStepper.hh"
 
class G4Mag_EqRhs;
 
class BDSIntegratorMag: public G4MagIntegratorStepper, public BDSIntegratorDrift, public BDSAuxiliaryNavigator
{
public:
  BDSIntegratorMag(G4Mag_EqRhs* eqOfMIn,
                   G4int        nVariablesIn);
  
  virtual ~BDSIntegratorMag();
  BDSIntegratorMag& operator=(const BDSIntegratorMag&) = delete;
  BDSIntegratorMag(BDSIntegratorMag&) = delete;
 
  inline virtual G4double DistChord() const {return distChordPrivate;}
 
  inline virtual G4int IntegratorOrder() const {return 2;}
 
  static G4double thinElementLength;
 
  static G4double nominalMatrixRelativeMomCut;
 
  static G4bool currentTrackIsPrimary;
 
protected:
  void ConvertToGlobal(const G4ThreeVector &localPos,
                       const G4ThreeVector &localMom,
                       G4double             yOut[],
                       G4double             yErr[],
                       const G4double       momScaling = 1.0);
 
  inline void SetDistChord(G4double distChordIn) {distChordPrivate = distChordIn;}
  
  G4Mag_EqRhs* eqOfM;
  
  const G4int nVariables;
 
  G4MagIntegratorStepper* backupStepper;
 
  G4bool zeroStrength;
 
  G4double backupStepperMomLimit;
 
private:
  BDSIntegratorMag() = delete;
 
  G4double distChordPrivate;
};
 
#endif