2022-12-13 18:24:26 -05:00
<!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;
}
2022-12-13 23:02:47 -05:00
2022-12-13 18:24:26 -05:00
< / style >
< body >
2022-12-14 14:31:54 -05:00
< h1 > HELIOSmatics< / h1 >
2022-12-14 16:25:37 -05:00
< h1 id = 'reactionName' style = "color: #1363A7" > 24F(d,p)25F@10MeV/u< / h1 >
2022-12-13 23:02:47 -05:00
< table >
< tr >
< td style = "text-align:right" > Beam (A):< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "24F" id = "beam" / > < / td >
2022-12-13 23:02:47 -05:00
< td style = "text-align:right" > Beam Ex:< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "0" id = "beamEx" / > < / td >
2022-12-13 23:02:47 -05:00
< td > MeV< / td >
< / tr >
< tr >
< td style = "text-align:right" > Target (a):< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "d" id = "target" / > < / td >
2022-12-13 23:02:47 -05:00
< / tr >
< tr >
< td style = "text-align:right" > Light (b):< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "p" id = "light" / > < / td >
2022-12-13 23:02:47 -05:00
< 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 >
2022-12-13 18:24:26 -05:00
< 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 >
2022-12-13 23:02:47 -05:00
< table >
2022-12-13 18:24:26 -05:00
< tr >
2022-12-13 23:02:47 -05:00
< td style = "text-align:right" > B-field:< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "-2" id = 'BField' / > < / td >
2022-12-13 23:02:47 -05:00
< td > T< / td >
< td > < Input type = "range" min = "-4" max = "4" step = "0.05" value = "-2" style = "width:400px" id = 'BRange' / > < / td >
2022-12-13 18:24:26 -05:00
< / tr >
2022-12-13 23:02:47 -05:00
2022-12-13 18:24:26 -05:00
< tr >
< td style = "text-align:right" > Beam Energy:< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "10" id = 'KEA' / > < / td >
2022-12-13 18:24:26 -05:00
< td > MeV/u< / td >
2022-12-13 23:02:47 -05:00
< td > < Input type = "range" min = "0" max = "20" step = "0.1" value = "10" style = "width:400px" id = 'KEARange' / > < / td >
2022-12-13 18:24:26 -05:00
< / tr >
< tr >
2022-12-13 23:02:47 -05:00
< td style = "text-align:right" > min Beam Energy:< / td >
2022-12-13 18:24:26 -05:00
< td id = 'minKEA' > -1.107< / td >
< td > MeV/u< / td >
< / tr >
< / table >
2022-12-14 14:31:54 -05:00
< p > < / p >
< table id = "ExTable" , style = "border:1px solid; text-align:center;" >
2022-12-13 18:24:26 -05:00
< tr >
2022-12-14 14:31:54 -05:00
< th style = "width:85px" > E< sub > x< / sub > [MeV]< / th >
< th style = "width:85px" > θ < sub > cm< / sub > [deg]< / th >
< th style = "width:70px" > E< sub > b< / sub > [MeV]< / th >
< th style = "width:70px" > Z< sub > b0< / sub > [mm]< / th >
< th style = "width:70px" > Z< sub > b< / sub > [mm]< / th >
< th style = "width:70px" > 2ρ < sub > b< / sub > [mm]< / th >
< th style = "width:70px" > θ Lab< sub > b< / sub > [deg]< / th >
< th style = "width:60px" > T< sub > b< / sub > [ns]< / th >
< th style = "width:70px" > E< sub > B< / sub > [MeV]< / th >
< th style = "width:90px" > θ Lab< sub > B< / sub > [deg]< / th >
< th style = "width:80px" > Z< sub > B0< / sub > /2[mm]< / th >
< th style = "width:70px" > 2*ρ < sub > B< / sub > [mm]< / th >
2022-12-13 18:24:26 -05:00
< / tr >
< / tr >
2022-12-14 14:31:54 -05:00
< td > < input type = "text" id = 'Ex1' name = "Ex" size = "8" value = "0" / > < / td >
< td > < input type = "text" id = 'theta1' name = "thetaCM" size = "8" value = "10" / > < / td >
2022-12-13 18:24:26 -05:00
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< / tr >
2022-12-14 16:59:17 -05:00
< / tr >
< td > < input type = "text" id = 'Ex2' name = "Ex" size = "8" value = "0" / > < / td >
< td > < input type = "text" id = 'theta2' name = "thetaCM" size = "8" value = "40" / > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< td > < / td >
< / tr >
2022-12-13 18:24:26 -05:00
< tr >
2022-12-13 23:02:47 -05:00
< td > < button type = "button" onclick = "addRow()" style = "width:85px" > Add E< sub > x< / sub > < / button > < / td >
2022-12-13 18:24:26 -05:00
< td > < button type = "button" onclick = "deleteRow()" > Remove E< sub > x< / sub > < / button > < / td >
< / tr >
< / table >
2022-12-13 23:02:47 -05:00
< p > < / p >
< table >
< tr >
< td style = "text-align:right" > θ < sub > CM< / sub > :< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "0" id = 'thetaCM' / > < / td >
2022-12-13 23:02:47 -05:00
< td > deg< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "range" min = "0" max = "50" step = "0.1" value = "0" style = "width:600px" id = 'thetaCMRange' / > < / td >
2022-12-13 23:02:47 -05:00
< / tr >
< td style = "text-align:right" > Array Pos:< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "text" size = "3" value = "-100" id = 'posArray' / > < / td >
2022-12-13 23:02:47 -05:00
< td > mm< / td >
2022-12-14 14:31:54 -05:00
< td > < Input type = "range" min = "-500" max = "1000" step = "1" value = "-100" style = "width:600px" id = 'posArrayRange' / > < / td >
2022-12-13 23:02:47 -05:00
< / tr >
2022-12-14 16:25:37 -05:00
< tr >
< td style = "text-align:right" > Recoil Pos:< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "text" size = "3" value = "1500" id = 'posRecoil' / > < / td >
2022-12-14 16:25:37 -05:00
< td > mm< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "range" min = "0" max = "2000" step = "1" value = "1500" style = "width:600px" id = 'posRecoilRange' / > < / td >
2022-12-14 16:25:37 -05:00
< / tr >
2022-12-13 23:02:47 -05:00
< / table >
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
< p > < / p >
2022-12-13 18:24:26 -05:00
2022-12-14 16:25:37 -05:00
< table >
< tr >
< td > < / td >
< td style = "text-align:left" >
< table >
< tr >
< td > zRange(Min):< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "text" size = "4" value = "-200" id = 'zRange1' / > < / td >
2022-12-14 16:25:37 -05:00
< td > mm< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "range" min = "-2000" max = "4000" step = "1" value = "-200" style = "width:300px" id = 'zRange1Slider' / > < / td >
2022-12-14 16:25:37 -05:00
< / tr >
< tr >
< td > zRange(Max):< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "text" size = "4" value = "2000" id = 'zRange2' / > < / td >
2022-12-14 16:25:37 -05:00
< td > mm< / td >
2022-12-14 16:59:17 -05:00
< td > < Input type = "range" min = "-2000" max = "4000" step = "1" value = "2000" style = "width:300px" id = 'zRange2Slider' / > < / td >
2022-12-14 16:25:37 -05:00
< / tr >
< tr >
< td > rRange:< / td >
< td > < Input type = "text" size = "4" value = "50" id = 'rRange' / > < / td >
< td > mm< / td >
< td > < Input type = "range" min = "1" max = "400" step = "1" value = "50" style = "width:300px" id = 'rRangeSlider' / > < / td >
< / tr >
< /table
< / td >
< tr >
< td >
< div id = "Plot_EZ" style = "width:650px;height:600px;" > < / div >
< / td >
< td >
< div id = "Plot_RZ" style = "width:650px;height:600px;" > < / div >
< / td >
< / tr >
< / table >
2022-12-13 18:24:26 -05:00
< p id = "n0" > < / p >
< p id = "n1" > < / p >
< p id = "n2" > < / p >
< p id = "n3" > < / p >
< / body >
2022-12-14 14:31:54 -05:00
<!-- ##################################### -->
2022-12-13 18:24:26 -05:00
< script >
function addRow() {
2022-12-14 16:59:17 -05:00
let table = document.getElementById("ExTable");
let nRow = table.rows.length;
let row = table.insertRow(nRow-1);
2022-12-14 16:25:37 -05:00
row.innerHTML = '< td > < input type = "text" name = "Ex" size = "8" value = \"' + ( nRow-2 ) + ' \ " id = \"Ex' + ( nRow-1 ) + ' \ " / > < / td > \
< td > < input type = "text" name = "thetaCM" size = "8" value = \"10\" id = \"theta' + ( nRow-1 ) + ' \ " / > < / td > \
2022-12-13 18:24:26 -05:00
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > \
< td > < / td > ';
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
2022-12-14 16:59:17 -05:00
AdjustRangeEZ(document.getElementById('posArray').value);
2022-12-13 18:24:26 -05:00
}
function deleteRow(){
2022-12-14 16:59:17 -05:00
let table = document.getElementById("ExTable");
let nRow = table.rows.length;
2022-12-13 18:24:26 -05:00
if ( nRow > 3){
table.deleteRow(nRow-2);
}
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
2022-12-13 18:24:26 -05:00
}
2022-12-13 23:02:47 -05:00
var beam = []; //A, Z, Mass, Name
var beamMass;
2022-12-13 18:24:26 -05:00
var target= [];
var light=[];
var heavy=[];
2022-12-13 23:02:47 -05:00
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
2022-12-14 14:31:54 -05:00
var arrayLen = 50 * 10 + 2 * 9 // SOLARIS
var detLen = 50; // SOLARIS
var nDet = 10; // SOLARIS
var detGap = 2; // SOLARIS
2022-12-13 23:02:47 -05:00
2022-12-14 14:31:54 -05:00
const c = 299.792468; // mm/ns
2022-12-13 23:02:47 -05:00
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 ;
2022-12-14 16:25:37 -05:00
var xList =[]; // for E-Z plot
var yList =[]; // for E-Z plot
2022-12-13 23:02:47 -05:00
var ExList = [];
2022-12-14 16:25:37 -05:00
var Zb0List = [];
var rhoList = [];
var ZB0List = [];
var rhoBList = [];
2022-12-14 14:31:54 -05:00
2022-12-14 16:59:17 -05:00
var xRange ;
2022-12-14 14:31:54 -05:00
var yRange = [0, 12];
2022-12-13 23:02:47 -05:00
2022-12-14 16:25:37 -05:00
var zList = []; // for R-Z plot
2022-12-14 16:59:17 -05:00
var rbList = []; // for R-Z plot
var rBList = []; // for R-Z plot
2022-12-14 16:25:37 -05:00
2022-12-13 23:02:47 -05:00
function GetMassFromSym(AZ, id){
let str = 'massProxy.py?AZ=' + AZ;
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
let client = new XMLHttpRequest();
2022-12-13 18:24:26 -05:00
client.onreadystatechange = function() {
2022-12-13 23:02:47 -05:00
let mass = client.responseText.split(",");
2022-12-13 18:24:26 -05:00
if( id == 0 ){
beam[0] = parseInt(mass[0]);
beam[1] = parseInt(mass[1]);
beam[2] = parseFloat(mass[2]);
2022-12-13 23:02:47 -05:00
beam[3] = AZ;
2022-12-13 18:24:26 -05:00
}
if( id == 1 ){
target[0] = parseInt(mass[0]);
target[1] = parseInt(mass[1]);
target[2] = parseFloat(mass[2]);
2022-12-13 23:02:47 -05:00
target[3] = AZ;
2022-12-13 18:24:26 -05:00
}
if( id == 2 ){
light[0] = parseInt(mass[0]);
light[1] = parseInt(mass[1]);
light[2] = parseFloat(mass[2]);
2022-12-13 23:02:47 -05:00
light[3] = AZ;
2022-12-13 18:24:26 -05:00
}
}
client.open('GET', str, false);
client.send();
}
2022-12-13 23:02:47 -05:00
function GetMassFromAZ(A,Z){
let str = 'massProxy.py?A=' + A + '& Z=' + Z;
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
let client = new XMLHttpRequest();
2022-12-13 18:24:26 -05:00
client.onreadystatechange = function() {
2022-12-13 23:02:47 -05:00
let mass = client.responseText.split(",");
2022-12-13 18:24:26 -05:00
heavy[2] = parseFloat(mass[2]);
heavy[3] = mass[3].replace(/\s/g, '');
}
client.open('GET', str, false);
client.send();
}
2022-12-13 23:02:47 -05:00
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);
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
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]
2022-12-13 18:24:26 -05:00
}
2022-12-13 23:02:47 -05:00
function CalConstants(){
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
beamEx = parseFloat(document.getElementById('beamEx').value);
beam[2] = beamMass + beamEx;
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
BField = parseFloat(document.getElementById('BField').value);
KEA = document.getElementById('KEA').value;
KE = KEA * beam[0];
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
reactionName = beam[3] +"(" + target[3] + "," + light[3] + ")" + heavy[3] + "@" + KEA + "MeV/u";
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
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;
2022-12-13 18:24:26 -05:00
document.getElementById('minKEA').innerHTML = minKEA.toFixed(3);
document.getElementById('Q-value').innerHTML = Qvalue.toFixed(3);
2022-12-13 23:02:47 -05:00
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);
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
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];
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
}
function SetSSType(){
let helios = document.getElementById('HELIOS').checked;
let solaris = document.getElementById('SOLARIS').checked;
let iss = document.getElementById('ISS').checked;
2022-12-14 14:31:54 -05:00
if ( helios == true ) {
perpDistant = 11.5;
detGap = 5;
detLen = 50;
nDet = 6;
}
if ( solaris == true ) {
perpDistant = 11.5;
detGap = 5;
detLen = 50;
nDet = 10;
}
if ( iss == true ) {
perpDistant = 20;
detGap = 5;
detLen = 100;
nDet = 4;
}
arrayLen = detLen * nDet + detGap * (nDet-1);
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
//document.getElementById('n1').innerHTML = perpDistant;
}
2022-12-13 18:24:26 -05:00
2022-12-14 16:25:37 -05:00
function CalculateEZ(){
2022-12-13 23:02:47 -05:00
let tableEx = document.getElementById("ExTable");
let nRow = tableEx.rows.length;
xList = [];
yList = [];
ExList = [];
2022-12-14 16:25:37 -05:00
Zb0List = [];
rhoList = [];
ZB0List = [];
rhoBList = [];
2022-12-14 16:59:17 -05:00
//alert("CalculateEZ called, nRow = " + nRow);
2022-12-13 18:24:26 -05:00
for( let i = 1; i < nRow-1 ; i + + ) {
let Ex = parseFloat(document.getElementById("Ex" + i).value);
let theta = parseFloat(document.getElementById("theta" + i).value);
2022-12-13 23:02:47 -05:00
ExList.push(Ex);
//alert( i, ", Ex : " + Ex);
2022-12-13 18:24:26 -05:00
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
2022-12-14 16:25:37 -05:00
Zb0List.push(Zb0);
2022-12-13 23:02:47 -05:00
let rho = k_cm * ss/c/light[1]/Math.abs(BField) * 1000; // mm
2022-12-14 16:25:37 -05:00
rhoList.push(rho);
2022-12-13 18:24:26 -05:00
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
2022-12-13 23:02:47 -05:00
let rhoB = k_cm * ss/c/heavy[1]/Math.abs(BField) * 1000; // mm
2022-12-14 16:25:37 -05:00
ZB0List.push(ZB0);
rhoBList.push(rhoB);
2022-12-13 18:24:26 -05:00
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);
2022-12-13 23:02:47 -05:00
let xxx = [];
let yyy = [];
2022-12-13 18:24:26 -05:00
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);
};
2022-12-14 16:59:17 -05:00
let haha = parseFloat(document.getElementById('posArrayRange').value);
let xStart = (haha < 0 ? haha - arrayLen - 100 : haha - 100 ) ;
let xEnd = (haha < 0 ? haha + 100: haha + arrayLen + 100 ) ;
xRange = [xStart, xEnd];
//document.getElementById('n0').innerHTML = xRange;
2022-12-14 16:25:37 -05:00
PlotEZ(xRange, yRange);
2022-12-13 18:24:26 -05:00
2022-12-14 14:31:54 -05:00
}
2022-12-14 16:25:37 -05:00
function PlotEZ(rangeX, rangeY){
2022-12-14 14:31:54 -05:00
SetSSType();
Plotly.purge("Plot_EZ");
2022-12-14 16:59:17 -05:00
2022-12-14 14:31:54 -05:00
let nEx = xList.length;
2022-12-14 16:25:37 -05:00
let data = [];
for( let i = 0; i < nEx ; i + + ) {
let newData = {
x : xList[i], y : yList[i], mode:"lines", type:"scatter", name:"Ex="+ExList[i]
}
data.push(newData);
}
2022-12-14 14:31:54 -05:00
let 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 );
2022-12-14 16:25:37 -05:00
2022-12-14 14:31:54 -05:00
PlotThetaCMLine(document.getElementById('thetaCM').value);
PlotBore();
2022-12-14 16:25:37 -05:00
AdjustRangeEZ(-100);
PlotRZ();
2022-12-13 23:02:47 -05:00
}
function PlotThetaCMLine(thetaCM){
let cs = Math.cos(thetaCM * Math.PI /180);
let ss = Math.sin(thetaCM * Math.PI /180);
2022-12-14 16:25:37 -05:00
let zzList = [];
2022-12-13 23:02:47 -05:00
let eList = [];
for( let z = -2000; z < 2000 ; z + = 30 ) {
2022-12-14 16:25:37 -05:00
zzList.push(z);
2022-12-13 23:02:47 -05:00
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);
2022-12-13 18:24:26 -05:00
}
2022-12-14 16:25:37 -05:00
Plotly.addTraces("Plot_EZ", {x : zzList, y: eList, name:"thetaCM=" + thetaCM }, 0);
2022-12-13 23:02:47 -05:00
}
function PlotBore(){
2022-12-14 16:25:37 -05:00
let zzList = [];
2022-12-13 23:02:47 -05:00
let eList = [];
for( let z = -2000; z < 2000 ; z + = 30 ) {
2022-12-14 16:25:37 -05:00
zzList.push(z);
2022-12-13 23:02:47 -05:00
let haha = Math.sqrt((z*z+Math.PI*Math.PI*bore*bore)*alpha*alpha + light[2]*light[2]) - light[2];
eList.push(haha);
2022-12-14 16:25:37 -05:00
}
Plotly.addTraces("Plot_EZ", {x : zzList, y: eList, name:"Bore" }, 0);
}
function CalculateRZ(){
// this rquire ZB0List and rhoBList from CalculateEZ();
zList = [];
2022-12-14 16:59:17 -05:00
rbList = [];
rBList = [];
2022-12-14 16:25:37 -05:00
for( let z = -2000; z < 4000 ; z + = 2 ) zList . push ( z ) ;
// cal the heavy recoil first
for( let i = 0; i < ZB0List.length ; i + + ) {
let rrr = []
for( let j = 0; j < zList.length ; j + + ) {
if( zList[j] < 0 ) {
rrr.push(NaN);
}else{
rrr.push(2*rhoBList[i] *Math.abs( Math.sin(zList[j] * Math.PI / ZB0List[i])));
}
}
2022-12-14 16:59:17 -05:00
rBList.push(rrr);
}
// cal the light recoil first
for( let i = 0; i < Zb0List.length ; i + + ) {
let rrr = []
for( let j = 0; j < zList.length ; j + + ) {
if( target[0] > light[0] ){
if( zList[j] > 0 ){
rrr.push(NaN);
}else{
rrr.push(2*rhoList[i] *Math.abs( Math.sin(zList[j] * Math.PI / Zb0List[i])));
}
}else{
if( zList[j] < 0 ) {
rrr.push(NaN);
}else{
rrr.push(2*rhoList[i] *Math.abs( Math.sin(zList[j] * Math.PI / Zb0List[i])));
}
}
}
rbList.push(rrr);
2022-12-13 23:02:47 -05:00
}
2022-12-14 16:25:37 -05:00
PlotRZ();
AdjustRecoilPos();
}
function PlotRZ(){
Plotly.purge("Plot_RZ");
let data = [];
2022-12-14 16:59:17 -05:00
let nEx = ExList.length;
2022-12-14 16:25:37 -05:00
for(let i = 0 ; i < nEx ; i + + ) {
let newData = {
x : zList,
2022-12-14 16:59:17 -05:00
y : rBList[i],
mode : "lines",
type : "scatter",
name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
}
data.push(newData);
}
for(let i = 0 ; i < nEx ; i + + ) {
let newData = {
x : zList,
y : rbList[i],
2022-12-14 16:25:37 -05:00
mode : "lines",
2022-12-14 16:59:17 -05:00
line : {dash : 'dashdot', width : 1 },
2022-12-14 16:25:37 -05:00
type : "scatter",
name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
}
data.push(newData);
}
let xxx = [parseInt(document.getElementById('zRange1').value), parseInt(document.getElementById('zRange2').value)];
let yyy = [0, parseInt(document.getElementById('rRange').value)];
let layout = {
xaxis: {range: xxx, title: "Z [mm]", mirror : "allticks", linewidth : "1"},
yaxis: {range: yyy, title: "R [mm]" , mirror : "allticks", linewidth : "1"},
title: reactionName,
showlegend : true,
legend: { yanchor:"top", xanchor:"left", x:"0.01",y:"0.99" }
};
Plotly.newPlot( "Plot_RZ", data, layout);
2022-12-13 23:02:47 -05:00
}
2022-12-13 18:24:26 -05:00
document.getElementById('beam').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
2022-12-13 23:02:47 -05:00
GetMass();
CalConstants();
2022-12-13 18:24:26 -05:00
}
}, false
);
2022-12-14 14:31:54 -05:00
document.getElementById('beamEx').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
CalConstants();
}
}, false
);
2022-12-13 18:24:26 -05:00
document.getElementById('target').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
2022-12-13 23:02:47 -05:00
GetMass();
CalConstants();
2022-12-13 18:24:26 -05:00
}
}, false
);
document.getElementById('light').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
2022-12-13 23:02:47 -05:00
GetMass();
CalConstants();
2022-12-13 18:24:26 -05:00
}
}, false
);
document.getElementById('BField').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
2022-12-13 23:02:47 -05:00
document.getElementById('BRange').value = this.value;
2022-12-14 16:25:37 -05:00
CalConstants();
CalculateEZ();
CalculateRZ();
2022-12-13 18:24:26 -05:00
}
}, false
);
2022-12-14 14:31:54 -05:00
document.getElementById('BRange').oninput = function(){
document.getElementById('BField').value = this.value;
CalConstants();
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
2022-12-14 14:31:54 -05:00
}
2022-12-13 18:24:26 -05:00
document.getElementById('KEA').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
2022-12-13 23:02:47 -05:00
document.getElementById('KEARange').value = this.value;
2022-12-14 16:25:37 -05:00
CalConstants();
CalculateEZ();
CalculateRZ();
2022-12-13 23:02:47 -05:00
}
}, false
);
2022-12-14 14:31:54 -05:00
document.getElementById('KEARange').oninput = function(){
document.getElementById('KEA').value = this.value;
CalConstants();
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
2022-12-14 14:31:54 -05:00
}
2022-12-13 23:02:47 -05:00
document.getElementById('thetaCM').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('thetaCMRange').value = this.value;
2022-12-14 16:25:37 -05:00
CalConstants();
CalculateEZ();
CalculateRZ()
2022-12-13 23:02:47 -05:00
}
}, false
);
2022-12-14 14:31:54 -05:00
document.getElementById('thetaCMRange').oninput = function(){
document.getElementById('thetaCM').value = this.value;
CalConstants();
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ()
2022-12-14 14:31:54 -05:00
}
2022-12-14 16:25:37 -05:00
function AdjustRangeEZ(value){
2022-12-14 14:31:54 -05:00
let haha = parseFloat(value);
let xStart = (haha < 0 ? haha - arrayLen - 100 : haha - 100 ) ;
let xEnd = (haha < 0 ? haha + 100: haha + arrayLen + 100 ) ;
xRange = [xStart, xEnd];
let pos = [];
for( let i = 0; i < nDet ; i + + ) {
let kaka = [];
if( haha < 0 ) {
kaka.push(haha - (i+1) * detLen - i * detGap);
kaka.push(haha - (i) * detLen - i * detGap);
}else{
kaka.push(haha + (i+1) * detLen + i * detGap);
kaka.push(haha + (i) * detLen + i * detGap);
}
pos.push(kaka);
}
//document.getElementById('n0').innerHTML = pos;
//document.getElementById('n2').innerHTML = arrayLen;
//document.getElementById('n1').innerHTML = xRange;
2022-12-14 16:25:37 -05:00
2022-12-14 14:31:54 -05:00
let shapeArray = [];
for( let i = 0; i < nDet ; i + + ) {
let newBlock = {
type: 'rect',
xref: 'x',
yref: 'paper',
x0 : pos[i][0],
x1 : pos[i][1],
y0 : 0,
y1 : 1,
fillcolor : '#9999FF',
opacity : 0.1,
line : { width : 0} }
shapeArray.push(newBlock);
}
var update = {
'xaxis.range' : xRange,
'yaxis.range' : yRange,
'shapes' : shapeArray
}
Plotly.relayout("Plot_EZ", update);
}
2022-12-13 23:02:47 -05:00
document.getElementById('posArray').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('posArrayRange').value = this.value;
2022-12-14 16:25:37 -05:00
AdjustRangeEZ(this.value);
AdjustRecoilPos();
2022-12-13 18:24:26 -05:00
}
}, false
);
2022-12-14 14:31:54 -05:00
document.getElementById('posArrayRange').oninput = function(){
document.getElementById('posArray').value = this.value;
2022-12-14 16:25:37 -05:00
AdjustRangeEZ(this.value);
AdjustRecoilPos();
2022-12-14 14:31:54 -05:00
}
2022-12-14 16:25:37 -05:00
function AdjustRecoilPos(){
let pos1 = parseInt(document.getElementById('posRecoil').value);
let pos2 = pos1+10;
//document.getElementById('n0').innerHTML = pos1;
let shapeRecoil = {
type: 'rect',
xref: 'x',
yref: 'y',
x0 : pos1,
x1 : pos2,
y0 : 10,
y1 : 45,
fillcolor : '#FF0000',
opacity : 0.4,
line : { width : 0}
};
let haha = parseFloat(document.getElementById('posArray').value);
let pos = [];
for( let i = 0; i < nDet ; i + + ) {
let kaka = [];
if( haha < 0 ) {
kaka.push(haha - (i+1) * detLen - i * detGap);
kaka.push(haha - (i) * detLen - i * detGap);
}else{
kaka.push(haha + (i+1) * detLen + i * detGap);
kaka.push(haha + (i) * detLen + i * detGap);
}
pos.push(kaka);
}
let totalShape = [];
for( let i = 0; i < nDet ; i + + ) {
let newhaha = {
type: 'rect',
xref: 'x',
yref: 'y',
x0 : pos[i][0],
x1 : pos[i][1],
y0 : 0,
y1 : perpDistant,
fillcolor : '#9999FF',
opacity : 0.1,
line : { width : 0} };
totalShape.push(newhaha);
}
totalShape.push(shapeRecoil);
var update = {
'shapes' : totalShape
}
Plotly.relayout("Plot_RZ", update);
}
document.getElementById('posRecoil').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('posRecoilRange').value = this.value;
AdjustRecoilPos();
}
}, false
);
document.getElementById('posRecoilRange').oninput = function(){
document.getElementById('posRecoil').value = this.value;
AdjustRecoilPos();
}
function AdjustRangeRZ(){
let zmin = parseInt(document.getElementById('zRange1').value);
let zmax = parseInt(document.getElementById('zRange2').value);
let rmax = parseInt(document.getElementById('rRange').value);
let update = {
'xaxis.range' : [zmin, zmax],
'yaxis.range' : [0, rmax]
}
Plotly.relayout("Plot_RZ", update);
}
document.getElementById('zRange1').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('zRange1Slider').value = this.value;
AdjustRangeRZ();
}
}, false
);
document.getElementById('zRange1Slider').oninput = function(){
document.getElementById('zRange1').value = this.value;
AdjustRangeRZ();
}
document.getElementById('zRange2').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
document.getElementById('zRange2Slider').value = this.value;
AdjustRangeRZ();
}
}, false
);
document.getElementById('zRange2Slider').oninput = function(){
document.getElementById('zRange2').value = this.value;
AdjustRangeRZ();
}
document.getElementById('rRange').addEventListener('keypress',
function(e){
if(e.keyCode == 13){
let rrrr = parseInt(this.value);
if ( rrrr < 1 ) {
document.getElementById('rRange').value = 1;
rrrr = 1;
}
document.getElementById('rRangeSlider').value = rrrr;
AdjustRangeRZ();
}
}, false
);
document.getElementById('rRangeSlider').oninput = function(){
document.getElementById('rRange').value = this.value;
AdjustRangeRZ();
}
2022-12-13 23:02:47 -05:00
window.logMeThis = function(){
SetSSType();
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
2022-12-13 18:24:26 -05:00
}
var FuncEx = window.logMeThis.bind(null, "Ex");
window.addEventListener('keypress', FuncEx);
var FuncThetaCM = window.logMeThis.bind(null, "thetaCM");
window.addEventListener('keypress', FuncThetaCM);
2022-12-13 23:02:47 -05:00
window.checkSSType = function(){
SetSSType();
2022-12-14 16:25:37 -05:00
CalculateEZ();
2022-12-13 23:02:47 -05:00
}
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
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);
2022-12-13 18:24:26 -05:00
2022-12-13 23:02:47 -05:00
GetMass();
CalConstants();
SetSSType();
2022-12-14 16:25:37 -05:00
CalculateEZ();
CalculateRZ();
AdjustRecoilPos();
2022-12-13 18:24:26 -05:00
< / script >
< html >
