forward muon selection added

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
index 3edb82e313a170e39021382ee7c223e984683a43..6056fe1abc1a4e8b1fac5fd4b0e6af6816e2f905 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
@@ -1153,7 +1153,7 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
   
   if (fTriggersFromDetector == 0)
     tracks = fAnalyseUE->GetAcceptedParticles(inputEvent, 0, kTRUE, fParticleSpeciesTrigger, kTRUE);
   
   if (fTriggersFromDetector == 0)
     tracks = fAnalyseUE->GetAcceptedParticles(inputEvent, 0, kTRUE, fParticleSpeciesTrigger, kTRUE);

-  else if (fTriggersFromDetector == 1 || fTriggersFromDetector == 2 || fTriggersFromDetector == 3)
+  else if (fTriggersFromDetector <= 4)

     tracks=GetParticlesFromDetector(inputEvent,fTriggersFromDetector);
   else
     AliFatal(Form("Invalid setting for fTriggersFromDetector: %d", fTriggersFromDetector));
     tracks=GetParticlesFromDetector(inputEvent,fTriggersFromDetector);
   else
     AliFatal(Form("Invalid setting for fTriggersFromDetector: %d", fTriggersFromDetector));


@@ -1208,7 +1208,7 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
     if (fParticleSpeciesAssociated != fParticleSpeciesTrigger || fTriggersFromDetector > 0 )
       tracksCorrelate = fAnalyseUE->GetAcceptedParticles(inputEvent, 0, kTRUE, fParticleSpeciesAssociated, kTRUE);
   }
     if (fParticleSpeciesAssociated != fParticleSpeciesTrigger || fTriggersFromDetector > 0 )
       tracksCorrelate = fAnalyseUE->GetAcceptedParticles(inputEvent, 0, kTRUE, fParticleSpeciesAssociated, kTRUE);
   }

-  else if (fAssociatedFromDetector == 1 || fAssociatedFromDetector == 2 || fAssociatedFromDetector == 3){
+  else if (fAssociatedFromDetector <= 4){

     if(fAssociatedFromDetector != fTriggersFromDetector)
       tracksCorrelate = GetParticlesFromDetector(inputEvent,fAssociatedFromDetector);
   }
     if(fAssociatedFromDetector != fTriggersFromDetector)
       tracksCorrelate = GetParticlesFromDetector(inputEvent,fAssociatedFromDetector);
   }


@@ -1507,6 +1507,27 @@ TObjArray* AliAnalysisTaskPhiCorrelations::GetParticlesFromDetector(AliVEvent* i
          obj->Add(particle);
                }      
     }
          obj->Add(particle);
                }      
     }

+  else if (idet == 4)
+    {
+      if(!fAOD)
+       AliFatal("Muon selection only implemented on AOD");//FIXME to be implemented also for ESDs as in AliAnalyseLeadingTrackUE::GetAcceptedPArticles
+      for (Int_t iTrack = 0; iTrack < fAOD->GetNTracks(); iTrack++) {
+       AliAODTrack* track = fAOD->GetTrack(iTrack);
+       if (!track->IsMuonTrack()) continue;
+       //Float_t dca    = track->DCA();
+       //Float_t chi2   = track->Chi2perNDF();
+       //Int_t   mask   = track->GetMatchTrigger();
+       Float_t eta    = track->Eta();
+       Float_t rabs   = track->GetRAtAbsorberEnd();
+       if (rabs < 17.6 || rabs > 89.5) continue;
+       if (eta < -4 || eta > -2.5) continue;
+       
+       AliDPhiBasicParticle* particle = new AliDPhiBasicParticle(eta,track->Phi(), track->Pt(), track->Charge()); 
+       particle->SetUniqueID(fAnalyseUE->GetEventCounter() * 100000 + iTrack);
+       
+       obj->Add(particle);
+      }
+    }      

   else
     AliFatal(Form("GetParticlesFromDetector: Invalid idet value: %d", idet));
   
   else
     AliFatal(Form("GetParticlesFromDetector: Invalid idet value: %d", idet));
   

diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
index c5882aa5ceadc717acd643df0e9401067a713d97..d08337508b3929185289657ac2f3089259b944cd 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
+++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h
@@ -239,8 +239,8 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     Bool_t fWeightPerEvent;       // weight with the number of trigger particles per event
     TString fCustomBinning;       // supersedes default binning if set, see AliUEHist::GetBinning or AliUEHistograms::AliUEHistograms for syntax and examples
     Bool_t fPtOrder;              // apply pT,a < pt,t condition; default: kTRUE
     Bool_t fWeightPerEvent;       // weight with the number of trigger particles per event
     TString fCustomBinning;       // supersedes default binning if set, see AliUEHist::GetBinning or AliUEHistograms::AliUEHistograms for syntax and examples
     Bool_t fPtOrder;              // apply pT,a < pt,t condition; default: kTRUE

-    Int_t fTriggersFromDetector;   // 0 = tracks (default); 1 = VZERO_A; 2 = VZERO_C; 3 = SPD tracklets
-    Int_t fAssociatedFromDetector;   // 0 = tracks (default); 1 = VZERO_A; 2 = VZERO_C; 3 = SPD tracklets
+    Int_t fTriggersFromDetector;   // 0 = tracks (default); 1 = VZERO_A; 2 = VZERO_C; 3 = SPD tracklets; 4 = forward muons
+    Int_t fAssociatedFromDetector;   // 0 = tracks (default); 1 = VZERO_A; 2 = VZERO_C; 3 = SPD tracklets; 4 = forward muons

     Bool_t fMCUseUncheckedCentrality; // use unchecked centrality (only applies to MC); default: kFALSE
     
     Bool_t fFillpT;                // fill sum pT instead of number density
     Bool_t fMCUseUncheckedCentrality; // use unchecked centrality (only applies to MC); default: kFALSE
     
     Bool_t fFillpT;                // fill sum pT instead of number density