// $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, string namepath = "Magnets/Solenoid/solenoid_1200_poisson_20091123_01");
		DMagneticFieldMapCalibDB(JCalibration *jcalib, string namepath = "Magnets/Solenoid/solenoid_1200_poisson_20091123_01");
		virtual ~DMagneticFieldMapCalibDB();
		
		int ReadMap(string namepath, int runnumber=1, string context="");
		
		void GetField(double x, double y, double z, double &Bx, double &By, double &Bz, int method=0) 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_