DMagneticFieldMapCalibDB.h

Go to the documentation of this file.

// $Id$
//
//    File: DMagneticFieldMapCalibDB.h
// Created: Thu Jul 19 13:58:21 EDT 2007
// Creator: davidl (on Darwin fwing-dhcp61.jlab.org 8.10.1 i386)
//

#ifndef _DMagneticFieldMapCalibDB_
#define _DMagneticFieldMapCalibDB_

#include <JANA/jerror.h>

#include <HDGEOMETRY/DMagneticFieldMap.h>

#include <vector>
#include <string>
using std::vector;
using std::string;

#include <JANA/JApplication.h>
#include <JANA/JCalibration.h>
using namespace jana;

class DMagneticFieldMapCalibDB:public DMagneticFieldMap{
   public:
      DMagneticFieldMapCalibDB(JApplication *japp, int32_t runnumber=1, string namepath = "Magnets/Solenoid/solenoid_1500_poisson_20090814_01");
      DMagneticFieldMapCalibDB(JCalibration *jcalib, string namepath = "Magnets/Solenoid/solenoid_1500_poisson_20090814_01");
      virtual ~DMagneticFieldMapCalibDB();
      
      int ReadMap(string namepath, int32_t runnumber=1, string context="");
      void GetField(const DVector3 &pos,DVector3 &Bout) const;
      void GetField(double x, double y, double z, double &Bx, double &By, double &Bz, int method=0) const;
      double GetBz(double x,double y, double z) const;
      void GetFieldGradient(double x, double y, double z,
                  double &dBxdx, double &dBxdy,
                  double &dBxdz,
                  double &dBydx, double &dBydy,
                  double &dBydz,    
                  double &dBzdx, double &dBzdy,
                  double &dBzdz) const;
      void GetFieldBicubic(double x,double y,double z,
                 double &Bx,double &By,double &Bz) const;
      void GetFieldAndGradient(double x,double y,double z,
                double &Bx,double &By,
                double &Bz,
                double &dBxdx, double &dBxdy,
                double &dBxdz,
                double &dBydx, double &dBydy,
                double &dBydz,
                double &dBzdx, double &dBzdy,
                double &dBzdz) const;

      typedef struct{
         float x,y,z,Bx,By,Bz;
         double dBxdx, dBxdy, dBxdz;
         double dBydx, dBydy, dBydz;
         double dBzdx, dBzdy, dBzdz;
        double dBxdxdy,dBxdxdz,dBxdydz;
        double dBydxdy,dBydxdz,dBydydz;
        double dBzdxdy,dBzdxdz,dBzdydz;
      }DBfieldPoint_t;

   protected:
      
      JCalibration *jcalib;
      vector< vector< vector<DBfieldPoint_t> > > Btable;

      float xmin, xmax, ymin, ymax, zmin, zmax;
      int Nx, Ny, Nz;
      double dx, dy,dz;
};

#endif // _DMagneticFieldMapCalibDB_