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