X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=HBTAN%2FAliHBTCorrFitFctn.cxx;h=b84d5e255436dba4274ecef7fd8ae9456483b695;hp=66d4edbc3f792ea0a0cbba55768cd3ef9f3dd96e;hb=0bb0710808d5c005df0875008431072121ca6177;hpb=12ba5510c60191c882c2e2299ade3cf43fcf552f diff --git a/HBTAN/AliHBTCorrFitFctn.cxx b/HBTAN/AliHBTCorrFitFctn.cxx index 66d4edbc3f7..b84d5e25543 100644 --- a/HBTAN/AliHBTCorrFitFctn.cxx +++ b/HBTAN/AliHBTCorrFitFctn.cxx @@ -31,10 +31,15 @@ void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTP //Fills the numerator using pair from the same event partpair = CheckPair(partpair); if(partpair == 0x0) return; + trackpair = CheckPair(trackpair); + if(trackpair == 0x0) return; Double_t q = trackpair->GetQInv(); Bool_t fill = kFALSE; + Double_t weight = AliHBTLLWeights::Instance()->GetWeight(partpair); + fNumerator->Fill(q,weight); + if ( (q < 0.15) && (fNPairsFitArea < 2.e+5)) { fNPairsFitArea++; @@ -49,12 +54,10 @@ void AliHBTCorrFitFctn::ProcessSameEventParticles(AliHBTPair* trackpair, AliHBTP if (fill) { - Double_t weight = AliHBTLLWeights::Instance()->GetWeight(partpair); - fNumerator->Fill(q,weight); const AliHBTParticle& p1 = *(trackpair->Particle1()); const AliHBTParticle& p2 = *(trackpair->Particle2()); fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.Energy(), - p2.Px(),p2.Py(),p2.Pz(),p2.Energy()); + p2.Px(),p2.Py(),p2.Pz(),p2.Energy()); } } /****************************************************************/