Proper filling of correlation function

[u/mrichter/AliRoot.git] / PWGCF / Correlations / DPhi / MuonHadron / AliAnalysisTaskDiMuonCorrelations.cxx

diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx
index b3abe5fdb59eb94c04a5debb4dcc902decb3d20f..07db0a7b2d5e6e9afda29b82fd0c5b4b0e55f725 100644 (file)
--- a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx
+++ b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx
@@ -330,7 +330,8 @@ void AliAnalysisTaskDiMuonCorrelations::UserExec(Option_t *) {
     fMuonTrack[0] = (AliAODTrack*) tracksMuonArm->At(iTrMuon1);
     fHistSingleMuonsPt[centBin]->Fill(fMuonTrack[0]->Pt());
     fHistSingleMuonsEtaPt[centBin]->Fill(fMuonTrack[0]->Pt(),fMuonTrack[0]->Eta());
-    for (Int_t iTrMuon2=iTrMuon1+1; iTrMuon2<tracksMuonArm->GetEntriesFast(); iTrMuon2++) {
+    for (Int_t iTrMuon2=0; iTrMuon2<tracksMuonArm->GetEntriesFast(); iTrMuon2++) {
+      if (iTrMuon2 == iTrMuon1) continue;
       fMuonTrack[1] = (AliAODTrack*) tracksMuonArm -> At(iTrMuon2);
       FillHistograms(centBin, kSingleEvent);
     }
@@ -376,6 +377,8 @@ void AliAnalysisTaskDiMuonCorrelations::FillHistograms(Int_t centrality, Int_t o
 
   if (ptBinTrackMuon1<1 || ptBinTrackMuon1>fNbinsPt || ptBinTrackMuon2<1 || ptBinTrackMuon2>fNbinsPt) return;
 
+  if (ptBinTrackMuon1 < ptBinTrackMuon2) return;
+
   Double_t deltaPhi = fMuonTrack[0]->Phi() - fMuonTrack[1]->Phi();
   if (deltaPhi >  1.5*TMath::Pi()) deltaPhi -= TMath::TwoPi();
   if (deltaPhi < -0.5*TMath::Pi()) deltaPhi += TMath::TwoPi();