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