new functions added
[u/mrichter/AliRoot.git] / HBTAN / AliHBTPIDPurityFctns.cxx
1 #include "AliHBTPIDPurityFctns.h"
2 //_______________________________________________________________________________
3 /////////////////////////////////////////////////////////////////////////////////
4 //
5 // class AliHBTMonPhiResolutionVsPtFctn;
6 // class AliHBTMonThetaResolutionVsPtFctn;
7 //
8 // file: AliHBTPIDPurityFctns.cxx AliHBTPIDPurityFctns.h
9 //
10 // Caution: On 2D plots on X axis in simulated values
11 // That is contrary to two-particle resolutions where it is reconstructed one
12 //
13 // added by Piotr.Skowronski@cern.ch
14 // 
15 //
16 //////////////////////////////////////////////////////////////////////////////////
17
18
19 /******************************************************************/
20 /******************************************************************/
21 /******************************************************************/
22 AliHBTQInvCorrelFctnPerfectPID f;
23
24 ClassImp(AliHBTMonPIDPurityVsPtFctn)
25
26 AliHBTMonPIDPurityVsPtFctn::AliHBTMonPIDPurityVsPtFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
27  AliHBTMonTwoParticleFctn1D(nbins,maxXval,minXval),
28  fGood(0x0),
29  fAll(0x0)
30 {
31   //ctor
32   Rename("pidpurityvspt","PIDPurityVsPt");
33   
34 }
35 /******************************************************************/
36
37 AliHBTMonPIDPurityVsPtFctn::~AliHBTMonPIDPurityVsPtFctn()
38 {
39  //dtor
40   delete fGood;
41   delete fAll;
42 }
43 /******************************************************************/
44 void AliHBTMonPIDPurityVsPtFctn::Write()
45 {
46  AliHBTMonitorFunction::Write();
47  fGood->Write();
48  fAll->Write();
49 }
50 /******************************************************************/
51
52 void AliHBTMonPIDPurityVsPtFctn::Init()
53 {
54 //Initializes fuction
55   if (AliHBTParticle::GetDebug()>0) Info("Init","%s",GetName());
56
57   if (fResult == 0x0)
58    {   
59       Warning("Init","Function has NULL result histogram!");
60       return;
61    }
62   
63   if (fGood == 0x0)
64    {
65      TString numstr = fName + " Good";  //title and name of the
66                                            //result histogram
67      TAxis* xax = fResult->GetXaxis();
68      fGood = new TH1D(numstr,numstr,xax->GetNbins(),xax->GetXmin(),xax->GetXmax());
69    }
70
71   if (fAll == 0x0)
72    {
73      TString numstr = fName + " All";  //title and name of the
74                                            //result histogram
75      TAxis* xax = fResult->GetXaxis();
76      fAll = new TH1D(numstr,numstr,xax->GetNbins(),xax->GetXmin(),xax->GetXmax());
77    }
78
79   fResult->Reset();
80   fResult->SetDirectory(0x0);
81   fResult->Sumw2();
82   fGood->Reset();
83   fGood->SetDirectory(0x0);
84   fGood->Sumw2();
85   fAll->Reset();
86   fAll->SetDirectory(0x0);
87   fAll->Sumw2();
88
89   if (AliHBTParticle::GetDebug()>0) Info("Init","%s Done.",GetName());
90 }
91
92 /******************************************************************/
93
94 void AliHBTMonPIDPurityVsPtFctn::Rename(const Char_t * name)
95
96   //Rename fuctions and all histograms belonging to it
97   SetName(name);
98   SetTitle(name);
99   
100   if (fResult)
101    {
102      TString numstr = fName + " Result";  //title and name of the result histogram
103      fResult->SetName(numstr);
104      fResult->SetTitle(numstr);
105    }
106    
107   if (fGood)
108    {
109      TString numstr = fName + " Good";
110      fGood->SetName(numstr);
111      fGood->SetTitle(numstr);
112    }
113    
114   if (fAll)
115    {
116      TString numstr = fName + " All";
117      fAll->SetName(numstr);
118      fAll->SetTitle(numstr);
119    }
120 }
121 /******************************************************************/
122
123 void AliHBTMonPIDPurityVsPtFctn::Rename(const Char_t * name, const Char_t * title)
124 {
125  //renames and retitle the function and histograms
126  
127   SetName(name);
128   SetTitle(title);
129   
130   if (fResult)
131    {
132      TString numstrn = fName + " Result";  //name of the result histogram
133      TString numstrt = fTitle + " Result";  //title of the result histogram
134      fResult->SetName(numstrn);
135      fResult->SetTitle(numstrt);
136    }
137   if (fGood)
138    {
139      TString numstrn = fName + " Good";  //name of the Good histogram
140      TString numstrt = fTitle + " Good";  //title of the Good histogram
141      fGood->SetName(numstrn);
142      fGood->SetTitle(numstrt);
143    }
144
145   if (fAll)
146    {
147      TString numstrn = fName + " All";  //name of the All histogram
148      TString numstrt = fTitle + " All";  //title of the All histogram
149      fAll->SetName(numstrn);
150      fAll->SetTitle(numstrt);
151    }
152
153 }
154 /******************************************************************/
155
156 TH1* AliHBTMonPIDPurityVsPtFctn::GetResult()
157 {
158   //Returns the result of the fuction
159   //that is histogram with effciency and contamination
160   
161   fResult->Divide(fGood,fAll);
162   return fResult;
163 }
164 /******************************************************************/
165
166 void AliHBTMonPIDPurityVsPtFctn::Process(AliHBTParticle * track,AliHBTParticle * part)
167 {
168  //process the particle/track
169  Double_t pt = part->Pt();
170  fAll->Fill(pt);
171  if (track->GetPdgCode() == part->GetPdgCode()) 
172   {
173     fGood->Fill(pt);
174   }
175 // else
176 //  {
177 //    Info("Process","Catched pid impurity ...");
178 //  }
179 }
180
181 /******************************************************************/
182 /******************************************************************/
183 /******************************************************************/
184
185 ClassImp(AliHBTMonPIDContaminationVsPtFctn)
186
187 AliHBTMonPIDContaminationVsPtFctn::AliHBTMonPIDContaminationVsPtFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
188  AliHBTMonTwoParticleFctn1D(nbins,maxXval,minXval),
189  fWrong(0x0),
190  fAll(0x0)
191 {
192   //ctor
193   Rename("pidcontaminationvspt","PIDContaminationVsPt");
194 }
195 /******************************************************************/
196
197 AliHBTMonPIDContaminationVsPtFctn::~AliHBTMonPIDContaminationVsPtFctn()
198 {
199  //dtor
200   delete fWrong;
201   delete fAll;
202 }
203 /******************************************************************/
204
205 void AliHBTMonPIDContaminationVsPtFctn::Write()
206 {
207  //Writes the function results
208  AliHBTMonitorFunction::Write();
209  fWrong->Write();
210  fAll->Write();
211 }
212 /******************************************************************/
213
214 void AliHBTMonPIDContaminationVsPtFctn::Init()
215 {
216 //Initializes fuction
217   if (AliHBTParticle::GetDebug()>0) Info("Init","%s",GetName());
218
219   if (fResult == 0x0)
220    {   
221       Warning("Init","Function has NULL result histogram!");
222       return;
223    }
224   
225   if (fWrong == 0x0)
226    {
227      TString numstr = fName + " Wrong";  //title and name of the
228                                            //result histogram
229      TAxis* xax = fResult->GetXaxis();
230      fWrong = new TH1D(numstr,numstr,xax->GetNbins(),xax->GetXmin(),xax->GetXmax());
231    }
232
233   if (fAll == 0x0)
234    {
235      TString numstr = fName + " All";  //title and name of the
236                                            //result histogram
237      TAxis* xax = fResult->GetXaxis();
238      fAll = new TH1D(numstr,numstr,xax->GetNbins(),xax->GetXmin(),xax->GetXmax());
239    }
240   fResult->Reset();
241   fResult->SetDirectory(0x0);
242   fResult->Sumw2();
243   fWrong->Reset();
244   fWrong->SetDirectory(0x0);
245   fWrong->Sumw2();
246   fAll->Reset();
247   fAll->SetDirectory(0x0);
248   fAll->Sumw2();
249   
250   if (AliHBTParticle::GetDebug()>0) Info("Init","%s Done.",GetName());
251 }
252
253 /******************************************************************/
254
255 void AliHBTMonPIDContaminationVsPtFctn::Rename(const Char_t * name)
256
257   //Rename fuctions and all histograms belonging to it
258   SetName(name);
259   SetTitle(name);
260   
261   if (fResult)
262    {
263      TString numstr = fName + " Result";  //title and name of the result histogram
264      fResult->SetName(numstr);
265      fResult->SetTitle(numstr);
266    }
267   if (fWrong)
268    {
269      TString numstr = fName + " Wrong";
270      fWrong->SetName(numstr);
271      fWrong->SetTitle(numstr);
272    }
273
274   if (fAll)
275    {
276      TString numstrn = fName + " All";  //name of the All histogram
277      TString numstrt = fTitle + " All";  //title of the All histogram
278      fAll->SetName(numstrn);
279      fAll->SetTitle(numstrt);
280    }
281 }
282 /******************************************************************/
283
284 void AliHBTMonPIDContaminationVsPtFctn::Rename(const Char_t * name, const Char_t * title)
285 {
286  //renames and retitle the function and histograms
287  
288   SetName(name);
289   SetTitle(title);
290   
291   if (fResult)
292    {
293      TString numstrn = fName + " Result";  //name of the result histogram
294      TString numstrt = fTitle + " Result";  //title of the result histogram
295      fResult->SetName(numstrn);
296      fResult->SetTitle(numstrt);
297    }
298   if (fWrong)
299    {
300      TString numstrn = fName + " Wrong";  //name of the Wrong histogram
301      TString numstrt = fTitle + " Wrong";  //title of the Wrong histogram
302      fWrong->SetName(numstrn);
303      fWrong->SetTitle(numstrt);
304    }
305    
306   if (fAll)
307    {
308      TString numstr = fName + " All";
309      fAll->SetName(numstr);
310      fAll->SetTitle(numstr);
311    }
312
313 }
314 /******************************************************************/
315
316 TH1* AliHBTMonPIDContaminationVsPtFctn::GetResult()
317 {
318   //Returns the result of the fuction
319   //that is histogram with effciency and contamination
320   
321   fResult->Divide(fWrong,fAll);
322   return fResult;
323 }
324 /******************************************************************/
325
326 void AliHBTMonPIDContaminationVsPtFctn::Process(AliHBTParticle * track, AliHBTParticle * part)
327 {
328  //process the particle/track
329  Double_t pt = part->Pt();
330  fAll->Fill(pt);
331  
332  if (track->GetPdgCode() != part->GetPdgCode()) 
333   {
334 //    Info("Process","Catched contamination");
335 //    track->Print();part->Print();
336     fWrong->Fill(pt);
337   }
338 }
339
340 /******************************************************************/
341 /******************************************************************/
342 /******************************************************************/
343 ClassImp(AliHBTQInvCorrelFctnPerfectPID)
344
345 AliHBTQInvCorrelFctnPerfectPID::AliHBTQInvCorrelFctnPerfectPID(Int_t nbins, Double_t maxXval, Double_t minXval):
346  AliHBTTwoPairFctn1D(nbins,maxXval,minXval)
347 {
348 //ctor
349  fWriteNumAndDen = kTRUE;//change default behaviour
350  Rename("qinvcfprfctpid","Q_{inv} Correlation Function Perfect PID");
351 }
352 /*************************************************************/
353
354 void AliHBTQInvCorrelFctnPerfectPID::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
355 {
356 //Fills numerator
357   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
358   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
359
360   trackpair  = CheckPair(trackpair);
361   if (trackpair == 0x0) return;
362   partpair  = CheckPair(partpair);
363   if (partpair == 0x0) return;
364
365   fNumerator->Fill(trackpair->GetQInv());
366 }
367
368 /*************************************************************/
369 void AliHBTQInvCorrelFctnPerfectPID::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
370 {
371 //Fills numerator
372   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
373   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
374
375   trackpair  = CheckPair(trackpair);
376   if (trackpair == 0x0) return;
377   partpair  = CheckPair(partpair);
378   if (partpair == 0x0) return;
379
380   fDenominator->Fill(trackpair->GetQInv());
381 }
382 /*************************************************************/
383
384 TH1* AliHBTQInvCorrelFctnPerfectPID::GetResult()
385 {
386  //returns the scaled ratio
387  
388  delete fRatio;
389  fRatio = GetRatio(Scale());
390  return fRatio;
391 }
392 /*************************************************************/
393 /******************************************************************/
394 /******************************************************************/
395 /******************************************************************/
396
397 ClassImp(AliHBTWeightQInvCorrelFctnPerfectPID)
398
399 AliHBTWeightQInvCorrelFctnPerfectPID::AliHBTWeightQInvCorrelFctnPerfectPID(Int_t nbins, Double_t maxXval, Double_t minXval):
400  AliHBTTwoPairFctn1D(nbins,maxXval,minXval)
401 {
402 //ctor
403  fWriteNumAndDen = kTRUE;//change default behaviour
404  Rename("wqinvcfprfctpid","Q_{inv} Weight Correlation Function Perfect PID");
405 }
406 /*************************************************************/
407
408 void AliHBTWeightQInvCorrelFctnPerfectPID::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
409 {
410 //Fills numerator
411   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
412   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
413
414   trackpair  = CheckPair(trackpair);
415   if (trackpair == 0x0) return;
416   partpair  = CheckPair(partpair);
417   if (partpair == 0x0) return;
418
419   fNumerator->Fill(trackpair->GetQInv(),partpair->GetWeight());
420 }
421
422 /*************************************************************/
423 void AliHBTWeightQInvCorrelFctnPerfectPID::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
424 {
425 //Fills numerator
426   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
427   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
428
429   trackpair  = CheckPair(trackpair);
430   if (trackpair == 0x0) return;
431   partpair  = CheckPair(partpair);
432   if (partpair == 0x0) return;
433   
434   fDenominator->Fill(trackpair->GetQInv());
435 }
436 /*************************************************************/
437
438 TH1* AliHBTWeightQInvCorrelFctnPerfectPID::GetResult()
439 {
440  //returns the scaled ratio
441  
442  delete fRatio;
443  fRatio = GetRatio(Scale());
444  return fRatio;
445 }
446
447
448 /******************************************************************/
449 /******************************************************************/
450 /******************************************************************/
451
452 ClassImp(AliHBTWeightQOutSQideQLongFctnPerfectPID)
453
454 AliHBTWeightQOutSQideQLongFctnPerfectPID::AliHBTWeightQOutSQideQLongFctnPerfectPID
455   (Int_t nXbins, Double_t maxXval, Double_t minXval,
456    Int_t nYbins, Double_t maxYval, Double_t minYval,
457    Int_t nZbins, Double_t maxZval, Double_t minZval):
458  AliHBTTwoPairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
459 {
460   //ctor
461   fWriteNumAndDen = kTRUE;//change default behaviour
462   Rename("wqoslprfctpid","Q_{out}-Q_{side}-Q_{long} Weight Fctn with perfect PID");
463 }
464 /*************************************************************/
465
466 void AliHBTWeightQOutSQideQLongFctnPerfectPID::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
467 {
468 //Fills numerator
469   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
470   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
471
472   trackpair  = CheckPair(trackpair);
473   if (trackpair == 0x0) return;
474   partpair  = CheckPair(partpair);
475   if (partpair == 0x0) return;
476
477   Double_t weight = partpair->GetWeight();
478   Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
479   Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
480   Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
481   fNumerator->Fill(out,side,lon,weight);
482 }
483 /*************************************************************/
484
485 void AliHBTWeightQOutSQideQLongFctnPerfectPID::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
486 {
487 //Fills numerator
488   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
489   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
490
491   trackpair  = CheckPair(trackpair);
492   if (trackpair == 0x0) return;
493   partpair  = CheckPair(partpair);
494   if (partpair == 0x0) return;
495
496   Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
497   Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
498   Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
499   fDenominator->Fill(out,side,lon);
500 }
501 /******************************************************************/
502
503 TH1* AliHBTWeightQOutSQideQLongFctnPerfectPID::GetResult()
504 {
505  //returns the scaled ratio
506  delete fRatio;
507  fRatio = GetRatio(Scale());
508  return fRatio;
509 }
510
511
512 /******************************************************************/
513 /******************************************************************/
514 /******************************************************************/
515
516 ClassImp(AliHBTQOutSQideQLongFctnPerfectPID)
517
518 AliHBTQOutSQideQLongFctnPerfectPID::AliHBTQOutSQideQLongFctnPerfectPID
519    (Int_t nXbins, Double_t maxXval, Double_t minXval,
520     Int_t nYbins, Double_t maxYval, Double_t minYval,
521     Int_t nZbins, Double_t maxZval, Double_t minZval):
522  AliHBTTwoPairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
523 {
524   //ctor
525   fWriteNumAndDen = kTRUE;//change default behaviour
526   Rename("qoslprfctpid","Q_{out}-Q_{side}-Q_{long} Fctn with perfect PID");
527 }
528 /*************************************************************/
529
530 void AliHBTQOutSQideQLongFctnPerfectPID::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
531 {
532 //Fills numerator
533   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
534   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
535
536   trackpair  = CheckPair(trackpair);
537   if (trackpair == 0x0) return;
538   partpair  = CheckPair(partpair);
539   if (partpair == 0x0) return;
540   
541   Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
542   Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
543   Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
544   fNumerator->Fill(out,side,lon);
545 }
546 /*************************************************************/
547
548 void AliHBTQOutSQideQLongFctnPerfectPID::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
549 {
550 //Fills numerator
551   if (trackpair->Particle1()->GetPid() != partpair->Particle1()->GetPid()) return;
552   if (trackpair->Particle2()->GetPid() != partpair->Particle2()->GetPid()) return;
553
554   trackpair  = CheckPair(trackpair);
555   if (trackpair == 0x0) return;
556   partpair  = CheckPair(partpair);
557   if (partpair == 0x0) return;
558   
559   Double_t out = TMath::Abs(trackpair->GetQOutCMSLC());
560   Double_t side = TMath::Abs(trackpair->GetQSideCMSLC());
561   Double_t lon = TMath::Abs(trackpair->GetQLongCMSLC());
562   fDenominator->Fill(out,side,lon);
563 }
564 /******************************************************************/
565
566 TH1* AliHBTQOutSQideQLongFctnPerfectPID::GetResult()
567 {
568  //returns the scaled ratio
569  delete fRatio;
570  fRatio = GetRatio(Scale());
571  return fRatio;
572 }
573
574 /******************************************************************/
575 /******************************************************************/
576 /******************************************************************/
577
578
579
580 ClassImp(AliHBTPairPIDProbVsQInvFctn)
581
582 AliHBTPairPIDProbVsQInvFctn::AliHBTPairPIDProbVsQInvFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
583  AliHBTOnePairFctn1D(nbins,maxXval,minXval)
584 {
585  fWriteNumAndDen = kTRUE;//change default behaviour
586  Rename("qinvpidpur","Q_{inv}  Function");
587 }
588 /******************************************************************/
589
590 TH1* AliHBTPairPIDProbVsQInvFctn::GetResult()
591 {
592  //returns the scaled ratio
593  delete fRatio;
594  fRatio = GetRatio(Scale());
595  return fRatio;
596 }
597 /******************************************************************/
598
599 void AliHBTPairPIDProbVsQInvFctn::ProcessSameEventParticles(AliHBTPair* pair)
600 {
601  //Fills the numerator using pair from the same event
602    pair = CheckPair(pair);
603    if(pair) fNumerator->Fill(pair->GetQInv(),pair->GetPIDProb());
604 }
605 /******************************************************************/
606
607 void AliHBTPairPIDProbVsQInvFctn::ProcessDiffEventParticles(AliHBTPair* pair)
608  {
609   //Fills the denumerator using mixed pairs
610    pair = CheckPair(pair);
611    if(pair) fDenominator->Fill(pair->GetQInv(),pair->GetPIDProb());
612   }
613
614 /******************************************************************/
615 /******************************************************************/
616 /******************************************************************/
617
618 ClassImp(AliHBTPairPIDProbVsQOutSQideQLongFctn)
619
620 AliHBTPairPIDProbVsQOutSQideQLongFctn::AliHBTPairPIDProbVsQOutSQideQLongFctn(Int_t nXbins, Double_t maxXval, Double_t minXval,
621                                                    Int_t nYbins, Double_t maxYval, Double_t minYval,
622                                                    Int_t nZbins, Double_t maxZval, Double_t minZval):
623  AliHBTOnePairFctn3D(nXbins,maxXval,minXval,nYbins,maxYval,minYval,nZbins,maxZval,minZval)
624 {
625   //ctor
626   fWriteNumAndDen = kTRUE;//change default behaviour
627   Rename("qoslpidpur","Pair PID Probablilty .vs. Q_{out}-Q_{side}-Q_{long} Fctn");
628 }
629 /*************************************************************/
630
631 void AliHBTPairPIDProbVsQOutSQideQLongFctn::ProcessSameEventParticles(AliHBTPair* pair)
632 {
633 //Fills numerator
634   pair  = CheckPair(pair);
635   if (pair == 0x0) return;
636   Double_t weight = pair->GetPIDProb();
637   Double_t out = TMath::Abs(pair->GetQOutCMSLC());
638   Double_t side = TMath::Abs(pair->GetQSideCMSLC());
639   Double_t lon = TMath::Abs(pair->GetQLongCMSLC());
640   fNumerator->Fill(out,side,lon,weight);
641 }
642 /*************************************************************/
643
644 void AliHBTPairPIDProbVsQOutSQideQLongFctn::ProcessDiffEventParticles(AliHBTPair* pair)
645 {
646 //Fills numerator
647   pair  = CheckPair(pair);
648   if (pair == 0x0) return;
649   Double_t weight = pair->GetPIDProb();
650   Double_t out = TMath::Abs(pair->GetQOutCMSLC());
651   Double_t side = TMath::Abs(pair->GetQSideCMSLC());
652   Double_t lon = TMath::Abs(pair->GetQLongCMSLC());
653   fDenominator->Fill(out,side,lon,weight);
654 }
655 /*************************************************************/
656
657 TH1* AliHBTPairPIDProbVsQOutSQideQLongFctn::GetResult()
658 {
659  //returns the scaled ratio
660  
661  delete fRatio;
662  fRatio = GetRatio(Scale());
663  return fRatio;
664 }