SimG4MagneticFieldTool.h

Go to the documentation of this file.

#ifndef SIMG4COMPONENTS_G4MAGNETICFIELDTOOL_H
#define SIMG4COMPONENTS_G4MAGNETICFIELDTOOL_H
 
// Gaudi
#include "GaudiKernel/AlgTool.h"
 
// k4FWCore
#include "k4Interface/ISimG4MagneticFieldTool.h"
#include "k4Interface/IGeoSvc.h"
 
// Geant4
#include "G4SystemOfUnits.hh"
#include "G4MagneticField.hh"
 
// Forward declarations:
// Geant4 classes
class G4MagIntegratorStepper;
 
 
class SimG4MagneticFieldTool : public AlgTool, virtual public ISimG4MagneticFieldTool {
public:
  SimG4MagneticFieldTool(const std::string& type,
                         const std::string& name,
                         const IInterface* parent);
 
  virtual ~SimG4MagneticFieldTool();
 
  virtual StatusCode initialize() final;
 
  virtual StatusCode finalize() final;
 
  virtual const G4MagneticField* field() const final;
 
  G4MagIntegratorStepper* stepper(const std::string&, G4MagneticField*) const;
 
private:
  ServiceHandle<IGeoSvc> m_geoSvc;
 
  G4MagneticField* m_field = nullptr;
 
  Gaudi::Property<bool> m_fieldOn{this, "FieldOn", false, "Switch to turn field off"};
 
  Gaudi::Property<double> m_minEps{this, "MinimumEpsilon", 0, "Minimum epsilon (see G4 documentation)"};
 
  Gaudi::Property<double> m_maxEps{this, "MaximumEpsilon", 0, "Maximum epsilon (see G4 documentation)"};
 
  Gaudi::Property<double> m_deltaChord{this, "DeltaChord", 0, "Missing distance for the chord finder"};
 
  Gaudi::Property<double> m_deltaOneStep{this, "DeltaOneStep", 0, "Delta(one-step)"};
 
  Gaudi::Property<double> m_maxStep{this, "MaximumStep", 1. * m, "Maximum step length in field (see G4 documentation)"};
 
  Gaudi::Property<double> m_minStep{this, "MinimumStep", 0.01 * mm, "Minimum step length in field (see G4 documentation)"};
 
  Gaudi::Property<std::string> m_integratorStepper{this, "IntegratorStepper", "NystromRK4", "Integrator stepper name"};
};
 
#endif /* SIMG4COMPONENTS_G4MAGNETICFIELDTOOL_H */