SOLARIS_Web_Simulation/heliosmatics.html
2022-12-13 23:02:47 -05:00

621 lines
16 KiB
HTML

<!DOCTYPE html>
<html>
<head>
<title>Heliosmatics</title>
<meta name="description" content="HELIOSMATICS was first build by Ben P. Kay in MS Excel around 2010. Later, it was modified by Ryan Tang. Now, it migrate to web.">
<script src="https://cdn.plot.ly/plotly-2.16.1.min.js"></script>
</head>
<style>
body{
font-family: Arial, Helvetica, sans-serif;
}
</style>
<body>
<h1>HELIOSmatics (experimental...)</h1>
<h2 id='reactionName'> 24F(d,p)25F@10MeV/u</h2>
<table>
<tr>
<td style="text-align:right"> Beam (A):</td>
<td><Input type="text" size="2" value="24F" id="beam"/></td>
<td style="text-align:right"> Beam Ex:</td>
<td><Input type="text" size="2" value="0" id="beamEx"/></td>
<td>MeV</td>
</tr>
<tr>
<td style="text-align:right"> Target (a):</td>
<td><Input type="text" size="2" value="d" id="target"/></td>
</tr>
<tr>
<td style="text-align:right"> Light (b):</td>
<td><Input type="text" size="2" value="p" id="light"/></td>
<td style="text-align:right"> Q-value:</td>
<td id='Q-value'>2.057</td>
<td>MeV</td>
</tr>
<tr>
<td style="text-align:right"> Heavy (B):</td>
<td id='heavyName'>25F</td>
</tr>
</table>
<p></p>
<table>
<tr>
<td>
<input type="radio" name="SSType" id='HELIOS' value="HELIOS"/>HELIOS
</td>
<td>
<input type="radio" name="SSType" id='SOLARIS' value="SOLARIS" checked="checked"/>SOLARIS
</td>
<td>
<input type="radio" name="SSType" id='ISS' value="ISS"/>ISS
</td>
</tr>
</table>
<p></p>
<table id="setting">
<table>
<tr>
<td style="text-align:right"> B-field:</td>
<td><Input type="text" size="1" value="-2" id='BField'/></td>
<td>T</td>
<td><Input type="range" min="-4" max="4" step="0.05" value="-2" style="width:400px" id='BRange'/> </td>
</tr>
<tr>
<td style="text-align:right"> Beam Energy:</td>
<td><Input type="text" size="1" value="10" id='KEA'/></td>
<td>MeV/u</td>
<td><Input type="range" min="0" max="20" step="0.1" value="10" style="width:400px" id='KEARange'/> </td>
</tr>
<tr>
<td style="text-align:right"> min Beam Energy:</td>
<td id='minKEA'>-1.107</td>
<td>MeV/u</td>
</tr>
</table>
<table id="ExTable", style="border:1px solid">
<tr>
<td> E<sub>x</sub> [MeV]</td>
<td> &theta;<sub>cm</sub>[deg]</td>
<td>E<sub>b</sub>[MeV]</td>
<td>Z<sub>b0</sub>[mm]</td>
<td>Z<sub>b</sub>[mm]</td>
<td>2*&rho;<sub>b</sub>[mm]</td>
<td>&theta;Lab<sub>b</sub>[deg]</td>
<td>T<sub>b</sub>[ns]</td>
<td>E<sub>B</sub>[MeV]</td>
<td>&theta;Lab<sub>B</sub>[deg]</td>
<td>Z<sub>B0</sub>/2[mm]</td>
<td>2*&rho;<sub>B</sub>[mm]</td>
</tr>
</tr>
<td><input type="text" id='Ex1' name="Ex" size="7" value="0"/></td>
<td><input type="text" id='theta1' name="thetaCM" size="7" value="10"/></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
<td></td>
</tr>
<tr>
<td><button type="button" onclick="addRow()" style="width:85px">Add E<sub>x</sub></button></td>
<td><button type="button" onclick="deleteRow()">Remove E<sub>x</sub></button></td>
</tr>
</table>
<p></p>
<table>
<tr>
<td style="text-align:right"> &theta;<sub>CM</sub>:</td>
<td><Input type="text" size="1" value="0" id='thetaCM'/></td>
<td>deg</td>
<td><Input type="range" min="0" max="40" step="0.1" value="0" style="width:600px" id='thetaCMRange'/> </td>
</tr>
<td style="text-align:right"> Array Pos:</td>
<td><Input type="text" size="1" value="0" id='posArray'/></td>
<td>mm</td>
<td><Input type="range" min="-500" max="1000" step="1" value="0" style="width:600px" id='posArrayRange'/> </td>
</tr>
</table>
<p></p>
<div id="Plot_EZ" style="width:650px;height:600px;"></div>
<p id="n0"></p>
<p id="n1"></p>
<p id="n2"></p>
<p id="n3"></p>
<!--
<p id='beam_k_lab'></p>
<p id='E_tot_cm'></p>
<p id='KE_cm'></p>
<p id='max_Ex'></p>
<p id='beta'></p>
<p id='gamma'></p>
<p id='msg'></p>
-->
</body>
<script>
function addRow() {
var table = document.getElementById("ExTable");
var nRow = table.rows.length;
var row = table.insertRow(nRow-1);
row.innerHTML = '<td><input type="text" name="Ex" size="7" value=\"' + (nRow-2) + '\" id=\"Ex' + (nRow-1) + '\"/></td> \
<td><input type="text" name="thetaCM" size="7" value=\"10\" id=\"theta' + (nRow-1) + '\"/></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td> \
<td></td>';
Calculation();
}
function deleteRow(){
var table = document.getElementById("ExTable");
var nRow = table.rows.length;
if ( nRow > 3){
table.deleteRow(nRow-2);
}
Calculation();
}
var beam = []; //A, Z, Mass, Name
var beamMass;
var target= [];
var light=[];
var heavy=[];
var beamEx;
var BField;
var KEA;
var KE;
var reactionName;
var Qvalue;
var minKEA;
var perpDistant = 11.5; //mm, detector prepdicular distance
var bore = 462.0; // mm
const c = 299.792468;
var beam_k_lab;
var E_tot_cm ;
var KE_cm ;
var max_Ex ;
var beta ;
var gamma ;
var ez_slope ; // MeV/mm
var alpha ;
var alpha_B ;
var ExList = [];
function GetMassFromSym(AZ, id){
let str = 'massProxy.py?AZ=' + AZ;
let client = new XMLHttpRequest();
client.onreadystatechange = function() {
let mass = client.responseText.split(",");
if( id == 0 ){
beam[0] = parseInt(mass[0]);
beam[1] = parseInt(mass[1]);
beam[2] = parseFloat(mass[2]);
beam[3] = AZ;
}
if( id == 1 ){
target[0] = parseInt(mass[0]);
target[1] = parseInt(mass[1]);
target[2] = parseFloat(mass[2]);
target[3] = AZ;
}
if( id == 2 ){
light[0] = parseInt(mass[0]);
light[1] = parseInt(mass[1]);
light[2] = parseFloat(mass[2]);
light[3] = AZ;
}
}
client.open('GET', str, false);
client.send();
}
function GetMassFromAZ(A,Z){
let str = 'massProxy.py?A=' + A + '&Z=' + Z;
let client = new XMLHttpRequest();
client.onreadystatechange = function() {
let mass = client.responseText.split(",");
heavy[2] = parseFloat(mass[2]);
heavy[3] = mass[3].replace(/\s/g, '');
}
client.open('GET', str, false);
client.send();
}
function GetMass(){
let beam_AZ = document.getElementById('beam').value;
GetMassFromSym(beam_AZ, 0);
let target_AZ = document.getElementById('target').value;
GetMassFromSym(target_AZ, 1);
let light_AZ = document.getElementById('light').value;
GetMassFromSym(light_AZ, 2);
beamMass = beam[2];
heavy[0] = beam[0]+target[0]-light[0];
heavy[1] = beam[1]+target[1]-light[1];
GetMassFromAZ(heavy[0], heavy[1]);
document.getElementById('heavyName').innerHTML = heavy[3];
//document.getElementById('n0').innerHTML = beam[0] + "," + beam[1] + "," + beam[2]
//document.getElementById('n1').innerHTML = target[0] + "," + target[1] + "," + target[2]
//document.getElementById('n2').innerHTML = light[0] + "," + light[1] + "," + light[2]
//document.getElementById('n3').innerHTML = heavy[0] + "," + heavy[1] + "," + heavy[2]
}
function CalConstants(){
beamEx = parseFloat(document.getElementById('beamEx').value);
beam[2] = beamMass + beamEx;
BField = parseFloat(document.getElementById('BField').value);
KEA = document.getElementById('KEA').value;
KE = KEA * beam[0];
reactionName = beam[3] +"(" + target[3] + "," + light[3] + ")" + heavy[3] + "@" + KEA + "MeV/u";
Qvalue = - heavy[2] - light[2] + beam[2] + target[2] ;
minKEA = (Math.pow(light[2] + heavy[2],2) - Math.pow(beam[2] + target[2],2))/2/target[2]/beam[0];
document.getElementById('reactionName').innerHTML = reactionName;
document.getElementById('minKEA').innerHTML = minKEA.toFixed(3);
document.getElementById('Q-value').innerHTML = Qvalue.toFixed(3);
beam_k_lab = Math.sqrt(Math.pow(beam[2] + KE,2) - Math.pow(beam[2],2));
E_tot_cm = Math.sqrt(Math.pow(target[2] + beam[2],2) + 2*target[2]*KE);
KE_cm = E_tot_cm - beam[2] - target[2];
max_Ex = KE_cm + Qvalue;
beta = beam_k_lab/(beam[2] + target[2] + KE);
gamma = 1./Math.sqrt(1-beta*beta);
ez_slope = Math.abs(BField) * c * light[1]*beta/2/Math.PI/1000; // MeV/mm
alpha = ez_slope/beta;
alpha_B = alpha * heavy[1]/light[1];
/*
document.getElementById('beam_k_lab').innerHTML = beam_k_lab;
document.getElementById('E_tot_cm').innerHTML = E_tot_cm;
document.getElementById('max_Ex').innerHTML = max_Ex;
document.getElementById('KE_cm').innerHTML = KE_cm;
document.getElementById('beta').innerHTML = beta;
document.getElementById('gamma').innerHTML = gamma;
*/
}
function SetSSType(){
let helios = document.getElementById('HELIOS').checked;
let solaris = document.getElementById('SOLARIS').checked;
let iss = document.getElementById('ISS').checked;
if ( helios == true ) perpDistant = 11.5;
if ( solaris == true ) perpDistant = 11.5;
if ( iss == true ) perpDistant = 20;
//document.getElementById('n1').innerHTML = perpDistant;
}
var xList =[];
var yList =[];
function Calculation(){
let tableEx = document.getElementById("ExTable");
let nRow = tableEx.rows.length;
xList = [];
yList = [];
ExList = [];
//alert("Calculation called");
let xmin ;
for( let i = 1; i < nRow-1; i++){
let Ex = parseFloat(document.getElementById("Ex" + i).value);
let theta = parseFloat(document.getElementById("theta" + i).value);
ExList.push(Ex);
//alert( i, ", Ex : " + Ex);
let haha1 = E_tot_cm*E_tot_cm - Math.pow(heavy[2] + Ex + light[2],2);
let haha2 = E_tot_cm*E_tot_cm - Math.pow(heavy[2] + Ex - light[2],2);
let k_cm = Math.sqrt(haha1*haha2)/2/E_tot_cm;
let cs = Math.cos(theta*Math.PI/180.);
let ss = Math.sin(theta*Math.PI/180.);
let qb = Math.sqrt(light[2]*light[2]+k_cm*k_cm);
let Eb = gamma * qb - gamma * beta * k_cm * cs - light[2];
let Zb0 = (gamma*beta* qb - gamma * k_cm * cs )/alpha; //mm
if( i == 1 ) xmin = Zb0;
if( i > 1 & Zb0 < xmin) xmin = Zb0;
let rho = k_cm * ss/c/light[1]/Math.abs(BField) * 1000; // mm
let Zb = Zb0 * (1- Math.asin(perpDistant/rho)/2/Math.PI);
let thetaLab = 180 - Math.atan2(k_cm * ss, gamma * (beta * qb - k_cm * cs)) * 180/Math.PI;
let Tcyc = (light[2] + Eb)/alpha/c;
tableEx.rows[i].cells[2].innerHTML = Eb.toFixed(3);
tableEx.rows[i].cells[3].innerHTML = Zb0.toFixed(1);
tableEx.rows[i].cells[4].innerHTML = Zb.toFixed(1);
tableEx.rows[i].cells[5].innerHTML = (2*rho).toFixed(1);
tableEx.rows[i].cells[6].innerHTML = thetaLab.toFixed(2);
tableEx.rows[i].cells[7].innerHTML = Tcyc.toFixed(2);
let qB = Math.sqrt(heavy[2]*heavy[2]+k_cm*k_cm);
let EB = gamma * qB + gamma * beta * k_cm * cs - heavy[2];
let thetaLab_B = Math.atan2(-k_cm * ss, gamma * (beta * qB + k_cm * cs)) * 180/Math.PI;
let ZB0 = (gamma*beta* qB + gamma * k_cm * cs )/alpha_B; //mm
let rhoB = k_cm * ss/c/heavy[1]/Math.abs(BField) * 1000; // mm
tableEx.rows[i].cells[8].innerHTML = EB.toFixed(3);
tableEx.rows[i].cells[9].innerHTML = thetaLab_B.toFixed(2);
tableEx.rows[i].cells[10].innerHTML = (ZB0/2).toFixed(2);
tableEx.rows[i].cells[11].innerHTML = (2*rhoB).toFixed(2);
let xxx = [];
let yyy = [];
for( let j = 0; j < 100 ; j++){
let thetaCM = Math.PI/Math.log10(100) * Math.log10(100/(100-j)) ;
let temp = 2* Math.PI * ez_slope / beta / k_cm * perpDistant / Math.sin(thetaCM);
if( !isFinite(temp) ) continue;
let pxTemp = beta /ez_slope * (gamma * beta * qb - gamma * k_cm * Math.cos(thetaCM)) * (1 - Math.asin(temp)/Math.PI/2) ;
let pyTemp = gamma * qb - light[2] - gamma * beta * k_cm * Math.cos(thetaCM);
if( !isNaN(pxTemp) && !isNaN(pyTemp) ){
xxx.push(pxTemp);
yyy.push(pyTemp);
}
};
xList.push(xxx);
yList.push(yyy);
};
//var xRange = xmin < 0 ? [xmin, 0] : [xmin, xmin + 600]
let xRange = [-700, 0];
Plot(xRange, [0, 12]);
}
function PlotThetaCMLine(thetaCM){
let cs = Math.cos(thetaCM * Math.PI /180);
let ss = Math.sin(thetaCM * Math.PI /180);
let zList = [];
let eList = [];
for( let z = -2000; z < 2000; z+=30){
zList.push(z);
let kaka = Math.pow(gamma * ss,2);
let a1 = light[2]*light[2]*(1-kaka);
let haha = (cs*Math.sqrt(alpha*alpha*z*z + a1) - kaka * ez_slope * z)/(1-kaka)- light[2];
eList.push(haha);
}
Plotly.addTraces("Plot_EZ", {x : zList, y: eList, name:"thetaCM=" + thetaCM }, 0);
}
function PlotBore(){
let zList = [];
let eList = [];
for( let z = -2000; z < 2000; z+=30){
zList.push(z);
let haha = Math.sqrt((z*z+Math.PI*Math.PI*bore*bore)*alpha*alpha + light[2]*light[2]) - light[2];
eList.push(haha);
}
Plotly.addTraces("Plot_EZ", {x : zList, y: eList, name:"Bore" }, 0);
}
function Plot(rangeX, rangeY){
let nEx = xList.length;
//document.getElementById('n0').innerHTML = nEx;
var data = [
{ x : xList[0], y : yList[0], mode:"lines", type:"scatter", name:"Ex="+ExList[0]}
];
var layout = {
xaxis: {range: rangeX, title: "Z [mm]", mirror : "allticks", linewidth : "1"},
yaxis: {range: rangeY, title: "Energy [MeV]" , mirror : "allticks", linewidth : "1"},
title: reactionName,
legend: {yanchor:"top", xanchor:"left", x:"0.01",y:"0.99" }
};
Plotly.newPlot( "Plot_EZ", data, layout );
for( let i = 1; i < nEx; i++){
Plotly.addTraces("Plot_EZ", {x : xList[i], y: yList[i], name:"Ex="+ExList[i]});
}
PlotThetaCMLine(document.getElementById('thetaCM').value);
PlotBore();
}
document.getElementById('beam').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
GetMass();
CalConstants();
}
}, false
);
document.getElementById('target').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
GetMass();
CalConstants();
}
}, false
);
document.getElementById('light').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
GetMass();
CalConstants();
}
}, false
);
document.getElementById('BField').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('BRange').value = this.value;
//CalConstants();
}
}, false
);
document.getElementById('beamEx').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
CalConstants();
}
}, false
);
document.getElementById('KEA').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('KEARange').value = this.value;
//CalConstants();
}
}, false
);
document.getElementById('thetaCM').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('thetaCMRange').value = this.value;
}
}, false
);
document.getElementById('posArray').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('posArrayRange').value = this.value;
}
}, false
);
window.logMeThis = function(){
SetSSType();
Calculation();
}
var FuncEx = window.logMeThis.bind(null, "Ex");
window.addEventListener('keypress', FuncEx);
var FuncThetaCM = window.logMeThis.bind(null, "thetaCM");
window.addEventListener('keypress', FuncThetaCM);
window.checkSSType = function(){
SetSSType();
Calculation();
}
var FuncCheckSSType = window.checkSSType.bind(null, "SSType");
document.getElementById('HELIOS').addEventListener('click', FuncCheckSSType);
document.getElementById('SOLARIS').addEventListener('click', FuncCheckSSType);
document.getElementById('ISS').addEventListener('click', FuncCheckSSType);
document.getElementById('BRange').oninput = function(){
document.getElementById('BField').value = this.value;
CalConstants();
Calculation();
}
document.getElementById('KEARange').oninput = function(){
document.getElementById('KEA').value = this.value;
CalConstants();
Calculation();
}
document.getElementById('thetaCMRange').oninput = function(){
document.getElementById('thetaCM').value = this.value;
CalConstants();
Calculation();
}
document.getElementById('posArrayRange').oninput = function(){
document.getElementById('posArray').value = this.value;
}
GetMass();
CalConstants();
SetSSType();
Calculation();
</script>
<html>