data added

root/mm.C

+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);
+}

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}