Fixed new version

[u/mrichter/AliRoot.git] / PWGGA / GammaConv / AliAnalysisTaskGammaConvV1.cxx

diff --git a/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx b/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx
index 1a48767f82f434bcde008a48081fc0e22d698461..a83eaf3cf33c63b181e6d0e66f8c743cbb6e4e06 100644 (file)
--- a/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx
+++ b/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx
@@ -86,6 +86,7 @@ AliAnalysisTaskGammaConvV1::AliAnalysisTaskGammaConvV1(): AliAnalysisTaskSE(),
        hESDConvGammaPt(NULL),
        hESDConvGammaR(NULL),
        hESDConvGammaEta(NULL),
+       hESDConvGammaPhi(NULL),
        tESDConvGammaPtDcazCat(NULL),
        fPtGamma(0),
        fDCAzPhoton(0),
@@ -156,9 +157,9 @@ AliAnalysisTaskGammaConvV1::AliAnalysisTaskGammaConvV1(): AliAnalysisTaskSE(),
        hESDTrueEtaPtOpenAngle(NULL),
        hESDTrueMotherDalitzInvMassPt(NULL),
        hESDTrueConvGammaPt(NULL),
-        hESDTrueConvGammaR(NULL),
+       hESDTrueConvGammaR(NULL),
        hESDTrueConvGammaPtMC(NULL),
-        hESDTrueConvGammaRMC(NULL),
+       hESDTrueConvGammaRMC(NULL),
        hESDTrueConvGammaEta(NULL),
        hESDCombinatorialPt(NULL),
        hESDTruePrimaryConvGammaPt(NULL),
@@ -235,6 +236,7 @@ AliAnalysisTaskGammaConvV1::AliAnalysisTaskGammaConvV1(const char *name):
        hESDConvGammaPt(NULL),
        hESDConvGammaR(NULL),
        hESDConvGammaEta(NULL),
+       hESDConvGammaPhi(NULL),
        tESDConvGammaPtDcazCat(NULL),
        fPtGamma(0),
        fDCAzPhoton(0),
@@ -475,6 +477,7 @@ void AliAnalysisTaskGammaConvV1::UserCreateOutputObjects(){
        if (fDoPhotonQA > 0){
                hESDConvGammaR = new TH1F*[fnCuts];
                hESDConvGammaEta = new TH1F*[fnCuts];
+               hESDConvGammaPhi = new TH1F*[fnCuts];
        } 
        
        if(fDoMesonAnalysis){
@@ -577,6 +580,8 @@ void AliAnalysisTaskGammaConvV1::UserCreateOutputObjects(){
                        fESDList[iCut]->Add(hESDConvGammaR[iCut]);
                        hESDConvGammaEta[iCut] = new TH1F("ESD_ConvGamma_Eta","ESD_ConvGamma_Eta",2000,-2,2);
                        fESDList[iCut]->Add(hESDConvGammaEta[iCut]);
+                       hESDConvGammaPhi[iCut] = new TH1F("ESD_ConvGamma_Phi","ESD_ConvGamma_Phi",360,0,2*TMath::Pi());
+                       fESDList[iCut]->Add(hESDConvGammaPhi[iCut]);
                }
                
                if(fDoMesonAnalysis){
@@ -1184,8 +1189,9 @@ void AliAnalysisTaskGammaConvV1::ProcessPhotonCandidates()
                        if(fIsFromMBHeader){
                                hESDConvGammaPt[fiCut]->Fill(PhotonCandidate->Pt());
                                if (fDoPhotonQA > 0){
-                               hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
-                               hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
+                                       hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
+                                       hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
+                                       hESDConvGammaPhi[fiCut]->Fill(PhotonCandidate->Phi());
                                }   
                        }
                        if(fIsMC){
@@ -1236,11 +1242,12 @@ void AliAnalysisTaskGammaConvV1::ProcessPhotonCandidates()
                        if(!((AliConversionPhotonCuts*)fCutArray->At(fiCut))->UseToCloseV0sCut()){ // To Colse v0s cut diabled, step two not needed
                                fGammaCandidates->Add(PhotonCandidate);
                                if(fIsFromMBHeader){
-                               hESDConvGammaPt[fiCut]->Fill(PhotonCandidate->Pt());
-                               if (fDoPhotonQA > 0){
-                                       hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
-                                       hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
-                               }
+                                       hESDConvGammaPt[fiCut]->Fill(PhotonCandidate->Pt());
+                                       if (fDoPhotonQA > 0){
+                                               hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
+                                               hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
+                                               hESDConvGammaPhi[fiCut]->Fill(PhotonCandidate->Phi());
+                                       }
                                }
                        }
                        if(fIsMC){
@@ -1286,8 +1293,9 @@ void AliAnalysisTaskGammaConvV1::ProcessPhotonCandidates()
                        if(fIsFromMBHeader){
                                hESDConvGammaPt[fiCut]->Fill(PhotonCandidate->Pt());
                                if (fDoPhotonQA > 0){
-                               hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
-                               hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
+                                       hESDConvGammaR[fiCut]->Fill(PhotonCandidate->GetConversionRadius());
+                                       hESDConvGammaEta[fiCut]->Fill(PhotonCandidate->Eta());
+                                       hESDConvGammaPhi[fiCut]->Fill(PhotonCandidate->Phi());
                                }
                        }
                        if(fIsMC){