STEER/AliMagF.h

   1 #ifndef ALIMAGF_H
   2 #define ALIMAGF_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 //
   7 // Interface between the TVirtualMagField and AliMagWrapCheb: wrapper to
   8 // the set of magnetic field data + Tosca parameterization by
   9 // Chebyshev polynomials
  10 //
  11 // Author: ruben.shahoyan@cern.ch
  12 //
  13
  14 //#include <TGeoGlobalMagField.h>
  15 #include <TVirtualMagField.h>
  16 class AliMagWrapCheb;
  17
  18 class AliMagF : public TVirtualMagField
  19 {
  20  public:
  21   enum BMap_t      {k2kG, k5kG, k5kGUniform};
  22   enum BeamType_t  {kBeamTypeAA, kBeamTypepp, kNoBeamField};
  23   enum PolarityConvention_t {kConvLHC,kConvDCS2008,kConvMap2005};
  24   enum             {kOverrideGRP=BIT(14)}; // don't recreate from GRP if set
  25   //
  26   AliMagF();
  27   AliMagF(const char *name, const char* title, Int_t integ,
  28           Double_t factorSol=1., Double_t factorDip=1.,
  29           Double_t fmax=15, BMap_t maptype = k5kG,
  30           const char* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root",
  31           BeamType_t btype=kBeamTypepp, Double_t benergy=7000.);
  32   AliMagF(const AliMagF& src);
  33   AliMagF& operator=(const AliMagF& src);
  34   virtual ~AliMagF();
  35   //
  36   virtual void Field(const Double_t *x, Double_t *b);
  37   void       GetTPCInt(const Double_t *xyz, Double_t *b)        const;
  38   void       GetTPCIntCyl(const Double_t *rphiz, Double_t *b)   const;
  39   Double_t   GetBz(const Double_t *xyz)                         const;
  40   //
  41   AliMagWrapCheb* GetMeasuredMap()                              const {return fMeasuredMap;}
  42   //
  43   // former AliMagF methods or their aliases
  44   void         SetFactorSol(Float_t fc=1.);
  45   void         SetFactorDip(Float_t fc=1.);
  46   Double_t     GetFactorSol()                                   const;
  47   Double_t     GetFactorDip()                                   const;
  48   Double_t     Factor()                                         const {return GetFactorSol();}
  49   Bool_t       IsUniform()                                      const {return fMapType == k5kGUniform;}
  50   //
  51   void         MachineField(const Double_t  *x, Double_t *b)    const;
  52   BMap_t       GetMapType()                                     const {return fMapType;}
  53   BeamType_t   GetBeamType()                                    const {return fBeamType;}
  54   Double_t     GetBeamEnergy()                                  const {return fBeamEnergy;}
  55   Double_t     Max()                                            const {return fMax;}
  56   Int_t        Integ()                                          const {return fInteg;}
  57   Int_t        PrecInteg()                                      const {return fPrecInteg;}
  58   Double_t     SolenoidField()                                  const {return fFactorSol*fSolenoid;}
  59   //
  60   Char_t*      GetDataFileName()                                const {return (Char_t*)fParNames.GetName();}
  61   Char_t*      GetParamName()                                   const {return (Char_t*)fParNames.GetTitle();}
  62   void         SetDataFileName(const Char_t* nm)                      {fParNames.SetName(nm);}
  63   void         SetParamName(const Char_t* nm)                         {fParNames.SetTitle(nm);}
  64   //
  65   Bool_t       LoadParameterization();
  66   static Int_t GetPolarityConvention()                                {return Int_t(fgkPolarityConvention);}
  67   //
  68  protected:
  69   // not supposed to be changed during the run, set only at the initialization via constructor
  70   void         InitMachineField(BeamType_t btype, Double_t benergy);
  71   void         SetBeamType(BeamType_t type)                           {fBeamType = type;}
  72   void         SetBeamEnergy(Float_t energy)                          {fBeamEnergy = energy;}
  73   //
  74  protected:
  75   AliMagWrapCheb*  fMeasuredMap;     //! Measured part of the field map
  76   BMap_t           fMapType;         // field map type
  77   Double_t         fSolenoid;        // Solenoid field setting
  78   BeamType_t       fBeamType;        // Beam type: A-A (fBeamType=0) or p-p (fBeamType=1)
  79   Double_t         fBeamEnergy;      // Beam energy in GeV
  80   //
  81   Int_t            fInteg;           // Default integration method as indicated in Geant
  82   Int_t            fPrecInteg;       // Alternative integration method, e.g. for higher precision
  83   Double_t         fFactorSol;       // Multiplicative factor for solenoid
  84   Double_t         fFactorDip;       // Multiplicative factor for dipole
  85   Double_t         fMax;             // Max Field as indicated in Geant
  86   Bool_t           fDipoleOFF;       // Dipole ON/OFF flag
  87   //
  88   Double_t         fQuadGradient;    // Gradient field for inner triplet quadrupoles
  89   Double_t         fDipoleField;     // Field value for D1 and D2 dipoles
  90   Double_t         fCCorrField;      // Side C 2nd compensator field
  91   Double_t         fACorr1Field;     // Side A 1st compensator field
  92   Double_t         fACorr2Field;     // Side A 2nd compensator field
  93   //
  94   TNamed           fParNames;        // file and parameterization loadad
  95   //
  96   static const Double_t  fgkSol2DipZ;    // conventional Z of transition from L3 to Dipole field
  97   static const UShort_t  fgkPolarityConvention; // convention for the mapping of the curr.sign on main component sign
  98   //
  99   ClassDef(AliMagF, 2)           // Class for all Alice MagField wrapper for measured data + Tosca parameterization
 100 };
 101
 102
 103 #endif