TrkReco/TrkReco-00-08-59-p04/TrkReco/Bfield.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Bfield class
//
// 27-Mar-1999 : KUNIYA Toshio
//   Enabled QCS compornent, phi compornent
//
// 21-Feb-1999 : KUNIYA Toshio
//   Keeping comatibility, Bfield class is modified.
//   No longer fortran common block is used for bfield map.
//   Access functions are prepared for fortran call.
//
// 22-Feb-1999 : KUNIYA Toshio
//   constructor, destructor are go to protected entries.
// 
 
#if !defined(BFIELD_H)
#define BFIELD_H
 
#include "CLHEP/Vector/ThreeVector.h"
using CLHEP::Hep3Vector;
 
#ifndef CLHEP_POINT3D_H
#include "CLHEP/Geometry/Point3D.h"
#endif
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
    typedef HepGeom::Point3D<double> HepPoint3D;
#endif
using HepGeom::Point3D;
 
//class HepPoint3D;
 
//class Bfield;
 
class Bfield {
  
public:
  /// returns Bfield object.
  static Bfield * getBfield(int);
  
protected:
  /// Constructor, Destructor
  Bfield(int);
//zsl  virtual ~Bfield(void) {};
  
public://selectors
  
  ///returns B field Map ID
  int fieldID(void) const {return _fieldID;}
  
  ///returns B field
  const Hep3Vector & fieldMap(float x, float y, float z) const;
  const Hep3Vector & fieldMap(const HepPoint3D &xyz) const;
  void fieldMap(float *position, float *bfield);
  
  ///returns an element of B field
  float bx(float x, float y, float z) const;
  float by(float x, float y, float z) const;
  float bz(float x, float y, float z) const;
  float bx(const HepPoint3D &xyz) const;
  float by(const HepPoint3D &xyz) const;
  float bz(const HepPoint3D &xyz) const;
  
private:
  static  Bfield *_field[200];
  float _Bz[175][399];
  float _Br[175][399];
  float _Bphi[175][399];
  float _BzQR[101][163];
  float _BrQR[101][163];
  float _BphiQR[101][163];
  float _BzQL[17][51][52];
  float _BrQL[17][51][52];
  float _BphiQL[17][51][52];
  int   _fieldID;
  mutable float m_Bx;
  mutable float m_By;
  mutable float m_Bz;
  mutable float m_x;
  mutable float m_y;
  mutable float m_z;
  mutable Hep3Vector m_Bfld;
  
private:
  void updateCache(float, float, float) const;
};
 
#endif /* BFIELD_H */