reverting to prev.

[u/mrichter/AliRoot.git] / HBTAN / AliHBTFunction.cxx

diff --git a/HBTAN/AliHBTFunction.cxx b/HBTAN/AliHBTFunction.cxx
index 8fb3f3ec2d22e133cd92c53f429318f6e92b2275..04be138f56b37deaf3b7d4ab3d1ab808c1740796 100644 (file)
--- a/HBTAN/AliHBTFunction.cxx
+++ b/HBTAN/AliHBTFunction.cxx
@@ -121,6 +121,7 @@ TH1* AliHBTFunction::GetRatio(Double_t normfactor)
     }
    TString str = fName + " ratio";
    TH1 *result = (TH1*)GetNumerator()->Clone(str.Data());
+   result->SetDirectory(0x0);
    
    result->SetTitle(str.Data());
    
@@ -234,6 +235,7 @@ AliHBTCorrelFunction& AliHBTCorrelFunction::operator=(const AliHBTCorrelFunction
   if (&in == this) return *this;
   delete fRatio;
   fRatio=(in.fRatio)?(TH1*)in.fRatio->Clone():0x0;
+  fRatio->SetDirectory(0x0);
   return *this;
 }
 
@@ -393,8 +395,8 @@ Double_t AliHBTFunction1D::Scale(TH1D* num,TH1D* den)
    {
     if ( num->GetBinContent(i) > 0.0 )
      {
-       densum = fDenominator->GetBinContent(i);
-       numsum = fNumerator->GetBinContent(i);
+       densum += fDenominator->GetBinContent(i);
+       numsum += fNumerator->GetBinContent(i);
      }
    }
   
@@ -587,8 +589,8 @@ Double_t AliHBTFunction2D::Scale()
      {
       if ( fNumerator->GetBinContent(i,j) > 0.0 )
        {
-           densum = fDenominator->GetBinContent(i,j);
-           numsum = fNumerator->GetBinContent(i,j);
+           densum += fDenominator->GetBinContent(i,j);
+           numsum += fNumerator->GetBinContent(i,j);
        }
      }
   
@@ -803,8 +805,8 @@ Double_t AliHBTFunction3D::Scale()
         if ( fNumerator->GetBinContent(i,j,k) > 0.0 )
          {
            
-           densum = fDenominator->GetBinContent(i,j,k);
-           numsum = fNumerator->GetBinContent(i,j,k);
+           densum += fDenominator->GetBinContent(i,j,k);
+           numsum += fNumerator->GetBinContent(i,j,k);
          }
        }