add one addtask for VZERO SP pid hadrons, other centralites

diff --git a/PWGCF/FLOW/macros/AddTaskPidVZEROSP_2.C b/PWGCF/FLOW/macros/AddTaskPidVZEROSP_2.C
new file mode 100644 (file)
index 0000000..ee478bf
--- /dev/null
+++ b/PWGCF/FLOW/macros/AddTaskPidVZEROSP_2.C
@@ -0,0 +1,78 @@
+AddTaskPidVZEROSP_2(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8,Int_t side=0){
+  gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C");
+
+  const Int_t ncentr = 4;
+  Int_t cmin[ncentr]={30,50,60,70};
+  Int_t cmax[ncentr]={40,60,70,80};
+
+  for(Int_t i=0;i < ncentr;i++){
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side); // no pid
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax,"",side);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side);
+  }
+}
+
+createSPres(){
+  gSystem->Load("libVMC.so");
+  gSystem->Load("libPhysics.so");
+  gSystem->Load("libTree.so");
+  gSystem->Load("libMinuit.so");
+  gSystem->Load("libSTEERBase.so");
+  gSystem->Load("libANALYSIS.so");
+  gSystem->Load("libAOD.so");
+  gSystem->Load("libESD.so");
+  gSystem->Load("libANALYSIS.so");
+  gSystem->Load("libANALYSISalice.so");
+  gSystem->Load("libCORRFW.so");
+  gSystem->Load("libNetx.so");
+  gSystem->Load("libPWGflowBase.so");
+  
+  char name[200];
+  char *spe[4]={"pion","kaon","antipr","proton"};
+  char *tech[2]={"TOF","TPC"};
+  const Int_t ncentr = 4;
+  Int_t cmin[ncentr]={30,50,60,70};
+  Int_t cmax[ncentr]={40,60,70,80};
+  
+  TFile *f = new TFile("AnalysisResults.root");
+  TFile *fo = new TFile("results.root","RECREATE");
+  TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputSPanalysisTPCstandalone"));
+  TList* listTemp = directory->GetListOfKeys();
+  for(Int_t i=0;i < ncentr;i++){
+    TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr)->GetName()));
+    AliFlowAnalysisWithScalarProduct* sp2 = new AliFlowAnalysisWithScalarProduct();
+    sp2->GetOutputHistograms(list2);
+    sp2->Finish();
+    AliFlowCommonHistResults* res2=sp2->GetCommonHistsRes();
+    TH1D *h2=res2->GetHistDiffFlowPtPOI();
+    sprintf(name,"v2SP_%s_%i_%i","AllCharged",cmin[i],cmax[i]);
+    h2->SetName(name);
+    fo->cd();
+    h2->Write();
+    for(Int_t j=0;j < 2;j++){
+      for(Int_t k=0;k < 4;k++){
+       TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr+j*4+k+1)->GetName()));
+       AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();
+       sp->GetOutputHistograms(list);
+       sp->Finish();
+       AliFlowCommonHistResults* res=sp->GetCommonHistsRes();
+       TH1D *h=res->GetHistDiffFlowPtPOI();
+       sprintf(name,"v2SP_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]);
+       h->SetName(name);
+       fo->cd();
+       h->Write();
+      }
+    }
+  }
+  fo->Close();
+}
+
+
+
+