‭{‬
	gROOT‭->‬Reset‭();‬
	gROOT‭->‬SetStyle‭("‬Plain‭");‬
‭   ‬gStyle‭->‬SetOptTitle(0‭);‬
‭   ‬gStyle.SetOptStat(0‭);‬
‭   ‬gStyle.SetFillColor(0‭);‬
‭   ‬gStyle.SetPadBorderMode(0‭);‬
‭   ‬m‭ = (‬TH1F‭*)‬gROOT‭->‬FindObject‭("‬h‭");‬
‭   ‬if‭ (‬m‭) ‬m‭->‬Delete‭();‬
‭   ‬TCanvas‭ *‬c1‭ = ‬new TCanvas‭("‬c1‭","‬Spincut parameter‭",‬600,600‭);‬
‭   ‬TH2F‭ *‬h‭ = ‬new TH2F‭("‬h‭"," ",‬10,0.0‭,  ‬15.07,50,0.001‭,    ‬7.4‭);‬
‭   ‬ifstream spincutfile‭("‬spincut.cnt‭");
	‬ifstream spin07("spin_07.cnt")‬;
‭   ‬float energy[70‭],‬spincut[70‭], spin[120], ener[120];‬
	float enererr[120] = {0.};
	float spinerr[120] = {0.};
‭   ‬int i‭ = ‬0‭;‬
‭   ‬float a0‭ =   ‬0.0392‭;‬
‭   ‬float a1‭ =   ‬0.2210‭;‬
‭   ‬float x,y,z‭;‬
‭   ‬while(spincutfile‭){‬
‭   ‬	spincutfile‭ >> ‬x‭;‬
	‭   ‬spincut[i‭]=‬x‭;‬
	‭   ‬energy[i‭]=‬a0+‭(‬a1‭*‬i‭);‬
	‭   ‬i‭++;‬
‭   }‬
	
	i=0;
	while(spin07){
	spin07 >> x;
	ener[i] = x;
	spin[i] =x;
	i++;
	}

‭   ‬TGraph‭ *‬spincutgraph‭ = ‬new TGraph(70‭,‬energy,spincut‭);‬
	TGrapherrors *spin_07 = new TGrapherrors(40, ener, spin, enererr, spinerr);
‭   ‬c1‭.‬SetLeftMargin(0.14‭);‬
‭   ‬h.GetXaxis‭().‬CenterTitle‭();‬
‭   ‬h.GetXaxis‭().‬SetTitle‭("‬Excitation energy E‭ (‬MeV‭)");‬
‭   ‬h.GetYaxis‭().‬CenterTitle‭();‬
‭   ‬h.GetYaxis‭().‬SetTitleOffset(1.4‭);‬
‭   ‬h.GetYaxis‭().‬SetTitle‭("‬Spin cutoff‭ #‬sigma‭");‬
‭   ‬h.Draw‭();‬
‭   ‬spincutgraph.Draw‭("‬L‭");‬
	spin.SetMarkerstyle(20);
	spin.Draw("P");
‭   ‬c1-‭>‬Update‭();‬
‭   ‬c1-‭>‬Print‭("‬spincut.eps‭");‬
‭   ‬c1-‭>‬Print‭("‬spincut.ps‭");‬
‭}‬