MagFieldScanner.h

Go to the documentation of this file.

#ifndef MAGFIELDSCANNER_H
#define MAGFIELDSCANNER_H
 
// Gaudi
#include "GaudiKernel/Service.h"
 
// k4FWCore
#include "k4Interface/IGeoSvc.h"
#include "k4Interface/ISimG4MagneticFieldTool.h"
#include "k4Interface/ISimG4Svc.h"
 
class MagFieldScanner : public Service {
public:
  explicit MagFieldScanner(const std::string& name, ISvcLocator* svcLoc);
 
  virtual StatusCode initialize();
  virtual StatusCode finalize();
  virtual ~MagFieldScanner() {};
 
private:
  ServiceHandle<IGeoSvc> m_geoSvc;
 
  ServiceHandle<ISimG4Svc> m_simG4Svc;
 
  Gaudi::Property<std::string> m_outFilePath{this, "outFilePath", "magFieldProbes.root", "Output file path"};
 
  Gaudi::Property<std::vector<std::vector<double>>> m_xyPlaneProbes{
      this, "xyPlaneProbes", {}, "xy-plane probe definitions"};
 
  Gaudi::Property<std::vector<std::vector<double>>> m_zPlaneProbes{
      this, "zPlaneProbes", {}, "z-plane probe definitions"};
 
  Gaudi::Property<std::vector<std::vector<double>>> m_tubeProbes{this, "tubeProbes", {}, "Tube probe definitions"};
 
  struct XYPlaneProbe {
    const double xMax;
    const double yMax;
    const double z;
  };
 
  struct ZPlaneProbe {
    const double zMin;
    const double zMax;
    const double rMax;
    const double phi;
  };
 
  struct TubeProbe {
    const double zMin;
    const double zMax;
    const double r;
  };
 
  friend std::ostream& operator<<(std::ostream& outStream, const XYPlaneProbe& probe);
  friend std::ostream& operator<<(std::ostream& outStream, const ZPlaneProbe& probe);
  friend std::ostream& operator<<(std::ostream& outStream, const TubeProbe& probe);
};
 
#endif /* MAGFIELDSCANNER_H */