0b3636b366c2835f6e0d75c7ff03636b159e1ebe
[u/mrichter/AliRoot.git] / STEER / AliMagFDM.h
1 #ifndef ALIMAGFDM_H
2 #define ALIMAGFDM_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 "AliMagFC.h"
9 //
10 class AliMagFDM : public AliMagFC
11 {
12 //Alice Magnetic Field:Magnetic field map from IP to muon filter for Muon arm
13
14 public:
15   AliMagFDM();
16   AliMagFDM(const char *name, const char *title, const Int_t integ,
17             const Float_t factor, const Float_t fmax);
18   virtual ~AliMagFDM(){} 
19   virtual void Field(Float_t *x, Float_t *b);
20   virtual void ReadField(); 
21   virtual void SetSolenoidField(Float_t field = 2.) {fSolenoid = field;}
22   virtual Float_t SolenoidField() const {return fSolenoid;}
23   Int_t FZ(Double_t u, Float_t *Ar, Float_t du, Int_t ki, Int_t nu) const;
24   Double_t Ba(Int_t kai, Double_t za1, Double_t za2, Double_t al1, Double_t al2, Double_t al3, Int_t ka, Int_t ma);
25   Double_t Bb(Double_t z1, Double_t z2, Double_t y1, Double_t y2, Double_t x1, Double_t x2, Int_t kvr, Int_t k, Int_t l, Int_t m); 
26
27
28 protected:
29
30 //
31   Float_t fSolenoid; // Solenoid Field Strength
32   Int_t   fInd;      // Character number of validity Map region
33
34   Float_t fZmin;  // Start of the cartesian  part  of MAP in z
35   Float_t fZmax;  // End of Map in z   
36   Float_t fYmax;  // Start of the cartesian  part  of MAP in y
37   Float_t fYmin;  // End  of the cartesian  part  of MAP in y
38   Float_t fZpmx;  // End of the polar  part  of MAP in z
39   Float_t fZpmn;  // Start of the polar  part  of MAP in z
40   Float_t fRmax;  // Maximal radius of the polar  part  of MAP 
41   Float_t fRmin;  // Minimal radius of the polar  part  of MAP  
42               
43
44   Float_t    fXdel;  //  step in x - cartesian  part  of MAP
45   Float_t    fYdel;  //  step in y - cartesian  part  of MAP
46   Float_t    fZdel;  //  step in z - cartesian  part  of MAP
47   
48   Float_t    fRdel;  //  step in r - polar  part  of MAP
49   Float_t    fPhid;  //  step in Phi - polar  part  of MAP
50   Float_t    fZpdl;  //  step in z - polar  part  of MAP 
51   
52   Float_t    fCx1;   // Field constant
53   Float_t    fCx2;   // Field constant
54   Float_t    fAx1;   // Field constant
55   Float_t    fAx2;   // Field constant
56    
57   Float_t fZc[81];  // z coordinates in cartesian  part
58   Float_t fY[81];   // y coordinates in cartesian  part 
59   Float_t fBcx[81][81][44]; // Bx array for cartesian  part
60   Float_t fBcy[81][81][44]; // By array for cartesian  part
61   Float_t fBcz[81][81][44]; // Bz array for cartesian  part
62
63   Float_t  fZp[51];  // z coordinates in polar  part
64   Float_t  fR[10];   // r coordinates in polar  part  
65   Float_t  fPhi[33]; // Phi coordinates in polar  part
66
67   Float_t  fBpx[51][10][33]; // Bx array for polar  part
68   Float_t  fBpy[51][10][33]; // By array for polar  part
69   Float_t  fBpz[51][10][33]; // Bx array for polar  part 
70   Float_t  fB[2][2][32];     // Limits of field
71   
72   Int_t      fXl;    // Number steps in x for cartesian  part
73   Int_t      fYl;    // Number steps in y  for cartesian  par
74   Int_t      fZl;    // Number steps in z  for cartesian  part
75     
76   Int_t      fRn;    // Number steps in r for polar  part
77   Int_t      fPhin;  // Number steps in Phi for polar  part
78   Int_t      fZpl;   // Number steps in z for polar  part 
79   
80
81   ClassDef(AliMagFDM,1) //Class Magnetic field map from IP till muon filter
82 };
83
84 #endif