SimG4ConstantMagneticFieldTool.h

Go to the documentation of this file.

#ifndef SIMG4COMPONENTS_G4CONSTANTMAGNETICFIELDTOOL_H
#define SIMG4COMPONENTS_G4CONSTANTMAGNETICFIELDTOOL_H
 
// Gaudi
#include "GaudiKernel/AlgTool.h"
 
// FCCSW
#include "SimG4Interface/ISimG4MagneticFieldTool.h"
 
// Geant4
#include "G4SystemOfUnits.hh"
 
// Forward declarations:
// Geant 4 classes
class G4MagIntegratorStepper;
 
// FCCSW
namespace sim {
class ConstantField;
}
 
class SimG4ConstantMagneticFieldTool : public AlgTool, virtual public ISimG4MagneticFieldTool {
public:
  SimG4ConstantMagneticFieldTool(const std::string& type, const std::string& name, const IInterface* parent);
 
  virtual ~SimG4ConstantMagneticFieldTool();
 
  virtual StatusCode initialize() final;
 
  virtual StatusCode finalize() final;
 
  virtual const G4MagneticField* field() const final;
 
  G4MagIntegratorStepper* stepper(const std::string&, G4MagneticField*) const;
 
private:
  sim::ConstantField* m_field;
  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"};
 
  Gaudi::Property<double> m_fieldComponentX{this, "FieldComponentX", 0, "Field X component"};
  Gaudi::Property<double> m_fieldComponentY{this, "FieldComponentY", 0, "Field Y component"};
  Gaudi::Property<double> m_fieldComponentZ{this, "FieldComponentZ", -4 * tesla, "Field Z component"};
  Gaudi::Property<double> m_fieldRadMax{this, "FieldRMax", 6 * m, "Field max radius"};
  Gaudi::Property<double> m_fieldZMax{this, "FieldZMax", 20. * m, "Field max Z"};
};
 
#endif