]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG1/AliComparisonRes.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Fixing the decoding of regional header bits.
[u/mrichter/AliRoot.git] / PWG1 / AliComparisonRes.cxx
CommitLineData
96b3109a 1//------------------------------------------------------------------------------
2// Implementation of AliComparisonRes class. It keeps information from
3// comparison of reconstructed and MC particle tracks. In addtion,
4// it keeps selection cuts used during comparison. The comparison
5// information is stored in the ROOT histograms. Analysis of these
6// histograms can be done by using Analyse() class function. The result of
3baa4bfd 7// the analysis (histograms/graphs) are stored in the folder which is
8// a data member of AliComparisonRes.
9//
96b3109a 10// Author: J.Otwinowski 04/02/2008
11//------------------------------------------------------------------------------
12
13/*
3baa4bfd 14
15 // after running comparison task, read the file, and get component
26e10df5 16 gROOT->LoadMacro("$ALICE_ROOT/PWG1/Macros/LoadMyLibs.C");
17 LoadMyLibs();
18
96b3109a 19 TFile f("Output.root");
3baa4bfd 20 AliComparisonRes * compObj = (AliComparisonRes*)f.Get("AliComparisonRes");
21
22 // analyse comparison data
23 compObj->Analyse();
24
25 // the output histograms/graphs will be stored in the folder "folderRes"
26 compObj->GetAnalysisFolder()->ls("*");
27
28 // user can save whole comparison object (or only folder with anlysed histograms)
29 // in the seperate output file (e.g.)
30 TFile fout("Analysed_Res.root","recreate");
31 compObj->Write(); // compObj->GetAnalysisFolder()->Write();
32 fout.Close();
96b3109a 33
96b3109a 34*/
35
36#include <iostream>
37
38#include "TFile.h"
39#include "TCint.h"
40#include "TH3F.h"
41#include "TH2F.h"
42#include "TF1.h"
43#include "TProfile.h"
44#include "TProfile2D.h"
45#include "TGraph2D.h"
46#include "TCanvas.h"
47#include "TGraph.h"
48
49#include "AliESDEvent.h"
50#include "AliESD.h"
51#include "AliESDfriend.h"
52#include "AliESDfriendTrack.h"
53#include "AliESDVertex.h"
54#include "AliRecInfoCuts.h"
55#include "AliMCInfoCuts.h"
56#include "AliLog.h"
57#include "AliTracker.h"
58
59#include "AliMathBase.h"
60#include "AliTreeDraw.h"
61
62#include "AliMCInfo.h"
63#include "AliESDRecInfo.h"
64#include "AliComparisonRes.h"
65
66using namespace std;
67
68ClassImp(AliComparisonRes)
69
70//_____________________________________________________________________________
71AliComparisonRes::AliComparisonRes():
3baa4bfd 72 AliComparisonObject("AliComparisonRes"),
96b3109a 73
74 // Resolution
75 fPtResolLPT(0), // pt resolution - low pt
76 fPtResolHPT(0), // pt resolution - high pt
77 fPtPullLPT(0), // pt resolution - low pt
78 fPtPullHPT(0), // pt resolution - high pt
decf0997 79 fPhiResolTan(0), //-> angular resolution
80 fTanResolTan(0), //-> angular resolution
81 fPhiPullTan(0), //-> angular resolution
82 fTanPullTan(0), //-> angular resolution
83
96b3109a 84 //
85 // Resolution constrained param
86 //
87 fCPhiResolTan(0), // angular resolution - constrained
88 fCTanResolTan(0), // angular resolution - constrained
89 fCPtResolTan(0), // pt resolution - constrained
90 fCPhiPullTan(0), // angular resolution - constrained
91 fCTanPullTan(0), // angular resolution - constrained
92 fCPtPullTan(0), // pt resolution - constrained
93
94 //
95 // Parametrisation histograms
96 //
97
3baa4bfd 98 f1Pt2ResolS1PtTPC(0),
99 f1Pt2ResolS1PtTPCITS(0),
100 fYResolS1PtTPC(0),
101 fYResolS1PtTPCITS(0),
102 fZResolS1PtTPC(0),
103 fZResolS1PtTPCITS(0),
104 fPhiResolS1PtTPC(0),
105 fPhiResolS1PtTPCITS(0),
106 fThetaResolS1PtTPC(0),
107 fThetaResolS1PtTPCITS(0),
96b3109a 108
109 // constrained
3baa4bfd 110 fC1Pt2ResolS1PtTPC(0),
111 fC1Pt2ResolS1PtTPCITS(0),
112 fCYResolS1PtTPC(0),
113 fCYResolS1PtTPCITS(0),
114 fCZResolS1PtTPC(0),
115 fCZResolS1PtTPCITS(0),
116 fCPhiResolS1PtTPC(0),
117 fCPhiResolS1PtTPCITS(0),
118 fCThetaResolS1PtTPC(0),
119 fCThetaResolS1PtTPCITS(0),
96b3109a 120
121 // vertex
122 fVertex(0),
123
124 // Cuts
125 fCutsRC(0),
3baa4bfd 126 fCutsMC(0),
127
128 // histogram folder
129 fAnalysisFolder(0)
96b3109a 130{
3baa4bfd 131 Init();
96b3109a 132
133 // vertex (0,0,0)
134 fVertex = new AliESDVertex();
135 fVertex->SetXv(0.0);
136 fVertex->SetYv(0.0);
137 fVertex->SetZv(0.0);
138}
139
140//_____________________________________________________________________________
141AliComparisonRes::~AliComparisonRes(){
142
143 // Resolution histograms
144 if(fPtResolLPT) delete fPtResolLPT; fPtResolLPT=0;
145 if(fPtResolHPT) delete fPtResolHPT; fPtResolHPT=0;
146 if(fPtPullLPT) delete fPtPullLPT; fPtPullLPT=0;
147 if(fPtPullHPT) delete fPtPullHPT; fPtPullHPT=0;
148
decf0997 149 if(fPhiResolTan) delete fPhiResolTan; fPhiResolTan=0;
150 if(fTanResolTan) delete fTanResolTan; fTanResolTan=0;
151 if(fPhiPullTan) delete fPhiPullTan; fPhiPullTan=0;
152 if(fTanPullTan) delete fTanPullTan; fTanPullTan=0;
153
96b3109a 154 // Resolution histograms (constrained param)
155 if(fCPhiResolTan) delete fCPhiResolTan; fCPhiResolTan=0;
156 if(fCTanResolTan) delete fCTanResolTan; fCTanResolTan=0;
157 if(fCPtResolTan) delete fCPtResolTan; fCPtResolTan=0;
158 if(fCPhiPullTan) delete fCPhiPullTan; fCPhiPullTan=0;
159 if(fCTanPullTan) delete fCTanPullTan; fCTanPullTan=0;
160 if(fCPtPullTan) delete fCPtPullTan; fCPtPullTan=0;
161
162 // Parametrisation histograms
163 //
3baa4bfd 164 if(f1Pt2ResolS1PtTPC) delete f1Pt2ResolS1PtTPC; f1Pt2ResolS1PtTPC=0;
165 if(f1Pt2ResolS1PtTPCITS) delete f1Pt2ResolS1PtTPCITS; f1Pt2ResolS1PtTPCITS=0;
166 if(fYResolS1PtTPC) delete fYResolS1PtTPC; fYResolS1PtTPC=0;
167 if(fYResolS1PtTPCITS) delete fYResolS1PtTPCITS; fYResolS1PtTPCITS=0;
168 if(fZResolS1PtTPC) delete fZResolS1PtTPC; fZResolS1PtTPC=0;
169 if(fZResolS1PtTPCITS) delete fZResolS1PtTPCITS; fZResolS1PtTPCITS=0;
170 if(fPhiResolS1PtTPC) delete fPhiResolS1PtTPC; fPhiResolS1PtTPC=0;
171 if(fPhiResolS1PtTPCITS) delete fPhiResolS1PtTPCITS; fPhiResolS1PtTPCITS=0;
172 if(fThetaResolS1PtTPC) delete fThetaResolS1PtTPC; fThetaResolS1PtTPC=0;
173 if(fThetaResolS1PtTPCITS) delete fThetaResolS1PtTPCITS; fThetaResolS1PtTPCITS=0;
96b3109a 174
175 // constrained
3baa4bfd 176 if(fC1Pt2ResolS1PtTPC) delete fC1Pt2ResolS1PtTPC; fC1Pt2ResolS1PtTPC=0;
177 if(fC1Pt2ResolS1PtTPCITS) delete fC1Pt2ResolS1PtTPCITS; fC1Pt2ResolS1PtTPCITS=0;
178 if(fCYResolS1PtTPC) delete fCYResolS1PtTPC; fCYResolS1PtTPC=0;
179 if(fCYResolS1PtTPCITS) delete fCYResolS1PtTPCITS; fCYResolS1PtTPCITS=0;
180 if(fCZResolS1PtTPC) delete fCZResolS1PtTPC; fCZResolS1PtTPC=0;
181 if(fCZResolS1PtTPCITS) delete fCZResolS1PtTPCITS; fCZResolS1PtTPCITS=0;
182 if(fCPhiResolS1PtTPC) delete fCPhiResolS1PtTPC; fCPhiResolS1PtTPC=0;
183 if(fCPhiResolS1PtTPCITS) delete fCPhiResolS1PtTPCITS; fCPhiResolS1PtTPCITS=0;
184 if(fCThetaResolS1PtTPC) delete fCThetaResolS1PtTPC; fCThetaResolS1PtTPC=0;
185 if(fCThetaResolS1PtTPCITS) delete fCThetaResolS1PtTPCITS; fCThetaResolS1PtTPCITS=0;
96b3109a 186
187 if(fVertex) delete fVertex; fVertex=0;
188
3baa4bfd 189 if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
96b3109a 190}
191
192//_____________________________________________________________________________
3baa4bfd 193void AliComparisonRes::Init(){
96b3109a 194
195 // Init histograms
196 fCPhiResolTan = new TH2F("CPhiResolTan","CPhiResolTan",50, -2,2,200,-0.025,0.025);
197 fCPhiResolTan->SetXTitle("tan(#theta)");
198 fCPhiResolTan->SetYTitle("#Delta#phi");
199
200 fCTanResolTan = new TH2F("CTanResolTan","CTanResolTan",50, -2,2,200,-0.025,0.025);
201 fCTanResolTan->SetXTitle("tan(#theta)");
202 fCTanResolTan->SetYTitle("#Delta#theta");
203
204 fCPtResolTan=new TH2F("CPtResol","CPtResol",50, -2,2,200,-0.2,0.2);
205 fCPtResolTan->SetXTitle("Tan(#theta)");
206 fCPtResolTan->SetYTitle("#Deltap_{t}/p_{t}");
207
208 fCPhiPullTan = new TH2F("CPhiPullTan","CPhiPullTan",50, -2,2,200,-5,5);
209 fCPhiPullTan->SetXTitle("Tan(#theta)");
210 fCPhiPullTan->SetYTitle("#Delta#phi/#Sigma");
211
212 fCTanPullTan = new TH2F("CTanPullTan","CTanPullTan",50, -2,2,200,-5,5);
213 fCTanPullTan->SetXTitle("Tan(#theta)");
214 fCTanPullTan->SetYTitle("#Delta#theta/#Sigma");
215
216 fCPtPullTan=new TH2F("CPtPull","CPtPull",50, -2,2,200,-5,5);
217 fCPtPullTan->SetXTitle("Tan(#theta)");
218 fCPtPullTan->SetYTitle("(1/mcp_{t}-1/p_{t})/#Sigma");
219
220 fPtResolLPT = new TH2F("Pt_resol_lpt","pt resol",10, 0.1,3,200,-0.2,0.2);
221 fPtResolLPT->SetXTitle("p_{t}");
222 fPtResolLPT->SetYTitle("#Deltap_{t}/p_{t}");
223
224 fPtResolHPT = new TH2F("Pt_resol_hpt","pt resol",10, 2,100,200,-0.3,0.3);
225 fPtResolHPT->SetXTitle("p_{t}");
226 fPtResolHPT->SetYTitle("#Deltap_{t}/p_{t}");
227
228 fPtPullLPT = new TH2F("Pt_pull_lpt","pt pull",10, 0.1,3,200,-6,6);
229 fPtPullLPT->SetXTitle("p_{t}");
230 fPtPullLPT->SetYTitle("#Deltap_{t}/#Sigma");
231
232 fPtPullHPT = new TH2F("Pt_pull_hpt","pt pull",10,2,100,200,-6,6);
233 fPtPullHPT->SetXTitle("p_{t}");
234 fPtPullHPT->SetYTitle("#Deltap_{t}/#Sigma");
decf0997 235
236 fPhiResolTan = new TH2F("PhiResolTan","PhiResolTan",50, -2,2,200,-0.025,0.025);
237 fPhiResolTan->SetXTitle("tan(#theta)");
238 fPhiResolTan->SetYTitle("#Delta#phi");
239
240 fTanResolTan = new TH2F("TanResolTan","TanResolTan",50, -2,2,200,-0.025,0.025);
241 fTanResolTan->SetXTitle("tan(#theta)");
242 fTanResolTan->SetYTitle("#Delta#theta");
243
244 fPhiPullTan = new TH2F("PhiPullTan","PhiPullTan",50, -2,2,200,-5,5);
245 fPhiPullTan->SetXTitle("Tan(#theta)");
246 fPhiPullTan->SetYTitle("#Delta#phi/#Sigma");
247
248 fTanPullTan = new TH2F("TanPullTan","TanPullTan",50, -2,2,200,-5,5);
249 fTanPullTan->SetXTitle("Tan(#theta)");
250 fTanPullTan->SetYTitle("#Delta#theta/#Sigma");
96b3109a 251
252 //
253 // Parametrisation histograms
254 //
255
3baa4bfd 256 f1Pt2ResolS1PtTPC = new TH2F("f1Pt2ResolS1PtTPC","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs sqrt(1/pt))",100,0,3,200,-0.010,0.010);
257 f1Pt2ResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
258 f1Pt2ResolS1PtTPC->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)");
96b3109a 259
3baa4bfd 260 f1Pt2ResolS1PtTPCITS = new TH2F("f1Pt2ResolS1PtTPCITS","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs sqrt(1/pt))",100,0,3,200,-0.010,0.010);
261 f1Pt2ResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
262 f1Pt2ResolS1PtTPCITS->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)");
96b3109a 263
3baa4bfd 264 fYResolS1PtTPC = new TH2F("fYResolS1PtTPC","fYResolS1PtTPC",100, 0,3,200,-1.0,1.0);
265 fYResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
266 fYResolS1PtTPC->SetYTitle("#DeltaY");
96b3109a 267
3baa4bfd 268 fYResolS1PtTPCITS = new TH2F("fYResolS1PtTPCITS","fYResolS1PtTPCITS",100, 0,3,200,-0.05,0.05);
269 fYResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
270 fYResolS1PtTPCITS->SetYTitle("#DeltaY");
96b3109a 271
3baa4bfd 272 fZResolS1PtTPC = new TH2F("fZResolS1PtTPC","fZResolS1PtTPC",100, 0,3,200,-1.0,1.0);
273 fZResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
274 fZResolS1PtTPC->SetYTitle("#DeltaZ");
96b3109a 275
3baa4bfd 276 fZResolS1PtTPCITS = new TH2F("fZResolS1PtTPCITS","fZResolS1PtTPCITS",100, 0,3,200,-0.05,0.05);
277 fZResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
278 fZResolS1PtTPCITS->SetYTitle("#DeltaZ");
96b3109a 279
3baa4bfd 280 fPhiResolS1PtTPC = new TH2F("fPhiResolS1PtTPC","fPhiResolS1PtTPC",100, 0,3,200,-0.025,0.025);
281 fPhiResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
282 fPhiResolS1PtTPC->SetYTitle("#Delta#phi");
96b3109a 283
3baa4bfd 284 fPhiResolS1PtTPCITS = new TH2F("fPhiResolS1PtTPCITS","fPhiResolS1PtTPCITS",100, 0,3,200,-0.01,0.01);
285 fPhiResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
286 fPhiResolS1PtTPCITS->SetYTitle("#Delta#phi");
96b3109a 287
3baa4bfd 288 fThetaResolS1PtTPC = new TH2F("fThetaResolS1PtTPC","fThetaResolS1PtTPC",100, 0,3,200,-0.025,0.025);
289 fThetaResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
290 fThetaResolS1PtTPC->SetYTitle("#Delta#theta");
96b3109a 291
3baa4bfd 292 fThetaResolS1PtTPCITS = new TH2F("fThetaResolS1PtTPCITS","fThetaResolS1PtTPCITS",100, 0,3,200,-0.01,0.01);
293 fThetaResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
294 fThetaResolS1PtTPCITS->SetYTitle("#Delta#theta");
96b3109a 295
296 // constrained
3baa4bfd 297 fC1Pt2ResolS1PtTPC = new TH2F("fC1Pt2ResolS1PtTPC","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs 1/pt)",100,0,3,200,-0.010,0.010);
298 fC1Pt2ResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
299 fC1Pt2ResolS1PtTPC->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)");
96b3109a 300
3baa4bfd 301 fC1Pt2ResolS1PtTPCITS = new TH2F("fC1Pt2ResolS1PtTPCITS","(1/mcpt-1/pt)/(1+1/mcpt)^2 vs 1/pt)",100,0,3,200,-0.010,0.010);
302 fC1Pt2ResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
303 fC1Pt2ResolS1PtTPCITS->SetYTitle("(1/mcp_{t}-1/p_{t})/(1+1/mcp_{t})^2)");
96b3109a 304
3baa4bfd 305 fCYResolS1PtTPC = new TH2F("fCYResolS1PtTPC","fCYResolS1PtTPC",100, 0,3,200,-1.0,1.0);
306 fCYResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
307 fCYResolS1PtTPC->SetYTitle("#DeltaY");
96b3109a 308
3baa4bfd 309 fCYResolS1PtTPCITS = new TH2F("fCYResolS1PtTPCITS","fCYResolS1PtTPCITS",100, 0,3,200,-0.01,0.01);
310 fCYResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
311 fCYResolS1PtTPCITS->SetYTitle("#DeltaY");
96b3109a 312
3baa4bfd 313 fCZResolS1PtTPC = new TH2F("fCZResolS1PtTPC","fCZResolS1PtTPC",100, 0,3,200,-1.0,1.0);
314 fCZResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
315 fCZResolS1PtTPC->SetYTitle("#DeltaZ");
96b3109a 316
3baa4bfd 317 fCZResolS1PtTPCITS = new TH2F("fCZResolS1PtTPCITS","fCZResolS1PtTPCITS",100, 0,3,200,-0.025,0.025);
318 fCZResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
319 fCZResolS1PtTPCITS->SetYTitle("#DeltaZ");
96b3109a 320
3baa4bfd 321 fCPhiResolS1PtTPC = new TH2F("fCPhiResolS1PtTPC","fCPhiResolS1PtTPC",100, 0,3,200,-0.025,0.025);
322 fCPhiResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
323 fCPhiResolS1PtTPC->SetYTitle("#Delta#phi");
96b3109a 324
3baa4bfd 325 fCPhiResolS1PtTPCITS = new TH2F("fCPhiResolS1PtTPCITS","fCPhiResolS1PtTPCITS",100, 0,3,200,-0.003,0.003);
326 fCPhiResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
327 fCPhiResolS1PtTPCITS->SetYTitle("#Delta#phi");
96b3109a 328
3baa4bfd 329 fCThetaResolS1PtTPC = new TH2F("fCThetaResolS1PtTPC","fCThetaResolS1PtTPC",100, 0,3,200,-0.025,0.025);
330 fCThetaResolS1PtTPC->SetXTitle("#sqrt{1/mcp_{t}}");
331 fCThetaResolS1PtTPC->SetYTitle("#Delta#theta");
96b3109a 332
3baa4bfd 333 fCThetaResolS1PtTPCITS = new TH2F("fCThetaResolS1PtTPCITS","fCThetaResolS1PtTPCITS",100, 0,3,200,-0.005,0.005);
334 fCThetaResolS1PtTPCITS->SetXTitle("#sqrt{1/mcp_{t}}");
335 fCThetaResolS1PtTPCITS->SetYTitle("#Delta#theta");
96b3109a 336
96b3109a 337 // Init cuts
338 if(!fCutsMC)
339 AliDebug(AliLog::kError, "ERROR: Cannot find AliMCInfoCuts object");
340 if(!fCutsRC)
341 AliDebug(AliLog::kError, "ERROR: Cannot find AliRecInfoCuts object");
3baa4bfd 342
343 // init folder
344 fAnalysisFolder = CreateFolder("folderRes","Analysis Resolution Folder");
96b3109a 345}
346
347//_____________________________________________________________________________
348void AliComparisonRes::Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC)
349{
350 // Fill resolution comparison information
351 AliExternalTrackParam *track = 0;
96b3109a 352 Double_t field = AliTracker::GetBz(); // nominal Bz field [kG]
353 Double_t kMaxD = 123456.0; // max distance
354
355 Int_t clusterITS[200];
356 Double_t dca[2], cov[3]; // dca_xy, dca_z, sigma_xy, sigma_xy_z, sigma_z
357
decf0997 358 Float_t deltaPt, pullPt, deltaPhi,pullPhi, deltaTan, pullTan, delta1Pt2, deltaY1Pt, deltaZ1Pt, deltaPhi1Pt, deltaTheta1Pt;
96b3109a 359 Float_t deltaPtTPC, pullPtTPC, deltaPhiTPC, deltaTanTPC, delta1Pt2TPC, deltaY1PtTPC, deltaZ1PtTPC, deltaPhi1PtTPC, deltaTheta1PtTPC;
360
361 Float_t mcpt = infoMC->GetParticle().Pt();
3baa4bfd 362 Float_t s1mcpt = TMath::Sqrt(1./infoMC->GetParticle().Pt());
decf0997 363 Float_t tantheta = TMath::Tan(infoMC->GetParticle().Theta()-TMath::Pi()*0.5);
96b3109a 364
365 // distance to Prim. vertex
366 const Double_t* dv = infoMC->GetVDist();
367 Bool_t isPrim = TMath::Sqrt(dv[0]*dv[0] + dv[1]*dv[1])<fCutsMC->GetMaxR() && TMath::Abs(dv[2])<fCutsMC->GetMaxVz();
368
369 // Check selection cuts
370 if (fCutsMC->IsPdgParticle(TMath::Abs(infoMC->GetParticle().GetPdgCode())) == kFALSE) return;
371 if (!isPrim) return;
96b3109a 372 if (infoRC->GetStatus(1)!=3) return; // TPC refit
373 if (!infoRC->GetESDtrack()) return;
374 if (infoRC->GetESDtrack()->GetTPCNcls()<fCutsRC->GetMinNClustersTPC()) return;
375
376 // calculate and set prim. vertex
377 fVertex->SetXv( infoMC->GetParticle().Vx() - dv[0] );
378 fVertex->SetYv( infoMC->GetParticle().Vy() - dv[1] );
379 fVertex->SetZv( infoMC->GetParticle().Vz() - dv[2] );
380
381 deltaPt= (mcpt-infoRC->GetESDtrack()->Pt())/mcpt;
382 pullPt= (1/mcpt-infoRC->GetESDtrack()->OneOverPt())/TMath::Sqrt(infoRC->GetESDtrack()->GetSigma1Pt2());
383 deltaPhi = TMath::ATan2(infoRC->GetESDtrack()->Py(),infoRC->GetESDtrack()->Px())-
384 TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px());
decf0997 385 pullPhi = deltaPhi/TMath::Sqrt(infoRC->GetESDtrack()->GetSigmaSnp2());
96b3109a 386
387 deltaTan = TMath::ATan2(infoRC->GetESDtrack()->Pz(),infoRC->GetESDtrack()->Pt())-
388 TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt());
decf0997 389 pullTan = deltaPhi/TMath::Sqrt(infoRC->GetESDtrack()->GetSigmaSnp2());
96b3109a 390
391 delta1Pt2 = (1/mcpt-infoRC->GetESDtrack()->OneOverPt())/TMath::Power(1+1/mcpt,2);
392 deltaY1Pt = (infoMC->GetParticle().Vy()-infoRC->GetESDtrack()->GetY()) / (0.2+1/mcpt);
393 deltaZ1Pt = (infoMC->GetParticle().Vz()-infoRC->GetESDtrack()->GetZ()) / (0.2+1/mcpt);
394 deltaPhi1Pt = deltaPhi / (0.1+1/mcpt);
395 deltaTheta1Pt = deltaTan / (0.1+1/mcpt);
396
96b3109a 397 // calculate track parameters at vertex
398 const AliExternalTrackParam *innerTPC = 0;
399 if ((innerTPC = infoRC->GetESDtrack()->GetTPCInnerParam()) != 0)
400 {
401 if ((track = new AliExternalTrackParam(*infoRC->GetESDtrack()->GetTPCInnerParam())) != 0 )
402 {
96b3109a 403 Bool_t bDCAStatus = track->PropagateToDCA(fVertex,field,kMaxD,dca,cov);
404
405 // Fill parametrisation histograms (only TPC track)
3baa4bfd 406 if(bDCAStatus)
96b3109a 407 {
408 deltaPtTPC= (mcpt-innerTPC->Pt())/mcpt;
409 pullPtTPC= (1/mcpt-innerTPC->OneOverPt())/TMath::Sqrt(innerTPC->GetSigma1Pt2());
410 deltaPhiTPC = TMath::ATan2(innerTPC->Py(),innerTPC->Px())-
411 TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px());
412
413 deltaTanTPC = TMath::ATan2(innerTPC->Pz(),innerTPC->Pt())-
414 TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt());
415
416 delta1Pt2TPC = (1/mcpt-innerTPC->OneOverPt())/TMath::Power(1+1/mcpt,2);
417 deltaY1PtTPC= (infoMC->GetParticle().Vy()-innerTPC->GetY()) / (0.2+1/mcpt);
418 deltaZ1PtTPC = (infoMC->GetParticle().Vz()-innerTPC->GetZ()) / (0.2+1/mcpt);
419 deltaPhi1PtTPC = deltaPhiTPC / (0.1+1/mcpt);
420 deltaTheta1PtTPC = deltaTanTPC / (0.1+1/mcpt);
421
3baa4bfd 422 f1Pt2ResolS1PtTPC->Fill(s1mcpt,delta1Pt2TPC);
423 fYResolS1PtTPC->Fill(s1mcpt,deltaY1PtTPC);
424 fZResolS1PtTPC->Fill(s1mcpt,deltaZ1PtTPC);
425 fPhiResolS1PtTPC->Fill(s1mcpt,deltaPhi1PtTPC);
426 fThetaResolS1PtTPC->Fill(s1mcpt,deltaTheta1PtTPC);
96b3109a 427 }
428 delete track;
429 }
430 }
431
432 // TPC and ITS (nb. of clusters >2) in the system
433 if(infoRC->GetESDtrack()->GetITSclusters(clusterITS)>2)
434 {
3baa4bfd 435 f1Pt2ResolS1PtTPCITS->Fill(s1mcpt,delta1Pt2);
436 fYResolS1PtTPCITS->Fill(s1mcpt,deltaY1Pt);
437 fZResolS1PtTPCITS->Fill(s1mcpt,deltaZ1Pt);
438 fPhiResolS1PtTPCITS->Fill(s1mcpt,deltaPhi1Pt);
439 fThetaResolS1PtTPCITS->Fill(s1mcpt,deltaTheta1Pt);
96b3109a 440 }
441
442 // Fill histograms
443 fPtResolLPT->Fill(mcpt,deltaPt);
444 fPtResolHPT->Fill(mcpt,deltaPt);
445 fPtPullLPT->Fill(mcpt,pullPt);
446 fPtPullHPT->Fill(mcpt,pullPt);
decf0997 447
448 fPhiResolTan->Fill(tantheta,deltaPhi);
449 fPhiPullTan->Fill(tantheta,pullPhi);
450 fTanResolTan->Fill(tantheta,deltaTan);
451 fTanPullTan->Fill(tantheta,pullTan);
452
96b3109a 453}
454
455//_____________________________________________________________________________
456void AliComparisonRes::ProcessConstrained(AliMCInfo* infoMC, AliESDRecInfo *infoRC)
457{
458 // Fill resolution comparison information (constarained parameters)
459 //
460 AliExternalTrackParam *track = 0;
96b3109a 461 Double_t field = AliTracker::GetBz(); // nominal Bz field [kG]
462 Double_t kMaxD = 123456.0; // max distance
463
464 Int_t clusterITS[200];
465 Double_t dca[2], cov[3]; // dca_xy, dca_z, sigma_xy, sigma_xy_z, sigma_z
466
467 Float_t deltaPt, pullPt, deltaPhi, pullPhi, deltaTan, pullTan, delta1Pt2, deltaY1Pt, deltaZ1Pt, deltaPhi1Pt, deltaTheta1Pt;
468 Float_t deltaPtTPC, pullPtTPC, deltaPhiTPC, deltaTanTPC, delta1Pt2TPC, deltaY1PtTPC, deltaZ1PtTPC, deltaPhi1PtTPC, deltaTheta1PtTPC;
469
470 Float_t mcpt = infoMC->GetParticle().Pt();
3baa4bfd 471 Float_t s1mcpt = TMath::Sqrt(1./infoMC->GetParticle().Pt());
96b3109a 472 Float_t tantheta = TMath::Tan(infoMC->GetParticle().Theta()-TMath::Pi()*0.5);
473
474 // distance to Prim. vertex
475 const Double_t* dv = infoMC->GetVDist();
476 Bool_t isPrim = TMath::Sqrt(dv[0]*dv[0] + dv[1]*dv[1])<fCutsMC->GetMaxR() && TMath::Abs(dv[2])<fCutsMC->GetMaxVz();
477
478 // Check selection cuts
479 if (fCutsMC->IsPdgParticle(TMath::Abs(infoMC->GetParticle().GetPdgCode())) == kFALSE) return;
480 if (!isPrim) return;
481 if (infoRC->GetStatus(1)!=3) return;
482 if (!infoRC->GetESDtrack()) return;
483 if (infoRC->GetESDtrack()->GetTPCNcls()<fCutsRC->GetMinNClustersTPC()) return;
484 if (!infoRC->GetESDtrack()->GetConstrainedParam()) return;
485
486// calculate and set prim. vertex
487 fVertex->SetXv( infoMC->GetParticle().Vx() - dv[0] );
488 fVertex->SetYv( infoMC->GetParticle().Vy() - dv[1] );
489 fVertex->SetZv( infoMC->GetParticle().Vz() - dv[2] );
490
491 // constrained parameters resolution
492 const AliExternalTrackParam * cparam = infoRC->GetESDtrack()->GetConstrainedParam();
493 deltaPt= (mcpt-cparam->Pt())/mcpt;
494 pullPt= (1/mcpt-cparam->OneOverPt())/TMath::Sqrt(cparam->GetSigma1Pt2());
495 deltaPhi = TMath::ATan2(cparam->Py(),cparam->Px())-
496 TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px());
497 pullPhi = deltaPhi/TMath::Sqrt(cparam->GetSigmaSnp2());
498 deltaTan = TMath::ATan2(cparam->Pz(),cparam->Pt())-TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt());
499 pullTan = deltaPhi/TMath::Sqrt(cparam->GetSigmaSnp2());
500
501
502 delta1Pt2 = (1/mcpt-cparam->OneOverPt())/TMath::Power(1+1/mcpt,2);
503
504 deltaY1Pt = (infoMC->GetParticle().Vy()-cparam->GetY()) / (0.2+1/mcpt);
505 deltaZ1Pt = (infoMC->GetParticle().Vz()-cparam->GetZ()) / (0.2+1/mcpt);
506 deltaPhi1Pt = deltaPhi / (0.1+1/mcpt);
507 deltaTheta1Pt = deltaTan / (0.1+1/mcpt);
508
509 // calculate track parameters at vertex
510 const AliExternalTrackParam *innerTPC = 0;
511 if ((innerTPC = infoRC->GetESDtrack()->GetTPCInnerParam()) != 0)
512 {
513 if ((track = new AliExternalTrackParam(*infoRC->GetESDtrack()->GetTPCInnerParam())) != 0 )
514 {
96b3109a 515 Bool_t bDCAStatus = track->PropagateToDCA(fVertex,field,kMaxD,dca,cov);
516
517 // Fill parametrisation histograms (only TPC track)
3baa4bfd 518 if(bDCAStatus)
96b3109a 519 {
520 deltaPtTPC= (mcpt-innerTPC->Pt())/mcpt;
521 pullPtTPC= (1/mcpt-innerTPC->OneOverPt())/TMath::Sqrt(innerTPC->GetSigma1Pt2());
522 deltaPhiTPC = TMath::ATan2(innerTPC->Py(),innerTPC->Px())-
523 TMath::ATan2(infoMC->GetParticle().Py(),infoMC->GetParticle().Px());
524
525 deltaTanTPC = TMath::ATan2(innerTPC->Pz(),innerTPC->Pt())-
526 TMath::ATan2(infoMC->GetParticle().Pz(),infoMC->GetParticle().Pt());
527
528 delta1Pt2TPC = (1/mcpt-innerTPC->OneOverPt())/TMath::Power(1+1/mcpt,2);
529 deltaY1PtTPC= (infoMC->GetParticle().Vy()-innerTPC->GetY()) / (0.2+1/mcpt);
530 deltaZ1PtTPC = (infoMC->GetParticle().Vz()-innerTPC->GetZ()) / (0.2+1/mcpt);
531 deltaPhi1PtTPC = deltaPhiTPC / (0.1+1/mcpt);
532 deltaTheta1PtTPC = deltaTanTPC / (0.1+1/mcpt);
533
3baa4bfd 534 fC1Pt2ResolS1PtTPC->Fill(s1mcpt,delta1Pt2TPC);
535 fCYResolS1PtTPC->Fill(s1mcpt,deltaY1PtTPC);
536 fCZResolS1PtTPC->Fill(s1mcpt,deltaZ1PtTPC);
537 fCPhiResolS1PtTPC->Fill(s1mcpt,deltaPhi1PtTPC);
538 fCThetaResolS1PtTPC->Fill(s1mcpt,deltaTheta1PtTPC);
96b3109a 539 }
540 delete track;
541 }
542 }
543
544 // TPC and ITS (nb. of clusters >2) in the system
545 if(infoRC->GetESDtrack()->GetITSclusters(clusterITS)>2)
546 {
3baa4bfd 547 fC1Pt2ResolS1PtTPCITS->Fill(s1mcpt,delta1Pt2);
548 fCYResolS1PtTPCITS->Fill(s1mcpt,deltaY1Pt);
549 fCZResolS1PtTPCITS->Fill(s1mcpt,deltaZ1Pt);
550 fCPhiResolS1PtTPCITS->Fill(s1mcpt,deltaPhi1Pt);
551 fCThetaResolS1PtTPCITS->Fill(s1mcpt,deltaTheta1Pt);
96b3109a 552 }
553
554 // Fill histograms
555 fCPtResolTan->Fill(tantheta,deltaPt);
556 fCPtPullTan->Fill(tantheta,pullPt);
557 fCPhiResolTan->Fill(tantheta,deltaPhi);
558 fCPhiPullTan->Fill(tantheta,pullPhi);
559 fCTanResolTan->Fill(tantheta,deltaTan);
560 fCTanPullTan->Fill(tantheta,pullTan);
561}
562
563//_____________________________________________________________________________
564void AliComparisonRes::Exec(AliMCInfo* infoMC, AliESDRecInfo *infoRC){
565
566 // Process comparison information
567 Process(infoMC,infoRC);
568 ProcessConstrained(infoMC,infoRC);
569}
570
571//_____________________________________________________________________________
572TH1F* AliComparisonRes::MakeResol(TH2F * his, Int_t integ, Bool_t type){
573 // Create resolution histograms
574
575 TH1F *hisr, *hism;
576 if (!gPad) new TCanvas;
577 hisr = AliTreeDraw::CreateResHistoI(his,&hism,integ);
578 if (type) return hism;
579 else
580 return hisr;
581}
582
583//_____________________________________________________________________________
584void AliComparisonRes::Analyse(){
3baa4bfd 585 // Analyse comparison information and store output histograms
586 // in the folder "folderRes"
587 //
588
589 TH1::AddDirectory(kFALSE);
590
96b3109a 591 AliComparisonRes * comp=this;
592 TH1F *hiss=0;
b4126c69 593 TObjArray *aFolderObj = new TObjArray;
3baa4bfd 594
595 // write results in the folder
96b3109a 596
597 TCanvas * c = new TCanvas("Phi resol Tan","Phi resol Tan");
598 c->cd();
599 //
600 hiss = comp->MakeResol(comp->fCPtResolTan,1,0);
601 hiss->SetXTitle("Tan(#theta)");
602 hiss->SetYTitle("#sigmap_{t}/p_{t}");
603 hiss->Draw();
3baa4bfd 604 hiss->SetName("CptResolTan");
605
b4126c69 606 aFolderObj->Add(hiss);
3baa4bfd 607
96b3109a 608 //
decf0997 609 hiss = comp->MakeResol(comp->fPtResolLPT,1,0);
610 hiss->SetXTitle("p_{t} (GeV)");
611 hiss->SetYTitle("#sigmap_{t}/p_{t}");
612 hiss->Draw();
613 hiss->SetName("PtResolPt");
614
615 aFolderObj->Add(hiss);
616
617 //
618 hiss = comp->MakeResol(comp->fPtResolLPT,1,1);
619 hiss->SetXTitle("p_{t} (GeV)");
620 hiss->SetYTitle("mean #Deltap_{t}/p_{t}");
621 hiss->Draw();
622 hiss->SetName("PtMeanPt");
623
624 aFolderObj->Add(hiss);
625
626 //
627 hiss = comp->MakeResol(comp->fPhiResolTan,1,0);
96b3109a 628 hiss->SetXTitle("Tan(#theta)");
629 hiss->SetYTitle("#sigma#phi (rad)");
630 hiss->Draw();
3baa4bfd 631 hiss->SetName("PhiResolTan");
632
b4126c69 633 aFolderObj->Add(hiss);
96b3109a 634 //
decf0997 635 hiss = comp->MakeResol(comp->fTanResolTan,1,0);
96b3109a 636 hiss->SetXTitle("Tan(#theta)");
637 hiss->SetYTitle("#sigma#theta (rad)");
638 hiss->Draw();
3baa4bfd 639 hiss->SetName("ThetaResolTan");
640
decf0997 641 aFolderObj->Add(hiss);
642
643 //
644 hiss = comp->MakeResol(comp->fPhiResolTan,1,1);
645 hiss->SetXTitle("Tan(#theta)");
646 hiss->SetYTitle("mean #Delta#phi (rad)");
647 hiss->Draw();
648 hiss->SetName("PhiMeanTan");
649
650 aFolderObj->Add(hiss);
651 //
652 hiss = comp->MakeResol(comp->fTanResolTan,1,1);
653 hiss->SetXTitle("Tan(#theta)");
654 hiss->SetYTitle("mean #Delta#theta (rad)");
655 hiss->Draw();
656 hiss->SetName("ThetaMeanTan");
657
658 aFolderObj->Add(hiss);
659
660 //
661 hiss = comp->MakeResol(comp->fCPhiResolTan,1,0);
662 hiss->SetXTitle("Tan(#theta)");
663 hiss->SetYTitle("#sigma#phi (rad)");
664 hiss->Draw();
665 hiss->SetName("CPhiResolTan");
666
667 aFolderObj->Add(hiss);
668 //
669 hiss = comp->MakeResol(comp->fCTanResolTan,1,0);
670 hiss->SetXTitle("Tan(#theta)");
671 hiss->SetYTitle("#sigma#theta (rad)");
672 hiss->Draw();
673 hiss->SetName("CThetaResolTan");
674
b4126c69 675 aFolderObj->Add(hiss);
96b3109a 676 //
677 hiss = comp->MakeResol(comp->fCPtPullTan,1,0);
678 hiss->SetXTitle("Tan(#theta)");
679 hiss->SetYTitle("1/mcp_{t}-1/p_{t}/#Sigma(1/p_{t})");
680 hiss->Draw();
3baa4bfd 681 hiss->SetName("CptPullTan");
682
b4126c69 683 aFolderObj->Add(hiss);
96b3109a 684 //
3baa4bfd 685 hiss = comp->MakeResol(comp->fC1Pt2ResolS1PtTPC,1,0);
b4126c69 686 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 687 hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2");
688 hiss->Draw();
3baa4bfd 689 hiss->SetName("C1Pt2ResolS1PtTPC");
690
b4126c69 691 aFolderObj->Add(hiss);
96b3109a 692
3baa4bfd 693 hiss = comp->MakeResol(comp->fC1Pt2ResolS1PtTPCITS,1,0);
b4126c69 694 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 695 hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2");
696 hiss->Draw();
3baa4bfd 697 hiss->SetName("C1Pt2ResolS1PtTPCITS");
698
b4126c69 699 aFolderObj->Add(hiss);
96b3109a 700 //
3baa4bfd 701 hiss = comp->MakeResol(comp->fCYResolS1PtTPC,1,0);
b4126c69 702 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 703 hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})");
704 hiss->Draw();
3baa4bfd 705 hiss->SetName("CYResolS1PtTPC");
706
b4126c69 707 aFolderObj->Add(hiss);
96b3109a 708
3baa4bfd 709 hiss = comp->MakeResol(comp->fCYResolS1PtTPCITS,1,0);
b4126c69 710 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 711 hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})");
712 hiss->Draw();
3baa4bfd 713 hiss->SetName("CYResolS1PtTPCITS");
714
b4126c69 715 aFolderObj->Add(hiss);
96b3109a 716 //
3baa4bfd 717 hiss = comp->MakeResol(comp->fCZResolS1PtTPC,1,0);
b4126c69 718 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 719 hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})");
720 hiss->Draw();
3baa4bfd 721 hiss->SetName("CZResolS1PtTPC");
722
b4126c69 723 aFolderObj->Add(hiss);
96b3109a 724
3baa4bfd 725 hiss = comp->MakeResol(comp->fCZResolS1PtTPCITS,1,0);
b4126c69 726 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 727 hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})");
728 hiss->Draw();
3baa4bfd 729 hiss->SetName("CZResolS1PtTPCITS");
730
b4126c69 731 aFolderObj->Add(hiss);
96b3109a 732 //
3baa4bfd 733 hiss = comp->MakeResol(comp->fCPhiResolS1PtTPC,1,0);
b4126c69 734 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 735 hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})");
736 hiss->Draw();
3baa4bfd 737 hiss->SetName("CPhiResolS1PtTPC");
738
b4126c69 739 aFolderObj->Add(hiss);
96b3109a 740
3baa4bfd 741 hiss = comp->MakeResol(comp->fCPhiResolS1PtTPCITS,1,0);
b4126c69 742 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 743 hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})");
744 hiss->Draw();
3baa4bfd 745 hiss->SetName("CPhiResolS1PtTPCITS");
746
b4126c69 747 aFolderObj->Add(hiss);
96b3109a 748 //
3baa4bfd 749 hiss = comp->MakeResol(comp->fCThetaResolS1PtTPC,1,0);
b4126c69 750 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 751 hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})");
752 hiss->Draw();
3baa4bfd 753 hiss->SetName("CThetaResolS1PtTPC");
754
b4126c69 755 aFolderObj->Add(hiss);
96b3109a 756
3baa4bfd 757 hiss = comp->MakeResol(comp->fCThetaResolS1PtTPCITS,1,0);
b4126c69 758 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 759 hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})");
760 hiss->Draw();
3baa4bfd 761 hiss->SetName("CThetaResolS1PtTPCITS");
762
b4126c69 763 aFolderObj->Add(hiss);
96b3109a 764
765 //
3baa4bfd 766 hiss = comp->MakeResol(comp->f1Pt2ResolS1PtTPC,1,0);
b4126c69 767 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 768 hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2");
769 hiss->Draw();
3baa4bfd 770 hiss->SetName("OnePt2ResolS1PtTPC");
771
b4126c69 772 aFolderObj->Add(hiss);
96b3109a 773
3baa4bfd 774 hiss = comp->MakeResol(comp->f1Pt2ResolS1PtTPCITS,1,0);
b4126c69 775 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 776 hiss->SetYTitle("1/mcp_{t}-1/p_{t}/(1+1/p_{t})^2");
777 hiss->Draw();
3baa4bfd 778 hiss->SetName("OnePt2ResolS1PtTPCITS");
779
b4126c69 780 aFolderObj->Add(hiss);
96b3109a 781 //
3baa4bfd 782 hiss = comp->MakeResol(comp->fYResolS1PtTPC,1,0);
b4126c69 783 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 784 hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})");
785 hiss->Draw();
3baa4bfd 786 hiss->SetName("YResolS1PtTPC");
787
b4126c69 788 aFolderObj->Add(hiss);
96b3109a 789
3baa4bfd 790 hiss = comp->MakeResol(comp->fYResolS1PtTPCITS,1,0);
b4126c69 791 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 792 hiss->SetYTitle("(mcy-y)/(0.2+1/mcp_{t})");
793 hiss->Draw();
3baa4bfd 794 hiss->SetName("YResolS1PtTPCITS");
795
b4126c69 796 aFolderObj->Add(hiss);
96b3109a 797 //
3baa4bfd 798 hiss = comp->MakeResol(comp->fZResolS1PtTPC,1,0);
b4126c69 799 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 800 hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})");
801 hiss->Draw();
3baa4bfd 802 hiss->SetName("ZResolS1PtTPC");
803
b4126c69 804 aFolderObj->Add(hiss);
96b3109a 805
3baa4bfd 806 hiss = comp->MakeResol(comp->fZResolS1PtTPCITS,1,0);
b4126c69 807 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 808 hiss->SetYTitle("(mcz-z)/(0.2+1/mcp_{t})");
809 hiss->Draw();
3baa4bfd 810 hiss->SetName("ZResolS1PtTPCITS");
811
b4126c69 812 aFolderObj->Add(hiss);
96b3109a 813 //
3baa4bfd 814 hiss = comp->MakeResol(comp->fPhiResolS1PtTPC,1,0);
b4126c69 815 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 816 hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})");
817 hiss->Draw();
3baa4bfd 818 hiss->SetName("PhiResolS1PtTPC");
819
b4126c69 820 aFolderObj->Add(hiss);
96b3109a 821
3baa4bfd 822 hiss = comp->MakeResol(comp->fPhiResolS1PtTPCITS,1,0);
b4126c69 823 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 824 hiss->SetYTitle("(mc#phi-#phi)/(0.1+1/mcp_{t})");
825 hiss->Draw();
3baa4bfd 826 hiss->SetName("PhiResolS1PtTPCITS");
827
b4126c69 828 aFolderObj->Add(hiss);
96b3109a 829 //
3baa4bfd 830 hiss = comp->MakeResol(comp->fThetaResolS1PtTPC,1,0);
b4126c69 831 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 832 hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})");
833 hiss->Draw();
3baa4bfd 834 hiss->SetName("ThetaResolS1PtTPC");
835
b4126c69 836 aFolderObj->Add(hiss);
96b3109a 837
3baa4bfd 838 hiss = comp->MakeResol(comp->fThetaResolS1PtTPCITS,1,0);
b4126c69 839 hiss->SetXTitle("#sqrt(1/mcp_{t})");
96b3109a 840 hiss->SetYTitle("(mc#theta-#theta)/(0.1+1/mcp_{t})");
841 hiss->Draw();
3baa4bfd 842 hiss->SetName("ThetaResolS1PtTPCITS");
843
b4126c69 844 aFolderObj->Add(hiss);
845
846 // export objects to analysis folder
847 fAnalysisFolder = ExportToFolder(aFolderObj);
848
849 // delete only TObjArray
850 if(aFolderObj) delete aFolderObj;
851}
852
853//_____________________________________________________________________________
854TFolder* AliComparisonRes::ExportToFolder(TObjArray * array)
855{
856 // recreate folder avery time and export objects to new one
857 //
858 AliComparisonRes * comp=this;
859 TFolder *folder = comp->GetAnalysisFolder();
860
861 TString name, title;
862 TFolder *newFolder = 0;
863 Int_t i = 0;
864 Int_t size = array->GetSize();
865
866 if(folder) {
867 // get name and title from old folder
868 name = folder->GetName();
869 title = folder->GetTitle();
870
871 // delete old one
872 delete folder;
873
874 // create new one
875 newFolder = CreateFolder(name.Data(),title.Data());
876 newFolder->SetOwner();
877
878 // add objects to folder
879 while(i < size) {
880 newFolder->Add(array->At(i));
881 i++;
882 }
883 }
96b3109a 884
b4126c69 885return newFolder;
96b3109a 886}
887
888//_____________________________________________________________________________
889Long64_t AliComparisonRes::Merge(TCollection* list)
890{
891 // Merge list of objects (needed by PROOF)
892
893 if (!list)
894 return 0;
895
896 if (list->IsEmpty())
897 return 1;
898
899 TIterator* iter = list->MakeIterator();
900 TObject* obj = 0;
901
902 // collection of generated histograms
903 Int_t count=0;
904 while((obj = iter->Next()) != 0)
905 {
906 AliComparisonRes* entry = dynamic_cast<AliComparisonRes*>(obj);
907 if (entry == 0) continue;
908
909 fPtResolLPT->Add(entry->fPtResolLPT);
910 fPtResolHPT->Add(entry->fPtResolHPT);
911 fPtPullLPT->Add(entry->fPtPullLPT);
912 fPtPullHPT->Add(entry->fPtPullHPT);
decf0997 913 fPhiResolTan->Add(entry->fPhiResolTan);
914 fTanResolTan->Add(entry->fTanResolTan);
915 fPhiPullTan->Add(entry->fPhiPullTan);
916 fTanPullTan->Add(entry->fTanPullTan);
96b3109a 917
918 // Histograms for 1/pt parameterisation
3baa4bfd 919 f1Pt2ResolS1PtTPC->Add(entry->f1Pt2ResolS1PtTPC);
920 fYResolS1PtTPC->Add(entry->fYResolS1PtTPC);
921 fZResolS1PtTPC->Add(entry->fZResolS1PtTPC);
922 fPhiResolS1PtTPC->Add(entry->fPhiResolS1PtTPC);
923 fThetaResolS1PtTPC->Add(entry->fThetaResolS1PtTPC);
924
925 f1Pt2ResolS1PtTPCITS->Add(entry->f1Pt2ResolS1PtTPCITS);
926 fYResolS1PtTPCITS->Add(entry->fYResolS1PtTPCITS);
927 fZResolS1PtTPCITS->Add(entry->fZResolS1PtTPCITS);
928 fPhiResolS1PtTPCITS->Add(entry->fPhiResolS1PtTPCITS);
929 fThetaResolS1PtTPCITS->Add(entry->fThetaResolS1PtTPCITS);
96b3109a 930
931 // Resolution histograms (constrained param)
932 fCPhiResolTan->Add(entry->fCPhiResolTan);
933 fCTanResolTan->Add(entry->fCTanResolTan);
934 fCPtResolTan->Add(entry->fCPtResolTan);
935 fCPhiPullTan->Add(entry->fCPhiPullTan);
936 fCTanPullTan->Add(entry->fCTanPullTan);
937 fCPtPullTan->Add(entry->fCPtPullTan);
938
939 // Histograms for 1/pt parameterisation (constrained)
3baa4bfd 940 fC1Pt2ResolS1PtTPC->Add(entry->fC1Pt2ResolS1PtTPC);
941 fCYResolS1PtTPC->Add(entry->fCYResolS1PtTPC);
942 fCZResolS1PtTPC->Add(entry->fCZResolS1PtTPC);
943 fCPhiResolS1PtTPC->Add(entry->fCPhiResolS1PtTPC);
944 fCThetaResolS1PtTPC->Add(entry->fCThetaResolS1PtTPC);
945
946 fC1Pt2ResolS1PtTPCITS->Add(entry->fC1Pt2ResolS1PtTPCITS);
947 fCYResolS1PtTPCITS->Add(entry->fCYResolS1PtTPCITS);
948 fCZResolS1PtTPCITS->Add(entry->fCZResolS1PtTPCITS);
949 fCPhiResolS1PtTPCITS->Add(entry->fCPhiResolS1PtTPCITS);
950 fCThetaResolS1PtTPCITS->Add(entry->fCThetaResolS1PtTPCITS);
96b3109a 951
952 count++;
953 }
954
955return count;
956}
3baa4bfd 957
958//_____________________________________________________________________________
959TFolder* AliComparisonRes::CreateFolder(TString name,TString title) {
960// create folder for analysed histograms
961//
962TFolder *folder = 0;
963 folder = new TFolder(name.Data(),title.Data());
964
965 return folder;
966}