SimG4MagneticFieldFromMapTool.h

Go to the documentation of this file.

#ifndef SIMG4COMPONENTS_G4MAGNETICFIELDFROMMAPTOOL_H
#define SIMG4COMPONENTS_G4MAGNETICFIELDFROMMAPTOOL_H
 
// Gaudi
#include "GaudiKernel/AlgTool.h"
 
// FCCSW
#include "SimG4Interface/ISimG4MagneticFieldTool.h"
 
// Geant4
#include "G4SystemOfUnits.hh"
#include "G4MagneticField.hh"
 
// Forward declarations:
// Geant 4 classes
class G4MagIntegratorStepper;
 
// FCCSW
/*
namespace sim {
  class MapField;
}
*/
 
class SimG4MagneticFieldFromMapTool : public AlgTool, virtual public ISimG4MagneticFieldTool {
public:
  SimG4MagneticFieldFromMapTool(const std::string& type, const std::string& name, const IInterface* parent);
 
  virtual ~SimG4MagneticFieldFromMapTool();
 
  virtual StatusCode initialize() final;
 
  virtual StatusCode finalize() final;
 
  virtual const G4MagneticField* field() const final;
 
  G4MagIntegratorStepper* stepper(const std::string&, G4MagneticField*) const;
 
private:
  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"};
  Gaudi::Property<std::string> m_mapFilePath{this, "MapFile", "", "Path to file containing fieldmap"};
  Gaudi::Property<double> m_addFieldBz{this, "AddFieldBz", 0., "Additional constant field, z component (default: 0.)"};
  Gaudi::Property<double> m_addFieldMaxR{this, "AddFieldMaxR", -1., "Maximum radius of additional constant field (default: no limit)"};
  Gaudi::Property<double> m_addFieldMaxZ{this, "AddFieldMaxZ", -1., "Maximum z coordinate of additional constant field (default: no limit)"};
  Gaudi::Property<double> m_fieldMaxR{this, "FieldMaxR", -1., "Field maximum radius (default: no limit)"};
  Gaudi::Property<double> m_fieldMaxZ{this, "FieldMaxZ", -1., "Field maximum z coordinate (default: no limit)"};
 
  StatusCode loadRootMap();
  StatusCode loadComsolMap();
};
 
#endif