]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
nan check (needed if class is used in standalone root)
authorjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 27 Oct 2010 09:55:03 +0000 (09:55 +0000)
committerjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 27 Oct 2010 09:55:03 +0000 (09:55 +0000)
PWG0/AliUnfolding.cxx

index b3dd89f93cddacc2fc3731736468fbb2f457b4e4..b35c6e2d394ac5e28029b007a3c1d1a5300838b5 100644 (file)
@@ -443,7 +443,11 @@ Int_t AliUnfolding::UnfoldWithMinuit(TH2* correlation, TH1* efficiency, TH1* mea
     results[i] = minuit->GetParameter(i);
     Double_t value = results[i] * results[i];
     // error is : (relError) * (value) = (minuit->GetParError(i) / minuit->GetParameter(i)) * (minuit->GetParameter(i) * minuit->GetParameter(i))
-    Double_t error = minuit->GetParError(i) * results[i];
+    Double_t error = 0;
+    if (TMath::IsNaN(minuit->GetParError(i)))
+      Printf("WARNING: Parameter %d error is nan", i);
+    else 
+      error = minuit->GetParError(i) * results[i];
     
     if (efficiency)
     {  
@@ -465,7 +469,7 @@ Int_t AliUnfolding::UnfoldWithMinuit(TH2* correlation, TH1* efficiency, TH1* mea
   
   fgCallCount = 0;
   Chi2Function(dummy, 0, chi2, results, 0);
-  printf("AliUnfolding::UnfoldWithMinuit: Chi2 of final parameters is = %f\n", chi2);
+  Printf("AliUnfolding::UnfoldWithMinuit: Chi2 of final parameters is = %f", chi2);
   
   if (fgDebug)
     DrawGuess(results);