updated for nonHFE reconstruction study

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
index cd1ae4d67bb003d12917c05a622da9309762d09a..11841dacd99f8181fcff0b75d99ad9b634c42091 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
@@ -127,6 +127,7 @@ ClassImp(AliAnalysisTaskFlowTPCEMCalQCSP)
 ,fCentralityMin(0)
 ,fCentralityMax(0)
 ,fInvmassCut(0)
+,fpTCut(0)
 ,fTrigger(0)
 ,fPhi(0)
 ,fEta(0)
@@ -210,6 +211,7 @@ AliAnalysisTaskFlowTPCEMCalQCSP::AliAnalysisTaskFlowTPCEMCalQCSP()
 ,fCentralityMin(0)
 ,fCentralityMax(0)
 ,fInvmassCut(0)
+,fpTCut(0)
 ,fTrigger(0)
 ,fPhi(0)
 ,fEta(0)
@@ -476,7 +478,7 @@ if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInput
    if(TMath::Abs(cluster->GetTrackDx()) > 0.05 || TMath::Abs(cluster->GetTrackDz()) > 0.05) continue;
 
    pt = track->Pt();         //pt track after cuts
-   if(pt<1) continue;
+   if(pt<fpTCut) continue;
    fClsE = cluster->E();
    p = track->P();
   // dEdx = track->GetTPCsignal();

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h
index 0801560ed098cd95a2b348992fd05375b16fe928..bc005b0620ce8eb975ee8a379a414f7599c5fcbc 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h
@@ -1,5 +1,5 @@
-#ifndef AlIANALYSISTASKFLOWTPCEMCALQCSP_H
-#define AlIANALYSISTASKFLOWTPCEMCALQCSP_H
+#ifndef ALIANALYSISTASKFLOWTPCEMCALQCSP_H
+#define ALIANALYSISTASKFLOWTPCEMCALQCSP_H
 
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
@@ -51,6 +51,7 @@ class AliAnalysisTaskFlowTPCEMCalQCSP : public AliAnalysisTaskSE {
     void                                 CheckCentrality(AliAODEvent *event,Bool_t &centralitypass); //to use only events with the correct centrality....
     void                                 SelectPhotonicElectron(Int_t itrack,const AliAODTrack *track, Bool_t &fFlagPhotonicElec);
     void                                 SetInvariantMassCut(Double_t invmass) {fInvmassCut = invmass;};
+    void                                 SetpTCuttrack(Double_t ptcut) {fpTCut = ptcut;};
     void                                 SetTrigger(Int_t trig) {fTrigger = trig;};
     void                                 SetFlowSideBands(Bool_t sidebandsflow){fSideBandsFlow = sidebandsflow;}
     void                                 Setphiminuspsi(Bool_t phipsi){fPhiminusPsi = phipsi;}
@@ -91,6 +92,7 @@ class AliAnalysisTaskFlowTPCEMCalQCSP : public AliAnalysisTaskSE {
     Double_t             fCentralityMin; // lower bound of cenrality bin
     Double_t             fCentralityMax; // upper bound of centrality bin
     Double_t            fInvmassCut;           //invariant mass cut value
+    Double_t             fpTCut;    //pt cut value
     Int_t                   fTrigger;          //invariant mass cut value
     TH1F                 *fPhi; //! QA plot of azimuthal distribution of tracks used for event plane estimation
     TH1F                 *fEta; //! QA plot of eta distribution of tracks used for event plane estimation

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
index d3192c46d2b8fa447ac8b940f991277b1dc8e892..2d9fecef0c6607fada2d1fb571380f54acbc6b19 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
@@ -123,6 +123,7 @@ ClassImp(AliAnalysisTaskFlowTPCTOFQCSP)
 ,fCentralityMin(0)
 ,fCentralityMax(0)
 ,fInvmassCut(0)
+,fpTCut(0)
 ,fTrigger(0)
 ,fPhi(0)
 ,fEta(0)
@@ -202,6 +203,7 @@ AliAnalysisTaskFlowTPCTOFQCSP::AliAnalysisTaskFlowTPCTOFQCSP()
 ,fCentralityMin(0)
 ,fCentralityMax(0)
 ,fInvmassCut(0)
+,fpTCut(0)
 ,fTrigger(0)
 ,fPhi(0)
 ,fEta(0)
@@ -455,7 +457,7 @@ if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInput
      Double_t eta = track->Eta();
      Double_t phi = track->Phi();
      Double_t pt = track->Pt();         //pt track after cuts
-     if(pt<1.) continue;
+     if(pt<fpTCut) continue;
 //==========================================================================================================
 //=========================================PID==============================================================
      if(track->GetTPCsignalN() < fTPCS) continue;
@@ -614,9 +616,7 @@ void AliAnalysisTaskFlowTPCTOFQCSP::SelectPhotonicElectron(Int_t itrack,const Al
     openingAngle = ge1.GetAngle(ge2);
     if(fFlagLS) fOpeningAngleLS->Fill(openingAngle);
     if(fFlagULS) fOpeningAngleULS->Fill(openingAngle);
-      
     if(openingAngle > fOpeningAngleCut) continue;
-
       }
       
       

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
index b7ae661febdf43d3252d1313c59c44b160978f65..8ca52c49e739d1e64597438987df05767387f41c 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
@@ -48,6 +48,7 @@ class AliAnalysisTaskFlowTPCTOFQCSP : public AliAnalysisTaskSE {
     void                                 CheckCentrality(AliAODEvent *event,Bool_t &centralitypass); //to use only events with the correct centrality....
     void                                 SelectPhotonicElectron(Int_t itrack,const AliAODTrack *track, Bool_t &fFlagPhotonicElec);
     void                                 SetInvariantMassCut(Double_t invmass) {fInvmassCut = invmass;};
+    void                                 SetpTCuttrack(Double_t ptcut) {fpTCut = ptcut;};
     void                                 SetTrigger(Int_t trig) {fTrigger = trig;};
     void                                 SetTPCS(Int_t sig) {fTPCS = sig;};
     void                                 SetVz(Int_t ver) {fVz = ver;};
@@ -90,6 +91,7 @@ class AliAnalysisTaskFlowTPCTOFQCSP : public AliAnalysisTaskSE {
     Double_t             fCentralityMin; // lower bound of cenrality bin
     Double_t             fCentralityMax; // upper bound of centrality bin
     Double_t            fInvmassCut;           //invariant mass cut value
+    Double_t             fpTCut;    //pt cut value
     Int_t                   fTrigger;          //invariant mass cut value
     TH1F                 *fPhi; //! QA plot of azimuthal distribution of tracks used for event plane estimation
     TH1F                 *fEta; //! QA plot of eta distribution of tracks used for event plane estimation

diff --git a/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C b/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
index 0e8ab75bef824df43371586f556a40bd08af00cb..fd941390611c13bef86781e20554573f37647a37 100644 (file)
--- a/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
+++ b/PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
@@ -17,6 +17,7 @@ AliAnalysisTaskFlowTPCEMCalQCSP*  AddTaskFlowTPCEMCalQCSP(
                                               Float_t centrMin ,\r
                                               Float_t centrMax ,\r
                                               Double_t InvmassCut,\r
+                                              Double_t pTCut,\r
                                               Int_t Trigger,\r
                                               Double_t minTPC,\r
                                               Double_t maxTPC,\r
@@ -78,7 +79,7 @@ AliAnalysisTaskFlowTPCEMCalQCSP*  AddTaskFlowTPCEMCalQCSP(
   taskHFE->SetFlowSideBands(SideBandsFlow);\r
   taskHFE->Setphiminuspsi(Phi_minus_psi);\r
   taskHFE->SetPurity(purity);\r
-    \r
+  taskHFE->SetpTCuttrack(pTCut);\r
     \r
 //set RP cuts for flow package analysis\r
   cutsRP = new AliFlowTrackCuts(Form("RFPcuts%s",uniqueID));\r

diff --git a/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C b/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
index cda526070b8653314145c9e9393a75152ee0599e..fbc6f85aa8820c3948b2e1b532aa2f71f5ab5219 100644 (file)
--- a/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
+++ b/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
@@ -17,7 +17,8 @@ AliAnalysisTaskFlowTPCTOFQCSP* AddTaskFlowTPCTOFQCSP(
                                               Float_t centrMin ,\r
                                               Float_t centrMax ,\r
                                               Double_t InvmassCut,\r
-                                              Bool_t op_ang = kTRUE,\r
+                                              Double_t pTCut,\r
+                                              Bool_t op_ang = kFALSE,\r
                                               Double_t op_angle_cut,\r
                                               Int_t Trigger,\r
                                               Double_t minTPCnsigma,\r
@@ -79,7 +80,8 @@ AliAnalysisTaskFlowTPCTOFQCSP* AddTaskFlowTPCTOFQCSP(
   taskHFE->SetPhiCut(PhiCut);\r
   taskHFE->SetOpeningAngleflag(op_ang);\r
   taskHFE->SetOpeningAngleCut(op_angle_cut);\r
-    \r
+  taskHFE->SetpTCuttrack(pTCut);\r
+\r
     \r
 //set RP cuts for flow package analysis\r
   cutsRP = new AliFlowTrackCuts(Form("RFPcuts%s",uniqueID));\r