Adding new file also for the new Makefile

[u/mrichter/AliRoot.git] / HBTAN / AliHBTPair.h

diff --git a/HBTAN/AliHBTPair.h b/HBTAN/AliHBTPair.h
index 1d87ec9b80cdfdca7ba48f1be09f1148791b4e52..b416788e8561c6549a2e2e3c6c0e8b2a996557de 100644 (file)
--- a/HBTAN/AliHBTPair.h
+++ b/HBTAN/AliHBTPair.h
@@ -24,6 +24,8 @@ class AliHBTPair: public TObject
    
    Double_t GetInvMass(); //returns invariant mass of the pair
    
+
+   Double_t GetQInv(); //returns Q invariant
    Double_t GetQSideCMSLC(); //returns Q Side CMS longitudionally co-moving
    Double_t GetQOutCMSLC(); //returns Q out CMS longitudionally co-moving
    Double_t GetQLongCMSLC(); //returns Q Long CMS longitudionally co-moving
@@ -32,11 +34,10 @@ class AliHBTPair: public TObject
    Double_t GetKt();  //returns K transverse
    Double_t GetKStar();
    
-   Double_t GetQInv(); //returns Q invariant
-   Double_t GetQSide(); //returns Q side
-   Double_t GetQLong(); //returns Q long 
-   Double_t GetQOut(); //returns Q out
-   
+   Double_t GetDeltaP(); //return difference of momenta
+   Double_t GetDeltaPx();
+   Double_t GetDeltaPy();
+   Double_t GetDeltaPz();
    
  protected:
    AliHBTParticle* fPart1;  //pointer to first particle
@@ -68,7 +69,7 @@ class AliHBTPair: public TObject
    Double_t fKt; //K == sum vector of particle's momenta. Kt transverse component
    Bool_t   fKtNotCalc;//flag indicating if fKt is already calculated for this pair
    
-   Double_t fKStar;
+   Double_t fKStar; //
    Bool_t   fKStarNotCalc;
    
    Double_t fPInv;  //invariant momentum
@@ -138,11 +139,6 @@ void AliHBTPair::SetParticles(AliHBTParticle* p1,AliHBTParticle* p2)
  //and do nothing until will be asked for
 }
 /****************************************************************/
-inline
-Double_t AliHBTPair::GetQSide()
-{
- return 0.0;
-}
 /****************************************************************/
 inline 
 void AliHBTPair::CalculateInvMassSqr()
@@ -162,6 +158,7 @@ void AliHBTPair::CalculateInvMassSqr()
 inline 
 void AliHBTPair::CalculateQInvL()
  {
+ //Calculates square root of Qinv
   if (fQInvLNotCalc)
   {
    CalculateDiffs();
@@ -196,4 +193,35 @@ void AliHBTPair::CalculateDiffs()
     }
  }
 
+/****************************************************************/
+inline 
+Double_t AliHBTPair::GetDeltaP() //return difference of momenta
+{
+ CalculateDiffs();
+ return TMath::Sqrt(fPxDiff*fPxDiff + fPyDiff*fPyDiff + fPzDiff*fPzDiff);
+}
+/****************************************************************/
+inline 
+Double_t AliHBTPair::GetDeltaPx()
+ {
+   CalculateDiffs();
+   return fPxDiff;
+ }
+/****************************************************************/
+inline 
+Double_t AliHBTPair::GetDeltaPy()
+ {
+   CalculateDiffs();
+   return fPyDiff;
+ }
+
+/****************************************************************/
+inline 
+Double_t AliHBTPair::GetDeltaPz()
+ {
+   CalculateDiffs();
+   return fPzDiff;
+ }
+
+
 #endif