X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGGA%2FGammaConv%2FAliAnalysisTaskGammaConvV1.cxx;h=a83eaf3cf33c63b181e6d0e66f8c743cbb6e4e06;hb=3fe4753a5fca46510ba73a98f18c12482098140f;hp=740fc77cf9e531023f44763559a8fc606ea56471;hpb=614a1866472dcf05cc90661ced7a5843accc1f0d;p=u%2Fmrichter%2FAliRoot.git

diff --git a/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx b/PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx
index 740fc77cf9e..a83eaf3cf33 100644
--- 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){
@@ -1511,7 +1519,7 @@ void AliAnalysisTaskGammaConvV1::ProcessTruePhotonCandidates(AliAODConversionPho
 				if (fMCStack->Particle(fMCStack->Particle(Photon->GetMother(0))->GetMother(0))->GetPdgCode() == 221){
 					iPhotonMCInfo = 3;
 				}
-		}	
+			}	
 		}
 	}
 }