// 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 "AliLog.h"
-#include "AliRun.h"
-#include "AliMagF.h"
+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),
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) AliFatalClass("Can't access the field map !");
- AliKalmanTrack::SetConvConst(1000/0.299792458/(map->SolenoidField()+1e-13));
+
+ 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;
+
}
//__________________________________________________________________________
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);
+}