]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
Updates.(D. Perrino)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 13 Dec 2009 20:32:02 +0000 (20:32 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 13 Dec 2009 20:32:02 +0000 (20:32 +0000)
STEER/AliAODDiJet.cxx
STEER/AliAODDiJet.h

index ba159f2e3f6584a65871332f976e2c73bc4768c2..3dcc516bd9cef6ef36a3c42ec7f2175257e7900b 100644 (file)
@@ -74,3 +74,30 @@ void AliAODDiJet::SetJetRefs(AliAODJet* jet1, AliAODJet* jet2)
     fJet2 = jet2;
 }
 
+
+//______________________________________________________________________________
+Float_t AliAODDiJet::DeltaPhi()
+{
+  // phi distance between the two jets
+  // the result is in the interval [0,pi]
+
+    Float_t phi1 = Jet(0)->Phi();
+    Float_t phi2 = Jet(1)->Phi();
+    Float_t dphi = TMath::Abs(phi1 - phi2);
+    dphi = (dphi > TMath::Pi()) ? (2*TMath::Pi()-dphi) : dphi;
+    return dphi;
+}
+
+
+//______________________________________________________________________________
+Float_t AliAODDiJet::PhiImbalance()
+{
+  // phi imbalance wrt back-to-back
+  // the result is in the interval [-pi,pi]
+
+    Float_t phi1 = Jet(0)->Phi();
+    Float_t phi2 = Jet(1)->Phi();
+    Float_t dphi = phi1 - phi2;
+    dphi = dphi - TMath::Sign((Float_t)TMath::Pi(),dphi);
+    return dphi;
+}
index 4cba120f98dff9116754a5eb030d121aaafcccf1..8bf257625d70043a6cd9354d2ece5c01acfbda79 100644 (file)
@@ -23,16 +23,21 @@ class AliAODDiJet : public AliAODJet {
     AliAODDiJet();
     AliAODDiJet(Double_t px, Double_t py, Double_t pz, Double_t e);
     AliAODDiJet(TLorentzVector & p);
+    virtual ~AliAODDiJet();
+
     void SetJetRefs(AliAODJet* jet1, AliAODJet* jet2);
     AliAODJet* Jet(Int_t i) {return ((AliAODJet*) (fJetR->At(i)));}
-    virtual ~AliAODDiJet();
+    Float_t    DeltaPhi();
+    Float_t    PhiImbalance();
+
  private:
-    AliAODDiJet(const AliAODDiJet& jet); 
+    AliAODDiJet(const AliAODDiJet& jet);
     AliAODDiJet& operator=(const AliAODDiJet& jet);
+
  private:
     TRefArray*  fJetR;
     TRef        fJet1;
-    TRef        fJet2;    
+    TRef        fJet2;
     ClassDef(AliAODDiJet, 1);
 };
 #endif