Update EMCal pPb

diff --git a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx
index bc1421709a5df03317d5035194fe2b792ba00570..f06fb4c5ec145592f93c6b616c31d2a7952d6645 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx
@@ -19,7 +19,7 @@
        //      Task for Heavy-flavour electron analysis in pPb collisions    //
        //      (+ Electron-Hadron Jetlike Azimuthal Correlation)             //
        //                                                                                                                                        //
-       //              version: June 18, 2014.                                                               //
+       //              version: July 28, 2014.                                                               //
        //                                                                    //
        //          Authors                                                                               //
        //              Elienos Pereira de Oliveira Filho (epereira@cern.ch)          //
@@ -304,7 +304,9 @@ AliAnalysisTaskEMCalHFEpA::AliAnalysisTaskEMCalHFEpA(const char *name)
 ,fM02CutMax(10)
 ,fAngleCut(999)
 ,fChi2Cut(3.5)
-,fDCAcut(999)
+,fDCAcut(999)//dca between two tracks
+,fDCAcutr(999)//dca to vertex
+,fDCAcutz(999)//dca to vertex
 ,fMassCutFlag(kTRUE)
 ,fAngleCutFlag(kFALSE)
 ,fChi2CutFlag(kFALSE)
@@ -590,7 +592,9 @@ AliAnalysisTaskEMCalHFEpA::AliAnalysisTaskEMCalHFEpA()
 ,fM02CutMax(10)
 ,fAngleCut(999)
 ,fChi2Cut(3.5)
-,fDCAcut(999)
+,fDCAcut(999)//dca between two tracks
+,fDCAcutr(999)//dca to vertex
+,fDCAcutz(999)//dca to vertex
 ,fMassCutFlag(kTRUE)
 ,fAngleCutFlag(kFALSE)
 ,fChi2CutFlag(kFALSE)
@@ -984,7 +988,7 @@ void AliAnalysisTaskEMCalHFEpA::UserCreateOutputObjects()
        
        for(Int_t i = 0; i < 3; i++)
        {
-               fEoverP_pt[i] = new TH2F(Form("fEoverP_pt%d",i),";p_{t} (GeV/c);E / p ",1000,0,30,500,0,2);
+               fEoverP_pt[i] = new TH2F(Form("fEoverP_pt%d",i),";p_{t} (GeV/c);E / p ",1000,0,30,2000,0,2);
                fTPC_p[i] = new TH2F(Form("fTPC_p%d",i),";pt (GeV/c);TPC dE/dx (a. u.)",1000,0.3,15,1000,-20,200);
                fTPCnsigma_p[i] = new TH2F(Form("fTPCnsigma_p%d",i),";p (GeV/c);TPC Electron N#sigma",1000,0.3,15,1000,-15,10);
                
@@ -1244,13 +1248,13 @@ void AliAnalysisTaskEMCalHFEpA::UserCreateOutputObjects()
        
        
        fNcells_pt=new TH2F("fNcells_pt","fNcells_pt",1000, 0,20,100,0,30);
-       fEoverP_pt_pions= new TH2F("fEoverP_pt_pions","fEoverP_pt_pions",1000,0,30,500,0,2);
+       fEoverP_pt_pions= new TH2F("fEoverP_pt_pions","fEoverP_pt_pions",1000,0,30,2000,0,2);
        
        ftpc_p_EoverPcut= new TH2F("ftpc_p_EoverPcut","ftpc_p_EoverPcut",1000,0,30,200,20,200);
        fnsigma_p_EoverPcut= new TH2F("fnsigma_p_EoverPcut","fnsigma_p_EoverPcut",1000,0,30,500,-15,15);
        
-       fEoverP_pt_pions2= new TH2F("fEoverP_pt_pions2","fEoverP_pt_pions2",1000,0,30,500,0,2);
-       fEoverP_pt_hadrons= new TH2F("fEoverP_pt_hadrons","fEoverP_pt_hadrons",1000,0,30,500,0,2);
+       fEoverP_pt_pions2= new TH2F("fEoverP_pt_pions2","fEoverP_pt_pions2",1000,0,30,2000,0,2);
+       fEoverP_pt_hadrons= new TH2F("fEoverP_pt_hadrons","fEoverP_pt_hadrons",1000,0,30,2000,0,2);
        
        
        fOutputList->Add(fTPCnsigma_eta);
@@ -2353,20 +2357,41 @@ void AliAnalysisTaskEMCalHFEpA::UserExec(Option_t *)
                } 
                
 //RecPrim
