From 60544aea18090d5f33f525d46c3725e254139259 Mon Sep 17 00:00:00 2001
From: ssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 15 Oct 2012 10:10:12 +0000
Subject: [PATCH] updated

---
 PWGHF/hfe/AliAnalysisTaskHFECal.cxx | 12 +++++++++---
 PWGHF/hfe/AliAnalysisTaskHFECal.h   |  1 +
 2 files changed, 10 insertions(+), 3 deletions(-)

diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
index 8da91cee081..2804acefe6e 100644
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
@@ -144,6 +144,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal(const char *name)
   ,fIncpTMChfeAll(0)	
   ,fIncpTMCM20hfe(0)	
   ,fIncpTMCM20hfeAll(0)	
+  ,fIncpTMCM20hfeCheck(0)	
   ,fIncpTMCpho(0)	
   ,fIncpTMCM20pho(0)	
   ,fPhoElecPtMC(0)
@@ -237,6 +238,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal()
   ,fIncpTMChfeAll(0)	
   ,fIncpTMCM20hfe(0)	
   ,fIncpTMCM20hfeAll(0)	
+  ,fIncpTMCM20hfeCheck(0)	
   ,fIncpTMCpho(0)	
   ,fIncpTMCM20pho(0)	
   ,fPhoElecPtMC(0)
@@ -650,6 +652,7 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
          {
           fIncpTMChfe->Fill(cent,pt);    
           if(m20>0.0 && m20<0.3)fIncpTMCM20hfe->Fill(cent,pt);    
+          if(m20>0.0 && m20<0.3)fIncpTMCM20hfeCheck->Fill(cent,mcpT);    
          }
       
        if(mcPho) // select photonic electrons
@@ -914,6 +917,8 @@ void AliAnalysisTaskHFECal::UserCreateOutputObjects()
   fIncpTMCM20hfeAll = new TH2F("fIncpTMCM20hfeAll","MC Alle pid electro vs. centrality with M20",200,0,100,100,0,50);
   fOutputList->Add(fIncpTMCM20hfeAll);
 
+  fIncpTMCM20hfeCheck = new TH2F("fIncpTMCM20hfeCheck","MC HFE pid electro vs. centrality with M20 Check",200,0,100,100,0,50);
+  fOutputList->Add(fIncpTMCM20hfeCheck);
 
   Int_t nBinspho2[5] =  { 200, 100,    7,    3, 100};
   Double_t minpho2[5] = {  0.,  0., -2.5, -0.5, 0.};   
@@ -1065,12 +1070,13 @@ void AliAnalysisTaskHFECal::SelectPhotonicElectron(Int_t itrack, Double_t cent,
     Double_t chi2recg = recg.GetChi2()/recg.GetNDF();
     if(TMath::Sqrt(TMath::Abs(chi2recg))>3.) continue;
     
+    /*  for v5-03-70-AN
     AliKFVertex primV(*pVtx);
     primV += recg;
     recg.SetProductionVertex(primV);
     
     recg.SetMassConstraint(0,0.0001);
-    
+    */
     openingAngle = ge1.GetAngle(ge2);
     if(fFlagLS) fOpeningAngleLS->Fill(openingAngle);
     if(fFlagULS) fOpeningAngleULS->Fill(openingAngle);
@@ -1102,8 +1108,8 @@ void AliAnalysisTaskHFECal::SelectPhotonicElectron(Int_t itrack, Double_t cent,
 
     if(openingAngle > fOpeningAngleCut) continue;
       
-    //if(mass<fInvmassCut && fFlagULS && !flagPhotonicElec){
-    if(mass<fInvmassCut && fFlagULS && !flagPhotonicElec && (p1==p2)){
+    if(mass<fInvmassCut && fFlagULS && !flagPhotonicElec){
+    //if(mass<fInvmassCut && fFlagULS && !flagPhotonicElec && (p1==p2)){ <--- only MC train (55,56) v5-03-68-AN & 69 for check
       flagPhotonicElec = kTRUE;
     }
     if(mass<fInvmassCut && fFlagLS && !flagConvinatElec){
diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.h b/PWGHF/hfe/AliAnalysisTaskHFECal.h
index ae09b2ead55..a40e0a769ff 100644
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.h
@@ -133,6 +133,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
   TH2F			*fIncpTMChfeAll;		//! MC HFE pid electron vs centrality
   TH2F			*fIncpTMCM20hfe;	//! MC HFE pid electron vs centrality
   TH2F			*fIncpTMCM20hfeAll;	//! MC HFE pid electron vs centrality
+  TH2F			*fIncpTMCM20hfeCheck;	//! MC HFE pid electron vs centrality
  THnSparseD		*fIncpTMCpho;		//! MC HFE pid electron vs centrality
  THnSparseD 		*fIncpTMCM20pho;	//! MC HFE pid electron vs centrality
  THnSparseD 		*fPhoElecPtMC;	        //! Pho inclusive ele pt
-- 
2.43.0