| 6fa5a1ef |
1 | void fig4() |
| 2 | { |
| 3 | |
| 4 | |
| 5 | myOptions(); |
| 6 | gROOT->ForceStyle(); |
| 7 | TDatime now; |
| 8 | int iDate = now.GetDate(); |
| 9 | int iYear=iDate/10000; |
| 10 | int iMonth=(iDate%10000)/100; |
| 11 | int iDay=iDate%100; |
| 12 | char* cMonth[12]={"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"}; |
| 13 | char cStamp1[25],cStamp2[25]; |
| 14 | sprintf(cStamp1,"%i %s %i",iDay,cMonth[iMonth-1],iYear); |
| 15 | sprintf(cStamp2,"%i/%.2d/%i",iDay,iMonth,iYear); |
| 16 | |
| 17 | |
| 18 | |
| 19 | //=========Macro generated from canvas: c/c |
| 20 | //========= (Wed Apr 20 21:53:35 2011) by ROOT version5.27/05 |
| 21 | TCanvas *c = new TCanvas("c", "correlation function",4,23,800,700); |
| 22 | gStyle->SetOptStat(0); |
| 23 | |
| 24 | // ------------>Primitives in pad: c_1 |
| 25 | TPad *c_1 = new TPad("c_1", "c_1",0,0.,1,1); |
| 26 | c_1->Draw(); |
| 27 | c_1->cd(); |
| 28 | c_1->Range(-3.246312,0.9897024,5.131268,1.019464); |
| 29 | c_1->SetFillColor(0); |
| 30 | c_1->SetTickx(1); |
| 31 | c_1->SetTicky(1); |
| 32 | c_1->SetLeftMargin(0.2); |
| 33 | c_1->SetRightMargin(0.05); |
| 34 | c_1->SetTopMargin(0.02); |
| 35 | c_1->SetBottomMargin(0.15); |
| 36 | |
| 37 | |
| 38 | Double_t xAxis5[37] = {-1.570796, -1.396263, -1.22173, -1.047198, -0.8726646, -0.6981317, -0.5235988, -0.3490659, -0.1745329, 0, 0.1745329, 0.3490659, 0.5235988, 0.6981317, 0.8726646, 1.047198, 1.22173, 1.396263, 1.570796, 1.745329, 1.919862, 2.094395, 2.268928, 2.443461, 2.617994, 2.792527, 2.96706, 3.141593, 3.316126, 3.490659, 3.665191, 3.839724, 4.014257, 4.18879, 4.363323, 4.537856, 4.712389}; |
| 39 | |
| 40 | TH1D *hist = new TH1D("hist","",36, xAxis5); |
| 41 | hist->SetBinContent(1,0.9969217); |
| 42 | hist->SetBinContent(2,0.994212); |
| 43 | hist->SetBinContent(3,0.9924679); |
| 44 | hist->SetBinContent(4,0.9927132); |
| 45 | hist->SetBinContent(5,0.9941483); |
| 46 | hist->SetBinContent(6,0.99839); |
| 47 | hist->SetBinContent(7,1.004135); |
| 48 | hist->SetBinContent(8,1.009233); |
| 49 | hist->SetBinContent(9,1.012218); |
| 50 | hist->SetBinContent(10,1.012218); |
| 51 | hist->SetBinContent(11,1.009233); |
| 52 | hist->SetBinContent(12,1.004135); |
| 53 | hist->SetBinContent(13,0.99839); |
| 54 | hist->SetBinContent(14,0.9941483); |
| 55 | hist->SetBinContent(15,0.9927132); |
| 56 | hist->SetBinContent(16,0.9924679); |
| 57 | hist->SetBinContent(17,0.994212); |
| 58 | hist->SetBinContent(18,0.9969217); |
| 59 | hist->SetBinContent(19,0.9993932); |
| 60 | hist->SetBinContent(20,1.000563); |
| 61 | hist->SetBinContent(21,1.001791); |
| 62 | hist->SetBinContent(22,1.001702); |
| 63 | hist->SetBinContent(23,1.001143); |
| 64 | hist->SetBinContent(24,1.000711); |
| 65 | hist->SetBinContent(25,1.000244); |
| 66 | hist->SetBinContent(26,0.9996513); |
| 67 | hist->SetBinContent(27,0.9998416); |
| 68 | hist->SetBinContent(28,0.9998416); |
| 69 | hist->SetBinContent(29,0.9996513); |
| 70 | hist->SetBinContent(30,1.000244); |
| 71 | hist->SetBinContent(31,1.000711); |
| 72 | hist->SetBinContent(32,1.001143); |
| 73 | hist->SetBinContent(33,1.001702); |
| 74 | hist->SetBinContent(34,1.001791); |
| 75 | hist->SetBinContent(35,1.000563); |
| 76 | hist->SetBinContent(36,0.9993932); |
| 77 | hist->SetBinError(1,0.0001929295); |
| 78 | hist->SetBinError(2,0.0001921493); |
| 79 | hist->SetBinError(3,0.0001918335); |
| 80 | hist->SetBinError(4,0.000192209); |
| 81 | hist->SetBinError(5,0.0001920719); |
| 82 | hist->SetBinError(6,0.0001927733); |
| 83 | hist->SetBinError(7,0.0001938162); |
| 84 | hist->SetBinError(8,0.0001943988); |
| 85 | hist->SetBinError(9,0.0001947595); |
| 86 | hist->SetBinError(10,0.0001947595); |
| 87 | hist->SetBinError(11,0.0001943988); |
| 88 | hist->SetBinError(12,0.0001938162); |
| 89 | hist->SetBinError(13,0.0001927733); |
| 90 | hist->SetBinError(14,0.0001920719); |
| 91 | hist->SetBinError(15,0.000192209); |
| 92 | hist->SetBinError(16,0.0001918335); |
| 93 | hist->SetBinError(17,0.0001921493); |
| 94 | hist->SetBinError(18,0.0001929295); |
| 95 | hist->SetBinError(19,0.0001931266); |
| 96 | hist->SetBinError(20,0.00019325); |
| 97 | hist->SetBinError(21,0.0001934904); |
| 98 | hist->SetBinError(22,0.0001933326); |
| 99 | hist->SetBinError(23,0.0001931764); |
| 100 | hist->SetBinError(24,0.0001931073); |
| 101 | hist->SetBinError(25,0.0001930693); |
| 102 | hist->SetBinError(26,0.0001929318); |
| 103 | hist->SetBinError(27,0.0001934329); |
| 104 | hist->SetBinError(28,0.0001934329); |
| 105 | hist->SetBinError(29,0.0001929318); |
| 106 | hist->SetBinError(30,0.0001930693); |
| 107 | hist->SetBinError(31,0.0001931073); |
| 108 | hist->SetBinError(32,0.0001931764); |
| 109 | hist->SetBinError(33,0.0001933326); |
| 110 | hist->SetBinError(34,0.0001934904); |
| 111 | hist->SetBinError(35,0.00019325); |
| 112 | hist->SetBinError(36,0.0001931266); |
| 113 | hist->SetMinimum(0.99); |
| 114 | hist->SetMaximum(1.015); |
| 115 | hist->SetEntries(4.826814e+08); |
| 116 | hist->SetStats(0); |
| 117 | hist->SetLineWidth(2); |
| 118 | hist->GetXaxis()->SetTitle("#Delta#phi (rad.)"); |
| 119 | hist->GetXaxis()->SetLabelSize(0.05); |
| 120 | hist->GetXaxis()->SetTitleSize(0.06); |
| 121 | hist->GetYaxis()->SetTitle("C(#Delta#phi)"); |
| 122 | hist->GetYaxis()->SetLabelOffset(0.01); |
| 123 | hist->GetYaxis()->SetLabelSize(0.05); |
| 124 | hist->GetYaxis()->SetTitleSize(0.06); |
| 125 | hist->GetYaxis()->SetTitleOffset(1.7); |
| 126 | hist->GetZaxis()->SetLabelSize(0.05); |
| 127 | hist->GetZaxis()->SetTitleSize(0.06); |
| 128 | |
| 129 | hist->SetMarkerStyle(kOpenCircle); |
| 130 | //hist->Draw(""); |
| 131 | /* |
| 132 | |
| 133 | Double_t xAxis6[37] = {-1.570796, -1.396263, -1.22173, -1.047198, -0.8726646, -0.6981317, -0.5235988, -0.3490659, -0.1745329, 0, 0.1745329, 0.3490659, 0.5235988, 0.6981317, 0.8726646, 1.047198, 1.22173, 1.396263, 1.570796, 1.745329, 1.919862, 2.094395, 2.268928, 2.443461, 2.617994, 2.792527, 2.96706, 3.141593, 3.316126, 3.490659, 3.665191, 3.839724, 4.014257, 4.18879, 4.363323, 4.537856, 4.712389}; |
| 134 | |
| 135 | TH1D *GetDistAndFlow_1__2 = new TH1D("GetDistAndFlow_1__2","",36, xAxis6); |
| 136 | GetDistAndFlow_1__2->SetBinContent(1,0.9974817); |
| 137 | GetDistAndFlow_1__2->SetBinContent(2,0.9945214); |
| 138 | GetDistAndFlow_1__2->SetBinContent(3,0.9927582); |
| 139 | GetDistAndFlow_1__2->SetBinContent(4,0.993248); |
| 140 | GetDistAndFlow_1__2->SetBinContent(5,0.9938889); |
| 141 | GetDistAndFlow_1__2->SetBinContent(6,0.9976373); |
| 142 | GetDistAndFlow_1__2->SetBinContent(7,1.003408); |
| 143 | GetDistAndFlow_1__2->SetBinContent(8,1.007627); |
| 144 | GetDistAndFlow_1__2->SetBinContent(9,1.010439); |
| 145 | GetDistAndFlow_1__2->SetBinContent(10,1.010439); |
| 146 | GetDistAndFlow_1__2->SetBinContent(11,1.007627); |
| 147 | GetDistAndFlow_1__2->SetBinContent(12,1.003408); |
| 148 | GetDistAndFlow_1__2->SetBinContent(13,0.9976373); |
| 149 | GetDistAndFlow_1__2->SetBinContent(14,0.9938889); |
| 150 | GetDistAndFlow_1__2->SetBinContent(15,0.993248); |
| 151 | GetDistAndFlow_1__2->SetBinContent(16,0.9927582); |
| 152 | GetDistAndFlow_1__2->SetBinContent(17,0.9945214); |
| 153 | GetDistAndFlow_1__2->SetBinContent(18,0.9974817); |
| 154 | GetDistAndFlow_1__2->SetBinContent(19,1.000223); |
| 155 | GetDistAndFlow_1__2->SetBinContent(20,1.000993); |
| 156 | GetDistAndFlow_1__2->SetBinContent(21,1.002592); |
| 157 | GetDistAndFlow_1__2->SetBinContent(22,1.002206); |
| 158 | GetDistAndFlow_1__2->SetBinContent(23,1.001343); |
| 159 | GetDistAndFlow_1__2->SetBinContent(24,1.000875); |
| 160 | GetDistAndFlow_1__2->SetBinContent(25,1.000312); |
| 161 | GetDistAndFlow_1__2->SetBinContent(26,0.9996157); |
| 162 | GetDistAndFlow_1__2->SetBinContent(27,1.000308); |
| 163 | GetDistAndFlow_1__2->SetBinContent(28,1.000308); |
| 164 | GetDistAndFlow_1__2->SetBinContent(29,0.9996157); |
| 165 | GetDistAndFlow_1__2->SetBinContent(30,1.000312); |
| 166 | GetDistAndFlow_1__2->SetBinContent(31,1.000875); |
| 167 | GetDistAndFlow_1__2->SetBinContent(32,1.001343); |
| 168 | GetDistAndFlow_1__2->SetBinContent(33,1.002206); |
| 169 | GetDistAndFlow_1__2->SetBinContent(34,1.002592); |
| 170 | GetDistAndFlow_1__2->SetBinContent(35,1.000993); |
| 171 | GetDistAndFlow_1__2->SetBinContent(36,1.000223); |
| 172 | GetDistAndFlow_1__2->SetBinError(1,0.0003562986); |
| 173 | GetDistAndFlow_1__2->SetBinError(2,0.000354732); |
| 174 | GetDistAndFlow_1__2->SetBinError(3,0.0003541372); |
| 175 | GetDistAndFlow_1__2->SetBinError(4,0.0003549467); |
| 176 | GetDistAndFlow_1__2->SetBinError(5,0.000354547); |
| 177 | GetDistAndFlow_1__2->SetBinError(6,0.0003558383); |
| 178 | GetDistAndFlow_1__2->SetBinError(7,0.0003578351); |
| 179 | GetDistAndFlow_1__2->SetBinError(8,0.0003588045); |
| 180 | GetDistAndFlow_1__2->SetBinError(9,0.0003594313); |
| 181 | GetDistAndFlow_1__2->SetBinError(10,0.0003594313); |
| 182 | GetDistAndFlow_1__2->SetBinError(11,0.0003588045); |
| 183 | GetDistAndFlow_1__2->SetBinError(12,0.0003578351); |
| 184 | GetDistAndFlow_1__2->SetBinError(13,0.0003558383); |
| 185 | GetDistAndFlow_1__2->SetBinError(14,0.000354547); |
| 186 | GetDistAndFlow_1__2->SetBinError(15,0.0003549467); |
| 187 | GetDistAndFlow_1__2->SetBinError(16,0.0003541372); |
| 188 | GetDistAndFlow_1__2->SetBinError(17,0.000354732); |
| 189 | GetDistAndFlow_1__2->SetBinError(18,0.0003562986); |
| 190 | GetDistAndFlow_1__2->SetBinError(19,0.000356631); |
| 191 | GetDistAndFlow_1__2->SetBinError(20,0.0003568194); |
| 192 | GetDistAndFlow_1__2->SetBinError(21,0.0003573007); |
| 193 | GetDistAndFlow_1__2->SetBinError(22,0.0003569495); |
| 194 | GetDistAndFlow_1__2->SetBinError(23,0.0003566265); |
| 195 | GetDistAndFlow_1__2->SetBinError(24,0.0003564958); |
| 196 | GetDistAndFlow_1__2->SetBinError(25,0.0003564227); |
| 197 | GetDistAndFlow_1__2->SetBinError(26,0.0003561458); |
| 198 | GetDistAndFlow_1__2->SetBinError(27,0.0003572533); |
| 199 | GetDistAndFlow_1__2->SetBinError(28,0.0003572533); |
| 200 | GetDistAndFlow_1__2->SetBinError(29,0.0003561458); |
| 201 | GetDistAndFlow_1__2->SetBinError(30,0.0003564227); |
| 202 | GetDistAndFlow_1__2->SetBinError(31,0.0003564958); |
| 203 | GetDistAndFlow_1__2->SetBinError(32,0.0003566265); |
| 204 | GetDistAndFlow_1__2->SetBinError(33,0.0003569495); |
| 205 | GetDistAndFlow_1__2->SetBinError(34,0.0003573007); |
| 206 | GetDistAndFlow_1__2->SetBinError(35,0.0003568194); |
| 207 | GetDistAndFlow_1__2->SetBinError(36,0.000356631); |
| 208 | GetDistAndFlow_1__2->SetEntries(1.416101e+08); |
| 209 | GetDistAndFlow_1__2->SetDirectory(0); |
| 210 | GetDistAndFlow_1__2->SetStats(0); |
| 211 | GetDistAndFlow_1__2->SetLineColor(4); |
| 212 | GetDistAndFlow_1__2->SetLineWidth(2); |
| 213 | GetDistAndFlow_1__2->GetXaxis()->SetTitle("#Delta#phi (rad.)"); |
| 214 | GetDistAndFlow_1__2->GetXaxis()->SetLabelSize(0.05); |
| 215 | GetDistAndFlow_1__2->GetXaxis()->SetTitleSize(0.06); |
| 216 | GetDistAndFlow_1__2->GetYaxis()->SetTitle("C(#Delta#phi)"); |
| 217 | GetDistAndFlow_1__2->GetYaxis()->SetLabelOffset(0.01); |
| 218 | GetDistAndFlow_1__2->GetYaxis()->SetLabelSize(0.05); |
| 219 | GetDistAndFlow_1__2->GetYaxis()->SetTitleSize(0.06); |
| 220 | GetDistAndFlow_1__2->GetYaxis()->SetTitleOffset(1.3); |
| 221 | GetDistAndFlow_1__2->GetZaxis()->SetLabelSize(0.05); |
| 222 | GetDistAndFlow_1__2->GetZaxis()->SetTitleSize(0.06); |
| 223 | GetDistAndFlow_1__2->SetMarkerStyle(kFullCircle); |
| 224 | |
| 225 | GetDistAndFlow_1__2->Draw(""); |
| 226 | */ |
| 227 | |
| 228 | |
| 229 | Double_t xAxis2[37] = {-1.570796, -1.396263, -1.22173, -1.047198, -0.8726646, -0.6981317, -0.5235988, -0.3490659, -0.1745329, 0, 0.1745329, 0.3490659, 0.5235988, 0.6981317, 0.8726646, 1.047198, 1.22173, 1.396263, 1.570796, 1.745329, 1.919862, 2.094395, 2.268928, 2.443461, 2.617994, 2.792527, 2.96706, 3.141593, 3.316126, 3.490659, 3.665191, 3.839724, 4.014257, 4.18879, 4.363323, 4.537856, 4.712389}; |
| 230 | |
| 231 | TH1D *GetDistAndFlow_1__1 = new TH1D("GetDistAndFlow_1__1","",36, xAxis2); |
| 232 | GetDistAndFlow_1__1->SetBinContent(1,0.9978028); |
| 233 | GetDistAndFlow_1__1->SetBinContent(2,0.9945293); |
| 234 | GetDistAndFlow_1__1->SetBinContent(3,0.9926859); |
| 235 | GetDistAndFlow_1__1->SetBinContent(4,0.99349); |
| 236 | GetDistAndFlow_1__1->SetBinContent(5,0.9934892); |
| 237 | GetDistAndFlow_1__1->SetBinContent(6,0.9974816); |
| 238 | GetDistAndFlow_1__1->SetBinContent(7,1.003301); |
| 239 | GetDistAndFlow_1__1->SetBinContent(8,1.007196); |
| 240 | GetDistAndFlow_1__1->SetBinContent(9,1.010189); |
| 241 | GetDistAndFlow_1__1->SetBinContent(10,1.010189); |
| 242 | GetDistAndFlow_1__1->SetBinContent(11,1.007196); |
| 243 | GetDistAndFlow_1__1->SetBinContent(12,1.003301); |
| 244 | GetDistAndFlow_1__1->SetBinContent(13,0.9974816); |
| 245 | GetDistAndFlow_1__1->SetBinContent(14,0.9934892); |
| 246 | GetDistAndFlow_1__1->SetBinContent(15,0.99349); |
| 247 | GetDistAndFlow_1__1->SetBinContent(16,0.9926859); |
| 248 | GetDistAndFlow_1__1->SetBinContent(17,0.9945293); |
| 249 | GetDistAndFlow_1__1->SetBinContent(18,0.9978028); |
| 250 | GetDistAndFlow_1__1->SetBinContent(19,1.000258); |
| 251 | GetDistAndFlow_1__1->SetBinContent(20,1.000897); |
| 252 | GetDistAndFlow_1__1->SetBinContent(21,1.002466); |
| 253 | GetDistAndFlow_1__1->SetBinContent(22,1.001989); |
| 254 | GetDistAndFlow_1__1->SetBinContent(23,1.001441); |
| 255 | GetDistAndFlow_1__1->SetBinContent(24,1.000879); |
| 256 | GetDistAndFlow_1__1->SetBinContent(25,1.000151); |
| 257 | GetDistAndFlow_1__1->SetBinContent(26,0.99953); |
| 258 | GetDistAndFlow_1__1->SetBinContent(27,1.000522); |
| 259 | GetDistAndFlow_1__1->SetBinContent(28,1.000522); |
| 260 | GetDistAndFlow_1__1->SetBinContent(29,0.99953); |
| 261 | GetDistAndFlow_1__1->SetBinContent(30,1.000151); |
| 262 | GetDistAndFlow_1__1->SetBinContent(31,1.000879); |
| 263 | GetDistAndFlow_1__1->SetBinContent(32,1.001441); |
| 264 | GetDistAndFlow_1__1->SetBinContent(33,1.001989); |
| 265 | GetDistAndFlow_1__1->SetBinContent(34,1.002466); |
| 266 | GetDistAndFlow_1__1->SetBinContent(35,1.000897); |
| 267 | GetDistAndFlow_1__1->SetBinContent(36,1.000258); |
| 268 | GetDistAndFlow_1__1->SetBinError(1,0.0004580698); |
| 269 | GetDistAndFlow_1__1->SetBinError(2,0.0004559522); |
| 270 | GetDistAndFlow_1__1->SetBinError(3,0.0004551655); |
| 271 | GetDistAndFlow_1__1->SetBinError(4,0.0004563055); |
| 272 | GetDistAndFlow_1__1->SetBinError(5,0.0004556771); |
| 273 | GetDistAndFlow_1__1->SetBinError(6,0.0004573828); |
| 274 | GetDistAndFlow_1__1->SetBinError(7,0.0004599905); |
| 275 | GetDistAndFlow_1__1->SetBinError(8,0.0004611871); |
| 276 | GetDistAndFlow_1__1->SetBinError(9,0.0004620077); |
| 277 | GetDistAndFlow_1__1->SetBinError(10,0.0004620077); |
| 278 | GetDistAndFlow_1__1->SetBinError(11,0.0004611871); |
| 279 | GetDistAndFlow_1__1->SetBinError(12,0.0004599905); |
| 280 | GetDistAndFlow_1__1->SetBinError(13,0.0004573828); |
| 281 | GetDistAndFlow_1__1->SetBinError(14,0.0004556771); |
| 282 | GetDistAndFlow_1__1->SetBinError(15,0.0004563055); |
| 283 | GetDistAndFlow_1__1->SetBinError(16,0.0004551655); |
| 284 | GetDistAndFlow_1__1->SetBinError(17,0.0004559522); |
| 285 | GetDistAndFlow_1__1->SetBinError(18,0.0004580698); |
| 286 | GetDistAndFlow_1__1->SetBinError(19,0.0004584189); |
| 287 | GetDistAndFlow_1__1->SetBinError(20,0.0004586533); |
| 288 | GetDistAndFlow_1__1->SetBinError(21,0.0004592669); |
| 289 | GetDistAndFlow_1__1->SetBinError(22,0.0004587856); |
| 290 | GetDistAndFlow_1__1->SetBinError(23,0.0004584035); |
| 291 | GetDistAndFlow_1__1->SetBinError(24,0.0004582231); |
| 292 | GetDistAndFlow_1__1->SetBinError(25,0.000458117); |
| 293 | GetDistAndFlow_1__1->SetBinError(26,0.0004577663); |
| 294 | GetDistAndFlow_1__1->SetBinError(27,0.000459302); |
| 295 | GetDistAndFlow_1__1->SetBinError(28,0.000459302); |
| 296 | GetDistAndFlow_1__1->SetBinError(29,0.0004577663); |
| 297 | GetDistAndFlow_1__1->SetBinError(30,0.000458117); |
| 298 | GetDistAndFlow_1__1->SetBinError(31,0.0004582231); |
| 299 | GetDistAndFlow_1__1->SetBinError(32,0.0004584035); |
| 300 | GetDistAndFlow_1__1->SetBinError(33,0.0004587856); |
| 301 | GetDistAndFlow_1__1->SetBinError(34,0.0004592669); |
| 302 | GetDistAndFlow_1__1->SetBinError(35,0.0004586533); |
| 303 | GetDistAndFlow_1__1->SetBinError(36,0.0004584189); |
| 304 | //GetDistAndFlow_1__1->SetMinimum(0.9430516); |
| 305 | //GetDistAndFlow_1__1->SetMaximum(1.060699); |
| 306 | GetDistAndFlow_1__1->SetEntries(8.569633e+07); |
| 307 | GetDistAndFlow_1__1->SetDirectory(0); |
| 308 | GetDistAndFlow_1__1->SetStats(0); |
| 309 | GetDistAndFlow_1__1->SetLineColor(4); |
| 310 | GetDistAndFlow_1__1->SetLineWidth(2); |
| 311 | GetDistAndFlow_1__1->GetXaxis()->SetTitle("#Delta#phi (rad.)"); |
| 312 | GetDistAndFlow_1__1->GetXaxis()->SetLabelSize(0.05); |
| 313 | GetDistAndFlow_1__1->GetXaxis()->SetTitleSize(0.06); |
| 314 | GetDistAndFlow_1__1->GetYaxis()->SetTitle("C(#Delta#phi)"); |
| 315 | GetDistAndFlow_1__1->GetYaxis()->SetLabelOffset(0.01); |
| 316 | GetDistAndFlow_1__1->GetYaxis()->SetLabelSize(0.05); |
| 317 | GetDistAndFlow_1__1->GetYaxis()->SetTitleSize(0.06); |
| 318 | GetDistAndFlow_1__1->GetYaxis()->SetTitleOffset(1.3); |
| 319 | GetDistAndFlow_1__1->GetZaxis()->SetLabelSize(0.05); |
| 320 | GetDistAndFlow_1__1->GetZaxis()->SetTitleSize(0.06); |
| 321 | GetDistAndFlow_1__1->SetMarkerStyle(kFullCircle); |
| 322 | GetDistAndFlow_1__1->SetMarkerSize(1.2) |
| 323 | GetDistAndFlow_1__1->SetMarkerColor(kBlue); |
| 324 | |
| 325 | |
| 326 | |
| 327 | GetDistAndFlow_1__1->Draw(""); |
| 328 | |
| 329 | TF1 *flowFunc = new TF1("flowFunc","[0]*(1+2*[1]*cos(2*x)+2*[2]*cos(3*x)+2*[3]*cos(4*x)+2*[4]*cos(5*x))",-1.570796,4.712389); |
| 330 | flowFunc->SetFillColor(19); |
| 331 | flowFunc->SetFillStyle(0); |
| 332 | flowFunc->SetLineColor(2); |
| 333 | flowFunc->SetLineWidth(3); |
| 334 | flowFunc->SetLineStyle(1); |
| 335 | flowFunc->SetChisquare(62.12704); |
| 336 | flowFunc->SetNDF(35); |
| 337 | flowFunc->GetXaxis()->SetLabelSize(0.05); |
| 338 | flowFunc->GetXaxis()->SetTitleSize(0.06); |
| 339 | flowFunc->GetYaxis()->SetLabelOffset(0.01); |
| 340 | flowFunc->GetYaxis()->SetLabelSize(0.05); |
| 341 | flowFunc->GetYaxis()->SetTitleSize(0.06); |
| 342 | flowFunc->GetYaxis()->SetTitleOffset(1.3); |
| 343 | flowFunc->SetParameter(0,0.999971); |
| 344 | flowFunc->SetParError(0,5.941878e-05); |
| 345 | flowFunc->SetParLimits(0,0,0); |
| 346 | flowFunc->SetParameter(1,0.001639691); |
| 347 | flowFunc->SetParError(1,0); |
| 348 | flowFunc->SetParLimits(1,0.001639691,0.001639691); |
| 349 | flowFunc->SetParameter(2,0.002702756); |
| 350 | flowFunc->SetParError(2,0); |
| 351 | flowFunc->SetParLimits(2,0.002702756,0.002702756); |
| 352 | flowFunc->SetParameter(3,0.001165077); |
| 353 | flowFunc->SetParError(3,0); |
| 354 | flowFunc->SetParLimits(3,0.001165077,0.001165077); |
| 355 | flowFunc->SetParameter(4,0.0001550361); |
| 356 | flowFunc->SetParError(4,0); |
| 357 | flowFunc->SetParLimits(4,0.0001550361,0.0001550361); |
| 358 | flowFunc->Draw("SAME"); |
| 359 | |
| 360 | TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389); |
| 361 | flowFuncPart->SetFillColor(19); |
| 362 | flowFuncPart->SetFillStyle(0); |
| 363 | flowFuncPart->SetLineWidth(1); |
| 364 | flowFuncPart->GetXaxis()->SetLabelSize(0.05); |
| 365 | flowFuncPart->GetXaxis()->SetTitleSize(0.06); |
| 366 | flowFuncPart->GetYaxis()->SetLabelOffset(0.01); |
| 367 | flowFuncPart->GetYaxis()->SetLabelSize(0.05); |
| 368 | flowFuncPart->GetYaxis()->SetTitleSize(0.06); |
| 369 | flowFuncPart->GetYaxis()->SetTitleOffset(1.3); |
| 370 | flowFuncPart->SetParameter(0,0.999971); |
| 371 | flowFuncPart->SetParError(0,0); |
| 372 | flowFuncPart->SetParLimits(0,0,0); |
| 373 | flowFuncPart->SetParameter(1,0); |
| 374 | flowFuncPart->SetParError(1,0); |
| 375 | flowFuncPart->SetParLimits(1,0,0); |
| 376 | flowFuncPart->SetParameter(2,1); |
| 377 | flowFuncPart->SetParError(2,0); |
| 378 | flowFuncPart->SetParLimits(2,0,0); |
| 379 | flowFuncPart->Draw("SAME"); |
| 380 | |
| 381 | TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389); |
| 382 | flowFuncPart->SetFillColor(19); |
| 383 | flowFuncPart->SetFillStyle(0); |
| 384 | flowFuncPart->SetLineWidth(1); |
| 385 | flowFuncPart->SetLineStyle(2); |
| 386 | flowFuncPart->GetXaxis()->SetLabelSize(0.05); |
| 387 | flowFuncPart->GetXaxis()->SetTitleSize(0.06); |
| 388 | flowFuncPart->GetYaxis()->SetLabelOffset(0.01); |
| 389 | flowFuncPart->GetYaxis()->SetLabelSize(0.05); |
| 390 | flowFuncPart->GetYaxis()->SetTitleSize(0.06); |
| 391 | flowFuncPart->GetYaxis()->SetTitleOffset(1.3); |
| 392 | flowFuncPart->SetParameter(0,0.999971); |
| 393 | flowFuncPart->SetParError(0,0); |
| 394 | flowFuncPart->SetParLimits(0,0,0); |
| 395 | flowFuncPart->SetParameter(1,0.001639691); |
| 396 | flowFuncPart->SetParError(1,0); |
| 397 | flowFuncPart->SetParLimits(1,0,0); |
| 398 | flowFuncPart->SetParameter(2,2); |
| 399 | flowFuncPart->SetParError(2,0); |
| 400 | flowFuncPart->SetParLimits(2,0,0); |
| 401 | flowFuncPart->Draw("SAME"); |
| 402 | |
| 403 | TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389); |
| 404 | flowFuncPart->SetFillColor(19); |
| 405 | flowFuncPart->SetFillStyle(0); |
| 406 | flowFuncPart->SetLineWidth(1); |
| 407 | flowFuncPart->SetLineStyle(3); |
| 408 | flowFuncPart->GetXaxis()->SetLabelSize(0.05); |
| 409 | flowFuncPart->GetXaxis()->SetTitleSize(0.06); |
| 410 | flowFuncPart->GetYaxis()->SetLabelOffset(0.01); |
| 411 | flowFuncPart->GetYaxis()->SetLabelSize(0.05); |
| 412 | flowFuncPart->GetYaxis()->SetTitleSize(0.06); |
| 413 | flowFuncPart->GetYaxis()->SetTitleOffset(1.3); |
| 414 | flowFuncPart->SetParameter(0,0.999971); |
| 415 | flowFuncPart->SetParError(0,0); |
| 416 | flowFuncPart->SetParLimits(0,0,0); |
| 417 | flowFuncPart->SetParameter(1,0.002702756); |
| 418 | flowFuncPart->SetParError(1,0); |
| 419 | flowFuncPart->SetParLimits(1,0,0); |
| 420 | flowFuncPart->SetParameter(2,3); |
| 421 | flowFuncPart->SetParError(2,0); |
| 422 | flowFuncPart->SetParLimits(2,0,0); |
| 423 | flowFuncPart->Draw("SAME"); |
| 424 | |
| 425 | TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389); |
| 426 | flowFuncPart->SetFillColor(19); |
| 427 | flowFuncPart->SetFillStyle(0); |
| 428 | flowFuncPart->SetLineWidth(1); |
| 429 | flowFuncPart->SetLineStyle(4); |
| 430 | flowFuncPart->GetXaxis()->SetLabelSize(0.05); |
| 431 | flowFuncPart->GetXaxis()->SetTitleSize(0.06); |
| 432 | flowFuncPart->GetYaxis()->SetLabelOffset(0.01); |
| 433 | flowFuncPart->GetYaxis()->SetLabelSize(0.05); |
| 434 | flowFuncPart->GetYaxis()->SetTitleSize(0.06); |
| 435 | flowFuncPart->GetYaxis()->SetTitleOffset(1.3); |
| 436 | flowFuncPart->SetParameter(0,0.999971); |
| 437 | flowFuncPart->SetParError(0,0); |
| 438 | flowFuncPart->SetParLimits(0,0,0); |
| 439 | flowFuncPart->SetParameter(1,0.001165077); |
| 440 | flowFuncPart->SetParError(1,0); |
| 441 | flowFuncPart->SetParLimits(1,0,0); |
| 442 | flowFuncPart->SetParameter(2,4); |
| 443 | flowFuncPart->SetParError(2,0); |
| 444 | flowFuncPart->SetParLimits(2,0,0); |
| 445 | flowFuncPart->Draw("SAME"); |
| 446 | |
| 447 | TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389); |
| 448 | flowFuncPart->SetFillColor(19); |
| 449 | flowFuncPart->SetFillStyle(0); |
| 450 | flowFuncPart->SetLineWidth(1); |
| 451 | flowFuncPart->SetLineStyle(5); |
| 452 | flowFuncPart->GetXaxis()->SetLabelSize(0.05); |
| 453 | flowFuncPart->GetXaxis()->SetTitleSize(0.06); |
| 454 | flowFuncPart->GetYaxis()->SetLabelOffset(0.01); |
| 455 | flowFuncPart->GetYaxis()->SetLabelSize(0.05); |
| 456 | flowFuncPart->GetYaxis()->SetTitleSize(0.06); |
| 457 | flowFuncPart->GetYaxis()->SetTitleOffset(1.3); |
| 458 | flowFuncPart->SetParameter(0,0.999971); |
| 459 | flowFuncPart->SetParError(0,0); |
| 460 | flowFuncPart->SetParLimits(0,0,0); |
| 461 | flowFuncPart->SetParameter(1,0.0001550361); |
| 462 | flowFuncPart->SetParError(1,0); |
| 463 | flowFuncPart->SetParLimits(1,0,0); |
| 464 | flowFuncPart->SetParameter(2,5); |
| 465 | flowFuncPart->SetParError(2,0); |
| 466 | flowFuncPart->SetParLimits(2,0,0); |
| 467 | flowFuncPart->Draw("SAME"); |
| 468 | |
| 469 | //hist->SetMarkerStyle(kOpenCircle); |
| 470 | //hist->Draw("SAME"); |
| 471 | |
| 472 | |
| 473 | GetDistAndFlow_1__1->Draw("SAME"); |
| 474 | TLatex * tex = new TLatex(0.6,0.27,"2.0 < p_{t,trig} < 3.0"); |
| 475 | tex->SetNDC(); |
| 476 | tex->SetTextSize(0.04); |
| 477 | tex->SetLineWidth(2); |
| 478 | tex->Draw(); |
| 479 | tex = new TLatex(0.6,0.22,"1.0 < p_{t,assoc} < 2.0"); |
| 480 | tex->SetNDC(); |
| 481 | tex->SetTextSize(0.04); |
| 482 | tex->SetLineWidth(2); |
| 483 | tex->Draw(); |
| 484 | |
| 485 | |
| 486 | //TLegend *leg = new TLegend(0.48,0.82,0.91,0.95,NULL,"brNDC"); |
| 487 | TLegend *leg = new TLegend(0.48,0.78,0.91,0.95,"Centrality 0-1%, |#eta| < 0.8"); |
| 488 | myLegendSetUp(leg,0.04); |
| 489 | //leg->SetTextAlign(22); |
| 490 | //leg->SetTextFont(22); |
| 491 | //leg->SetTextSize(0.034); |
| 492 | //leg->SetFillColor(10); |
| 493 | //leg->SetFillStyle(0); |
| 494 | //leg->AddEntry("hist","All #Delta#eta","p"); |
| 495 | leg->AddEntry("GetDistAndFlow_1__1","|#Delta#eta| > 1","p"); |
| 496 | leg->AddEntry("flowFuncPart","v_{2,3,4,5}{2, |#Delta#eta| > 1}","l"); |
| 497 | leg->Draw(); |
| 498 | } |
| 499 | |
| 500 | |
| 501 | void myLegendSetUp(TLegend *currentLegend=0,float currentTextSize=0.07){ |
| 502 | currentLegend->SetTextFont(42); |
| 503 | currentLegend->SetBorderSize(0); |
| 504 | currentLegend->SetFillStyle(0); |
| 505 | currentLegend->SetFillColor(0); |
| 506 | currentLegend->SetMargin(0.25); |
| 507 | currentLegend->SetTextSize(currentTextSize); |
| 508 | currentLegend->SetEntrySeparation(0.5); |
| 509 | return; |
| 510 | } |
| 511 | |
| 512 | void myPadSetUp(TPad *currentPad, float currentLeft=0.11, float currentTop=0.04, float currentRight=0.04, float currentBottom=0.15){ |
| 513 | currentPad->SetLeftMargin(currentLeft); |
| 514 | currentPad->SetTopMargin(currentTop); |
| 515 | currentPad->SetRightMargin(currentRight); |
| 516 | currentPad->SetBottomMargin(currentBottom); |
| 517 | return; |
| 518 | } |
| 519 | |
| 520 | void myGraphSetUp(TGraphErrors *currentGraph=0, Float_t currentMarkerSize = 1.0, |
| 521 | int currentMarkerStyle=20, int currentMarkerColor=0, |
| 522 | int currentLineStyle=1, int currentLineColor=0, int currentFillColor=0, int currentFillStyle=0) |
| 523 | { |
| 524 | currentGraph->SetMarkerSize(currentMarkerSize); |
| 525 | currentGraph->SetMarkerStyle(currentMarkerStyle); |
| 526 | currentGraph->SetMarkerColor(currentMarkerColor); |
| 527 | currentGraph->SetLineStyle(currentLineStyle); |
| 528 | currentGraph->SetLineWidth(2); |
| 529 | currentGraph->SetLineColor(currentLineColor); |
| 530 | currentGraph->SetFillColor(currentFillColor); |
| 531 | currentGraph->SetFillStyle(currentFillStyle); |
| 532 | return; |
| 533 | } |
| 534 | |
| 535 | void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift) |
| 536 | { |
| 537 | // Shift original TGraphErrors along x-axis by Double_t shift. |
| 538 | |
| 539 | if(!ge){cout<<"!ge"<<endl;exit(0);} |
| 540 | |
| 541 | Int_t nPoints = ge->GetN(); |
| 542 | Double_t x = 0.; |
| 543 | Double_t y = 0.; |
| 544 | for(Int_t p=0;p<nPoints;p++) |
| 545 | { |
| 546 | ge->GetPoint(p,x,y); |
| 547 | x+=shift; |
| 548 | ge->SetPoint(p,x,y); |
| 549 | } // end of for(Int_t p=0;p<nPoints;p++) |
| 550 | |
| 551 | } // end of void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift) |
| 552 | |
| 553 | void myOptions(Int_t lStat=0){ |
| 554 | // Set gStyle |
| 555 | int font = 42; |
| 556 | // From plain |
| 557 | gStyle->SetFrameBorderMode(0); |
| 558 | gStyle->SetFrameFillColor(0); |
| 559 | gStyle->SetCanvasBorderMode(0); |
| 560 | gStyle->SetPadBorderMode(0); |
| 561 | gStyle->SetPadColor(10); |
| 562 | gStyle->SetCanvasColor(10); |
| 563 | gStyle->SetTitleFillColor(10); |
| 564 | gStyle->SetTitleBorderSize(1); |
| 565 | gStyle->SetStatColor(10); |
| 566 | gStyle->SetStatBorderSize(1); |
| 567 | gStyle->SetLegendBorderSize(1); |
| 568 | // |
| 569 | gStyle->SetDrawBorder(0); |
| 570 | gStyle->SetTextFont(font); |
| 571 | gStyle->SetStatFont(font); |
| 572 | gStyle->SetStatFontSize(0.05); |
| 573 | gStyle->SetStatX(0.97); |
| 574 | gStyle->SetStatY(0.98); |
| 575 | gStyle->SetStatH(0.03); |
| 576 | gStyle->SetStatW(0.3); |
| 577 | gStyle->SetTickLength(0.02,"y"); |
| 578 | gStyle->SetEndErrorSize(3); |
| 579 | gStyle->SetLabelSize(0.05,"xyz"); |
| 580 | gStyle->SetLabelFont(font,"xyz"); |
| 581 | gStyle->SetLabelOffset(0.01,"xyz"); |
| 582 | gStyle->SetTitleFont(font,"xyz"); |
| 583 | gStyle->SetTitleOffset(1.0,"xyz"); |
| 584 | gStyle->SetTitleSize(0.06,"xyz"); |
| 585 | gStyle->SetMarkerSize(1); |
| 586 | gStyle->SetPalette(1,0); |
| 587 | if (lStat){ |
| 588 | gStyle->SetOptTitle(1); |
| 589 | gStyle->SetOptStat(1111); |
| 590 | gStyle->SetOptFit(1111); |
| 591 | } |
| 592 | else { |
| 593 | gStyle->SetOptTitle(0); |
| 594 | gStyle->SetOptStat(0); |
| 595 | gStyle->SetOptFit(0); |
| 596 | } |
| 597 | } |
| 598 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame