]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliMagFCM.h
Modified AddTracks. Should be backward compatible
[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 #include "AliMagF.h"
9 class TVector;
10
11 class AliMagFCM : public AliMagF
12 {
13   //Alice Magnetic Field with constan mesh
14
15 public:
16   AliMagFCM(){fB=0;}
17   AliMagFCM(const char *name, const char *title, const Int_t integ,
18            const Float_t factor, const Float_t fmax);
19   AliMagFCM(const AliMagFCM &mag);
20   virtual ~AliMagFCM() {delete fB;}
21   virtual void Field(Float_t *x, Float_t *b);
22   virtual void ReadField();
23   virtual void SetSolenoidField(Float_t field = 2.) {fSolenoid = field;}
24   
25   void Copy(AliMagFCM &magf) const;
26   virtual AliMagFCM & operator=(const AliMagFCM &magf);
27
28   Float_t Bx(const Int_t ix, const Int_t iy, const Int_t iz) {
29     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix));
30   }
31   Float_t By(const Int_t ix, const Int_t iy, const Int_t iz) {
32     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix)+1);
33   }
34   Float_t Bz(const Int_t ix, const Int_t iy, const Int_t iz) {
35     return (*fB)(3*(iz*(fXn*fYn)+iy*fXn+ix)+2);
36   }
37
38 protected:
39
40   Float_t    fXbeg;     // Start of mesh in x
41   Float_t    fYbeg;     // Start of mesh in y
42   Float_t    fZbeg;     // Start of mesh in z
43   Float_t    fXdel;     // Mesh step in x
44   Float_t    fYdel;     // Mesh step in y
45   Float_t    fZdel;     // Mesh step in z
46   Float_t    fSolenoid; // Solenoid Field Strength
47   Double_t   fXdeli;    // Inverse of Mesh step in x
48   Double_t   fYdeli;    // Inverse of Mesh step in y
49   Double_t   fZdeli;    // Inverse of Mesh step in z
50   Int_t      fXn;       // Number of mesh points in x
51   Int_t      fYn;       // Number of mesh points in y
52   Int_t      fZn;       // Number of mesh points in z
53   TVector   *fB;        // Field map
54   
55   ClassDef(AliMagFCM,1)  //Class for all Alice MagField with Constant Mesh
56 };
57
58 #endif