-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot_dm2.C
97 lines (79 loc) · 3.2 KB
/
plot_dm2.C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
void plot_dm2(){
TCanvas *c1 = new TCanvas("c1","c1",800,1000);
c1->SetFillColor(10);
// Neutrino 2016
// double idx[7]={1,7,3,2,4,5,6};
// double value[7]={2.45,2.57,2.42,2.545,2.67,2.50,2.50};
// double valueUncLow[7]={0.08,0.264,0.09,0.082,0.12,0.2,0.2};
// double valueUncHigh[7]={0.08,0.242,0.09,0.084,0.12,0.2,0.19};
//Neutrino 2018
double idx[8]={2,7,8,1,3,5,4,6};
double value[8]={2.47,2.63,2.48,2.434,2.42,2.51,2.50,2.31};
double valueUncLow[8]={0.07,0.14,0.28,0.064,0.085,0.08,0.12,0.13};
double valueUncHigh[8]={0.07,0.14,0.32,0.064,0.085,0.12,0.05,0.11};
TGraphAsymmErrors *oscGraph = new TGraphAsymmErrors(3, value, idx, valueUncLow, valueUncHigh, 0, 0);
oscGraph->SetTitle("");
oscGraph->SetLineColor(2);
oscGraph->SetLineWidth(1);
oscGraph->SetMarkerColor(2);
oscGraph->SetMarkerSize(1.5);
oscGraph->SetMarkerStyle(20);
oscGraph->GetXaxis()->SetLimits(2.1,3.3);
oscGraph->GetYaxis()->SetRangeUser(0,9.);
oscGraph->GetYaxis()->SetNdivisions(0);
oscGraph->GetXaxis()->SetTitle("#Delta m^{2}_{32} (#times 10^{-3} eV^{2})");
oscGraph->SetTitle("Normal Hierarchy Assumed");
oscGraph->Draw("AP");
TLatex latex;
latex.SetTextSize(0.05);
latex.SetTextColor(2);
latex.DrawLatex(2.83,idx[0]-0.1,"Daya Bay");
latex.SetTextColor(2);
latex.DrawLatex(2.83,idx[1]-0.1,"RENO (nGd)");
latex.DrawLatex(2.83,idx[2]-0.1,"RENO (nH)");
Double_t dyb_update = value[3];
Double_t dyb_update_low = valueUncLow[3];
Double_t dyb_update_high = valueUncHigh[3];
Double_t dyb_update_x = -4.0;
TBox box;
box.SetFillColor(17);
box.DrawBox(dyb_update-dyb_update_low, 9.0,
dyb_update+dyb_update_high, 0.0);
oscGraph->Draw("Psame");
TGraphAsymmErrors *oscGraph1 = new TGraphAsymmErrors(3, &value[3], &idx[3], &valueUncLow[3], &valueUncHigh[3], 0, 0);
oscGraph1->SetTitle("");
oscGraph1->SetLineColor(4);
oscGraph1->SetLineWidth(1);
oscGraph1->SetMarkerColor(4);
oscGraph1->SetMarkerSize(1.5);
oscGraph1->SetMarkerStyle(20);
oscGraph1->GetXaxis()->SetLimits(2.2,3.2);
oscGraph1->GetYaxis()->SetRangeUser(0,8.);
oscGraph1->GetYaxis()->SetNdivisions(0);
oscGraph1->GetXaxis()->SetTitle("#Delta m^{2}_{32} (#times 10^{3} eV^{2})");
oscGraph1->SetTitle("Normal Hierarchy Assumed");
oscGraph1->Draw("Psame");
latex.SetTextColor(4);
latex.DrawLatex(2.83,idx[4]-0.1,"MINOS(+)");
latex.SetTextColor(4);
latex.DrawLatex(2.83,idx[3]-0.1,"T2K");
latex.SetTextColor(4);
latex.DrawLatex(2.83,idx[5]-0.1,"NO#nuA");
TGraphAsymmErrors *oscGraph2 = new TGraphAsymmErrors(2, &value[6], &idx[6], &valueUncLow[6], &valueUncHigh[6], 0, 0);
oscGraph2->SetTitle("");
oscGraph2->SetLineColor(8);
oscGraph2->SetLineWidth(1);
oscGraph2->SetMarkerColor(8);
oscGraph2->SetMarkerSize(1.5);
oscGraph2->SetMarkerStyle(20);
oscGraph2->GetXaxis()->SetLimits(2.2,3.2);
oscGraph2->GetYaxis()->SetRangeUser(0,8.);
oscGraph2->GetYaxis()->SetNdivisions(0);
oscGraph2->GetXaxis()->SetTitle("#Delta m^{2}_{32} (#times 10^{3} eV^{2})");
oscGraph2->SetTitle("Normal Hierarchy Assumed");
oscGraph2->Draw("Psame");
latex.SetTextColor(8);
latex.DrawLatex(2.83,idx[6]-0.1,"Super-K");
latex.SetTextColor(8);
latex.DrawLatex(2.83,idx[7]-0.1,"IceCube");
}