G4FieldManagerToolBase.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef G4ATLASTOOLS_G4FieldManagerToolBase_H
#define G4ATLASTOOLS_G4FieldManagerToolBase_H
 
// Framework includes
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
#include "AthenaBaseComps/AthAlgTool.h"
 
// G4Atlas includes
#include "G4AtlasInterfaces/IEquationOfMotionTool.h"
#include "G4AtlasInterfaces/IFieldManagerTool.h"
#include "G4AtlasInterfaces/IG4FieldSvc.h"
 
// Geant4 includes
#include "G4Version.hh"
 
// Forward declarations
class G4VIntegrationDriver;
class G4MagIntegratorStepper;
 

class G4FieldManagerToolBase : public extends<AthAlgTool, IFieldManagerTool>
{
 
  public:
    G4FieldManagerToolBase(const std::string& type, const std::string& name,
                            const IInterface* parent);

    virtual ~G4FieldManagerToolBase(){}

    virtual StatusCode initialize( ) override;
 
  protected:
 
#if G4VERSION_NUMBER >= 1040
    G4VIntegrationDriver*
    createDriverAndStepper(const std::string& stepperType, G4MagneticField* field) const;
#endif
 
#if G4VERSION_NUMBER < 1040
    G4MagIntegratorStepper*
    getStepper(const std::string& stepperType, G4MagneticField* field) const;
#endif

    StatusCode setFieldParameters(G4FieldManager* fieldMgr) const;

    ServiceHandle<IG4FieldSvc> m_fieldSvc{this, "FieldSvc", "G4FieldSvc", "Service providing a G4MagneticField"};

    Gaudi::Property<bool> m_fieldOn{this, "FieldOn", true, "Toggles field on/off"};

    Gaudi::Property<std::string> m_integratorStepper{this, "IntegratorStepper", "AtlasRK4", "Integrator stepper name"};

    ToolHandle<IEquationOfMotionTool> m_equationOfMotion{this, "EquationOfMotion", "", ""};
    //std::string m_equationOfMotion;

    Gaudi::Property<double> m_minEps{this, "MinimumEpsilonStep", -1.0, "Minimum epsilon (see G4 documentation)"};
    Gaudi::Property<double> m_maxEps{this, "MaximumEpsilonStep", -1.0, "Maximum epsilon (see G4 documentation)"};
    Gaudi::Property<double> m_deltaChord{this, "DeltaChord", -1.0, "Missing distance for the chord finder"};
    Gaudi::Property<double> m_deltaOneStep{this, "DeltaOneStep", -1.0, "Delta(one-step)"};
    Gaudi::Property<double> m_deltaIntersection{this, "DeltaIntersection", -1.0, "Accuracy for boundary intersection"};
    Gaudi::Property<double> m_maxStep{this, "MaximumStep", -1.0, "Maximum step length in field (see G4 documentation)"};
    Gaudi::Property<double> m_minStep{this, "MinimumStep",1e-2, "Minimum step length in field (see G4 documentation)"};
 
};
 
#endif // G4ATLASTOOLS_G4FieldManagerToolBase_H