added histogram for cross check of delta phi vs centrality
authorpluettig <philipp.luettig@cern.ch>
Mon, 8 Sep 2014 13:52:23 +0000 (15:52 +0200)
committerpluettig <philipp.luettig@cern.ch>
Mon, 8 Sep 2014 13:52:23 +0000 (15:52 +0200)
PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.cxx
PWGLF/SPECTRA/ChargedHadrons/dNdPt/AlidNdPtAnalysisPbPbAOD.h

index 0c14648..91084fd 100644 (file)
@@ -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
   
index ef6f8f7..17a9708 100644 (file)
@@ -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