/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkExtAlg/TrkExtAlg-00-00-64/src/ExtBesMagneticField.cxx

Go to the documentation of this file.

#include "TrkExtAlg/ExtBesMagneticField.h"

#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SvcFactory.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/Bootstrap.h"

#include "G4FieldManager.hh"
#include "G4MagneticField.hh"
#include "globals.hh"

#include "CLHEP/Geometry/Vector3D.h"
#include "CLHEP/Geometry/Point3D.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();
}

---