1 #include "AliHBTCorrFitFctn.h"
2 //____________________________________________________________
3 ///////////////////////////////////////////////////////////////
5 // class AliHBTCorrFitFctn //
8 ///////////////////////////////////////////////////////////////
10 ClassImp(AliHBTCorrFitFctn)
12 /*****************************************************************/
14 AliHBTCorrFitFctn::AliHBTCorrFitFctn(Int_t nbins, Double_t maxXval, Double_t minXval):
15 AliHBTTwoPairFctn1D(nbins,maxXval,minXval),
16 fNtuple(new TNtuple("pair", "pair", "px1:py1:pz1:e1:px2:py2:pz2:e2")),
21 fWriteNumAndDen = kTRUE;//change default behaviour
22 Rename("wqinvcfCorrFit","Lednicky Weught Theoretical Q_{inv} Correlation Function");
24 /*****************************************************************/
26 void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
28 //Fills the numerator using pair from the same event
29 partpair = CheckPair(partpair);
30 if(partpair == 0x0) return;
31 trackpair = CheckPair(trackpair);
32 if(trackpair == 0x0) return;
34 Double_t q = trackpair->GetQInv();
37 Double_t weight = partpair->GetLLWeight();
38 fNumerator->Fill(q,weight);
40 if ( (q < 0.15) && (fNPairsFitArea < 2.e+5))
46 if ( (q > 0.15) && (q < 0.3) && (fNPairsFitArea < 1.e+5))
54 const AliHBTParticle& p1 = *(trackpair->Particle1());
55 const AliHBTParticle& p2 = *(trackpair->Particle2());
56 fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.Energy(),
57 p2.Px(),p2.Py(),p2.Pz(),p2.Energy());
60 /****************************************************************/
62 void AliHBTCorrFitFctn::ProcessDiffEventParticles(AliHBTPair* trackpair, AliHBTPair* partpair)
64 // Fills the denominator using mixed pairs
65 trackpair = CheckPair(trackpair);
66 partpair = CheckPair(partpair);
67 if ( trackpair && partpair)
69 fDenominator->Fill(trackpair->GetQInv());
72 /*****************************************************************/
74 TH1* AliHBTCorrFitFctn::GetResult()
76 //returns ratio of numerator and denominator
77 return GetRatio(Scale());
79 /**************************************************************/
81 void AliHBTCorrFitFctn::WriteFunction()
84 AliHBTFunction::WriteFunction();
85 fNtuple->Write(0,TObject::kOverwrite);