double pnt[3]={0.,0.,0.}; //Position where to get field strength
  double Bf[3]; //result goes here
  // retrieve the field from the framework
  FairRunAna::Instance()->GetField()->GetFieldValue ( pnt, Bf );  //[kGs]
  return Bf[2]; //OK, Bf[2] is your Bz