GFField.cxx

Go to the documentation of this file.

#include "GFField.h"
 
#include <iostream>
 
#include <CLHEP/Units/SystemOfUnits.h>
 
namespace SQGenFit
{
GFField::GFField(const PHField* field): 
  _field(field),
  _scale(1.),
  _offset(0.) 
{
  _scale = 1.;
  _disable = false;
}
 
TVector3 GFField::get(const TVector3& pos) const
{
  double x = pos.X();
  double y = pos.Y();
  double z = pos.Z();
 
  double Bx, By, Bz;
  get(x, y, z, Bx, By, Bz);
 
  return TVector3(Bx, By, Bz);
}
 
void GFField::get(const double& x, const double& y, const double& z, double& Bx, double& By, double& Bz) const
{
  if(_disable)
  {
    Bx = 0.;
    By = 0.;
    Bz = 0.;
    return;
  }
 
  const double Point[] = {x*CLHEP::cm, y*CLHEP::cm, (z + _offset)*CLHEP::cm, 0.};
  double Bfield[6];
  for(int i = 0; i < 6; ++i)
  {
    Bfield[i] = 0.;
  }
 
  _field->GetFieldValue(Point, Bfield);
  Bx = _scale*Bfield[0]/CLHEP::kilogauss;
  By = _scale*Bfield[1]/CLHEP::kilogauss;
  Bz = _scale*Bfield[2]/CLHEP::kilogauss;
}
 
}