+#include "AliHBTQResolutionFctns.h"
+
//__________________________________________________________________
////////////////////////////////////////////////////////////////////
// //
// General Remark: //
-// CMSLC means //
+// LCMS means //
// Center of Mass System Longitudially Co-moving //
// //
// //
// //
////////////////////////////////////////////////////////////////////
-#include "AliHBTQResolutionFctns.h"
/******************************************************************/
Rename("KtResolVsQInv","K_{t} Resolution vs. Q_{Inv}");
}
/******************************************************************/
+
+void AliHBTKtResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+//returns values of the functiion
+ y = partpair->GetKt() - trackpair->GetKt();
+ x = partpair->GetQInv();
+}
+/******************************************************************/
/******************************************************************/
/******************************************************************/
ClassImp( AliHBTQInvResolVsQInvFctn )
Rename("QInvResolVsQInv","Q_{Inv} Resolution vs. Q_{Inv}");
}
/******************************************************************/
+void AliHBTQInvResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+//returns values of the functiion
+ y = partpair->GetQInv() - trackpair->GetQInv();
+ x = partpair->GetQInv();
+}
+/******************************************************************/
/******************************************************************/
/******************************************************************/
//ctor
Rename("QOutResolVsQInv","Q_{Out} Resolution vs. Q_{Inv}");
}
+/******************************************************************/
+
+void AliHBTQOutResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qoutsim-Qoutrec for y
+ //returns Qinv for x
+ Double_t tqout = trackpair->GetQOutLCMS();
+ y = partpair->GetQOutLCMS() - tqout;
+ if (tqout < 0.0) y = -y;
+ x = partpair->GetQInv();
+}
+
/******************************************************************/
/******************************************************************/
/******************************************************************/
//ctor
Rename("QSideResolVsQInv","Q_{Side} Resolution vs. Q_{Inv}");
}
+/******************************************************************/
+
+void AliHBTQSideResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qsidesim-Qsiderec for y
+ //returns Qinv for x
+ y = partpair->GetQSideLCMS() - trackpair->GetQSideLCMS();
+ if (trackpair->GetQSideLCMS() < 0.0) y = -y;
+ x = partpair->GetQInv();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QLongResolVsQInv","Q_{Long} Resolution vs. Q_{Inv}");
}
+/******************************************************************/
+
+void AliHBTQLongResolVsQInvFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qlongsim-Qlongrec for y
+ //returns Qinv for x
+ y = partpair->GetQLongLCMS() - trackpair->GetQLongLCMS();
+ if (trackpair->GetQLongLCMS() < 0.0) y = -y;
+ x = partpair->GetQInv();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QInvResolVsKt","Q_{Inv} Resolution vs. K_{t}");
}
+/******************************************************************/
+
+void AliHBTQInvResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns values of the function
+ y = partpair->GetQInv() - trackpair->GetQInv();
+ x = partpair->GetKt();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QOutResolVsKt","Q_{Out} Resolution vs. K_{t} ");
}
+/******************************************************************/
+void AliHBTQOutResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qoutsim-Qoutrec for y
+ //returns Kt for x
+ y = partpair->GetQOutLCMS() - trackpair->GetQOutLCMS();
+ if (trackpair->GetQOutLCMS() < 0.0) y = -y;
+ x = partpair->GetKt();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QSideResolVsKt","Q_{Side} Resolution vs. K_{t} ");
}
+/******************************************************************/
+
+void AliHBTQSideResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qsidesim-Qsiderec for y
+ //returns Kt for x
+ y = partpair->GetQSideLCMS() - trackpair->GetQSideLCMS();
+ if (trackpair->GetQSideLCMS() < 0.0) y = -y;
+ x = partpair->GetKt();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QLongResolVsKt","Q_{Long} Resolution vs. K_{t} ");
}
+/******************************************************************/
+void AliHBTQLongResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Qlongsim-Qlongrec for y
+ //returns Kt for x
+ y = partpair->GetQLongLCMS() - trackpair->GetQLongLCMS();
+ if (trackpair->GetQLongLCMS() < 0.0) y = -y;
+ x = partpair->GetKt();
+}
/******************************************************************/
/******************************************************************/
/******************************************************************/
//ctor
Rename("QOutResolVsQOut","Q_{Out} Resolution vs. Q_{Out} ");
}
+/******************************************************************/
+
+void AliHBTQOutResolVsQOutFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+//returns values of the function
+ x = partpair->GetQOutLCMS();
+ y = x - trackpair->GetQOutLCMS();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QSideResolVsQSide","Q_{Side} Resolution vs. Q_{Side} ");
}
+/******************************************************************/
+
+void AliHBTQSideResolVsQSideFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+//returns values of the function
+ x = partpair->GetQSideLCMS();
+ y = x - trackpair->GetQSideLCMS();
+}
/******************************************************************/
/******************************************************************/
//ctor
Rename("QLongResolVsQLong","Q_{Long} Resolution vs. Q_{Long} ");
}
+/******************************************************************/
-
+void AliHBTQLongResolVsQLongFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+//returns values of the function
+ x = partpair->GetQLongLCMS();
+ y = x - trackpair->GetQLongLCMS();
+}
/******************************************************************/
/******************************************************************/
Rename("PairThetaResolVsQInv","Pair Theta Angle Resolution vs. Q_{Inv} ");
}
/******************************************************************/
+
+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();
+ }
+/******************************************************************/
/******************************************************************/
/******************************************************************/
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;
+}
+/******************************************************************/
/******************************************************************/
/******************************************************************/
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();
+ }
+
+/******************************************************************/
/******************************************************************/
/******************************************************************/
Rename("PairThetaResolVsKt","Pair Theta Angle Resolution vs. K_{t} ");
}
/******************************************************************/
+
+void AliHBTPairThetaResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Pair Theta sim - Pair Theta rec for y
+ //returns Kt sim 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->GetKt();
+}
+/******************************************************************/
/******************************************************************/
/******************************************************************/
//ctor
Rename("PairPhiResolVsKt","Pair Phi Angle Resolution vs. K_{t} ");
}
+/******************************************************************/
+
+void AliHBTPairPhiResolVsKtFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Pair Phi sim - Pair Phi rec for y
+ //returns Kt 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->GetKt();
+}
+
/******************************************************************/
/******************************************************************/
/******************************************************************/
//ctor
Rename("PairPhiResolVsPairPhi","Pair Phi Angle Resolution vs. Pair Phi ");
}
+/******************************************************************/
+
+void AliHBTPairPhiResolVsPairPhiFctn::GetValues(AliHBTPair* trackpair, AliHBTPair* partpair, Double_t& x, Double_t& y) const
+{
+ //returns Pair Phi sim - Pair Phi rec for y
+ //returns Pair Phi 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 = trackPhi;
+}
/******************************************************************/
/******************************************************************/