1 GeV cut on gamma and e+/- in HALL outer concrete, air, steel

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

diff --git a/STEER/AliTracker.cxx b/STEER/AliTracker.cxx
index 75d16017ad040f577d2c690b1156d1d6790aa833..e0262c5125fa01f54e600624bc7c7b7ee81cc40d 100644 (file)
--- a/STEER/AliTracker.cxx
+++ b/STEER/AliTracker.cxx
@@ -23,12 +23,13 @@
 #include <TClass.h>
 #include <TMath.h>
 
+#include "AliMagF.h"
 #include "AliTracker.h"
 #include "AliCluster.h"
 #include "AliKalmanTrack.h"
 
 Bool_t AliTracker::fgUniformField=kTRUE;
-Double_t AliTracker::fgBz=0.;
+Double_t AliTracker::fgBz=kAlmost0Field;
 const AliMagF *AliTracker::fgkFieldMap=0;
 
 ClassImp(AliTracker)
@@ -83,7 +84,7 @@ void AliTracker::SetFieldMap(const AliMagF* map, Bool_t uni) {
   //Double_t bz=-b[2];
  
   Double_t bz=-map->SolenoidField();
-  fgBz=TMath::Sign(1e-13,bz) + bz;
+  fgBz=TMath::Sign(kAlmost0Field,bz) + bz;
 
 }
 
@@ -147,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);
+}