void ExpManager::DrawMap(TString NameTitle, double xmin, double xmax, double ymin, double ymax, double zmin, double zmax) {
TMultiGraph *mg = new TMultiGraph();
//Draw a cross
TGraphErrors *frame = new TGraphErrors(); //= new TGraph2DErrors(np, x_array, y_array, bz_array, ex, ey, ez);
frame->SetPoint(0,+0,+100);
frame->SetPoint(1,+0,-100);
frame->SetPoint(2,0,0);
frame->SetPoint(3,-100,0);
frame->SetPoint(4,+100,0);
frame->SetPoint(5,+0,+0);
frame->SetPoint(6,+0,-100);
frame->SetMarkerColor(kWhite);
frame->SetDrawOption("ap");
//Draw the map
TGraphErrors *fGraph = new TGraphErrors(); //= new TGraph2DErrors(np, x_array, y_array, bz_array, ex, ey, ez);
fGraph->SetMarkerSize(1.2);
fGraph->SetMarkerStyle(20);
fGraph->SetMarkerColor(kBlue);
fGraph->SetLineColor(kBlue);
fGraph->SetLineWidth(1);
//fGraph->SetDrawOption("ap");
int graph_counter = 0 ;
for (unsigned i=0; i< fExpY.size(); i++) {
if( (fExpX.at(i) >= xmin && fExpX.at(i) <= xmax) && (fExpY.at(i) >= ymin && fExpY.at(i) <= ymax) && (fExpZ.at(i) >= zmin && fExpZ.at(i) <= zmax) ) {
fGraph->SetPoint(graph_counter,fExpX.at(i),fExpY.at(i));
fGraph->SetPointError(graph_counter,fExpXErr.at(i),fExpYErr.at(i));
graph_counter++;
}
}
//fGraph->Write();
mg->SetTitle(NameTitle+Form(" Map : %.2f < X < %.2f mm %.2f < Y < %.2f mm %.2f < Z < %.2f mm;X (mm);Y (mm)",xmin,xmax,ymin,ymax,zmin,zmax));
mg->SetName(NameTitle+Form("_Map__X_%.2f_%.2fmm__Y_%.2f_%.2fmm__Z_%.2f_%.2fmm",xmin,xmax,ymin,ymax,zmin,zmax));
mg->Add(frame);
mg->Add(fGraph);
mg->Write();
}
void result_JES_akPu4PF_(refpt> 75 && refpt < 120)&&(abs(refeta)<2)_Centrality()
{
//=========Macro generated from canvas: Can_result_0/
//========= (Mon Apr 18 02:39:41 2016) by ROOT version6.02/13
TCanvas *Can_result_0 = new TCanvas("Can_result_0", "",18,40,700,500);
gStyle->SetOptFit(1);
Can_result_0->Range(-18.75,0.9236134,118.75,1.043845);
Can_result_0->SetFillColor(0);
Can_result_0->SetBorderMode(0);
Can_result_0->SetBorderSize(2);
Can_result_0->SetFrameBorderMode(0);
Can_result_0->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("name");
multigraph->SetTitle("JES_akPu4PF");
Double_t Graph_fx1001[4] = {
5,
20,
40,
75};
Double_t Graph_fy1001[4] = {
1.025461,
0.982262,
0.9644553,
0.970611};
Double_t Graph_fex1001[4] = {
5,
10,
10,
25};
Double_t Graph_fey1001[4] = {
0.001598024,
0.0009318739,
0.0008002418,
0.0004502591};
TGraphErrors *gre = new TGraphErrors(4,Graph_fx1001,Graph_fy1001,Graph_fex1001,Graph_fey1001);
gre->SetName("Graph");
gre->SetTitle("some title_0");
gre->SetFillColor(1);
gre->SetFillStyle(0);
TH1F *Graph_Graph1001 = new TH1F("Graph_Graph1001","some title_0",100,0,110);
Graph_Graph1001->SetMinimum(0.9573146);
Graph_Graph1001->SetMaximum(1.0334);
Graph_Graph1001->SetDirectory(0);
Graph_Graph1001->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
Graph_Graph1001->SetLineColor(ci);
Graph_Graph1001->GetXaxis()->SetLabelFont(42);
Graph_Graph1001->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1001->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1001->GetXaxis()->SetTitleFont(42);
Graph_Graph1001->GetYaxis()->SetLabelFont(42);
Graph_Graph1001->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1001->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1001->GetYaxis()->SetTitleFont(42);
Graph_Graph1001->GetZaxis()->SetLabelFont(42);
Graph_Graph1001->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1001->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1001->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1001);
multigraph->Add(gre,"");
Double_t Graph_fx1002[4] = {
5,
20,
40,
75};
Double_t Graph_fy1002[4] = {
1.025808,
0.9823451,
0.964104,
0.9707841};
Double_t Graph_fex1002[4] = {
5,
10,
10,
25};
Double_t Graph_fey1002[4] = {
0.001641175,
0.0009556419,
0.0008206184,
0.000461058};
gre = new TGraphErrors(4,Graph_fx1002,Graph_fy1002,Graph_fex1002,Graph_fey1002);
gre->SetName("Graph");
gre->SetTitle("some title_1");
gre->SetFillColor(1);
gre->SetFillStyle(0);
gre->SetLineColor(2);
gre->SetMarkerColor(2);
TH1F *Graph_Graph1002 = new TH1F("Graph_Graph1002","some title_1",100,0,110);
Graph_Graph1002->SetMinimum(0.9568668);
Graph_Graph1002->SetMaximum(1.033866);
Graph_Graph1002->SetDirectory(0);
Graph_Graph1002->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1002->SetLineColor(ci);
Graph_Graph1002->GetXaxis()->SetLabelFont(42);
Graph_Graph1002->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1002->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1002->GetXaxis()->SetTitleFont(42);
Graph_Graph1002->GetYaxis()->SetLabelFont(42);
Graph_Graph1002->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1002->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1002->GetYaxis()->SetTitleFont(42);
Graph_Graph1002->GetZaxis()->SetLabelFont(42);
Graph_Graph1002->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1002->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1002->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1002);
multigraph->Add(gre,"");
Double_t Graph_fx1003[4] = {
5,
20,
40,
75};
Double_t Graph_fy1003[4] = {
0.9944711,
0.9453074,
0.9619222,
0.9556338};
Double_t Graph_fex1003[4] = {
5,
10,
10,
25};
Double_t Graph_fey1003[4] = {
0.01053149,
0.005298861,
0.004412705,
0.002663901};
gre = new TGraphErrors(4,Graph_fx1003,Graph_fy1003,Graph_fex1003,Graph_fey1003);
gre->SetName("Graph");
gre->SetTitle("some title_2");
gre->SetFillColor(1);
gre->SetFillStyle(0);
gre->SetLineColor(3);
gre->SetMarkerColor(3);
TH1F *Graph_Graph1003 = new TH1F("Graph_Graph1003","some title_2",100,0,110);
Graph_Graph1003->SetMinimum(0.9335092);
Graph_Graph1003->SetMaximum(1.011502);
Graph_Graph1003->SetDirectory(0);
Graph_Graph1003->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1003->SetLineColor(ci);
Graph_Graph1003->GetXaxis()->SetLabelFont(42);
Graph_Graph1003->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1003->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1003->GetXaxis()->SetTitleFont(42);
Graph_Graph1003->GetYaxis()->SetLabelFont(42);
Graph_Graph1003->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1003->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1003->GetYaxis()->SetTitleFont(42);
Graph_Graph1003->GetZaxis()->SetLabelFont(42);
Graph_Graph1003->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1003->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1003->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1003);
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("Centrality");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("#mu_{Reco./Gen.}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLegend *leg = new TLegend(0.68,0.7,0.88,0.85,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("Graph","allJets","l");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","lightJets","l");
entry->SetLineColor(2);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","bJets","l");
entry->SetLineColor(3);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.25,0.7,"#bf{#color[2]{|#eta_{jet}|<2.0 &75<refpt<120}}");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(24);
tex->SetLineWidth(2);
tex->Draw();
TPaveText *pt = new TPaveText(0.3750575,0.9365254,0.6249425,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("JES_akPu4PF");
pt->Draw();
Can_result_0->Modified();
Can_result_0->cd();
Can_result_0->SetSelected(Can_result_0);
}
void Final_av_Lambda_VarCoupling_40Percent()
{
//=========Macro generated from canvas: c/c
//========= (Sat Feb 27 17:06:04 2016) by ROOT version6.04/06
TCanvas *c = new TCanvas("c", "c",0,0,800,700);
c->SetHighLightColor(2);
c->Range(-0.4813593,1.36281,3.148656,6.672724);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetLogx();
c->SetLogy();
c->SetTickx(1);
c->SetTicky(1);
c->SetLeftMargin(0.12);
c->SetRightMargin(0.04);
c->SetTopMargin(0.07);
c->SetBottomMargin(0.12);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t Graph_fx1[12] = {
1,
10,
100,
400,
700,
1000,
1000,
700,
400,
100,
10,
1};
Double_t Graph_fy1[12] = {
881.525,
896.549,
871.907,
667.806,
471.744,
308.777,
382.832,
584.437,
826.374,
1079.16,
1110.85,
1092.23};
TGraph *graph = new TGraph(12,Graph_fx1,Graph_fy1);
graph->SetName("Graph");
graph->SetTitle("");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffcc00");
graph->SetFillColor(ci);
ci = TColor::GetColor("#ff6666");
graph->SetLineColor(ci);
graph->SetLineWidth(3);
TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","",100,0.9,1099.9);
Graph_Graph1->SetMinimum(228.5697);
Graph_Graph1->SetMaximum(1191.057);
Graph_Graph1->SetDirectory(0);
Graph_Graph1->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1->SetLineColor(ci);
Graph_Graph1->GetXaxis()->SetLabelFont(42);
Graph_Graph1->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1->GetXaxis()->SetTitleFont(42);
Graph_Graph1->GetYaxis()->SetLabelFont(42);
Graph_Graph1->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1->GetYaxis()->SetTitleFont(42);
Graph_Graph1->GetZaxis()->SetLabelFont(42);
Graph_Graph1->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph1);
multigraph->Add(graph,"");
multigraph->Draw("af");
multigraph->GetXaxis()->SetTitle("M_{#chi} (GeV)");
multigraph->GetXaxis()->SetRange(0,96);
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetTitleSize(0.05);
multigraph->GetXaxis()->SetTitleOffset(1.1);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("#Lambda (GeV)");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetTitleSize(0.05);
multigraph->GetYaxis()->SetTitleOffset(1.05);
multigraph->GetYaxis()->SetTitleFont(42);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t Graph_fx2[6] = {
1,
10,
100,
400,
700,
1000};
Double_t Graph_fy2[6] = {
1092.23,
1110.85,
1079.16,
826.374,
584.437,
382.832};
graph = new TGraph(6,Graph_fx2,Graph_fy2);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff6666");
graph->SetLineColor(ci);
TH1F *Graph_Graph2 = new TH1F("Graph_Graph2","Graph",100,0.9,1099.9);
Graph_Graph2->SetMinimum(310.0302);
Graph_Graph2->SetMaximum(1183.652);
Graph_Graph2->SetDirectory(0);
Graph_Graph2->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph2->SetLineColor(ci);
Graph_Graph2->GetXaxis()->SetLabelFont(42);
Graph_Graph2->GetXaxis()->SetLabelSize(0.035);
Graph_Graph2->GetXaxis()->SetTitleSize(0.035);
Graph_Graph2->GetXaxis()->SetTitleFont(42);
Graph_Graph2->GetYaxis()->SetLabelFont(42);
Graph_Graph2->GetYaxis()->SetLabelSize(0.035);
Graph_Graph2->GetYaxis()->SetTitleSize(0.035);
Graph_Graph2->GetYaxis()->SetTitleFont(42);
Graph_Graph2->GetZaxis()->SetLabelFont(42);
Graph_Graph2->GetZaxis()->SetLabelSize(0.035);
Graph_Graph2->GetZaxis()->SetTitleSize(0.035);
Graph_Graph2->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph2);
multigraph->Add(graph,"");
Double_t Graph_fx3[6] = {
1,
10,
100,
400,
700,
1000};
Double_t Graph_fy3[6] = {
980.038,
996.742,
968.943,
742.067,
524.44,
343.371};
graph = new TGraph(6,Graph_fx3,Graph_fy3);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff6666");
graph->SetLineColor(ci);
graph->SetLineWidth(3);
TH1F *Graph_Graph3 = new TH1F("Graph_Graph3","Graph",100,0.9,1099.9);
Graph_Graph3->SetMinimum(278.0339);
Graph_Graph3->SetMaximum(1062.079);
Graph_Graph3->SetDirectory(0);
Graph_Graph3->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph3->SetLineColor(ci);
Graph_Graph3->GetXaxis()->SetLabelFont(42);
Graph_Graph3->GetXaxis()->SetLabelSize(0.035);
Graph_Graph3->GetXaxis()->SetTitleSize(0.035);
Graph_Graph3->GetXaxis()->SetTitleFont(42);
Graph_Graph3->GetYaxis()->SetLabelFont(42);
Graph_Graph3->GetYaxis()->SetLabelSize(0.035);
Graph_Graph3->GetYaxis()->SetTitleSize(0.035);
Graph_Graph3->GetYaxis()->SetTitleFont(42);
Graph_Graph3->GetZaxis()->SetLabelFont(42);
Graph_Graph3->GetZaxis()->SetLabelSize(0.035);
Graph_Graph3->GetZaxis()->SetTitleSize(0.035);
Graph_Graph3->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph3);
multigraph->Add(graph,"");
Double_t Graph_fx4[6] = {
1,
10,
100,
400,
700,
1000};
Double_t Graph_fy4[6] = {
881.525,
896.549,
871.907,
667.806,
471.744,
308.777};
graph = new TGraph(6,Graph_fx4,Graph_fy4);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff6666");
graph->SetLineColor(ci);
TH1F *Graph_Graph4 = new TH1F("Graph_Graph4","Graph",100,0.9,1099.9);
Graph_Graph4->SetMinimum(249.9998);
Graph_Graph4->SetMaximum(955.3262);
Graph_Graph4->SetDirectory(0);
Graph_Graph4->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph4->SetLineColor(ci);
Graph_Graph4->GetXaxis()->SetLabelFont(42);
Graph_Graph4->GetXaxis()->SetLabelSize(0.035);
Graph_Graph4->GetXaxis()->SetTitleSize(0.035);
Graph_Graph4->GetXaxis()->SetTitleFont(42);
Graph_Graph4->GetYaxis()->SetLabelFont(42);
Graph_Graph4->GetYaxis()->SetLabelSize(0.035);
Graph_Graph4->GetYaxis()->SetTitleSize(0.035);
Graph_Graph4->GetYaxis()->SetTitleFont(42);
Graph_Graph4->GetZaxis()->SetLabelFont(42);
Graph_Graph4->GetZaxis()->SetLabelSize(0.035);
Graph_Graph4->GetZaxis()->SetTitleSize(0.035);
Graph_Graph4->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph4);
multigraph->Add(graph,"");
Double_t Graph_fx5[6] = {
1,
10,
100,
400,
700,
1000};
Double_t Graph_fy5[6] = {
1028.6,
1012.32,
1016.7,
751.989,
524.44,
360.443};
graph = new TGraph(6,Graph_fx5,Graph_fy5);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#0000ff");
graph->SetLineColor(ci);
graph->SetLineWidth(3);
TH1F *Graph_Graph5 = new TH1F("Graph_Graph5","Graph",100,0.9,1099.9);
Graph_Graph5->SetMinimum(293.6273);
Graph_Graph5->SetMaximum(1095.416);
Graph_Graph5->SetDirectory(0);
Graph_Graph5->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph5->SetLineColor(ci);
Graph_Graph5->GetXaxis()->SetLabelFont(42);
Graph_Graph5->GetXaxis()->SetLabelSize(0.035);
Graph_Graph5->GetXaxis()->SetTitleSize(0.035);
Graph_Graph5->GetXaxis()->SetTitleFont(42);
Graph_Graph5->GetYaxis()->SetLabelFont(42);
Graph_Graph5->GetYaxis()->SetLabelSize(0.035);
Graph_Graph5->GetYaxis()->SetTitleSize(0.035);
Graph_Graph5->GetYaxis()->SetTitleFont(42);
Graph_Graph5->GetZaxis()->SetLabelFont(42);
Graph_Graph5->GetZaxis()->SetLabelSize(0.035);
Graph_Graph5->GetZaxis()->SetTitleSize(0.035);
Graph_Graph5->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph5);
multigraph->Add(graph,"");
Double_t Graph_fx6[35] = {
3.981,
4.365,
4.786,
5.248,
5.754,
6.31,
6.918,
7.586,
8.318,
9.12,
10,
10.965,
12.023,
13.183,
14.454,
15.849,
17.378,
19.055,
20.893,
22.909,
25.119,
31.623,
39.811,
50.119,
63.096,
79.433,
100,
125.893,
158.489,
199.526,
251.189,
316.228,
1000,
3162.278,
10000};
Double_t Graph_fy6[35] = {
19.76533,
32.58338,
47.96959,
63.10164,
79.00995,
95.55754,
111.9766,
129.0118,
146.2376,
163.5577,
179.3658,
194.9584,
209.5672,
223.2999,
236.0683,
247.7404,
258.366,
268.0662,
277.1383,
285.1146,
292.3901,
306.3527,
314.441,
317.1059,
314.8684,
308.9039,
300.2537,
289.8022,
278.2241,
266.1082,
253.7422,
241.4337,
184.6404,
139.2985,
104.656};
graph = new TGraph(35,Graph_fx6,Graph_fy6);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#6666ff");
graph->SetLineColor(ci);
graph->SetLineStyle(6);
graph->SetLineWidth(3);
TH1F *Graph_Graph6 = new TH1F("Graph_Graph6","Graph",100,3.5829,10999.6);
Graph_Graph6->SetMinimum(17.78879);
Graph_Graph6->SetMaximum(346.84);
Graph_Graph6->SetDirectory(0);
Graph_Graph6->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph6->SetLineColor(ci);
Graph_Graph6->GetXaxis()->SetLabelFont(42);
Graph_Graph6->GetXaxis()->SetLabelSize(0.035);
Graph_Graph6->GetXaxis()->SetTitleSize(0.035);
Graph_Graph6->GetXaxis()->SetTitleFont(42);
Graph_Graph6->GetYaxis()->SetLabelFont(42);
Graph_Graph6->GetYaxis()->SetLabelSize(0.035);
Graph_Graph6->GetYaxis()->SetTitleSize(0.035);
Graph_Graph6->GetYaxis()->SetTitleFont(42);
Graph_Graph6->GetZaxis()->SetLabelFont(42);
Graph_Graph6->GetZaxis()->SetLabelSize(0.035);
Graph_Graph6->GetZaxis()->SetTitleSize(0.035);
Graph_Graph6->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph6);
multigraph->Add(graph,"");
Double_t Graph_fx7[33] = {
5.050109,
5.570498,
6.092375,
6.79836,
7.629504,
8.933253,
10.64203,
12.05014,
13.68495,
14.93901,
17.1131,
19.54765,
22.0792,
25.65551,
29.39341,
34.84152,
41.41884,
48.68491,
59.2071,
75.34994,
97.26741,
131.0246,
163.9459,
202.8261,
243.2248,
309.5828,
383.0166,
473.8718,
622.2993,
833.6,
1028.448,
1283.339,
2026.806};
Double_t Graph_fy7[33] = {
71.77067,
88.82209,
110.1401,
130.7052,
154.2557,
185.7482,
213.1504,
232.735,
248.5277,
258.1632,
268.4587,
280.1818,
286.5229,
293.3758,
296.7995,
299.6924,
299.6947,
297.0875,
293.2892,
285.6511,
275.8293,
262.4571,
252.7364,
243.3346,
234.001,
222.0121,
212.1495,
202.5032,
189.9106,
177.908,
168.9168,
159.7122,
142.9221};
graph = new TGraph(33,Graph_fx7,Graph_fy7);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff9900");
graph->SetLineColor(ci);
graph->SetLineStyle(2);
graph->SetLineWidth(3);
TH1F *Graph_Graph7 = new TH1F("Graph_Graph7","Graph",100,4.545098,2228.982);
Graph_Graph7->SetMinimum(48.97827);
Graph_Graph7->SetMaximum(322.4871);
Graph_Graph7->SetDirectory(0);
Graph_Graph7->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph7->SetLineColor(ci);
Graph_Graph7->GetXaxis()->SetLabelFont(42);
Graph_Graph7->GetXaxis()->SetLabelSize(0.035);
Graph_Graph7->GetXaxis()->SetTitleSize(0.035);
Graph_Graph7->GetXaxis()->SetTitleFont(42);
Graph_Graph7->GetYaxis()->SetLabelFont(42);
Graph_Graph7->GetYaxis()->SetLabelSize(0.035);
Graph_Graph7->GetYaxis()->SetTitleSize(0.035);
Graph_Graph7->GetYaxis()->SetTitleFont(42);
Graph_Graph7->GetZaxis()->SetLabelFont(42);
Graph_Graph7->GetZaxis()->SetLabelSize(0.035);
Graph_Graph7->GetZaxis()->SetTitleSize(0.035);
Graph_Graph7->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph7);
multigraph->Add(graph,"");
Double_t Graph_fx8[90] = {
18.74425,
19.90368,
20.73627,
21.82644,
23.00752,
24.07546,
25.0827,
25.71454,
26.71212,
27.99356,
28.53122,
29.59491,
30.69819,
31.93543,
33.71249,
35.53626,
36.69926,
39.48566,
41.74415,
44.52103,
47.27444,
49.54151,
51.99283,
54.56543,
56.59926,
59.57426,
61.97591,
64.94736,
67.7638,
72.165,
77.52983,
78.67276,
79.25885,
79.50015,
79.8577,
79.98559,
80.10781,
80.93361,
82.97299,
85.0636,
87.3346,
89.6657,
93.68976,
97.60831,
101.0972,
103.4923,
107.3498,
111.0254,
114.4909,
118.2376,
120.3333,
131.3772,
140.9375,
145.9754,
148.1289,
150.9739,
154.3249,
158.6761,
165.3121,
170.7206,
172.2259,
173.4824,
174.2374,
175.2489,
175.7591,
179.6599,
186.3538,
196.4347,
207.6675,
217.3037,
228.0541,
241.4457,
256.7473,
271.0309,
292.4573,
319.2919,
352.6953,
387.3116,
421.6106,
459.6163,
505.4719,
563.2603,
635.9718,
692.2802,
745.897,
789.6904,
831.1771,
886.4295,
950.901,
983.4362};
Double_t Graph_fy8[90] = {
190.042,
192.3104,
192.9376,
195.179,
197.2136,
198.0799,
199.3876,
199.4802,
201.01,
204.4632,
204.529,
206.7971,
208.8441,
208.9712,
209.8662,
210.2649,
211.3333,
212.7612,
214.3898,
216.0449,
217.1898,
218.8107,
219.1693,
219.524,
220.873,
221.9988,
223.3628,
223.7069,
225.3419,
225.4478,
225.5616,
225.584,
239.8787,
256.2561,
271.2293,
293.7999,
304.8973,
307.0445,
308.5203,
309.6434,
310.7722,
310.8221,
310.5429,
310.2559,
309.9569,
309.996,
311.8566,
313.3588,
314.1336,
315.2763,
316.7668,
317.2607,
317.3559,
317.4016,
319.264,
319.2881,
319.6855,
319.7195,
318.6591,
318.3273,
317.9688,
309.2591,
300.7849,
289.5104,
286.1793,
286.2011,
285.5743,
285.6238,
285.3433,
284.0635,
283.1184,
281.2031,
278.656,
277.7288,
275.8579,
273.0541,
270.9095,
265.6796,
261.7535,
256.9895,
253.7797,
247.4406,
240.1473,
234.4071,
229.3303,
226.7081,
224.3734,
220.2731,
215.7486,
213.7682};
graph = new TGraph(90,Graph_fx8,Graph_fy8);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#cc33ff");
graph->SetLineColor(ci);
graph->SetLineStyle(6);
graph->SetLineWidth(3);
TH1F *Graph_Graph8 = new TH1F("Graph_Graph8","Graph",100,16.86982,1079.905);
Graph_Graph8->SetMinimum(177.0742);
Graph_Graph8->SetMaximum(332.6873);
Graph_Graph8->SetDirectory(0);
Graph_Graph8->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph8->SetLineColor(ci);
Graph_Graph8->GetXaxis()->SetLabelFont(42);
Graph_Graph8->GetXaxis()->SetLabelSize(0.035);
Graph_Graph8->GetXaxis()->SetTitleSize(0.035);
Graph_Graph8->GetXaxis()->SetTitleFont(42);
Graph_Graph8->GetYaxis()->SetLabelFont(42);
Graph_Graph8->GetYaxis()->SetLabelSize(0.035);
Graph_Graph8->GetYaxis()->SetTitleSize(0.035);
Graph_Graph8->GetYaxis()->SetTitleFont(42);
Graph_Graph8->GetZaxis()->SetLabelFont(42);
Graph_Graph8->GetZaxis()->SetLabelSize(0.035);
Graph_Graph8->GetZaxis()->SetTitleSize(0.035);
Graph_Graph8->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph8);
multigraph->Add(graph,"");
Double_t Graph_fx9[5] = {
250,
500,
1000,
3000,
5000};
Double_t Graph_fy9[5] = {
597.2217,
587.5585,
450.4588,
245.3363,
177.6733};
graph = new TGraph(5,Graph_fx9,Graph_fy9);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#3399ff");
graph->SetLineColor(ci);
graph->SetLineStyle(4);
graph->SetLineWidth(3);
TH1F *Graph_Graph9 = new TH1F("Graph_Graph9","Graph",100,225,5475);
Graph_Graph9->SetMinimum(135.7185);
Graph_Graph9->SetMaximum(639.1765);
Graph_Graph9->SetDirectory(0);
Graph_Graph9->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph9->SetLineColor(ci);
Graph_Graph9->GetXaxis()->SetLabelFont(42);
Graph_Graph9->GetXaxis()->SetLabelSize(0.035);
Graph_Graph9->GetXaxis()->SetTitleSize(0.035);
Graph_Graph9->GetXaxis()->SetTitleFont(42);
Graph_Graph9->GetYaxis()->SetLabelFont(42);
Graph_Graph9->GetYaxis()->SetLabelSize(0.035);
Graph_Graph9->GetYaxis()->SetTitleSize(0.035);
Graph_Graph9->GetYaxis()->SetTitleFont(42);
Graph_Graph9->GetZaxis()->SetLabelFont(42);
Graph_Graph9->GetZaxis()->SetLabelSize(0.035);
Graph_Graph9->GetZaxis()->SetTitleSize(0.035);
Graph_Graph9->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph9);
multigraph->Add(graph,"");
Double_t Graph_fx10[518] = {
6.04,
6.04,
6.04,
6.04,
6.05,
6.08,
6.11,
6.12,
6.12,
6.13,
6.16,
6.16,
6.16,
6.19,
6.2,
6.2,
6.2,
6.21,
6.24,
6.27,
6.28,
6.28,
6.28,
6.28,
6.28,
6.29,
6.34,
6.37,
6.37,
6.37,
6.37,
6.37,
6.38,
6.43,
6.45,
6.45,
6.45,
6.45,
6.45,
6.45,
6.46,
6.53,
6.54,
6.54,
6.54,
6.54,
6.54,
6.54,
6.62,
6.62,
6.62,
6.62,
6.62,
6.62,
6.64,
6.69,
6.71,
6.71,
6.71,
6.71,
6.71,
6.72,
6.78,
6.8,
6.8,
6.8,
6.8,
6.8,
6.81,
6.87,
6.89,
6.89,
6.89,
6.89,
6.89,
6.91,
6.94,
6.97,
6.98,
6.98,
7,
7.05,
7.08,
7.09,
7.15,
7.17,
7.19,
7.24,
7.27,
7.28,
7.34,
7.37,
7.37,
7.45,
7.46,
7.47,
7.55,
7.56,
7.57,
7.65,
7.66,
7.67,
7.76,
7.77,
7.77,
7.86,
7.88,
7.95,
7.99,
8.05,
8.1,
8.15,
8.19,
8.24,
8.27,
8.31,
8.38,
8.43,
8.48,
8.52,
8.58,
8.61,
8.65,
8.71,
8.75,
8.81,
8.84,
8.88,
8.93,
8.97,
9,
9.05,
9.12,
9.17,
9.24,
9.29,
9.36,
9.42,
9.49,
9.54,
9.62,
9.67,
9.74,
9.8,
9.87,
9.93,
10,
10.2,
10.3,
10.5,
10.6,
10.8,
10.9,
10.9,
11,
11.2,
11.2,
11.3,
11.5,
11.6,
11.8,
11.9,
12.1,
12.2,
12.4,
12.6,
12.7,
12.9,
13.1,
13.3,
13.5,
13.5,
13.6,
13.8,
14,
14.2,
14.4,
14.5,
14.5,
14.8,
15,
15.2,
15.3,
15.3,
15.6,
15.8,
15.9,
16,
16.2,
16.5,
16.6,
16.9,
17.2,
17.3,
17.6,
17.9,
18,
18.3,
18.6,
18.7,
19,
19.3,
19.4,
19.5,
19.8,
20.6,
21.4,
21.8,
22,
22.2,
22.6,
22.9,
23.2,
23.5,
23.8,
24.1,
24.5,
24.8,
25.1,
25.5,
25.8,
26.1,
26.5,
26.8,
27.2,
27.6,
28.1,
28.7,
29.1,
29.4,
29.8,
30.8,
31.7,
32.4,
32.9,
34.2,
35.3,
35.8,
36.4,
36.9,
37.4,
38.3,
38.9,
39.4,
40.4,
41,
41.6,
42.6,
43.2,
43.8,
44.4,
45,
45.6,
46.2,
46.8,
47.5,
48.1,
48.7,
49.4,
50.1,
50.7,
51.4,
52.1,
53.3,
54.8,
55.7,
56.4,
57.2,
57.9,
58.7,
59.5,
60.3,
61.1,
61.9,
62.7,
63.5,
64.4,
65.3,
66.3,
67.6,
68.8,
69.7,
70.7,
71.6,
72.8,
74.2,
75.5,
76.5,
77.5,
78.8,
80.3,
81.8,
82.8,
84,
85.1,
86.4,
88.2,
89.7,
90.9,
92.1,
93.6,
95.5,
97.1,
98.4,
100,
102,
104,
105,
108,
110,
112,
114,
115,
118,
120,
122,
124,
127,
130,
132,
134,
137,
139,
142,
145,
147,
150,
152,
155,
157,
161,
164,
167,
170,
173,
176,
179,
182,
187,
190,
194,
197,
200,
205,
208,
213,
216,
222,
225,
231,
234,
240,
244,
253,
262,
267,
272,
278,
283,
289,
295,
301,
307,
312,
317,
323,
329,
334,
341,
347,
353,
358,
367,
374,
381,
387,
395,
403,
411,
419,
427,
436,
445,
454,
463,
475,
485,
494,
504,
514,
528,
545,
564,
583,
603,
615,
628,
639,
649,
662,
675,
689,
701,
712,
726,
741,
756,
771,
786,
802,
818,
833,
846,
863,
878,
892,
910,
935,
964,
992,
1030,
1050,
1070,
1090,
1110,
1150,
1190,
1210,
1230,
1260,
1280,
1310,
1340,
1370,
1400,
1430,
1460,
1490,
1520,
1550,
1600,
1650,
1690,
1720,
1750,
1790,
1850,
1910,
1950,
1990,
2030,
2070,
2130,
2170,
2210,
2260,
2300,
2370,
2410,
2460,
2540,
2630,
2680,
2730,
2820,
2920,
3020,
3120,
3190,
3240,
3320,
3370,
3490,
3610,
3690,
3750,
3880,
4020,
4100,
4180,
4270,
4350,
4440,
4530,
4620,
4690,
4810,
4880,
5000,
5080,
5200,
5310,
5420,
5520,
5640,
5750,
5860,
5980,
6100,
6220,
6330,
6440,
6560,
6700,
6830,
7050,
7300,
7450,
7590,
7840,
8110,
8380,
8670,
8850,
8990,
9210,
9390,
9610};
Double_t Graph_fy10[518] = {
10.38676,
10.50231,
10.59957,
10.7015,
10.89364,
11.07546,
11.20647,
11.30821,
11.41342,
11.62109,
11.82995,
11.96215,
12.07818,
12.17911,
12.30743,
12.414,
12.52534,
12.76538,
12.97118,
13.12865,
13.25128,
13.37442,
13.49524,
13.62174,
13.74534,
13.99711,
14.24146,
14.40622,
14.54013,
14.67427,
14.8083,
14.94868,
15.22224,
15.47959,
15.6669,
15.80805,
15.95581,
16.10083,
16.2527,
16.40114,
16.65049,
17.00656,
17.18103,
17.34867,
17.50979,
17.663,
17.82315,
18.09534,
18.48112,
18.68225,
18.85134,
19.02836,
19.21403,
19.38419,
19.72512,
20.07365,
20.32073,
20.51337,
20.68111,
20.89181,
21.0759,
21.47197,
21.82747,
22.10836,
22.30316,
22.50687,
22.72023,
22.88709,
23.32871,
23.73763,
24.00677,
24.23435,
24.45738,
24.68282,
24.91041,
25.36578,
25.78503,
26.09729,
26.34398,
26.58842,
27.06598,
27.53896,
27.86668,
28.37133,
28.86478,
29.20029,
29.71859,
30.24769,
30.60227,
31.15273,
31.68307,
32.06927,
32.5362,
33.23597,
33.60661,
34.0793,
34.82676,
35.19988,
35.69491,
36.46783,
36.87798,
37.44003,
38.20997,
38.65276,
39.20006,
40.07932,
40.8414,
41.55343,
42.4036,
43.14178,
44.02907,
44.76045,
45.25133,
46.05417,
46.59723,
47.07559,
47.84267,
48.80958,
49.64402,
50.16063,
51.06094,
51.65136,
52.17979,
52.98822,
53.58735,
54.51024,
55.15295,
55.7454,
56.61769,
57.29215,
57.83828,
58.59123,
59.58957,
60.25911,
61.2608,
61.9507,
63.01658,
63.72388,
64.8542,
65.56339,
66.64773,
67.45883,
68.57788,
69.38042,
70.52167,
71.44279,
73.10832,
74.86664,
76.70132,
78.48531,
80.39639,
82.27539,
83.96915,
84.81317,
86.19446,
88.02168,
88.90821,
90.11315,
92.12023,
94.32802,
96.50886,
98.84577,
100.8332,
102.8086,
104.7116,
106.793,
108.8109,
110.7898,
112.9515,
115.1677,
117.0762,
118.1521,
119.1161,
120.7506,
122.937,
125.5722,
127.4196,
128.5312,
129.6638,
131.542,
134.0308,
136.1693,
137.3979,
138.4952,
140.691,
142.7079,
144.0108,
145.2215,
146.8251,
149.5702,
150.7827,
152.4733,
155.2857,
156.5944,
158.387,
161.2726,
162.6707,
164.4505,
167.5191,
168.8856,
170.7937,
172.6887,
174.1507,
175.4199,
177.3078,
182.9573,
187.7626,
189.8365,
190.4369,
191.0454,
193.294,
195.0114,
196.8028,
198.6746,
200.6337,
202.4317,
204.8624,
205.4538,
206.0518,
208.4005,
210.5569,
212.4997,
214.5452,
216.3183,
218.547,
220.4993,
222.5493,
224.7033,
226.9325,
229.2592,
231.2259,
233.3437,
236.0739,
236.692,
237.8515,
240.2765,
242.669,
242.8986,
243.7,
244.4396,
245.1237,
245.324,
246.1552,
246.9239,
247.738,
248.6028,
249.34,
250.183,
251.0113,
251.7779,
252.4054,
252.4402,
252.4741,
252.5072,
252.5394,
252.576,
252.6065,
251.8236,
251.565,
251.5977,
251.9175,
252.7618,
252.7922,
253.0667,
254.0334,
254.1442,
254.1704,
254.1995,
254.2243,
254.2519,
254.2788,
254.305,
254.3305,
254.3554,
254.3796,
254.4032,
254.4291,
254.4543,
254.2526,
253.4561,
253.3367,
253.3588,
252.5675,
252.295,
252.1031,
251.2675,
251.223,
250.3929,
250.0624,
249.8782,
249.0786,
249.0366,
248.174,
247.9263,
247.944,
247.6981,
246.9351,
246.8915,
246.9084,
246.0196,
245.5296,
244.7362,
244.6935,
243.8445,
243.6187,
242.4377,
241.8683,
241.295,
240.1783,
239.0772,
238.001,
237.4805,
236.957,
236.4553,
236.9969,
237.012,
234.4527,
232.994,
232.0529,
231.124,
230.6725,
229.7764,
228.8918,
227.5962,
226.7575,
226.3463,
225.53,
224.7235,
224.3345,
223.5489,
222.4,
221.6514,
220.9151,
220.5568,
219.4781,
218.4252,
218.4349,
217.4066,
216.4078,
215.7535,
214.7888,
214.4764,
213.537,
212.6227,
211.7229,
210.5545,
209.6973,
208.5851,
207.767,
206.7044,
205.9222,
204.9054,
204.1581,
202.7048,
200.6119,
199.4946,
198.6241,
197.7734,
196.9394,
195.9209,
195.126,
194.1531,
193.393,
192.6467,
192.2811,
191.5566,
190.8456,
190.4965,
189.8054,
189.1259,
188.7928,
187.9657,
186.9986,
186.3681,
185.748,
185.29,
184.688,
183.8008,
183.0782,
182.0883,
181.3984,
180.4531,
179.6626,
178.8891,
178.132,
176.4223,
175.5989,
174.7941,
174.0077,
173.2387,
171.644,
170.3208,
168.5686,
167.267,
165.4898,
164.7208,
163.9695,
163.8889,
163.1556,
162.5181,
161.6606,
160.976,
160.9779,
160.2332,
159.5059,
158.7245,
158.0995,
157.2843,
156.6211,
155.8431,
155.2095,
155.1483,
154.4668,
153.8004,
153.7418,
153.0319,
152.3957,
150.9962,
149.8166,
148.4291,
146.8649,
146.1595,
145.5618,
144.8424,
144.2273,
143.077,
141.5716,
140.9059,
140.2935,
139.6204,
138.999,
138.3208,
137.7278,
136.5779,
135.9239,
135.2852,
134.6612,
134.0513,
133.455,
132.8717,
131.7421,
130.3946,
129.6185,
129.1138,
128.6188,
127.8935,
126.9576,
125.6149,
124.9697,
124.5488,
123.9302,
123.3267,
122.1621,
121.6,
121.0507,
120.5136,
119.9882,
118.9705,
118.3154,
117.8359,
116.905,
115.7187,
115.1476,
114.5904,
113.7792,
112.4881,
111.6267,
110.4508,
109.8856,
109.4437,
108.7971,
108.3764,
107.6595,
106.4832,
106.0117,
105.4593,
104.657,
103.6329,
103.0581,
102.6571,
102.1086,
101.7257,
101.2013,
100.7625,
100.2622,
99.77413,
99.29779,
98.89833,
98.37802,
97.9966,
97.43805,
97.13459,
96.65859,
96.25151,
95.79658,
95.35219,
94.91791,
94.49337,
94.02686,
93.6217,
93.27425,
93.07861,
92.78892,
92.31587,
91.94598,
91.2278,
90.24254,
89.78968,
89.42747,
88.72436,
87.79184,
87.11485,
86.19591,
85.77034,
85.35504,
84.91885,
84.61402,
84.25532};
graph = new TGraph(518,Graph_fx10,Graph_fy10);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#66ccff");
graph->SetLineColor(ci);
graph->SetLineStyle(5);
graph->SetLineWidth(3);
TH1F *Graph_Graph10 = new TH1F("Graph_Graph10","Graph",518,5.436,10570.4);
Graph_Graph10->SetMinimum(9.348083);
Graph_Graph10->SetMaximum(278.861);
Graph_Graph10->SetDirectory(0);
Graph_Graph10->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph10->SetLineColor(ci);
Graph_Graph10->GetXaxis()->SetLabelFont(42);
Graph_Graph10->GetXaxis()->SetLabelSize(0.035);
Graph_Graph10->GetXaxis()->SetTitleSize(0.035);
Graph_Graph10->GetXaxis()->SetTitleFont(42);
Graph_Graph10->GetYaxis()->SetLabelFont(42);
Graph_Graph10->GetYaxis()->SetLabelSize(0.035);
Graph_Graph10->GetYaxis()->SetTitleSize(0.035);
Graph_Graph10->GetYaxis()->SetTitleFont(42);
Graph_Graph10->GetZaxis()->SetLabelFont(42);
Graph_Graph10->GetZaxis()->SetLabelSize(0.035);
Graph_Graph10->GetZaxis()->SetTitleSize(0.035);
Graph_Graph10->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph10);
multigraph->Add(graph,"");
Double_t Graph_fx11[35] = {
3.98,
4.37,
4.79,
5.25,
5.75,
6.31,
6.92,
7.59,
8.32,
9.12,
10,
10.96,
12.02,
13.18,
14.45,
15.85,
17.38,
19.05,
20.89,
22.91,
25.12,
31.62,
39.81,
50.12,
63.1,
79.43,
100,
125.89,
158.49,
199.53,
251.19,
316.23,
1000,
3162.3,
10000};
Double_t Graph_fy11[35] = {
64.47703,
92.96188,
120.9916,
144.7661,
171.1873,
194.8983,
219.9867,
244.5213,
267.5613,
289.6608,
311.7509,
330.2783,
345.9273,
362.605,
376.5173,
387.5872,
398.0475,
406.9399,
414.3522,
421.5604,
426.3119,
433.1181,
435.0014,
429.503,
421.4079,
411.5566,
398.3919,
383.7713,
369.2768,
352.1727,
335.9504,
319.8066,
245.9105,
185.1913,
139.1868};
graph = new TGraph(35,Graph_fx11,Graph_fy11);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#66cc66");
graph->SetLineColor(ci);
graph->SetLineStyle(2);
graph->SetLineWidth(3);
TH1F *Graph_Graph11 = new TH1F("Graph_Graph11","Graph",100,3.582,10999.6);
Graph_Graph11->SetMinimum(27.4246);
Graph_Graph11->SetMaximum(472.0538);
Graph_Graph11->SetDirectory(0);
Graph_Graph11->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph11->SetLineColor(ci);
Graph_Graph11->GetXaxis()->SetLabelFont(42);
Graph_Graph11->GetXaxis()->SetLabelSize(0.035);
Graph_Graph11->GetXaxis()->SetTitleSize(0.035);
Graph_Graph11->GetXaxis()->SetTitleFont(42);
Graph_Graph11->GetYaxis()->SetLabelFont(42);
Graph_Graph11->GetYaxis()->SetLabelSize(0.035);
Graph_Graph11->GetYaxis()->SetTitleSize(0.035);
Graph_Graph11->GetYaxis()->SetTitleFont(42);
Graph_Graph11->GetZaxis()->SetLabelFont(42);
Graph_Graph11->GetZaxis()->SetLabelSize(0.035);
Graph_Graph11->GetZaxis()->SetTitleSize(0.035);
Graph_Graph11->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph11);
multigraph->Add(graph,"");
Double_t Graph_fx12[35] = {
4,
5,
6,
7,
8,
9,
10,
12,
14,
16,
18,
20,
24,
28,
32,
36,
40,
44,
48,
52,
56,
60,
65,
70,
75,
80,
85,
90,
95,
100,
150,
200,
300,
400,
500};
Double_t Graph_fy12[35] = {
91.75045,
115.9467,
132.4599,
144.0692,
152.856,
159.9043,
165.4462,
172.653,
177.2831,
179.5654,
180.8904,
181.5836,
181.5776,
181.2419,
178.737,
176.8406,
175.2212,
173.4303,
171.2294,
169.3498,
167.5194,
165.3283,
163.1921,
161.2337,
159.198,
157.2876,
155.7156,
153.813,
152.3262,
150.9621,
138.5165,
129.9907,
118.4184,
110.7202,
104.9163};
graph = new TGraph(35,Graph_fx12,Graph_fy12);
graph->SetName("Graph");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#333399");
graph->SetLineColor(ci);
graph->SetLineStyle(2);
graph->SetLineWidth(3);
TH1F *Graph_Graph12 = new TH1F("Graph_Graph12","Graph",100,3.6,549.6);
Graph_Graph12->SetMinimum(82.76714);
Graph_Graph12->SetMaximum(190.5669);
Graph_Graph12->SetDirectory(0);
Graph_Graph12->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph12->SetLineColor(ci);
Graph_Graph12->GetXaxis()->SetLabelFont(42);
Graph_Graph12->GetXaxis()->SetLabelSize(0.035);
Graph_Graph12->GetXaxis()->SetTitleSize(0.035);
Graph_Graph12->GetXaxis()->SetTitleFont(42);
Graph_Graph12->GetYaxis()->SetLabelFont(42);
Graph_Graph12->GetYaxis()->SetLabelSize(0.035);
Graph_Graph12->GetYaxis()->SetTitleSize(0.035);
Graph_Graph12->GetYaxis()->SetTitleFont(42);
Graph_Graph12->GetZaxis()->SetLabelFont(42);
Graph_Graph12->GetZaxis()->SetLabelSize(0.035);
Graph_Graph12->GetZaxis()->SetTitleSize(0.035);
Graph_Graph12->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph12);
multigraph->Add(graph,"");
multigraph->Draw("L");
TLegend *leg = new TLegend(0.16,0.61,0.5,0.85,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.03);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("NULL","Razor-0#mu 90% CL limit: AV EFT operator","h");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","Expected limit, with uncertainty","lf");
ci = TColor::GetColor("#ffcc00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
ci = TColor::GetColor("#ff6666");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","Observed limit","l");
ci = TColor::GetColor("#0000ff");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
leg = new TLegend(0.67,0.61,0.92,0.9,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.03);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
entry=leg->AddEntry("Graph","IceCube W^{+}W^{-}","l");
ci = TColor::GetColor("#3399ff");
entry->SetLineColor(ci);
entry->SetLineStyle(4);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","COUPP 2012","l");
ci = TColor::GetColor("#6666ff");
entry->SetLineColor(ci);
entry->SetLineStyle(6);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","SIMPLE 2012","l");
ci = TColor::GetColor("#ff9900");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","Super-K W^{+}W^{-}","l");
ci = TColor::GetColor("#cc33ff");
entry->SetLineColor(ci);
entry->SetLineStyle(6);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","XENON 100","l");
ci = TColor::GetColor("#66ccff");
entry->SetLineColor(ci);
entry->SetLineStyle(5);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","PICO","l");
ci = TColor::GetColor("#66cc66");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","PICASSO","l");
ci = TColor::GetColor("#333399");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(3);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.955,0.945,"18.8 fb^{-1} (8 TeV)");
tex->SetNDC();
tex->SetTextAlign(31);
tex->SetTextFont(42);
tex->SetTextSize(0.045);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.245,0.865,"CMS");
tex->SetNDC();
tex->SetTextAlign(31);
tex->SetTextFont(61);
tex->SetTextSize(0.045);
tex->SetLineWidth(2);
tex->Draw();
TH1F *_copy__1 = new TH1F("_copy__1","",100,0.9,1049.95);
_copy__1->SetMinimum(100);
_copy__1->SetMaximum(2000000);
_copy__1->SetDirectory(0);
_copy__1->SetStats(0);
ci = TColor::GetColor("#000099");
_copy__1->SetLineColor(ci);
_copy__1->GetXaxis()->SetTitle("M_{#chi} (GeV)");
_copy__1->GetXaxis()->SetRange(0,96);
_copy__1->GetXaxis()->SetLabelFont(42);
_copy__1->GetXaxis()->SetTitleSize(0.05);
_copy__1->GetXaxis()->SetTitleOffset(1.1);
_copy__1->GetXaxis()->SetTitleFont(42);
_copy__1->GetYaxis()->SetTitle("#Lambda (GeV)");
_copy__1->GetYaxis()->SetLabelFont(42);
_copy__1->GetYaxis()->SetTitleSize(0.05);
_copy__1->GetYaxis()->SetTitleOffset(1.05);
_copy__1->GetYaxis()->SetTitleFont(42);
_copy__1->GetZaxis()->SetLabelFont(42);
_copy__1->GetZaxis()->SetLabelSize(0.035);
_copy__1->GetZaxis()->SetTitleSize(0.035);
_copy__1->GetZaxis()->SetTitleFont(42);
_copy__1->Draw("sameaxis");
c->Modified();
c->cd();
c->SetSelected(c);
}
void BSA_InCoherent_t()
{
//=========Macro generated from canvas: c17/Incoherent channe
//========= (Tue Aug 21 12:00:26 2018) by ROOT version6.09/01
TCanvas *c17 = new TCanvas("c17", "Incoherent channe",0,0,1200,400);
gStyle->SetOptFit(1);
c17->Range(0,0,1,1);
c17->SetFillColor(0);
c17->SetBorderMode(0);
c17->SetBorderSize(2);
c17->SetLeftMargin(0.025);
c17->SetRightMargin(0.02325581);
c17->SetFrameBorderMode(0);
// ------------>Primitives in pad: c17_1
TPad *c17_1 = new TPad("c17_1", "Incoherent channe_1",0,0,0.262936,0.9);
c17_1->Draw();
c17_1->cd();
c17_1->Range(-32.86883,-0.4888889,353.6809,0.4);
c17_1->SetFillColor(0);
c17_1->SetBorderMode(0);
c17_1->SetBorderSize(2);
c17_1->SetGridx();
c17_1->SetGridy();
c17_1->SetRightMargin(0);
c17_1->SetTopMargin(0);
c17_1->SetFrameBorderMode(0);
c17_1->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.05<-t<0.21 [GeV^{2}/c^{2}]");
Double_t Graph_fx1154[9] = {
21.59954,
59.03178,
96.06391,
141.7703,
179.2242,
222.8541,
259.5793,
303.38,
337.8675};
Double_t Graph_fy1154[9] = {
0.1020231,
0.03114899,
0.2309077,
0.1341245,
0.01684177,
-0.1267998,
-0.2824572,
-0.2260703,
-0.0186068};
Double_t Graph_fex1154[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1154[9] = {
0.04242532,
0.05165487,
0.05483741,
0.07646672,
0.1006626,
0.09089962,
0.05905575,
0.04462244,
0.05521722};
TGraphErrors *gre = new TGraphErrors(9,Graph_fx1154,Graph_fy1154,Graph_fex1154,Graph_fey1154);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1154 = new TH1F("Graph_Graph1154","Graph",100,0,369.4943);
Graph_Graph1154->SetMinimum(-0.4042387);
Graph_Graph1154->SetMaximum(0.3484709);
Graph_Graph1154->SetDirectory(0);
Graph_Graph1154->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1154->SetLineColor(ci);
Graph_Graph1154->GetXaxis()->SetLabelFont(42);
Graph_Graph1154->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetXaxis()->SetTitleFont(42);
Graph_Graph1154->GetYaxis()->SetLabelFont(42);
Graph_Graph1154->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetYaxis()->SetTitleFont(42);
Graph_Graph1154->GetZaxis()->SetLabelFont(42);
Graph_Graph1154->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1154);
TF1 *myfit1155 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,369.4943);
myfit1155->SetFillColor(19);
myfit1155->SetFillStyle(0);
myfit1155->SetLineColor(2);
myfit1155->SetLineWidth(2);
myfit1155->SetChisquare(12.03234);
myfit1155->SetNDF(7);
myfit1155->GetXaxis()->SetLabelFont(42);
myfit1155->GetXaxis()->SetLabelSize(0.035);
myfit1155->GetXaxis()->SetTitleSize(0.035);
myfit1155->GetXaxis()->SetTitleFont(42);
myfit1155->GetYaxis()->SetLabelFont(42);
myfit1155->GetYaxis()->SetLabelSize(0.035);
myfit1155->GetYaxis()->SetTitleSize(0.035);
myfit1155->GetYaxis()->SetTitleFont(42);
myfit1155->SetParameter(0,0.2160608);
myfit1155->SetParError(0,0.02670852);
myfit1155->SetParLimits(0,0,0);
myfit1155->SetParameter(1,0.2378828);
myfit1155->SetParError(1,0.2292802);
myfit1155->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1155);
TPaveStats *ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("#chi^{2} / ndf = 12.03 / 7");
AText = ptstats->AddText("p0 = 0.2161 #pm 0.02671 ");
AText = ptstats->AddText("p1 = 0.2379 #pm 0.2293 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("A_{LU}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
TPaveText *pt = new TPaveText(0.2292451,0.9268996,0.7707549,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.05<-t<0.21 [GeV^{2}/c^{2}]");
pt->Draw();
c17_1->Modified();
c17->cd();
// ------------>Primitives in pad: c17_2
TPad *c17_2 = new TPad("c17_2", "Incoherent channe_2",0.262936,0,0.5008721,0.9);
c17_2->Draw();
c17_2->cd();
c17_2->Range(5.895345,-0.4888889,352.6198,0.4);
c17_2->SetFillColor(0);
c17_2->SetBorderMode(0);
c17_2->SetBorderSize(2);
c17_2->SetGridx();
c17_2->SetGridy();
c17_2->SetLeftMargin(0);
c17_2->SetRightMargin(0);
c17_2->SetTopMargin(0);
c17_2->SetFrameBorderMode(0);
c17_2->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.21<-t<0.36 [GeV^{2}/c^{2}]");
Double_t Graph_fx1157[9] = {
21.65555,
61.41754,
96.1904,
139.625,
182.9173,
222.6445,
259.1526,
302.149,
336.8596};
Double_t Graph_fy1157[9] = {
0.04372413,
0.1019325,
0.1737435,
0.1708402,
0.2316287,
-0.1435087,
-0.1227365,
-0.1969296,
-0.07714485};
Double_t Graph_fex1157[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1157[9] = {
0.0517677,
0.05885852,
0.0646232,
0.09789667,
0.1185262,
0.09557657,
0.06953456,
0.04921871,
0.06282415};
gre = new TGraphErrors(9,Graph_fx1157,Graph_fy1157,Graph_fex1157,Graph_fey1157);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1157 = new TH1F("Graph_Graph1157","Graph",100,0,368.38);
Graph_Graph1157->SetMinimum(-0.3057786);
Graph_Graph1157->SetMaximum(0.4097851);
Graph_Graph1157->SetDirectory(0);
Graph_Graph1157->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1157->SetLineColor(ci);
Graph_Graph1157->GetXaxis()->SetLabelFont(42);
Graph_Graph1157->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetXaxis()->SetTitleFont(42);
Graph_Graph1157->GetYaxis()->SetLabelFont(42);
Graph_Graph1157->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetYaxis()->SetTitleFont(42);
Graph_Graph1157->GetZaxis()->SetLabelFont(42);
Graph_Graph1157->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1157);
TF1 *myfit1158 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,368.38);
myfit1158->SetFillColor(19);
myfit1158->SetFillStyle(0);
myfit1158->SetLineColor(2);
myfit1158->SetLineWidth(2);
myfit1158->SetChisquare(6.955655);
myfit1158->SetNDF(7);
myfit1158->GetXaxis()->SetLabelFont(42);
myfit1158->GetXaxis()->SetLabelSize(0.035);
myfit1158->GetXaxis()->SetTitleSize(0.035);
myfit1158->GetXaxis()->SetTitleFont(42);
myfit1158->GetYaxis()->SetLabelFont(42);
myfit1158->GetYaxis()->SetLabelSize(0.035);
myfit1158->GetYaxis()->SetTitleSize(0.035);
myfit1158->GetYaxis()->SetTitleFont(42);
myfit1158->SetParameter(0,0.1721707);
myfit1158->SetParError(0,0.0324291);
myfit1158->SetParLimits(0,0,0);
myfit1158->SetParameter(1,0.001708824);
myfit1158->SetParError(1,0.4101565);
myfit1158->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1158);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 6.956 / 7");
AText = ptstats->AddText("p0 = 0.1722 #pm 0.03243 ");
AText = ptstats->AddText("p1 = 0.001709 #pm 0.4102 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.21<-t<0.36 [GeV^{2}/c^{2}]");
pt->Draw();
c17_2->Modified();
c17->cd();
// ------------>Primitives in pad: c17_3
TPad *c17_3 = new TPad("c17_3", "Incoherent channe_3",0.5008721,0,0.7388081,0.9);
c17_3->Draw();
c17_3->cd();
c17_3->Range(7.121854,-0.4888889,349.9594,0.4);
c17_3->SetFillColor(0);
c17_3->SetBorderMode(0);
c17_3->SetBorderSize(2);
c17_3->SetGridx();
c17_3->SetGridy();
c17_3->SetLeftMargin(0);
c17_3->SetRightMargin(0);
c17_3->SetTopMargin(0);
c17_3->SetFrameBorderMode(0);
c17_3->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.36<-t<0.67 [GeV^{2}/c^{2}]");
Double_t Graph_fx1160[9] = {
22.70538,
60.97743,
94.23011,
138.1892,
183.0975,
222.79,
259.8951,
302.6808,
334.3759};
Double_t Graph_fy1160[9] = {
0.1761498,
0.1945204,
0.1034809,
-0.09380414,
0.2967375,
-0.1094475,
-0.1583282,
-0.1819771,
-0.09543453};
Double_t Graph_fex1160[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1160[9] = {
0.06358895,
0.0584754,
0.07622601,
0.115077,
0.1268505,
0.1278475,
0.07389284,
0.04701977,
0.07415904};
gre = new TGraphErrors(9,Graph_fx1160,Graph_fy1160,Graph_fex1160,Graph_fey1160);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1160 = new TH1F("Graph_Graph1160","Graph",100,0,365.5429);
Graph_Graph1160->SetMinimum(-0.3033832);
Graph_Graph1160->SetMaximum(0.4896763);
Graph_Graph1160->SetDirectory(0);
Graph_Graph1160->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1160->SetLineColor(ci);
Graph_Graph1160->GetXaxis()->SetLabelFont(42);
Graph_Graph1160->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetXaxis()->SetTitleFont(42);
Graph_Graph1160->GetYaxis()->SetLabelFont(42);
Graph_Graph1160->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetYaxis()->SetTitleFont(42);
Graph_Graph1160->GetZaxis()->SetLabelFont(42);
Graph_Graph1160->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1160);
TF1 *myfit1161 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.5429);
myfit1161->SetFillColor(19);
myfit1161->SetFillStyle(0);
myfit1161->SetLineColor(2);
myfit1161->SetLineWidth(2);
myfit1161->SetChisquare(8.926216);
myfit1161->SetNDF(7);
myfit1161->GetXaxis()->SetLabelFont(42);
myfit1161->GetXaxis()->SetLabelSize(0.035);
myfit1161->GetXaxis()->SetTitleSize(0.035);
myfit1161->GetXaxis()->SetTitleFont(42);
myfit1161->GetYaxis()->SetLabelFont(42);
myfit1161->GetYaxis()->SetLabelSize(0.035);
myfit1161->GetYaxis()->SetTitleSize(0.035);
myfit1161->GetYaxis()->SetTitleFont(42);
myfit1161->SetParameter(0,0.1371093);
myfit1161->SetParError(0,0.03383483);
myfit1161->SetParLimits(0,0,0);
myfit1161->SetParameter(1,-0.6800102);
myfit1161->SetParError(1,0.1832383);
myfit1161->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1161);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.926 / 7");
AText = ptstats->AddText("p0 = 0.1371 #pm 0.03383 ");
AText = ptstats->AddText("p1 = -0.68 #pm 0.1832 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.36<-t<0.67 [GeV^{2}/c^{2}]");
pt->Draw();
c17_3->Modified();
c17->cd();
// ------------>Primitives in pad: c17_4
TPad *c17_4 = new TPad("c17_4", "Incoherent channe_4",0.7388081,0,0.9767442,0.9);
c17_4->Draw();
c17_4->cd();
c17_4->Range(7.951117,-0.4888889,350.3641,0.4);
c17_4->SetFillColor(0);
c17_4->SetBorderMode(0);
c17_4->SetBorderSize(2);
c17_4->SetGridx();
c17_4->SetGridy();
c17_4->SetLeftMargin(0);
c17_4->SetRightMargin(0);
c17_4->SetTopMargin(0);
c17_4->SetFrameBorderMode(0);
c17_4->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.67<-t<0.90 [GeV^{2}/c^{2}]");
Double_t Graph_fx1163[9] = {
23.51534,
59.3759,
90.55341,
131.0631,
176.9476,
229.0094,
260.8643,
304.7922,
334.7999};
Double_t Graph_fy1163[9] = {
0.0123906,
0.195403,
0.18963,
-0.1677207,
-0.8232306,
-0.7760588,
-0.2579708,
-0.2579435,
-0.08614626};
Double_t Graph_fex1163[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1163[9] = {
0.09527872,
0.1168373,
0.1651642,
0.3137499,
0.3909388,
0.4596559,
0.1665094,
0.09339681,
0.09961735};
gre = new TGraphErrors(9,Graph_fx1163,Graph_fy1163,Graph_fex1163,Graph_fey1163);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1163 = new TH1F("Graph_Graph1163","Graph",100,0,365.9283);
Graph_Graph1163->SetMinimum(-1.394766);
Graph_Graph1163->SetMaximum(0.5138451);
Graph_Graph1163->SetDirectory(0);
Graph_Graph1163->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1163->SetLineColor(ci);
Graph_Graph1163->GetXaxis()->SetLabelFont(42);
Graph_Graph1163->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetXaxis()->SetTitleFont(42);
Graph_Graph1163->GetYaxis()->SetLabelFont(42);
Graph_Graph1163->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetYaxis()->SetTitleFont(42);
Graph_Graph1163->GetZaxis()->SetLabelFont(42);
Graph_Graph1163->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1163);
TF1 *myfit1164 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.9283);
myfit1164->SetFillColor(19);
myfit1164->SetFillStyle(0);
myfit1164->SetLineColor(2);
myfit1164->SetLineWidth(2);
myfit1164->SetChisquare(8.810638);
myfit1164->SetNDF(7);
myfit1164->GetXaxis()->SetLabelFont(42);
myfit1164->GetXaxis()->SetLabelSize(0.035);
myfit1164->GetXaxis()->SetTitleSize(0.035);
myfit1164->GetXaxis()->SetTitleFont(42);
myfit1164->GetYaxis()->SetLabelFont(42);
myfit1164->GetYaxis()->SetLabelSize(0.035);
myfit1164->GetYaxis()->SetTitleSize(0.035);
myfit1164->GetYaxis()->SetTitleFont(42);
myfit1164->SetParameter(0,0.2352442);
myfit1164->SetParError(0,0.0807391);
myfit1164->SetParLimits(0,0,0);
myfit1164->SetParameter(1,0.003684647);
myfit1164->SetParError(1,0.5673774);
myfit1164->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1164);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.811 / 7");
AText = ptstats->AddText("p0 = 0.2352 #pm 0.08074 ");
AText = ptstats->AddText("p1 = 0.003685 #pm 0.5674 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.67<-t<0.90 [GeV^{2}/c^{2}]");
pt->Draw();
c17_4->Modified();
c17->cd();
c17->Modified();
c17->cd();
c17->SetSelected(c17);
}
void Pythia8_Dijet30_result_JER_akPu4PF_(abs(refeta)<2)&&(hiBin>60)_refpt()
{
//=========Macro generated from canvas: Can_result_3/
//========= (Mon Apr 18 03:09:33 2016) by ROOT version6.02/13
TCanvas *Can_result_3 = new TCanvas("Can_result_3", "",18,40,700,500);
gStyle->SetOptFit(1);
Can_result_3->Range(-2.301397,-0.01150186,237.5545,0.3828568);
Can_result_3->SetFillColor(0);
Can_result_3->SetBorderMode(0);
Can_result_3->SetBorderSize(2);
Can_result_3->SetFrameBorderMode(0);
Can_result_3->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("name");
multigraph->SetTitle("JER_akPu4PF");
Double_t Graph_fx1007[11] = {
45,
55,
65,
75,
85,
95,
105,
115,
130,
150,
180};
Double_t Graph_fy1007[11] = {
0.1976721,
0.1639212,
0.1562002,
0.1480066,
0.1479035,
0.1333827,
0.1160421,
0.1205967,
0.1105508,
0.08141581,
0.09210746};
Double_t Graph_fex1007[11] = {
5,
5,
5,
5,
5,
5,
5,
5,
10,
10,
20};
Double_t Graph_fey1007[11] = {
0.001003849,
0.001317066,
0.001658995,
0.002251738,
0.003095601,
0.004152981,
0.005688263,
0.008277583,
0.006320737,
0.01258781,
0.01599127};
TGraphErrors *gre = new TGraphErrors(11,Graph_fx1007,Graph_fy1007,Graph_fex1007,Graph_fey1007);
gre->SetName("Graph");
gre->SetTitle("some title_0");
gre->SetFillColor(1);
gre->SetFillStyle(0);
TH1F *Graph_Graph1007 = new TH1F("Graph_Graph1007","some title_0",100,24,216);
Graph_Graph1007->SetMinimum(0.0558432);
Graph_Graph1007->SetMaximum(0.2116608);
Graph_Graph1007->SetDirectory(0);
Graph_Graph1007->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
Graph_Graph1007->SetLineColor(ci);
Graph_Graph1007->GetXaxis()->SetLabelFont(42);
Graph_Graph1007->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetXaxis()->SetTitleFont(42);
Graph_Graph1007->GetYaxis()->SetLabelFont(42);
Graph_Graph1007->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetYaxis()->SetTitleFont(42);
Graph_Graph1007->GetZaxis()->SetLabelFont(42);
Graph_Graph1007->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1007);
multigraph->Add(gre,"");
Double_t Graph_fx1008[10] = {
33.60148,
40.66353,
49.20982,
59.5523,
72.06846,
87.21514,
105.5452,
127.7278,
154.5724,
187.059};
Double_t Graph_fy1008[10] = {
0.3165438,
0.2200061,
0.1916064,
0.1736762,
0.1550799,
0.1417295,
0.1352239,
0.1191096,
0.1233824,
0.06176231};
Double_t Graph_fex1008[10] = {
3.19525,
3.8668,
4.67949,
5.66299,
6.85317,
8.29351,
10.03657,
12.14595,
14.6987,
17.7879};
Double_t Graph_fey1008[10] = {
0.01253682,
0.005719637,
0.004019816,
0.005043033,
0.006710735,
0.008656837,
0.01379847,
0.02334296,
0.03556038,
0.01948799};
gre = new TGraphErrors(10,Graph_fx1008,Graph_fy1008,Graph_fex1008,Graph_fey1008);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
gre->SetFillStyle(0);
gre->SetLineColor(2);
gre->SetMarkerColor(2);
TH1F *Graph_Graph1008 = new TH1F("Graph_Graph1008","Graph",100,12.96216,222.291);
Graph_Graph1008->SetMinimum(0.01359369);
Graph_Graph1008->SetMaximum(0.3577612);
Graph_Graph1008->SetDirectory(0);
Graph_Graph1008->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1008->SetLineColor(ci);
Graph_Graph1008->GetXaxis()->SetLabelFont(42);
Graph_Graph1008->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetXaxis()->SetTitleFont(42);
Graph_Graph1008->GetYaxis()->SetLabelFont(42);
Graph_Graph1008->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetYaxis()->SetTitleFont(42);
Graph_Graph1008->GetZaxis()->SetLabelFont(42);
Graph_Graph1008->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1008);
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("refpt");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("#sigma_{Reco./Gen.}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLegend *leg = new TLegend(0.65,0.7,0.85,0.85,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("Graph","allJets","l");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","from chris","l");
entry->SetLineColor(2);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.35,0.8,"#bf{#color[2]{|#eta_{jet}|<2.0}}");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(26);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.17,0.25,"Centrality:30-100%");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(26);
tex->SetLineWidth(2);
tex->Draw();
TPaveText *pt = new TPaveText(0.3743391,0.9365254,0.6256609,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("JER_akPu4PF");
pt->Draw();
Can_result_3->Modified();
Can_result_3->cd();
Can_result_3->SetSelected(Can_result_3);
}
void Pythia8_Dijet30_result_JER_akPu4PF_((abs(refeta)<2)&&(hiBin>60))&&(refpt>30)_refeta()
{
//=========Macro generated from canvas: Can_result_3/
//========= (Mon Apr 18 04:08:04 2016) by ROOT version6.02/13
TCanvas *Can_result_3 = new TCanvas("Can_result_3", "",18,40,700,500);
gStyle->SetOptFit(1);
Can_result_3->Range(-2.75,0.1715425,2.75,0.2926317);
Can_result_3->SetFillColor(0);
Can_result_3->SetBorderMode(0);
Can_result_3->SetBorderSize(2);
Can_result_3->SetFrameBorderMode(0);
Can_result_3->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("name");
multigraph->SetTitle("JER_akPu4PF");
Double_t Graph_fx1007[9] = {
-1.8,
-1.4,
-1,
-0.6,
-0.2,
0.2,
0.6,
1,
1.6};
Double_t Graph_fy1007[9] = {
0.2213775,
0.2313726,
0.212349,
0.1992318,
0.2015521,
0.2047929,
0.2065286,
0.2141753,
0.2298097};
Double_t Graph_fex1007[9] = {
0.2,
0.2,
0.2,
0.2,
0.2,
0.2,
0.2,
0.2,
0.4};
Double_t Graph_fey1007[9] = {
0.001529502,
0.001378741,
0.001274153,
0.001286487,
0.00132694,
0.001354738,
0.001397823,
0.001305814,
0.001034891};
TGraphErrors *gre = new TGraphErrors(9,Graph_fx1007,Graph_fy1007,Graph_fex1007,Graph_fey1007);
gre->SetName("Graph");
gre->SetTitle("some title_0");
gre->SetFillColor(1);
gre->SetFillStyle(0);
TH1F *Graph_Graph1007 = new TH1F("Graph_Graph1007","some title_0",100,-2.4,2.4);
Graph_Graph1007->SetMinimum(0.1944648);
Graph_Graph1007->SetMaximum(0.236232);
Graph_Graph1007->SetDirectory(0);
Graph_Graph1007->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
Graph_Graph1007->SetLineColor(ci);
Graph_Graph1007->GetXaxis()->SetLabelFont(42);
Graph_Graph1007->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetXaxis()->SetTitleFont(42);
Graph_Graph1007->GetYaxis()->SetLabelFont(42);
Graph_Graph1007->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetYaxis()->SetTitleFont(42);
Graph_Graph1007->GetZaxis()->SetLabelFont(42);
Graph_Graph1007->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1007->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1007->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1007);
multigraph->Add(gre,"");
Double_t Graph_fx1008[20] = {
-1.9,
-1.7,
-1.5,
-1.2,
-0.9,
-0.7,
-0.5,
-0.3,
-0.1,
0.1,
0.3,
0.5,
0.7,
0.9,
1.1,
1.3,
1.5,
1.7,
1.9,
1};
Double_t Graph_fy1008[20] = {
0.2311591,
0.2477474,
0.2655955,
0.2576075,
0.2220231,
0.206768,
0.2064983,
0.2020811,
0.19917,
0.2115871,
0.1952014,
0.1965805,
0.2158946,
0.2003286,
0.213619,
0.213695,
0.2466632,
0.2597152,
0.2553281,
0.2344014};
Double_t Graph_fex1008[20] = {
0.1,
0.1,
0.1,
0.2,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
0.1,
-1};
Double_t Graph_fey1008[20] = {
0.009894431,
0.01000488,
0.01052403,
0.009913006,
0.008527985,
0.007570562,
0.007309829,
0.007320613,
0.00699798,
0.00776784,
0.007146697,
0.00713059,
0.007530691,
0.007364154,
0.007896077,
0.007855752,
0.009496669,
0.01034782,
0.01046585,
0.009252961};
gre = new TGraphErrors(20,Graph_fx1008,Graph_fy1008,Graph_fex1008,Graph_fey1008);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
gre->SetFillStyle(0);
gre->SetLineColor(2);
gre->SetMarkerColor(2);
TH1F *Graph_Graph1008 = new TH1F("Graph_Graph1008","Graph",100,-2.4,2.4);
Graph_Graph1008->SetMinimum(0.1792482);
Graph_Graph1008->SetMaximum(0.284926);
Graph_Graph1008->SetDirectory(0);
Graph_Graph1008->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1008->SetLineColor(ci);
Graph_Graph1008->GetXaxis()->SetLabelFont(42);
Graph_Graph1008->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetXaxis()->SetTitleFont(42);
Graph_Graph1008->GetYaxis()->SetLabelFont(42);
Graph_Graph1008->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetYaxis()->SetTitleFont(42);
Graph_Graph1008->GetZaxis()->SetLabelFont(42);
Graph_Graph1008->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1008->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1008->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1008);
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("refeta");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("#sigma_{Reco./Gen.}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLegend *leg = new TLegend(0.68,0.7,0.88,0.85,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("Graph","allJets","l");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph","from chris","l");
entry->SetLineColor(2);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.35,0.8,"#bf{#color[2]{|#eta_{jet}|<2.0 & refpt>30}}");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(24);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.17,0.17,"Centrality:30-100%");
tex->SetNDC();
tex->SetTextFont(43);
tex->SetTextSize(24);
tex->SetLineWidth(2);
tex->Draw();
TPaveText *pt = new TPaveText(0.3743391,0.9365254,0.6256609,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("JER_akPu4PF");
pt->Draw();
Can_result_3->Modified();
Can_result_3->cd();
Can_result_3->SetSelected(Can_result_3);
}
void EfficiencyByEta_vs_Lz()
{
//=========Macro generated from canvas: eta_lz/eta_lz
//========= (Tue Apr 5 09:30:27 2016) by ROOT version6.02/05
TCanvas *eta_lz = new TCanvas("eta_lz", "eta_lz",0,0,700,500);
eta_lz->Range(-6.5625,-0.125,41.5625,1.125);
eta_lz->SetFillColor(0);
eta_lz->SetBorderMode(0);
eta_lz->SetBorderSize(2);
eta_lz->SetFrameBorderMode(0);
eta_lz->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001[10] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001[10] = {
0.9255743,
0.887184,
0.8789174,
0.8852085,
0.8950226,
0.8479381,
0.7993119,
0.7869822,
0.7105263,
0.5384615};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001[10] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001[10] = {
0.001481033,
0.00433367,
0.006472957,
0.00813138,
0.01004846,
0.01400104,
0.0144885,
0.03622352,
0.06131092,
0.1720041};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001[10] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001[10] = {
0.00145469,
0.00419584,
0.006195321,
0.007675046,
0.009298784,
0.01306018,
0.01377348,
0.03252222,
0.05512412,
0.1648473};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(10,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001);
grae->SetName("divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ");
grae->SetTitle("num_LEp8_1D_A0_LZ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001 = new TH1F("Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001","num_LEp8_1D_A0_LZ",100,0,27.5);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMinimum(0.3104003);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMaximum(0.9830861);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetDirectory(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001);
multigraph->Add(grae,"");
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002[12] = {
0.9411031,
0.9015597,
0.8916002,
0.8833977,
0.8820225,
0.8774584,
0.8594848,
0.8327138,
0.8206522,
0.8528529,
0.8434783,
0.8540146};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002[12] = {
0.001557248,
0.004997005,
0.00715926,
0.009612966,
0.01182494,
0.0141127,
0.01291605,
0.02592262,
0.03284668,
0.02203922,
0.02772573,
0.03663361};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002[12] = {
0.001520079,
0.004786173,
0.006779026,
0.008995481,
0.01092244,
0.01290395,
0.01203332,
0.02325302,
0.02903958,
0.01976026,
0.02447678,
0.03086749};
grae = new TGraphAsymmErrors(12,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002);
grae->SetName("divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ");
grae->SetTitle("num_LEp8_1D_A1_LZ");
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002 = new TH1F("Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002","num_LEp8_1D_A1_LZ",100,0,38.5);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMinimum(0.7723237);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMaximum(0.958105);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetDirectory(0);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003[12] = {
0.9502501,
0.935305,
0.9068493,
0.8950382,
0.9219144,
0.921875,
0.89,
0.8737374,
0.8074866,
0.8093645,
0.8384615,
0.8478261};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003[12] = {
0.003894584,
0.008370582,
0.01203572,
0.01515001,
0.01589315,
0.01802717,
0.01582595,
0.02818923,
0.03323571,
0.02546897,
0.03893712,
0.0470388};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003[12] = {
0.003634742,
0.007516641,
0.01085101,
0.01353658,
0.01357227,
0.01512129,
0.01415643,
0.02402349,
0.02965425,
0.0232519,
0.03313846,
0.03844226};
grae = new TGraphAsymmErrors(12,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003);
grae->SetName("divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ");
grae->SetTitle("num_Gp8_1D_A0_LZ");
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
grae->SetLineStyle(2);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003 = new TH1F("Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003","num_Gp8_1D_A0_LZ",100,0,38.5);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMinimum(0.7562875);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMaximum(0.9718483);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetDirectory(0);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004[12] = {
0.9382671,
0.9028476,
0.8930886,
0.8943534,
0.8732719,
0.8631179,
0.8828829,
0.7391304,
0.8125,
0.7848101,
0.9032258,
0.9473684};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004[12] = {
0.002390155,
0.007511495,
0.01114096,
0.01479634,
0.01803571,
0.02457885,
0.020362,
0.04739985,
0.06185409,
0.05634454,
0.08525764,
0.1108357};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004[12] = {
0.002308279,
0.007043833,
0.01024736,
0.01326161,
0.016194,
0.02157686,
0.01787098,
0.04282533,
0.05064369,
0.04821565,
0.05219903,
0.04358048};
grae = new TGraphAsymmErrors(12,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004);
grae->SetName("divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ");
grae->SetTitle("num_Gp8_1D_A1_LZ");
grae->SetLineStyle(2);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004 = new TH1F("Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004","num_Gp8_1D_A1_LZ",100,0,38.5);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMinimum(0.6618088);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMaximum(1.020871);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetDirectory(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004);
multigraph->Add(grae,"");
multigraph->Draw("ALP");
multigraph->GetXaxis()->SetTitle("L_{Z} [cm]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("Efficiency");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLine *line = new TLine(34.5,0,34.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
line = new TLine(46.5,0,46.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
eta_lz->Modified();
eta_lz->cd();
eta_lz->SetSelected(eta_lz);
}
void EfficiencyByEta_vs_Lz()
{
//=========Macro generated from canvas: eta_lz/eta_lz
//========= (Thu Oct 29 10:19:05 2015) by ROOT version6.02/05
TCanvas *eta_lz = new TCanvas("eta_lz", "eta_lz",1123,69,700,500);
eta_lz->Range(-15,-0.125,95,1.125);
eta_lz->SetFillColor(0);
eta_lz->SetBorderMode(0);
eta_lz->SetBorderSize(2);
eta_lz->SetFrameBorderMode(0);
eta_lz->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001[17] = {
0.8207681,
0.7183406,
0.6823529,
0.6296296,
0.6016949,
0.5540541,
0.4752187,
0.3962264,
0.4306569,
0.4033898,
0.2211982,
0.2352941,
0.1679389,
0.1630435,
0,
0,
0};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001[17] = {
0.01564874,
0.02258889,
0.02724798,
0.03011546,
0.03436919,
0.03578772,
0.02830245,
0.04109258,
0.04522796,
0.02985828,
0.0290375,
0.03208442,
0.03356904,
0.03963993,
0,
0,
0};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001[17] = {
0.01469211,
0.02158727,
0.02611493,
0.02919205,
0.03345254,
0.03526936,
0.02845448,
0.04249806,
0.04632244,
0.03055065,
0.03180918,
0.03514912,
0.03920329,
0.04793882,
0.02750881,
0.02093884,
0.05124115};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(17,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001);
grae->SetName("divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ");
grae->SetTitle("num_LEp8_1D_A0_LZ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001 = new TH1F("Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001","num_LEp8_1D_A0_LZ",100,0,88);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMinimum(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMaximum(0.9190063);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetDirectory(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001);
multigraph->Add(grae,"");
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002[17] = {
0.6232558,
0.5583756,
0.5841584,
0.5437956,
0.496,
0.4319249,
0.3764045,
0.3773585,
0.4468085,
0.3511905,
0.2879377,
0.2727273,
0.2324324,
0.2269504,
0.1729323,
0.139738,
0.1242424};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002[17] = {
0.02003841,
0.02640712,
0.03020659,
0.03202259,
0.03353345,
0.03585729,
0.02667804,
0.04061892,
0.04483449,
0.02702119,
0.02925602,
0.0311239,
0.0321134,
0.03658286,
0.03373231,
0.02335216,
0.01843313};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002[17] = {
0.01964108,
0.02609626,
0.02961994,
0.03168431,
0.03356739,
0.03653866,
0.02740908,
0.04228106,
0.04564915,
0.02795583,
0.03101567,
0.03334003,
0.035247,
0.04082589,
0.03918793,
0.02678903,
0.02089594};
grae = new TGraphAsymmErrors(17,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002);
grae->SetName("divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ");
grae->SetTitle("num_LEp8_1D_A1_LZ");
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002 = new TH1F("Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002","num_LEp8_1D_A1_LZ",100,0,88);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMinimum(0.05210053);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMaximum(0.6966057);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetDirectory(0);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003[17] = {
0.6829268,
0.6,
0.6851852,
0.6,
0.6060606,
0.5588235,
0.5571429,
0.5853659,
0.6071429,
0.5454545,
0.3875,
0.3333333,
0.2727273,
0.4090909,
0.2957746,
0.1875,
0.1058201};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003[17] = {
0.08888679,
0.1079526,
0.07564299,
0.09880854,
0.1021968,
0.1001643,
0.06689406,
0.09028307,
0.1123558,
0.0589631,
0.05875506,
0.06081322,
0.05358554,
0.06599562,
0.05755151,
0.0355912,
0.02268587};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003[17] = {
0.07863764,
0.1002229,
0.06787653,
0.09223778,
0.0947793,
0.09618714,
0.06508215,
0.08553781,
0.1034446,
0.05782793,
0.06186464,
0.06620896,
0.06018711,
0.06906434,
0.06399194,
0.04099953,
0.02733468};
grae = new TGraphAsymmErrors(17,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003);
grae->SetName("divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ");
grae->SetTitle("num_Gp8_1D_A0_LZ");
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
grae->SetLineStyle(2);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003 = new TH1F("Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003","num_Gp8_1D_A0_LZ",100,0,88);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMinimum(0.01529121);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMaximum(0.8294075);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetDirectory(0);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004[17] = {
0.7601626,
0.7016575,
0.6848485,
0.6185567,
0.648,
0.5301205,
0.6028369,
0.5737705,
0.4423077,
0.3730159,
0.3243243,
0.2906977,
0.1764706,
0.3265306,
0.1176471,
0.09677419,
0.03703704};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004[17] = {
0.03033946,
0.03787332,
0.04029006,
0.05551651,
0.04781274,
0.06080294,
0.04537447,
0.07221676,
0.07639601,
0.04576098,
0.04698051,
0.05176228,
0.05528972,
0.07233421,
0.05543815,
0.03771076,
0.02388645};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004[17] = {
0.02807219,
0.03547253,
0.03787114,
0.05283465,
0.04523635,
0.06000368,
0.04379748,
0.06951453,
0.07889241,
0.047951,
0.05044641,
0.05716632,
0.07010552,
0.08038169,
0.08329135,
0.05332885,
0.04677086};
grae = new TGraphAsymmErrors(17,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004);
grae->SetName("divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ");
grae->SetTitle("num_Gp8_1D_A1_LZ");
grae->SetLineStyle(2);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004 = new TH1F("Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004","num_Gp8_1D_A1_LZ",100,0,88);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMinimum(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMaximum(0.8657432);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetDirectory(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004);
multigraph->Add(grae,"");
multigraph->Draw("ALP");
multigraph->GetXaxis()->SetTitle("L_{Z} [cm]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("Efficiency");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLine *line = new TLine(34.5,0,34.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
line = new TLine(46.5,0,46.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
TLegend *leg = new TLegend(0.6982759,0.6701903,0.8979885,0.8794926,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.03171247);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("divide_num_LEp8_1D_A0_LXY_by_den_LEp8_1D_A0_LXY","|#eta_{#gamma D}| #leq 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A0_LXY_by_den_Gp8_1D_A0_LXY","|#eta_{#gamma D}| > 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_LEp8_1D_A1_LXY_by_den_LEp8_1D_A1_LXY","|#eta_{#gamma D}| #leq 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A1_LXY_by_den_Gp8_1D_A1_LXY","|#eta_{#gamma D}| > 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
eta_lz->Modified();
eta_lz->cd();
eta_lz->SetSelected(eta_lz);
}
void BSA_InCoherent_xB()
{
//=========Macro generated from canvas: c16/Incoherent channe
//========= (Tue Aug 21 12:00:25 2018) by ROOT version6.09/01
TCanvas *c16 = new TCanvas("c16", "Incoherent channe",0,0,1200,400);
gStyle->SetOptFit(1);
c16->Range(0,0,1,1);
c16->SetFillColor(0);
c16->SetBorderMode(0);
c16->SetBorderSize(2);
c16->SetLeftMargin(0.025);
c16->SetRightMargin(0.02325581);
c16->SetFrameBorderMode(0);
// ------------>Primitives in pad: c16_1
TPad *c16_1 = new TPad("c16_1", "Incoherent channe_1",0,0,0.262936,0.9);
c16_1->Draw();
c16_1->cd();
c16_1->Range(-29.34641,-0.4888889,350.2581,0.4);
c16_1->SetFillColor(0);
c16_1->SetBorderMode(0);
c16_1->SetBorderSize(2);
c16_1->SetGridx();
c16_1->SetGridy();
c16_1->SetRightMargin(0);
c16_1->SetTopMargin(0);
c16_1->SetFrameBorderMode(0);
c16_1->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.12<x_{B}<0.21 ");
Double_t Graph_fx1134[9] = {
24.14331,
62.70082,
96.14623,
141.7289,
180.8857,
220.5862,
258.7856,
303.3336,
334.7288};
Double_t Graph_fy1134[9] = {
0.1330996,
-0.03883743,
0.1388413,
0.2120595,
0.03749243,
-0.09129921,
-0.1273289,
-0.1643618,
0.03103783};
Double_t Graph_fex1134[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1134[9] = {
0.07077824,
0.06445297,
0.05939535,
0.06833722,
0.08559024,
0.08401363,
0.06136237,
0.05845869,
0.08718727};
TGraphErrors *gre = new TGraphErrors(9,Graph_fx1134,Graph_fy1134,Graph_fex1134,Graph_fey1134);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1134 = new TH1F("Graph_Graph1134","Graph",100,0,365.7873);
Graph_Graph1134->SetMinimum(-0.2731422);
Graph_Graph1134->SetMaximum(0.3307184);
Graph_Graph1134->SetDirectory(0);
Graph_Graph1134->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1134->SetLineColor(ci);
Graph_Graph1134->GetXaxis()->SetLabelFont(42);
Graph_Graph1134->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetXaxis()->SetTitleFont(42);
Graph_Graph1134->GetYaxis()->SetLabelFont(42);
Graph_Graph1134->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetYaxis()->SetTitleFont(42);
Graph_Graph1134->GetZaxis()->SetLabelFont(42);
Graph_Graph1134->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1134);
TF1 *myfit1135 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.7873);
myfit1135->SetFillColor(19);
myfit1135->SetFillStyle(0);
myfit1135->SetLineColor(2);
myfit1135->SetLineWidth(2);
myfit1135->SetChisquare(10.13595);
myfit1135->SetNDF(7);
myfit1135->GetXaxis()->SetLabelFont(42);
myfit1135->GetXaxis()->SetLabelSize(0.035);
myfit1135->GetXaxis()->SetTitleSize(0.035);
myfit1135->GetXaxis()->SetTitleFont(42);
myfit1135->GetYaxis()->SetLabelFont(42);
myfit1135->GetYaxis()->SetLabelSize(0.035);
myfit1135->GetYaxis()->SetTitleSize(0.035);
myfit1135->GetYaxis()->SetTitleFont(42);
myfit1135->SetParameter(0,0.1205343);
myfit1135->SetParError(0,0.03387604);
myfit1135->SetParLimits(0,0,0);
myfit1135->SetParameter(1,0.6858633);
myfit1135->SetParError(1,0.2646153);
myfit1135->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1135);
TPaveStats *ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("#chi^{2} / ndf = 10.14 / 7");
AText = ptstats->AddText("p0 = 0.1205 #pm 0.03388 ");
AText = ptstats->AddText("p1 = 0.6859 #pm 0.2646 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("A_{LU}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_1->Modified();
c16->cd();
// ------------>Primitives in pad: c16_2
TPad *c16_2 = new TPad("c16_2", "Incoherent channe_2",0.262936,0,0.5008721,0.9);
c16_2->Draw();
c16_2->cd();
c16_2->Range(5.85529,-0.4888889,352.0975,0.4);
c16_2->SetFillColor(0);
c16_2->SetBorderMode(0);
c16_2->SetBorderSize(2);
c16_2->SetGridx();
c16_2->SetGridy();
c16_2->SetLeftMargin(0);
c16_2->SetRightMargin(0);
c16_2->SetTopMargin(0);
c16_2->SetFrameBorderMode(0);
c16_2->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.21<x_{B}<0.31 ");
Double_t Graph_fx1137[9] = {
21.59357,
60.08662,
95.00395,
140.4897,
181.0443,
224.9978,
260.1264,
302.891,
336.3592};
Double_t Graph_fy1137[9] = {
0.0984356,
0.1891659,
0.2419886,
0.08897874,
0.2077051,
-0.3071227,
-0.1933684,
-0.2875592,
-0.08265727};
Double_t Graph_fex1137[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1137[9] = {
0.04446748,
0.058618,
0.06898028,
0.1082398,
0.135301,
0.1152745,
0.0716065,
0.04680018,
0.05547391};
gre = new TGraphErrors(9,Graph_fx1137,Graph_fy1137,Graph_fex1137,Graph_fey1137);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1137 = new TH1F("Graph_Graph1137","Graph",100,0,367.8358);
Graph_Graph1137->SetMinimum(-0.4989376);
Graph_Graph1137->SetMaximum(0.4195464);
Graph_Graph1137->SetDirectory(0);
Graph_Graph1137->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1137->SetLineColor(ci);
Graph_Graph1137->GetXaxis()->SetLabelFont(42);
Graph_Graph1137->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetXaxis()->SetTitleFont(42);
Graph_Graph1137->GetYaxis()->SetLabelFont(42);
Graph_Graph1137->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetYaxis()->SetTitleFont(42);
Graph_Graph1137->GetZaxis()->SetLabelFont(42);
Graph_Graph1137->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1137);
TF1 *myfit1138 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,367.8358);
myfit1138->SetFillColor(19);
myfit1138->SetFillStyle(0);
myfit1138->SetLineColor(2);
myfit1138->SetLineWidth(2);
myfit1138->SetChisquare(7.336971);
myfit1138->SetNDF(7);
myfit1138->GetXaxis()->SetLabelFont(42);
myfit1138->GetXaxis()->SetLabelSize(0.035);
myfit1138->GetXaxis()->SetTitleSize(0.035);
myfit1138->GetXaxis()->SetTitleFont(42);
myfit1138->GetYaxis()->SetLabelFont(42);
myfit1138->GetYaxis()->SetLabelSize(0.035);
myfit1138->GetYaxis()->SetTitleSize(0.035);
myfit1138->GetYaxis()->SetTitleFont(42);
myfit1138->SetParameter(0,0.2508538);
myfit1138->SetParError(0,0.03583101);
myfit1138->SetParLimits(0,0,0);
myfit1138->SetParameter(1,-0.1422748);
myfit1138->SetParError(1,0.2426713);
myfit1138->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1138);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 7.337 / 7");
AText = ptstats->AddText("p0 = 0.2509 #pm 0.03583 ");
AText = ptstats->AddText("p1 = -0.1423 #pm 0.2427 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_2->Modified();
c16->cd();
// ------------>Primitives in pad: c16_3
TPad *c16_3 = new TPad("c16_3", "Incoherent channe_3",0.5008721,0,0.7388081,0.9);
c16_3->Draw();
c16_3->cd();
c16_3->Range(4.159976,-0.4888889,357.6921,0.4);
c16_3->SetFillColor(0);
c16_3->SetBorderMode(0);
c16_3->SetBorderSize(2);
c16_3->SetGridx();
c16_3->SetGridy();
c16_3->SetLeftMargin(0);
c16_3->SetRightMargin(0);
c16_3->SetTopMargin(0);
c16_3->SetFrameBorderMode(0);
c16_3->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.31<x_{B}<0.39 ");
Double_t Graph_fx1140[9] = {
20.22962,
57.37981,
95.86069,
135.3486,
182.9227,
229.5531,
263.4275,
301.2495,
341.6225};
Double_t Graph_fy1140[9] = {
0.09662004,
0.06965272,
0.269744,
-0.170935,
-0.1509905,
0.6309487,
0.06270791,
-0.2582676,
-0.152689};
Double_t Graph_fex1140[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1140[9] = {
0.05297617,
0.09107484,
0.121636,
0.2287053,
0.3066164,
0.3879093,
0.1400804,
0.07623737,
0.06739512};
gre = new TGraphErrors(9,Graph_fx1140,Graph_fy1140,Graph_fex1140,Graph_fey1140);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1140 = new TH1F("Graph_Graph1140","Graph",100,0,373.7618);
Graph_Graph1140->SetMinimum(-0.6052533);
Graph_Graph1140->SetMaximum(1.166504);
Graph_Graph1140->SetDirectory(0);
Graph_Graph1140->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1140->SetLineColor(ci);
Graph_Graph1140->GetXaxis()->SetLabelFont(42);
Graph_Graph1140->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetXaxis()->SetTitleFont(42);
Graph_Graph1140->GetYaxis()->SetLabelFont(42);
Graph_Graph1140->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetYaxis()->SetTitleFont(42);
Graph_Graph1140->GetZaxis()->SetLabelFont(42);
Graph_Graph1140->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1140);
TF1 *myfit1141 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,373.7618);
myfit1141->SetFillColor(19);
myfit1141->SetFillStyle(0);
myfit1141->SetLineColor(2);
myfit1141->SetLineWidth(2);
myfit1141->SetChisquare(10.84482);
myfit1141->SetNDF(7);
myfit1141->GetXaxis()->SetLabelFont(42);
myfit1141->GetXaxis()->SetLabelSize(0.035);
myfit1141->GetXaxis()->SetTitleSize(0.035);
myfit1141->GetXaxis()->SetTitleFont(42);
myfit1141->GetYaxis()->SetLabelFont(42);
myfit1141->GetYaxis()->SetLabelSize(0.035);
myfit1141->GetYaxis()->SetTitleSize(0.035);
myfit1141->GetYaxis()->SetTitleFont(42);
myfit1141->SetParameter(0,0.124049);
myfit1141->SetParError(0,0.05069177);
myfit1141->SetParLimits(0,0,0);
myfit1141->SetParameter(1,-0.6977989);
myfit1141->SetParError(1,0.2049165);
myfit1141->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1141);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 10.84 / 7");
AText = ptstats->AddText("p0 = 0.124 #pm 0.05069 ");
AText = ptstats->AddText("p1 = -0.6978 #pm 0.2049 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_3->Modified();
c16->cd();
// ------------>Primitives in pad: c16_4
TPad *c16_4 = new TPad("c16_4", "Incoherent channe_4",0.7388081,0,0.9767442,0.9);
c16_4->Draw();
c16_4->cd();
c16_4->Range(-17.15326,-0.4888889,360.2184,0.4);
c16_4->SetFillColor(0);
c16_4->SetBorderMode(0);
c16_4->SetBorderSize(2);
c16_4->SetGridx();
c16_4->SetGridy();
c16_4->SetLeftMargin(0);
c16_4->SetRightMargin(0);
c16_4->SetTopMargin(0);
c16_4->SetFrameBorderMode(0);
c16_4->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.39<x_{B}<0.60 ");
Double_t Graph_fx1143[9] = {
19.13311,
55.36321,
90.36243,
127.5684,
0,
232.109,
268.7816,
302.6639,
343.0651};
Double_t Graph_fy1143[9] = {
0.1307559,
0.1996612,
0.006608881,
-2.830736,
0,
-2.034047,
0.1276631,
-0.2950426,
-0.003784223};
Double_t Graph_fex1143[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1143[9] = {
0.05218982,
0.09871893,
0.1820512,
4.457049,
0,
1.502801,
0.207584,
0.0825667,
0.06276702};
gre = new TGraphErrors(9,Graph_fx1143,Graph_fy1143,Graph_fex1143,Graph_fey1143);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1143 = new TH1F("Graph_Graph1143","Graph",100,0,377.3716);
Graph_Graph1143->SetMinimum(-8.179195);
Graph_Graph1143->SetMaximum(2.517722);
Graph_Graph1143->SetDirectory(0);
Graph_Graph1143->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1143->SetLineColor(ci);
Graph_Graph1143->GetXaxis()->SetLabelFont(42);
Graph_Graph1143->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetXaxis()->SetTitleFont(42);
Graph_Graph1143->GetYaxis()->SetLabelFont(42);
Graph_Graph1143->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetYaxis()->SetTitleFont(42);
Graph_Graph1143->GetZaxis()->SetLabelFont(42);
Graph_Graph1143->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1143);
TF1 *myfit1144 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,377.3716);
myfit1144->SetFillColor(19);
myfit1144->SetFillStyle(0);
myfit1144->SetLineColor(2);
myfit1144->SetLineWidth(2);
myfit1144->SetChisquare(8.6938);
myfit1144->SetNDF(7);
myfit1144->GetXaxis()->SetLabelFont(42);
myfit1144->GetXaxis()->SetLabelSize(0.035);
myfit1144->GetXaxis()->SetTitleSize(0.035);
myfit1144->GetXaxis()->SetTitleFont(42);
myfit1144->GetYaxis()->SetLabelFont(42);
myfit1144->GetYaxis()->SetLabelSize(0.035);
myfit1144->GetYaxis()->SetTitleSize(0.035);
myfit1144->GetYaxis()->SetTitleFont(42);
myfit1144->SetParameter(0,0.1663497);
myfit1144->SetParError(0,0.06722025);
myfit1144->SetParLimits(0,0,0);
myfit1144->SetParameter(1,-0.5093054);
myfit1144->SetParError(1,0.2952021);
myfit1144->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1144);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.694 / 7");
AText = ptstats->AddText("p0 = 0.1663 #pm 0.06722 ");
AText = ptstats->AddText("p1 = -0.5093 #pm 0.2952 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_4->Modified();
c16->cd();
c16->Modified();
c16->cd();
c16->SetSelected(c16);
}
int main(int argc, char *argv[])
{
setenv("TZ","UTC",1);
tzset();
gROOT->SetStyle("Plain");
printf("[email protected] (build %s %s)\n", __DATE__, __TIME__);
char *RunFileName="$HOME/analyisis/default_00.root";
char *OutputFileName="$HOME/analysis/RunFlag_out.root";
bool debug = false;
int c;
opterr = 0;
while ((c = getopt (argc, argv, "hf:i:")) != -1) //: after var ->Waits for input
switch (c)
{
case 'f':
OutputFileName = optarg ;
break;
case 'i':
RunFileName = optarg ;
break;
case '?':
if (isprint (optopt))
fprintf (stderr, "Unknown option `-%c'.\nTry '-h'\n", optopt);
else
fprintf (stderr,
"Unknown option character `\\x%x'.\nTry '-h'\n",
optopt);
return 1;
case 'h':
cout << endl;
cout << " Overview: optimal window size\n";
cout << "===============================\n\n";
cout << " -h HELP -- you already tried :-)\n";
cout << " -f output file name (default: " << OutputFileName << ")\n";
cout << " -i input file name in STAR format (default: " << RunFileName << ")\n";
cout << "\n\n\nmailto: [email protected]\n";
exit(1);
break;
default:
abort ();
}
TChain *Tevent = new TChain("Tevent");
CreateTChain(Tevent, RunFileName);
TChain *Theader = new TChain("Theader");
Theader->Add(RunFileName);
TChain *Tchannel_profile = new TChain("Tchannel_profile");
Tchannel_profile->Add(RunFileName);
struct event event;
struct header header;
struct channel_profiles channel_profiles;
TFile *output = new TFile (OutputFileName,"Recreate","RunFlagOut");
CreateTevent(Tevent, &event);
CreateTheader(Theader, &header);
CreateTchannel_profiles(Tchannel_profile, &channel_profiles);
Tevent->GetEntry(0);
Theader->GetEntry(0);
//no of channels per event in this run
int NoChannels = Tchannel_profile->GetEntries();
if(debug) cout << "NoChannels = " << NoChannels << endl;
///////////////////////////////
int START = 0;
unsigned int event_id_start = 1553544;
for(int w=0; w<Tevent->GetEntries(); w++){
Tevent->GetEntry(w);
START = w;
if(event.event_id == event_id_start) break;
if(event.event_id > event_id_start) break;
}
Tevent->GetEntry(0);
if(debug) cout << START << " START" << endl;
int STOP = START+NoChannels*1; //Tevent->GetEntries();
///////////////////////////////
if(STOP> Tevent->GetEntries()){
cerr << "entry doesn't exist" << endl;
STOP = Tevent->GetEntries();
cerr << "changed end number." << endl;
// exit(1);
}
if(debug) cout << "STOP " << STOP << endl;
/***************************************/
int stepsize = int(1*1024);
int startbin = 8*1024;
/***************************************/
int tmp_start=0;
int tmp_stop=0;
//calculate how many window sizes to prozess
Tevent->GetEntry(0);
int kwindow_size = event.window_size / stepsize;
//ZeroPadding
int NoZeros = 7;
int New_window_size = (NoZeros+1)*event.window_size;
float *TraceZP[(STOP-START)*(kwindow_size+1)];
for(int i=0; i<(STOP-START)*(kwindow_size+1); i++) TraceZP[i] = new float [New_window_size];
float TimeFFTZP[New_window_size];
for(int i=0; i<New_window_size; i++) TimeFFTZP[i]=(float)i/(NoZeros+1);
//spectra
float spectrum[(STOP-START)*(kwindow_size+1)][event.window_size/2];
float freq[event.window_size/2];
for(unsigned int i=0; i<event.window_size/2; i++)
freq[i]=40.0+40.0/(event.window_size/2)*2+(40.0/(event.window_size/2.0))*(i);
float phase[event.window_size/2], rawfft[2*event.window_size];
//graph data
float d_ws[kwindow_size+1];
float d_mean[NoChannels][kwindow_size+1];
float d_rms[NoChannels][kwindow_size+1];
float d_RMS_with[NoChannels][kwindow_size+1];
float d_RMS_without[NoChannels][kwindow_size+1];
float tmp1_mean, tmp2_mean, tmp1_rms, tmp2_rms;
unsigned int min_channel_id=999;
//stopwatch
TStopwatch timer[NoChannels][kwindow_size+1];
float real_time[NoChannels][kwindow_size+1];
float cpu_time[NoChannels][kwindow_size+1];
//define output TGraphs and canvases
TGraph **gWin = new TGraph *[(STOP-START)*(kwindow_size+1)];
TCanvas **canWin = new TCanvas *[(STOP-START)*(kwindow_size+1)];
TGraph **gSpec = new TGraph *[(STOP-START)*(kwindow_size+1)];
TCanvas **canSpec = new TCanvas *[(STOP-START)*(kwindow_size+1)];
TGraph **gd_mean = new TGraph *[NoChannels];
TGraph **gd_rms = new TGraph *[NoChannels];
TGraph **gd_RMS_with = new TGraph *[NoChannels];
TGraph **gd_RMS_without = new TGraph *[NoChannels];
TGraph **gtime_cpu = new TGraph *[NoChannels];
TGraph **gtime_real = new TGraph *[NoChannels];
float tmpTrace[event.window_size];
float tmpTraceSpec[event.window_size];
char name[1024];
char title[1024];
int ChId = 0;
int window_size = 0;
float tmpTraceStat[New_window_size];
for(int i=0; i<Tchannel_profile->GetEntries(); i++){
Tchannel_profile->GetEntry(i);
if(channel_profiles.ch_id < min_channel_id)
min_channel_id = channel_profiles.ch_id;
}
cout << "min channel id = " << min_channel_id << endl;
for(int i=0; i<(kwindow_size+1); i++){
printf("\nwithout RFISuppression\t\t\t\t\twith RFISuppression\n");
printf("Channel\t\twindow_size\tmean\trms\t\tChannel\t\twindow_size\tmean\trms\n");
if(i==1) window_size = event.window_size+stepsize;
if(i!=0){
window_size -= stepsize;
New_window_size = (NoZeros+1)*window_size;
for(int c=0; c<New_window_size; c++) TimeFFTZP[c]=(float)c/(NoZeros+1);
for(int c=0; c<window_size/2; c++)
freq[c]=40.0+40.0/(window_size/2)*2+(40.0/(window_size/2.0))*(c);
}
else window_size = event.window_size;
d_ws[i] = window_size;
if(debug){
cout << "kwindow_size = " << kwindow_size << " - i = " << i << endl;
cout << "window_size = " << window_size << endl;
cout << "New_window_size = " << New_window_size << endl;
}
for(int a=START;a<STOP; a++){
Tevent->GetEntry(a);
ChId=0;
Tchannel_profile->GetEntry(ChId);
while( event.channel_profile_id != channel_profiles.ch_id ){
ChId++;
if(Tchannel_profile->GetEntry(ChId)==0){
cerr << "error - Tchannel_profile->GetEntry()" << endl;
exit(0);
}
}
if(i!=0)sprintf(title,"D%i_%s_%s, %s, run=%i, event_id=%i, ws=%i",header.DAQ_id, channel_profiles.antenna_pos,
channel_profiles.polarisation, event.timestamp, header.run_id, event.event_id, window_size);
else sprintf(title,"D%i_%s_%s, %s, run=%i, event_id=%i, ws=%i without RFI suppression",header.DAQ_id, channel_profiles.antenna_pos,
channel_profiles.polarisation, event.timestamp, header.run_id, event.event_id, window_size);
if(debug) cout << title << " - " << a-START+(i)*NoChannels << " " << a << " " << i << endl;
//define readout range
tmp_start = startbin-window_size/2;
if( (tmp_start)<0 ) tmp_start = 0;
tmp_stop = tmp_start + window_size;
if((unsigned int) tmp_stop > event.window_size) tmp_stop = event.window_size;
for(int j=tmp_start; j<tmp_stop; j++) tmpTrace[j-tmp_start] = (event.trace[j] & 0x0fff);
if(debug) cout << "tmp: start:stop = " << tmp_start << ":" << tmp_stop << " - ws=" << tmp_stop-tmp_start << endl;
for(int v=0; v<2; v++){
if(v==1)RFISuppressionMedian(window_size, tmpTrace, tmpTrace);
ZeroPaddingFFT(window_size, tmpTrace, NoZeros, TraceZP[a-START+(i)*NoChannels]);
for(int w=0; w<window_size; w++) tmpTraceStat[w] = Abs(TraceZP[a-START+(i)*NoChannels][w]);
if(v==0){ tmp1_mean = Mean(window_size,tmpTraceStat); tmp1_rms = RMS(window_size,tmpTraceStat);}
if(v==1){ tmp2_mean = Mean(window_size,tmpTraceStat); tmp2_rms = RMS(window_size,tmpTraceStat);}
printf("%s:%s\t%d\t\t%.2f\t%.2f",channel_profiles.antenna_pos, channel_profiles.polarisation, window_size,
Mean(window_size,tmpTraceStat), RMS(window_size,tmpTraceStat));
if(v==0) cout << "\t\t";
}
cout << endl;
d_mean[event.channel_profile_id-min_channel_id][i] = (tmp1_mean-tmp2_mean)/tmp1_mean;
d_rms[event.channel_profile_id-min_channel_id][i] = (tmp1_rms-tmp2_rms)/tmp1_rms;
d_RMS_with[event.channel_profile_id-min_channel_id][i] = tmp2_rms;
d_RMS_without[event.channel_profile_id-min_channel_id][i] = tmp1_rms;
//be sure to have the original data
for(int j=tmp_start; j<tmp_stop; j++) tmpTrace[j-tmp_start] = (event.trace[j] & 0x0fff);
timer[event.channel_profile_id-min_channel_id][i].Start();
if(i!=0)RFISuppressionMedian(window_size, tmpTrace, tmpTrace);
for(int j=0; j<window_size; j++) tmpTraceSpec[j] = tmpTrace[j];
TraceFFT(window_size , tmpTraceSpec , spectrum[a-START+(i)*NoChannels] , phase, rawfft);
ZeroPaddingFFT(window_size, tmpTrace, NoZeros, TraceZP[a-START+(i)*NoChannels]);
timer[event.channel_profile_id-min_channel_id][i].Stop();
sprintf(name,"can%03i_ws%i",a-START+(i)*NoChannels,window_size);
gWin[a-START+(i)*NoChannels] = new TGraph(New_window_size, TimeFFTZP, TraceZP[a-START+(i)*NoChannels]);
gWin[a-START+(i)*NoChannels]->SetTitle(title);
gWin[a-START+(i)*NoChannels]->GetXaxis()->SetTitle(" time #upoint 12.5 / ns");
gWin[a-START+(i)*NoChannels]->GetYaxis()->SetTitle(" ADC value ");
canWin[a-START+(i)*NoChannels] = new TCanvas();
canWin[a-START+(i)*NoChannels]->SetTitle(name);
canWin[a-START+(i)*NoChannels]->SetName(name);
sprintf(name,"spec%03i_ws%i",a-START+(i)*NoChannels,window_size);
gSpec[a-START+(i)*NoChannels] = new TGraph(window_size/2, freq, spectrum[a-START+(i)*NoChannels]);
gSpec[a-START+(i)*NoChannels]->SetTitle(title);
gSpec[a-START+(i)*NoChannels]->GetXaxis()->SetTitle(" frequency / MHz ");
gSpec[a-START+(i)*NoChannels]->GetYaxis()->SetTitle(" Amp ");
canSpec[a-START+(i)*NoChannels] = new TCanvas();
canSpec[a-START+(i)*NoChannels]->SetTitle(name);
canSpec[a-START+(i)*NoChannels]->SetName(name);
} //end over START STOP
}//end for over kwindwo_size
if(debug) cout << "end of loop over kwindow_size" << endl;
TMultiGraph *mg = new TMultiGraph();
mg->SetTitle("optimal window size");
mg->SetName("mg");
TLegend leg = TLegend(0.87,0.7,0.99,0.99);
TMultiGraph *mg3 = new TMultiGraph();
mg3->SetTitle("optimal window size");
mg3->SetName("mg3");
TLegend leg3 = TLegend(0.655779, 0.807692, 0.981156, 0.972028);
for(int i=0; i<NoChannels; i++){
sprintf(name,"gd_mean%02i",i+min_channel_id);
gd_mean[i] = new TGraph(kwindow_size+1, d_ws, d_mean[i]);
gd_mean[i]->SetTitle(name);
gd_mean[i]->SetName(name);
gd_mean[i]->SetLineColor(i+1);
gd_mean[i]->SetLineWidth(2);
sprintf(title,"mean: channel %02i",i+min_channel_id);
mg->Add(gd_mean[i]);
leg.AddEntry(gd_mean[i],title,"L");
gd_mean[i]->Write();
sprintf(name,"gd_rms%02i",i+min_channel_id);
gd_rms[i] = new TGraph(kwindow_size+1, d_ws, d_rms[i]);
gd_rms[i]->SetTitle(name);
gd_rms[i]->SetName(name);
gd_rms[i]->SetLineColor(i+1);
gd_rms[i]->SetLineStyle(2);
gd_rms[i]->SetLineWidth(4);
sprintf(title,"rms: channel %02i",i+min_channel_id);
mg->Add(gd_rms[i]);
leg.AddEntry(gd_rms[i],title,"L");
gd_rms[i]->Write();
sprintf(name,"gd_rms_wiht%02i",i+min_channel_id);
gd_RMS_with[i] = new TGraph(kwindow_size+1, d_ws, d_RMS_with[i]);
gd_RMS_with[i]->SetTitle(name);
gd_RMS_with[i]->SetName(name);
gd_RMS_with[i]->SetLineColor(i+1);
gd_RMS_with[i]->SetLineStyle(2);
gd_RMS_with[i]->SetLineWidth(4);
sprintf(title,"RMS with RFISupp: channel %02i",i+min_channel_id);
mg3->Add(gd_RMS_with[i]);
// leg3.AddEntry(gd_RMS_with[i],title,"L");
gd_RMS_with[i]->Write();
sprintf(name,"gd_rms_wiht%02i",i+min_channel_id);
gd_RMS_without[i] = new TGraph(kwindow_size+1, d_ws, d_RMS_without[i]);
gd_RMS_without[i]->SetTitle(name);
gd_RMS_without[i]->SetName(name);
gd_RMS_without[i]->SetLineColor(i+1);
gd_RMS_without[i]->SetLineWidth(2);
sprintf(title,"RMS without RFISupp: channel %02i",i+min_channel_id);
mg3->Add(gd_RMS_without[i]);
// leg3.AddEntry(gd_RMS_without[i],title,"L");
gd_RMS_without[i]->Write();
if(i==0){
sprintf(title,"without RFI suppression");
leg3.AddEntry(gd_RMS_without[0],title,"L");
sprintf(title,"with RFI suppression");
leg3.AddEntry(gd_RMS_with[0],title,"L");
}
}
mg->Write();
mg3->Write();
TCanvas *gdcan = new TCanvas();
gdcan->SetName("canmg");
gdcan->SetTitle("canmg");
gdcan->cd();
mg->Draw("AL");
leg.Draw("SAME");
gdcan->Write();
TCanvas *gdcan3 = new TCanvas();
gdcan3->SetName("canmg3");
gdcan3->SetTitle("canmg3");
gdcan3->cd();
mg3->Draw("AL");
mg3->GetXaxis()->SetTitle("window size");
mg3->GetYaxis()->SetTitle("RMS / ADC unit");
leg3.Draw("SAME");
gdcan3->Write();
//stopwatch
for(int j=0; j<NoChannels; j++)for(int i=0; i<kwindow_size+1; i++) {
real_time[j][i] = timer[j][i].RealTime();
cpu_time[j][i] = timer[j][i].CpuTime();
// cout << "WS = " << i << " - no channel = " << j << endl;
// cout << "Real = " << timer[j[i].RealTime();
// cout << "CPU = " << timer[j][i].CpuTime() << endl;
}
TMultiGraph *mg2 = new TMultiGraph();
mg2->SetTitle("processing time");
mg2->SetName("mg2");
TLegend leg2 = TLegend(0.655779, 0.807692, 0.981156, 0.972028);
for(int i=0; i<NoChannels; i++){
cpu_time[i][0] = cpu_time[i][1];
sprintf(name,"gtime_cpu%02i",i+min_channel_id);
gtime_cpu[i] = new TGraph(kwindow_size+1, d_ws, cpu_time[i]);
gtime_cpu[i]->SetTitle(name);
gtime_cpu[i]->SetName(name);
gtime_cpu[i]->SetLineColor(i+1);
gtime_cpu[i]->SetLineWidth(4);
gtime_cpu[i]->SetLineStyle(2);
sprintf(title,"cpu time: channel %02i",i+min_channel_id);
mg2->Add(gtime_cpu[i]);
// leg2.AddEntry(gtime_cpu[i],title,"L");
gd_mean[i]->Write();
real_time[i][0] = real_time[i][1];
sprintf(name,"gtime_real%02i",i+min_channel_id);
gtime_real[i] = new TGraph(kwindow_size+1, d_ws, real_time[i]);
gtime_real[i]->SetTitle(name);
gtime_real[i]->SetName(name);
gtime_real[i]->SetLineColor(i+1);
gtime_real[i]->SetLineWidth(2);
sprintf(title,"real time: channel %02i",i+min_channel_id);
mg2->Add(gtime_real[i]);
// leg2.AddEntry(gtime_real[i],title,"L");
gd_rms[i]->Write();
if(i==0){
leg2.AddEntry(gtime_cpu[0],"cpu time","L");
leg2.AddEntry(gtime_real[0],"real time","L");
}
}
mg2->Write();
TCanvas *gtimecan = new TCanvas();
gtimecan->SetName("canmg2");
gtimecan->SetTitle("canmg2");
gtimecan->cd();
mg2->Draw("AL");
mg2->GetYaxis()->SetTitle("t / s");
mg2->GetYaxis()->CenterTitle();
mg2->GetXaxis()->SetTitle("window size");
leg2.Draw("SAME");
gtimecan->Write();
char AntPos[NoChannels][50];
char AntPol[NoChannels][50];
char position[5];
char pol[50];
for(int i=0; i<NoChannels; i++){
if(Tchannel_profile->GetEntry(i)==0){
cerr << "Tchannel_profile->GetEntry() error!" << endl;
exit(0);
}
strcpy(AntPos[i],channel_profiles.antenna_pos);
strcpy(AntPol[i],channel_profiles.polarisation);
strcat(AntPol[i],",");
//cout << AntPos[i] << " " << AntPol[i] << endl;
}
//fill and write the output file
output->cd();
int division=(int)NoChannels/2;
if((int)NoChannels%2 != 0 ) division++;
for(int j=0; j<kwindow_size+1; j++){
canWin[j*NoChannels]->cd();
canWin[j*NoChannels]->Divide(2,division);
canSpec[j*NoChannels]->cd();
canSpec[j*NoChannels]->Divide(2,division);
for(int i=0; i<NoChannels; i++){
sscanf(gWin[i+j*NoChannels]->GetTitle(),"%*c%*i_%c%c%c_%s CH0",&position[0], &position[1], &position[2], pol);
for(int a=0; a<NoChannels; a++){
if(!strcmp(position,AntPos[a]) && !strcmp(pol,AntPol[a])){
//cout << a << " " << AntPos[a] << " " << AntPol[a] << " " << position << " " << pol << endl;
canWin[j*NoChannels]->cd(a+1);
gWin[i+j*NoChannels]->Draw("AL");
//define user range
// gWin[i+j*NoChannels]->GetXaxis()->SetLimits(gWin[i+j*NoChannels]->GetN()/2-stepsize/2,gWin[i+j*NoChannels]->GetN()/2+stepsize/2);
canSpec[j*NoChannels]->cd(a+1);
gPad->SetLogy();
gSpec[i+j*NoChannels]->Draw("AL");
}
}
}
canWin[j*NoChannels]->Write();
canSpec[j*NoChannels]->Write();
}
output->Close();
//do not forget to clean up
delete output;
delete Tevent;
delete Theader;
delete Tchannel_profile;
}
void EfficiencyByEta_vs_Lz()
{
//=========Macro generated from canvas: eta_lz/eta_lz
//========= (Mon Mar 7 08:02:08 2016) by ROOT version6.02/05
TCanvas *eta_lz = new TCanvas("eta_lz", "eta_lz",321,1245,700,500);
eta_lz->Range(-15,-0.125,95,1.125);
eta_lz->SetFillColor(0);
eta_lz->SetBorderMode(0);
eta_lz->SetBorderSize(2);
eta_lz->SetFrameBorderMode(0);
eta_lz->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001[17] = {
0.9080298,
0.8202334,
0.7553908,
0.7090909,
0.670787,
0.6224337,
0.5515311,
0.4741744,
0.4378014,
0.3948586,
0.3029491,
0.2437673,
0.1630435,
0.1233333,
0.05417607,
0.004065041,
0.005617978};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001[17] = {
0.001541312,
0.004453759,
0.006612443,
0.008263867,
0.009814857,
0.01115846,
0.009204997,
0.01492435,
0.01581867,
0.0112827,
0.01209253,
0.01325854,
0.01301159,
0.01358835,
0.01078914,
0.002625212,
0.004647921};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001[17] = {
0.001518644,
0.004370433,
0.006495727,
0.008130245,
0.009671114,
0.01103306,
0.00917045,
0.01496934,
0.01594231,
0.01139329,
0.012364,
0.0137482,
0.01386317,
0.01492055,
0.01301025,
0.005336259,
0.01279935};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(17,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001);
grae->SetName("divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ");
grae->SetTitle("num_LEp8_1D_A0_LZ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001 = new TH1F("Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001","num_LEp8_1D_A0_LZ",100,0,88);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMinimum(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMaximum(1.000406);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetDirectory(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001);
multigraph->Add(grae,"");
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002[17] = {
0.9043252,
0.790366,
0.7452243,
0.6891468,
0.647705,
0.6043902,
0.5508524,
0.5215094,
0.4822203,
0.4753036,
0.3975904,
0.3444882,
0.2953822,
0.2689721,
0.2352321,
0.2118467,
0.1684266};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002[17] = {
0.001634743,
0.004828313,
0.006544022,
0.008413468,
0.009714908,
0.01109307,
0.008688844,
0.01411191,
0.01512672,
0.01023838,
0.01140515,
0.01239205,
0.01309829,
0.01398831,
0.01400723,
0.01092521,
0.008681733};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002[17] = {
0.001610364,
0.004749066,
0.006436997,
0.008292559,
0.00959664,
0.01098896,
0.008658439,
0.01407857,
0.01515845,
0.01025878,
0.01151492,
0.01260155,
0.01343403,
0.01444742,
0.01458683,
0.01133962,
0.009042502};
grae = new TGraphAsymmErrors(17,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002);
grae->SetName("divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ");
grae->SetTitle("num_LEp8_1D_A1_LZ");
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002 = new TH1F("Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002","num_LEp8_1D_A1_LZ",100,0,88);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMinimum(0.08512581);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMaximum(0.9805547);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetDirectory(0);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003[17] = {
0.9137676,
0.8726961,
0.8518519,
0.7936772,
0.7884406,
0.7881706,
0.7514079,
0.7118644,
0.7027027,
0.6429259,
0.5347222,
0.4498607,
0.4194915,
0.4225092,
0.3761468,
0.2380952,
0.1374374};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003[17] = {
0.003377242,
0.007279377,
0.009358388,
0.01233732,
0.01435381,
0.01624815,
0.01286564,
0.02099558,
0.02230088,
0.01548052,
0.01757905,
0.01917595,
0.01912036,
0.02196986,
0.02141066,
0.01425232,
0.009299189};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003[17] = {
0.003263802,
0.006950033,
0.008908972,
0.01183151,
0.01369971,
0.0154201,
0.01244805,
0.02016033,
0.02141713,
0.01519666,
0.01749599,
0.01932073,
0.01935624,
0.0222679,
0.02188472,
0.01484072,
0.009838194};
grae = new TGraphAsymmErrors(17,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003);
grae->SetName("divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ");
grae->SetTitle("num_Gp8_1D_A0_LZ");
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
grae->SetLineStyle(2);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003 = new TH1F("Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003","num_Gp8_1D_A0_LZ",100,0,88);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMinimum(0.04924885);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMaximum(0.9959208);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetDirectory(0);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004[17] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5,
37.5,
42.5,
47.5,
55,
70};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004[17] = {
0.9168931,
0.827886,
0.7737686,
0.7452504,
0.6772031,
0.6503145,
0.6317539,
0.5810526,
0.5179856,
0.4847775,
0.3793103,
0.3253968,
0.2425,
0.2758621,
0.2034483,
0.1534653,
0.08866995};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004[17] = {
0.002468044,
0.007328789,
0.01045828,
0.01345009,
0.01511763,
0.01772701,
0.01442878,
0.02384346,
0.02568126,
0.01765808,
0.01852453,
0.02147849,
0.02197595,
0.0258088,
0.02423641,
0.01824262,
0.01424295};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004[17] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5,
2.5,
2.5,
2.5,
5,
10};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004[17] = {
0.002404119,
0.007095188,
0.01013583,
0.01301246,
0.01476713,
0.0173351,
0.01420331,
0.02348679,
0.0255909,
0.01769492,
0.01886992,
0.0222034,
0.02333645,
0.02729915,
0.02642606,
0.02007379,
0.01644289};
grae = new TGraphAsymmErrors(17,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004);
grae->SetName("divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ");
grae->SetTitle("num_Gp8_1D_A1_LZ");
grae->SetLineStyle(2);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004 = new TH1F("Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004","num_Gp8_1D_A1_LZ",100,0,88);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMinimum(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMaximum(1.003784);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetDirectory(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004);
multigraph->Add(grae,"");
multigraph->Draw("ALP");
multigraph->GetXaxis()->SetTitle("L_{Z} [cm]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("Efficiency");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLine *line = new TLine(34.5,0,34.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
line = new TLine(46.5,0,46.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
TLegend *leg = new TLegend(0.6982759,0.6701903,0.8979885,0.8794926,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.03171247);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("divide_num_LEp8_1D_A0_LXY_by_den_LEp8_1D_A0_LXY","|#eta_{#gamma D}| #leq 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A0_LXY_by_den_Gp8_1D_A0_LXY","|#eta_{#gamma D}| > 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_LEp8_1D_A1_LXY_by_den_LEp8_1D_A1_LXY","|#eta_{#gamma D}| #leq 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A1_LXY_by_den_Gp8_1D_A1_LXY","|#eta_{#gamma D}| > 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
eta_lz->Modified();
eta_lz->cd();
eta_lz->SetSelected(eta_lz);
}
void EfficiencyByEta_vs_Lz()
{
//=========Macro generated from canvas: eta_lz/eta_lz
//========= (Tue Apr 5 07:29:20 2016) by ROOT version6.02/05
TCanvas *eta_lz = new TCanvas("eta_lz", "eta_lz",0,0,700,500);
eta_lz->Range(-6.5625,-0.125,41.5625,1.125);
eta_lz->SetFillColor(0);
eta_lz->SetBorderMode(0);
eta_lz->SetBorderSize(2);
eta_lz->SetFrameBorderMode(0);
eta_lz->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("");
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001[10] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001[10] = {
0.9255811,
0.887184,
0.8789174,
0.8852085,
0.8950226,
0.8479381,
0.7993119,
0.7869822,
0.7105263,
0.5384615};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001[10] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001[10] = {
0.001480903,
0.00433367,
0.006472957,
0.00813138,
0.01004846,
0.01400104,
0.0144885,
0.03622352,
0.06131092,
0.1720041};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001[10] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5};
Double_t divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001[10] = {
0.001454562,
0.00419584,
0.006195321,
0.007675046,
0.009298784,
0.01306018,
0.01377348,
0.03252222,
0.05512412,
0.1648473};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(10,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fy3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_felx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehx3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fely3001,divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ_fehy3001);
grae->SetName("divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ");
grae->SetTitle("num_LEp8_1D_A0_LZ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001 = new TH1F("Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001","num_LEp8_1D_A0_LZ",100,0,27.5);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMinimum(0.3103996);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetMaximum(0.9830935);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetDirectory(0);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A0_LZ_by_den_LEp8_1D_A0_LZ3001);
multigraph->Add(grae,"");
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002[12] = {
0.9399291,
0.9001285,
0.8902202,
0.8792473,
0.8573519,
0.8445122,
0.8166895,
0.7649485,
0.7375,
0.7456576,
0.6967871,
0.7176871};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002[12] = {
0.001405873,
0.004229466,
0.005917511,
0.007916015,
0.01009448,
0.01242246,
0.01068294,
0.02084312,
0.0238508,
0.01626845,
0.02199496,
0.02869559};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002[12] = {
0.001376087,
0.004079244,
0.005657549,
0.007505729,
0.00955166,
0.01169187,
0.01023525,
0.01969351,
0.02259777,
0.01563471,
0.02116716,
0.02712215};
grae = new TGraphAsymmErrors(12,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fy3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_felx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehx3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fely3002,divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ_fehy3002);
grae->SetName("divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ");
grae->SetTitle("num_LEp8_1D_A1_LZ");
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002 = new TH1F("Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002","num_LEp8_1D_A1_LZ",100,0,38.5);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMinimum(0.6481409);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetMaximum(0.9679565);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetDirectory(0);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->SetLineColor(ci);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetXaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetYaxis()->SetTitleFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelFont(42);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_LEp8_1D_A1_LZ_by_den_LEp8_1D_A1_LZ3002);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003[12] = {
0.9226259,
0.8929907,
0.878169,
0.835122,
0.8467742,
0.8478632,
0.8158954,
0.7928389,
0.7811634,
0.735119,
0.7272727,
0.6801802};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003[12] = {
0.003328733,
0.007105987,
0.009308612,
0.01241554,
0.01436289,
0.0163258,
0.01312262,
0.02254821,
0.02394814,
0.01811984,
0.0258094,
0.03434908};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003[12] = {
0.003205148,
0.006725961,
0.008754592,
0.01173451,
0.01338404,
0.01506936,
0.0124614,
0.02095935,
0.02229414,
0.01739016,
0.02444323,
0.032614};
grae = new TGraphAsymmErrors(12,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fy3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_felx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehx3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fely3003,divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ_fehy3003);
grae->SetName("divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ");
grae->SetTitle("num_Gp8_1D_A0_LZ");
ci = TColor::GetColor("#ff0000");
grae->SetLineColor(ci);
grae->SetLineStyle(2);
ci = TColor::GetColor("#ff0000");
grae->SetMarkerColor(ci);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003 = new TH1F("Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003","num_Gp8_1D_A0_LZ",100,0,38.5);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMinimum(0.6178311);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetMaximum(0.953831);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetDirectory(0);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A0_LZ_by_den_Gp8_1D_A0_LZ3003);
multigraph->Add(grae,"");
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004[12] = {
1,
3,
5,
7,
9,
11,
14,
17,
19,
22.5,
27.5,
32.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004[12] = {
0.9385068,
0.9029813,
0.8920233,
0.8859935,
0.8653061,
0.8677419,
0.8693467,
0.7536232,
0.7717391,
0.7142857,
0.7833333,
0.7941176};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004[12] = {
0.002317847,
0.007117628,
0.01056605,
0.01430681,
0.01723106,
0.02212471,
0.0191375,
0.0421842,
0.05225554,
0.0459041,
0.06646625,
0.09317978};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004[12] = {
1,
1,
1,
1,
1,
1,
2,
1,
1,
2.5,
2.5,
2.5};
Double_t divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004[12] = {
0.002240421,
0.00669502,
0.009765471,
0.01297094,
0.0156433,
0.01955835,
0.01714648,
0.03816421,
0.04566168,
0.04218646,
0.05566409,
0.07263947};
grae = new TGraphAsymmErrors(12,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fy3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_felx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehx3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fely3004,divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ_fehy3004);
grae->SetName("divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ");
grae->SetTitle("num_Gp8_1D_A1_LZ");
grae->SetLineStyle(2);
grae->SetMarkerStyle(2);
TH1F *Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004 = new TH1F("Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004","num_Gp8_1D_A1_LZ",100,0,38.5);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMinimum(0.6411451);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetMaximum(0.9679838);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetDirectory(0);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->SetLineColor(ci);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetXaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetYaxis()->SetTitleFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelFont(42);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetLabelSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleSize(0.035);
Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_divide_num_Gp8_1D_A1_LZ_by_den_Gp8_1D_A1_LZ3004);
multigraph->Add(grae,"");
multigraph->Draw("ALP");
multigraph->GetXaxis()->SetTitle("L_{Z} [cm]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("Efficiency");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
TLine *line = new TLine(34.5,0,34.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
line = new TLine(46.5,0,46.5,1);
ci = TColor::GetColor("#0000ff");
line->SetLineColor(ci);
line->SetLineStyle(2);
line->Draw();
TLegend *leg = new TLegend(0.6982759,0.6701903,0.8979885,0.8794926,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.03171247);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("divide_num_LEp8_1D_A0_LXY_by_den_LEp8_1D_A0_LXY","|#eta_{#gamma D}| #leq 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A0_LXY_by_den_Gp8_1D_A0_LXY","|#eta_{#gamma D}| > 0.8 A0","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_LEp8_1D_A1_LXY_by_den_LEp8_1D_A1_LXY","|#eta_{#gamma D}| #leq 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("divide_num_Gp8_1D_A1_LXY_by_den_Gp8_1D_A1_LXY","|#eta_{#gamma D}| > 0.8 A1","L");
entry->SetLineColor(1);
entry->SetLineStyle(2);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
eta_lz->Modified();
eta_lz->cd();
eta_lz->SetSelected(eta_lz);
}