Bug corrected - pairs must be filled from mixed events
authorskowron <skowron@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 22 Nov 2004 18:11:11 +0000 (18:11 +0000)
committerskowron <skowron@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 22 Nov 2004 18:11:11 +0000 (18:11 +0000)
HBTAN/AliHBTCorrFitFctn.cxx

index ab5f9a64a3993bfb65d2c0ea54772a9b0535a19b..5b77f060a8e09561543081150bcb20f05fefdd11 100644 (file)
@@ -32,30 +32,11 @@ void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTP
    if(trackpair == 0x0) return; 
    
    Double_t q = trackpair->GetQInv();
    if(trackpair == 0x0) return; 
    
    Double_t q = trackpair->GetQInv();
-   Bool_t fill = kFALSE;
    
    Double_t weight = partpair->GetWeight();
    fNumerator->Fill(q,weight);
    
    
    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)
   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);
+   }
 }
 /*****************************************************************/
 
 }
 /*****************************************************************/