Q transverse added
[u/mrichter/AliRoot.git] / ANALYSIS / AliAODPair.cxx
index 31de69be0aaa4e9b905983835351d0b6d5a2f408..acd30e75f4b52a738c85360ca8e0e280cae4146a 100644 (file)
@@ -27,6 +27,8 @@ AliAODPair::AliAODPair(Bool_t rev):
  fQOutLCMSNotCalc(kTRUE),
  fQLongLCMS(0.0),
  fQLongLCMSNotCalc(kTRUE),
+ fQtLCMS(0.0),
+ fQtLCMSNotCalc(kTRUE),
  fQInv(0.0),
  fQInvNotCalc(kTRUE),
  fInvMass(0.0),
@@ -80,6 +82,8 @@ AliAODPair::AliAODPair(AliVAODParticle* part1, AliVAODParticle* part2, Bool_t re
  fQOutLCMSNotCalc(kTRUE),
  fQLongLCMS(0.0),
  fQLongLCMSNotCalc(kTRUE),
+ fQtLCMS(0.0),
+ fQtLCMSNotCalc(kTRUE),
  fQInv(0.0),
  fQInvNotCalc(kTRUE),
  fInvMass(0.0),
@@ -133,6 +137,8 @@ AliAODPair::AliAODPair(const AliAODPair& in):
  fQOutLCMSNotCalc(kTRUE),
  fQLongLCMS(0.0),
  fQLongLCMSNotCalc(kTRUE),
+ fQtLCMS(0.0),
+ fQtLCMSNotCalc(kTRUE),
  fQInv(0.0),
  fQInvNotCalc(kTRUE),
  fInvMass(0.0),
@@ -298,6 +304,17 @@ Double_t AliAODPair::GetQLongLCMS()
 }
 /************************************************************************/
 
+Double_t AliAODPair::GetQtLCMS()
+{
+ //returns Q transverse CMS longitudionally co-moving
+ if (fQtLCMSNotCalc)
+  {
+    fQtLCMS = TMath::Hypot(GetQOutLCMS(),GetQSideLCMS());
+    fQtLCMSNotCalc = kFALSE;
+  }
+ return fQtLCMS; 
+}
+
 Double_t AliAODPair::GetKt()
 {
  //calculates the evarage momentum of the pair