fix in the average separation calculation

diff --git a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx
index f1374d05e0b4813caf4a7ed26f82f97108617023..d340cbb4ed0aff4b243f3010c35b6c4ca383c343 100644 (file)
--- a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx
+++ b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx
@@ -910,7 +910,17 @@ void AliFemtoEventReaderAOD::CopyAODtoFemtoTrack(AliAODTrack *tAodTrack,
   tFemtoTrack->SetTPCSharedMap(tAodTrack->GetTPCSharedMap());
 
   float globalPositionsAtRadii[9][3];
+
+  float s = fEvent->GetMagneticField();
+  if(s>0)
+    fMagFieldSign = 1;
+  else if(s<0)
+    fMagFieldSign = -1;
+  else
+    fMagFieldSign = 0;
+
   float bfield = 5*fMagFieldSign;
+
   GetGlobalPositionAtGlobalRadiiThroughTPC(tAodTrack,bfield,globalPositionsAtRadii);
   double tpcEntrance[3]={globalPositionsAtRadii[0][0],globalPositionsAtRadii[0][1],globalPositionsAtRadii[0][2]};
   double **tpcPositions;