Parameters of miscalibration in the head of the file (Marian)

[u/mrichter/AliRoot.git] / STEER / AliTracker.cxx

diff --git a/STEER/AliTracker.cxx b/STEER/AliTracker.cxx
index e8da02f0adea5a4d79f05555ab9e0f860c7241fc..e0262c5125fa01f54e600624bc7c7b7ee81cc40d 100644 (file)
--- a/STEER/AliTracker.cxx
+++ b/STEER/AliTracker.cxx
@@ -20,24 +20,22 @@
 //  that is the base for AliTPCtracker, AliITStrackerV2 and AliTRDtracker    
 //        Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
 //-------------------------------------------------------------------------
-
+#include <TClass.h>
 #include <TMath.h>
 
+#include "AliMagF.h"
 #include "AliTracker.h"
 #include "AliCluster.h"
 #include "AliKalmanTrack.h"
-#include "AliRun.h"
-#include "AliMagF.h"
-
-extern AliRun* gAlice;
 
+Bool_t AliTracker::fgUniformField=kTRUE;
+Double_t AliTracker::fgBz=kAlmost0Field;
 const AliMagF *AliTracker::fgkFieldMap=0;
 
 ClassImp(AliTracker)
 
 AliTracker::AliTracker():
-  fEventN(0),
-  fStoreBarrel(1),
+  TObject(),
   fX(0),
   fY(0),
   fZ(0),
@@ -48,18 +46,46 @@ AliTracker::AliTracker():
   //--------------------------------------------------------------------
   // The default constructor.
   //--------------------------------------------------------------------
- AliMagF *field=gAlice->Field();
- if (field==0) Fatal("AliTracker()","Can't access the field map !");
- SetFieldMap(field);
+  if (!fgkFieldMap) AliWarning("Field map is not set. Call AliTracker::SetFieldMap before creating a tracker!");
 }
 
-void AliTracker::SetFieldMap(const AliMagF* map) {
+//__________________________________________________________________________
+AliTracker::AliTracker(const AliTracker &atr):
+  TObject(atr),
+  fX(atr.fX),
+  fY(atr.fY),
+  fZ(atr.fZ),
+  fSigmaX(atr.fSigmaX),
+  fSigmaY(atr.fSigmaY),
+  fSigmaZ(atr.fSigmaZ)
+{
+  //--------------------------------------------------------------------
+  // The default constructor.
+  //--------------------------------------------------------------------
+  if (!fgkFieldMap) AliWarning("Field map is not set. Call AliTracker::SetFieldMap before creating a tracker!");
+}
+
+//__________________________________________________________________________
+void AliTracker::SetFieldMap(const AliMagF* map, Bool_t uni) {
   //--------------------------------------------------------------------
   //This passes the field map to the reconstruction.
   //--------------------------------------------------------------------
-  if (map==0) ::Fatal("SetFieldMap","Can't access the field map !");
-  AliKalmanTrack::SetConvConst(1000/0.299792458/map->SolenoidField());
+  if (map==0) AliFatalClass("Can't access the field map !");
+
+  if (fgkFieldMap) {
+     AliWarningClass("The magnetic field map has been already set !");
+     return;
+  }
+
+  fgUniformField=uni;
   fgkFieldMap=map;
+
+  //Float_t r[3]={0.,0.,0.},b[3]; map->Field(r,b);
+  //Double_t bz=-b[2];
+ 
+  Double_t bz=-map->SolenoidField();
+  fgBz=TMath::Sign(kAlmost0Field,bz) + bz;
+
 }
 
 //__________________________________________________________________________
@@ -122,3 +148,12 @@ void AliTracker::UseClusters(const AliKalmanTrack *t, Int_t from) const {
      c->Use();   
   }
 }
+
+Double_t AliTracker::GetBz(Float_t *r) {
+  //------------------------------------------------------------------
+  // Returns Bz (kG) at the point "r" .
+  //------------------------------------------------------------------
+    Float_t b[3]; fgkFieldMap->Field(r,b);
+    Double_t bz=-Double_t(b[2]);
+    return  (TMath::Sign(kAlmost0Field,bz) + bz);
+}