diff --git a/add.tex b/add.tex new file mode 100644 index 0000000000000000000000000000000000000000..71861e120cdb5d9d4b1fac27ddec6834e1d0caab --- /dev/null +++ b/add.tex @@ -0,0 +1,47 @@ +\documentclass[12pt,a4paper]{article} + +\usepackage[utf8]{inputenc} % размер + +\usepackage[T2A]{fontenc} % указывает внутреннюю кодировку TeX + +\usepackage[ngerman,english]{babel} + +\usepackage{epsfig} +\usepackage{xcolor} + +\usepackage[labelformat=simple]{subcaption} +\usepackage{array,graphicx,caption} +\usepackage{xcolor,color} + +\begin{document} + +%------------------------------------------------------------------------------- +\begin{figure}[t] + \begin{center} + \includegraphics[width=0.7\textwidth]{figures/single_angular.png} + \end{center} + % + \caption{ + Angular correlation. Red circles - events satisfied the kinematical triangle + } + % + \label{fig:s_angular} +\end{figure} +%------------------------------------------------------------------------------- + +%------------------------------------------------------------------------------- +\begin{figure}[t] + \begin{center} + \includegraphics[width=1.\textwidth]{figures/angular.png} + \end{center} + % + \caption{ + Distributions for different angular cuts. Full statistics, $<40$ deg, $<30$ deg. + } + % + \label{fig:angular_diff} +\end{figure} +%------------------------------------------------------------------------------ + + +\end{document} \ No newline at end of file diff --git a/figures/angular.png b/figures/angular.png new file mode 100644 index 0000000000000000000000000000000000000000..a76e96b673619f958891e9b91e1edd7241bd8ab1 Binary files /dev/null and b/figures/angular.png differ diff --git a/figures/single_angular.png b/figures/single_angular.png new file mode 100644 index 0000000000000000000000000000000000000000..a6d8620b85e703fddc1096e47c811433077f1037 Binary files /dev/null and b/figures/single_angular.png differ diff --git a/root/angular.C b/root/angular.C new file mode 100644 index 0000000000000000000000000000000000000000..cd29cb32149d913b4b4f425551c01a57307967ed --- /dev/null +++ b/root/angular.C @@ -0,0 +1,2019 @@ + void angular() +{ +//=========Macro generated from canvas: c_he7/title +//========= (Mon Oct 30 12:23:13 2023) by ROOT version6.10/08 + TCanvas *c_he7 = new TCanvas("c_he7", "title",67,55,1800,1000); + gStyle->SetOptStat(0); + c_he7->Range(0,0,1,1); + c_he7->SetFillColor(0); + c_he7->SetBorderMode(0); + c_he7->SetBorderSize(2); + c_he7->SetFrameBorderMode(0); + +// ------------>Primitives in pad: c_he7_1 + TPad *c_he7_1 = new TPad("c_he7_1", "c_he7_1",0.01,0.6766667,0.49,0.99); + c_he7_1->Draw(); + c_he7_1->cd(); + c_he7_1->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_1->SetFillColor(0); + c_he7_1->SetBorderMode(0); + c_he7_1->SetBorderSize(2); + c_he7_1->SetFrameBorderMode(0); + c_he7_1->SetFrameBorderMode(0); + + TH2F *htemp__2 = new TH2F("htemp__2","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__2->SetDirectory(0); + htemp__2->SetStats(0); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#000099"); + htemp__2->SetLineColor(ci); + htemp__2->SetMarkerStyle(20); + htemp__2->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__2->GetXaxis()->SetRange(1,40); + htemp__2->GetXaxis()->SetLabelFont(42); + htemp__2->GetXaxis()->SetLabelSize(0.035); + htemp__2->GetXaxis()->SetTitleSize(0.035); + htemp__2->GetXaxis()->SetTitleFont(42); + htemp__2->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__2->GetYaxis()->SetRange(1,40); + htemp__2->GetYaxis()->SetLabelFont(42); + htemp__2->GetYaxis()->SetLabelSize(0.035); + htemp__2->GetYaxis()->SetTitleSize(0.035); + htemp__2->GetYaxis()->SetTitleOffset(0); + htemp__2->GetYaxis()->SetTitleFont(42); + htemp__2->GetZaxis()->SetLabelFont(42); + htemp__2->GetZaxis()->SetLabelSize(0.035); + htemp__2->GetZaxis()->SetTitleSize(0.035); + htemp__2->GetZaxis()->SetTitleFont(42); + htemp__2->Draw(""); + + TPaveText *pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + TText *pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx3[136] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 9.533222, + 19.7312, + 2.250434, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 11.88648, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 1.294231, + 16.38851, + 12.18283, + 3.692083, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 14.27111, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy3[136] = { + 9.163973, + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 4.734602, + 1.949477, + 0.9119229, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 7.070025, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 17.6368, + 5.868251, + 16.82547, + 7.0717, + 17.36037, + 5.493919, + 1.762917, + 8.676415, + 19.98228, + 11.87412, + 20.94647, + 9.686372, + 21.24376, + 14.36828, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 21.01835, + 3.011658, + 3.843498, + 3.301521, + 3.249645, + 9.771441, + 8.041018, + 2.327877, + 19.21397, + 7.417117, + 1.852323, + 19.29377, + 2.539789, + 2.364044, + 15.25004, + 2.678995, + 0.9787855, + 7.585063, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 0.7892501, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 2.623029, + 4.777777, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 7.656126, + 3.515664, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 0.8742778, + 1.442736, + 6.330581, + 4.266314, + 6.523858, + 7.44585, + 10.23886, + 10.57767, + 3.315528, + 5.560736, + 19.32783, + 2.242471, + 1.46184, + 2.058325, + 8.574683, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + TGraph *graph = new TGraph(136,Graph0_fx3,Graph0_fy3); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph3 = new TH1F("Graph_Graph3","Graph",136,-3.580114,28.89493); + Graph_Graph3->SetMinimum(0.03102773); + Graph_Graph3->SetMaximum(27.16594); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + Double_t Graph1_fx4[108] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 19.7312, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 14.27111, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy4[108] = { + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 4.734602, + 1.949477, + 0.9119229, + 5.710012, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 5.868251, + 5.493919, + 1.762917, + 8.676415, + 11.87412, + 9.686372, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 3.301521, + 3.249645, + 8.041018, + 2.327877, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 0.7892501, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 2.623029, + 4.777777, + 1.353373, + 2.366684, + 7.656126, + 3.515664, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 0.8742778, + 1.442736, + 6.330581, + 6.523858, + 7.44585, + 10.23886, + 10.57767, + 3.315528, + 5.560736, + 2.242471, + 1.46184, + 2.058325, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(108,Graph1_fx4,Graph1_fy4); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph4 = new TH1F("Graph_Graph4","Graph",108,0.4340556,28.53001); + Graph_Graph4->SetMinimum(0.03102773); + Graph_Graph4->SetMaximum(16.33892); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + TF1 *fa0 = new TF1("fa","x*3/4",0,30); + fa0->SetFillColor(19); + fa0->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa0->SetLineColor(ci); + fa0->SetLineWidth(2); + fa0->GetXaxis()->SetLabelFont(42); + fa0->GetXaxis()->SetLabelSize(0.035); + fa0->GetXaxis()->SetTitleSize(0.035); + fa0->GetXaxis()->SetTitleFont(42); + fa0->GetYaxis()->SetLabelFont(42); + fa0->GetYaxis()->SetLabelSize(0.035); + fa0->GetYaxis()->SetTitleSize(0.035); + fa0->GetYaxis()->SetTitleOffset(0); + fa0->GetYaxis()->SetTitleFont(42); + fa0->Draw("same"); + + auto legend1 = new TLegend(0.1,0.9,0.9,1.); + legend1->SetHeader("All angles","C"); // option "C" allows to center the heade + legend1->SetTextColor(kRed); + legend1->SetTextSize(0.07); + legend1->Draw("same"); + + + c_he7_1->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_2 + TPad *c_he7_2 = new TPad("c_he7_2", "c_he7_2",0.51,0.6766667,0.99,0.99); + c_he7_2->Draw(); + c_he7_2->cd(); + c_he7_2->Range(-3.75,-1.575,33.75,14.175); + c_he7_2->SetFillColor(0); + c_he7_2->SetBorderMode(0); + c_he7_2->SetBorderSize(2); + c_he7_2->SetFrameBorderMode(0); + c_he7_2->SetFrameBorderMode(0); + + TH1F *h1__3 = new TH1F("h1__3","lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0}",30,0,30); + h1__3->SetBinContent(1,2); + h1__3->SetBinContent(2,1); + h1__3->SetBinContent(3,5); + h1__3->SetBinContent(4,5); + h1__3->SetBinContent(5,5); + h1__3->SetBinContent(6,2); + h1__3->SetBinContent(7,3); + h1__3->SetBinContent(8,4); + h1__3->SetBinContent(9,3); + h1__3->SetBinContent(10,7); + h1__3->SetBinContent(11,3); + h1__3->SetBinContent(12,5); + h1__3->SetBinContent(13,5); + h1__3->SetBinContent(14,6); + h1__3->SetBinContent(15,11); + h1__3->SetBinContent(16,12); + h1__3->SetBinContent(17,7); + h1__3->SetBinContent(18,9); + h1__3->SetBinContent(19,10); + h1__3->SetBinContent(20,5); + h1__3->SetBinContent(21,5); + h1__3->SetBinContent(22,2); + h1__3->SetBinContent(23,3); + h1__3->SetBinContent(24,4); + h1__3->SetBinContent(25,2); + h1__3->SetBinContent(26,9); + h1__3->SetBinContent(27,1); + h1__3->SetEntries(136); + h1__3->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h1__3->SetMarkerColor(ci); + h1__3->SetMarkerStyle(20); + h1__3->GetXaxis()->SetLabelFont(42); + h1__3->GetXaxis()->SetLabelSize(0.035); + h1__3->GetXaxis()->SetTitleSize(0.035); + h1__3->GetXaxis()->SetTitleFont(42); + h1__3->GetYaxis()->SetLabelFont(42); + h1__3->GetYaxis()->SetLabelSize(0.035); + h1__3->GetYaxis()->SetTitleSize(0.035); + h1__3->GetYaxis()->SetTitleOffset(0); + h1__3->GetYaxis()->SetTitleFont(42); + h1__3->GetZaxis()->SetLabelFont(42); + h1__3->GetZaxis()->SetLabelSize(0.035); + h1__3->GetZaxis()->SetTitleSize(0.035); + h1__3->GetZaxis()->SetTitleFont(42); + h1__3->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0}"); + pt->Draw(); + + TH1F *htemp__4 = new TH1F("htemp__4","lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__4->SetBinContent(3,1); + htemp__4->SetBinContent(4,2); + htemp__4->SetBinContent(5,2); + htemp__4->SetBinContent(6,1); + htemp__4->SetBinContent(7,1); + htemp__4->SetBinContent(8,1); + htemp__4->SetBinContent(9,3); + htemp__4->SetBinContent(10,4); + htemp__4->SetBinContent(11,3); + htemp__4->SetBinContent(12,4); + htemp__4->SetBinContent(13,4); + htemp__4->SetBinContent(14,6); + htemp__4->SetBinContent(15,10); + htemp__4->SetBinContent(16,12); + htemp__4->SetBinContent(17,6); + htemp__4->SetBinContent(18,9); + htemp__4->SetBinContent(19,10); + htemp__4->SetBinContent(20,4); + htemp__4->SetBinContent(21,5); + htemp__4->SetBinContent(22,2); + htemp__4->SetBinContent(23,3); + htemp__4->SetBinContent(24,4); + htemp__4->SetBinContent(25,2); + htemp__4->SetBinContent(26,8); + htemp__4->SetBinContent(27,1); + htemp__4->SetEntries(108); + htemp__4->SetDirectory(0); + htemp__4->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__4->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__4->SetMarkerColor(ci); + htemp__4->SetMarkerStyle(20); + htemp__4->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__4->GetXaxis()->SetLabelFont(42); + htemp__4->GetXaxis()->SetLabelSize(0.035); + htemp__4->GetXaxis()->SetTitleSize(0.035); + htemp__4->GetXaxis()->SetTitleFont(42); + htemp__4->GetYaxis()->SetLabelFont(42); + htemp__4->GetYaxis()->SetLabelSize(0.035); + htemp__4->GetYaxis()->SetTitleSize(0.035); + htemp__4->GetYaxis()->SetTitleOffset(0); + htemp__4->GetYaxis()->SetTitleFont(42); + htemp__4->GetZaxis()->SetLabelFont(42); + htemp__4->GetZaxis()->SetLabelSize(0.035); + htemp__4->GetZaxis()->SetTitleSize(0.035); + htemp__4->GetZaxis()->SetTitleFont(42); + htemp__4->Draw("same"); + c_he7_2->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_3 + TPad *c_he7_3 = new TPad("c_he7_3", "c_he7_3",0.01,0.3433333,0.49,0.6566667); + c_he7_3->Draw(); + c_he7_3->cd(); + c_he7_3->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_3->SetFillColor(0); + c_he7_3->SetBorderMode(0); + c_he7_3->SetBorderSize(2); + c_he7_3->SetFrameBorderMode(0); + c_he7_3->SetFrameBorderMode(0); + + TH2F *htemp__5 = new TH2F("htemp__5","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__5->SetDirectory(0); + htemp__5->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__5->SetLineColor(ci); + htemp__5->SetMarkerStyle(20); + htemp__5->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__5->GetXaxis()->SetRange(1,40); + htemp__5->GetXaxis()->SetLabelFont(42); + htemp__5->GetXaxis()->SetLabelSize(0.035); + htemp__5->GetXaxis()->SetTitleSize(0.035); + htemp__5->GetXaxis()->SetTitleFont(42); + htemp__5->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__5->GetYaxis()->SetRange(1,40); + htemp__5->GetYaxis()->SetLabelFont(42); + htemp__5->GetYaxis()->SetLabelSize(0.035); + htemp__5->GetYaxis()->SetTitleSize(0.035); + htemp__5->GetYaxis()->SetTitleOffset(0); + htemp__5->GetYaxis()->SetTitleFont(42); + htemp__5->GetZaxis()->SetLabelFont(42); + htemp__5->GetZaxis()->SetLabelSize(0.035); + htemp__5->GetZaxis()->SetTitleSize(0.035); + htemp__5->GetZaxis()->SetTitleFont(42); + htemp__5->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx5[124] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy5[124] = { + 9.163973, + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 4.734602, + 1.949477, + 0.9119229, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 7.070025, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 7.0717, + 17.36037, + 5.493919, + 1.762917, + 8.676415, + 19.98228, + 11.87412, + 20.94647, + 9.686372, + 21.24376, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 21.01835, + 3.011658, + 3.843498, + 3.249645, + 9.771441, + 8.041018, + 2.327877, + 19.21397, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 7.585063, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 4.777777, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 7.656126, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 0.8742778, + 6.330581, + 4.266314, + 6.523858, + 7.44585, + 10.23886, + 3.315528, + 5.560736, + 19.32783, + 2.242471, + 1.46184, + 2.058325, + 8.574683, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(124,Graph0_fx5,Graph0_fy5); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph5 = new TH1F("Graph_Graph5","Graph",124,-3.580114,28.89493); + Graph_Graph5->SetMinimum(0.03102773); + Graph_Graph5->SetMaximum(27.16594); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + Double_t Graph1_fx6[101] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy6[101] = { + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 4.734602, + 1.949477, + 0.9119229, + 5.710012, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 5.493919, + 1.762917, + 8.676415, + 11.87412, + 9.686372, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 3.249645, + 8.041018, + 2.327877, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 4.777777, + 1.353373, + 2.366684, + 7.656126, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 0.8742778, + 6.330581, + 6.523858, + 7.44585, + 10.23886, + 3.315528, + 5.560736, + 2.242471, + 1.46184, + 2.058325, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(101,Graph1_fx6,Graph1_fy6); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph6 = new TH1F("Graph_Graph6","Graph",101,0.4340556,28.53001); + Graph_Graph6->SetMinimum(0.03102773); + Graph_Graph6->SetMaximum(16.33892); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + TF1 *fa1 = new TF1("fa","x*3/4",0,30); + fa1->SetFillColor(19); + fa1->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa1->SetLineColor(ci); + fa1->SetLineWidth(2); + fa1->GetXaxis()->SetLabelFont(42); + fa1->GetXaxis()->SetLabelSize(0.035); + fa1->GetXaxis()->SetTitleSize(0.035); + fa1->GetXaxis()->SetTitleFont(42); + fa1->GetYaxis()->SetLabelFont(42); + fa1->GetYaxis()->SetLabelSize(0.035); + fa1->GetYaxis()->SetTitleSize(0.035); + fa1->GetYaxis()->SetTitleOffset(0); + fa1->GetYaxis()->SetTitleFont(42); + fa1->Draw("same"); + + auto legend2 = new TLegend(0.1,0.9,0.9,1.); + legend2->SetHeader("ThetaCM<40","C"); // option "C" allows to center the heade + legend2->SetTextColor(kRed); + legend2->SetTextSize(0.07); + legend2->Draw("same"); + + c_he7_3->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_4 + TPad *c_he7_4 = new TPad("c_he7_4", "c_he7_4",0.51,0.3433333,0.99,0.6566667); + c_he7_4->Draw(); + c_he7_4->cd(); + c_he7_4->Range(-3.75,-1.575,33.75,14.175); + c_he7_4->SetFillColor(0); + c_he7_4->SetBorderMode(0); + c_he7_4->SetBorderSize(2); + c_he7_4->SetFrameBorderMode(0); + c_he7_4->SetFrameBorderMode(0); + + TH1F *h2__6 = new TH1F("h2__6","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40}",30,0,30); + h2__6->SetBinContent(0,108); + h2__6->SetBinContent(1,2); + h2__6->SetBinContent(3,4); + h2__6->SetBinContent(4,3); + h2__6->SetBinContent(5,5); + h2__6->SetBinContent(6,2); + h2__6->SetBinContent(7,3); + h2__6->SetBinContent(8,4); + h2__6->SetBinContent(9,3); + h2__6->SetBinContent(10,6); + h2__6->SetBinContent(11,3); + h2__6->SetBinContent(12,4); + h2__6->SetBinContent(13,5); + h2__6->SetBinContent(14,6); + h2__6->SetBinContent(15,9); + h2__6->SetBinContent(16,12); + h2__6->SetBinContent(17,6); + h2__6->SetBinContent(18,7); + h2__6->SetBinContent(19,10); + h2__6->SetBinContent(20,4); + h2__6->SetBinContent(21,5); + h2__6->SetBinContent(22,2); + h2__6->SetBinContent(23,3); + h2__6->SetBinContent(24,4); + h2__6->SetBinContent(25,2); + h2__6->SetBinContent(26,9); + h2__6->SetBinContent(27,1); + h2__6->SetEntries(232); + h2__6->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h2__6->SetMarkerColor(ci); + h2__6->SetMarkerStyle(20); + h2__6->GetXaxis()->SetLabelFont(42); + h2__6->GetXaxis()->SetLabelSize(0.035); + h2__6->GetXaxis()->SetTitleSize(0.035); + h2__6->GetXaxis()->SetTitleFont(42); + h2__6->GetYaxis()->SetLabelFont(42); + h2__6->GetYaxis()->SetLabelSize(0.035); + h2__6->GetYaxis()->SetTitleSize(0.035); + h2__6->GetYaxis()->SetTitleOffset(0); + h2__6->GetYaxis()->SetTitleFont(42); + h2__6->GetZaxis()->SetLabelFont(42); + h2__6->GetZaxis()->SetLabelSize(0.035); + h2__6->GetZaxis()->SetTitleSize(0.035); + h2__6->GetZaxis()->SetTitleFont(42); + h2__6->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40}"); + pt->Draw(); + + TH1F *htemp__7 = new TH1F("htemp__7","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__7->SetBinContent(3,1); + htemp__7->SetBinContent(4,1); + htemp__7->SetBinContent(5,2); + htemp__7->SetBinContent(6,1); + htemp__7->SetBinContent(7,1); + htemp__7->SetBinContent(8,1); + htemp__7->SetBinContent(9,3); + htemp__7->SetBinContent(10,4); + htemp__7->SetBinContent(11,3); + htemp__7->SetBinContent(12,4); + htemp__7->SetBinContent(13,4); + htemp__7->SetBinContent(14,6); + htemp__7->SetBinContent(15,8); + htemp__7->SetBinContent(16,12); + htemp__7->SetBinContent(17,5); + htemp__7->SetBinContent(18,7); + htemp__7->SetBinContent(19,10); + htemp__7->SetBinContent(20,3); + htemp__7->SetBinContent(21,5); + htemp__7->SetBinContent(22,2); + htemp__7->SetBinContent(23,3); + htemp__7->SetBinContent(24,4); + htemp__7->SetBinContent(25,2); + htemp__7->SetBinContent(26,8); + htemp__7->SetBinContent(27,1); + htemp__7->SetEntries(101); + htemp__7->SetDirectory(0); + htemp__7->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__7->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__7->SetMarkerColor(ci); + htemp__7->SetMarkerStyle(20); + htemp__7->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__7->GetXaxis()->SetLabelFont(42); + htemp__7->GetXaxis()->SetLabelSize(0.035); + htemp__7->GetXaxis()->SetTitleSize(0.035); + htemp__7->GetXaxis()->SetTitleFont(42); + htemp__7->GetYaxis()->SetLabelFont(42); + htemp__7->GetYaxis()->SetLabelSize(0.035); + htemp__7->GetYaxis()->SetTitleSize(0.035); + htemp__7->GetYaxis()->SetTitleOffset(0); + htemp__7->GetYaxis()->SetTitleFont(42); + htemp__7->GetZaxis()->SetLabelFont(42); + htemp__7->GetZaxis()->SetLabelSize(0.035); + htemp__7->GetZaxis()->SetTitleSize(0.035); + htemp__7->GetZaxis()->SetTitleFont(42); + htemp__7->Draw("same"); + c_he7_4->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_5 + TPad *c_he7_5 = new TPad("c_he7_5", "c_he7_5",0.01,0.01,0.49,0.3233333); + c_he7_5->Draw(); + c_he7_5->cd(); + c_he7_5->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_5->SetFillColor(0); + c_he7_5->SetBorderMode(0); + c_he7_5->SetBorderSize(2); + c_he7_5->SetFrameBorderMode(0); + c_he7_5->SetFrameBorderMode(0); + + TH2F *htemp__8 = new TH2F("htemp__8","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__8->SetDirectory(0); + htemp__8->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__8->SetLineColor(ci); + htemp__8->SetMarkerStyle(20); + htemp__8->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__8->GetXaxis()->SetRange(1,40); + htemp__8->GetXaxis()->SetLabelFont(42); + htemp__8->GetXaxis()->SetLabelSize(0.035); + htemp__8->GetXaxis()->SetTitleSize(0.035); + htemp__8->GetXaxis()->SetTitleFont(42); + htemp__8->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__8->GetYaxis()->SetRange(1,40); + htemp__8->GetYaxis()->SetLabelFont(42); + htemp__8->GetYaxis()->SetLabelSize(0.035); + htemp__8->GetYaxis()->SetTitleSize(0.035); + htemp__8->GetYaxis()->SetTitleOffset(0); + htemp__8->GetYaxis()->SetTitleFont(42); + htemp__8->GetZaxis()->SetLabelFont(42); + htemp__8->GetZaxis()->SetLabelSize(0.035); + htemp__8->GetZaxis()->SetTitleSize(0.035); + htemp__8->GetZaxis()->SetTitleFont(42); + htemp__8->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx7[72] = { + 2.293378, + 2.775385, + 25.618, + 25.618, + 15.84522, + 17.87187, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.44866, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 25.77312, + 14.66246, + 21.70794, + 15.92302, + 23.08188, + 11.01585, + -0.8738604, + 18.01523, + 19.00919, + 14.77225, + 19.56621, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 15.78947, + 14.9445, + 5.574539, + 16.38851, + 20.43209, + 13.05403, + 11.00986, + 25.02393, + 18.07501, + 9.807625, + 8.678942, + 25.70097, + 26.18868, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 16.66522, + 18.17658, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 2.64591, + 15.26418, + 24.80567, + 0.2177003, + 3.090493, + 23.81961, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy7[72] = { + 9.163973, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 0.03447526, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 1.949477, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.3992, + 4.866271, + 5.708108, + 1.499539, + 14.85669, + 2.782254, + 17.36037, + 1.762917, + 11.87412, + 9.686372, + 21.24376, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 2.327877, + 19.21397, + 1.852323, + 2.539789, + 2.678995, + 0.9787855, + 0.2521605, + 3.357897, + 4.577259, + 6.320226, + 2.202726, + 4.411285, + 3.224871, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 8.187073, + 2.693062, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 4.266314, + 3.315528, + 1.46184, + 8.574683, + 1.994378, + 5.511144, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(72,Graph0_fx7,Graph0_fy7); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph7 = new TH1F("Graph_Graph7","Graph",100,-3.580114,28.89493); + Graph_Graph7->SetMinimum(0.03102773); + Graph_Graph7->SetMaximum(27.16594); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + Double_t Graph1_fx8[57] = { + 2.775385, + 25.618, + 25.618, + 15.84522, + 17.87187, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.44866, + 18.51257, + 12.52081, + 25.77312, + 14.66246, + 21.70794, + 15.92302, + 23.08188, + 11.01585, + 18.01523, + 19.00919, + 14.77225, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 15.78947, + 5.574539, + 16.38851, + 20.43209, + 13.05403, + 11.00986, + 25.02393, + 18.07501, + 9.807625, + 8.678942, + 25.70097, + 26.18868, + 14.87749, + 14.87749, + 16.66522, + 18.17658, + 4.07879, + 17.75641, + 15.26418, + 24.80567, + 3.090493, + 23.81961, + 18.85983, + 17.55044, + 15.24921, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy8[57] = { + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 0.03447526, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 1.949477, + 5.710012, + 7.403358, + 3.3992, + 4.866271, + 5.708108, + 1.499539, + 14.85669, + 2.782254, + 1.762917, + 11.87412, + 9.686372, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 2.327877, + 1.852323, + 2.539789, + 2.678995, + 0.9787855, + 0.2521605, + 3.357897, + 4.577259, + 6.320226, + 2.202726, + 4.411285, + 3.224871, + 1.353373, + 2.366684, + 8.187073, + 2.693062, + 1.414101, + 5.831221, + 3.315528, + 1.46184, + 1.994378, + 5.511144, + 1.438863, + 6.547084, + 4.416309, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(57,Graph1_fx8,Graph1_fy8); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph8 = new TH1F("Graph_Graph8","Graph",100,0.4340556,28.53001); + Graph_Graph8->SetMinimum(0.03102773); + Graph_Graph8->SetMaximum(16.33892); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + TF1 *fa2 = new TF1("fa","x*3/4",0,30); + fa2->SetFillColor(19); + fa2->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa2->SetLineColor(ci); + fa2->SetLineWidth(2); + fa2->GetXaxis()->SetLabelFont(42); + fa2->GetXaxis()->SetLabelSize(0.035); + fa2->GetXaxis()->SetTitleSize(0.035); + fa2->GetXaxis()->SetTitleFont(42); + fa2->GetYaxis()->SetLabelFont(42); + fa2->GetYaxis()->SetLabelSize(0.035); + fa2->GetYaxis()->SetTitleSize(0.035); + fa2->GetYaxis()->SetTitleOffset(0); + fa2->GetYaxis()->SetTitleFont(42); + fa2->Draw("same"); + + auto legend3 = new TLegend(0.1,0.9,0.9,1.); + legend3->SetHeader("ThetaCM<30","C"); // option "C" allows to center the heade + legend3->SetTextColor(kRed); + legend3->SetTextSize(0.07); + legend3->Draw("same"); + + c_he7_5->Modified(); + c_he7->cd(); + + + +// ------------>Primitives in pad: c_he7_6 + TPad *c_he7_6 = new TPad("c_he7_6", "c_he7_6",0.51,0.01,0.99,0.3233333); + c_he7_6->Draw(); + c_he7_6->cd(); + c_he7_6->Range(-3.75,-1.05,33.75,9.45); + c_he7_6->SetFillColor(0); + c_he7_6->SetBorderMode(0); + c_he7_6->SetBorderSize(2); + c_he7_6->SetFrameBorderMode(0); + c_he7_6->SetFrameBorderMode(0); + + TH1F *h3__9 = new TH1F("h3__9","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30}",30,0,30); + h3__9->SetBinContent(0,99); + h3__9->SetBinContent(1,2); + h3__9->SetBinContent(3,4); + h3__9->SetBinContent(4,2); + h3__9->SetBinContent(5,4); + h3__9->SetBinContent(6,2); + h3__9->SetBinContent(7,2); + h3__9->SetBinContent(8,2); + h3__9->SetBinContent(9,3); + h3__9->SetBinContent(10,3); + h3__9->SetBinContent(12,2); + h3__9->SetBinContent(13,1); + h3__9->SetBinContent(14,3); + h3__9->SetBinContent(15,6); + h3__9->SetBinContent(16,8); + h3__9->SetBinContent(17,3); + h3__9->SetBinContent(18,3); + h3__9->SetBinContent(19,6); + h3__9->SetBinContent(20,2); + h3__9->SetBinContent(21,2); + h3__9->SetBinContent(22,1); + h3__9->SetBinContent(24,2); + h3__9->SetBinContent(25,1); + h3__9->SetBinContent(26,7); + h3__9->SetBinContent(27,1); + h3__9->SetEntries(171); + h3__9->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h3__9->SetMarkerColor(ci); + h3__9->SetMarkerStyle(20); + h3__9->GetXaxis()->SetLabelFont(42); + h3__9->GetXaxis()->SetLabelSize(0.035); + h3__9->GetXaxis()->SetTitleSize(0.035); + h3__9->GetXaxis()->SetTitleFont(42); + h3__9->GetYaxis()->SetLabelFont(42); + h3__9->GetYaxis()->SetLabelSize(0.035); + h3__9->GetYaxis()->SetTitleSize(0.035); + h3__9->GetYaxis()->SetTitleOffset(0); + h3__9->GetYaxis()->SetTitleFont(42); + h3__9->GetZaxis()->SetLabelFont(42); + h3__9->GetZaxis()->SetLabelSize(0.035); + h3__9->GetZaxis()->SetTitleSize(0.035); + h3__9->GetZaxis()->SetTitleFont(42); + h3__9->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30}"); + pt->Draw(); + + TH1F *htemp__10 = new TH1F("htemp__10","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__10->SetBinContent(3,1); + htemp__10->SetBinContent(4,1); + htemp__10->SetBinContent(5,1); + htemp__10->SetBinContent(6,1); + htemp__10->SetBinContent(7,1); + htemp__10->SetBinContent(8,1); + htemp__10->SetBinContent(9,3); + htemp__10->SetBinContent(10,2); + htemp__10->SetBinContent(12,2); + htemp__10->SetBinContent(13,1); + htemp__10->SetBinContent(14,3); + htemp__10->SetBinContent(15,5); + htemp__10->SetBinContent(16,8); + htemp__10->SetBinContent(17,3); + htemp__10->SetBinContent(18,3); + htemp__10->SetBinContent(19,6); + htemp__10->SetBinContent(20,1); + htemp__10->SetBinContent(21,2); + htemp__10->SetBinContent(22,1); + htemp__10->SetBinContent(24,2); + htemp__10->SetBinContent(25,1); + htemp__10->SetBinContent(26,7); + htemp__10->SetBinContent(27,1); + htemp__10->SetEntries(57); + htemp__10->SetDirectory(0); + htemp__10->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__10->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__10->SetMarkerColor(ci); + htemp__10->SetMarkerStyle(20); + htemp__10->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__10->GetXaxis()->SetLabelFont(42); + htemp__10->GetXaxis()->SetLabelSize(0.035); + htemp__10->GetXaxis()->SetTitleSize(0.035); + htemp__10->GetXaxis()->SetTitleFont(42); + htemp__10->GetYaxis()->SetLabelFont(42); + htemp__10->GetYaxis()->SetLabelSize(0.035); + htemp__10->GetYaxis()->SetTitleSize(0.035); + htemp__10->GetYaxis()->SetTitleOffset(0); + htemp__10->GetYaxis()->SetTitleFont(42); + htemp__10->GetZaxis()->SetLabelFont(42); + htemp__10->GetZaxis()->SetLabelSize(0.035); + htemp__10->GetZaxis()->SetTitleSize(0.035); + htemp__10->GetZaxis()->SetTitleFont(42); + htemp__10->Draw("same"); + c_he7_6->Modified(); + c_he7->cd(); + c_he7->Modified(); + c_he7->cd(); + c_he7->SetSelected(c_he7); +} diff --git a/root/mm.C b/root/mm.C index b772c84257c9d619bf1f35d439f75937ba4ed93b..ad7e7a92f12872b389e5cb3579fb180ed7994f48 100644 --- a/root/mm.C +++ b/root/mm.C @@ -1,7 +1,7 @@ void mm() { //=========Macro generated from canvas: c_angn4/title -//========= (Sun Oct 29 18:31:05 2023) by ROOT version6.10/08 +//========= (Mon Oct 30 11:37:23 2023) by ROOT version6.10/08 TCanvas *c_angn4 = new TCanvas("c_angn4", "title",117,57,1000,1000); gStyle->SetOptStat(0); c_angn4->Range(-4.558442,-1.575,30.50649,14.175); @@ -274,7 +274,7 @@ void mm() Graph_Graph1->GetZaxis()->SetTitleFont(42); graph->SetHistogram(Graph_Graph1); - graph->Draw(""); + graph->Draw("same"); c_angn4->Modified(); c_angn4->cd(); c_angn4->SetSelected(c_angn4); diff --git a/root/single_angular.C b/root/single_angular.C new file mode 100644 index 0000000000000000000000000000000000000000..3cd1c2b42b73e87fafffa0e792f1646ed3149a58 --- /dev/null +++ b/root/single_angular.C @@ -0,0 +1,616 @@ +void single_angular() +{ +//=========Macro generated from canvas: c_angn4/title +//========= (Mon Oct 30 12:28:11 2023) by ROOT version6.10/08 + TCanvas *c_angn4 = new TCanvas("c_angn4", "title",67,55,1800,1000); + gStyle->SetOptStat(0); + c_angn4->Range(-6.875,-5.000001,31.875,95); + c_angn4->SetFillColor(0); + c_angn4->SetBorderMode(0); + c_angn4->SetBorderSize(2); + c_angn4->SetFrameBorderMode(0); + c_angn4->SetFrameBorderMode(0); + + TH2F *htemp__1 = new TH2F("htemp__1","lv_4n_mm_CMR.Theta()*TMath::RadToDeg():lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1}",40,-3,28,40,5,85); + htemp__1->SetDirectory(0); + htemp__1->SetStats(0); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#000099"); + htemp__1->SetLineColor(ci); + htemp__1->SetMarkerStyle(8); + htemp__1->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__1->GetXaxis()->SetRange(1,40); + htemp__1->GetXaxis()->SetLabelFont(42); + htemp__1->GetXaxis()->SetLabelSize(0.035); + htemp__1->GetXaxis()->SetTitleSize(0.035); + htemp__1->GetXaxis()->SetTitleFont(42); + htemp__1->GetYaxis()->SetTitle("lv_4n_mm_CMR.Theta()*TMath::RadToDeg()"); + htemp__1->GetYaxis()->SetRange(1,40); + htemp__1->GetYaxis()->SetLabelFont(42); + htemp__1->GetYaxis()->SetLabelSize(0.035); + htemp__1->GetYaxis()->SetTitleSize(0.035); + htemp__1->GetYaxis()->SetTitleOffset(0); + htemp__1->GetYaxis()->SetTitleFont(42); + htemp__1->GetZaxis()->SetLabelFont(42); + htemp__1->GetZaxis()->SetLabelSize(0.035); + htemp__1->GetZaxis()->SetTitleSize(0.035); + htemp__1->GetZaxis()->SetTitleFont(42); + htemp__1->Draw(""); + + TPaveText *pt = new TPaveText(0.15,0.9344439,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + TText *pt_LaTex = pt->AddText("lv_4n_mm_CMR.Theta()*TMath::RadToDeg():lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1}"); + pt->Draw(); + + Double_t Graph0_fx1[137] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 4.118373, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 9.533222, + 19.7312, + 2.250434, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 11.88648, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 1.294231, + 16.38851, + 12.18283, + 3.692083, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 14.27111, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy1[137] = { + 22.96098, + 32.37453, + 20.42714, + 23.75573, + 23.75573, + 28.97319, + 33.52173, + 28.46252, + 38.64605, + 29.95049, + 24.18883, + 17.91793, + 22.82172, + 24.2163, + 25.36298, + 29.18495, + 17.86598, + 34.39316, + 17.95625, + 31.63668, + 26.6676, + 14.39871, + 21.97711, + 28.91384, + 35.50192, + 21.15304, + 39.67074, + 38.94612, + 15.47257, + 26.06359, + 28.86099, + 32.18519, + 31.19928, + 25.90003, + 25.30423, + 33.05285, + 18.97107, + 78.38055, + 43.17075, + 46.28093, + 33.2535, + 25.42674, + 39.67903, + 28.68732, + 30.61999, + 30.61999, + 25.82045, + 33.3565, + 19.97845, + 25.77351, + 43.49009, + 36.59207, + 27.43145, + 27.43145, + 20.90391, + 32.97262, + 18.1204, + 21.22906, + 47.32078, + 37.22993, + 37.22993, + 37.85647, + 21.61354, + 19.62796, + 34.49465, + 26.35171, + 48.47679, + 26.015, + 38.44679, + 43.79028, + 27.86231, + 26.99938, + 33.08223, + 37.60542, + 28.98927, + 29.76159, + 31.48392, + 32.51186, + 28.55604, + 31.75072, + 32.60382, + 35.39864, + 41.52173, + 21.38447, + 30.2175, + 39.88721, + 14.35446, + 28.19643, + 30.12882, + 29.8008, + 33.42579, + 46.05483, + 33.75885, + 21.86448, + 14.45356, + 14.45356, + 18.33214, + 23.83462, + 39.67193, + 42.60805, + 19.64076, + 31.97388, + 26.98648, + 36.02568, + 26.63101, + 23.97905, + 10.75056, + 20.97361, + 37.34012, + 40.78341, + 30.53927, + 17.6882, + 39.23266, + 32.33729, + 34.5398, + 44.46054, + 27.14791, + 34.37672, + 34.71008, + 38.29793, + 28.51064, + 30.15736, + 18.32953, + 20.052, + 37.68594, + 30.23097, + 29.02963, + 38.83401, + 19.82276, + 29.83649, + 29.06612, + 27.19644, + 38.30944, + 39.88242, + 26.82334, + 24.67148, + 28.65588}; + TGraph *graph = new TGraph(137,Graph0_fx1,Graph0_fy1); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(8); + + TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",137,-3.580114,28.89493); + Graph_Graph1->SetMinimum(3.987557); + Graph_Graph1->SetMaximum(85.14355); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + + Double_t Graph1_fx2[108] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 19.7312, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 14.27111, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy2[108] = { + 32.37453, + 20.42714, + 23.75573, + 23.75573, + 28.97319, + 33.52173, + 28.46252, + 38.64605, + 29.95049, + 24.18883, + 17.91793, + 22.82172, + 24.2163, + 25.36298, + 29.18495, + 34.39316, + 17.95625, + 31.63668, + 14.39871, + 28.91384, + 35.50192, + 21.15304, + 39.67074, + 38.94612, + 26.06359, + 28.86099, + 31.19928, + 25.90003, + 25.30423, + 33.05285, + 18.97107, + 43.17075, + 39.67903, + 28.68732, + 30.61999, + 25.82045, + 19.97845, + 36.59207, + 27.43145, + 27.43145, + 20.90391, + 18.1204, + 21.22906, + 47.32078, + 37.22993, + 37.85647, + 21.61354, + 34.49465, + 26.35171, + 26.015, + 38.44679, + 27.86231, + 26.99938, + 37.60542, + 28.98927, + 29.76159, + 31.48392, + 32.51186, + 28.55604, + 31.75072, + 32.60382, + 35.39864, + 41.52173, + 21.38447, + 30.2175, + 39.88721, + 14.35446, + 28.19643, + 30.12882, + 29.8008, + 33.42579, + 46.05483, + 33.75885, + 14.45356, + 14.45356, + 39.67193, + 42.60805, + 19.64076, + 31.97388, + 26.98648, + 36.02568, + 26.63101, + 23.97905, + 37.34012, + 40.78341, + 30.53927, + 39.23266, + 32.33729, + 34.5398, + 44.46054, + 27.14791, + 34.37672, + 38.29793, + 28.51064, + 30.15736, + 20.052, + 37.68594, + 30.23097, + 29.02963, + 38.83401, + 19.82276, + 29.83649, + 29.06612, + 38.30944, + 39.88242, + 26.82334, + 24.67148, + 28.65588}; + graph = new TGraph(108,Graph1_fx2,Graph1_fy2); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(8); + + TH1F *Graph_Graph2 = new TH1F("Graph_Graph2","Graph",108,0.4340556,28.53001); + Graph_Graph2->SetMinimum(11.05782); + Graph_Graph2->SetMaximum(50.61741); + 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()->SetTitleOffset(0); + 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); + + graph->Draw("p"); + c_angn4->Modified(); + c_angn4->cd(); + c_angn4->SetSelected(c_angn4); +} diff --git a/tetra.tex b/tetra.tex index ba87e816021ed95e65fc13479ec316e173329067..b46c62b9f24119f5c9ead5a95615bd05399b8a63 100644 --- a/tetra.tex +++ b/tetra.tex @@ -1,7 +1,7 @@ \documentclass[12pt,a4paper]{article} \usepackage[utf8]{inputenc} % размер - +\usepackage[final]{pdfpages} \usepackage[T2A]{fontenc} % указывает внутреннюю кодировку TeX \usepackage[ngerman,english]{babel} @@ -105,12 +105,35 @@ The missing mass (MM) spectrum of 4n reconstructed after kinematic selection is %------------------------------------------------------------------------------- Obviously, this spectrum shape can not be described just as a contribution of the 4-body phase space volume, see the orange dotted curve in Fig.\ \ref{fig:mm_4n} (b), indicating the presence of some resonance states of 4n lying at 3-4\,MeV above the 4n decay threshold. +Its important to notice that the observed peak is consistent with the reported in Ref.\ \cite{Duer:2022} and previous Ref.\ \cite{Kisamori:2016} results. +Its important evidence for the tetraneutron resonance at low energies above the 4n decay threshold obtained at forward angles in the reaction of alpha core knockout from the $^8$He beam. + \section{Conclusion} +The $^2$H($^8$He,$^6$Li)4n reaction was studied at forward angles. +The observed peak at 3-4 MeV above the 4n threshold is in agreement with Ref.\ \cite{Duer:2022,Kisamori:2016}. +The obtained results are pushing the investigation of the neutron clusters. +The configuration of the detector setup specifically for such reaction and improved neutron detection efficiency is demanded in the future experiments. +The latter will not also increase statistics but also will allow to study neutron cluster decays in full kinematics reconstructed by the invariant mass methods. + + +\section{Acknowledgement} + +%Institute of Physics of Silesian University in Opava +We also acknowledge the principal support of this work by the Russian Science +Foundation grant No.\ 22-12-00054. The authors are grateful to +Profs.\ Yu.Ts.\ Oganessian and S.N.\ Dmitriev for the long-term support and +development of this activity. We acknowledge important contribution of Prof.\ +M.S.\ Golovkov to the development of the experimental methods and useful +discussions. Also, authors express their gratitude to the +acceleration team for the stable work of U-400M cyclotron during all runs. + \bibliographystyle{unsrt} \clearpage \bibliography{all1} +\includepdf[pages=-]{add.pdf} + \end{document} \ No newline at end of file