ExtBesMagneticField.cxx

Go to the documentation of this file.

#include "ExtBesMagneticField.h"
 
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
 
#include "G4FieldManager.hh"
#include "G4MagneticField.hh"
#include "globals.hh"
 
#include "CLHEP/Geometry/Point3D.h"
#include "CLHEP/Geometry/Vector3D.h"
#include "CLHEP/Units/PhysicalConstants.h"
 
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Vector3D<double> HepVector3D;
#endif
 
using namespace CLHEP;
 
ExtBesMagneticField::ExtBesMagneticField() {
  GetFieldManager()->SetDetectorField( this );
  GetFieldManager()->CreateChordFinder( this );
  //    GetFieldManager()->SetDeltaOneStep(0.5*mm);
  //    GetFieldManager()->SetDeltaIntersection(0.5*mm);
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  StatusCode   sc         = svcLocator->service( "MagneticFieldSvc", m_pIMF );
  if ( sc != StatusCode::SUCCESS )
  { G4cout << "Unable to open Magnetic field service" << G4endl; }
}
 
ExtBesMagneticField::~ExtBesMagneticField() {}
 
void ExtBesMagneticField::GetFieldValue( const double Point[3], double* Bfield ) const {
  double x = Point[0];
  double y = Point[1];
  double z = Point[2];
 
  HepPoint3D  r( x, y, z );
  HepVector3D b;
 
  if ( m_pIMF->ifRealField() == true ) m_pIMF->fieldVector( r, b );
  else m_pIMF->uniFieldVector( r, b );
 
  Bfield[0] = b.x();
  Bfield[1] = b.y();
  Bfield[2] = b.z();
}
 
#include "G4TransportationManager.hh"
 
G4FieldManager* ExtBesMagneticField::GetFieldManager() {
  return G4TransportationManager::GetTransportationManager()->GetFieldManager();
}