#ifndef HFIELD_H
#define HFIELD_H
#include "TNamed.h"
#include "TString.h"

//Read the field
//Input: 3 dimensional Position vector 
//Output:3 dimensional Field vector
//created 16/01/2002/ 
//********************************************
//  nfz      number of steps in z-direction 
//  nfr      number of steps in rho-direction 
//  nfp      number of steps in phi-direction 
//  zflmin   minimum value of z-coordinate
//  zflmax   maximum value of z-coordinate
//  zfldel   step value of z-coordinate
//  rflmin   minimum value of rho-coordinate
//  rflmax   maximum value of rho-coordinate
//  rfldel   step value of rho-coordinate
//  pflmin   minimum value of phi-coordinate
//  pflmax   maximum value of phi-coordinate
//  pfldel   step value of phi-coordinate


// coordinate system:
// rho = sqr(x2 + y2)
// phi = atan y/x
// z   = z (beam)

class HField : public TNamed {
protected:
    Int_t nfz,nfr,nfp;
    Double_t zflmin,zflmax,zfldel;
    Double_t rflmin,rflmax,rfldel;
    Double_t pflmin,pflmax,pfldel;
    Double_t tzfl[161][176][25];
    Double_t trfl[161][176][25];
    Double_t tpfl[161][176][25];
public:
    HField(const Char_t* name="",const Char_t* title="");
    ~HField();
    void  init(TString);
    void  calcField(Double_t*,Double_t*,Int_t);
    void  initVariables();
    ClassDef(HField,1) // Class for Field Map
};

#endif  /*!HFIELD_H*/