Zhong-Bao Yin: Updated Cascade Flow task

diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.cxx b/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.cxx
index e68faad5bf85c4b32cbf4ba4536e3d082197c547..4f5c8d52b094cd67b79718f2e89be673136e94a9 100644 (file)
--- a/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.cxx
+++ b/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.cxx
@@ -15,9 +15,9 @@
 \r
 /////////////////////////////////////////////////////\r
 // AliAnalysisTaskFlowCascade:\r
-// Analysis task to select K0/Lambda candidates for flow analysis.\r
+// Analysis task to select Xi and Omega candidates for flow analysis.\r
 //\r
-// Author: \r
+// Author: Zhong-Bao.Yin@cern.ch\r
 //////////////////////////////////////////////////////\r
 \r
 #include "TChain.h"\r
@@ -151,7 +151,7 @@ void AliAnalysisTaskFlowCascade::UserCreateOutputObjects()
   fQAList->SetOwner();\r
   AddQAEvents();\r
   AddQACandidates();\r
-  PostData(3,fQAList);\r
+  //  PostData(3,fQAList);\r
 \r
   AliFlowCommonConstants* cc = AliFlowCommonConstants::GetMaster();\r
   cc->SetNbinsMult(1);\r
@@ -470,9 +470,9 @@ void AliAnalysisTaskFlowCascade::UserExec(Option_t *)
   ((TH1D*)((TList*)fQAList->FindObject("Events"))->FindObject("RFPVZE"))\r
     ->Fill( mult );\r
 \r
-  if(fDebug) printf("TPCevent %d | VZEevent %d\n",\r
-                    fFlowEventTPC->NumberOfTracks(),\r
-                    fFlowEventVZE->NumberOfTracks() );\r
+  //  if(fDebug) printf("TPCevent %d | VZEevent %d\n",\r
+  //                  fFlowEventTPC->NumberOfTracks(),\r
+  //                  fFlowEventVZE->NumberOfTracks() );\r
   AddCandidates();\r
 \r
   PostData(1,fFlowEventTPC);\r
