Additional background suppression in the V0 finder (M.Ivanov)

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

diff --git a/STEER/AliMagFC.cxx b/STEER/AliMagFC.cxx
index 4f665ede5e0321b60409f68b6216979357445adf..0233b0296865c9c8f21ef2f3f1aa6fb0adcaca31 100644 (file)
--- a/STEER/AliMagFC.cxx
+++ b/STEER/AliMagFC.cxx
@@ -23,13 +23,14 @@
 
 #include <stdlib.h>
 
+#include "AliLog.h"
 #include "AliMagFC.h"
 
 ClassImp(AliMagFC)
 
 //________________________________________
-AliMagFC::AliMagFC(const char *name, const char *title, const Int_t integ, 
-                  const Float_t factor, const Float_t fmax)
+AliMagFC::AliMagFC(const char *name, const char *title, Int_t integ, 
+                  Float_t factor, Float_t fmax)
   : AliMagF(name,title,integ,factor,fmax)
 {
   // 
@@ -40,7 +41,7 @@ AliMagFC::AliMagFC(const char *name, const char *title, const Int_t integ,
 }
 
 //________________________________________
-void AliMagFC::Field(Float_t *x, Float_t *b)
+void AliMagFC::Field(Float_t *x, Float_t *b) const
 {
   //
   // Method to return the field in a point
@@ -64,13 +65,12 @@ void AliMagFC::Field(Float_t *x, Float_t *b)
        b[2]*=fFactor;
     }
   } else {
-      printf("Invalid field map for constant field %d\n",fMap);
-      exit(1);
+      AliFatal(Form("Invalid field map for constant field %d",fMap));
   }
 }
 
 
-void AliMagFC::ZDCField(Float_t *x, Float_t *b)
+void AliMagFC::ZDCField(Float_t *x, Float_t *b) const
 {
 //This is the ZDC part
     Float_t rad2=x[0]*x[0]+x[1]*x[1];
@@ -83,25 +83,25 @@ void AliMagFC::ZDCField(Float_t *x, Float_t *b)
     else if(x[2] < kZ1BEG && x[2] > kZ1END){  
        if(rad2<kZ1RA2){
            b[0] =  kG1*x[1];
-           b[1] = -kG1*x[0];
+           b[1] =  kG1*x[0];
        }
     }
     else if(x[2] < kZ2BEG && x[2] > kZ2END){  
        if(rad2<kZ2RA2){
            b[0] = -kG1*x[1];
-           b[1] =  kG1*x[0];
+           b[1] = -kG1*x[0];
        }
     }
     else if(x[2] < kZ3BEG && x[2] > kZ3END){  
        if(rad2<kZ3RA2){
            b[0] = -kG1*x[1];
-           b[1] =  kG1*x[0];
+           b[1] = -kG1*x[0];
        }
     }
     else if(x[2] < kZ4BEG && x[2] > kZ4END){  
        if(rad2<kZ4RA2){
            b[0] =  kG1*x[1];
-           b[1] = -kG1*x[0];
+           b[1] =  kG1*x[0];
        }
     }
     else if(x[2] < kD1BEG && x[2] > kD1END){