Remove lgamma and replace by TMath::LnGamma which should work across platforms.

[u/mrichter/AliRoot.git] / PWG2 / EVCHAR / AliCentralityGlauberFit.cxx

diff --git a/PWG2/EVCHAR/AliCentralityGlauberFit.cxx b/PWG2/EVCHAR/AliCentralityGlauberFit.cxx
index ee25481fe9cbcd735f46f9d64021d3de64fa12e0..f0086d20816dcb26f01af99fd89bcc50cee7f21b 100644 (file)
--- a/PWG2/EVCHAR/AliCentralityGlauberFit.cxx
+++ b/PWG2/EVCHAR/AliCentralityGlauberFit.cxx
@@ -490,7 +490,8 @@ void AliCentralityGlauberFit::SaveHisto(TH1D *hist1, TH1D *hist2, TH1D *heffi, T
 Double_t AliCentralityGlauberFit::NBD(Int_t n, Double_t mu, Double_t k)
 {
   // Compute NBD.
-  Double_t ret = exp( lgamma(n+k) - lgamma(k) - lgamma(n+1) ) * TMath::Power(mu/(mu+k),n) * TMath::Power(1-mu/(mu+k),k);
+  Double_t ret = TMath::Exp( TMath::LnGamma(n+k) - TMath::LnGamma(k) - TMath::LnGamma(n+1) ) * 
+                 TMath::Power(mu/(mu+k),n) * TMath::Power(1-mu/(mu+k),k);
   return ret;
 }