MucMagneticField.h

Go to the documentation of this file.

#ifndef MUCMAGNETICFIELD_H
#define MUCMAGNETICFIELD_H 1
 
#include <string>
 
#include "CLHEP/Geometry/Point3D.h"
#include "CLHEP/Geometry/Vector3D.h"
 
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
//  backwards compatibility will be enabled ONLY in CLHEP 1.9
typedef HepGeom::Point3D<double>  HepPoint3D;
typedef HepGeom::Vector3D<double> HepVector3D;
#endif
 
class MucMagneticField {
public:
#ifndef BEAN
  MucMagneticField();
#else
  MucMagneticField( const std::string Path );
#endif
  ~MucMagneticField();
 
  void getMucField( int part, int layer, int mat, HepPoint3D& r, HepVector3D& b );
  void readPar();
 
  std::string getPath() { return path; }
 
private:
  static const int npar = 18;
  double bipx0[npar], bipy0[npar], bipz0[npar]; // barrel iron par bx,by,bz in each layer
  double bipx1[npar], bipy1[npar], bipz1[npar];
  double bipx2[npar], bipy2[npar], bipz2[npar];
  double bipx3[npar], bipy3[npar], bipz3[npar];
  double bipx4[npar], bipy4[npar], bipz4[npar];
  double bipx5[npar], bipy5[npar], bipz5[npar];
  double bipx6[npar], bipy6[npar], bipz6[npar];
  double bipx7[npar], bipy7[npar], bipz7[npar];
  double bipx8[npar], bipy8[npar], bipz8[npar];
 
  double bapx0[npar], bapy0[npar], bapz0[npar]; // barrel air par bx,by,bz in each layer
  double bapx1[npar], bapy1[npar], bapz1[npar];
  double bapx2[npar], bapy2[npar], bapz2[npar];
  double bapx3[npar], bapy3[npar], bapz3[npar];
  double bapx4[npar], bapy4[npar], bapz4[npar];
  double bapx5[npar], bapy5[npar], bapz5[npar];
  double bapx6[npar], bapy6[npar], bapz6[npar];
  double bapx7[npar], bapy7[npar], bapz7[npar];
 
  static const int npar1 = 17;
  double aipx0[npar1], aipy0[npar1], aipz0[npar1]; // endcap iron par bx,by,bz in each layer
  double aipx1[npar1], aipy1[npar1], aipz1[npar1];
  double aipx2[npar1], aipy2[npar1], aipz2[npar1];
  double aipx3[npar1], aipy3[npar1], aipz3[npar1];
  double aipx4[npar1], aipy4[npar1], aipz4[npar1];
  double aipx5[npar1], aipy5[npar1], aipz5[npar1];
  double aipx6[npar1], aipy6[npar1], aipz6[npar1];
  double aipx7[npar1], aipy7[npar1], aipz7[npar1];
  double aipx8[npar1], aipy8[npar1], aipz8[npar1];
 
  double aapx0[npar1], aapy0[npar1], aapz0[npar1]; // endcap air par bx,by,bz in each layer
  double aapx1[npar1], aapy1[npar1], aapz1[npar1];
  double aapx2[npar1], aapy2[npar1], aapz2[npar1];
  double aapx3[npar1], aapy3[npar1], aapz3[npar1];
  double aapx4[npar1], aapy4[npar1], aapz4[npar1];
  double aapx5[npar1], aapy5[npar1], aapz5[npar1];
  double aapx6[npar1], aapy6[npar1], aapz6[npar1];
  double aapx7[npar1], aapy7[npar1], aapz7[npar1];
 
  double bp[npar];  // barrel par
  double ep[npar1]; // endcap par
 
  std::string filename;
  std::string path; // path to data files
};
#endif