Merging correctly the changes with the previous code clean-up

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

diff --git a/HBTAN/AliHBTLLWeightFctn.cxx b/HBTAN/AliHBTLLWeightFctn.cxx
index d2dc5a4aac6a73fc9ea719a177a04fe0b906346a..192dde6ad724ecc46d7fb57e562a0b5c0cd02fce 100644 (file)
--- a/HBTAN/AliHBTLLWeightFctn.cxx
+++ b/HBTAN/AliHBTLLWeightFctn.cxx
@@ -8,7 +8,6 @@
 //One needs both pairs 
 //(simulated and recontructed), thus function is of class AliHBTTwoPairFctn1D.
 //Author: Ludmila Malinina, JINR (malinina@sunhe.jinr.ru)
-
 #include "AliHBTLLWeightFctn.h"
 #include "AliHBTLLWeights.h"
 
@@ -32,7 +31,7 @@ void  AliHBTLLWeightQInvFctn::ProcessSameEventParticles(AliHBTPair* trackpair, A
     Double_t weightPID=1.;
     Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
     Double_t weight=weightHBT*weightPID;
-    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQInv(),weight);
+    if(TMath::Abs(weight)<=10.)fNumerator->Fill(trackpair->GetQInv(),weight);
   }
 } 
 /****************************************************************/
@@ -55,3 +54,156 @@ TH1* AliHBTLLWeightQInvFctn::GetResult()
  return GetRatio(Scale());                                                  
 }                    
                                                               
+/**************************************************************************************/
+/**************************************************************************************/
+/**************************************************************************************/
+/**************************************************************************************/
+
+ClassImp(AliHBTLLWeightQOutFctn)
+/****************************************************************/
+void AliHBTLLWeightQOutFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)     
+  {
+    Double_t weightPID=1.;
+    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
+    Double_t weight=weightHBT*weightPID;
+    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQOutCMSLC(),weight);
+  }
+} 
+/****************************************************************/
+
+void AliHBTLLWeightQOutFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)  
+  {
+     fDenominator->Fill(trackpair->GetQOutCMSLC());
+  }
+}
+/**************************************************************/
+TH1* AliHBTLLWeightQOutFctn::GetResult() 
+                                                                               
+{ 
+//returns ratio of numerator and denominator                                    
+ return GetRatio(Scale());                                                  
+}                    
+                                                              
+/*************************************************************************************/ 
+/*************************************************************************************/ 
+/*************************************************************************************/ 
+
+ClassImp(AliHBTLLWeightQLongFctn)
+/****************************************************************/
+void  AliHBTLLWeightQLongFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+//Processes Particles and tracks Same different even
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)     
+  {
+    Double_t weightPID=1.;
+    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
+    Double_t weight=weightHBT*weightPID;
+    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQLongCMSLC(),weight);
+  }
+} 
+/****************************************************************/
+
+void  AliHBTLLWeightQLongFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)  
+  {
+     fDenominator->Fill(trackpair->GetQLongCMSLC());
+  }
+}
+/**************************************************************/
+TH1* AliHBTLLWeightQLongFctn::GetResult() 
+                                                                               
+{ 
+//returns ratio of numerator and denominator                                    
+ return GetRatio(Scale());                                                  
+}                    
+                                                              
+/*************************************************************************************/ 
+/*************************************************************************************/ 
+/*************************************************************************************/ 
+
+ClassImp(AliHBTLLWeightQSideFctn)
+/****************************************************************/
+void  AliHBTLLWeightQSideFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+//Processes Particles and tracks Same different even
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)     
+  {
+    Double_t weightPID=1.;
+    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
+    Double_t weight=weightHBT*weightPID;
+    if(TMath::Abs(weight)<=10.) fNumerator->Fill(trackpair->GetQSideCMSLC(),weight);
+  }
+} 
+/****************************************************************/
+
+void  AliHBTLLWeightQSideFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)  
+  {
+     fDenominator->Fill(trackpair->GetQSideCMSLC());
+  }
+}
+/**************************************************************/
+TH1* AliHBTLLWeightQSideFctn::GetResult() 
+                                                                               
+{ 
+//returns ratio of numerator and denominator                                    
+ return GetRatio(Scale());                                                  
+}                    
+                                                              
+/*************************************************************************************/ 
+/*************************************************************************************/ 
+
+ClassImp(AliHBTLLWeightTwoKStarFctn)
+/****************************************************************/
+void  AliHBTLLWeightTwoKStarFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+//Processes Particles and tracks Same different even
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)     
+  {
+    Double_t weightPID=1.;
+    Double_t weightHBT=AliHBTLLWeights::Instance()->GetWeight(partpair);
+    Double_t weight=weightHBT*weightPID;
+    if(TMath::Abs(weight)<=10.) fNumerator->Fill(2.0*(trackpair->GetKStar()),weight);
+  }
+} 
+/****************************************************************/
+
+void  AliHBTLLWeightTwoKStarFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
+{
+  trackpair = CheckPair(trackpair);
+  partpair  = CheckPair(partpair);
+  if ( trackpair && partpair)  
+  {
+     fDenominator->Fill(2.0*(trackpair->GetKStar()));
+  }
+}
+/**************************************************************/
+TH1* AliHBTLLWeightTwoKStarFctn::GetResult() 
+                                                                               
+{ 
+//returns ratio of numerator and denominator                                    
+ return GetRatio(Scale());                                                  
+}                    
+                                                              
+/*************************************************************************************/ 
+