updated to check fake issue

diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
index fac67977d2397b7690014fe353d47c760195c5ba..ed124cb2d7eacd039e20e282ea59cdc38757f258 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
@@ -212,6 +212,8 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal(const char *name)
   ,fIncRecoMaxE(0)
   ,fPhoRecoMaxE(0)
   ,fSamRecoMaxE(0) 
+  ,fPhoVertexReco0(0)
+  ,fPhoVertexReco1(0)
   //,fnSigEtaCorr(NULL)
 {
   //Named constructor
@@ -356,6 +358,8 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal()
   ,fIncRecoMaxE(0)
   ,fPhoRecoMaxE(0)
   ,fSamRecoMaxE(0)
+  ,fPhoVertexReco0(0)
+  ,fPhoVertexReco1(0)
   //,fnSigEtaCorr(NULL)
 {
        //Default constructor
@@ -558,6 +562,7 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
     double mcMompT = 0.0;
     //double mcGrandMompT = 0.0;
     double mcWeight = -10.0;
+    double conv_proR = -1.0;
 
     int iHijing = 1;
     int mcLabel = -1;
@@ -576,6 +581,7 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
               TParticle* particle = stack->Particle(label);
               int mcpid = particle->GetPdgCode();
               mcpT = particle->Pt();
+               conv_proR = particle->R();
               //printf("MCpid = %d",mcpid);
               if(particle->GetFirstMother()>-1)
               {
@@ -977,13 +983,16 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
                   if(fFlagPhotonicElec) fPhoElecPtMCM20_pi0e->Fill(phoval,mcWeight);
                   if(fFlagConvinatElec) fSameElecPtMCM20_pi0e->Fill(phoval,mcWeight);
                  }
+               // --- eta
                if(mcOrgEta)
                  {
                   fIncpTMCM20pho_eta->Fill(phoval,mcWeight);    
                   if(fFlagPhotonicElec) fPhoElecPtMCM20_eta->Fill(phoval,mcWeight);
                   if(fFlagConvinatElec) fSameElecPtMCM20_eta->Fill(phoval,mcWeight);
                  }
-               // --- eta
+                // check production vertex
+                fPhoVertexReco0->Fill(pt,conv_proR);
+                if(fFlagPhotonicElec) fPhoVertexReco1->Fill(pt,conv_proR);
               }
            }
         } 
@@ -1451,6 +1460,12 @@ void AliAnalysisTaskHFECal::UserCreateOutputObjects()
   fSamRecoMaxE = new TH2D("fSamRecoMaxE","Same",10,0,100,100,0,500);
   fOutputList->Add(fSamRecoMaxE);
 
+  fPhoVertexReco0 = new TH2D("fPhoVertexReco0","photon production Vertex in ele sample",40,0,20,200,0,40);
+  fOutputList->Add(fPhoVertexReco0);
+
+  fPhoVertexReco1 = new TH2D("fPhoVertexReco1","photon production Vertex mass selection",40,0,20,200,0,40);
+  fOutputList->Add(fPhoVertexReco1);
+
   PostData(1,fOutputList);
 }
 

diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.h b/PWGHF/hfe/AliAnalysisTaskHFECal.h
index 88953c469eda1ef11c71c48308874ad98ddf8ad2..469e756dbc1d8a88b6944e1563096a31674767b0 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.h
@@ -208,6 +208,8 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
  TH2D                   *fIncRecoMaxE;
  TH2D                   *fPhoRecoMaxE;
  TH2D                   *fSamRecoMaxE; 
+ TH2D                   *fPhoVertexReco0;
+ TH2D                   *fPhoVertexReco1;
 
  //<----- correction
  TGraphErrors           *fnSigEtaCorr[7];