Proper filling of correlation function
authorcvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 7 Mar 2013 16:26:19 +0000 (16:26 +0000)
committercvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 7 Mar 2013 16:26:19 +0000 (16:26 +0000)
PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx

index b3abe5f..07db0a7 100644 (file)
@@ -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();