diff --git a/heliosmatics.html b/heliosmatics.html
index 772039a..f0772c6 100644
--- a/heliosmatics.html
+++ b/heliosmatics.html
@@ -413,1102 +413,8 @@ Experimental 3D model
-
-function addRow() {
- let table = document.getElementById("ExTable");
- let nRow = table.rows.length;
- let row = table.insertRow(nRow-1);
-
- let energy = Math.random()*Math.min(heavy[4], heavy[5], heavy[6]);
- let angle = Math.floor(Math.random() * 30) + 10;
-
- row.innerHTML = ' | \
- | \
- | \
- | \
- | \
- | \
- | \
- | \
- | \
- | \
- | \
- | ';
-
- CalculateEZ();
- CalculateRZ();
-
- //let table2 = document.getElementById("DWBATable");
- //row = table2.insertRow(nRow-1);
- //row.innerHTML = '' + (nRow-2) + ' | \
- // | \
- // | ';
-}
-
-function deleteRow(){
- let table = document.getElementById("ExTable");
- let nRow = table.rows.length;
- if ( nRow > 3){
- table.deleteRow(nRow-2);
- }
- CalculateEZ();
- CalculateRZ();
-
- //let table2 = document.getElementById("DWBATable");
- //if ( nRow > 3){
- // table2.deleteRow(nRow-2);
- //}
-}
-
-var beam = []; //A, Z, Mass, Name, Sn, Sp, Sa
-var beamMass;
-var target= [];
-var light=[];
-var heavy=[]; //A, Z, Mass, Name, Sn, Sp, Sa
-
-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
-var arrayLen = 50 * 10 + 2 * 9 // SOLARIS
-var detLen = 50; // SOLARIS
-var nDet = 10; // SOLARIS
-var detGap = 2; // SOLARIS
-var arrayPos = [];
-
-const c = 299.792468; // mm/ns
-
-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 xList =[]; // for E-Z plot
-var yList =[]; // for E-Z plot
-var ExList = [];
-var Zb0List = [];
-var rhoList = [];
-
-var ZB0List = [];
-var rhoBList = [];
-
-var xRange ;
-var yRange = [0, 12];
-
-var zList = []; // for R-Z plot
-var rbList = []; // for R-Z plot
-var rBList = []; // for R-Z plot
-
-var color = ['rgb(31,119,180)', // muted blue
- 'rgb(255,127,14)', // safety orange
- 'rgb(44,160,44)', // cooked asparagus green
- 'rgb(214,39,40)', // brick red
- 'rgb(148,103,189)', // muted purple
- 'rgb(140,86,75)', // chestnut brown
- 'rgb(277,119,194)', // raspberry yogurt pink
- 'rgb(127,127,127)', // middle gray
- 'rgb(118,189,34)', // curry yellow-green
- 'rgb(23,190,207)']; // blue-teal
-
-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;
- beam[4] = parseFloat(mass[4]);
- beam[5] = parseFloat(mass[5]);
- beam[6] = parseFloat(mass[6]);
- }
- 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]?.trim();
- heavy[4] = parseFloat(mass[4]);
- heavy[5] = parseFloat(mass[5]);
- heavy[6] = parseFloat(mass[6]);
- }
- client.open('GET', str, false);
- client.send();
-
-}
-
-function GetMass(){
-
- GetMassFromSym(document.getElementById('beam').value, 0);
- GetMassFromSym(document.getElementById('target').value, 1);
- GetMassFromSym(document.getElementById('light').value, 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('heavySp').innerHTML = 'Sn: ' + heavy[4] + ' MeV, Sp: ' + heavy[5] + ' MeV, Sa : ' + heavy[6] + ' MeV';
-
- //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 = "min Beam Energy: " + minKEA.toFixed(3) + " MeV/u";
- 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];
-
-}
-
-function CalArrayPos(){
- let haha = parseFloat(document.getElementById('posArray').value);
- arrayPos = [];
- 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) * detLen + i * detGap);
- kaka.push(haha + (i+1) * detLen + i * detGap);
- }
- arrayPos.push(kaka);
- }
-}
-
-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;
- 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);
-
- CalArrayPos();
- //document.getElementById('n1').innerHTML = perpDistant;
-}
-
-function CalculateEZ(){
-
- let tableEx = document.getElementById("ExTable");
- let nRow = tableEx.rows.length;
-
- xList = [];
- yList = [];
- ExList = [];
- Zb0List = [];
- rhoList = [];
- ZB0List = [];
- rhoBList = [];
-
- //alert("CalculateEZ called, nRow = " + nRow);
-
- 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
- Zb0List.push(Zb0);
-
- let rho = k_cm * ss/c/light[1]/Math.abs(BField) * 1000; // mm
- rhoList.push(rho);
- 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
- ZB0List.push(ZB0);
- rhoBList.push(rhoB);
-
- 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 = Math.PI * ez_slope / beta / k_cm * perpDistant / Math.sin(thetaCM); // perpDistant / 2/ rho(thetaCM)
- if( !isFinite(temp) ) continue;
- let pxTemp = 1. /alpha * (gamma * beta * qb - gamma * k_cm * Math.cos(thetaCM)) * (1 - Math.asin(temp)/Math.PI) ;
- 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);
-
- };
-
- PlotEZ();
- AdjustRangeEZ();
-
-}
-
-function PlotEZ(){
-
- SetSSType();
- Plotly.purge("Plot_EZ");
-
- let nEx = xList.length;
- let data = [];
-
- for( let i = 0; i < nEx; i++){
- let kaka = color[i%10];
- let newData = {
- x : xList[i],
- y : yList[i],
- mode:"lines",
- type:"scatter",
- name:"Ex="+ExList[i],
- marker : { color : kaka}
- }
- data.push(newData);
- }
-
- 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;
-
- let haha2 = parseFloat(document.getElementById('eRangeSlider').value);
- yRange = [0, haha2];
-
- let layout = {
- xaxis: {range: xRange, title: { text : "Z [mm]", standoff : 1}, mirror : "allticks", linewidth : "1"},
- yaxis: {range: yRange, title: "Energy [MeV]" , mirror : "allticks", linewidth : "1"},
- title: reactionName,
- dragmode : "pan",
- margin: { l: 40, r: 40, b : 60, t : 40},
- legend: {yanchor:"top", xanchor:"left", x:"0.01",y:"0.99" }
- };
-
-
- Plotly.newPlot( "Plot_EZ", data, layout, {responsive: true});
-
- PlotThetaCMLine(document.getElementById('thetaCM').value);
- PlotBore();
-
- PlotRZ();
-}
-
-function PlotThetaCMLine(thetaCM){
- let cs = Math.cos(thetaCM * Math.PI /180);
- let ss = Math.sin(thetaCM * Math.PI /180);
-
- let zzList = [];
- let eList = [];
-
- for( let z = -2000; z < 2000; z+=30){
- zzList.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 : zzList,
- y: eList,
- name:"thetaCM=" + thetaCM ,
- marker : { color : 'rgb(100,100,100)'},
- line : {dash : 'solid', width : 1 }
- }
- , 0);
-}
-
-function PlotBore(){
- let zzList = [];
- let eList = [];
-
- for( let z = -2000; z < 2000; z+=30){
- zzList.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 : zzList, y: eList, name:"Bore", marker : { color : 'rgb(200,200, 200)'} }, 0);
-}
-
-function CalculateRZ(){
- // this rquire ZB0List and rhoBList from CalculateEZ();
-
- zList = [];
- rbList = [];
- rBList = [];
-
- 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])));
- }
- }
- 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);
- }
-
- PlotRZ();
- AdjustRecoilPos();
-}
-
-function PlotRZ(){
- Plotly.purge("Plot_RZ");
-
- let data = [];
- let nEx = ExList.length;
-
- for(let i = 0 ; i < nEx; i++ ){
- let kaka = color[i%10];
- let newData = {
- x : zList,
- y : rBList[i],
- mode : "lines",
- type : "scatter",
- name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
- marker : { color : kaka}
- }
- data.push(newData);
- }
- for(let i = 0 ; i < nEx; i++ ){
- let kaka = color[i%10];
- let newData = {
- x : zList,
- y : rbList[i],
- mode : "lines",
- line : {dash : 'dashdot', width : 1 },
- type : "scatter",
- name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
- marker : { color : kaka}
- }
- 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: { text : "Z [mm]", standoff : 1}, mirror : "allticks", linewidth : "1"},
- yaxis: {range: yyy, title: "R [mm]" , mirror : "allticks", linewidth : "1"},
- title: reactionName,
- showlegend : true,
- dragmode : "pan",
- margin: { l: 40, r: 40, b : 60, t : 40},
- legend: { yanchor:"top", xanchor:"left", x:"0.01",y:"0.99"}
- };
-
- Plotly.newPlot( "Plot_RZ", data, layout, {responsive: true});
-
-}
-
-function CalZ(theta, k_cm, qb){
- let cs = Math.cos(theta*Math.PI/180.);
- let ss = Math.sin(theta*Math.PI/180.);
-
- let Zb0 = (gamma*beta* qb - gamma * k_cm * cs )/alpha; //mm
- let rho = k_cm * ss / c/ light[1]/ Math.abs(BField) * 1000; // mm
-
- return Zb0 * ( 1 - Math.sin( perpDistant/2/rho ) / Math.PI );
-}
-
-function SearchThetaCM(Z, thetaMin, k_cm, qb){
-
- let step = 0.01;
-
- for( let yy = thetaMin; yy < 180; yy += step){
- if( yy > 12 ) step = 0.1;
- if( yy > 20 ) step = 0.5;
- if( yy > 40 ) step = 1.0;
- let z1 = CalZ(yy, k_cm, qb);
- let z2 = CalZ(yy+step, k_cm, qb);
- //console.log(Z, ", ", yy.toFixed(3), ", ", z1.toFixed(3), ", ", z2.toFixed(3), ", ", (z2-z1).toFixed(3) );
- if( z1 > z2 ){ /// in the bending range
- continue;
- }else{
- if( Z < z1 ){
- // return yy;
- break;
- }
- if( z1<= Z && Z <= z2 ){ // do a linear approximation
- return (Z-z1)/(z2-z1)*step + yy;
- }
- }
- }
- return NaN;
-}
-
-function CalThetaCM(){
- let ex = parseFloat(document.getElementById('Ex0').value);
- let angGate = parseFloat(document.getElementById('thetaCMGate').value);
- let xGate = parseFloat(document.getElementById('XGate').value);
-
- CalArrayPos();
-
- 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 qb = Math.sqrt(light[2]*light[2]+k_cm*k_cm);
-
- let table = document.getElementById('thetaCMTable');
- let nRow = table.rows.length;
- if( nRow > 1 ){
- for( let i = nRow-1; i > 0; i -- ){
- table.deleteRow(i);
- }
- }
-
- for( let i = 0; i < arrayPos.length; i++){
- let row = table.insertRow(i+1);
- row.insertCell().innerHTML = i;
-
- let p1 = (arrayPos[i][0] + detLen * (100. - xGate)/200.);
- let p2 = (arrayPos[i][1] - detLen * (100. - xGate)/200.);
-
- row.insertCell().innerHTML = arrayPos[i][0] + "(" + p1 + ")";
- row.insertCell().innerHTML = arrayPos[i][1] + "(" + p2 + ")";
-
- ///Search thetaCM for Z
- let a1 = SearchThetaCM(p1, angGate, k_cm, qb);
- let a2 = SearchThetaCM(p2, angGate, k_cm, qb);
- row.insertCell().innerHTML = a1.toFixed(2);
- row.insertCell().innerHTML = a2.toFixed(2);
-
- let am = (a2+a1)/2;
- let da = Math.abs(a2-a1);
-
- row.insertCell().innerHTML = am.toFixed(2);
- row.insertCell().innerHTML = da.toFixed(2);
- row.insertCell().innerHTML = (Math.sin(am*Math.PI/180) * da * Math.PI/180).toExponential(2);
-
- }
-
-}
-
-document.getElementById('beam').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- GetMass();
- CalConstants();
- }
- }, false
-);
-
-document.getElementById('beamEx').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- 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){
- //document.getElementById('n0').innerHTML = e.keyCode;
- if(e.keyCode == 13){
- document.getElementById('BRange').value = this.value;
- CalConstants();
- CalculateEZ();
- CalculateRZ();
- CalThetaCM();
- }
- }, false
-);
-
-document.getElementById('BRange').oninput = function(){
- document.getElementById('BField').value = this.value;
- CalConstants();
- CalculateEZ();
- CalculateRZ();
- AdjustRangeEZ();
- CalThetaCM();
-}
-
-document.getElementById('KEA').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- document.getElementById('KEARange').value = this.value;
- CalConstants();
- CalculateEZ();
- CalculateRZ();
- CalThetaCM();
- }
- }, false
-);
-
-document.getElementById('KEARange').oninput = function(){
- document.getElementById('KEA').value = this.value;
- CalConstants();
- CalculateEZ();
- CalculateRZ();
- CalThetaCM();
-}
-
-document.getElementById('thetaCM').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- document.getElementById('thetaCMRange').value = this.value;
- CalConstants();
- CalculateEZ();
- }
- }, false
-);
-
-document.getElementById('thetaCMRange').oninput = function(){
- document.getElementById('thetaCM').value = this.value;
- CalConstants();
- CalculateEZ();
-}
-
-function AdjustRangeEZ(){
- let haha = parseFloat(document.getElementById('posArray').value);
- let xStart = (haha < 0 ? haha - arrayLen - 100 : haha - 100);
- let xEnd = (haha < 0 ? haha + 100: haha + arrayLen + 100);
-
- CalArrayPos();
-
- xRange = [xStart, xEnd];
- yRange = [0, parseFloat(document.getElementById('eRangeSlider').value)];
-
- //document.getElementById('n0').innerHTML = pos;
- //document.getElementById('n2').innerHTML = arrayLen;
- //document.getElementById('n1').innerHTML = xRange;
-
- let shapeArray = [];
- for( let i = 0; i < nDet; i++){
- let newBlock = {
- type: 'rect',
- xref: 'x',
- yref: 'paper',
- x0 : arrayPos[i][0],
- x1 : arrayPos[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);
-}
-
-document.getElementById('posArray').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- document.getElementById('posArrayRange').value = this.value;
- AdjustRangeEZ();
- AdjustRecoilPos();
- CalThetaCM();
- }
- }, false
-);
-
-document.getElementById('posArrayRange').oninput = function(){
- document.getElementById('posArray').value = this.value;
- AdjustRangeEZ();
- AdjustRecoilPos();
- CalThetaCM();
-}
-
-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 totalShape = [];
- for( let i = 0; i < nDet; i++){
- let newhaha = {
- type: 'rect',
- xref: 'x',
- yref: 'y',
- x0 : arrayPos[i][0],
- x1 : arrayPos[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();
-}
-
-document.getElementById('eRange').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- let rrrr = parseInt(this.value);
- if ( rrrr < 1 ){
- document.getElementById('eRange').value = 1;
- rrrr = 1;
- }
- document.getElementById('eRangeSlider').value = rrrr;
- AdjustRangeEZ();
- }
- }, false
-);
-
-document.getElementById('eRangeSlider').oninput = function(){
- document.getElementById('eRange').value = this.value;
- AdjustRangeEZ();
-}
-
-window.logMeThis = function(){
- SetSSType();
- CalculateEZ();
- CalculateRZ();
-}
-
-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();
- CalculateEZ();
- CalThetaCM();
-}
-
-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('Ex0').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- CalThetaCM();
- }
- }, false
-);
-
-document.getElementById('thetaCMGate').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- CalThetaCM();
- }
- }, false
-);
-
-document.getElementById('XGate').addEventListener('keypress',
- function(e){
- if(e.keyCode == 13){
- CalThetaCM();
- }
- }, false
-);
-
-
-function HalfCylinder(r, y, d, up){
- let a = [];
- let b = [];
- let c = [];
-
- let nStep = 60;
- let step = 2*r/nStep;
-
- for( let i = 0; i <= nStep ; i++){
- let x = -r + i * step;
- a.push(x);
- c.push(Math.sqrt(r*r - x*x) * up);
- b.push(y);
- }
-
-
- for( let i = 0; i <= nStep; i ++){
- let x = r - i * step;
- a.push(x);
- c.push(Math.sqrt(r*r - x*x) * up);
- b.push(y+d);
- }
-
-
- let haha = [];
- haha.push(a);
- haha.push(b);
- haha.push(c);
-
- return haha;
-}
-
-function detMesh(startPos, phi){
- let aaa = 110; // prepdist
- let len = 500; // detLen
- let www = 100; //det width
-
- let a =[];
- let b =[];
- let c =[];
-
- let cs = Math.cos(phi * Math.PI / 180);
- let ss = Math.sin(phi * Math.PI / 180);
-
- let A0 = aaa * cs - www / 2 * ss;
- let A1 = aaa * ss + www / 2 * cs;
- let B0 = aaa * cs + www / 2 * ss;
- let B1 = aaa * ss - www / 2 * cs;
-
- a.push(A0); b.push(A1); c.push(startPos);
- a.push(B0); b.push(B1); c.push(startPos);
- a.push(B0); b.push(B1); c.push(startPos + len);
- a.push(A0); b.push(A1); c.push(startPos + len);
-
- let haha = [];
- haha.push(a);
- haha.push(b);
- haha.push(c);
-
- return haha;
-
-}
-
-function Helix(theta, phi, rho, sign, nCyc){
- // sign = B-field
- // nCyc < 0 = updatream
-
- let a = [];
- let b = [];
- let c = [];
-
- let deg = Math.PI/180;
- let nStep = 100;
- let ts = Math.tan(theta * deg);
- let zRange = nCyc * rho/ts * Math.PI * 2;
- let zStep = zRange/nStep;
-
-
- for( let i = 0; i < nStep; i++){
- let zzz;
- zzz = i * zStep;
- b.push(zzz);
- a.push( rho * ( Math.sin( ts * zzz / rho - phi * deg ) + Math.sin(phi * deg) ) );
- c.push( rho * sign * (Math.cos( ts * zzz / rho - phi * deg ) - Math.cos(phi * deg)) );
- }
-
- let haha = [];
- haha.push(a);
- haha.push(b);
- haha.push(c);
-
- return haha;
-
-}
-
-
-function Plot3D(){
-
- let u1 = HalfCylinder(bore, -1000, 2500., 1);
- let d1 = HalfCylinder(bore, -1000, 2500., -1);
-
- let line1 = Helix(40, 0, 80, -1, -1);
- let line2 = Helix(2, 180, 10, -1, 0.5);
-
- let data = [
- {type : 'mesh3d', x : u1[0], y : u1[1], z : u1[2], hoverinfo: 'none', opacity : 0.2, color : "#DDDDDD"}
- ,{type : 'mesh3d', x : d1[0], y : d1[1], z : d1[2], hoverinfo: 'none', opacity : 0.2, color : "#DDDDDD"}
- ,{type : 'scatter3d', mode : 'lines', x : line1[0], y : line1[1], z : line1[2], name : "haha"}
- ,{type : 'scatter3d', mode : 'lines', x : line2[0], y : line2[1], z : line2[2], name : "kaka"}
- ];
-
-
- /*
- let data = [];
-
- for( let i = 0; i < 1; i++){
- let haha = detMesh(-200, 360/6*i + 90);
-
- console.log(haha[0], "| ", haha[1], "| ", haha[2]);
- data.push(
- {type : 'mesh3d', x : haha[0], y : haha[1], z : haha[2], hoverinfo: 'none', opacity : 1.0, color : "#FF0000"}
- );
- }*/
-
- let layout = {
- aspectmode: 'manual',
- aspectratio:{ x : 0.5, y : 0.5, z : 2.5 },
- margin: { l: 40, r: 40, b : 60, t : 40},
- };
-
- Plotly.newPlot('Plot_3D', data, layout);
-
-}
-
-
-GetMass();
-CalConstants();
-SetSSType();
-CalculateEZ();
-CalculateRZ();
-AdjustRecoilPos();
-CalThetaCM();
-
-//Plot3D();
-
-
diff --git a/heliosmatics.js b/heliosmatics.js
new file mode 100644
index 0000000..f9e052b
--- /dev/null
+++ b/heliosmatics.js
@@ -0,0 +1,1095 @@
+
+function addRow() {
+ let table = document.getElementById("ExTable");
+ let nRow = table.rows.length;
+ let row = table.insertRow(nRow-1);
+
+ let energy = Math.random()*Math.min(heavy[4], heavy[5], heavy[6]);
+ let angle = Math.floor(Math.random() * 30) + 10;
+
+ row.innerHTML = ' | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | \
+ | ';
+
+ CalculateEZ();
+ CalculateRZ();
+
+ //let table2 = document.getElementById("DWBATable");
+ //row = table2.insertRow(nRow-1);
+ //row.innerHTML = '' + (nRow-2) + ' | \
+ // | \
+ // | ';
+}
+
+function deleteRow(){
+ let table = document.getElementById("ExTable");
+ let nRow = table.rows.length;
+ if ( nRow > 3){
+ table.deleteRow(nRow-2);
+ }
+ CalculateEZ();
+ CalculateRZ();
+
+ //let table2 = document.getElementById("DWBATable");
+ //if ( nRow > 3){
+ // table2.deleteRow(nRow-2);
+ //}
+}
+
+var beam = []; //A, Z, Mass, Name, Sn, Sp, Sa
+var beamMass;
+var target= [];
+var light=[];
+var heavy=[]; //A, Z, Mass, Name, Sn, Sp, Sa
+
+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
+var arrayLen = 50 * 10 + 2 * 9 // SOLARIS
+var detLen = 50; // SOLARIS
+var nDet = 10; // SOLARIS
+var detGap = 2; // SOLARIS
+var arrayPos = [];
+
+const c = 299.792468; // mm/ns
+
+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 xList =[]; // for E-Z plot
+var yList =[]; // for E-Z plot
+var ExList = [];
+var Zb0List = [];
+var rhoList = [];
+
+var ZB0List = [];
+var rhoBList = [];
+
+var xRange ;
+var yRange = [0, 12];
+
+var zList = []; // for R-Z plot
+var rbList = []; // for R-Z plot
+var rBList = []; // for R-Z plot
+
+var color = ['rgb(31,119,180)', // muted blue
+ 'rgb(255,127,14)', // safety orange
+ 'rgb(44,160,44)', // cooked asparagus green
+ 'rgb(214,39,40)', // brick red
+ 'rgb(148,103,189)', // muted purple
+ 'rgb(140,86,75)', // chestnut brown
+ 'rgb(277,119,194)', // raspberry yogurt pink
+ 'rgb(127,127,127)', // middle gray
+ 'rgb(118,189,34)', // curry yellow-green
+ 'rgb(23,190,207)']; // blue-teal
+
+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;
+ beam[4] = parseFloat(mass[4]);
+ beam[5] = parseFloat(mass[5]);
+ beam[6] = parseFloat(mass[6]);
+ }
+ 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]?.trim();
+ heavy[4] = parseFloat(mass[4]);
+ heavy[5] = parseFloat(mass[5]);
+ heavy[6] = parseFloat(mass[6]);
+ }
+ client.open('GET', str, false);
+ client.send();
+
+}
+
+function GetMass(){
+
+ GetMassFromSym(document.getElementById('beam').value, 0);
+ GetMassFromSym(document.getElementById('target').value, 1);
+ GetMassFromSym(document.getElementById('light').value, 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('heavySp').innerHTML = 'Sn: ' + heavy[4] + ' MeV, Sp: ' + heavy[5] + ' MeV, Sa : ' + heavy[6] + ' MeV';
+
+ //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 = "min Beam Energy: " + minKEA.toFixed(3) + " MeV/u";
+ 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];
+
+}
+
+function CalArrayPos(){
+ let haha = parseFloat(document.getElementById('posArray').value);
+ arrayPos = [];
+ 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) * detLen + i * detGap);
+ kaka.push(haha + (i+1) * detLen + i * detGap);
+ }
+ arrayPos.push(kaka);
+ }
+}
+
+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;
+ 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);
+
+ CalArrayPos();
+ //document.getElementById('n1').innerHTML = perpDistant;
+}
+
+function CalculateEZ(){
+
+ let tableEx = document.getElementById("ExTable");
+ let nRow = tableEx.rows.length;
+
+ xList = [];
+ yList = [];
+ ExList = [];
+ Zb0List = [];
+ rhoList = [];
+ ZB0List = [];
+ rhoBList = [];
+
+ //alert("CalculateEZ called, nRow = " + nRow);
+
+ 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
+ Zb0List.push(Zb0);
+
+ let rho = k_cm * ss/c/light[1]/Math.abs(BField) * 1000; // mm
+ rhoList.push(rho);
+ 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
+ ZB0List.push(ZB0);
+ rhoBList.push(rhoB);
+
+ 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 = Math.PI * ez_slope / beta / k_cm * perpDistant / Math.sin(thetaCM); // perpDistant / 2/ rho(thetaCM)
+ if( !isFinite(temp) ) continue;
+ let pxTemp = 1. /alpha * (gamma * beta * qb - gamma * k_cm * Math.cos(thetaCM)) * (1 - Math.asin(temp)/Math.PI) ;
+ 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);
+
+ };
+
+ PlotEZ();
+ AdjustRangeEZ();
+
+}
+
+function PlotEZ(){
+
+ SetSSType();
+ Plotly.purge("Plot_EZ");
+
+ let nEx = xList.length;
+ let data = [];
+
+ for( let i = 0; i < nEx; i++){
+ let kaka = color[i%10];
+ let newData = {
+ x : xList[i],
+ y : yList[i],
+ mode:"lines",
+ type:"scatter",
+ name:"Ex="+ExList[i],
+ marker : { color : kaka}
+ }
+ data.push(newData);
+ }
+
+ 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;
+
+ let haha2 = parseFloat(document.getElementById('eRangeSlider').value);
+ yRange = [0, haha2];
+
+ let layout = {
+ xaxis: {range: xRange, title: { text : "Z [mm]", standoff : 1}, mirror : "allticks", linewidth : "1"},
+ yaxis: {range: yRange, title: "Energy [MeV]" , mirror : "allticks", linewidth : "1"},
+ title: reactionName,
+ dragmode : "pan",
+ margin: { l: 40, r: 40, b : 60, t : 40},
+ legend: {yanchor:"top", xanchor:"left", x:"0.01",y:"0.99" }
+ };
+
+
+ Plotly.newPlot( "Plot_EZ", data, layout, {responsive: true});
+
+ PlotThetaCMLine(document.getElementById('thetaCM').value);
+ PlotBore();
+
+ PlotRZ();
+}
+
+function PlotThetaCMLine(thetaCM){
+ let cs = Math.cos(thetaCM * Math.PI /180);
+ let ss = Math.sin(thetaCM * Math.PI /180);
+
+ let zzList = [];
+ let eList = [];
+
+ for( let z = -2000; z < 2000; z+=30){
+ zzList.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 : zzList,
+ y: eList,
+ name:"thetaCM=" + thetaCM ,
+ marker : { color : 'rgb(100,100,100)'},
+ line : {dash : 'solid', width : 1 }
+ }
+ , 0);
+}
+
+function PlotBore(){
+ let zzList = [];
+ let eList = [];
+
+ for( let z = -2000; z < 2000; z+=30){
+ zzList.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 : zzList, y: eList, name:"Bore", marker : { color : 'rgb(200,200, 200)'} }, 0);
+}
+
+function CalculateRZ(){
+ // this rquire ZB0List and rhoBList from CalculateEZ();
+
+ zList = [];
+ rbList = [];
+ rBList = [];
+
+ 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])));
+ }
+ }
+ 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);
+ }
+
+ PlotRZ();
+ AdjustRecoilPos();
+}
+
+function PlotRZ(){
+ Plotly.purge("Plot_RZ");
+
+ let data = [];
+ let nEx = ExList.length;
+
+ for(let i = 0 ; i < nEx; i++ ){
+ let kaka = color[i%10];
+ let newData = {
+ x : zList,
+ y : rBList[i],
+ mode : "lines",
+ type : "scatter",
+ name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
+ marker : { color : kaka}
+ }
+ data.push(newData);
+ }
+ for(let i = 0 ; i < nEx; i++ ){
+ let kaka = color[i%10];
+ let newData = {
+ x : zList,
+ y : rbList[i],
+ mode : "lines",
+ line : {dash : 'dashdot', width : 1 },
+ type : "scatter",
+ name : "Ex="+ExList[i] + ",theta=" + document.getElementById('theta'+(i+1)).value,
+ marker : { color : kaka}
+ }
+ 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: { text : "Z [mm]", standoff : 1}, mirror : "allticks", linewidth : "1"},
+ yaxis: {range: yyy, title: "R [mm]" , mirror : "allticks", linewidth : "1"},
+ title: reactionName,
+ showlegend : true,
+ dragmode : "pan",
+ margin: { l: 40, r: 40, b : 60, t : 40},
+ legend: { yanchor:"top", xanchor:"left", x:"0.01",y:"0.99"}
+ };
+
+ Plotly.newPlot( "Plot_RZ", data, layout, {responsive: true});
+
+}
+
+function CalZ(theta, k_cm, qb){
+ let cs = Math.cos(theta*Math.PI/180.);
+ let ss = Math.sin(theta*Math.PI/180.);
+
+ let Zb0 = (gamma*beta* qb - gamma * k_cm * cs )/alpha; //mm
+ let rho = k_cm * ss / c/ light[1]/ Math.abs(BField) * 1000; // mm
+
+ return Zb0 * ( 1 - Math.sin( perpDistant/2/rho ) / Math.PI );
+}
+
+function SearchThetaCM(Z, thetaMin, k_cm, qb){
+
+ let step = 0.01;
+
+ for( let yy = thetaMin; yy < 180; yy += step){
+ if( yy > 12 ) step = 0.1;
+ if( yy > 20 ) step = 0.5;
+ if( yy > 40 ) step = 1.0;
+ let z1 = CalZ(yy, k_cm, qb);
+ let z2 = CalZ(yy+step, k_cm, qb);
+ //console.log(Z, ", ", yy.toFixed(3), ", ", z1.toFixed(3), ", ", z2.toFixed(3), ", ", (z2-z1).toFixed(3) );
+ if( z1 > z2 ){ /// in the bending range
+ continue;
+ }else{
+ if( Z < z1 ){
+ // return yy;
+ break;
+ }
+ if( z1<= Z && Z <= z2 ){ // do a linear approximation
+ return (Z-z1)/(z2-z1)*step + yy;
+ }
+ }
+ }
+ return NaN;
+}
+
+function CalThetaCM(){
+ let ex = parseFloat(document.getElementById('Ex0').value);
+ let angGate = parseFloat(document.getElementById('thetaCMGate').value);
+ let xGate = parseFloat(document.getElementById('XGate').value);
+
+ CalArrayPos();
+
+ 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 qb = Math.sqrt(light[2]*light[2]+k_cm*k_cm);
+
+ let table = document.getElementById('thetaCMTable');
+ let nRow = table.rows.length;
+ if( nRow > 1 ){
+ for( let i = nRow-1; i > 0; i -- ){
+ table.deleteRow(i);
+ }
+ }
+
+ for( let i = 0; i < arrayPos.length; i++){
+ let row = table.insertRow(i+1);
+ row.insertCell().innerHTML = i;
+
+ let p1 = (arrayPos[i][0] + detLen * (100. - xGate)/200.);
+ let p2 = (arrayPos[i][1] - detLen * (100. - xGate)/200.);
+
+ row.insertCell().innerHTML = arrayPos[i][0] + "(" + p1 + ")";
+ row.insertCell().innerHTML = arrayPos[i][1] + "(" + p2 + ")";
+
+ ///Search thetaCM for Z
+ let a1 = SearchThetaCM(p1, angGate, k_cm, qb);
+ let a2 = SearchThetaCM(p2, angGate, k_cm, qb);
+ row.insertCell().innerHTML = a1.toFixed(2);
+ row.insertCell().innerHTML = a2.toFixed(2);
+
+ let am = (a2+a1)/2;
+ let da = Math.abs(a2-a1);
+
+ row.insertCell().innerHTML = am.toFixed(2);
+ row.insertCell().innerHTML = da.toFixed(2);
+ row.insertCell().innerHTML = (Math.sin(am*Math.PI/180) * da * Math.PI/180).toExponential(2);
+
+ }
+
+}
+
+document.getElementById('beam').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ GetMass();
+ CalConstants();
+ }
+ }, false
+);
+
+document.getElementById('beamEx').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ 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){
+ //document.getElementById('n0').innerHTML = e.keyCode;
+ if(e.keyCode == 13){
+ document.getElementById('BRange').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ CalculateRZ();
+ CalThetaCM();
+ }
+ }, false
+);
+
+document.getElementById('BRange').oninput = function(){
+ document.getElementById('BField').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ CalculateRZ();
+ AdjustRangeEZ();
+ CalThetaCM();
+}
+
+document.getElementById('KEA').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ document.getElementById('KEARange').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ CalculateRZ();
+ CalThetaCM();
+ }
+ }, false
+);
+
+document.getElementById('KEARange').oninput = function(){
+ document.getElementById('KEA').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ CalculateRZ();
+ CalThetaCM();
+}
+
+document.getElementById('thetaCM').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ document.getElementById('thetaCMRange').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ AdjustRecoilPos();
+ }
+ }, false
+);
+
+document.getElementById('thetaCMRange').oninput = function(){
+ document.getElementById('thetaCM').value = this.value;
+ CalConstants();
+ CalculateEZ();
+ AdjustRecoilPos();
+}
+
+function AdjustRangeEZ(){
+ let haha = parseFloat(document.getElementById('posArray').value);
+ let xStart = (haha < 0 ? haha - arrayLen - 100 : haha - 100);
+ let xEnd = (haha < 0 ? haha + 100: haha + arrayLen + 100);
+
+ CalArrayPos();
+
+ xRange = [xStart, xEnd];
+ yRange = [0, parseFloat(document.getElementById('eRangeSlider').value)];
+
+ //document.getElementById('n0').innerHTML = pos;
+ //document.getElementById('n2').innerHTML = arrayLen;
+ //document.getElementById('n1').innerHTML = xRange;
+
+ let shapeArray = [];
+ for( let i = 0; i < nDet; i++){
+ let newBlock = {
+ type: 'rect',
+ xref: 'x',
+ yref: 'paper',
+ x0 : arrayPos[i][0],
+ x1 : arrayPos[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);
+}
+
+document.getElementById('posArray').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ document.getElementById('posArrayRange').value = this.value;
+ AdjustRangeEZ();
+ AdjustRecoilPos();
+ CalThetaCM();
+ }
+ }, false
+);
+
+document.getElementById('posArrayRange').oninput = function(){
+ document.getElementById('posArray').value = this.value;
+ AdjustRangeEZ();
+ AdjustRecoilPos();
+ CalThetaCM();
+}
+
+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 totalShape = [];
+ for( let i = 0; i < nDet; i++){
+ let newhaha = {
+ type: 'rect',
+ xref: 'x',
+ yref: 'y',
+ x0 : arrayPos[i][0],
+ x1 : arrayPos[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();
+}
+
+document.getElementById('eRange').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ let rrrr = parseInt(this.value);
+ if ( rrrr < 1 ){
+ document.getElementById('eRange').value = 1;
+ rrrr = 1;
+ }
+ document.getElementById('eRangeSlider').value = rrrr;
+ AdjustRangeEZ();
+ }
+ }, false
+);
+
+document.getElementById('eRangeSlider').oninput = function(){
+ document.getElementById('eRange').value = this.value;
+ AdjustRangeEZ();
+}
+
+window.logMeThis = function(){
+ SetSSType();
+ CalculateEZ();
+ CalculateRZ();
+}
+
+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();
+ CalculateEZ();
+ CalThetaCM();
+}
+
+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('Ex0').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ CalThetaCM();
+ }
+ }, false
+);
+
+document.getElementById('thetaCMGate').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ CalThetaCM();
+ }
+ }, false
+);
+
+document.getElementById('XGate').addEventListener('keypress',
+ function(e){
+ if(e.keyCode == 13){
+ CalThetaCM();
+ }
+ }, false
+);
+
+
+function HalfCylinder(r, y, d, up){
+ let a = [];
+ let b = [];
+ let c = [];
+
+ let nStep = 60;
+ let step = 2*r/nStep;
+
+ for( let i = 0; i <= nStep ; i++){
+ let x = -r + i * step;
+ a.push(x);
+ c.push(Math.sqrt(r*r - x*x) * up);
+ b.push(y);
+ }
+
+
+ for( let i = 0; i <= nStep; i ++){
+ let x = r - i * step;
+ a.push(x);
+ c.push(Math.sqrt(r*r - x*x) * up);
+ b.push(y+d);
+ }
+
+
+ let haha = [];
+ haha.push(a);
+ haha.push(b);
+ haha.push(c);
+
+ return haha;
+}
+
+function detMesh(startPos, phi){
+ let aaa = 110; // prepdist
+ let len = 500; // detLen
+ let www = 100; //det width
+
+ let a =[];
+ let b =[];
+ let c =[];
+
+ let cs = Math.cos(phi * Math.PI / 180);
+ let ss = Math.sin(phi * Math.PI / 180);
+
+ let A0 = aaa * cs - www / 2 * ss;
+ let A1 = aaa * ss + www / 2 * cs;
+ let B0 = aaa * cs + www / 2 * ss;
+ let B1 = aaa * ss - www / 2 * cs;
+
+ a.push(A0); b.push(A1); c.push(startPos);
+ a.push(B0); b.push(B1); c.push(startPos);
+ a.push(B0); b.push(B1); c.push(startPos + len);
+ a.push(A0); b.push(A1); c.push(startPos + len);
+
+ let haha = [];
+ haha.push(a);
+ haha.push(b);
+ haha.push(c);
+
+ return haha;
+
+}
+
+function Helix(theta, phi, rho, sign, nCyc){
+ // sign = B-field
+ // nCyc < 0 = updatream
+
+ let a = [];
+ let b = [];
+ let c = [];
+
+ let deg = Math.PI/180;
+ let nStep = 100;
+ let ts = Math.tan(theta * deg);
+ let zRange = nCyc * rho/ts * Math.PI * 2;
+ let zStep = zRange/nStep;
+
+
+ for( let i = 0; i < nStep; i++){
+ let zzz;
+ zzz = i * zStep;
+ b.push(zzz);
+ a.push( rho * ( Math.sin( ts * zzz / rho - phi * deg ) + Math.sin(phi * deg) ) );
+ c.push( rho * sign * (Math.cos( ts * zzz / rho - phi * deg ) - Math.cos(phi * deg)) );
+ }
+
+ let haha = [];
+ haha.push(a);
+ haha.push(b);
+ haha.push(c);
+
+ return haha;
+
+}
+
+
+function Plot3D(){
+
+ let u1 = HalfCylinder(bore, -1000, 2500., 1);
+ let d1 = HalfCylinder(bore, -1000, 2500., -1);
+
+ let line1 = Helix(40, 0, 80, -1, -1);
+ let line2 = Helix(2, 180, 10, -1, 0.5);
+
+ let data = [
+ {type : 'mesh3d', x : u1[0], y : u1[1], z : u1[2], hoverinfo: 'none', opacity : 0.2, color : "#DDDDDD"}
+ ,{type : 'mesh3d', x : d1[0], y : d1[1], z : d1[2], hoverinfo: 'none', opacity : 0.2, color : "#DDDDDD"}
+ ,{type : 'scatter3d', mode : 'lines', x : line1[0], y : line1[1], z : line1[2], name : "haha"}
+ ,{type : 'scatter3d', mode : 'lines', x : line2[0], y : line2[1], z : line2[2], name : "kaka"}
+ ];
+
+
+ /*
+ let data = [];
+
+ for( let i = 0; i < 1; i++){
+ let haha = detMesh(-200, 360/6*i + 90);
+
+ console.log(haha[0], "| ", haha[1], "| ", haha[2]);
+ data.push(
+ {type : 'mesh3d', x : haha[0], y : haha[1], z : haha[2], hoverinfo: 'none', opacity : 1.0, color : "#FF0000"}
+ );
+ }*/
+
+ let layout = {
+ aspectmode: 'manual',
+ aspectratio:{ x : 0.5, y : 0.5, z : 2.5 },
+ margin: { l: 40, r: 40, b : 60, t : 40},
+ };
+
+ Plotly.newPlot('Plot_3D', data, layout);
+
+}
+
+
+GetMass();
+CalConstants();
+SetSSType();
+CalculateEZ();
+CalculateRZ();
+AdjustRecoilPos();
+CalThetaCM();
+
+//Plot3D();
diff --git a/index.html b/index.html
index 477bd6e..6d435f1 100644
--- a/index.html
+++ b/index.html
@@ -124,15 +124,15 @@