1 #include "AliHBTTwoTrackEffFctn.h"
2 //____________________________________________________________________
3 //////////////////////////////////////////////////////////////////////
5 // class AliHBTTwoTrackEffFctn //
7 // classes for calculating two track efficiency of the tracking //
8 // binning is done using value of simulated pair montum difference //
9 // pair must be recontructed, that is why we need both pairs //
10 // (simulated and recontructed), thus functions are "two pair" //
11 // Piotr.Skowronski@cern.ch //
13 //////////////////////////////////////////////////////////////////////
16 ClassImp(AliHBTTwoTrackEffFctn)
17 /******************************************************************/
19 AliHBTTwoTrackEffFctn::AliHBTTwoTrackEffFctn()
23 /******************************************************************/
25 AliHBTTwoTrackEffFctn::AliHBTTwoTrackEffFctn(Int_t nbins, Double_t maxval, Double_t minval):
26 AliHBTOnePairFctn1D("TwoTrackEff","Two Track Efficiency",nbins,maxval,minval)
29 //nbins - numner of bins of the function
30 //maxval - max X of the fctn
31 //minval - min X of the fctn
32 GetNumerator()->GetXaxis()->SetTitle("dP [GeV]");
33 GetDenominator()->GetXaxis()->SetTitle("dP [GeV]");
35 /******************************************************************/
37 TH1* AliHBTTwoTrackEffFctn::GetResult()
39 //returns ratio of numerator and denominator
41 fRatio = GetRatio(Scale());
44 fRatio->GetXaxis()->SetTitle("dP [GeV]");
45 fRatio->GetYaxis()->SetTitle("C(dP)");
46 fRatio->SetTitle("Double Track Resolution: dP Correlation Fctn.");
50 /******************************************************************/
51 /******************************************************************/
52 /******************************************************************/
53 ClassImp(AliHBTTwoTrackEffFctnPxPyPz)
55 AliHBTTwoTrackEffFctnPxPyPz::AliHBTTwoTrackEffFctnPxPyPz(Int_t nXbins, Double_t maxXval, Double_t minXval,
56 Int_t nYbins, Double_t maxYval, Double_t minYval,
57 Int_t nZbins, Double_t maxZval, Double_t minZval):
58 AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
62 fWriteNumAndDen = kTRUE;
63 Rename("tteffpxpypz","P_{x} P_{y} P_{z} Two Track Efficiency Function");
66 fNumerator->GetXaxis()->SetTitle("\\Delta P_{x} [GeV]");
67 fNumerator->GetYaxis()->SetTitle("\\Delta P_{y} [GeV]");
68 fNumerator->GetZaxis()->SetTitle("\\Delta P_{z} [GeV]");
73 fDenominator->GetXaxis()->SetTitle("\\Delta P_{x} [GeV]");
74 fDenominator->GetYaxis()->SetTitle("\\Delta P_{y} [GeV]");
75 fDenominator->GetZaxis()->SetTitle("\\Delta P_{z} [GeV]");
79 /******************************************************************/
81 void AliHBTTwoTrackEffFctnPxPyPz::GetValues(AliHBTPair* pair, Double_t& x, Double_t&y ,Double_t& z) const
83 //Returns values to be histogrammed
85 x = pair->GetDeltaPx();
86 y = pair->GetDeltaPy();
87 z = pair->GetDeltaPz();
96 /******************************************************************/
98 TH1* AliHBTTwoTrackEffFctnPxPyPz::GetResult()
100 //returns ratio of numerator and denominator
102 fRatio = GetRatio(Scale());
106 /******************************************************************/
107 /******************************************************************/
108 /******************************************************************/
109 ClassImp(AliHBTTwoTrackEffFctnPtThetaPhi)
111 AliHBTTwoTrackEffFctnPtThetaPhi::AliHBTTwoTrackEffFctnPtThetaPhi(Int_t nXbins, Double_t maxXval, Double_t minXval,
112 Int_t nYbins, Double_t maxYval, Double_t minYval,
113 Int_t nZbins, Double_t maxZval, Double_t minZval):
114 AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
118 fWriteNumAndDen = kTRUE;
119 Rename("tteffptthetaphi","P_{t} \\theta \\phi Two Track Efficiency Function");
122 fNumerator->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
123 fNumerator->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
124 fNumerator->GetZaxis()->SetTitle("\\Delta \\phi [rad]");
129 fDenominator->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
130 fDenominator->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
131 fDenominator->GetZaxis()->SetTitle("\\Delta \\phi [rad]");
134 /******************************************************************/
136 void AliHBTTwoTrackEffFctnPtThetaPhi::GetValues(AliHBTPair* pair, Double_t& x, Double_t&y ,Double_t& z) const
138 //Returns values to be histogrammed
140 x = pair->GetDeltaPt();
141 y = pair->GetDeltaTheta();
142 z = pair->GetDeltaPhi();
152 /******************************************************************/
154 TH1* AliHBTTwoTrackEffFctnPtThetaPhi::GetResult()
156 //returns ratio of numerator and denominator
158 fRatio = GetRatio(Scale());
161 fRatio->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
162 fRatio->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
163 fRatio->GetZaxis()->SetTitle("\\Delta \\phi [rad]");