@@ -485,31 +485,32 @@ void AliAnalysisTaskFlowCascade::UserExec(Option_t *)
 //_____________________________________________________________________________\r
 void AliAnalysisTaskFlowCascade::AddCandidates(){\r
   \r
-  if(fDebug) printf("I received %d candidates\n",\r
-                   fCandidates->GetEntriesFast());\r
+  //  if(fDebug) printf("I received %d candidates\n",\r
+  //               fCandidates->GetEntriesFast());\r
   for(int iCand=0; iCand!=fCandidates->GetEntriesFast(); ++iCand ) {\r
     AliFlowCandidateTrack *cand \r
       = dynamic_cast<AliFlowCandidateTrack*>(fCandidates->At(iCand));\r
     if(!cand) continue;\r
-    if(fDebug) \r
-      printf(" >Checking at candidate %d with %d daughters: mass %f\n",\r
-            iCand, cand->GetNDaughters(), cand->Mass());\r
+    // if(fDebug) \r
+    //  printf(" >Checking at candidate %d with %d daughters: mass %f\n",\r
+    //      iCand, cand->GetNDaughters(), cand->Mass());\r
+    \r
     // untagging ===>                      \r
     for(int iDau=0; iDau != cand->GetNDaughters(); ++iDau) {\r
-      if(fDebug) \r
-       printf(" >Daughter %d with fID %d", iDau, cand->GetIDDaughter(iDau));\r
+      // if(fDebug) \r
+      // printf(" >Daughter %d with fID %d", iDau, cand->GetIDDaughter(iDau));\r
       for(int iRPs=0; iRPs != fFlowEventTPC->NumberOfTracks(); ++iRPs ) {\r
         AliFlowTrack *iRP \r
          = dynamic_cast<AliFlowTrack*>(fFlowEventTPC->GetTrack( iRPs ));\r
         if (!iRP) continue;\r
         if( !iRP->InRPSelection() ) continue;\r
         if( cand->GetIDDaughter(iDau) == iRP->GetID() ) {\r
-          if(fDebug) printf(" was in RP set");\r
+          //if(fDebug) printf(" was in RP set");\r
           iRP->SetForRPSelection(kFALSE);\r
           fFlowEventTPC->SetNumberOfRPs( fFlowEventTPC->GetNumberOfRPs() -1 );\r
         }\r
       }\r
-      if(fDebug) printf("\n");\r
+      //if(fDebug) printf("\n");\r
     }\r
     // <=== untagging\r
     cand->SetForPOISelection(kTRUE);\r
@@ -517,9 +518,9 @@ void AliAnalysisTaskFlowCascade::AddCandidates(){
     fFlowEventVZE->InsertTrack( ((AliFlowTrack*) cand) );\r
   }\r
 \r
-  if(fDebug) printf("TPCevent %d | VZEevent %d\n",\r
-                    fFlowEventTPC->NumberOfTracks(),\r
-                    fFlowEventVZE->NumberOfTracks() );\r
+  //  if(fDebug) printf("TPCevent %d | VZEevent %d\n",\r
+  //                fFlowEventTPC->NumberOfTracks(),\r
+  //                fFlowEventVZE->NumberOfTracks() );\r
 \r
 }\r
 \r

diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.h b/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.h
index 6e22540469b093663fcecf1a139d96629e5e889b..cd24925805bdc64ecfe8e0cb4e65d9788d0893b9 100644 (file)
--- a/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.h
+++ b/PWG/FLOW/Tasks/AliAnalysisTaskFlowCascade.h
@@ -69,7 +69,7 @@ class AliAnalysisTaskFlowCascade : public AliAnalysisTaskSE {
                             AliFlowTrackCuts *cutsRPTPC,\r
                             AliFlowTrackCuts *cutsRPVZE,\r
                             /* AliESDtrackCuts */AliFlowTrackCuts *cutsDau);\r
-  void SetDebug() {fDebug = true;}\r
+  //void SetDebug() {fDebug = true;}\r
   void SetSpecie(int specie){fSpecie = specie;}\r
   void SetCommonConstants(Int_t massBins, Double_t minMass, Double_t maxMass);\r
   void SetCuts2010(int setOfCuts);\r
@@ -79,7 +79,7 @@ class AliAnalysisTaskFlowCascade : public AliAnalysisTaskSE {
   virtual void Terminate(Option_t *);\r
   virtual void NotifyRun();\r
   \r
-  ClassDef(AliAnalysisTaskFlowCascade, 1);\r
+  ClassDef(AliAnalysisTaskFlowCascade, 2);\r
 };\r
 \r
 #endif\r

diff --git a/PWG/FLOW/Tasks/AddTaskFlowCascade.C b/PWGCF/FLOW/macros/AddTaskFlowCascade.C
similarity index 79%
rename from PWG/FLOW/Tasks/AddTaskFlowCascade.C
rename to PWGCF/FLOW/macros/AddTaskFlowCascade.C
index 4eaa158e7864e586343c43a10240bdc4f6072e8e..9b6837e2588b6b08d0a161d9ef58101e69715471 100644 (file)
--- a/PWG/FLOW/Tasks/AddTaskFlowCascade.C
+++ b/PWGCF/FLOW/macros/AddTaskFlowCascade.C
@@ -59,8 +59,8 @@ void AddTaskFlowCascade(int trigger, float centrMin, float centrMax,
   //taskSel->SetDebug();\r
   taskSel->SetCuts2010(cuts);\r
   taskSel->SetSpecie(specie);\r
-  //printf( "CMM %d %f %f\n", MassBins(specie), MinMass(specie), MaxMass(specie) );\r
-  taskSel->SetCommonConstants( MassBins(specie), MinMass(specie), MaxMass(specie) );\r
+  //printf( "CMM %d %f %f\n", MassBins(specie), MinMass(specie), MSFT_MaxMass(specie) );\r
+  taskSel->SetCommonConstants( MSFT_MassBins(specie), MSFT_MinMass(specie), MSFT_MaxMass(specie) );\r
   AliAnalysisDataContainer *cOutHist \r
     = mgr->CreateContainer(Form("OutHistos_%s",suffixName.Data()),\r
                           TList::Class(),\r
@@ -83,47 +83,48 @@ void AddTaskFlowCascade(int trigger, float centrMin, float centrMax,
 \r
   //-------------------FLOW TASKS----------------------------\r
   AliFlowTrackSimpleCuts *filter[15], *filterhf[15][2]; // MASS BANDS\r
-  for(int mb=0; mb!=MassBands(0); ++mb) {\r
+  for(int mb=0; mb!=MSFT_MassBands(0); ++mb) {\r
     filter[mb] = new AliFlowTrackSimpleCuts( Form("Filter_MB%d",mb) );\r
     filter[mb]->SetEtaMin( -0.8 ); \r
     filter[mb]->SetEtaMax( +0.8 );\r
-    filter[mb]->SetMassMin( MassBandLowEdge(specie, mb) ); \r
-    filter[mb]->SetMassMax( MassBandLowEdge(specie, mb+1) );\r
+    filter[mb]->SetMassMin( MSFT_MassBandLowEdge(specie, mb) ); \r
+    filter[mb]->SetMassMax( MSFT_MassBandLowEdge(specie, mb+1) );\r
 \r
     filterhf[mb][0] = new AliFlowTrackSimpleCuts( Form("Filterhf0_MB%d",mb) );\r
     filterhf[mb][0]->SetEtaMin( 0.0 ); \r
     filterhf[mb][0]->SetEtaMax( +0.8 );\r
-    filterhf[mb][0]->SetMassMin( MassBandLowEdge(specie, mb) ); \r
-    filterhf[mb][0]->SetMassMax( MassBandLowEdge(specie, mb+1) );\r
+    filterhf[mb][0]->SetMassMin( MSFT_MassBandLowEdge(specie, mb) ); \r
+    filterhf[mb][0]->SetMassMax( MSFT_MassBandLowEdge(specie, mb+1) );\r
 \r
     filterhf[mb][1] = new AliFlowTrackSimpleCuts( Form("Filterhf1_MB%d",mb) );\r
     filterhf[mb][1]->SetEtaMin( -0.8 ); \r
     filterhf[mb][1]->SetEtaMax( 0.0 );\r
-    filterhf[mb][1]->SetMassMin( MassBandLowEdge(specie, mb) ); \r
-    filterhf[mb][1]->SetMassMax( MassBandLowEdge(specie, mb+1) );\r
+    filterhf[mb][1]->SetMassMin( MSFT_MassBandLowEdge(specie, mb) ); \r
+    filterhf[mb][1]->SetMassMax( MSFT_MassBandLowEdge(specie, mb+1) );\r
 \r
     if(method.Contains("QC"))\r
-      AddQCmethod( Form("QCTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
+      MSFT_AddQCmethod( Form("QCTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
                   exc_TPC, filter[mb]); // QC TPC\r
     if(method.Contains("TPCSP")){\r
-      AddSPmethod( Form("SPTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
+      MSFT_AddSPmethod( Form("SPTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
                   exc_TPC, filterhf[mb][0], "Qa" ); // SP TPC Qa\r
-      AddSPmethod( Form("SPTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
+      MSFT_AddSPmethod( Form("SPTPCMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
                 exc_TPC, filterhf[mb][1], "Qb" ); // SP TPC Qb\r
     }\r
     if(method.Contains("VZESP")){\r
-      AddSPmethod( Form("SPVZEMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
+      MSFT_AddSPmethod( Form("SPVZEMB%d", mb), folderName.Data(), suffixName.Data(), harmonic, \r
                   exc_VZE, filter[mb], "Qa" ); // SP VZE Qa\r
-      AddSPmethod( Form("SPVZEMB%d",mb), folderName.Data(), suffixName.Data(), harmonic, \r
+      MSFT_AddSPmethod( Form("SPVZEMB%d",mb), folderName.Data(), suffixName.Data(), harmonic, \r
                   exc_VZE, filter[mb], "Qb" ); // SP VZE Qa\r
-      AddSPmethod( Form("SPVZEMB%d",mb), folderName.Data(), suffixName.Data(), harmonic,\r
+      MSFT_AddSPmethod( Form("SPVZEMB%d",mb), folderName.Data(), suffixName.Data(), harmonic,\r
                   exc_VZE, filter[mb], "QaQb" ); // SP VZE QaQb\r
     }\r
   }\r
 }\r
 \r
-void AddQCmethod(char *name, TString myFolder, char *thecuts, int harmonic, \r
-                AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI=NULL) {\r
+void MSFT_AddQCmethod(char *name, TString myFolder, char *thecuts, int harmonic, \r
+                     AliAnalysisDataContainer *flowEvent, \r
+                     AliFlowTrackSimpleCuts *cutsPOI=NULL) {\r
   TString fileName = AliAnalysisManager::GetCommonFileName();\r
   myFolder.Append( Form("v%d",harmonic) );\r
   TString myName = Form("%sv%d_%s", name, harmonic, thecuts);\r
@@ -153,9 +154,9 @@ void AddQCmethod(char *name, TString myFolder, char *thecuts, int harmonic,
 }\r
 \r
 \r
-void AddSPmethod(char *name, TString myFolder, char *thecuts, int harmonic,\r
-                AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI=NULL,\r
-                 char *Qvector) {\r
+void MSFT_AddSPmethod(char *name, TString myFolder, char *thecuts, int harmonic,\r
+                     AliAnalysisDataContainer *flowEvent, AliFlowTrackSimpleCuts *cutsPOI=NULL,\r
+                     char *Qvector) {\r
   TString fileName = AliAnalysisManager::GetCommonFileName();\r
   myFolder.Append( Form("v%d",harmonic) );\r
   TString myNameSP = Form("%sv%d%s_%s", name, harmonic, Qvector, thecuts);\r
@@ -167,7 +168,7 @@ void AddSPmethod(char *name, TString myFolder, char *thecuts, int harmonic,
   AliAnalysisTaskFilterFE *tskFilter \r
     = new AliAnalysisTaskFilterFE( Form("TaskFilter_%s",myNameSP.Data()),\r
                                   NULL, cutsPOI);\r
-  tskFilter->SetSubeventEtaRange( -5.0, 0.0, 0.0, +5.0 );\r
+  tskFilter->SetSubeventEtaRange( -5.0, -1.0, 1.0, +5.0 );\r
   mgr->AddTask(tskFilter);\r
   mgr->ConnectInput( tskFilter, 0, flowEvent);\r
   mgr->ConnectOutput(tskFilter, 1, flowEvent2);\r
@@ -185,7 +186,7 @@ void AddSPmethod(char *name, TString myFolder, char *thecuts, int harmonic,
   mgr->ConnectOutput(tskSP,1,outSP);\r
 }\r
 \r
-double MassBandLowEdge( int nv0, int mb ) {\r
+double MSFT_MassBandLowEdge( int nv0, int mb ) {\r
   if(nv0==0) {\r
     double lowEdge[15+1]={ 1.282, 1.292, 1.302, 1.312, 1.317, 1.319, 1.321, 1.322, \r
                           1.323, 1.325, 1.327, 1.332, 1.342, 1.352, 1.362, 1.372};\r
@@ -196,7 +197,7 @@ double MassBandLowEdge( int nv0, int mb ) {
   return lowEdge[mb];\r
 }\r
 \r
-int MassBands( int nv0 ) {\r
+int MSFT_MassBands( int nv0 ) {\r
   if(nv0==0) {\r
     return 15;\r
   } else if(nv0==1) {\r
@@ -204,7 +205,7 @@ int MassBands( int nv0 ) {
   }\r
 }\r
 \r
-int MassBins( int nv0 ) {\r
+int MSFT_MassBins( int nv0 ) {\r
   if(nv0==0) {\r
     return 90;\r
   } else if(nv0==1) {\r
@@ -212,10 +213,10 @@ int MassBins( int nv0 ) {
   }\r
 }\r
 \r
-double MinMass( int nv0 ) {\r
-  return MassBandLowEdge( nv0, 0 );\r
+double MSFT_MinMass( int nv0 ) {\r
+  return MSFT_MassBandLowEdge( nv0, 0 );\r
 }\r
 \r
-double MaxMass( int nv0 ) {\r
-  return MassBandLowEdge( nv0, MassBands(nv0) );\r
+double MSFT_MaxMass( int nv0 ) {\r
+  return MSFT_MassBandLowEdge( nv0, MSFT_MassBands(nv0) );\r
 }\r