From 387dd2d0b5d3b5578d83536798aa148ab622f423 Mon Sep 17 00:00:00 2001
From: skowron <skowron@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 22 Nov 2004 18:11:11 +0000
Subject: [PATCH] Bug corrected - pairs must be filled from mixed events

---
 HBTAN/AliHBTCorrFitFctn.cxx | 49 ++++++++++++++++++++-----------------
 1 file changed, 27 insertions(+), 22 deletions(-)

diff --git a/HBTAN/AliHBTCorrFitFctn.cxx b/HBTAN/AliHBTCorrFitFctn.cxx
index ab5f9a64a39..5b77f060a8e 100644
--- a/HBTAN/AliHBTCorrFitFctn.cxx
+++ b/HBTAN/AliHBTCorrFitFctn.cxx
@@ -32,30 +32,11 @@ void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTP
    if(trackpair == 0x0) return; 
    
    Double_t q = trackpair->GetQInv();
-   Bool_t fill = kFALSE;
    
    Double_t weight = partpair->GetWeight();
    fNumerator->Fill(q,weight);
    
-   if ( (q < 0.15) && (fNPairsFitArea < 2.e+5))
-     {
-       fNPairsFitArea++;
-       fill = kTRUE;
-     }
 
-   if ( (q > 0.15) && (q < 0.3) && (fNPairsFitArea < 1.e+5))
-     {
-       fNPairsNormArea++;
-       fill = kTRUE;
-     }
-   
-   if (fill)
-    {  
-      const AliVAODParticle& p1 = *(trackpair->Particle1());
-      const AliVAODParticle& p2 = *(trackpair->Particle2());
-      fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.E(),
-                    p2.Px(),p2.Py(),p2.Pz(),p2.E());
-    }
 }
 /****************************************************************/
 
@@ -65,9 +46,33 @@ void  AliHBTCorrFitFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBT
   trackpair = CheckPair(trackpair);
 //  partpair  = CheckPair(partpair);
   if ( trackpair && partpair)
-  {
-     fDenominator->Fill(trackpair->GetQInv());
-  }
+   {
+     Double_t q = trackpair->GetQInv();
+
+     Bool_t fill = kFALSE;
+    
+     if ( (q < 0.15) && (fNPairsFitArea < 2.e+5))
+       {
+         fNPairsFitArea++;
+         fill = kTRUE;
+       }
+
+     if ( (q > 0.15) && (q < 0.3) && (fNPairsFitArea < 1.e+5))
+       {
+         fNPairsNormArea++;
+         fill = kTRUE;
+       }
+
+     if (fill)
+      {  
+        const AliVAODParticle& p1 = *(trackpair->Particle1());
+        const AliVAODParticle& p2 = *(trackpair->Particle2());
+        fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.E(),
+                      p2.Px(),p2.Py(),p2.Pz(),p2.E());
+      }
+
+       fDenominator->Fill(q);
+   }
 }
 /*****************************************************************/
 
-- 
2.39.3