updated for nonHFE reconstruction study
authorssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 3 Jun 2013 08:27:59 +0000 (08:27 +0000)
committerssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 3 Jun 2013 08:27:59 +0000 (08:27 +0000)
PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.cxx
PWGHF/hfe/AliAnalysisTaskFlowTPCEMCalQCSP.h
PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
PWGHF/hfe/macros/AddTaskFlowTPCEMCalQCSP.C
PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C

index cd1ae4d..11841da 100644 (file)
@@ -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();
index 0801560..bc005b0 100644 (file)
@@ -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
index d3192c4..2d9fece 100644 (file)
@@ -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;
-
       }
       
       
index b7ae661..8ca52c4 100644 (file)
@@ -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
index 0e8ab75..fd94139 100644 (file)
@@ -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
index cda5260..fbc6f85 100644 (file)
@@ -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