This commit was generated by cvs2svn to compensate for changes in r4472,
[u/mrichter/AliRoot.git] / HBTAN / AliHBTPair.cxx
1 #include "AliHBTPair.h"
2 #include "AliHBTParticle.h"
3
4 ClassImp(AliHBTPair)
5
6 //value of rev 
7 /************************************************************************/
8 AliHBTPair::AliHBTPair(Bool_t rev)
9  {
10 //value of rev defines if it is Swaped
11 //if you pass kTRUE swpaped pair will NOT be created
12 //though you wont be able to get the swaped pair from this pair
13
14   if(!rev) fSwapedPair = new AliHBTPair(kTRUE); //if false create swaped pair
15   else fSwapedPair = 0x0; //if true set the pointer to NULL
16   
17  }
18 /************************************************************************/
19 Double_t AliHBTPair::GetInvMass()
20 {
21   if(fInvMassNotCalc)
22    {
23      CalculateInvMassSqr(); //method is inline so we not waste th time for jumping into method 
24      
25      if(fInvMassSqr<0)  fInvMass = TMath::Sqrt(-fInvMassSqr);
26      else fInvMass = TMath::Sqrt(fInvMassSqr); 
27      
28      fInvMassNotCalc = kFALSE;
29    }
30   return fInvMass;
31 }
32 /************************************************************************/
33 Double_t AliHBTPair::GetQSideCMSLC()
34 {
35   //return Q Side in Central Of Mass System in Longitudialy Comoving Frame
36  
37   if (fQSideCMSLCNotCalc)
38    {
39     fQSideCMSLC = (fPart1->Px()*fPart2->Py()-fPart2->Px()*fPart1->Py())/GetKt();
40     fQSideCMSLCNotCalc = kFALSE;
41    }
42   return fQSideCMSLC;
43 }
44 /************************************************************************/
45 Double_t AliHBTPair::GetQOutCMSLC()
46 {
47  if(fQOutCMSLCNotCalc)
48   {
49    CalculateSums();
50    CalculateDiffs();
51    Double_t k2 = fPxSum*fPxDiff+fPySum*fPyDiff;
52    fQOutCMSLC = 0.5*k2/GetKt();
53    fQOutCMSLCNotCalc = kFALSE;
54   }
55  return fQOutCMSLC;
56 }
57 /************************************************************************/
58 Double_t AliHBTPair::GetQLongCMSLC()
59 {
60  if (fQLongCMSLCNotCalc)
61   {
62     CalculateSums();
63     CalculateDiffs();
64     Double_t beta = fPzSum/fESum;
65     Double_t gamma = 1.0/TMath::Sqrt(1.0 - beta*beta);
66     fQLongCMSLC = gamma*( fPzDiff - beta*fEDiff );
67     fQLongCMSLCNotCalc = kFALSE;
68   }
69  return fQLongCMSLC; 
70 }
71 /************************************************************************/
72 Double_t AliHBTPair::GetKt()
73 {
74   if(fKtNotCalc)
75    { 
76      CalculateSums();
77      fKt =  0.5*TMath::Hypot(fPxSum,fPySum);
78      fKtNotCalc = kFALSE;
79    }
80   return fKt;
81 }
82 /************************************************************************/
83
84 Double_t AliHBTPair::GetKStar()
85 {
86   if (fKStarNotCalc)
87    { 
88     
89     CalculateSums();
90
91     Double_t Ptrans = fPxSum*fPxSum + fPySum*fPySum;
92     Double_t Mtrans = fESum*fESum - fPzSum*fPzSum;
93     Double_t Pinv =   TMath::Sqrt(Mtrans - Ptrans);
94
95     Double_t Q = ( fPart1->GetMass()*fPart1->GetMass() - fPart2->GetMass()*fPart2->GetMass())/Pinv;
96     
97     CalculateQInvL();
98     
99     Q = TMath::Sqrt( Q*Q - fQInvL);
100     fKStar = Q/2;
101     fKStarNotCalc = kFALSE;
102    }
103   return fKStar;
104 }
105 /************************************************************************/
106
107 Double_t AliHBTPair::GetQInv()
108 {
109   if(fQInvNotCalc)
110    {
111     CalculateQInvL();
112
113     if (fQInvL < 0.) fQInv =  TMath::Sqrt(-fQInvL);
114     else             fQInv = -TMath::Sqrt( fQInvL);
115
116     fQInvNotCalc = kFALSE;
117    }
118    
119   return fQInv;
120 }
121
122
123
124
125
126
127
128