diff --git a/figures/mm_spec.png b/figures/mm_spec.png new file mode 100644 index 0000000000000000000000000000000000000000..6dbfbba59ab20b0b4abcdbbe4b297b45947bd991 Binary files /dev/null and b/figures/mm_spec.png differ diff --git a/root/mm.C b/root/mm.C new file mode 100644 index 0000000000000000000000000000000000000000..b772c84257c9d619bf1f35d439f75937ba4ed93b --- /dev/null +++ b/root/mm.C @@ -0,0 +1,281 @@ +void mm() +{ +//=========Macro generated from canvas: c_angn4/title +//========= (Sun Oct 29 18:31:05 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); + c_angn4->SetFillColor(0); + c_angn4->SetBorderMode(0); + c_angn4->SetBorderSize(2); + c_angn4->SetLeftMargin(0.13); + c_angn4->SetFrameBorderMode(0); + c_angn4->SetFrameBorderMode(0); + + TH1F *h1__1 = new TH1F("h1__1","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && lv_4n_mm.Mag()-4*939.565>-1 && 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)}",27,0,27); + h1__1->SetBinContent(3,1); + h1__1->SetBinContent(4,2); + h1__1->SetBinContent(5,2); + h1__1->SetBinContent(6,1); + h1__1->SetBinContent(7,1); + h1__1->SetBinContent(8,1); + h1__1->SetBinContent(9,3); + h1__1->SetBinContent(10,4); + h1__1->SetBinContent(11,3); + h1__1->SetBinContent(12,4); + h1__1->SetBinContent(13,4); + h1__1->SetBinContent(14,6); + h1__1->SetBinContent(15,10); + h1__1->SetBinContent(16,12); + h1__1->SetBinContent(17,6); + h1__1->SetBinContent(18,9); + h1__1->SetBinContent(19,10); + h1__1->SetBinContent(20,4); + h1__1->SetBinContent(21,5); + h1__1->SetBinContent(22,2); + h1__1->SetBinContent(23,3); + h1__1->SetBinContent(24,4); + h1__1->SetBinContent(25,2); + h1__1->SetBinContent(26,8); + h1__1->SetBinContent(27,1); + h1__1->SetEntries(108); + h1__1->SetStats(0); + h1__1->GetXaxis()->SetTitle("E_{T}, MeV"); + h1__1->GetXaxis()->CenterTitle(true); + h1__1->GetXaxis()->SetLabelFont(42); + h1__1->GetXaxis()->SetLabelSize(0.035); + h1__1->GetXaxis()->SetTitleSize(0.05); + h1__1->GetXaxis()->SetTitleOffset(0.8); + h1__1->GetXaxis()->SetTitleFont(42); + h1__1->GetYaxis()->SetTitle("Events/1 MeV"); + h1__1->GetYaxis()->CenterTitle(true); + h1__1->GetYaxis()->SetLabelFont(42); + h1__1->GetYaxis()->SetLabelSize(0.035); + h1__1->GetYaxis()->SetTitleSize(0.05); + h1__1->GetYaxis()->SetTitleOffset(0.9); + h1__1->GetYaxis()->SetTitleFont(42); + h1__1->GetZaxis()->SetLabelFont(42); + h1__1->GetZaxis()->SetLabelSize(0.035); + h1__1->GetZaxis()->SetTitleSize(0.035); + h1__1->GetZaxis()->SetTitleFont(42); + h1__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.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && lv_4n_mm.Mag()-4*939.565>-1 && 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)}"); + pt->Draw(); + + Double_t Graph0_fx1[84] = { + -1.75, + -1.25, + -0.75, + -0.25, + 0.25, + 0.75, + 1.25, + 1.75, + 2.25, + 2.75, + 3.25, + 3.75, + 4.25, + 4.75, + 5.25, + 5.75, + 6.25, + 6.75, + 7.25, + 7.75, + 8.25, + 8.75, + 9.25, + 9.75, + 10.25, + 10.75, + 11.25, + 11.75, + 12.25, + 12.75, + 13.25, + 13.75, + 14.25, + 14.75, + 15.25, + 15.75, + 16.25, + 16.75, + 17.25, + 17.75, + 18.25, + 18.75, + 19.25, + 19.75, + 20.25, + 20.75, + 21.25, + 21.75, + 22.25, + 22.75, + 23.25, + 23.75, + 24.25, + 24.75, + 25.25, + 25.75, + 26.25, + 26.75, + 27.25, + 27.75, + 28.25, + 28.75, + 29.25, + 29.75, + 30.25, + 30.75, + 31.25, + 31.75, + 32.25, + 32.75, + 33.25, + 33.75, + 34.25, + 34.75, + 35.25, + 35.75, + 36.25, + 36.75, + 37.25, + 37.75, + 38.25, + 38.75, + 39.25, + 39.75}; + Double_t Graph0_fy1[84] = { + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0.0181, + 0.10863, + 0.14484, + 0.23536, + 0.30778, + 0.59745, + 0.61556, + 0.92334, + 1.14059, + 1.33974, + 1.93719, + 1.9553, + 2.7338, + 3.34936, + 3.56661, + 4.54426, + 5.37707, + 6.08315, + 6.93407, + 8.52728, + 10.08428, + 11.58696, + 14.37507, + 15.2803, + 16.85541, + 19.44437, + 23.78948, + 26.50517, + 30.19851, + 33.2582, + 36.60755, + 42.72692, + 46.42026, + 51.41714, + 58.315, + 63.8188, + 72.09262, + 76.90845, + 84.91069, + 92.27927, + 102.0015, + 110.6374, + 117.1731, + 125.5013, + 133.3949, + 142.1394, + 149.852, + 154.9756, + 161.9821, + 168.4092, + 172.9173, + 176.502, + 175.7959, + 174.782, + 167.504, + 159.6647, + 145.0905, + 132.8698, + 115.598, + 96.44334, + 76.92655, + 52.4491, + 32.57022, + 15.96828, + 6.75303, + 1.46647, + 0.16294, + 0.0181, + 0, + 0, + 0, + 0, + 0, + 0}; + TGraph *graph = new TGraph(84,Graph0_fx1,Graph0_fy1); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + graph->SetFillColor(1); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#ff9933"); + graph->SetLineColor(ci); + graph->SetLineStyle(2); + graph->SetLineWidth(6); + + TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",100,-5.9,43.9); + Graph_Graph1->SetMinimum(0); + Graph_Graph1->SetMaximum(194.1522); + 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(""); + c_angn4->Modified(); + c_angn4->cd(); + c_angn4->SetSelected(c_angn4); +} diff --git a/tetra.tex b/tetra.tex index 1fa0dbbddec63255d31fd2979fc3e1292c5d18f4..ba87e816021ed95e65fc13479ec316e173329067 100644 --- a/tetra.tex +++ b/tetra.tex @@ -91,19 +91,21 @@ The missing mass (MM) spectrum of 4n reconstructed after kinematic selection is \begin{tabular}{cc} \includegraphics[width=0.5\linewidth]{figures/triangle} & - \includegraphics[width=0.5\linewidth]{figures/mm} \\ + \includegraphics[width=0.5\linewidth]{figures/mm_spec} \\ (a) & (b) \end{tabular} \caption{ - (a) The dE-TAC correlation providing clear separation of gammas from neutrons. - (b) The ToF distribution. - The set of three peaks on the left corresponds to the gammas produced in the diaphragm installed 20 cm upstream the target plane, target frame, and CsI(Tl) array. - Green ad red-line histograms formed by the events, identified as gamma/neutron by the dE-TAC method. + (a) Correlation between the neutron energy in the 4n frame and the 4n decay energy. + The kinematical border showing the $E_{n} < 3/4 E_T (4$n) relation is shown with the red line. + (b) The 4n MM spectrum projected from red events in (a) by using the kinematical condition. + The orange dotted curve illustrates the 4-body phase volume $\sim E_T^{7/2}$. } \label{fig:mm_4n} \end{figure} %------------------------------------------------------------------------------- +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. + \section{Conclusion} \bibliographystyle{unsrt}