]>
Commit | Line | Data |
---|---|---|
24e5dc5c | 1 | AddTaskPidVZEROSP(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8,Int_t side=0,Int_t filterbit=1,Bool_t TOFbeta=kFALSE){ |
24373b38 | 2 | gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C"); |
3 | ||
b671ca76 | 4 | const Int_t ncentr = 5; |
5 | Int_t cmin[ncentr]={0,5,10,20,40}; | |
6 | Int_t cmax[ncentr]={5,10,20,30,50}; | |
24373b38 | 7 | |
8 | for(Int_t i=0;i < ncentr;i++){ | |
24e5dc5c | 9 | if(!TOFbeta){ |
10 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit); // no pid | |
11 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit); | |
12 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit); | |
13 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax,"",side,filterbit); | |
14 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit); | |
15 | } | |
16 | else{ | |
17 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit); // no pid | |
18 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit); | |
19 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit); | |
20 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbeta,-1,2,0,etamin,etamax,"",side,filterbit); | |
21 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit); | |
22 | } | |
e9d352d8 | 23 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit); |
24 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit); | |
25 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax,"",side,filterbit); | |
26 | AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit); | |
24373b38 | 27 | } |
28 | } | |
29 | ||
30 | createSPres(){ | |
31 | gSystem->Load("libVMC.so"); | |
32 | gSystem->Load("libPhysics.so"); | |
33 | gSystem->Load("libTree.so"); | |
34 | gSystem->Load("libMinuit.so"); | |
35 | gSystem->Load("libSTEERBase.so"); | |
36 | gSystem->Load("libANALYSIS.so"); | |
37 | gSystem->Load("libAOD.so"); | |
38 | gSystem->Load("libESD.so"); | |
39 | gSystem->Load("libANALYSIS.so"); | |
40 | gSystem->Load("libANALYSISalice.so"); | |
41 | gSystem->Load("libCORRFW.so"); | |
42 | gSystem->Load("libNetx.so"); | |
43 | gSystem->Load("libPWGflowBase.so"); | |
44 | ||
45 | char name[200]; | |
46 | char *spe[4]={"pion","kaon","antipr","proton"}; | |
47 | char *tech[2]={"TOF","TPC"}; | |
b671ca76 | 48 | const Int_t ncentr = 5; |
49 | Int_t cmin[ncentr]={0,5,10,20,40}; | |
50 | Int_t cmax[ncentr]={5,10,20,30,50}; | |
24373b38 | 51 | |
52 | TFile *f = new TFile("AnalysisResults.root"); | |
53 | TFile *fo = new TFile("results.root","RECREATE"); | |
54 | TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputSPanalysisTPCstandalone")); | |
55 | TList* listTemp = directory->GetListOfKeys(); | |
b671ca76 | 56 | for(Int_t i=0;i < ncentr;i++){ |
57 | TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr)->GetName())); | |
24373b38 | 58 | AliFlowAnalysisWithScalarProduct* sp2 = new AliFlowAnalysisWithScalarProduct(); |
59 | sp2->GetOutputHistograms(list2); | |
60 | sp2->Finish(); | |
61 | AliFlowCommonHistResults* res2=sp2->GetCommonHistsRes(); | |
62 | TH1D *h2=res2->GetHistDiffFlowPtPOI(); | |
63 | sprintf(name,"v2SP_%s_%i_%i","AllCharged",cmin[i],cmax[i]); | |
64 | h2->SetName(name); | |
65 | fo->cd(); | |
66 | h2->Write(); | |
67 | for(Int_t j=0;j < 2;j++){ | |
68 | for(Int_t k=0;k < 4;k++){ | |
b671ca76 | 69 | TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr+j*4+k+1)->GetName())); |
24373b38 | 70 | AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct(); |
71 | sp->GetOutputHistograms(list); | |
72 | sp->Finish(); | |
73 | AliFlowCommonHistResults* res=sp->GetCommonHistsRes(); | |
74 | TH1D *h=res->GetHistDiffFlowPtPOI(); | |
75 | sprintf(name,"v2SP_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]); | |
76 | h->SetName(name); | |
77 | fo->cd(); | |
78 | h->Write(); | |
79 | } | |
80 | } | |
81 | } | |
82 | fo->Close(); | |
83 | } | |
84 | ||
85 | ||
86 | ||
87 |