GetBZ (using in real tracking) was returning hardwired 5kGauss for +30kA etc. Changed...

diff --git a/STEER/AliMagF.cxx b/STEER/AliMagF.cxx
index 71852138f942e4ea152012d2fd9afe36de9cbbd7..c99e898fd1686e2e6b14bce686bbb625e1f2673f 100644 (file)
--- a/STEER/AliMagF.cxx
+++ b/STEER/AliMagF.cxx
@@ -92,26 +92,20 @@ AliMagF::AliMagF(const char *name, const char* title, Int_t integ,
   //
   const char* parname = 0;
   //  
-  if (fMapType == k2kG) {
-    fSolenoid = 2.;
-    parname = fDipoleOFF ? "Sol12_Dip0_Hole":"Sol12_Dip6_Hole";
-  } else if (fMapType == k5kG) {
-    fSolenoid = 5.;
-    parname = fDipoleOFF ? "Sol30_Dip0_Hole":"Sol30_Dip6_Hole";
-  } else if (fMapType == k5kGUniform) {
-    fSolenoid = 5.;
-    parname = "Sol30_Dip6_Uniform";
-  } else {
-    AliFatal(Form("Unknown field identifier %d is requested\n",fMapType)); 
-  }
+  if      (fMapType == k2kG) parname = fDipoleOFF ? "Sol12_Dip0_Hole":"Sol12_Dip6_Hole";
+  else if (fMapType == k5kG) parname = fDipoleOFF ? "Sol30_Dip0_Hole":"Sol30_Dip6_Hole";
+  else if (fMapType == k5kGUniform) parname = "Sol30_Dip6_Uniform";
+  else AliFatal(Form("Unknown field identifier %d is requested\n",fMapType));
   //
   SetDataFileName(path);
   SetParamName(parname);
   //
-  SetFactorSol(factorSol);
-  SetFactorDip(factorDip);
   LoadParameterization();
   InitMachineField(fBeamType,fBeamEnergy);
+  double xyz[3]={0.,0.,0.};
+  fSolenoid = GetBz(xyz);
+  SetFactorSol(factorSol);
+  SetFactorDip(factorDip);
 }
 
 //_______________________________________________________________________

diff --git a/STEER/AliMagF.h b/STEER/AliMagF.h
index 7b7df4348705fc00c2c749e4c852fce00f459dae..b7694b3fdf7e959546a487db53f20f65a726d21c 100644 (file)
--- a/STEER/AliMagF.h
+++ b/STEER/AliMagF.h
@@ -53,7 +53,7 @@ class AliMagF : public TVirtualMagField
   Double_t     Max()                                            const {return fMax;}
   Int_t        Integ()                                          const {return fInteg;}
   Int_t        PrecInteg()                                      const {return fPrecInteg;}  
-  Double_t     SolenoidField()                                  const {return -fFactorSol*fSolenoid;}
+  Double_t     SolenoidField()                                  const {return fFactorSol*fSolenoid;}
   //
   Char_t*      GetDataFileName()                                const {return (Char_t*)fParNames.GetName();}
   Char_t*      GetParamName()                                   const {return (Char_t*)fParNames.GetTitle();}