+++ /dev/null
-void fig4()
-{
-
-
- myOptions();
- gROOT->ForceStyle();
- TDatime now;
- int iDate = now.GetDate();
- int iYear=iDate/10000;
- int iMonth=(iDate%10000)/100;
- int iDay=iDate%100;
- char* cMonth[12]={"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"};
- char cStamp1[25],cStamp2[25];
- sprintf(cStamp1,"%i %s %i",iDay,cMonth[iMonth-1],iYear);
- sprintf(cStamp2,"%i/%.2d/%i",iDay,iMonth,iYear);
-
-
-
- //=========Macro generated from canvas: c/c
- //========= (Wed Apr 20 21:53:35 2011) by ROOT version5.27/05
- TCanvas *c = new TCanvas("c", "correlation function",4,23,800,700);
- gStyle->SetOptStat(0);
-
- // ------------>Primitives in pad: c_1
- TPad *c_1 = new TPad("c_1", "c_1",0,0.,1,1);
- c_1->Draw();
- c_1->cd();
- c_1->Range(-3.246312,0.9897024,5.131268,1.019464);
- c_1->SetFillColor(0);
- c_1->SetTickx(1);
- c_1->SetTicky(1);
- c_1->SetLeftMargin(0.2);
- c_1->SetRightMargin(0.05);
- c_1->SetTopMargin(0.02);
- c_1->SetBottomMargin(0.15);
-
-
- 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};
-
- TH1D *hist = new TH1D("hist","",36, xAxis5);
- hist->SetBinContent(1,0.9969217);
- hist->SetBinContent(2,0.994212);
- hist->SetBinContent(3,0.9924679);
- hist->SetBinContent(4,0.9927132);
- hist->SetBinContent(5,0.9941483);
- hist->SetBinContent(6,0.99839);
- hist->SetBinContent(7,1.004135);
- hist->SetBinContent(8,1.009233);
- hist->SetBinContent(9,1.012218);
- hist->SetBinContent(10,1.012218);
- hist->SetBinContent(11,1.009233);
- hist->SetBinContent(12,1.004135);
- hist->SetBinContent(13,0.99839);
- hist->SetBinContent(14,0.9941483);
- hist->SetBinContent(15,0.9927132);
- hist->SetBinContent(16,0.9924679);
- hist->SetBinContent(17,0.994212);
- hist->SetBinContent(18,0.9969217);
- hist->SetBinContent(19,0.9993932);
- hist->SetBinContent(20,1.000563);
- hist->SetBinContent(21,1.001791);
- hist->SetBinContent(22,1.001702);
- hist->SetBinContent(23,1.001143);
- hist->SetBinContent(24,1.000711);
- hist->SetBinContent(25,1.000244);
- hist->SetBinContent(26,0.9996513);
- hist->SetBinContent(27,0.9998416);
- hist->SetBinContent(28,0.9998416);
- hist->SetBinContent(29,0.9996513);
- hist->SetBinContent(30,1.000244);
- hist->SetBinContent(31,1.000711);
- hist->SetBinContent(32,1.001143);
- hist->SetBinContent(33,1.001702);
- hist->SetBinContent(34,1.001791);
- hist->SetBinContent(35,1.000563);
- hist->SetBinContent(36,0.9993932);
- hist->SetBinError(1,0.0001929295);
- hist->SetBinError(2,0.0001921493);
- hist->SetBinError(3,0.0001918335);
- hist->SetBinError(4,0.000192209);
- hist->SetBinError(5,0.0001920719);
- hist->SetBinError(6,0.0001927733);
- hist->SetBinError(7,0.0001938162);
- hist->SetBinError(8,0.0001943988);
- hist->SetBinError(9,0.0001947595);
- hist->SetBinError(10,0.0001947595);
- hist->SetBinError(11,0.0001943988);
- hist->SetBinError(12,0.0001938162);
- hist->SetBinError(13,0.0001927733);
- hist->SetBinError(14,0.0001920719);
- hist->SetBinError(15,0.000192209);
- hist->SetBinError(16,0.0001918335);
- hist->SetBinError(17,0.0001921493);
- hist->SetBinError(18,0.0001929295);
- hist->SetBinError(19,0.0001931266);
- hist->SetBinError(20,0.00019325);
- hist->SetBinError(21,0.0001934904);
- hist->SetBinError(22,0.0001933326);
- hist->SetBinError(23,0.0001931764);
- hist->SetBinError(24,0.0001931073);
- hist->SetBinError(25,0.0001930693);
- hist->SetBinError(26,0.0001929318);
- hist->SetBinError(27,0.0001934329);
- hist->SetBinError(28,0.0001934329);
- hist->SetBinError(29,0.0001929318);
- hist->SetBinError(30,0.0001930693);
- hist->SetBinError(31,0.0001931073);
- hist->SetBinError(32,0.0001931764);
- hist->SetBinError(33,0.0001933326);
- hist->SetBinError(34,0.0001934904);
- hist->SetBinError(35,0.00019325);
- hist->SetBinError(36,0.0001931266);
- hist->SetMinimum(0.99);
- hist->SetMaximum(1.015);
- hist->SetEntries(4.826814e+08);
- hist->SetStats(0);
- hist->SetLineWidth(2);
- hist->GetXaxis()->SetTitle("#Delta#phi (rad.)");
- hist->GetXaxis()->SetLabelSize(0.05);
- hist->GetXaxis()->SetTitleSize(0.06);
- hist->GetYaxis()->SetTitle("C(#Delta#phi)");
- hist->GetYaxis()->SetLabelOffset(0.01);
- hist->GetYaxis()->SetLabelSize(0.05);
- hist->GetYaxis()->SetTitleSize(0.06);
- hist->GetYaxis()->SetTitleOffset(1.7);
- hist->GetZaxis()->SetLabelSize(0.05);
- hist->GetZaxis()->SetTitleSize(0.06);
-
- hist->SetMarkerStyle(kOpenCircle);
- //hist->Draw("");
- /*
-
- 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};
-
- TH1D *GetDistAndFlow_1__2 = new TH1D("GetDistAndFlow_1__2","",36, xAxis6);
- GetDistAndFlow_1__2->SetBinContent(1,0.9974817);
- GetDistAndFlow_1__2->SetBinContent(2,0.9945214);
- GetDistAndFlow_1__2->SetBinContent(3,0.9927582);
- GetDistAndFlow_1__2->SetBinContent(4,0.993248);
- GetDistAndFlow_1__2->SetBinContent(5,0.9938889);
- GetDistAndFlow_1__2->SetBinContent(6,0.9976373);
- GetDistAndFlow_1__2->SetBinContent(7,1.003408);
- GetDistAndFlow_1__2->SetBinContent(8,1.007627);
- GetDistAndFlow_1__2->SetBinContent(9,1.010439);
- GetDistAndFlow_1__2->SetBinContent(10,1.010439);
- GetDistAndFlow_1__2->SetBinContent(11,1.007627);
- GetDistAndFlow_1__2->SetBinContent(12,1.003408);
- GetDistAndFlow_1__2->SetBinContent(13,0.9976373);
- GetDistAndFlow_1__2->SetBinContent(14,0.9938889);
- GetDistAndFlow_1__2->SetBinContent(15,0.993248);
- GetDistAndFlow_1__2->SetBinContent(16,0.9927582);
- GetDistAndFlow_1__2->SetBinContent(17,0.9945214);
- GetDistAndFlow_1__2->SetBinContent(18,0.9974817);
- GetDistAndFlow_1__2->SetBinContent(19,1.000223);
- GetDistAndFlow_1__2->SetBinContent(20,1.000993);
- GetDistAndFlow_1__2->SetBinContent(21,1.002592);
- GetDistAndFlow_1__2->SetBinContent(22,1.002206);
- GetDistAndFlow_1__2->SetBinContent(23,1.001343);
- GetDistAndFlow_1__2->SetBinContent(24,1.000875);
- GetDistAndFlow_1__2->SetBinContent(25,1.000312);
- GetDistAndFlow_1__2->SetBinContent(26,0.9996157);
- GetDistAndFlow_1__2->SetBinContent(27,1.000308);
- GetDistAndFlow_1__2->SetBinContent(28,1.000308);
- GetDistAndFlow_1__2->SetBinContent(29,0.9996157);
- GetDistAndFlow_1__2->SetBinContent(30,1.000312);
- GetDistAndFlow_1__2->SetBinContent(31,1.000875);
- GetDistAndFlow_1__2->SetBinContent(32,1.001343);
- GetDistAndFlow_1__2->SetBinContent(33,1.002206);
- GetDistAndFlow_1__2->SetBinContent(34,1.002592);
- GetDistAndFlow_1__2->SetBinContent(35,1.000993);
- GetDistAndFlow_1__2->SetBinContent(36,1.000223);
- GetDistAndFlow_1__2->SetBinError(1,0.0003562986);
- GetDistAndFlow_1__2->SetBinError(2,0.000354732);
- GetDistAndFlow_1__2->SetBinError(3,0.0003541372);
- GetDistAndFlow_1__2->SetBinError(4,0.0003549467);
- GetDistAndFlow_1__2->SetBinError(5,0.000354547);
- GetDistAndFlow_1__2->SetBinError(6,0.0003558383);
- GetDistAndFlow_1__2->SetBinError(7,0.0003578351);
- GetDistAndFlow_1__2->SetBinError(8,0.0003588045);
- GetDistAndFlow_1__2->SetBinError(9,0.0003594313);
- GetDistAndFlow_1__2->SetBinError(10,0.0003594313);
- GetDistAndFlow_1__2->SetBinError(11,0.0003588045);
- GetDistAndFlow_1__2->SetBinError(12,0.0003578351);
- GetDistAndFlow_1__2->SetBinError(13,0.0003558383);
- GetDistAndFlow_1__2->SetBinError(14,0.000354547);
- GetDistAndFlow_1__2->SetBinError(15,0.0003549467);
- GetDistAndFlow_1__2->SetBinError(16,0.0003541372);
- GetDistAndFlow_1__2->SetBinError(17,0.000354732);
- GetDistAndFlow_1__2->SetBinError(18,0.0003562986);
- GetDistAndFlow_1__2->SetBinError(19,0.000356631);
- GetDistAndFlow_1__2->SetBinError(20,0.0003568194);
- GetDistAndFlow_1__2->SetBinError(21,0.0003573007);
- GetDistAndFlow_1__2->SetBinError(22,0.0003569495);
- GetDistAndFlow_1__2->SetBinError(23,0.0003566265);
- GetDistAndFlow_1__2->SetBinError(24,0.0003564958);
- GetDistAndFlow_1__2->SetBinError(25,0.0003564227);
- GetDistAndFlow_1__2->SetBinError(26,0.0003561458);
- GetDistAndFlow_1__2->SetBinError(27,0.0003572533);
- GetDistAndFlow_1__2->SetBinError(28,0.0003572533);
- GetDistAndFlow_1__2->SetBinError(29,0.0003561458);
- GetDistAndFlow_1__2->SetBinError(30,0.0003564227);
- GetDistAndFlow_1__2->SetBinError(31,0.0003564958);
- GetDistAndFlow_1__2->SetBinError(32,0.0003566265);
- GetDistAndFlow_1__2->SetBinError(33,0.0003569495);
- GetDistAndFlow_1__2->SetBinError(34,0.0003573007);
- GetDistAndFlow_1__2->SetBinError(35,0.0003568194);
- GetDistAndFlow_1__2->SetBinError(36,0.000356631);
- GetDistAndFlow_1__2->SetEntries(1.416101e+08);
- GetDistAndFlow_1__2->SetDirectory(0);
- GetDistAndFlow_1__2->SetStats(0);
- GetDistAndFlow_1__2->SetLineColor(4);
- GetDistAndFlow_1__2->SetLineWidth(2);
- GetDistAndFlow_1__2->GetXaxis()->SetTitle("#Delta#phi (rad.)");
- GetDistAndFlow_1__2->GetXaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__2->GetXaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__2->GetYaxis()->SetTitle("C(#Delta#phi)");
- GetDistAndFlow_1__2->GetYaxis()->SetLabelOffset(0.01);
- GetDistAndFlow_1__2->GetYaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__2->GetYaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__2->GetYaxis()->SetTitleOffset(1.3);
- GetDistAndFlow_1__2->GetZaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__2->GetZaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__2->SetMarkerStyle(kFullCircle);
-
- GetDistAndFlow_1__2->Draw("");
- */
-
-
- 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};
-
- TH1D *GetDistAndFlow_1__1 = new TH1D("GetDistAndFlow_1__1","",36, xAxis2);
- GetDistAndFlow_1__1->SetBinContent(1,0.9978028);
- GetDistAndFlow_1__1->SetBinContent(2,0.9945293);
- GetDistAndFlow_1__1->SetBinContent(3,0.9926859);
- GetDistAndFlow_1__1->SetBinContent(4,0.99349);
- GetDistAndFlow_1__1->SetBinContent(5,0.9934892);
- GetDistAndFlow_1__1->SetBinContent(6,0.9974816);
- GetDistAndFlow_1__1->SetBinContent(7,1.003301);
- GetDistAndFlow_1__1->SetBinContent(8,1.007196);
- GetDistAndFlow_1__1->SetBinContent(9,1.010189);
- GetDistAndFlow_1__1->SetBinContent(10,1.010189);
- GetDistAndFlow_1__1->SetBinContent(11,1.007196);
- GetDistAndFlow_1__1->SetBinContent(12,1.003301);
- GetDistAndFlow_1__1->SetBinContent(13,0.9974816);
- GetDistAndFlow_1__1->SetBinContent(14,0.9934892);
- GetDistAndFlow_1__1->SetBinContent(15,0.99349);
- GetDistAndFlow_1__1->SetBinContent(16,0.9926859);
- GetDistAndFlow_1__1->SetBinContent(17,0.9945293);
- GetDistAndFlow_1__1->SetBinContent(18,0.9978028);
- GetDistAndFlow_1__1->SetBinContent(19,1.000258);
- GetDistAndFlow_1__1->SetBinContent(20,1.000897);
- GetDistAndFlow_1__1->SetBinContent(21,1.002466);
- GetDistAndFlow_1__1->SetBinContent(22,1.001989);
- GetDistAndFlow_1__1->SetBinContent(23,1.001441);
- GetDistAndFlow_1__1->SetBinContent(24,1.000879);
- GetDistAndFlow_1__1->SetBinContent(25,1.000151);
- GetDistAndFlow_1__1->SetBinContent(26,0.99953);
- GetDistAndFlow_1__1->SetBinContent(27,1.000522);
- GetDistAndFlow_1__1->SetBinContent(28,1.000522);
- GetDistAndFlow_1__1->SetBinContent(29,0.99953);
- GetDistAndFlow_1__1->SetBinContent(30,1.000151);
- GetDistAndFlow_1__1->SetBinContent(31,1.000879);
- GetDistAndFlow_1__1->SetBinContent(32,1.001441);
- GetDistAndFlow_1__1->SetBinContent(33,1.001989);
- GetDistAndFlow_1__1->SetBinContent(34,1.002466);
- GetDistAndFlow_1__1->SetBinContent(35,1.000897);
- GetDistAndFlow_1__1->SetBinContent(36,1.000258);
- GetDistAndFlow_1__1->SetBinError(1,0.0004580698);
- GetDistAndFlow_1__1->SetBinError(2,0.0004559522);
- GetDistAndFlow_1__1->SetBinError(3,0.0004551655);
- GetDistAndFlow_1__1->SetBinError(4,0.0004563055);
- GetDistAndFlow_1__1->SetBinError(5,0.0004556771);
- GetDistAndFlow_1__1->SetBinError(6,0.0004573828);
- GetDistAndFlow_1__1->SetBinError(7,0.0004599905);
- GetDistAndFlow_1__1->SetBinError(8,0.0004611871);
- GetDistAndFlow_1__1->SetBinError(9,0.0004620077);
- GetDistAndFlow_1__1->SetBinError(10,0.0004620077);
- GetDistAndFlow_1__1->SetBinError(11,0.0004611871);
- GetDistAndFlow_1__1->SetBinError(12,0.0004599905);
- GetDistAndFlow_1__1->SetBinError(13,0.0004573828);
- GetDistAndFlow_1__1->SetBinError(14,0.0004556771);
- GetDistAndFlow_1__1->SetBinError(15,0.0004563055);
- GetDistAndFlow_1__1->SetBinError(16,0.0004551655);
- GetDistAndFlow_1__1->SetBinError(17,0.0004559522);
- GetDistAndFlow_1__1->SetBinError(18,0.0004580698);
- GetDistAndFlow_1__1->SetBinError(19,0.0004584189);
- GetDistAndFlow_1__1->SetBinError(20,0.0004586533);
- GetDistAndFlow_1__1->SetBinError(21,0.0004592669);
- GetDistAndFlow_1__1->SetBinError(22,0.0004587856);
- GetDistAndFlow_1__1->SetBinError(23,0.0004584035);
- GetDistAndFlow_1__1->SetBinError(24,0.0004582231);
- GetDistAndFlow_1__1->SetBinError(25,0.000458117);
- GetDistAndFlow_1__1->SetBinError(26,0.0004577663);
- GetDistAndFlow_1__1->SetBinError(27,0.000459302);
- GetDistAndFlow_1__1->SetBinError(28,0.000459302);
- GetDistAndFlow_1__1->SetBinError(29,0.0004577663);
- GetDistAndFlow_1__1->SetBinError(30,0.000458117);
- GetDistAndFlow_1__1->SetBinError(31,0.0004582231);
- GetDistAndFlow_1__1->SetBinError(32,0.0004584035);
- GetDistAndFlow_1__1->SetBinError(33,0.0004587856);
- GetDistAndFlow_1__1->SetBinError(34,0.0004592669);
- GetDistAndFlow_1__1->SetBinError(35,0.0004586533);
- GetDistAndFlow_1__1->SetBinError(36,0.0004584189);
- //GetDistAndFlow_1__1->SetMinimum(0.9430516);
- //GetDistAndFlow_1__1->SetMaximum(1.060699);
- GetDistAndFlow_1__1->SetEntries(8.569633e+07);
- GetDistAndFlow_1__1->SetDirectory(0);
- GetDistAndFlow_1__1->SetStats(0);
- GetDistAndFlow_1__1->SetLineColor(4);
- GetDistAndFlow_1__1->SetLineWidth(2);
- GetDistAndFlow_1__1->GetXaxis()->SetTitle("#Delta#phi (rad.)");
- GetDistAndFlow_1__1->GetXaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__1->GetXaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__1->GetYaxis()->SetTitle("C(#Delta#phi)");
- GetDistAndFlow_1__1->GetYaxis()->SetLabelOffset(0.01);
- GetDistAndFlow_1__1->GetYaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__1->GetYaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__1->GetYaxis()->SetTitleOffset(1.3);
- GetDistAndFlow_1__1->GetZaxis()->SetLabelSize(0.05);
- GetDistAndFlow_1__1->GetZaxis()->SetTitleSize(0.06);
- GetDistAndFlow_1__1->SetMarkerStyle(kFullCircle);
- GetDistAndFlow_1__1->SetMarkerSize(1.2)
- GetDistAndFlow_1__1->SetMarkerColor(kBlue);
-
-
-
- GetDistAndFlow_1__1->Draw("");
-
- 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);
- flowFunc->SetFillColor(19);
- flowFunc->SetFillStyle(0);
- flowFunc->SetLineColor(2);
- flowFunc->SetLineWidth(3);
- flowFunc->SetLineStyle(1);
- flowFunc->SetChisquare(62.12704);
- flowFunc->SetNDF(35);
- flowFunc->GetXaxis()->SetLabelSize(0.05);
- flowFunc->GetXaxis()->SetTitleSize(0.06);
- flowFunc->GetYaxis()->SetLabelOffset(0.01);
- flowFunc->GetYaxis()->SetLabelSize(0.05);
- flowFunc->GetYaxis()->SetTitleSize(0.06);
- flowFunc->GetYaxis()->SetTitleOffset(1.3);
- flowFunc->SetParameter(0,0.999971);
- flowFunc->SetParError(0,5.941878e-05);
- flowFunc->SetParLimits(0,0,0);
- flowFunc->SetParameter(1,0.001639691);
- flowFunc->SetParError(1,0);
- flowFunc->SetParLimits(1,0.001639691,0.001639691);
- flowFunc->SetParameter(2,0.002702756);
- flowFunc->SetParError(2,0);
- flowFunc->SetParLimits(2,0.002702756,0.002702756);
- flowFunc->SetParameter(3,0.001165077);
- flowFunc->SetParError(3,0);
- flowFunc->SetParLimits(3,0.001165077,0.001165077);
- flowFunc->SetParameter(4,0.0001550361);
- flowFunc->SetParError(4,0);
- flowFunc->SetParLimits(4,0.0001550361,0.0001550361);
- flowFunc->Draw("SAME");
-
- TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389);
- flowFuncPart->SetFillColor(19);
- flowFuncPart->SetFillStyle(0);
- flowFuncPart->SetLineWidth(1);
- flowFuncPart->GetXaxis()->SetLabelSize(0.05);
- flowFuncPart->GetXaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetLabelOffset(0.01);
- flowFuncPart->GetYaxis()->SetLabelSize(0.05);
- flowFuncPart->GetYaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetTitleOffset(1.3);
- flowFuncPart->SetParameter(0,0.999971);
- flowFuncPart->SetParError(0,0);
- flowFuncPart->SetParLimits(0,0,0);
- flowFuncPart->SetParameter(1,0);
- flowFuncPart->SetParError(1,0);
- flowFuncPart->SetParLimits(1,0,0);
- flowFuncPart->SetParameter(2,1);
- flowFuncPart->SetParError(2,0);
- flowFuncPart->SetParLimits(2,0,0);
- flowFuncPart->Draw("SAME");
-
- TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389);
- flowFuncPart->SetFillColor(19);
- flowFuncPart->SetFillStyle(0);
- flowFuncPart->SetLineWidth(1);
- flowFuncPart->SetLineStyle(2);
- flowFuncPart->GetXaxis()->SetLabelSize(0.05);
- flowFuncPart->GetXaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetLabelOffset(0.01);
- flowFuncPart->GetYaxis()->SetLabelSize(0.05);
- flowFuncPart->GetYaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetTitleOffset(1.3);
- flowFuncPart->SetParameter(0,0.999971);
- flowFuncPart->SetParError(0,0);
- flowFuncPart->SetParLimits(0,0,0);
- flowFuncPart->SetParameter(1,0.001639691);
- flowFuncPart->SetParError(1,0);
- flowFuncPart->SetParLimits(1,0,0);
- flowFuncPart->SetParameter(2,2);
- flowFuncPart->SetParError(2,0);
- flowFuncPart->SetParLimits(2,0,0);
- flowFuncPart->Draw("SAME");
-
- TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389);
- flowFuncPart->SetFillColor(19);
- flowFuncPart->SetFillStyle(0);
- flowFuncPart->SetLineWidth(1);
- flowFuncPart->SetLineStyle(3);
- flowFuncPart->GetXaxis()->SetLabelSize(0.05);
- flowFuncPart->GetXaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetLabelOffset(0.01);
- flowFuncPart->GetYaxis()->SetLabelSize(0.05);
- flowFuncPart->GetYaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetTitleOffset(1.3);
- flowFuncPart->SetParameter(0,0.999971);
- flowFuncPart->SetParError(0,0);
- flowFuncPart->SetParLimits(0,0,0);
- flowFuncPart->SetParameter(1,0.002702756);
- flowFuncPart->SetParError(1,0);
- flowFuncPart->SetParLimits(1,0,0);
- flowFuncPart->SetParameter(2,3);
- flowFuncPart->SetParError(2,0);
- flowFuncPart->SetParLimits(2,0,0);
- flowFuncPart->Draw("SAME");
-
- TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389);
- flowFuncPart->SetFillColor(19);
- flowFuncPart->SetFillStyle(0);
- flowFuncPart->SetLineWidth(1);
- flowFuncPart->SetLineStyle(4);
- flowFuncPart->GetXaxis()->SetLabelSize(0.05);
- flowFuncPart->GetXaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetLabelOffset(0.01);
- flowFuncPart->GetYaxis()->SetLabelSize(0.05);
- flowFuncPart->GetYaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetTitleOffset(1.3);
- flowFuncPart->SetParameter(0,0.999971);
- flowFuncPart->SetParError(0,0);
- flowFuncPart->SetParLimits(0,0,0);
- flowFuncPart->SetParameter(1,0.001165077);
- flowFuncPart->SetParError(1,0);
- flowFuncPart->SetParLimits(1,0,0);
- flowFuncPart->SetParameter(2,4);
- flowFuncPart->SetParError(2,0);
- flowFuncPart->SetParLimits(2,0,0);
- flowFuncPart->Draw("SAME");
-
- TF1 *flowFuncPart = new TF1("flowFuncPart","[0]*(1+2*[1]*cos([2]*x))",-1.570796,4.712389);
- flowFuncPart->SetFillColor(19);
- flowFuncPart->SetFillStyle(0);
- flowFuncPart->SetLineWidth(1);
- flowFuncPart->SetLineStyle(5);
- flowFuncPart->GetXaxis()->SetLabelSize(0.05);
- flowFuncPart->GetXaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetLabelOffset(0.01);
- flowFuncPart->GetYaxis()->SetLabelSize(0.05);
- flowFuncPart->GetYaxis()->SetTitleSize(0.06);
- flowFuncPart->GetYaxis()->SetTitleOffset(1.3);
- flowFuncPart->SetParameter(0,0.999971);
- flowFuncPart->SetParError(0,0);
- flowFuncPart->SetParLimits(0,0,0);
- flowFuncPart->SetParameter(1,0.0001550361);
- flowFuncPart->SetParError(1,0);
- flowFuncPart->SetParLimits(1,0,0);
- flowFuncPart->SetParameter(2,5);
- flowFuncPart->SetParError(2,0);
- flowFuncPart->SetParLimits(2,0,0);
- flowFuncPart->Draw("SAME");
-
- //hist->SetMarkerStyle(kOpenCircle);
- //hist->Draw("SAME");
-
-
- GetDistAndFlow_1__1->Draw("SAME");
- TLatex * tex = new TLatex(0.6,0.27,"2.0 < p_{t,trig} < 3.0");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
- tex = new TLatex(0.6,0.22,"1.0 < p_{t,assoc} < 2.0");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
-
-
- //TLegend *leg = new TLegend(0.48,0.82,0.91,0.95,NULL,"brNDC");
- TLegend *leg = new TLegend(0.48,0.78,0.91,0.95,"Centrality 0-1%, |#eta| < 0.8");
- myLegendSetUp(leg,0.04);
- //leg->SetTextAlign(22);
- //leg->SetTextFont(22);
- //leg->SetTextSize(0.034);
- //leg->SetFillColor(10);
- //leg->SetFillStyle(0);
- //leg->AddEntry("hist","All #Delta#eta","p");
- leg->AddEntry("GetDistAndFlow_1__1","|#Delta#eta| > 1","p");
- leg->AddEntry("flowFuncPart","v_{2,3,4,5}{2, |#Delta#eta| > 1}","l");
- leg->Draw();
-}
-
-
-void myLegendSetUp(TLegend *currentLegend=0,float currentTextSize=0.07){
- currentLegend->SetTextFont(42);
- currentLegend->SetBorderSize(0);
- currentLegend->SetFillStyle(0);
- currentLegend->SetFillColor(0);
- currentLegend->SetMargin(0.25);
- currentLegend->SetTextSize(currentTextSize);
- currentLegend->SetEntrySeparation(0.5);
- return;
-}
-
-void myPadSetUp(TPad *currentPad, float currentLeft=0.11, float currentTop=0.04, float currentRight=0.04, float currentBottom=0.15){
- currentPad->SetLeftMargin(currentLeft);
- currentPad->SetTopMargin(currentTop);
- currentPad->SetRightMargin(currentRight);
- currentPad->SetBottomMargin(currentBottom);
- return;
-}
-
-void myGraphSetUp(TGraphErrors *currentGraph=0, Float_t currentMarkerSize = 1.0,
- int currentMarkerStyle=20, int currentMarkerColor=0,
- int currentLineStyle=1, int currentLineColor=0, int currentFillColor=0, int currentFillStyle=0)
-{
- currentGraph->SetMarkerSize(currentMarkerSize);
- currentGraph->SetMarkerStyle(currentMarkerStyle);
- currentGraph->SetMarkerColor(currentMarkerColor);
- currentGraph->SetLineStyle(currentLineStyle);
- currentGraph->SetLineWidth(2);
- currentGraph->SetLineColor(currentLineColor);
- currentGraph->SetFillColor(currentFillColor);
- currentGraph->SetFillStyle(currentFillStyle);
- return;
-}
-
-void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift)
-{
- // Shift original TGraphErrors along x-axis by Double_t shift.
-
- if(!ge){cout<<"!ge"<<endl;exit(0);}
-
- Int_t nPoints = ge->GetN();
- Double_t x = 0.;
- Double_t y = 0.;
- for(Int_t p=0;p<nPoints;p++)
- {
- ge->GetPoint(p,x,y);
- x+=shift;
- ge->SetPoint(p,x,y);
- } // end of for(Int_t p=0;p<nPoints;p++)
-
-} // end of void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift)
-
-void myOptions(Int_t lStat=0){
- // Set gStyle
- int font = 42;
- // From plain
- gStyle->SetFrameBorderMode(0);
- gStyle->SetFrameFillColor(0);
- gStyle->SetCanvasBorderMode(0);
- gStyle->SetPadBorderMode(0);
- gStyle->SetPadColor(10);
- gStyle->SetCanvasColor(10);
- gStyle->SetTitleFillColor(10);
- gStyle->SetTitleBorderSize(1);
- gStyle->SetStatColor(10);
- gStyle->SetStatBorderSize(1);
- gStyle->SetLegendBorderSize(1);
- //
- gStyle->SetDrawBorder(0);
- gStyle->SetTextFont(font);
- gStyle->SetStatFont(font);
- gStyle->SetStatFontSize(0.05);
- gStyle->SetStatX(0.97);
- gStyle->SetStatY(0.98);
- gStyle->SetStatH(0.03);
- gStyle->SetStatW(0.3);
- gStyle->SetTickLength(0.02,"y");
- gStyle->SetEndErrorSize(3);
- gStyle->SetLabelSize(0.05,"xyz");
- gStyle->SetLabelFont(font,"xyz");
- gStyle->SetLabelOffset(0.01,"xyz");
- gStyle->SetTitleFont(font,"xyz");
- gStyle->SetTitleOffset(1.0,"xyz");
- gStyle->SetTitleSize(0.06,"xyz");
- gStyle->SetMarkerSize(1);
- gStyle->SetPalette(1,0);
- if (lStat){
- gStyle->SetOptTitle(1);
- gStyle->SetOptStat(1111);
- gStyle->SetOptFit(1111);
- }
- else {
- gStyle->SetOptTitle(0);
- gStyle->SetOptStat(0);
- gStyle->SetOptFit(0);
- }
-}
-