]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliTracker.cxx
Support for automatic pedestal estimation per pad (Marian and Matevz).
[u/mrichter/AliRoot.git] / STEER / AliTracker.cxx
index e8da02f0adea5a4d79f05555ab9e0f860c7241fc..68f67faaa10647d0ca39999907a6d1c624f3c76b 100644 (file)
 //  that is the base for AliTPCtracker, AliITStrackerV2 and AliTRDtracker    
 //        Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
 //-------------------------------------------------------------------------
-
+#include <TClass.h>
 #include <TMath.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=0.;
 const AliMagF *AliTracker::fgkFieldMap=0;
 
 ClassImp(AliTracker)
 
 AliTracker::AliTracker():
-  fEventN(0),
-  fStoreBarrel(1),
   fX(0),
   fY(0),
   fZ(0),
@@ -48,18 +44,29 @@ 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) {
+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(1e-13,bz) + bz;
+
 }
 
 //__________________________________________________________________________