//---------------------------------------------------------------------------
function drawPot() {
	var pts=200,i,j,k,l,min,max;
	// OM
	/*
	var pot_om=new CountPotential();
	X=new Array();
	Y=new Array();
	pot_om.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_om.setCoulombParameters(R0_COUL.getValue());
	switch (POT_OM.getPotType()) {
		case "volume"		: pot_om.setNuclearParameters("Volume",POT_OM.getValue("_VOL_V"),POT_OM.getValue("_VOL_R"),POT_OM.getValue("_VOL_A"),0.0,0.0,0.0); break;
		case "surface"		: pot_om.setNuclearParameters("Surface",POT_OM.getValue("_SURF_V"),POT_OM.getValue("_SURF_R"),POT_OM.getValue("_SURF_A"),0.0,0.0,0.0); break;
		case "superposition": pot_om.setNuclearParameters("Superposition",POT_OM.getValue("_SUP_VOL_V"),POT_OM.getValue("_SUP_VOL_R"),POT_OM.getValue("_SUP_VOL_A"),
																	   POT_OM.getValue("_SUP_SURF_V"),POT_OM.getValue("_SUP_SURF_R"),POT_OM.getValue("_SUP_SURF_A")); break;
		case "gaussian"		: pot_om.setNuclearParameters("Gaussian",POT_OM.getValue("_GAUS_V"),POT_OM.getValue("_GAUS_R"),0.0,0.0,0.0,0.0); break;
	}
	pot_om.Calculate(0.001,15,pts,X,Y);
	var gf_om=createArray(pts,4);
	min=1.e5; max=-1.e5;
	for (i=0;i<pts;i++) { 
		gf_om[i][0]=X[i]; gf_om[i][1]=Y[i]; gf_om[i][2]=0; gf_om[i][3]=0;	
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}
	// AB
	var pot_ab=new CountPotential();
	pot_ab.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	switch (POT_AB.getPotType()) {
		case "volume"		: pot_ab.setNuclearParameters("Volume",POT_AB.getValue("_VOL_V"),POT_AB.getValue("_VOL_R"),POT_AB.getValue("_VOL_A"),0.0,0.0,0.0); break;
		case "surface"		: pot_ab.setNuclearParameters("Surface",POT_AB.getValue("_SURF_V"),POT_AB.getValue("_SURF_R"),POT_AB.getValue("_SURF_A"),0.0,0.0,0.0); break;
		case "superposition": pot_ab.setNuclearParameters("Superposition",POT_AB.getValue("_SUP_VOL_V"),POT_AB.getValue("_SUP_VOL_R"),POT_AB.getValue("_SUP_VOL_A"),
																	   POT_AB.getValue("_SUP_SURF_V"),POT_AB.getValue("_SUP_SURF_R"),POT_AB.getValue("_SUP_SURF_A")); break;
	}
	pot_ab.Calculate(0.001,15,pts,X,Y);
	*/
	var pot_om=new FRESCO_Potential();
	X=new Array();
	Y=new Array();
	pot_om.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_om.setNTerms(POT_OM.Potentials);
	pot_om.setCoulombParameters(R0_COUL.getValue());
	for (i=0;i<POT_OM.Potentials;i++) {		
		if (POT_OM.Type[i]<2) pot_om.setOneTermInNuclearPotential(i,1*POT_OM.Type[i]+1,POT_OM.Shape[i],POT_OM.V0[i],POT_OM.r0[i],POT_OM.a[i]);
		else pot_om.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_om.Calculate(0.001,15,pts,X,Y);
	var gf_om=createArray(pts,4);
	min=1.e5; max=-1.e5;
	for (i=0;i<pts;i++) { 
		gf_om[i][0]=X[i]; gf_om[i][1]=Y[i]; gf_om[i][2]=0; gf_om[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}	
	//
	var pot_ab=new FRESCO_Potential();
	pot_ab.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_ab.setNTerms(POT_AB.Potentials);
	for (i=0;i<POT_AB.Potentials;i++) {
		if (POT_AB.Type[i]<2) pot_ab.setOneTermInNuclearPotential(i,1*POT_AB.Type[i]+1,POT_AB.Shape[i],POT_AB.V0[i],POT_AB.r0[i],POT_AB.a[i]);
		else pot_ab.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_ab.Calculate(0.001,15,pts,X,Y);
	var gf_ab=createArray(pts,4);
	for (i=0;i<pts;i++) { 
		gf_ab[i][0]=X[i]; gf_ab[i][1]=Y[i]; gf_ab[i][2]=0; gf_ab[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}	
	// Draw OM and AB
	var myChart=new Chart('canvas_pot','12px Calibri','black','14px Calibri','black',"r (fm)","V, W (MeV)",1,true,"gray","lightgoldenrodyellow",6,false,false,false);
	myChart.addData(new ChartData(pts,4,gf_om,getBasicColor(0),0,1,1,1,"",false));
	myChart.addData(new ChartData(pts,4,gf_ab,getBasicColor(1),0,1,1,1,"",false));
	myChart.setXLimits(0,15.);
	myChart.draw(); 
	//-------------------------------------------------------
	// Spin-orbit and tensor
	//-------------------------------------------------------
	var pot_om_so=new FRESCO_Potential();
	pot_om_so.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_om_so.setNTerms(POT_OM.Potentials);
	for (i=0;i<POT_OM.Potentials;i++) {
		if ((POT_OM.Type[i]==2)||(POT_OM.Type[i]==3)) pot_om_so.setOneTermInNuclearPotential(i,1*POT_OM.Type[i]+1,POT_OM.Shape[i],POT_OM.V0[i],POT_OM.r0[i],POT_OM.a[i]);
		else pot_om_so.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_om_so.Calculate(0.001,15,pts,X,Y);
	var gf_om_so=createArray(pts,4);
	min=1.e5; max=-1.e5;
	for (i=0;i<pts;i++) { 
		gf_om_so[i][0]=X[i]; gf_om_so[i][1]=Y[i]; gf_om_so[i][2]=0; gf_om_so[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}	
	//
	var pot_ab_so=new FRESCO_Potential();
	pot_ab_so.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_ab_so.setNTerms(POT_AB.Potentials);
	for (i=0;i<POT_AB.Potentials;i++) {
		if ((POT_AB.Type[i]==2)||(POT_AB.Type[i]==3)) pot_ab_so.setOneTermInNuclearPotential(i,1*POT_AB.Type[i]+1,POT_AB.Shape[i],POT_AB.V0[i],POT_AB.r0[i],POT_AB.a[i]);
		else pot_om_so.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_ab_so.Calculate(0.001,15,pts,X,Y);
	var gf_ab_so=createArray(pts,4);
	for (i=0;i<pts;i++) { 
		gf_ab_so[i][0]=X[i]; gf_ab_so[i][1]=Y[i]; gf_ab_so[i][2]=0; gf_ab_so[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}
	//
	var pot_om_t=new FRESCO_Potential();
	pot_om_t.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_om_t.setNTerms(POT_OM.Potentials);
	for (i=0;i<POT_OM.Potentials;i++) {
		if ((POT_OM.Type[i]==4)||(POT_OM.Type[i]==5)) pot_om_t.setOneTermInNuclearPotential(i,1*POT_OM.Type[i]+1,POT_OM.Shape[i],POT_OM.V0[i],POT_OM.r0[i],POT_OM.a[i]);
		else pot_om_t.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_om_t.Calculate(0.001,15,pts,X,Y);
	var gf_om_t=createArray(pts,4);
	for (i=0;i<pts;i++) { 
		gf_om_t[i][0]=X[i]; gf_om_t[i][1]=Y[i]; gf_om_t[i][2]=0; gf_om_t[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}	
	//
	var pot_ab_t=new FRESCO_Potential();
	pot_ab_t.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_ab_t.setNTerms(POT_AB.Potentials);
	for (i=0;i<POT_AB.Potentials;i++) {
		if ((POT_AB.Type[i]==4)||(POT_AB.Type[i]==5)) pot_ab_t.setOneTermInNuclearPotential(i,1*POT_AB.Type[i]+1,POT_AB.Shape[i],POT_AB.V0[i],POT_AB.r0[i],POT_AB.a[i]);
		else pot_om_t.setOneTermInNuclearPotential(i,1,0,0,1.2,0.5);
	}
	pot_ab_t.Calculate(0.001,15,pts,X,Y);
	var gf_ab_t=createArray(pts,4);
	for (i=0;i<pts;i++) { 
		gf_ab_t[i][0]=X[i]; gf_ab_t[i][1]=Y[i]; gf_ab_t[i][2]=0; gf_ab_t[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}	
	//
	var myChart2=new Chart('canvas_pot_so','12px Calibri','black','14px Calibri','black',"r (fm)","Vso (MeV)",1,true,"gray","lightgoldenrodyellow",6,false,false,false);
	myChart2.addData(new ChartData(pts,4,gf_om_so,getBasicColor(0),0,1,1,1,"",false));
	myChart2.addData(new ChartData(pts,4,gf_ab_so,getBasicColor(1),0,1,1,1,"",false));
	myChart2.addData(new ChartData(pts,4,gf_om_t,getBasicColor(2),0,1,1,1,"",false));
	myChart2.addData(new ChartData(pts,4,gf_ab_t,getBasicColor(3),0,1,1,1,"",false));
	myChart2.setXLimits(0,15.);
	myChart2.draw(); 
	
	
	/*
	var 
		pot_so_re=new CountPotential(),
		pot_so_im=new CountPotential();
	pot_so_re.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_so_im.setReaction(PROJ.getZ(),PROJ.getA(),TARG.getZ(),TARG.getA());
	pot_so_re.setNuclearParameters("Thomas",POT_SO.getValue("_SO_V"),POT_SO.getValue("_SO_R"),POT_SO.getValue("_SO_A"),0.0,0.0,0.0);
	pot_so_im.setNuclearParameters("Thomas",POT_SO.getValue("_SO_W"),POT_SO.getValue("_SO_R"),POT_SO.getValue("_SO_A"),0.0,0.0,0.0);
	
	pot_so_re.Calculate(0.001,15,pts,X,Y);
	var gf_so_re=createArray(pts,4);
	min=Y[0]; max=min;
	for (i=0;i<pts;i++) { 
		gf_so_re[i][0]=X[i]; 
		gf_so_re[i][1]=Y[i]; 
		gf_so_re[i][2]=0; 
		gf_so_re[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}
	pot_so_im.Calculate(0.001,15,pts,X,Y);	
	var gf_so_im=createArray(pts,4);
	for (i=0;i<pts;i++) { 
		gf_so_im[i][0]=X[i]; 
		gf_so_im[i][1]=Y[i]; 
		gf_so_im[i][2]=0; 
		gf_so_im[i][3]=0;
		if (Y[i]>max) max=Y[i];
		if (Y[i]<min) min=Y[i];
	}
	
	// Draw spin-orbit
	var myChart2=new Chart('canvas_pot_so','12px Calibri','black','14px Calibri','black',"r (fm)","VSO (MeV)",1,true,"gray","lightgoldenrodyellow",6,false,false,false);
	myChart2.addData(new Data(pts,4,gf_so_re,getBasicColor(0),0,1,1,1,"",false));
	myChart2.addData(new Data(pts,4,gf_so_im,getBasicColor(1),0,1,1,1,"",false));
	myChart2.setXLimits(0,15.);
	myChart2.draw(); 
	*/
	//
	must_be_recalculate();
}
//---------------------------------------------------------------------------
function must_be_recalculate() {
	var canvas = document.getElementById('canvas_result'); 
    var context = canvas.getContext('2d');
	context.font = '30px Calibri';
	context.fillStyle = "#AAAAAA";
	context.textBaseline = "center";
	context.textAlign = "center";
    context.fillText("Must be recalculated", 250, 200);
}
//---------------------------------------------------------------------------
function check_L_max() {
var A1,A2,Ecm,Elab,EA,m,k,Lmax;
A1=document.getElementById("ID_PROJ_A").value;
A2=document.getElementById("ID_TARG_A").value;
m=1.0*A1*A2/(1.0*A1+1.0*A2);
sys=document.getElementById("ID_ENV_SYS").value;
if (sys=="lab") Ecm=1.0*document.getElementById("ID_ENV_E").value*A2/(1.0*A1+1.0*A2);
else if (sys=="E/A") Ecm=1.0*document.getElementById("ID_ENV_E").value*A1*A2/(1.0*A1+1.0*A2);
else if (sys=="cm") Ecm=1.0*document.getElementById("ID_ENV_E").value;
k=Math.sqrt(2.0*m*Ecm/41.47);
Lmax=document.getElementById("ID_LMAX").value;
if (isNaN(Lmax)) Lmax=1.0*k*document.getElementById("ID_RMAX").value;
if (1*Lmax<0) Lmax=0;
if (1*Lmax>500) Lmax=500;
document.getElementById("ID_LMAX").value=Math.floor(Lmax);
}
//---------------------------------------------------------------------------
function check_R_max() {
var A1,A2,Z1,Z2,m,R_turn,delta,RC_match,amax=0,a,r0=0,r,Rmax;
A1=document.getElementById("ID_PROJ_A").value;
A2=document.getElementById("ID_TARG_A").value;
Z1=document.getElementById("ID_PROJ_Z").value;
Z2=document.getElementById("ID_TARG_Z").value;

m=1.0*A1*A2/(1.0*A1+1.0*A2);
sys=document.getElementById("ID_ENV_SYS").value;
if (sys=="lab") Ecm=1.0*document.getElementById("ID_ENV_E").value*A2/(1.0*A1+1.0*A2);
else if (sys=="E/A") Ecm=1.0*document.getElementById("ID_ENV_E").value*A1*A2/(1.0*A1+1.0*A2);
else if (sys=="cm") Ecm=1.0*document.getElementById("ID_ENV_E").value;
k = Math.sqrt(2.0*m*Ecm/41.47);

// (Coulomb turning point + 2pi/k) - min distance in the case of energy under Coulomb barrier
R_turn = 1.43986*Z1*Z2/Ecm;
delta  = 2*3.1415926/k;
RC_match = R_turn + delta;

//alert(R_turn+"|"+RC_match);

for (i=0;i<POT_OM.Potentials;i++) {
	if (1*POT_OM.r0[i]>r0) r0=POT_OM.r0[i];
	if (1*POT_OM.a[i]>amax) amax=POT_OM.a[i];
}
for (i=0;i<POT_AB.Potentials;i++) {
	if (1*POT_AB.r0[i]>r0) r0=POT_AB.r0[i];
	if (1*POT_AB.a[i]>amax) amax=POT_AB.a[i];
}
Rmax=document.getElementById("ID_RMAX").value;
if (isNaN(Number(Rmax))) Rmax=Number(2*r*(Math.pow(A1,1./3.)+Math.pow(A2,1./3.))+5*amax);
if (1*Rmax < 15) Rmax = 15;
if (1*Rmax > 75) Rmax = 75;
if (1*Rmax < RC_match) Rmax = RC_match;

Rmax = Number(Rmax).toFixed(1);
document.getElementById("ID_RMAX").value=Rmax;
}
//---------------------------------------------------------------------------
function check_int_step() {
var A1,A2,Ecm,Elab,EA,E,m,k,Lmax;
A1=1*PROJ.getA();
A2=1*TARG.getA();
m=1*A1*A2/(1*A1+1*A2);
sys=ENV.getSystem();
E=1*ENV.getE();
if (sys=="lab") Ecm=E*A2/(A1+A2);
else if (sys=="E/A") Ecm=E*A1*A2/(A1+A2);
else if (sys=="cm") Ecm=E;
k=Math.sqrt(2*m*Ecm/41.47);
if (1*INTEGRATION_STEP.getValue()>=0.2/k) {
	m=0.2/k;
	m=m.toFixed(3);
	INTEGRATION_STEP.setValue(m);
}
}
//---------------------------------------------------------------------------
function def_int_params() {
document.getElementById("ID_INITIAL_ANGLE").value=0;
document.getElementById("ID_MAXIMAL_ANGLE").value=180;
document.getElementById("ID_STEP_ANGLE").value=1;
document.getElementById("ID_INTEGRATION_STEP").value=0.05;
RMAX.setValue("none"); // Портим значение, чтобы оно пересчиталось
check_R_max(); // Сперва  R потом L, потому что L зависит от R
LMAX.setValue("none"); // Портим значение, чтобы оно пересчиталось
check_L_max();
check_int_step();
}
//---------------------------------------------------------------------------
function update_expdata_list() {
	$.ajax({
		url: "../save/show_exp_data.php",
		type: "GET",
		data: { task: "inelastic"	},
		success: function( data ) {
			var V,s="&nbsp;"; 
			V=eval('(' +data+ ')');
			for (i=0;i<V.length;i++) s=s+"\n"+V[i].name;
			EXC_EXP1.updateList(s);
			EXC_EXP2.updateList(s);
			EXC_EXP3.updateList(s);
			EXC_EXP4.updateList(s);
			EXC_EXP5.updateList(s);
		}
	});
	$.ajax({
		url: "../save/show_exp_data.php",
		type: "GET",
		data: { task: "elastic"	},
		success: function( data ) {
			var V,s="&nbsp;"; 
			V=eval('(' +data+ ')');
			for (i=0;i<V.length;i++) s=s+"\n"+V[i].name;
			EXC_EXP0.updateList(s);
		}
	});
}
//---------------------------------------------------------------------------
function initValues() {
// Nucleus
PROJ.setZ("2");
PROJ.setA("4");
PROJ.setSpin("0");
PROJ.setParity("+1");
TARG.setZ("12");
TARG.setA("24");
TARG.setSpin("0");
TARG.setParity("+1");
ENV.setE("30.3");
ENV.setSystem("lab");
// Excitations
K.Spin();
K.setValue("0");
BETA2.setValue(0.374);
BETA4.setValue(-0.05);
COULOMB_EXCITATION.setValue("true");
EXCITATION_OF.setValue("target");
EXC_CHECK1.setValue("true");
EXC_CHECK2.setValue("false");
EXC_CHECK3.setValue("false");
EXC_CHECK4.setValue("false");
EXC_CHECK5.setValue("false");
EXC_SPIN1.Spin(); EXC_SPIN1.setValue("2");
EXC_SPIN2.Spin(); 
EXC_SPIN3.Spin(); 
EXC_SPIN4.Spin();
EXC_SPIN5.Spin();
EXC_E1.setValue(1.369);
EXC_E2.setValue(0.0);
EXC_E3.setValue(0.0);
EXC_E4.setValue(0.0);
EXC_E5.setValue(0.0);
EXC_MUL1.setValue(0.1);
EXC_MUL2.setValue(0.01);
EXC_MUL3.setValue(0.001);
EXC_MUL4.setValue(0.0001);
EXC_MUL5.setValue(0.00001);
// Potential
$("#f_pot_tabs").tabs();
R0_COUL.setValue("0.838");
// Integration parameters
def_int_params();
RMAX.setValue("18.0");
LMAX.setValue("30");
drawPot(); // Рисуем потенциалы
VARIANT.list(); // Показываем список вариантов
update_expdata_list(); // Обновляем список экспериментальных данных
excitations(); // Делает доступными только выбранные элементы
setTimeout('EXC_EXP1.setValue("(sample) 4He+24Mg, E=30.3MeV");',1000);
}
//---------------------------------------------------------------------------
function autoname() {
var s;
s=PROJ.getA()+PROJ.getElem()+" + "+TARG.getA()+TARG.getElem()+", E = "+ENV.getE()+" MeV";
VARIANT.set(s);
}
//---------------------------------------------------------------------------
function save() {
var jsonString='\
{\
"Z_PROJ": "'+PROJ.getZ()+'",\
"A_PROJ": "'+PROJ.getA()+'",\
"SPIN_PROJ": "'+PROJ.getSpin()+'",\
"PARITY_PROJ": "'+PROJ.getParity()+'",\
"Z_TARG": "'+TARG.getZ()+'",\
"A_TARG": "'+TARG.getA()+'",\
"SPIN_TARG": "'+TARG.getSpin()+'",\
"PARITY_TARG": "'+TARG.getParity()+'",\
"E": "'+ENV.getE()+'",\
"E_SYS": "'+ENV.getSystem()+'",\
"EXCITATION_OF": "'+EXCITATION_OF.getValue()+'",\
"K": "'+K.getValue()+'",\
"BETA2": "'+BETA2.getValue()+'",\
"BETA4": "'+BETA4.getValue()+'",\
"COULOMB_EXCITATION": "'+COULOMB_EXCITATION.getValue()+'",\
"EXC_CHECK1": "'+EXC_CHECK1.getValue()+'",\
"EXC_CHECK2": "'+EXC_CHECK2.getValue()+'",\
"EXC_CHECK3": "'+EXC_CHECK3.getValue()+'",\
"EXC_CHECK4": "'+EXC_CHECK4.getValue()+'",\
"EXC_CHECK5": "'+EXC_CHECK5.getValue()+'",\
"EXC_SPIN1": "'+EXC_SPIN1.getValue()+'",\
"EXC_SPIN2": "'+EXC_SPIN2.getValue()+'",\
"EXC_SPIN3": "'+EXC_SPIN3.getValue()+'",\
"EXC_SPIN4": "'+EXC_SPIN4.getValue()+'",\
"EXC_SPIN5": "'+EXC_SPIN5.getValue()+'",\
"EXC_EXP0": "'+EXC_EXP0.getValue()+'",\
"EXC_EXP1": "'+EXC_EXP1.getValue()+'",\
"EXC_EXP2": "'+EXC_EXP2.getValue()+'",\
"EXC_EXP3": "'+EXC_EXP3.getValue()+'",\
"EXC_EXP4": "'+EXC_EXP4.getValue()+'",\
"EXC_EXP5": "'+EXC_EXP5.getValue()+'",\
"EXC_E1": "'+EXC_E1.getValue()+'",\
"EXC_E2": "'+EXC_E2.getValue()+'",\
"EXC_E3": "'+EXC_E3.getValue()+'",\
"EXC_E4": "'+EXC_E4.getValue()+'",\
"EXC_E5": "'+EXC_E5.getValue()+'",\
"EXC_MUL1": "'+EXC_MUL1.getValue()+'",\
"EXC_MUL2": "'+EXC_MUL2.getValue()+'",\
"EXC_MUL3": "'+EXC_MUL3.getValue()+'",\
"EXC_MUL4": "'+EXC_MUL4.getValue()+'",\
"EXC_MUL5": "'+EXC_MUL5.getValue()+'",\
"OM": {\
"POTENTIALS": '+POT_OM.Potentials+',\
"POT1": { "TYPE": "'+POT_OM.getValue(0,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(0,"SHAPE")+'", "V0": "'+POT_OM.getValue(0,"V0")+'", "R0": "'+POT_OM.getValue(0,"R0")+'", "A": "'+POT_OM.getValue(0,"A")+'" },\
"POT2": { "TYPE": "'+POT_OM.getValue(1,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(1,"SHAPE")+'", "V0": "'+POT_OM.getValue(1,"V0")+'", "R0": "'+POT_OM.getValue(1,"R0")+'", "A": "'+POT_OM.getValue(1,"A")+'" },\
"POT3": { "TYPE": "'+POT_OM.getValue(2,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(2,"SHAPE")+'", "V0": "'+POT_OM.getValue(2,"V0")+'", "R0": "'+POT_OM.getValue(2,"R0")+'", "A": "'+POT_OM.getValue(2,"A")+'" },\
"POT4": { "TYPE": "'+POT_OM.getValue(3,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(3,"SHAPE")+'", "V0": "'+POT_OM.getValue(3,"V0")+'", "R0": "'+POT_OM.getValue(3,"R0")+'", "A": "'+POT_OM.getValue(3,"A")+'" },\
"POT5": { "TYPE": "'+POT_OM.getValue(4,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(4,"SHAPE")+'", "V0": "'+POT_OM.getValue(4,"V0")+'", "R0": "'+POT_OM.getValue(4,"R0")+'", "A": "'+POT_OM.getValue(4,"A")+'" },\
"POT6": { "TYPE": "'+POT_OM.getValue(5,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(5,"SHAPE")+'", "V0": "'+POT_OM.getValue(5,"V0")+'", "R0": "'+POT_OM.getValue(5,"R0")+'", "A": "'+POT_OM.getValue(5,"A")+'" },\
"POT7": { "TYPE": "'+POT_OM.getValue(6,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(6,"SHAPE")+'", "V0": "'+POT_OM.getValue(6,"V0")+'", "R0": "'+POT_OM.getValue(6,"R0")+'", "A": "'+POT_OM.getValue(6,"A")+'" },\
"POT8": { "TYPE": "'+POT_OM.getValue(7,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(7,"SHAPE")+'", "V0": "'+POT_OM.getValue(7,"V0")+'", "R0": "'+POT_OM.getValue(7,"R0")+'", "A": "'+POT_OM.getValue(7,"A")+'" },\
"POT9": { "TYPE": "'+POT_OM.getValue(8,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(8,"SHAPE")+'", "V0": "'+POT_OM.getValue(8,"V0")+'", "R0": "'+POT_OM.getValue(8,"R0")+'", "A": "'+POT_OM.getValue(8,"A")+'" },\
"POT10": { "TYPE": "'+POT_OM.getValue(9,"TYPE")+'", "SHAPE": "'+POT_OM.getValue(9,"SHAPE")+'", "V0": "'+POT_OM.getValue(9,"V0")+'", "R0": "'+POT_OM.getValue(9,"R0")+'", "A": "'+POT_OM.getValue(9,"A")+'" }\
},\
"AB": {\
"POTENTIALS": '+POT_AB.Potentials+',\
"POT1": { "TYPE": "'+POT_AB.getValue(0,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(0,"SHAPE")+'", "V0": "'+POT_AB.getValue(0,"V0")+'", "R0": "'+POT_AB.getValue(0,"R0")+'", "A": "'+POT_AB.getValue(0,"A")+'" },\
"POT2": { "TYPE": "'+POT_AB.getValue(1,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(1,"SHAPE")+'", "V0": "'+POT_AB.getValue(1,"V0")+'", "R0": "'+POT_AB.getValue(1,"R0")+'", "A": "'+POT_AB.getValue(1,"A")+'" },\
"POT3": { "TYPE": "'+POT_AB.getValue(2,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(2,"SHAPE")+'", "V0": "'+POT_AB.getValue(2,"V0")+'", "R0": "'+POT_AB.getValue(2,"R0")+'", "A": "'+POT_AB.getValue(2,"A")+'" },\
"POT4": { "TYPE": "'+POT_AB.getValue(3,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(3,"SHAPE")+'", "V0": "'+POT_AB.getValue(3,"V0")+'", "R0": "'+POT_AB.getValue(3,"R0")+'", "A": "'+POT_AB.getValue(3,"A")+'" },\
"POT5": { "TYPE": "'+POT_AB.getValue(4,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(4,"SHAPE")+'", "V0": "'+POT_AB.getValue(4,"V0")+'", "R0": "'+POT_AB.getValue(4,"R0")+'", "A": "'+POT_AB.getValue(4,"A")+'" },\
"POT6": { "TYPE": "'+POT_AB.getValue(5,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(5,"SHAPE")+'", "V0": "'+POT_AB.getValue(5,"V0")+'", "R0": "'+POT_AB.getValue(5,"R0")+'", "A": "'+POT_AB.getValue(5,"A")+'" },\
"POT7": { "TYPE": "'+POT_AB.getValue(6,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(6,"SHAPE")+'", "V0": "'+POT_AB.getValue(6,"V0")+'", "R0": "'+POT_AB.getValue(6,"R0")+'", "A": "'+POT_AB.getValue(6,"A")+'" },\
"POT8": { "TYPE": "'+POT_AB.getValue(7,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(7,"SHAPE")+'", "V0": "'+POT_AB.getValue(7,"V0")+'", "R0": "'+POT_AB.getValue(7,"R0")+'", "A": "'+POT_AB.getValue(7,"A")+'" },\
"POT9": { "TYPE": "'+POT_AB.getValue(8,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(8,"SHAPE")+'", "V0": "'+POT_AB.getValue(8,"V0")+'", "R0": "'+POT_AB.getValue(8,"R0")+'", "A": "'+POT_AB.getValue(8,"A")+'" },\
"POT10": { "TYPE": "'+POT_AB.getValue(9,"TYPE")+'", "SHAPE": "'+POT_AB.getValue(9,"SHAPE")+'", "V0": "'+POT_AB.getValue(9,"V0")+'", "R0": "'+POT_AB.getValue(9,"R0")+'", "A": "'+POT_AB.getValue(9,"A")+'" }\
},\
"R0_COUL": "'+R0_COUL.getValue()+'",\
"INITIAL_ANGLE": "'+INITIAL_ANGLE.getValue()+'",\
"MAXIMAL_ANGLE": "'+MAXIMAL_ANGLE.getValue()+'",\
"STEP_ANGLE": "'+STEP_ANGLE.getValue()+'",\
"INTEGRATION_STEP": "'+INTEGRATION_STEP.getValue()+'",\
"LMAX": "'+LMAX.getValue()+'",\
"RMAX": "'+RMAX.getValue()+'"\
}\
';
return(jsonString);
}
//---------------------------------------------------------------------------
function load(data) {
var obj=JSON.parse(data);
PROJ.setZ(obj.Z_PROJ);
PROJ.setA(obj.A_PROJ);
PROJ.setSpin(obj.SPIN_PROJ);
PROJ.setParity(obj.PARITY_PROJ);
TARG.setZ(obj.Z_TARG);
TARG.setA(obj.A_TARG);
TARG.setSpin(obj.SPIN_TARG);
TARG.setParity(obj.PARITY_TARG);
ENV.setE(obj.E);
ENV.setSystem(obj.E_SYS);
EXCITATION_OF.setValue(obj.EXCITATION_OF);
K.setValue(obj.K);
BETA2.setValue(obj.BETA2);
BETA4.setValue(obj.BETA4);
COULOMB_EXCITATION.setValue(obj.COULOMB_EXCITATION);
EXC_CHECK1.setValue(obj.EXC_CHECK1);
EXC_CHECK2.setValue(obj.EXC_CHECK2);
EXC_CHECK3.setValue(obj.EXC_CHECK3);
EXC_CHECK4.setValue(obj.EXC_CHECK4);
EXC_CHECK5.setValue(obj.EXC_CHECK5);
EXC_SPIN1.setValue(obj.EXC_SPIN1);
EXC_SPIN2.setValue(obj.EXC_SPIN2);
EXC_SPIN3.setValue(obj.EXC_SPIN3);
EXC_SPIN4.setValue(obj.EXC_SPIN4);
EXC_SPIN5.setValue(obj.EXC_SPIN5);
EXC_EXP0.setValue(obj.EXC_EXP0);
EXC_EXP1.setValue(obj.EXC_EXP1);
EXC_EXP2.setValue(obj.EXC_EXP2);
EXC_EXP3.setValue(obj.EXC_EXP3);
EXC_EXP4.setValue(obj.EXC_EXP4);
EXC_EXP5.setValue(obj.EXC_EXP5);
EXC_E1.setValue(obj.EXC_E1);
EXC_E2.setValue(obj.EXC_E2);
EXC_E3.setValue(obj.EXC_E3);
EXC_E4.setValue(obj.EXC_E4);
EXC_E5.setValue(obj.EXC_E5);
EXC_MUL1.setValue(obj.EXC_MUL1);
EXC_MUL2.setValue(obj.EXC_MUL2);
EXC_MUL3.setValue(obj.EXC_MUL3);
EXC_MUL4.setValue(obj.EXC_MUL4);
EXC_MUL5.setValue(obj.EXC_MUL5);

POT_OM.Potentials=obj.OM.POTENTIALS;
POT_OM.setValue(0,"TYPE", Number(obj.OM.POT1.TYPE));
POT_OM.setValue(0,"SHAPE",Number(obj.OM.POT1.SHAPE));
POT_OM.setValue(0,"V0",   Number(obj.OM.POT1.V0));
POT_OM.setValue(0,"R0",   Number(obj.OM.POT1.R0));
POT_OM.setValue(0,"A",    Number(obj.OM.POT1.A));
POT_OM.setValue(1,"TYPE", Number(obj.OM.POT2.TYPE));
POT_OM.setValue(1,"SHAPE",Number(obj.OM.POT2.SHAPE));
POT_OM.setValue(1,"V0",   Number(obj.OM.POT2.V0));
POT_OM.setValue(1,"R0",   Number(obj.OM.POT2.R0));
POT_OM.setValue(1,"A",    Number(obj.OM.POT2.A));
POT_OM.setValue(2,"TYPE", Number(obj.OM.POT3.TYPE));
POT_OM.setValue(2,"SHAPE",Number(obj.OM.POT3.SHAPE));
POT_OM.setValue(2,"V0",   Number(obj.OM.POT3.V0));
POT_OM.setValue(2,"R0",   Number(obj.OM.POT3.R0));
POT_OM.setValue(2,"A",    Number(obj.OM.POT3.A));
POT_OM.setValue(3,"TYPE", Number(obj.OM.POT4.TYPE));
POT_OM.setValue(3,"SHAPE",Number(obj.OM.POT4.SHAPE));
POT_OM.setValue(3,"V0",   Number(obj.OM.POT4.V0));
POT_OM.setValue(3,"R0",   Number(obj.OM.POT4.R0));
POT_OM.setValue(3,"A",    Number(obj.OM.POT4.A));
POT_OM.setValue(4,"TYPE", Number(obj.OM.POT5.TYPE));
POT_OM.setValue(4,"SHAPE",Number(obj.OM.POT5.SHAPE));
POT_OM.setValue(4,"V0",   Number(obj.OM.POT5.V0));
POT_OM.setValue(4,"R0",   Number(obj.OM.POT5.R0));
POT_OM.setValue(4,"A",    Number(obj.OM.POT5.A));
POT_OM.setValue(5,"TYPE", obj.OM.POT6.TYPE);
POT_OM.setValue(5,"SHAPE",obj.OM.POT6.SHAPE);
POT_OM.setValue(5,"V0",   obj.OM.POT6.V0);
POT_OM.setValue(5,"R0",   obj.OM.POT6.R0);
POT_OM.setValue(5,"A",    obj.OM.POT6.A);
POT_OM.setValue(6,"TYPE", obj.OM.POT7.TYPE);
POT_OM.setValue(6,"SHAPE",obj.OM.POT7.SHAPE);
POT_OM.setValue(6,"V0",   obj.OM.POT7.V0);
POT_OM.setValue(6,"R0",   obj.OM.POT7.R0);
POT_OM.setValue(6,"A",    obj.OM.POT7.A);
POT_OM.setValue(7,"TYPE", obj.OM.POT8.TYPE);
POT_OM.setValue(7,"SHAPE",obj.OM.POT8.SHAPE);
POT_OM.setValue(7,"V0",   obj.OM.POT8.V0);
POT_OM.setValue(7,"R0",   obj.OM.POT8.R0);
POT_OM.setValue(7,"A",    obj.OM.POT8.A);
POT_OM.setValue(8,"TYPE", obj.OM.POT9.TYPE);
POT_OM.setValue(8,"SHAPE",obj.OM.POT9.SHAPE);
POT_OM.setValue(8,"V0",   obj.OM.POT9.V0);
POT_OM.setValue(8,"R0",   obj.OM.POT9.R0);
POT_OM.setValue(8,"A",    obj.OM.POT9.A);
POT_OM.setValue(9,"TYPE", obj.OM.POT10.TYPE);
POT_OM.setValue(9,"SHAPE",obj.OM.POT10.SHAPE);
POT_OM.setValue(9,"V0",   obj.OM.POT10.V0);
POT_OM.setValue(9,"R0",   obj.OM.POT10.R0);
POT_OM.setValue(9,"A",    obj.OM.POT10.A);
POT_OM.Show();
POT_AB.Potentials=obj.AB.POTENTIALS;
POT_AB.setValue(0,"TYPE", Number(obj.AB.POT1.TYPE));
POT_AB.setValue(0,"SHAPE",Number(obj.AB.POT1.SHAPE));
POT_AB.setValue(0,"V0",   Number(obj.AB.POT1.V0));
POT_AB.setValue(0,"R0",   Number(obj.AB.POT1.R0));
POT_AB.setValue(0,"A",    Number(obj.AB.POT1.A));
POT_AB.setValue(1,"TYPE", Number(obj.AB.POT2.TYPE));
POT_AB.setValue(1,"SHAPE",Number(obj.AB.POT2.SHAPE));
POT_AB.setValue(1,"V0",   Number(obj.AB.POT2.V0));
POT_AB.setValue(1,"R0",   Number(obj.AB.POT2.R0));
POT_AB.setValue(1,"A",    Number(obj.AB.POT2.A));
POT_AB.setValue(2,"TYPE", Number(obj.AB.POT3.TYPE));
POT_AB.setValue(2,"SHAPE",Number(obj.AB.POT3.SHAPE));
POT_AB.setValue(2,"V0",   Number(obj.AB.POT3.V0));
POT_AB.setValue(2,"R0",   Number(obj.AB.POT3.R0));
POT_AB.setValue(2,"A",    Number(obj.AB.POT3.A));
POT_AB.setValue(3,"TYPE", Number(obj.AB.POT4.TYPE));
POT_AB.setValue(3,"SHAPE",Number(obj.AB.POT4.SHAPE));
POT_AB.setValue(3,"V0",   Number(obj.AB.POT4.V0));
POT_AB.setValue(3,"R0",   Number(obj.AB.POT4.R0));
POT_AB.setValue(3,"A",    Number(obj.AB.POT4.A));
POT_AB.setValue(4,"TYPE", Number(obj.AB.POT5.TYPE));
POT_AB.setValue(4,"SHAPE",Number(obj.AB.POT5.SHAPE));
POT_AB.setValue(4,"V0",   Number(obj.AB.POT5.V0));
POT_AB.setValue(4,"R0",   Number(obj.AB.POT5.R0));
POT_AB.setValue(4,"A",    Number(obj.AB.POT5.A));
POT_AB.setValue(5,"TYPE", obj.AB.POT6.TYPE);
POT_AB.setValue(5,"SHAPE",obj.AB.POT6.SHAPE);
POT_AB.setValue(5,"V0",   obj.AB.POT6.V0);
POT_AB.setValue(5,"R0",   obj.AB.POT6.R0);
POT_AB.setValue(5,"A",    obj.AB.POT6.A);
POT_AB.setValue(6,"TYPE", obj.AB.POT7.TYPE);
POT_AB.setValue(6,"SHAPE",obj.AB.POT7.SHAPE);
POT_AB.setValue(6,"V0",   obj.AB.POT7.V0);
POT_AB.setValue(6,"R0",   obj.AB.POT7.R0);
POT_AB.setValue(6,"A",    obj.AB.POT7.A);
POT_AB.setValue(7,"TYPE", obj.AB.POT8.TYPE);
POT_AB.setValue(7,"SHAPE",obj.AB.POT8.SHAPE);
POT_AB.setValue(7,"V0",   obj.AB.POT8.V0);
POT_AB.setValue(7,"R0",   obj.AB.POT8.R0);
POT_AB.setValue(7,"A",    obj.AB.POT8.A);
POT_AB.setValue(8,"TYPE", obj.AB.POT9.TYPE);
POT_AB.setValue(8,"SHAPE",obj.AB.POT9.SHAPE);
POT_AB.setValue(8,"V0",   obj.AB.POT9.V0);
POT_AB.setValue(8,"R0",   obj.AB.POT9.R0);
POT_AB.setValue(8,"A",    obj.AB.POT9.A);
POT_AB.setValue(9,"TYPE", obj.AB.POT10.TYPE);
POT_AB.setValue(9,"SHAPE",obj.AB.POT10.SHAPE);
POT_AB.setValue(9,"V0",   obj.AB.POT10.V0);
POT_AB.setValue(9,"R0",   obj.AB.POT10.R0);
POT_AB.setValue(9,"A",    obj.AB.POT10.A);
POT_AB.Show();
R0_COUL.setValue(obj.R0_COUL);
INITIAL_ANGLE.setValue(obj.INITIAL_ANGLE);
MAXIMAL_ANGLE.setValue(obj.MAXIMAL_ANGLE);
STEP_ANGLE.setValue(obj.STEP_ANGLE);
INTEGRATION_STEP.setValue(obj.INTEGRATION_STEP);
LMAX.setValue(obj.LMAX);
RMAX.setValue(obj.RMAX);
must_be_recalculate();
excitations();
drawPot();
$("#tabs").tabs("option","active",1);
}
//---------------------------------------------------------------------------
function calculate() {
	// Выводим сообщение
	var canvas = document.getElementById('canvas_result'); 
    var context = canvas.getContext('2d');
	context.font = '30px Calibri';
	context.fillStyle = "#AA0000";
	context.textBaseline = "center";
	context.textAlign = "center";
    context.fillText("Calculation",250,250);
	// Готовим данные
	var s=save();
	// Отправляем данные
	$.ajax({
		url: "calculate.php",
		type: "POST",
		data: { data: s	},
		success: result
	});	
}
//---------------------------------------------------------------------------
function excitations() {
if (EXC_CHECK1.getValue()==true) { EXC_SPIN1.enable(); EXC_E1.enable(); EXC_EXP1.enable(); EXC_MUL1.enable(); } else { EXC_SPIN1.disable(); EXC_E1.disable(); EXC_EXP1.disable(); EXC_MUL1.disable(); }
if (EXC_CHECK2.getValue()==true) { EXC_SPIN2.enable(); EXC_E2.enable(); EXC_EXP2.enable(); EXC_MUL2.enable(); } else { EXC_SPIN2.disable(); EXC_E2.disable(); EXC_EXP2.disable(); EXC_MUL2.disable(); }
if (EXC_CHECK3.getValue()==true) { EXC_SPIN3.enable(); EXC_E3.enable(); EXC_EXP3.enable(); EXC_MUL3.enable(); } else { EXC_SPIN3.disable(); EXC_E3.disable(); EXC_EXP3.disable(); EXC_MUL3.disable(); }
if (EXC_CHECK4.getValue()==true) { EXC_SPIN4.enable(); EXC_E4.enable(); EXC_EXP4.enable(); EXC_MUL4.enable(); } else { EXC_SPIN4.disable(); EXC_E4.disable(); EXC_EXP4.disable(); EXC_MUL4.disable(); }
if (EXC_CHECK5.getValue()==true) { EXC_SPIN5.enable(); EXC_E5.enable(); EXC_EXP5.enable(); EXC_MUL5.enable(); } else { EXC_SPIN5.disable(); EXC_E5.disable(); EXC_EXP5.disable(); EXC_MUL5.disable(); }
}
//---------------------------------------------------------------------------
function result(data) {
var i,j,pts,legend,sigmaR,scale;
// Меняем заголовок
if (EXCITATION_OF.getValue()=="projectile") { p1="*"; p2=""; } else { p1=""; p2="*"; }
document.getElementById("title_reaction").innerHTML="<sup>"+TARG.getA()+"</sup>"+Nuclei.Name[TARG.getZ()]+"(<sup>"+PROJ.getA()+"</sup>"+Nuclei.Name[PROJ.getZ()]+",<sup>"+PROJ.getA()+"</sup>"+Nuclei.Name[PROJ.getZ()]+p1+")<sup>"+TARG.getA()+"</sup>"+Nuclei.Name[TARG.getZ()]+p2+", E<sub>"+ENV.getSystem()+"</sub>="+ENV.getE()+" MeV";
// Создаём Chart
var myChart=new Chart('canvas_result','14px Arial','black','16px Calibri','black',"Θcm (deg)","dσ/dΩ (mb/sr)",1,true,"gray","lightgoldenrodyellow",6,true,false,true);
// Находим Ecm
var A1,A2,Ecm,Elab,EA,E,m;
A1=1*PROJ.getA();
A2=1*TARG.getA();
m=1*A1*A2/(1*A1+1*A2);
sys=ENV.getSystem();
E=1*ENV.getE();
if (sys=="lab") Ecm=E*A2/(A1+A2);
else if (sys=="E/A") Ecm=E*A1*A2/(A1+A2);
else if (sys=="cm") Ecm=E;
// Строим насчитанные графики
D=data.split("\n");	
for (i=0;i<6;i++) { // Цикл по файлам результатам
	if (D[i]!="") {
		// Находим кол-во точек
		C=D[i].split("\\n");
		pts=C.length;
		var gf=createArray(pts,2);
		// Определяем значения
		for (j=0;j<pts;j++) {
			s=getValuesFromString(C[j]);
			gf[j][0]=1*s[0]; 
			switch (i) {
				case 0:	
					gf[j][1]=s[1]*10.*Math.pow(PROJ.getZ()*TARG.getZ()*1.43986/4/Ecm,2)/Math.pow(Math.sin(0.5*gf[j][0]/180*3.1415927),4); 
					break;
				case 1:	gf[j][1]=s[1]*EXC_MUL1.getValue();  break;
				case 2:	gf[j][1]=s[1]*EXC_MUL2.getValue();  break;
				case 3:	gf[j][1]=s[1]*EXC_MUL3.getValue();  break;
				case 4:	gf[j][1]=s[1]*EXC_MUL4.getValue();  break;
				case 5:	gf[j][1]=s[1]*EXC_MUL5.getValue();  break;
			}
			//gf[j][2]=0; gf[j][3]=0;
		}
		if (EXCITATION_OF.getValue()=="projectile") parity=PROJ.getParity(); else parity=TARG.getParity();
		parity=parity.substr(0,1);
		switch (i) {
			case 0: legend="g.s."; break;
			case 1: legend=EXC_SPIN1.getValue()+parity+", "+EXC_E1.getValue()+" MeV"; break;
			case 2: legend=EXC_SPIN2.getValue()+parity+", "+EXC_E2.getValue()+" MeV"; break;
			case 3: legend=EXC_SPIN3.getValue()+parity+", "+EXC_E3.getValue()+" MeV"; break;
			case 4: legend=EXC_SPIN4.getValue()+parity+", "+EXC_E4.getValue()+" MeV"; break;
			case 5: legend=EXC_SPIN5.getValue()+parity+", "+EXC_E5.getValue()+" MeV"; break;
		}
		myChart.addData(new ChartData(pts,2,gf,getBasicColor(i),0,1,1,1,legend,true));
	}
}
// Строим графики с экспериментальными точками
// Начнём с G.S. - elastic
if (D[6]!="") {
	var obj=JSON.parse(D[i]);
	if (obj.DATA!="") {
		// Находим кол-во точек
		C=obj.DATA.split("\n");
		pts=C.length;
		var gf=createArray(pts,4);
		var Angle=new Array();
		var CS=new Array();
		var dCS=new Array();
		// Определяем значения
		for (j=0;j<pts;j++) {
			s=getValuesFromString(C[j]);
			Angle[j]=1*s[0]; 
			// Единицы измерения
			kk=1;
			if (obj.CROSS_SECTION=="b") kk=1e3; 
			else if (obj.CROSS_SECTION=="mb") kk=1; 
			else if (obj.CROSS_SECTION=="μb") kk=1e-3; 
			else if (obj.CROSS_SECTION=="nb") kk=1e-6; 
			else if (obj.CROSS_SECTION=="pb") kk=1e-9; 
			CS[j]=kk*s[1]; 
			// Ошибки
			if (obj.ERROR=="abs") dCS[j]=kk*s[2]; 
			else 
			{
				if(s[2]>99) s[2]=99;
				dCS[j]=1*s[2];
			}
		}
		Ecm=energy_recalculate2(obj.E,obj.A_PROJ,obj.A_TARG,obj.E_SYS,"cm");
		if (obj.ANGLE=="lab") 
			CS_lab_to_cm(obj.A_PROJ,obj.A_TARG,obj.A_PROJ,obj.A_TARG,Ecm,0,pts,Angle,CS,dCS,obj.CROSS_SECTION,obj.ERROR);
		if (obj.CROSS_SECTION=="σ/σ(Ruth)") 
			CS_Ratio_to_mb(obj.Z_PROJ,obj.A_PROJ,obj.Z_TARG,obj.A_TARG,Ecm,pts,Angle,CS,dCS,"cm",obj.CROSS_SECTION,obj.ERROR);
		for (i=0;i<pts;i++) 
		{ 
			gf[i][0]=Angle[i]; gf[i][1]=CS[i]; 
			if(obj.ERROR != "abs") dCS[i] = 0.01*dCS[i]*CS[i]; 
			gf[i][2]=dCS[i]; gf[i][3]=gf[i][2]; 
		}
		myChart.addData(new ChartData(pts,4,gf,getBasicColor(0),1,2,0,1,"",true));
	}
}
// Теперь другие данные
for (i=7;i<12;i++) { // Цикл по файлам результатам
	// Scale factor
	switch (i) {
		case  7: scale=1*EXC_MUL1.getValue(); break;
		case  8: scale=1*EXC_MUL2.getValue(); break;
		case  9: scale=1*EXC_MUL3.getValue(); break;
		case 10: scale=1*EXC_MUL4.getValue(); break;
		case 11: scale=1*EXC_MUL5.getValue(); break;
	}
	if (D[i]!="") {
		var obj=JSON.parse(D[i]);
		//alert(D[i]);
		if (obj.DATA!="") {
			// Находим кол-во точек
			C=obj.DATA.split("\n");
			pts=C.length;
			var gf=createArray(pts,4);
			// Определяем значения
			for (j=0;j<pts;j++) {
				s=getValuesFromString(C[j]);
				gf[j][0]=1*s[0]; 
				switch (i) {
					case  6: gf[j][1]=s[1]*scale; break;
					case  7: gf[j][1]=s[1]*scale; break;
					case  8: gf[j][1]=s[1]*scale; break;
					case  9: gf[j][1]=s[1]*scale; break;
					case 10: gf[j][1]=s[1]*scale; break;
				}
				// Единицы измерения
				kk=1;
				if (obj.CROSS_SECTION=="b") kk=1e3; 
				else if (obj.CROSS_SECTION=="mb") kk=1; 
				else if (obj.CROSS_SECTION=="μb") kk=1e-3; 
				else if (obj.CROSS_SECTION=="nb") kk=1e-6; 
				else if (obj.CROSS_SECTION=="pb") kk=1e-9; 
				gf[j][1]=kk*gf[j][1]; 
				// Ошибки
				if (obj.ERROR=="abs") gf[j][2]=scale*kk*s[2]; else gf[j][2]=0.01*gf[j][1]*s[2]; 
				gf[j][3]=gf[j][2];
			}
			myChart.addData(new ChartData(pts,4,gf,getBasicColor(i-6),1,2,0,1,"",true));
		}
	}
}
myChart.setXNLimits(0.0,180.,6);
myChart.setYmax(4);
myChart.draw(); 
}
//---------------------------------------------------------------------------
function getValuesFromString(s) {
w=s.trim();
w=w.replace(new RegExp("    ",'g')," ");
w=w.replace(new RegExp("   ",'g')," ");
w=w.replace(new RegExp("  ",'g')," ");
w=w.replace(new RegExp("  ",'g')," ");
r=w.split(" ");
return(r);
}
//---------------------------------------------------------------------------
//---------------------------------------------------------------------------
// constants
var
   e2   = 1.4399764,   // MeV*fm
   h2_mn = 41.8015876123,
   degTorad = 0.01745329252,
   radTodeg = 57.2957795131;
//-------------------------------------------------------------------------------
// Zp,Ap - projectile charge and mass
// Zt,At - target charge and mass
// Ecm   - c.m. energy
// Npoint- numebr of points
// Angle - scattering angles in degree (array [Npoint])
// CS    - cross section (array [Npoint])
// dCS   - errors for cross section (array [Npoint])
// cm_or_lab - scattering angle measured in 'cm' or 'lab' system
// cs_unit   - cross section units ('b','mb','μb','nb','pb','σ/σ(Ruth)')
// dcs_unit  - cr. section errors units ('abs' or '%')
function CS_mb_to_Ratio(Zp,Ap,Zt,At,Ecm,Npoint,Angle,CS,dCS,cm_or_lab,cs_unit,dcs_unit)
{
	// poschitaem coeff. dlya perescheta v mb
	coeff = 1.0;
	if(cs_unit == 'b')      coeff = 1.0e-3;
	else if(cs_unit == 'mb')  coeff = 1.0;
	else if(cs_unit == 'μb') coeff = 1.0e3;
	else if(cs_unit == 'nb')  coeff = 1.0e6;
	else if(cs_unit == 'pb')  coeff = 1.0e9;
	
	if(1*Zp*Zt!=0)
	{
		for(i=0;i<Npoint;i++)
		{
			if(dcs_unit == 'abs')// snachala oshibki perevedem v %
			{
				s = 100./CS[i]; dCS[i] *= s;
			}
			 
			CS[i]=CS[i]/(getRutherford(Zp,Ap,Zt,At,Ecm,Angle[i],cm_or_lab) *coeff);
			
			if(dcs_unit == 'abs')// teper' oshibki perevedem v abs values
			{
				s = 0.01*CS[i]; dCS[i] *= s;
			}
		}
	}
}
//-------------------------------------------------------------------------------
function CS_Ratio_to_mb(Zp,Ap,Zt,At,Ecm,Npoint,Angle,CS,dCS,cm_or_lab,cs_unit,dcs_unit)
{
	// poschitaem coeff. dlya perescheta v mb
	coeff = 1.0;
	if(cs_unit == 'b')        coeff = 1.0e-3;
	else if(cs_unit == 'mb')  coeff = 1.0;
	else if(cs_unit == 'μb')  coeff = 1.0e3;
	else if(cs_unit == 'nb')  coeff = 1.0e6;
	else if(cs_unit == 'pb')  coeff = 1.0e9;

	if(1*Zp*Zt!=0)
	{
		for(i=0;i<Npoint;i++)
		{
			if(dcs_unit == 'abs')// snachala oshibki perevedem v %
			{
				dCS[i] = 100.0*dCS[i]/CS[i];
			}

			CS[i] = CS[i]*coeff*getRutherford(Zp,Ap,Zt,At,Ecm,Angle[i],cm_or_lab);

			if(dcs_unit == 'abs')// teper' oshibki perevedem v abs values
			{
				dCS[i] = 0.01*dCS[i]*CS[i];
			}
		}
	}
}
//-------------------------------------------------------------------------------
function getRutherford(Zp,Ap,Zt,At,Ecm,th_deg,cm_or_lab) 
{
	x  = 0.25*Zp*Zt*e2/Ecm;
	th = degTorad*th_deg;
	s  = Math.sin(0.5*th);//sin(th/2)
	Jacobian=1.;
	mu = 1.*Ap*At/(1.*Ap + 1.*At);
	eta = 1.*Zp*Zt*e2*mu/h2_mn/Math.sqrt(2.0*Ecm*mu/h2_mn);	
	//---------------------------------------------------------------------------
	// 1. v sluchae esli sechenie nuzhno v lab systeme, nuzhno perevesti Ruth. cr. section to the lab. system
	// 2. ugol th_deg v lab sisteme, => nuzhno ego pereschitat' v cm. dlya vychisleniya J*ds/dO
	if(cm_or_lab == 'lab')
	{
		cos2_lab = Math.cos(th);
		sin_lab = Math.sin(th);
		sin2_lab = 1.*sin_lab*sin_lab,
		cos_cm,	xx, m12 = 1.*Ap/At;
		cos2_lab = 1.*cos2_lab*cos2_lab;

		xx=1.*m12*m12*(sin2_lab-1.)*sin2_lab + 1.*cos2_lab;
		
		if(th_deg<90) cos_cm = (-1.*m12*sin2_lab + 1.*Math.sqrt(xx));
		else cos_cm = (-1.*m12*sin2_lab - 1.*Math.sqrt(xx));
		s = Math.sqrt(0.5*(1.-cos_cm));// sin(th_cm/2)

		xx=1.*1+1.*m12*m12+2.*m12*cos_cm;
		Jacobian = 1.*xx*Math.sqrt(xx)/Math.abs(1.+1.*m12*cos_cm);
	}
	//---------------------------------------------------------------------------
	if(Ap==At && Zp==Zt) 
	{
		sin2_th = 1.*s*s;
		cos2_th = 1 - 1.*sin2_th;
		ln_sin2_th = Math.log(1./sin2_th);
		ln_cos2_th = Math.log(1./cos2_th);
		re = 1.*Math.cos(eta*ln_sin2_th)/sin2_th + 1.*Math.cos(eta*ln_cos2_th)/cos2_th;
		im = 1.*Math.sin(eta*ln_sin2_th)/sin2_th + 1.*Math.sin(eta*ln_cos2_th)/cos2_th;
		cs = 10.*x*x*(re*re + im*im);
	}
	else cs = 10.*x*x/(s*s*s*s);
	//---------------------------------------------------------------------------
	return cs*Jacobian;
}
//-----------------------------------------------------------------------
// Aplf - mass of projectile in exit channel
// Atlf - mass of target in exit channel
// Q    - Q-value: Q = Ecm(exit) - Ecm(entrance)
function CS_lab_to_cm(Ap,At,Aplf,Atlf,Ecm,Q,Npoint,Angle,CS,dCS,cs_unit,dcs_unit)
{
	m12=1.*Ap/At,
	mu  = 1.*Ap*At/(1.*Ap + 1.*At);
	_mu = 1.*Aplf*Atlf/(1.*Aplf + 1.*Atlf);
	alpha = Math.sqrt(_mu*(1.0 + Q/Ecm)/mu);
	
	for(i=0;i<Npoint;i++)
	{
		s = Math.sin(Angle[i]*degTorad);
		c = Math.cos(Angle[i]*degTorad);
		
		chi = (1.*c*Math.sqrt(alpha*alpha - m12*m12*s*s) - 1.*m12*s*s)/alpha;
		
		Angle[i]=Math.acos(chi)*radTodeg;
		
		x = 1.*alpha*alpha + 1.*m12*m12 + 2.*alpha*m12*chi; 
		x = alpha*alpha*Math.abs(1.*alpha+1.*m12*chi)/(x*Math.sqrt(x));

		if(cs_unit != 'σ/σ(Ruth)') // esli sechenie v absolutnyh ed.
		{
			CS[i]=CS[i]*x;
			if(dcs_unit == 'abs') dCS[i]=dCS[i]*x;
		}
	}
}
function CS_cm_to_lab(Ap,At,Aplf,Atlf,Ecm,Q,Npoint,Angle,CS,dCS,cs_unit,dcs_unit)
{
	m12=1.*Ap/At;
	mu  = 1.*Ap*At/(1.*Ap + 1.*At);
	_mu = 1.*Aplf*Atlf/(1.*Aplf + 1.*Atlf);	
	alpha = Math.sqrt(_mu*(1.0+ Q/Ecm)/mu);

	for(i=0;i<Npoint;i++)
	{
		c = Math.cos(Angle[i]*degTorad);
		
		x = 1.*alpha*alpha + 1.*m12*m12 + 2.*alpha*m12*c; 
		thlab = (1.*m12+1.*alpha*c)/Math.sqrt(x);
		thlab = Math.acos(thlab);
		
		Angle[i] = thlab*radTodeg;
		x = x*Math.sqrt(x)/Math.abs(1.*alpha+1.*m12*c)/(alpha*alpha);

		if(cs_unit != 'σ/σ(Ruth)') // esli sechenie v absolutnyh ed.
		{
			CS[i]=CS[i]*x;
			if(dcs_unit == 'abs') dCS[i]=dCS[i]*x;
		}
	}
}
//-----------------------------------------------------------------------