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)
83 //Returns values to be histogrammed
85 x = pair->GetDeltaPx();
86 y = pair->GetDeltaPy();
87 z = pair->GetDeltaPz();
89 /******************************************************************/
91 TH1* AliHBTTwoTrackEffFctnPxPyPz::GetResult()
93 //returns ratio of numerator and denominator
95 fRatio = GetRatio(Scale());
99 /******************************************************************/
100 /******************************************************************/
101 /******************************************************************/
102 ClassImp(AliHBTTwoTrackEffFctnPtThetaPhi)
104 AliHBTTwoTrackEffFctnPtThetaPhi::AliHBTTwoTrackEffFctnPtThetaPhi(Int_t nXbins, Double_t maxXval, Double_t minXval,
105 Int_t nYbins, Double_t maxYval, Double_t minYval,
106 Int_t nZbins, Double_t maxZval, Double_t minZval):
107 AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
111 fWriteNumAndDen = kTRUE;
112 Rename("tteffptthetaphi","P_{t} \\theta \\phi Two Track Efficiency Function");
115 fNumerator->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
116 fNumerator->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
117 fNumerator->GetZaxis()->SetTitle("\\Delta \\phi [rad]");
122 fDenominator->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
123 fDenominator->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
124 fDenominator->GetZaxis()->SetTitle("\\Delta \\phi [rad]");
127 /******************************************************************/
129 void AliHBTTwoTrackEffFctnPtThetaPhi::GetValues(AliHBTPair* pair, Double_t& x, Double_t&y ,Double_t& z)
131 //Returns values to be histogrammed
133 x = pair->GetDeltaPt();
134 y = pair->GetDeltaTheta();
135 z = pair->GetDeltaPhi();
137 /******************************************************************/
139 TH1* AliHBTTwoTrackEffFctnPtThetaPhi::GetResult()
141 //returns ratio of numerator and denominator
143 fRatio = GetRatio(Scale());
146 fRatio->GetXaxis()->SetTitle("\\Delta P_{t} [GeV]");
147 fRatio->GetYaxis()->SetTitle("\\Delta \\theta [rad]");
148 fRatio->GetZaxis()->SetTitle("\\Delta \\phi [rad]");