PWGCF/FLOW/papers/PRL107_032301/fig4.C

   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