ITSTPCTOF flow analysis in the acceptance of EMCal

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.cxx
index c838d2a5686c19d95a0581d852fa6dcca7591cbf..2e6db7d04c1eb9e79c36d77d1312fc4a1aac2355 100755 (executable)
--- a/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.cxx
@@ -207,6 +207,7 @@ AliAnalysisTaskFlowITSTPCTOFQCSP::AliAnalysisTaskFlowITSTPCTOFQCSP(const char *n
 ,fSparseMassLS(0)
 ,fAssoTPCCluster(0)
 ,fAssoITSRefit(0)
+,fPhiCut(0)
 {
     //Named constructor
     
@@ -325,6 +326,7 @@ AliAnalysisTaskFlowITSTPCTOFQCSP::AliAnalysisTaskFlowITSTPCTOFQCSP()
 ,fSparseMassLS(0)
 ,fAssoTPCCluster(0)
 ,fAssoITSRefit(0)
+,fPhiCut(0)
 {
     //Default constructor
     fPID = new AliHFEpid("hfePid");
@@ -588,6 +590,13 @@ void AliAnalysisTaskFlowITSTPCTOFQCSP::UserExec(Option_t*)
         //==========================================================================================================
         Double_t eta = track->Eta();
         Double_t phi = track->Phi();
+        
+        if(fPhiCut){
+            if(phi<1.4 || phi >3.14)continue; //to have same EMCal phi acceptance
+        }
+
+        
+        
         Double_t pt = track->Pt();         //pt track after cuts
         if(pt<fpTCutmin || pt>fpTCutmax) continue;
         //==========================================================================================================
@@ -1314,7 +1323,7 @@ void AliAnalysisTaskFlowITSTPCTOFQCSP::CheckCentrality(AliAODEvent* event, Bool_
     }//...after centrality selectrion
     //============================================================================================================================
     if(fMultCut){
-        if(fTrigger==1 || fTrigger==4){
+        if(fTrigger==1 || fTrigger==4 || fTrigger==5){
             if(! (multTPC > (-36.73 + 1.48*multGlob) && multTPC < (62.87 + 1.78*multGlob))){
                 //   cout <<" Trigger ==" <<fTrigger<< endl;
                 centralitypass = kFALSE;

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.h b/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.h
index c84242966046dfb0f3c62e2381078658caf540b8..b39149e6bf3e5911dd073077d1ed5e4635f345c0 100755 (executable)
--- a/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.h
+++ b/PWGHF/hfe/AliAnalysisTaskFlowITSTPCTOFQCSP.h
@@ -6,7 +6,7 @@
 
 ////////////////////////////////////////////////////////////////////////
 //                                                                    //
-//  Task for Heavy Flavour Electron Flow  TPC plus TOF                //
+//  Task for Heavy Flavour Electron Flow ITS TPC TOF                  //
 //                                                                    //
 //  Author: Andrea Dubla (Utrecht University)                         //
 //                                                                    //
@@ -74,6 +74,7 @@ public:
     Bool_t                               IsEventSelectedForCentrFlattening(Float_t centvalue);
     void                                 SetAssoITSRefit(Bool_t itsref) {fAssoITSRefit = itsref;};
     void                                 SetAssoTPCCluster(Int_t tpc_clust) {fAssoTPCCluster = tpc_clust;};
+    void                                 SetPhiCut(Bool_t phicut){fPhiCut = phicut;};
 
     
     AliHFEpid *GetPID() const { return fPID; };
@@ -183,7 +184,8 @@ private:
     THnSparseF           *fSparseMassLS;//!ssssss
     Int_t                 fAssoTPCCluster;//asso tpc cluster
     Bool_t                fAssoITSRefit;//asso its refit
-    
+    Bool_t                fPhiCut;//Phi cut to simulate emcal acc
+
     
     AliAnalysisTaskFlowITSTPCTOFQCSP(const AliAnalysisTaskFlowITSTPCTOFQCSP&); // not implemented
     AliAnalysisTaskFlowITSTPCTOFQCSP& operator=(const AliAnalysisTaskFlowITSTPCTOFQCSP&); // not implemented

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
index 1a57e7a0b7899e80b37112c6dc68f8228635be45..307af57e7422ce5bebf37fc2c3d965a8f9a75eb1 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
@@ -1304,7 +1304,7 @@ void AliAnalysisTaskFlowTPCEMCalQCSP::CheckCentrality(AliAODEvent* event, Bool_t
     }//...after centrality selectrion
     //============================================================================================================================
     if(fMultCut){
-        if(fTrigger==1 || fTrigger==4){
+        if(fTrigger==1 || fTrigger==4 || fTrigger==5){
             if(! (multTPC > (-36.73 + 1.48*multGlob) && multTPC < (62.87 + 1.78*multGlob))){
                 //   cout <<" Trigger ==" <<fTrigger<< endl;
                 centralitypass = kFALSE;

diff --git a/PWGHF/hfe/macros/AddTaskFlowITSTPCTOFQCSP.C b/PWGHF/hfe/macros/AddTaskFlowITSTPCTOFQCSP.C
index f3019c9a1702ad0abafc522b0ceeeba1bd34ac45..5e95af78318899b827fd69e82ceb48395ec668ca 100755 (executable)
--- a/PWGHF/hfe/macros/AddTaskFlowITSTPCTOFQCSP.C
+++ b/PWGHF/hfe/macros/AddTaskFlowITSTPCTOFQCSP.C
@@ -34,9 +34,10 @@ AliAnalysisTaskFlowITSTPCTOFQCSP* AddTaskFlowITSTPCTOFQCSP(
                                                      Int_t minTPCCluster,
                                                      Int_t TPCS,
                                                      AliHFEextraCuts::ITSPixel_t pixel,
-                                                    Int_t TPCClusterforAsso = 80,
+                                                     Int_t TPCClusterforAsso = 80,
                                                      Bool_t AssoITSref = kTRUE,
-                                                     Double_t ptminassocut = 0.25,
+                                                     Double_t ptminassocut = 0.0,
+                                                     Bool_t PhiCut = kFALSE,
                                                      Bool_t PhotonicElectronDCA = kFALSE,
                                                     // Bool_t QaPidSparse = kFALSE,
                                                      const char *Cent = "V0M",
@@ -103,6 +104,11 @@ AliAnalysisTaskFlowITSTPCTOFQCSP* AddTaskFlowITSTPCTOFQCSP(
     taskHFE->SetPtMinAssoCut(ptminassocut);
     taskHFE->SetAssoTPCCluster(TPCClusterforAsso);
     taskHFE->SetAssoITSRefit(AssoITSref);
+    taskHFE->SetPhiCut(PhiCut);
+
+    
+    
+    
     //set RP cuts for flow package analysis
     cutsRP = new AliFlowTrackCuts(Form("RFPcuts%s",uniqueID));
     if(!cutsRP) {