3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 // Interface between the TVirtualMagField and AliMagWrapCheb: wrapper to
8 // the set of magnetic field data + Tosca parameterization by
9 // Chebyshev polynomials
11 // Author: ruben.shahoyan@cern.ch
14 //#include <TGeoGlobalMagField.h>
15 #include <TVirtualMagField.h>
18 class AliMagF : public TVirtualMagField
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
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);
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;
41 AliMagWrapCheb* GetMeasuredMap() const {return fMeasuredMap;}
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;}
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;}
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);}
65 Bool_t LoadParameterization();
66 static Int_t GetPolarityConvention() {return Int_t(fgkPolarityConvention);}
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;}
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
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
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
94 TNamed fParNames; // file and parameterization loadad
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
99 ClassDef(AliMagF, 2) // Class for all Alice MagField wrapper for measured data + Tosca parameterization