]>
Commit | Line | Data |
---|---|---|
2803ac99 | 1 | //------------------------------------------------------------------------------ |
2 | // makeCompNSD.C | |
3 | // | |
4 | // create figure with ALICE comparison to CMS and ATLAS | |
5 | //------------------------------------------------------------------------------ | |
6 | ||
7 | ||
8 | ||
9 | void makeCompNSD() | |
10 | { | |
11 | ||
12 | gROOT->LoadMacro("divide.C"); | |
13 | ||
14 | TGraphErrors *graphCms = new TGraphErrors(binsCms,ptCms,nsdCms,errPtCms,err2NsdCms); | |
15 | graphCms->SetMarkerColor(colorCms); | |
16 | graphCms->SetLineColor(colorCms); | |
17 | graphCms->SetMarkerStyle(26); | |
18 | ||
19 | TGraphAsymmErrors *graphAtlas = new TGraphAsymmErrors(binsAtlas,ptAtlas,nsdAtlas,lowErrPtAtlas,highErrPtAtlas,err2NsdAtlas,err2NsdAtlas); | |
20 | graphAtlas->SetMarkerColor(colorAtlas); | |
21 | graphAtlas->SetLineColor(colorAtlas); | |
22 | graphAtlas->SetMarkerStyle(25); | |
23 | ||
24 | TGraphErrors *graphAlice = new TGraphErrors(binsNsdAlice,ptNsdAlice,nsdAlice,0,err2NsdAlice); | |
25 | graphAlice->SetMarkerColor(colorAliceFit); | |
26 | graphAlice->SetLineColor(colorAliceFit); | |
27 | graphAlice->SetMarkerStyle(20); | |
28 | ||
29 | TMultiGraph *all = new TMultiGraph("all",""); | |
30 | ||
31 | ||
32 | all->Add(graphAtlas); | |
33 | all->Add(graphCms); | |
34 | //all->Add(graphAlice); | |
35 | all->Add(graphAlice); | |
36 | //all->Draw("AP"); | |
37 | //setAttrib(all); | |
38 | ||
39 | ||
40 | TGraphErrors *ratioAtlasAlice = divide(graphAtlas,graphAlice); | |
41 | ratioAtlasAlice->SetFillColor(colorAtlas); | |
42 | ratioAtlasAlice->SetLineColor(colorAtlas); | |
43 | ratioAtlasAlice->SetFillStyle(3354); | |
44 | ||
45 | TGraphErrors *ratioCmsAlice = divide(graphCms,graphAlice); | |
46 | ratioCmsAlice->SetFillColor(colorCms); | |
47 | ratioCmsAlice->SetLineColor(colorCms); | |
48 | ratioCmsAlice->SetFillStyle(3345); | |
49 | ||
50 | Double_t one[binsNsdAlice]; | |
51 | for (Int_t i=0; i < binsNsdAlice; i++) { | |
52 | one[i] = 1.0; | |
53 | } | |
54 | ||
55 | TGraphErrors *ratioAliceAlice = new TGraphErrors(binsNsdAlice,ptNsdAlice,one,0,relErr2NsdAlice); | |
56 | ratioAliceAlice->SetFillColor(colorAliceErrors); | |
57 | ratioAliceAlice->SetLineColor(20); | |
58 | ||
59 | TCanvas *can3 = new TCanvas("can3","CompNSD",520,700); | |
60 | ||
61 | TPad *pad3_1 = new TPad("pad3_1","pad3_1",0.0,0.35,1.0,1.0); | |
62 | setAttrib(pad3_1); | |
63 | ||
64 | TPad *pad3_2 = new TPad("pad3_2","pad3_2",0.0,0.0,1.0,0.35); | |
65 | setAttrib(pad3_2); | |
66 | // | |
67 | can3->cd(); | |
68 | ||
69 | pad3_1->Draw(); | |
70 | pad3_1->cd(); | |
71 | pad3_1->SetLogx(); | |
72 | pad3_1->SetLogy(); | |
73 | //all->Draw("AE4"); | |
74 | all->Draw("APZ"); | |
75 | all->GetXaxis()->SetTitle("p_{T} [GeV/c]"); | |
76 | all->GetXaxis()->SetTitleOffset(1.6); | |
77 | all->GetYaxis()->SetTitleOffset(1.6); | |
78 | all->GetYaxis()->SetTitle("1/N_{evt} 1/(2#pi p_{T}) (d^{2}N_{ch})/(d#eta dp_{T}) (GeV/c)^{-2}"); | |
79 | all->GetXaxis()->SetLimits(minPt,maxPt); | |
80 | setAttrib(all); | |
81 | /* | |
82 | all->SetTitle(""); | |
83 | all->UseCurrentStyle(); | |
84 | all->GetXaxis()->SetRangeUser(minPt,maxPt); | |
85 | all->GetXaxis()->SetLabelOffset(0.0); | |
86 | all->GetXaxis()->SetTitleOffset(1.0); | |
87 | all->GetYaxis()->SetTitleOffset(1.5); | |
88 | all->GetYaxis()->SetTitleSize(0.06); | |
89 | all->GetYaxis()->SetLabelSize(0.06); | |
90 | all->GetXaxis()->SetNdivisions(405); | |
91 | all->GetXaxis()->SetLabelSize(0.0); | |
92 | */ | |
93 | ||
94 | ||
95 | ||
96 | graphAlice->SetFillColor(2); | |
97 | TLegend *l3_1 = new TLegend(0.2,0.03,0.65,0.261); | |
98 | l3_1->SetHeader("pp, NSD, #sqrt{s} = 900 GeV"); | |
99 | l3_1->AddEntry(graphAlice,"ALICE | #eta | < 0.8","LP"); | |
100 | l3_1->AddEntry(graphAtlas,"ATLAS | #eta | < 2.5","LP"); | |
101 | l3_1->AddEntry(graphCms,"CMS | #eta | < 2.4","LP"); | |
102 | //la->AddEntry(graphUA1,"UA1","p"); | |
103 | //la->AddEntry(graphCMS,"CMS","p"); | |
104 | //la->AddEntry(graphCMSred,"CMS, #eta<0.8","p"); | |
105 | //la->AddEntry(graphATLAS,"ATLAS","p"); | |
106 | l3_1->SetTextSize(legendTextSize); | |
107 | l3_1->SetFillColor(0); | |
108 | l3_1->SetLineColor(0); | |
109 | l3_1->Draw(); | |
110 | ||
111 | ||
112 | //histo1->GetXaxis()->SetNdivisions(405); | |
113 | ||
114 | can3->cd(); | |
115 | pad3_2->Draw(); | |
116 | pad3_2->cd(); | |
117 | pad3_2->SetLogx(); | |
118 | ||
119 | TMultiGraph *ratios = new TMultiGraph("ratios",""); | |
120 | ||
121 | ||
122 | ratios->Add(ratioAliceAlice); | |
123 | ratios->Add(ratioAtlasAlice); | |
124 | ratios->Add(ratioCmsAlice); | |
125 | ||
126 | ; | |
127 | //ratio->SetTitle("ratio ALICE/ATLAS (different #eta)"); | |
128 | ratios->Draw("AE3"); | |
129 | ratios->GetXaxis()->SetLimits(minPt,maxPt); | |
130 | ratios->GetXaxis()->SetTitle("p_{T} (GeV/c)"); | |
131 | ratios->GetXaxis()->SetTitleOffset(1.4); | |
132 | ratios->GetYaxis()->SetTitleOffset(0.9); | |
133 | ratios->GetYaxis()->SetTitle("ratio"); | |
134 | ratios->GetYaxis()->SetLabelSize(0.07); | |
135 | ratios->GetXaxis()->SetLabelSize(0.07); | |
136 | ratios->GetXaxis()->SetTitleSize(0.07); | |
137 | ratios->GetYaxis()->SetTitleSize(0.07); | |
138 | ratios->GetYaxis()->CenterTitle(); | |
139 | setAttrib(ratios); | |
140 | ||
141 | ratios->GetYaxis()->SetLabelOffset(0.023); | |
142 | ||
143 | //ratios->GetYaxis()->SetRangeUser(0.2,1.7); | |
144 | //ratios->GetYaxis()->SetLimits(0.2,1.7); | |
145 | ||
146 | /* | |
147 | ratios->UseCurrentStyle(); | |
148 | ratios->GetXaxis()->SetTitleOffset(1.2); | |
149 | ratios->GetYaxis()->SetTitleOffset(0.6); | |
150 | ratios->GetXaxis()->SetTitleSize(0.135); | |
151 | ratios->GetYaxis()->SetTitleSize(0.135); | |
152 | ratios->GetXaxis()->SetLabelSize(0.135); | |
153 | ratios->GetYaxis()->SetLabelSize(0.135); | |
154 | ratios->GetXaxis()->SetLabelOffset(0.025); | |
155 | ratios->GetYaxis()->SetLabelOffset(0.06); | |
156 | ratios->GetXaxis()->SetTickLength(0.09); | |
157 | ratios->GetYaxis()->SetLabelOffset(0.025); | |
158 | ratios->GetXaxis()->SetNdivisions(405); | |
159 | ratios->GetYaxis()->SetNdivisions(402); | |
160 | ratios->GetXaxis()->SetRangeUser(minPt,maxPt); | |
161 | */ | |
162 | ||
163 | graphAlice->SetFillColor(2); | |
164 | TLegend *l3_2 = new TLegend(0.2,0.323,0.65,0.637); | |
165 | l3_2->AddEntry(ratioAliceAlice,"ALICE uncertainties","F"); | |
166 | l3_2->AddEntry(ratioAtlasAlice,"ATLAS / ALICE","F"); | |
167 | l3_2->AddEntry(ratioCmsAlice,"CMS / ALICE","F"); | |
168 | l3_2->SetFillColor(0); | |
169 | l3_2->SetLineColor(0); | |
170 | l3_2->SetTextSize(legendTextSize); | |
171 | l3_2->Draw(); | |
172 | ||
173 | ||
174 | //aliceData->Draw("AE4"); | |
175 | TF1 *fOne = new TF1("one","1",0.1,20); | |
176 | fOne->SetRange(0.1,20); | |
177 | fOne->SetLineWidth(1); | |
178 | fOne->Draw("SAME"); | |
179 | //ratio->GetYaxis()->SetTitle("ratio ALICE/ATLAS"); | |
180 | /* | |
181 | histo1c->Draw(); | |
182 | histo1c->GetXaxis()->SetLabelSize(0.08); | |
183 | histo1c->GetYaxis()->SetLabelSize(0.08); | |
184 | histo1c->GetXaxis()->SetNdivisions(405); | |
185 | histo1c->GetYaxis()->SetNdivisions(405); | |
186 | */ | |
187 | ||
188 | ratios->Draw("E3"); | |
189 | fOne->Draw("SAME"); | |
190 | logoPrelim(can3); | |
191 | } |