]>
Commit | Line | Data |
---|---|---|
24373b38 | 1 | AddTaskPidQC(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){ |
2 | gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDQC.C"); | |
3 | ||
4 | const Int_t ncentr = 9; | |
5 | Int_t cmin[ncentr]={0,5,10,20,30,40,50,60,70}; | |
6 | Int_t cmax[ncentr]={5,10,20,30,40,50,60,70,80}; | |
7 | ||
8 | for(Int_t i=0;i < ncentr;i++){ | |
9 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,1,etamin,etamax); // no pid + QA | |
10 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); | |
11 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); | |
12 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax); | |
13 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); | |
14 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax); | |
15 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax); | |
16 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax); | |
17 | AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax); | |
18 | } | |
19 | } | |
20 | ||
21 | createQC4res(){ | |
4070f709 | 22 | gSystem->Load("libVMC"); |
23 | gSystem->Load("libPhysics"); | |
24 | gSystem->Load("libTree"); | |
25 | gSystem->Load("libMinuit"); | |
26 | gSystem->Load("libSTEERBase"); | |
27 | gSystem->Load("libANALYSIS"); | |
28 | gSystem->Load("libAOD"); | |
29 | gSystem->Load("libESD"); | |
30 | gSystem->Load("libANALYSIS"); | |
31 | gSystem->Load("libANALYSISalice"); | |
32 | gSystem->Load("libCORRFW"); | |
33 | gSystem->Load("libNetx"); | |
34 | gSystem->Load("libPWGflowBase"); | |
24373b38 | 35 | |
36 | char name[200]; | |
37 | char *spe[4]={"pion","kaon","antipr","proton"}; | |
38 | char *tech[2]={"TOF","TPC"}; | |
39 | Int_t cmin[9]={0,5,10,20,30,40,50,60,70}; | |
40 | Int_t cmax[9]={5,10,20,30,40,50,60,70,80}; | |
41 | ||
42 | TFile *f = new TFile("AnalysisResults.root"); | |
43 | TFile *fo = new TFile("resultsQC4.root","RECREATE"); | |
44 | TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputQCanalysisTPCstandalone")); | |
45 | TList* listTemp = directory->GetListOfKeys(); | |
46 | for(Int_t i=0;i < 9;i++){ | |
47 | TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->GetName())); | |
48 | AliFlowAnalysisWithQCumulants* sp2 = new AliFlowAnalysisWithQCumulants(); | |
49 | sp2->GetOutputHistograms(list2); | |
50 | sp2->Finish(); | |
51 | AliFlowCommonHistResults* res2=sp2->GetCommonHistsResults4th(); | |
52 | TH1D *h2=res2->GetHistDiffFlowPtPOI(); | |
53 | sprintf(name,"v2QC4_%s_%i_%i","AllCharged",cmin[i],cmax[i]); | |
54 | h2->SetName(name); | |
55 | fo->cd(); | |
56 | h2->Write(); | |
57 | for(Int_t j=0;j < 2;j++){ | |
58 | for(Int_t k=0;k < 4;k++){ | |
59 | TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9+j*4+k+1)->GetName())); | |
60 | AliFlowAnalysisWithQCumulants* sp = new AliFlowAnalysisWithQCumulants(); | |
61 | sp->GetOutputHistograms(list); | |
62 | sp->Finish(); | |
63 | AliFlowCommonHistResults* res=sp->GetCommonHistsResults4th(); | |
64 | TH1D *h=res->GetHistDiffFlowPtPOI(); | |
65 | sprintf(name,"v2QC4_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]); | |
66 | h->SetName(name); | |
67 | fo->cd(); | |
68 | h->Write(); | |
69 | } | |
70 | } | |
71 | } | |
72 | fo->Close(); | |
73 | } | |
74 | ||
75 | ||
76 | ||
77 | createQC2res(){ | |
4070f709 | 78 | gSystem->Load("libVMC"); |
79 | gSystem->Load("libPhysics"); | |
80 | gSystem->Load("libTree"); | |
81 | gSystem->Load("libMinuit"); | |
82 | gSystem->Load("libSTEERBase"); | |
83 | gSystem->Load("libANALYSIS"); | |
84 | gSystem->Load("libAOD"); | |
85 | gSystem->Load("libESD"); | |
86 | gSystem->Load("libANALYSIS"); | |
87 | gSystem->Load("libANALYSISalice"); | |
88 | gSystem->Load("libCORRFW"); | |
89 | gSystem->Load("libNetx"); | |
90 | gSystem->Load("libPWGflowBase"); | |
24373b38 | 91 | |
92 | char name[200]; | |
93 | char *spe[4]={"pion","kaon","antipr","proton"}; | |
94 | char *tech[2]={"TOF","TPC"}; | |
95 | Int_t cmin[9]={0,5,10,20,30,40,50,60,70}; | |
96 | Int_t cmax[9]={5,10,20,30,40,50,60,70,80}; | |
97 | ||
98 | TFile *f = new TFile("AnalysisResults.root"); | |
99 | TFile *fo = new TFile("resultsQC2.root","RECREATE"); | |
100 | TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputQCanalysisTPCstandalone")); | |
101 | TList* listTemp = directory->GetListOfKeys(); | |
102 | for(Int_t i=0;i < 9;i++){ | |
103 | TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->GetName())); | |
104 | AliFlowAnalysisWithQCumulants* sp2 = new AliFlowAnalysisWithQCumulants(); | |
105 | sp2->GetOutputHistograms(list2); | |
106 | sp2->Finish(); | |
107 | AliFlowCommonHistResults* res2=sp2->GetCommonHistsResults2nd(); | |
108 | TH1D *h2=res2->GetHistDiffFlowPtPOI(); | |
109 | sprintf(name,"v2QC2_%s_%i_%i","AllCharged",cmin[i],cmax[i]); | |
110 | h2->SetName(name); | |
111 | fo->cd(); | |
112 | h2->Write(); | |
113 | for(Int_t j=0;j < 2;j++){ | |
114 | for(Int_t k=0;k < 4;k++){ | |
115 | TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9+j*4+k+1)->GetName())); | |
116 | AliFlowAnalysisWithQCumulants* sp = new AliFlowAnalysisWithQCumulants(); | |
117 | sp->GetOutputHistograms(list); | |
118 | sp->Finish(); | |
119 | AliFlowCommonHistResults* res=sp->GetCommonHistsResults2nd(); | |
120 | TH1D *h=res->GetHistDiffFlowPtPOI(); | |
121 | sprintf(name,"v2QC2_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]); | |
122 | h->SetName(name); | |
123 | fo->cd(); | |
124 | h->Write(); | |
125 | } | |
126 | } | |
127 | } | |
128 | fo->Close(); | |
129 | } | |
130 | ||
131 | ||
132 |