Moved old AliMagFCheb to AliMagF, small fixes/optimizations in field classes
[u/mrichter/AliRoot.git] / STEER / AliMagFCM.h
1 #ifndef ALIMAGFCM_H
2 #define ALIMAGFCM_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id$ */
7
8 //-----------------------------------------------------------------------
9 //  Class for Alice magnetic field with constant mesh
10 //  Used in the configuration macros (macros/Config.C, etc.)
11 //  Author:
12 //-----------------------------------------------------------------------
13
14 #include "AliMagFC.h"
15 #include <TVector.h>
16
17 class AliMagFCM : public AliMagFC
18 {
19   //Alice Magnetic Field with constant mesh
20
21 public:
22   AliMagFCM();
23   AliMagFCM(const char *name, const char *title, Int_t integ,
24            Float_t factor, Float_t fmax);
25   AliMagFCM(const AliMagFCM &mag);
26   virtual ~AliMagFCM() {delete fB;}
27   virtual void Field(const float *x, float *b) const;
28   virtual void Field(const double *x, double *b) const;
29   virtual void ReadField();
30   virtual void    SetSolenoidField(Float_t field = 2.) {fSolenoid = field;}
31   virtual Float_t SolenoidField() const {
32     return -Factor()*fSolenoid;
33   }
34   
35   void Copy(TObject &magf) const;
36   virtual AliMagFCM & operator=(const AliMagFCM &magf)
37     {magf.Copy(*this); return *this;}
38
39   Float_t Bx(Int_t ix, Int_t iy, Int_t iz) const {
40     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix));
41   }
42   Float_t By(Int_t ix, Int_t iy, Int_t iz) const {
43     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix)+1);
44   }
45   Float_t Bz(Int_t ix, Int_t iy, Int_t iz) const {
46     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix)+2);
47   }
48
49 protected:
50
51   Float_t    fXbeg;     // Start of mesh in x
52   Float_t    fYbeg;     // Start of mesh in y
53   Float_t    fZbeg;     // Start of mesh in z
54   Float_t    fXdel;     // Mesh step in x
55   Float_t    fYdel;     // Mesh step in y
56   Float_t    fZdel;     // Mesh step in z
57   Float_t    fSolenoid; // Solenoid Field Strength
58   Double_t   fXdeli;    // Inverse of Mesh step in x
59   Double_t   fYdeli;    // Inverse of Mesh step in y
60   Double_t   fZdeli;    // Inverse of Mesh step in z
61   Int_t      fXn;       // Number of mesh points in x
62   Int_t      fYn;       // Number of mesh points in y
63   Int_t      fZn;       // Number of mesh points in z
64   TVector   *fB;        // Field map
65   
66   ClassDef(AliMagFCM,1)  //Class for all Alice MagField with Constant Mesh
67 };
68
69 #endif