1 #include "AliHBTCorrFitFctn.h"
2 //____________________________________________________________
3 ///////////////////////////////////////////////////////////////
5 // class AliHBTCorrFitFctn //
8 ///////////////////////////////////////////////////////////////
10 ClassImp(AliHBTCorrFitFctn)
12 /*****************************************************************/
14 AliHBTCorrFitFctn::AliHBTCorrFitFctn():
15 AliHBTOnePairFctn1D(),
18 fNMaxPairsFitArea(3000000),
21 fNMaxPairsNormArea(1000000),
26 fWriteNumAndDen = kTRUE;//change default behaviour
27 Rename("wqinvcfCorrFit","Lednicky Weught Theoretical Q_{inv} Correlation Function");
29 /*****************************************************************/
31 AliHBTCorrFitFctn::AliHBTCorrFitFctn(Int_t fit, Int_t norm):
32 AliHBTOnePairFctn1D(100,0.1,0.0),
33 fNtuple(new TNtuple("pair", "pair", "px1:py1:pz1:e1:px2:py2:pz2:e2")),
35 fNMaxPairsFitArea(fit),
38 fNMaxPairsNormArea(norm),
43 fWriteNumAndDen = kTRUE;//change default behaviour
44 Rename("wqinvcfCorrFit","Lednicky Weught Theoretical Q_{inv} Correlation Function");
46 /*****************************************************************/
48 void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* /*trackpair*/)
50 //Fills the numerator using pair from the same event
51 // partpair = CheckPair(partpair);
55 /****************************************************************/
57 void AliHBTCorrFitFctn::ProcessDiffEventParticles(AliHBTPair* trackpair)
59 // Fills the denominator using mixed pairs
60 trackpair = CheckPair(trackpair);
61 if ( trackpair == 0x0) return;
63 Double_t q = 2.* trackpair->GetKStar();
67 if ( q < fFitRangeMax )
69 if (fNPairsFitArea < fNMaxPairsFitArea)
76 Info("ProcessDiffEventParticles","Fit area already full");
80 if ( (q > fNormRangeMin) && (q < fNormRangeMax) )
82 if ( fNPairsNormArea < fNMaxPairsNormArea)
90 const AliVAODParticle& p1 = *(trackpair->Particle1());
91 const AliVAODParticle& p2 = *(trackpair->Particle2());
92 fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.E(),
93 p2.Px(),p2.Py(),p2.Pz(),p2.E());
95 fDenominator->Fill(q);
100 /*****************************************************************/
102 TH1* AliHBTCorrFitFctn::GetResult()
104 //returns denominator
107 /**************************************************************/
109 Int_t AliHBTCorrFitFctn::WriteFunction()
113 retval += AliHBTFunction::WriteFunction();
114 retval += fNtuple->Write(0,TObject::kOverwrite);