Bug Correction (lacking Abs)

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

diff --git a/HBTAN/AliHBTWeightFctn.cxx b/HBTAN/AliHBTWeightFctn.cxx
index f0538a9367a97d7450d1674c96a33e1b2baff2ac..ea10f41a0b7e3d00f39fe9035084b3c4ab69890c 100644 (file)
--- a/HBTAN/AliHBTWeightFctn.cxx
+++ b/HBTAN/AliHBTWeightFctn.cxx
@@ -495,6 +495,17 @@ void AliHBTWeightQOutSQideQLongFctn::ProcessSameEventParticles(AliHBTPair* track
     Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
     Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
     Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
+    
+    if (out < 0.005)
+     if (side < 0.005)
+       if (lon < 0.005)
+        {
+          trackpair->Particle1()->Print();
+          trackpair->Particle2()->Print();
+          Info("","Delta Theta %f, Delta Phi %f",
+              trackpair->GetDeltaTheta(),trackpair->GetDeltaPhi());
+        }
+    
     fNumerator->Fill(out,side,lon,weight);//here we fill in q's corresponding to track pair 
                                           //weight calculated for the simulated one
   }
@@ -508,7 +519,11 @@ void AliHBTWeightQOutSQideQLongFctn::ProcessDiffEventParticles(AliHBTPair* track
   partpair  = CheckPair(partpair);
   if ( trackpair && partpair)  
    {
-     fDenominator->Fill(trackpair->GetQOutCMSLC(),trackpair->GetQSideCMSLC(),trackpair->GetQLongCMSLC());
+     Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
+     Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
+     Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
+   
+     fDenominator->Fill(out,side,lon);
    }
 }
 /*************************************************************/