-               if(!fIsAOD)
-               {
-                       if(!ProcessCutStep(AliHFEcuts::kStepRecPrim, track)) continue;
-               }
+               
+                       //it was not working on aod... testing again
+                       //July 29th, 2014: aparently working again
+               
+///if(!fIsAOD)
+//{
+
+               if(!ProcessCutStep(AliHFEcuts::kStepRecPrim, track)) continue;
+//}
                
 //HFEcuts: ITS layers cuts
                if(!ProcessCutStep(AliHFEcuts::kStepHFEcutsITS, track)) continue;
                
 //HFE cuts: TPC PID cleanup
                if(!ProcessCutStep(AliHFEcuts::kStepHFEcutsTPC, track)) continue;
+               
+//DCA cut done by hand -- July 29th, 2014
+               if(fIsAOD){
+                       Double_t d0z0[2], cov[3];
+                       AliAODVertex *prim_vtx = fAOD->GetPrimaryVertex();
+                       track->PropagateToDCA(prim_vtx, fAOD->GetMagneticField(), 20., d0z0, cov);
+                       Double_t DCAxy = d0z0[0];
+                       Double_t DCAz = d0z0[1];
+                       
+                       if(DCAxy >= fDCAcutr || DCAz>=fDCAcutz) continue;
+               }
+               
+               
+               
+               
+               
 //______________________________________________________________________________________
                
                if(fIsAOD){     
-                               //AOD test -- Fancesco suggestion
+                               //AOD test -- Francesco suggestion
                                //aod test -- Francesco suggestion
                        AliAODTrack *aod_track=fAOD->GetTrack(iTracks);
 

diff --git a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.h b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.h
index af168bfc94d83f07a6955398d591801032db3bea..4c50654594e69ba366714775095d92aac3d2838c 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.h
+++ b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.h
@@ -9,7 +9,7 @@
        //      Task for Heavy-flavour electron analysis in pPb collisions    //
        //      (+ Electron-Hadron Jetlike Azimuthal Correlation)             //
        //                                                                                                                                        //
-       //              version: June 16, 2014.                                                               //
+       //              version: July 28, 2014.                                                               //
        //                                                                    //
        //          Authors                                                                               //
        //              Elienos Pereira de Oliveira Filho (epereira@cern.ch)          //
@@ -106,6 +106,10 @@ public:
        void SetNonHFEangleCut(Double_t AngleCut) { fAngleCut = AngleCut; fAngleCutFlag = kTRUE;};
        void SetNonHFEchi2Cut(Double_t Chi2Cut) { fChi2Cut = Chi2Cut; fChi2CutFlag = kTRUE;};
        void SetNonHFEdcaCut(Double_t DCAcut) { fDCAcut = DCAcut; fDCAcutFlag = kTRUE;};
+       
+       //DCA cut main particle
+       void SetdcaCut(Double_t DCAcutr, Double_t DCAcutz) { fDCAcutr = DCAcutr; fDCAcutz = DCAcutz;};
+       
        void SetUseEMCal() { fUseEMCal=kTRUE;};
        void SetUseTrigger() { fUseTrigger=kTRUE;};
        void SetUseShowerShapeCut(Bool_t UseShowerShapeCut=kFALSE) { fUseShowerShapeCut=UseShowerShapeCut;};
@@ -353,6 +357,9 @@ private:
        Double_t                        fAngleCut;
        Double_t                        fChi2Cut;
        Double_t                        fDCAcut;
+       Double_t                        fDCAcutr;
+       Double_t                        fDCAcutz;
+       
        Bool_t                          fMassCutFlag;
        Bool_t                          fAngleCutFlag;
        Bool_t                          fChi2CutFlag;

diff --git a/PWGHF/hfe/macros/configs/pPb/ConfigEMCalHFEpA.C b/PWGHF/hfe/macros/configs/pPb/ConfigEMCalHFEpA.C
index 2226c829f46b51ea7014004ce541db7cfdd1ad37..4bb6ec9fb683e69ee873fb6f81e854aae08b6967 100644 (file)
--- a/PWGHF/hfe/macros/configs/pPb/ConfigEMCalHFEpA.C
+++ b/PWGHF/hfe/macros/configs/pPb/ConfigEMCalHFEpA.C
@@ -64,8 +64,8 @@ Int_t EMCalThreshould = 0 //0 == EG1, 1 == EG2
        //testing this line for the DCA cut
        //hfecuts->SetRequireDCAToVertex();
        //DCA cut included in the analysis 12 March 2014
-       if(configIndex==18) hfecuts->SetMaxImpactParam(2,5); 
-       else if(configIndex==19) hfecuts->SetMaxImpactParam(0.5,1);
+       if(configIndex==18) hfecuts->SetMaxImpactParam(2,3); //changed z to 3
+       else if(configIndex==19) hfecuts->SetMaxImpactParam(0.5,1);//r,z   default on AOD 2.4, 3.2
        else if(configIndex==83) hfecuts->SetMaxImpactParam(0.1,0.2);
        else hfecuts->SetMaxImpactParam(1,2);                                                                   //DCA to vertex
        
@@ -96,6 +96,12 @@ Int_t EMCalThreshould = 0 //0 == EG1, 1 == EG2
                //Bool_t isTrigger = kFALSE;
        if(isTrigger) task->SetUseTrigger();
        
+       if(configIndex==104)  task->SetdcaCut(2,3);//r,z
+       else if(configIndex==105)  task->SetdcaCut(1.5,2.5);//r,z
+       else if(configIndex==106)  task->SetdcaCut(0.5,1);//r,z
+       else if(configIndex==107)  task->SetdcaCut(0.5,1.5);//r,z
+       else if(configIndex==108)  task->SetdcaCut(0.1,0.2);//r,z
+       //else task->SetdcaCut(1,2);//r,z
        
                
        if(configIndex==100){
@@ -145,8 +151,6 @@ Int_t EMCalThreshould = 0 //0 == EG1, 1 == EG2
        else if (configIndex==44) task->SetEtaCut(-0.5,0.5);
        else if (configIndex==45) task->SetEtaCut(-0.4,0.4);
        else if (configIndex==46) task->SetEtaCut(-0.3,0.3);
-               //else if (configIndex==47) task->SetEtaCut(-0.4,0.4);
-               //else if (configIndex==48) task->SetEtaCut(-0.3,0.3);
        else task->SetEtaCut(-0.6,0.6);
        
        //track matching cuts
@@ -157,21 +161,27 @@ Int_t EMCalThreshould = 0 //0 == EG1, 1 == EG2
 
        //E/p Cuts
        if (configIndex==60) task->SetEoverPCut(0.85,1.2);
-       else if (configIndex==61) task->SetEoverPCut(0.75,1.25);        
-       else if (configIndex==88) task->SetEoverPCut(0.76,1.2);
-       else if (configIndex==89) task->SetEoverPCut(0.78,1.2);
-       else if (configIndex==90) task->SetEoverPCut(0.79,1.2);
-       else if (configIndex==90) task->SetEoverPCut(0.81,1.2);
-       else if (configIndex==90) task->SetEoverPCut(0.82,1.2);
-       else if (configIndex==91) task->SetEoverPCut(0.84,1.2);
-       else if (configIndex==92) task->SetEoverPCut(0.86,1.2);
-       else if (configIndex==93) task->SetEoverPCut(0.85,1.3);
-       else if (configIndex==62) task->SetEoverPCut(0.70,1.2);
-       else if (configIndex==63) task->SetEoverPCut(0.80,1.25);
-       else if (configIndex==64) task->SetEoverPCut(0.9,1.3);
-       else if (configIndex==65) task->SetEoverPCut(0.95,1.3);
-       else if (configIndex==66) task->SetEoverPCut(0.75,1.2);
-       else task->SetEoverPCut(0.8,1.2);
+       else if (configIndex==61) task->SetEoverPCut(0.70,1.2);
+       else if (configIndex==62) task->SetEoverPCut(0.75,1.2);
+       else if (configIndex==63) task->SetEoverPCut(0.76,1.2);
+       else if (configIndex==64) task->SetEoverPCut(0.77,1.2);
+       else if (configIndex==65) task->SetEoverPCut(0.78,1.2);
+       else if (configIndex==66) task->SetEoverPCut(0.79,1.2);
+       
+       else if (configIndex==67) task->SetEoverPCut(0.81,1.2);
+       else if (configIndex==68) task->SetEoverPCut(0.82,1.2);
+       else if (configIndex==69) task->SetEoverPCut(0.83,1.2);
+       
+       else if (configIndex==88) task->SetEoverPCut(0.84,1.2);
+       else if (configIndex==89) task->SetEoverPCut(0.86,1.2);
+       
+       else if (configIndex==90) task->SetEoverPCut(0.85,1.3);
+       else if (configIndex==91) task->SetEoverPCut(0.80,1.3);
+       else if (configIndex==92) task->SetEoverPCut(0.81,1.3);
+       else if (configIndex==93) task->SetEoverPCut(0.78,1.3);
+       
+       
+       else task->SetEoverPCut(0.80,1.2);