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;
+}
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