added histogram for cross check of delta phi vs centrality

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
index 0c14648da3ad5d2c3aa9e7c46e3e0d23af57ac9b..91084fd5e3dcc461963fec261b21cf41bf52e5c1 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
@@ -83,6 +83,8 @@ fPcosEPCent(0),
 fPsinEPCent(0),
 fPcosPhiCent(0),
 fPsinPhiCent(0),
+// cross check for event plane determination
+fDeltaPhiCent(0),
 //global
 fIsMonteCarlo(0),
 fEPselector("Q"),
@@ -605,6 +607,11 @@ void AlidNdPtAnalysisPbPbAOD::UserCreateOutputObjects()
   fPsinPhiCent->GetYaxis()->SetTitle("#LT sin 2 #phi #GT");
   fPsinPhiCent->Sumw2();
   
+  fDeltaPhiCent = new TH2F("fDeltaPhiCent","fDeltaPhiCent",200, -2.*TMath::Pi(), 2.*TMath::Pi(), fCentralityNbins-1, fBinsCentrality);
+  fDeltaPhiCent->GetXaxis()->SetTitle("#Delta #phi");
+  fDeltaPhiCent->GetYaxis()->SetTitle("Centrality");
+  fDeltaPhiCent->Sumw2();
+  
   // Add Histos, Profiles etc to List
   fOutputList->Add(fZvPtEtaCent);
   fOutputList->Add(fDeltaphiPtEtaCent);
@@ -656,6 +663,8 @@ void AlidNdPtAnalysisPbPbAOD::UserCreateOutputObjects()
   fOutputList->Add(fPcosPhiCent);
   fOutputList->Add(fPsinPhiCent);
   
+  fOutputList->Add(fDeltaPhiCent);
+  
   StoreCutSettingsToHistogram();
   
   PostData(1, fOutputList);
@@ -986,6 +995,7 @@ void AlidNdPtAnalysisPbPbAOD::UserExec(Option_t *option)
        dTrackPhiPtEtaCent[2] = track->Eta();
        dTrackPhiPtEtaCent[3] = dCentrality;
        
+       
        if( fIsMonteCarlo )
        {
          mcPart = (AliAODMCParticle*)stack->At(TMath::Abs(track->GetLabel()));
@@ -1073,6 +1083,11 @@ void AlidNdPtAnalysisPbPbAOD::UserExec(Option_t *option)
          fPhiCent->Fill(track->Phi(), dCentrality);
          fPcosPhiCent->Fill(dCentrality, TMath::Cos(2.*track->Phi()));
          fPsinPhiCent->Fill(dCentrality, TMath::Sin(2.*track->Phi()));
+         
+         Double_t deltaphi = track->Phi() - dEventplaneAngleCorrected;
+         if(deltaphi > TMath::Pi()) deltaphi -= 2.*TMath::Pi();
+         
+         fDeltaPhiCent->Fill(deltaphi, dCentrality);
        }
   } // end track loop
   

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
index ef6f8f7a92ca629966da57dd713ce898fa7d2de5..17a9708daf567967fda10e92b182d8a13a7315eb 100644 (file)
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h
@@ -182,7 +182,7 @@ class AlidNdPtAnalysisPbPbAOD : public AliAnalysisTaskSE {
     TH1F           *fPt; // simple pT histogramm
     TH1F           *fMCPt; // simple pT truth histogramm
     THnSparseF         *fZvPtEtaCent; //-> Zv:Pt:Eta:Cent
-    THnSparseF         *fDeltaphiPtEtaCent; //-> Phi:Pt:Eta:Cent
+    THnSparseF         *fDeltaphiPtEtaCent; //-> DeltaPhi:Pt:Eta:Cent
     THnSparseF         *fPtResptCent; //-> 1/pt:ResolutionPt:Cent
     THnSparseF         *fMCRecPrimZvPtEtaCent; //-> MC Zv:Pt:Eta:Cent
     THnSparseF         *fMCGenZvPtEtaCent; //-> MC Zv:Pt:Eta:Cent
@@ -228,6 +228,9 @@ class AlidNdPtAnalysisPbPbAOD : public AliAnalysisTaskSE {
        TProfile        *fPcosPhiCent; // < cos 2 phi > vs centrality
        TProfile        *fPsinPhiCent; // < sin 2 phi > vs centrality
 
+       // cross check for event plane determination
+       TH2F            *fDeltaPhiCent; // DeltaPhi:Cent - DeltaPhi in the range from -pi to pi
+
        // global variables
     Bool_t fIsMonteCarlo;
        
@@ -292,7 +295,7 @@ class AlidNdPtAnalysisPbPbAOD : public AliAnalysisTaskSE {
     AlidNdPtAnalysisPbPbAOD(const AlidNdPtAnalysisPbPbAOD&); // not implemented
     AlidNdPtAnalysisPbPbAOD& operator=(const AlidNdPtAnalysisPbPbAOD&); // not implemented  
     
-    ClassDef(AlidNdPtAnalysisPbPbAOD,10); // has to be at least 1, otherwise not streamable...
+    ClassDef(AlidNdPtAnalysisPbPbAOD,11); // has to be at least 1, otherwise not streamable...
 };
 
 #endif