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/ISimG4Svc.h"
#include "k4Interface/ISimG4MagneticFieldTool.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 */