+void AliHBTPairThetaResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+ {
+ //returns Pair Theta sim - Pair Theta rec for y
+ //returns Qinv for x
+ Double_t partTheta = partpair->Particle1()->Theta() - partpair->Particle2()->Theta();
+ Double_t trackTheta = trackpair->Particle1()->Theta() - trackpair->Particle2()->Theta();
+ y = partTheta - trackTheta;
+ x = partpair->GetQInv();
+ }
+/******************************************************************/
+/******************************************************************/
+/******************************************************************/
+
+ClassImp( AliHBTPairThetaResolVsPairThetaFctn )
+
+AliHBTPairThetaResolVsPairThetaFctn::
+AliHBTPairThetaResolVsPairThetaFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
+ Int_t nYbins, Double_t maxYval, Double_t minYval):
+ AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
+{
+//ctor
+ Rename("PairThetaResolVsPairTheta","Pair Theta Angle Resolution vs. Pair Theta ");
+}
+/******************************************************************/
+void AliHBTPairThetaResolVsPairThetaFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Pair Theta sim - Pair Theta rec for y
+ //returns Pair Theta sim for x
+ Double_t partTheta = partpair->Particle1()->Theta() - partpair->Particle2()->Theta();
+ Double_t trackTheta = trackpair->Particle1()->Theta() - trackpair->Particle2()->Theta();
+ y = partTheta - trackTheta;
+ if (trackTheta < 0.0) y = -y;
+ x = trackTheta;
+}
+/******************************************************************/
+/******************************************************************/
+/******************************************************************/
+
+ClassImp( AliHBTPairPhiResolVsQInvFctn )
+
+AliHBTPairPhiResolVsQInvFctn::
+AliHBTPairPhiResolVsQInvFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
+ Int_t nYbins, Double_t maxYval, Double_t minYval):
+ AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
+{
+//ctor
+ Rename("PairPhiResolVsQInv","Pair Phi Angle Resolution vs. Q_{Inv} ");
+}
+/******************************************************************/
+
+void AliHBTPairPhiResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+ {
+ //returns Pair Phi sim - Pair Phi rec for y
+ //returns QInv sim for x
+ Double_t partPhi = partpair->Particle1()->Phi() - partpair->Particle2()->Phi();
+ Double_t trackPhi = trackpair->Particle1()->Phi() - trackpair->Particle2()->Phi();
+ y = partPhi - trackPhi;
+ x = partpair->GetQInv();
+ }
+
+/******************************************************************/
+/******************************************************************/
+/******************************************************************/
+
+
+ClassImp( AliHBTPairThetaResolVsKtFctn )
+
+AliHBTPairThetaResolVsKtFctn::
+AliHBTPairThetaResolVsKtFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
+ Int_t nYbins, Double_t maxYval, Double_t minYval):
+ AliHBTTwoPairFctn2D(nXbins,maxXval,minXval,nYbins,maxYval,minYval)
+{
+//ctor
+ Rename("PairThetaResolVsKt","Pair Theta Angle Resolution vs. K_{t} ");
+}
+/******************************************************************/