updated for invariant mass study of TPC+TOF

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
index c8b7548c6cd631c471483691a3837ac5f7a5648e..d3192c46d2b8fa447ac8b940f991277b1dc8e892 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.cxx
@@ -164,6 +164,10 @@ ClassImp(AliAnalysisTaskFlowTPCTOFQCSP)
 ,fTPCS(0)
 ,fVz(0)
 ,fPhiCut(0)
+,fOpeningAngleCut(0.1)
+,fOP_angle(0)
+,fOpeningAngleLS(0)
+,fOpeningAngleULS(0)
 {
   //Named constructor
 
@@ -239,6 +243,10 @@ AliAnalysisTaskFlowTPCTOFQCSP::AliAnalysisTaskFlowTPCTOFQCSP()
 ,fTPCS(0)
 ,fVz(0)
 ,fPhiCut(0)
+,fOpeningAngleCut(0.1)
+,fOP_angle(0)
+,fOpeningAngleLS(0)
+,fOpeningAngleULS(0)
 {
   //Default constructor
 //  fPID = new AliHFEpid("hfePid");
@@ -571,6 +579,7 @@ void AliAnalysisTaskFlowTPCTOFQCSP::SelectPhotonicElectron(Int_t itrack,const Al
     if(jTracks == itrack) continue;
     Double_t ptAsso=-999., nsigma=-999.0;
     Double_t mass=-999., width = -999;
+      Double_t openingAngle = -999.;
     Bool_t fFlagLS=kFALSE, fFlagULS=kFALSE;
 
   
@@ -601,6 +610,16 @@ void AliAnalysisTaskFlowTPCTOFQCSP::SelectPhotonicElectron(Int_t itrack,const Al
     if(TMath::Sqrt(TMath::Abs(chi2recg))>3.) continue;
     recg.GetMass(mass,width);
 
+      if(fOP_angle){
+    openingAngle = ge1.GetAngle(ge2);
+    if(fFlagLS) fOpeningAngleLS->Fill(openingAngle);
+    if(fFlagULS) fOpeningAngleULS->Fill(openingAngle);
+      
+    if(openingAngle > fOpeningAngleCut) continue;
+
+      }
+      
+      
     if(fFlagLS) fInvmassLS1->Fill(mass);
     if(fFlagULS) fInvmassULS1->Fill(mass);
            
@@ -764,6 +783,14 @@ void AliAnalysisTaskFlowTPCTOFQCSP::UserCreateOutputObjects()
   fMultvsCentr = new TH2F("fMultvsCentr", "Multiplicity vs centrality; centrality; Multiplicity", 100, 0., 100, 100, 0, 3000);
   fOutputList->Add(fMultvsCentr);
     
+  fOpeningAngleLS = new TH1F("fOpeningAngleLS","Opening angle for LS pairs",100,0,1);
+  fOutputList->Add(fOpeningAngleLS);
+    
+  fOpeningAngleULS = new TH1F("fOpeningAngleULS","Opening angle for ULS pairs",100,0,1);
+  fOutputList->Add(fOpeningAngleULS);
+  
+    
+    
 //----------------------------------------------------------------------------
     EPVzA = new TH1D("EPVzA", "EPVzA", 80, -2, 2);
     fOutputList->Add(EPVzA);

diff --git a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
index f6200e9d329525531464e832d9c15ee702f35546..b7ae661febdf43d3252d1313c59c44b160978f65 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
+++ b/PWGHF/hfe/AliAnalysisTaskFlowTPCTOFQCSP.h
@@ -53,6 +53,8 @@ class AliAnalysisTaskFlowTPCTOFQCSP : public AliAnalysisTaskSE {
     void                                 SetVz(Int_t ver) {fVz = ver;};
     void                                 SetQAPIDSparse(Bool_t qapidsparse) {fQAPIDSparse = qapidsparse;};
     void                                 SetPhiCut(Bool_t phicut){fPhiCut = phicut;};
+    void                                 SetOpeningAngleflag(Bool_t opang){fOP_angle = opang;};
+    void                                 SetOpeningAngleCut(Double_t opanglecut) {fOpeningAngleCut = opanglecut;};
     template <typename T> void           PlotVZeroMultiplcities(const T* event) const;
     template <typename T> void           SetNullCuts(T* aod);
     void                                 PrepareFlowEvent(Int_t iMulti, AliFlowEvent *FlowEv) const;
@@ -129,7 +131,11 @@ class AliAnalysisTaskFlowTPCTOFQCSP : public AliAnalysisTaskSE {
     Int_t                 fTPCS; //tpc signal cluster
     Int_t                 fVz; //vertex range
     Bool_t                           fPhiCut;       //Phi cut to simulate emcal acc
-
+    Double_t                 fOpeningAngleCut; //openingAngle cut value
+    Bool_t                fOP_angle; //to shitch on and off the op_angle cut
+    TH1F                            *fOpeningAngleLS;  //opening angle for LS pairs
+    TH1F                            *fOpeningAngleULS; //opening angle for ULS pairs
+    
     AliAnalysisTaskFlowTPCTOFQCSP(const AliAnalysisTaskFlowTPCTOFQCSP&); // not implemented
     AliAnalysisTaskFlowTPCTOFQCSP& operator=(const AliAnalysisTaskFlowTPCTOFQCSP&); // not implemented
 

diff --git a/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C b/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
index 25ec79fb077bdbc0998787816e4b92110a6f1139..cda526070b8653314145c9e9393a75152ee0599e 100644 (file)
--- a/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
+++ b/PWGHF/hfe/macros/AddTaskFlowTPCTOFQCSP.C
@@ -17,6 +17,8 @@ AliAnalysisTaskFlowTPCTOFQCSP* AddTaskFlowTPCTOFQCSP(
                                               Float_t centrMin ,\r
                                               Float_t centrMax ,\r
                                               Double_t InvmassCut,\r
+                                              Bool_t op_ang = kTRUE,\r
+                                              Double_t op_angle_cut,\r
                                               Int_t Trigger,\r
                                               Double_t minTPCnsigma,\r
                                               Double_t maxTPCnsigma,\r
@@ -75,6 +77,8 @@ AliAnalysisTaskFlowTPCTOFQCSP* AddTaskFlowTPCTOFQCSP(
   taskHFE->SetIDCuts(minTPCnsigma, maxTPCnsigma, minTOFnSigma, maxTOFnSigma);\r
   taskHFE->SetQAPIDSparse(QaPidSparse);\r
   taskHFE->SetPhiCut(PhiCut);\r
+  taskHFE->SetOpeningAngleflag(op_ang);\r
+  taskHFE->SetOpeningAngleCut(op_angle_cut);\r
     \r
     \r
 //set RP cuts for flow package analysis\r