]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoAvgSepCorrFctn.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoAvgSepCorrFctn.cxx
CommitLineData
07816065 1///////////////////////////////////////////////////////////////////////////
2// //
3// AliFemtoAvgSepCorrFctn: //
4// a simple Q-invariant correlation function //
5// //
6///////////////////////////////////////////////////////////////////////////
7
8#include "AliFemtoAvgSepCorrFctn.h"
9//#include "AliFemtoHisto.h"
10#include <cstdio>
11
12#ifdef __ROOT__
13ClassImp(AliFemtoAvgSepCorrFctn)
14#endif
15
16//____________________________
17AliFemtoAvgSepCorrFctn::AliFemtoAvgSepCorrFctn(char* title, const int& nbins, const float& Low, const float& High):
18fNumerator(0), //2 tracks
19 fDenominator(0),
20 fNumeratorPos(0), //track + V0
21 fDenominatorPos(0),
22 fNumeratorNeg(0),
23 fDenominatorNeg(0),
24 fNumeratorPosPos(0), //2 V0s
25 fDenominatorPosPos(0),
26 fNumeratorPosNeg(0),
27 fDenominatorPosNeg(0),
28 fNumeratorNegPos(0),
29 fDenominatorNegPos(0),
30 fNumeratorNegNeg(0),
31 fDenominatorNegNeg(0),
32 fRatio(0),
33 fPairType(kTracks)
34{
35 // set up numerator
36 // title = "Num Qinv (MeV/c)";
37 char tTitNum[101] = "Num";
38 strncat(tTitNum,title, 100);
39 fNumerator = new TH1D(tTitNum,title,nbins,Low,High);
40 // set up denominator
41 //title = "Den Qinv (MeV/c)";
42 char tTitDen[101] = "Den";
43 strncat(tTitDen,title, 100);
44 fDenominator = new TH1D(tTitDen,title,nbins,Low,High);
45 // set up ratio
46 //title = "Ratio Qinv (MeV/c)";
47 char tTitRat[101] = "Rat";
48 strncat(tTitRat,title, 100);
49 fRatio = new TH1D(tTitRat,title,nbins,Low,High);
50
51
52 char tTitNumPos[101] = "NumV0TrackPos";
53 strncat(tTitNumPos,title, 100);
54 fNumeratorPos = new TH1D(tTitNumPos,title,nbins,Low,High);
55 char tTitDenPos[101] = "DenV0TrackPos";
56 strncat(tTitDenPos,title, 100);
57 fDenominatorPos = new TH1D(tTitDenPos,title,nbins,Low,High);
58 char tTitNumNeg[101] = "NumV0TrackNeg";
59 strncat(tTitNumNeg,title, 100);
60 fNumeratorNeg = new TH1D(tTitNumNeg,title,nbins,Low,High);
61 char tTitDenNeg[101] = "DenV0TrackNeg";
62 strncat(tTitDenNeg,title, 100);
63 fDenominatorNeg = new TH1D(tTitDenNeg,title,nbins,Low,High);
64
65
66 char tTitNumPosPos[101] = "NumV0sPosPos";
67 strncat(tTitNumPosPos,title, 100);
68 fNumeratorPosPos = new TH1D(tTitNumPosPos,title,nbins,Low,High);
69 char tTitDenPosPos[101] = "DenV0sPosPos";
70 strncat(tTitDenPosPos,title, 100);
71 fDenominatorPosPos = new TH1D(tTitDenPosPos,title,nbins,Low,High);
72 char tTitNumPosNeg[101] = "NumV0sPosNeg";
73 strncat(tTitNumPosNeg,title, 100);
74 fNumeratorPosNeg = new TH1D(tTitNumPosNeg,title,nbins,Low,High);
75 char tTitDenPosNeg[101] = "DenV0sPosNeg";
76 strncat(tTitDenPosNeg,title, 100);
77 fDenominatorPosNeg = new TH1D(tTitDenPosNeg,title,nbins,Low,High);
78 char tTitNumNegPos[101] = "NumV0sNegPos";
79 strncat(tTitNumNegPos,title, 100);
80 fNumeratorNegPos = new TH1D(tTitNumNegPos,title,nbins,Low,High);
81 char tTitDenNegPos[101] = "DenV0sNegPos";
82 strncat(tTitDenNegPos,title, 100);
83 fDenominatorNegPos = new TH1D(tTitDenNegPos,title,nbins,Low,High);
84 char tTitNumNegNeg[101] = "NumV0sNegNeg";
85 strncat(tTitNumNegNeg,title, 100);
86 fNumeratorNegNeg = new TH1D(tTitNumNegNeg,title,nbins,Low,High);
87 char tTitDenNegNeg[101] = "DenV0sNegNeg";
88 strncat(tTitDenNegNeg,title, 100);
89 fDenominatorNegNeg = new TH1D(tTitDenNegNeg,title,nbins,Low,High);
90
91
92 // to enable error bar calculation...
93 fNumerator->Sumw2();
94 fDenominator->Sumw2();
95
96 fNumeratorPos->Sumw2();
97 fDenominatorPos->Sumw2();
98 fNumeratorNeg->Sumw2();
99 fDenominatorNeg->Sumw2();
100
101 fNumeratorPosPos->Sumw2();
102 fDenominatorPosPos->Sumw2();
103 fNumeratorPosNeg->Sumw2();
104 fDenominatorPosNeg->Sumw2();
105 fNumeratorNegPos->Sumw2();
106 fDenominatorNegPos->Sumw2();
107 fNumeratorNegNeg->Sumw2();
108 fDenominatorNegNeg->Sumw2();
109
110 fRatio->Sumw2();
111}
112
113//____________________________
114AliFemtoAvgSepCorrFctn::AliFemtoAvgSepCorrFctn(const AliFemtoAvgSepCorrFctn& aCorrFctn) :
115 AliFemtoCorrFctn(),
116 fNumerator(0),
117 fDenominator(0),
118 fNumeratorPos(0),
119 fDenominatorPos(0),
120 fNumeratorNeg(0),
121 fDenominatorNeg(0),
122 fNumeratorPosPos(0),
123 fDenominatorPosPos(0),
124 fNumeratorPosNeg(0),
125 fDenominatorPosNeg(0),
126 fNumeratorNegPos(0),
127 fDenominatorNegPos(0),
128 fNumeratorNegNeg(0),
129 fDenominatorNegNeg(0),
130 fRatio(0),
131 fPairType(kTracks)
132{
133 // copy constructor
134 fNumerator = new TH1D(*aCorrFctn.fNumerator);
135 fDenominator = new TH1D(*aCorrFctn.fDenominator);
136
137 fNumeratorPos = new TH1D(*aCorrFctn.fNumerator);
138 fDenominatorPos = new TH1D(*aCorrFctn.fDenominator);
139 fNumeratorNeg = new TH1D(*aCorrFctn.fNumerator);
140 fDenominatorNeg = new TH1D(*aCorrFctn.fDenominator);
141
142 fNumeratorPosPos = new TH1D(*aCorrFctn.fNumerator);
143 fDenominatorPosPos = new TH1D(*aCorrFctn.fDenominator);
144 fNumeratorPosNeg = new TH1D(*aCorrFctn.fNumerator);
145 fDenominatorPosNeg = new TH1D(*aCorrFctn.fDenominator);
146 fNumeratorNegPos = new TH1D(*aCorrFctn.fNumerator);
147 fDenominatorNegPos = new TH1D(*aCorrFctn.fDenominator);
148 fNumeratorNegNeg = new TH1D(*aCorrFctn.fNumerator);
149 fDenominatorNegNeg = new TH1D(*aCorrFctn.fDenominator);
150
151 fRatio = new TH1D(*aCorrFctn.fRatio);
152}
153//____________________________
154AliFemtoAvgSepCorrFctn::~AliFemtoAvgSepCorrFctn(){
155 // destructor
156 delete fNumerator;
157 delete fDenominator;
158
159 delete fNumeratorPos;
160 delete fDenominatorPos;
161 delete fNumeratorNeg;
162 delete fDenominatorNeg;
163
164 delete fNumeratorPosPos;
165 delete fDenominatorPosPos;
166 delete fNumeratorPosNeg;
167 delete fDenominatorPosNeg;
168 delete fNumeratorNegPos;
169 delete fDenominatorNegPos;
170 delete fNumeratorNegNeg;
171 delete fDenominatorNegNeg;
172
173 delete fRatio;
174}
175//_________________________
176AliFemtoAvgSepCorrFctn& AliFemtoAvgSepCorrFctn::operator=(const AliFemtoAvgSepCorrFctn& aCorrFctn)
177{
178 // assignment operator
179 if (this == &aCorrFctn)
180 return *this;
181
182 if (fNumerator) delete fNumerator;
183 fNumerator = new TH1D(*aCorrFctn.fNumerator);
184 if (fDenominator) delete fDenominator;
185 fDenominator = new TH1D(*aCorrFctn.fDenominator);
186
187 if (fNumeratorPos) delete fNumeratorPos;
188 fNumeratorPos = new TH1D(*aCorrFctn.fNumeratorPos);
189 if (fDenominatorPos) delete fDenominatorPos;
190 fDenominatorPos = new TH1D(*aCorrFctn.fDenominatorPos);
191 if (fNumeratorNeg) delete fNumeratorNeg;
192 fNumeratorNeg = new TH1D(*aCorrFctn.fNumeratorNeg);
193 if (fDenominatorNeg) delete fDenominatorNeg;
194 fDenominatorNeg = new TH1D(*aCorrFctn.fDenominatorNeg);
195
196 if (fNumeratorPosPos) delete fNumeratorPosPos;
197 fNumeratorPosPos = new TH1D(*aCorrFctn.fNumeratorPosPos);
198 if (fDenominatorPosPos) delete fDenominatorPosPos;
199 fDenominatorPosPos = new TH1D(*aCorrFctn.fDenominatorPosPos);
200 if (fNumeratorPosNeg) delete fNumeratorPosNeg;
201 fNumeratorPosNeg = new TH1D(*aCorrFctn.fNumeratorPosNeg);
202 if (fDenominatorPosNeg) delete fDenominatorPosNeg;
203 fDenominatorPosNeg = new TH1D(*aCorrFctn.fDenominatorPosNeg);
204 if (fNumeratorNegPos) delete fNumeratorNegPos;
205 fNumeratorNegPos = new TH1D(*aCorrFctn.fNumeratorNegPos);
206 if (fDenominatorNegPos) delete fDenominatorNegPos;
207 fDenominatorNegPos = new TH1D(*aCorrFctn.fDenominatorNegPos);
208 if (fNumeratorNegNeg) delete fNumeratorNegNeg;
209 fNumeratorNegNeg = new TH1D(*aCorrFctn.fNumeratorNegNeg);
210 if (fDenominatorNegNeg) delete fDenominatorNegNeg;
211 fDenominatorNegNeg = new TH1D(*aCorrFctn.fDenominatorNegNeg);
212
213 if (fRatio) delete fRatio;
214 fRatio = new TH1D(*aCorrFctn.fRatio);
215
216 return *this;
217}
218
219//_________________________
220void AliFemtoAvgSepCorrFctn::Finish(){
221 // here is where we should normalize, fit, etc...
222 // we should NOT Draw() the histos (as I had done it below),
223 // since we want to insulate ourselves from root at this level
224 // of the code. Do it instead at root command line with browser.
225 // fNumerator->Draw();
226 //fDenominator->Draw();
227 //fRatio->Draw();
228 fRatio->Divide(fNumerator,fDenominator,1.0,1.0);
229
230}
231
232//____________________________
233AliFemtoString AliFemtoAvgSepCorrFctn::Report(){
234 // construct report
235 string stemp = "Qinv Correlation Function Report:\n";
236 char ctemp[100];
237 snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fNumerator->GetEntries());
238 stemp += ctemp;
239 snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fDenominator->GetEntries());
240 stemp += ctemp;
241 snprintf(ctemp , 100, "Number of entries in ratio:\t%E\n",fRatio->GetEntries());
242 stemp += ctemp;
243 // stemp += mCoulombWeight->Report();
244 AliFemtoString returnThis = stemp;
245 return returnThis;
246}
247//____________________________
248void AliFemtoAvgSepCorrFctn::AddRealPair(AliFemtoPair* pair){
249 // add true pair
250 if (fPairCut)
251 if (!fPairCut->Pass(pair)) return;
252
253 double avgSep=0;
254 AliFemtoThreeVector first, second, tmp;
255
256
257 if(fPairType == 0) //2 tracks
258 {
259 for(int i=0; i<8 ;i++)
260 {
261 tmp = pair->Track1()->Track()->NominalTpcPoint(i);
262 first.SetX((double)(tmp.x()));
263 first.SetY((double)tmp.y());
264 first.SetZ((double)tmp.z());
265
266 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
267 second.SetX((double)tmp.x());
268 second.SetY((double)tmp.y());
269 second.SetZ((double)tmp.z());
270
271 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
272 }
273 avgSep /= 8;
274 fNumerator->Fill(avgSep);
275 }
276 else if(fPairType == 1) // track + V0
277 {
278 for(int i=0; i<8 ;i++)
279 {
280 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
281 first.SetX((double)(tmp.x()));
282 first.SetY((double)tmp.y());
283 first.SetZ((double)tmp.z());
284 //cout<<"V0pos x: "<<tmp.x()<<", y: "<<tmp.y()<<", z: "<<tmp.z()<<endl;
285 //cout<<"V0pos x: "<<first.x()<<", y: "<<first.y()<<", z: "<<first.z()<<endl;
286
287 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
288 second.SetX((double)tmp.x());
289 second.SetY((double)tmp.y());
290 second.SetZ((double)tmp.z());
291 //cout<<"Track x: "<<tmp.x()<<", y: "<<tmp.y()<<", z: "<<tmp.z()<<endl;
292 //cout<<"Track x: "<<second.x()<<", y: "<<second.y()<<", z: "<<second.z()<<endl;
293
294 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
295 }
296 //cout<<"****************************************************"<<endl;
297 avgSep /= 8;
298 fNumeratorPos->Fill(avgSep);
299 //cout<<"Track + Pos V0 Avg Sep: "<<avgSep<<endl;
300 avgSep = 0;
301
302 for(int i=0; i<8 ;i++)
303 {
304 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
305 first.SetX((double)(tmp.x()));
306 first.SetY((double)tmp.y());
307 first.SetZ((double)tmp.z());
308 //cout<<"V0 X: "<<tmp.x()<<endl;
309
310 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
311 second.SetX((double)tmp.x());
312 second.SetY((double)tmp.y());
313 second.SetZ((double)tmp.z());
314 //cout<<"Track X: "<<tmp.x()<<endl;
315
316 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
317 }
318 fNumeratorNeg->Fill(avgSep);
319
320 }
321 else if(fPairType == 2)
322 {
323 for(int i=0; i<8 ;i++)
324 {
325 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
326 //cout<<"X pos: "<<tmp.x()<<endl;
327 first.SetX((double)(tmp.x()));
328 first.SetY((double)tmp.y());
329 first.SetZ((double)tmp.z());
330
331 tmp = pair->Track2()->V0()->NominalTpcPointPos(i);
332 second.SetX((double)tmp.x());
333 second.SetY((double)tmp.y());
334 second.SetZ((double)tmp.z());
335
336 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
337 }
338 avgSep /= 8;
339 fNumeratorPosPos->Fill(avgSep);
340 //cout<<"PovV0 + PosV0 Avg Sep: "<<avgSep<<endl;
341
342 avgSep = 0;
343
344 for(int i=0; i<8 ;i++)
345 {
346 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
347 first.SetX((double)(tmp.x()));
348 first.SetY((double)tmp.y());
349 first.SetZ((double)tmp.z());
350
351 tmp = pair->Track2()->V0()->NominalTpcPointNeg(i);
352 //cout<<"X neg: "<<tmp.x()<<endl;
353 second.SetX((double)tmp.x());
354 second.SetY((double)tmp.y());
355 second.SetZ((double)tmp.z());
356
357 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
358 }
359 avgSep /= 8;
360 fNumeratorPosNeg->Fill(avgSep);
361
362 avgSep = 0;
363
364 for(int i=0; i<8 ;i++)
365 {
366 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
367 first.SetX((double)(tmp.x()));
368 first.SetY((double)tmp.y());
369 first.SetZ((double)tmp.z());
370
371 tmp = pair->Track2()->V0()->NominalTpcPointPos(i);
372 second.SetX((double)tmp.x());
373 second.SetY((double)tmp.y());
374 second.SetZ((double)tmp.z());
375
376 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
377 }
378 avgSep /= 8;
379 fNumeratorNegPos->Fill(avgSep);
380
381
382 avgSep = 0;
383
384 for(int i=0; i<8 ;i++)
385 {
386 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
387 first.SetX((double)(tmp.x()));
388 first.SetY((double)tmp.y());
389 first.SetZ((double)tmp.z());
390
391 tmp = pair->Track2()->V0()->NominalTpcPointNeg(i);
392 second.SetX((double)tmp.x());
393 second.SetY((double)tmp.y());
394 second.SetZ((double)tmp.z());
395
396 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
397 }
398 avgSep /= 8;
399 fNumeratorNegNeg->Fill(avgSep);
400
401
402 }
403
404
405 //2 daughters of first V0
406
407
408
409
410
411}
412//____________________________
413void AliFemtoAvgSepCorrFctn::AddMixedPair(AliFemtoPair* pair){
414 // add mixed (background) pair
415 if (fPairCut)
416 if (!fPairCut->Pass(pair)) return;
417
418 double avgSep=0;
419 AliFemtoThreeVector first, second, tmp;
420
421
422 if(fPairType == 0) //2 tracks
423 {
424 for(int i=0; i<8 ;i++)
425 {
426 tmp = pair->Track1()->Track()->NominalTpcPoint(i);
427 first.SetX((double)(tmp.x()));
428 first.SetY((double)tmp.y());
429 first.SetZ((double)tmp.z());
430
431 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
432 second.SetX((double)tmp.x());
433 second.SetY((double)tmp.y());
434 second.SetZ((double)tmp.z());
435
436 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
437 }
438 avgSep /= 8;
439 fDenominator->Fill(avgSep);
440 }
441 else if(fPairType == 1) // track + V0
442 {
443 for(int i=0; i<8 ;i++)
444 {
445 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
446 first.SetX((double)(tmp.x()));
447 first.SetY((double)tmp.y());
448 first.SetZ((double)tmp.z());
449
450 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
451 second.SetX((double)tmp.x());
452 second.SetY((double)tmp.y());
453 second.SetZ((double)tmp.z());
454
455 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
456 }
457 avgSep /= 8;
458 fDenominatorPos->Fill(avgSep);
459
460 avgSep = 0;
461
462 for(int i=0; i<8 ;i++)
463 {
464 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
465 first.SetX((double)(tmp.x()));
466 first.SetY((double)tmp.y());
467 first.SetZ((double)tmp.z());
468
469 tmp = pair->Track2()->Track()->NominalTpcPoint(i);
470 second.SetX((double)tmp.x());
471 second.SetY((double)tmp.y());
472 second.SetZ((double)tmp.z());
473
474 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
475 }
476 fDenominatorNeg->Fill(avgSep);
477
478 }
479 else if(fPairType == 2)
480 {
481 for(int i=0; i<8 ;i++)
482 {
483
484 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
485 first.SetX((double)(tmp.x()));
486 first.SetY((double)tmp.y());
487 first.SetZ((double)tmp.z());
488
489 tmp = pair->Track2()->V0()->NominalTpcPointPos(i);
490 second.SetX((double)tmp.x());
491 second.SetY((double)tmp.y());
492 second.SetZ((double)tmp.z());
493
494 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
495 }
496 avgSep /= 8;
497 fDenominatorPosPos->Fill(avgSep);
498
499 avgSep = 0;
500
501 for(int i=0; i<8 ;i++)
502 {
503 tmp = pair->Track1()->V0()->NominalTpcPointPos(i);
504 first.SetX((double)(tmp.x()));
505 first.SetY((double)tmp.y());
506 first.SetZ((double)tmp.z());
507
508 tmp = pair->Track2()->V0()->NominalTpcPointNeg(i);
509 second.SetX((double)tmp.x());
510 second.SetY((double)tmp.y());
511 second.SetZ((double)tmp.z());
512
513 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
514 }
515 avgSep /= 8;
516 fDenominatorPosNeg->Fill(avgSep);
517
518 avgSep = 0;
519
520 for(int i=0; i<8 ;i++)
521 {
522 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
523 first.SetX((double)(tmp.x()));
524 first.SetY((double)tmp.y());
525 first.SetZ((double)tmp.z());
526
527 tmp = pair->Track2()->V0()->NominalTpcPointPos(i);
528 second.SetX((double)tmp.x());
529 second.SetY((double)tmp.y());
530 second.SetZ((double)tmp.z());
531
532 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
533 }
534 avgSep /= 8;
535 fDenominatorNegPos->Fill(avgSep);
536
537 avgSep = 0;
538
539 for(int i=0; i<8 ;i++)
540 {
541 tmp = pair->Track1()->V0()->NominalTpcPointNeg(i);
542 first.SetX((double)(tmp.x()));
543 first.SetY((double)tmp.y());
544 first.SetZ((double)tmp.z());
545
546 tmp = pair->Track2()->V0()->NominalTpcPointNeg(i);
547 second.SetX((double)tmp.x());
548 second.SetY((double)tmp.y());
549 second.SetZ((double)tmp.z());
550
551 avgSep += TMath::Sqrt(((double)first.x()-(double)second.x())*((double)first.x()-(double)second.x())+((double)first.y()-(double)second.y())*((double)first.y()-second.y())+((double)first.z()-(double)second.z())*((double)first.z()-(double)second.z()));
552 }
553 avgSep /= 8;
554 fDenominatorNegNeg->Fill(avgSep);
555
556
557 }
558
559}
560//____________________________
561void AliFemtoAvgSepCorrFctn::Write(){
562 // Write out neccessary objects
563 if(fPairType == kTracks)
564 {
565 fNumerator->Write();
566 fDenominator->Write();
567 fRatio->Write();
568 }
569 else if(fPairType == kTrackV0)
570 {
571 fNumeratorPos->Write();
572 fDenominatorPos->Write();
573 fNumeratorNeg->Write();
574 fDenominatorNeg->Write();
575 }
576 else if(fPairType == kV0s)
577 {
578 fNumeratorPosPos->Write();
579 fDenominatorPosPos->Write();
580 fNumeratorPosNeg->Write();
581 fDenominatorPosNeg->Write();
582 fNumeratorNegPos->Write();
583 fDenominatorNegPos->Write();
584 fNumeratorNegNeg->Write();
585 fDenominatorNegNeg->Write();
586 }
587
588
589}
590//______________________________
591TList* AliFemtoAvgSepCorrFctn::GetOutputList()
592{
593 // Prepare the list of objects to be written to the output
594 TList *tOutputList = new TList();
595
596 if(fPairType == kTracks)
597 {
598 tOutputList->Add(fNumerator);
599 tOutputList->Add(fDenominator);
600 tOutputList->Add(fRatio);
601 }
602 else if(fPairType == kTrackV0)
603 {
604 tOutputList->Add(fNumeratorPos);
605 tOutputList->Add(fDenominatorPos);
606 tOutputList->Add(fNumeratorNeg);
607 tOutputList->Add(fDenominatorNeg);
608 }
609 else if(fPairType == kV0s)
610 {
611 tOutputList->Add(fNumeratorPosPos);
612 tOutputList->Add(fDenominatorPosPos);
613 tOutputList->Add(fNumeratorPosNeg);
614 tOutputList->Add(fDenominatorPosNeg);
615 tOutputList->Add(fNumeratorNegPos);
616 tOutputList->Add(fDenominatorNegPos);
617 tOutputList->Add(fNumeratorNegNeg);
618 tOutputList->Add(fDenominatorNegNeg);
619 }
620 return tOutputList;
621}
622
623
624void AliFemtoAvgSepCorrFctn::SetPairType(AliFemtoPairType pairtype)
625{
626 fPairType = pairtype;
627}