]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGCF/FLOW/macros/AddTaskPidVZEROSP.C
Fix the CRLF new line to unix type via dos2unix
[u/mrichter/AliRoot.git] / PWGCF / FLOW / macros / AddTaskPidVZEROSP.C
index 5699e98bc4b9403a4ffadba054bc71f9e3888826..5136ed129d8e0e47398e1cc3ac854229178cf0d8 100644 (file)
@@ -1,77 +1,77 @@
-AddTaskPidVZEROSP(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){\r
-  gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C");\r
-\r
-  const Int_t ncentr = 9;\r
-  Int_t cmin[ncentr]={0,5,10,20,30,40,50,60,70};\r
-  Int_t cmax[ncentr]={5,10,20,30,40,50,60,70,80};\r
-\r
-  for(Int_t i=0;i < ncentr;i++){\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); // no pid\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-  }\r
-}\r
-\r
-createSPres(){\r
-  gSystem->Load("libVMC.so");\r
-  gSystem->Load("libPhysics.so");\r
-  gSystem->Load("libTree.so");\r
-  gSystem->Load("libMinuit.so");\r
-  gSystem->Load("libSTEERBase.so");\r
-  gSystem->Load("libANALYSIS.so");\r
-  gSystem->Load("libAOD.so");\r
-  gSystem->Load("libESD.so");\r
-  gSystem->Load("libANALYSIS.so");\r
-  gSystem->Load("libANALYSISalice.so");\r
-  gSystem->Load("libCORRFW.so");\r
-  gSystem->Load("libNetx.so");\r
-  gSystem->Load("libPWGflowBase.so");\r
-  \r
-  char name[200];\r
-  char *spe[4]={"pion","kaon","antipr","proton"};\r
-  char *tech[2]={"TOF","TPC"};\r
-  Int_t cmin[9]={0,5,10,20,30,40,50,60,70};\r
-  Int_t cmax[9]={5,10,20,30,40,50,60,70,80};\r
-  \r
-  TFile *f = new TFile("AnalysisResults.root");\r
-  TFile *fo = new TFile("results.root","RECREATE");\r
-  TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputSPanalysisTPCstandalone"));\r
-  TList* listTemp = directory->GetListOfKeys();\r
-  for(Int_t i=0;i < 9;i++){\r
-    TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->GetName()));\r
-    AliFlowAnalysisWithScalarProduct* sp2 = new AliFlowAnalysisWithScalarProduct();\r
-    sp2->GetOutputHistograms(list2);\r
-    sp2->Finish();\r
-    AliFlowCommonHistResults* res2=sp2->GetCommonHistsRes();\r
-    TH1D *h2=res2->GetHistDiffFlowPtPOI();\r
-    sprintf(name,"v2SP_%s_%i_%i","AllCharged",cmin[i],cmax[i]);\r
-    h2->SetName(name);\r
-    fo->cd();\r
-    h2->Write();\r
-    for(Int_t j=0;j < 2;j++){\r
-      for(Int_t k=0;k < 4;k++){\r
-       TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9+j*4+k+1)->GetName()));\r
-       AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();\r
-       sp->GetOutputHistograms(list);\r
-       sp->Finish();\r
-       AliFlowCommonHistResults* res=sp->GetCommonHistsRes();\r
-       TH1D *h=res->GetHistDiffFlowPtPOI();\r
-       sprintf(name,"v2SP_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]);\r
-       h->SetName(name);\r
-       fo->cd();\r
-       h->Write();\r
-      }\r
-    }\r
-  }\r
-  fo->Close();\r
-}\r
-\r
-\r
-\r
-\r
+AddTaskPidVZEROSP(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){
+  gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C");
+
+  const Int_t ncentr = 9;
+  Int_t cmin[ncentr]={0,5,10,20,30,40,50,60,70};
+  Int_t cmax[ncentr]={5,10,20,30,40,50,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); // no pid
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
+  }
+}
+
+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"};
+  Int_t cmin[9]={0,5,10,20,30,40,50,60,70};
+  Int_t cmax[9]={5,10,20,30,40,50,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 < 9;i++){
+    TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->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*9+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();
+}
+
+
+
+