diff --git a/.gitignore b/.gitignore
index 19ca352..459d9a5 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,2 +1,3 @@
 files/*
 !files/.gitkeep
+*.txt
diff --git a/index.html b/index.html
index 976b47c..35beedf 100644
--- a/index.html
+++ b/index.html
@@ -171,7 +171,7 @@
 
   <!-- ////////////////////////////////////////// -->
   <article>
-      <iframe name="uploaded" style="border:none;" width="100%" height= "100%"></iframe>
+      <iframe name="uploaded" style="border:none;" width="100%" height= "100%" src="instruction.html"></iframe>
   </article>
 </section>
 
diff --git a/instruction.html b/instruction.html
index 370e8b7..5fe8119 100644
--- a/instruction.html
+++ b/instruction.html
@@ -26,8 +26,8 @@
   </ul>
 </p>
 
-<h4>This page is created and hosted by FSU. Contact Ryan Tang (rtang at fsu.edu)</h4>
-<img src="FSU_logo_640.png" width="300">
+<h4>This page is created and hosted by FSU. Contact: Ryan Tang (rtang at fsu.edu)</h4>
+<img src="FSU_logo_640.png" width="150">
 
 <h2> Todo : </h2>
 <ul>
diff --git a/simpleInput.py b/simpleInput.py
new file mode 100755
index 0000000..5d59e50
--- /dev/null
+++ b/simpleInput.py
@@ -0,0 +1,296 @@
+#!/usr/bin/env /usr/bin/python3
+import os
+import cgi, cgitb
+import subprocess
+
+
+def findName(A, Z):
+  result = subprocess.run(["../Cleopatra/Isotope", str(A), str(Z)], stdout=subprocess.PIPE).stdout.decode('utf-8')
+  ppp1 = result.find('31m')
+  ppp2 = result.find('[m')
+  return result[ppp1+3:ppp2-1]
+
+
+form = cgi.FieldStorage()
+
+beam_A = form.getvalue('beam_A')
+beam_Z = form.getvalue('beam_Z')
+beam_Jpi = form.getvalue('beam_Jpi')
+beam_Ex = form.getvalue('beam_Ex')
+
+target_A = form.getvalue('target_A')
+target_Z = form.getvalue('target_Z')
+
+lRecoil_A = form.getvalue('lRecoil_A')
+lRecoil_Z = form.getvalue('lRecoil_Z')
+
+KEA = form.getvalue('KEA')
+BField = form.getvalue('BField')
+posArray = form.getvalue('posArray')
+posRecoil = form.getvalue('posRecoil')
+
+isDWBA = form.getvalue('DWBA');
+if isDWBA == "On" :
+  isDWBA = True
+else:
+  isDWBA = False
+
+SSType = form.getvalue('SSType');
+
+ExList = form.getvalue('Ex');
+JpiList = form.getvalue('Jpi');
+OrbList = form.getvalue('Orb');
+
+op1 = form.getvalue('op1')
+op2 = form.getvalue('op2')
+
+plot = form.getvalue('plot')
+gate = form.getvalue('gate')
+
+#====== ch directory
+os.chdir(r"files")
+#================== write reactionConfig.txt
+
+fn1 = "sim_reactionConfig.txt"
+
+f = open(fn1, "w")
+f.write("%s //beam_A\n" % beam_A)
+f.write("%s //beam_Z\n" % beam_Z)
+f.write("%s //target_A\n" % target_A)
+f.write("%s //target_Z\n" % target_Z)
+f.write("%s //recoil_light_A\n" % lRecoil_A)
+f.write("%s //recoil_light_Z\n" % lRecoil_Z)
+f.write("%.3f     //beam_energy_in_MeV/u\n" % float(KEA))
+f.write("0.000  //beam_energy_sigma_in_MeV/u\n")
+f.write("0.000  //beam_angle_in_mrad\n")
+f.write("0.000  //beam_emittance_in_mrad\n")
+f.write("0.000  //beam_x-offset_in_mm\n")
+f.write("0.000  //beam_y-offset_in_mm\n")
+f.write("100000  //number_of_event_generated\n")
+f.write("false  //isTargetScattering\n")
+f.write("0.913  //Target_density_in_g/cm3\n")
+f.write("2.2e-4  //Target_thickness_in_cm\n")
+f.write("../SRIM/20F_in_CD2.txt  //stopping_power_for_beam\n")
+f.write("../SRIM/3H_in_CD2.txt   //stopping_power_for_light_recoil\n")
+f.write("../SRIM19F_in_CD2.txt  //stopping_power_for_heavy_recoil\n")
+f.write("false     //isDacay\n")
+f.write("32        //decayNucleus_A\n")
+f.write("14        //decayNucleus_Z\n")
+f.write("false     //isReDo\n")
+f.write("%s        //excitation_energy_of_A[MeV]\n" % beam_Ex)
+f.write("#===== end of file")
+f.close()
+
+#================== write detectorGeo.txt
+bore="462.5"
+prepDistant="11.5"
+detWidth="10.0"
+detLength="50.0"
+nSide="4"
+pos=['58.6', '117.9', '176.8', '235.8', '294.0']
+
+if SSType == "SOLARIS" :
+  nSide="6"
+  pos=['58.6', '117.9', '176.8', '235.8', '294.0', '353.2', '412.1', '471.0', '530.0']
+
+
+fn2 = "sim_detectorGeo.txt"
+f = open(fn2, "w")
+f.write("%s    //Bfield_[T]\n" % BField)
+f.write("0.00  //Bfield_direction_to_z-axis_[deg]_should_not_use\n")
+f.write("%s    //bore_[mm]\n" % bore)
+f.write("%s    //distance_from_axis_[mm]\n" % prepDistant)
+f.write("%s     //width_of_detector_[mm]\n" % detWidth)
+f.write("%s       //length_of_detector_[mm]\n" % detLength)
+f.write("%s      //recoil_position_+_for_downstream_[mm]\n" % posRecoil)
+f.write("10.0     //inner_radius_of_recoil_detector_[mm]\n")
+f.write("40.2     //outter_radius_of_recoil_detector_[mm]\n")
+f.write("false    //is_coincident_with_recoil\n")
+f.write("0        //Recoil_1_position_[mm]_when_0_disable_tree_branch\n")
+f.write("0        //Recoil_2_position_[mm]\n")
+f.write("0.00   //Elum_1_position_[mm]_(just_another_recoil_detector_but_for_light_recoil)\n")
+f.write("0.00   //Elum_2_position_[mm]_when_Elum=0_disable_tree_branch\n")
+f.write("0       //support_length_[mm]\n")
+f.write("%s     //first_position_-_for_upstream_[mm]\n" % posArray)
+f.write("0.03     //energy_resolution_of_PSD_array_[MeV]\n")
+f.write("1.00     //position_resolution_of_PSD_array_[mm]\n")
+f.write("Out           //detector_facing_Out_or_In\n")
+f.write("%s        //number_of_detector_as_same_side\n" % nSide)
+f.write("0.00     //1st_detector_near_position_in_reference_to_det6_[mm]\n")
+for a in pos:
+  f.write("%s\n" % a)
+f.write("#============= end of file")
+f.close()
+
+#================== if DWBA, write DWBA.txt, else write Ex.txt
+fn4 = "sim_DWBA.txt"
+beamName=""
+hRecoil_A = int(beam_A) + int(target_A) - int(lRecoil_A)
+hRecoil_Z = int(beam_Z) + int(target_Z) - int(lRecoil_Z)
+hRecoilName=""
+
+targetName=""
+lRecoilName=""
+
+if target_A == "1" and target_Z == "0":
+   targetName="n"
+if target_A == "1" and target_Z == "1":
+   targetName="p"
+if target_A == "2" and target_Z == "1":
+   targetName="d"
+if target_A == "3" and target_Z == "1":
+   targetName="t"
+if target_A == "3" and target_Z == "2":
+   targetName="3He"
+if target_A == "4" and target_Z == "2":
+   targetName="a"
+
+
+if lRecoil_A == "1" and lRecoil_Z == "0":
+   lRecoilName="n"
+if lRecoil_A == "1" and lRecoil_Z == "1":
+   lRecoilName="p"
+if lRecoil_A == "2" and lRecoil_Z == "1":
+   lRecoilName="d"
+if lRecoil_A == "3" and lRecoil_Z == "1":
+   lRecoilName="t"
+if lRecoil_A == "3" and lRecoil_Z == "2":
+   lRecoilName="3He"
+if lRecoil_A == "4" and lRecoil_Z == "2":
+   lRecoilName="a"
+
+if isDWBA :
+  #===== check if the reaction DWBA support
+  if targetName=="" or lRecoilName=="" :
+    isDWBA = False
+  else:
+    #===== find the symbol of Beam
+    beamName=findName(beam_A, beam_Z)
+    hRecoilName=findName(hRecoil_A, hRecoil_Z)
+    reactionName = beamName + "(" + targetName +","+lRecoilName+")" + hRecoilName;
+
+    f = open(fn4, "w")
+
+    if isinstance(ExList, list):
+       nState = len(ExList)
+       for i in range(0, nState):
+         f.write("%s  %s   %s   %s    %.3f   %sMeV/u   %s%s\n" % (reactionName, beam_Jpi, OrbList[i], JpiList[i], float(ExList[i]), KEA, op1, op2))
+    else :
+       f.write("%s  %s   %s   %s    %.3f   %sMeV/u   %s%s\n" % (reactionName, beam_Jpi, OrbList, JpiList, float(ExList), KEA, op1, op2))
+
+    f.close()
+
+
+fnEx="sim_Ex.txt"
+if isDWBA == False:
+  f = open(fnEx, "w")
+  f.write("//Ex   rel_xsec   SF   sigma_in_MeV\n")
+  if isinstance(ExList, list):
+    for i in range(0, len(ExList)):
+      f.write("%.3f  1.0   1.0  0.000\n" % float(ExList[i]))
+  else:
+    f.write("%.3f 1.0  1.0  0.000\n" % float(ExList))
+  f.write("#==============_End_of_file")
+  f.close()
+
+#================== write PlotConfig.txt (to be done)
+fn5=""
+
+#==================== Run Simulation
+dwba_1=""
+outFile=""
+dwba_3=""
+dwba_4=""
+isPtolemy=False
+dwbaPNG=""
+dwbaRootFile=""
+dwbaExFile=""
+pngName=""
+result=""
+haha=""
+
+if isDWBA :
+  dwba_1 = subprocess.run(["../Cleopatra/InFileCreator", fn4, "0", "180", "0.1"] , stdout=subprocess.PIPE).stdout.decode('utf-8')
+  cmd = "../Cleopatra/ptolemy <%s.in> %s.out" % (fn4, fn4)
+  outFile = fn4+ ".out"
+  os.system(cmd)
+  if os.path.exists(outFile) :
+    dwba_3 = subprocess.run(["../Cleopatra/ExtractXSec", outFile, "1"], stdout=subprocess.PIPE).stdout.decode('utf-8')
+    dwbaRootFile = fn4 + ".root"
+    dwbaExFile = fn4 + ".Ex.txt"
+    dwba_4 = subprocess.run(["../Cleopatra/PlotTGraphTObjArray", dwbaRootFile, "1"], stdout=subprocess.PIPE).stdout.decode('utf-8')
+    dwbaPNG = dwba_4.splitlines().pop()
+
+
+if isDWBA :
+  result = subprocess.run(["../Cleopatra/Transfer", fn1, fn2, dwbaExFile, dwbaRootFile, "transfer.root", "reaction.dat", "0"], stdout=subprocess.PIPE).stdout.decode('utf-8')    
+else:
+  result = subprocess.run(["../Cleopatra/Transfer", fn1, fn2, fnEx, "", "transfer.root", "reaction.dat", "0"], stdout=subprocess.PIPE).stdout.decode('utf-8')
+
+
+if fn5 =="":
+   haha = subprocess.run(["../Cleopatra/PlotSimulation", "transfer.root"],  stdout=subprocess.PIPE).stdout.decode('utf-8')
+else:
+   haha = subprocess.run(["../Cleopatra/PlotSimulation", "transfer.root", fn5],  stdout=subprocess.PIPE).stdout.decode('utf-8')
+pngName=haha.splitlines().pop()
+
+
+
+print ("Content-Type: text/html\r\n\r\n")
+print ("<html>")
+print ("<style> img {max-width: 100%} </style>")
+print ("<style> body { font-family: courier, courier new, serif;} </style>")
+print ("<style> a { color: #F7CF3C;} </style>")
+print ("<body>")
+
+'''
+print("<br> %s" % beam_A)
+print("<br> %s" % beam_Z)
+print("<br> %s" % beam_Jpi)
+print("<br> %s" % beam_Ex)
+print("<br> %s" % target_A)
+print("<br> %s" % target_Z)
+print("<br> %s" % lRecoil_A)
+print("<br> %s" % lRecoil_Z)
+print("<br> %s" % KEA)
+print("<br> %s" % BField)
+print("<br> %s" % posArray)
+print("<br> %s" % posRecoil)
+print("<br> %s" % SSType)
+print("<br> do DWBA : |%s|" % isDWBA)
+print("<br> %s" % ExList)
+print("<br> %s" % JpiList)
+print("<br> %s" % OrbList)
+print("<br> %s" % op1)
+print("<br> %s" % op2)
+print("<br> %s" % plot)
+print("<br> %s" % gate)
+'''
+
+
+if dwbaPNG != "" :
+  print ("<img src=\"files/%s\">" % dwbaPNG)
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">%s</a>" % (fn4, fn4, fn4))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">%s</a>" % (fn4+".in", fn4+".in", fn4+".in"))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">%s</a>" % (fn4+".out", fn4+".out", fn4+".out"))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">%s</a>" % (fn4+".Xsec.txt", fn4+".Xsec.txt", fn4+".Xec.txt"))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">%s</a>" % (fn4+".root", fn4+".root", fn4+".root"))
+  print ("<br>")
+
+if pngName != "" :
+  print ("<img src=\"files/%s\">" % pngName)
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">reactionConfig file</a>" % (fn1, fn1))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">detectorGeo file</a>" % (fn2, fn2))
+  print ("<br> Downlaod the <a href=\"files/%s\" download=\"%s\">Ex file</a>" % (fnEx, fnEx))
+  print ("<br> Downlaod the simulation <a href=\"files/transfer.root\" download=\"transfer.root\">root file</a>")
+
+print ("<br><br><br>#################################### programs text output")
+print ("<br><span style=\"white-space: pre-line\"> %s </span>" % dwba_1)
+print ("<br><span style=\"white-space: pre-line\"> %s </span>" % dwba_3)
+print ("<br><span style=\"white-space: pre-line\"> %s </span>" % dwba_4)
+print ("<br><span style=\"white-space: pre-line\"> %s </span>" % result)
+print ("<br><span style=\"white-space: pre-line\"> %s </span>" % haha)
+print ("</body>")
+print ("</html>")
+
+
diff --git a/simpleSim.html b/simpleSim.html
index f9e0374..86f9e21 100644
--- a/simpleSim.html
+++ b/simpleSim.html
@@ -1,3 +1,4 @@
+
 <!DOCTYPE html>
 <html>
 <head>
@@ -18,37 +19,37 @@
   <h3>Reaction:</h3>
 
     <table>
+      <tr>
+        <td>Beam Energy</td>
+        <td></td>
+        <td><input type = "text" name = "KEA" size="5" value="10"/></td>
+        <td style="text-align:left">MeV/u</td>
+      </tr>
       <tr>
         <td>Beam </td>
         <td style="text-align:right">A:</td>
-        <td><input type = "text" name = "beam_A" size="5"/></td>
+        <td><input type = "text" name = "beam_A" size="5" value="16"/></td>
         <td style="text-align:right">Z:</td>
-        <td><input type = "text" name = "beam_Z" size="5"/></td>
+        <td><input type = "text" name = "beam_Z" size="5" value="8"/></td>
         <td style="text-align:right">J-pi</td>
-        <td><input type = "text" name = "beam_Jpi" size="5"/><td>
+        <td><input type = "text" name = "beam_Jpi" size="5" value="0+"/><td>
         <td style="text-align:right">Ex:</td>
-        <td><input type = "text" name = "beam_Ex" size="5"/></td>
+        <td><input type = "text" name = "beam_Ex" size="5" value="0.00"/></td>
         <td style="text-align:left">MeV</td>
       </tr>
       <tr>
         <td>Target</td>
         <td style="text-align:right">A:</td>
-        <td><input type = "text" name = "target_a" size="5"/></td>
+        <td><input type = "text" name = "target_A" size="5" value="2"/></td>
         <td style="text-align:right">Z:</td>
-        <td><input type = "text" name = "target_z" size="5"/></td>
+        <td><input type = "text" name = "target_Z" size="5" value="1"/></td>
       </tr>
       <tr>
         <td>Light recoil</td>
         <td style="text-align:right">A:</td>
-        <td><input type = "text" name = "lRecoil_a" size="5"/></td>
+        <td><input type = "text" name = "lRecoil_A" size="5" value="1"/></td>
         <td style="text-align:right">Z:</td>
-        <td><input type = "text" name = "lRecoil_z" size="5"/></td>
-      </tr>
-      <tr>
-        <td>Beam Energy</td>
-        <td></td>
-        <td><input type = "text" name = "KEA" size="5"/></td>
-        <td style="text-align:left">MeV/u</td>
+        <td><input type = "text" name = "lRecoil_Z" size="5" value="1"/></td>
       </tr>
     </table>
   
@@ -56,29 +57,30 @@
   
     <table>
       <tr>
-        <input type="radio" id="helios" name="SSType" value="HELIOS"/><label for="helios">HELIOS</label>
-        <br><input type="radio" id="solaris" name="SSType" value="SOLARIS"/><label for="solaris">SOLARIS</label>
+        <input type="radio" name="SSType" value="HELIOS"/>HELIOS
+        <br>
+        <input type="radio" name="SSType" value="SOLARIS" checked="checked"/>SOLARIS
       </tr>
       <tr>
         <td>B-field </td>
-        <td><input type = "text" name = "Bfield" size="5"/></td>
-        <td style="text-align:left">T</td>
+        <td><input type = "text" name = "BField" size="5" value="-2.5"/></td>
+        <td style="text-align:left">T (minus sign = field point to upstream)</td>
       </tr>
       <tr>
         <td>Array Pos.</td>
-        <td><input type = "text" name = "posArray" size="5"/></td>
+        <td><input type = "text" name = "posArray" size="5" value="-100"/></td>
         <td style="text-align:left">mm (negative for upstream)</td>
       </tr>
       <tr>
         <td>Recoil Pos.</td>
-        <td><input type = "text" name = "posRecoil" size="5"/></td>
+        <td><input type = "text" name = "posRecoil" size="5" value="500"/></td>
         <td style="text-align:left">mm (negative for upstream)</td>
       </tr>
     </table>
 
 <h3>DWBA and Ex List:</h3>
 
-  <input type="checkbox" id="DWBA" name="DWBA" value="On"/><label for="DWBA">Cal. DWBA</label>
+  <input type="checkbox" name="DWBA" value="On"/>Cal. DWBA
   
   <p></p>
 
@@ -89,9 +91,9 @@
       <td> Orbital </td>
     </tr>
     </tr>
-     <td><input type="text" size="5" value="0"/></td>
-     <td><input type="text" size="5" value="0+"/></td>
-     <td><input type="text" size="5" value="0d3/2"/></td>
+     <td><input type="text" name="Ex" size="5" value="0"/></td>
+     <td><input type="text" name="Jpi" size="5" value="3/2+"/></td>
+     <td><input type="text" name="Orb" size="5" value="0d3/2"/></td>
     </tr>
     <tr>
       <td><button type="button" onclick="addRow()">Add Ex</button></td>
@@ -106,25 +108,58 @@
     <tr>
       <td>Incoming Channel</td>
       <td>
-        <select name="op1" id="op1">
-          <option value="An">An.</option>
-          <option value="Kai">Kai.</option>
+        <select name="op1">
+          <option value="A" selected>D | An & Cai (2006), E < 183, 12 < A < 238</option>
+          <option value="H">D | Han, Shi, & Shen, (2006), E < 200 | 12 < A < 209</option>
+          <option value="B">D | Bojowald et al., (1988), 50 < E < 80 | 27 < A < 208</option>
+          <option value="K">P | Koning & Delaroche, (2009), E < 200 | 24 < A < 209 | Iso.Dep.</option>
+          <option value="V">P | Varner et al., (1991), 16 < E < 65 | 4 < A < 209</option>
+          <option value="M">P | Menet et al., (1971), 30 < E < 60 | 40 < A </option>
+          <option value="x">A=3 | Xu, Guo, Han, & Shen, (2011), E < 250 | 20 < A < 209 </option>
+          <option value="l">A=3 | Liang, Li, & Cai, (2009), E < 270 | All masses </option>
+          <option value="x">A=4 | Su & Han, (2015), E < 398 | 20 < A < 209 </option>
         </select> 
       <td>
     </tr>
     <tr>
       <td>Outgoing Channel</td>
-      <td> <td>
+      <td>
+        <select name="op2">
+          <option value="A">D | An & Cai (2006), E < 183, 12 < A < 238</option>
+          <option value="H">D | Han, Shi, & Shen, (2006), E < 200 | 12 < A < 209</option>
+          <option value="B">D | Bojowald et al., (1988), 50 < E < 80 | 27 < A < 208</option>
+          <option value="K" selected>P | Koning & Delaroche, (2009), E < 200 | 24 < A < 209 | Iso.Dep.</option>
+          <option value="V">P | Varner et al., (1991), 16 < E < 65 | 4 < A < 209</option>
+          <option value="M">P | Menet et al., (1971), 30 < E < 60 | 40 < A </option>
+          <option value="x">A=3 | Xu, Guo, Han, & Shen, (2011), E < 250 | 20 < A < 209 </option>
+          <option value="l">A=3 | Liang, Li, & Cai, (2009), E < 270 | All masses </option>
+          <option value="x">A=4 | Su & Han, (2015), E < 398 | 20 < A < 209 </option>
+        </select> 
+      <td>
     </tr>
   </table>
 
-<h3> Plot config:</h3>
+<h3> Plot config (still working on):</h3>
+
+  <input type="checkbox" name="plot" value="pEZ" checked/>E vs Z<br>
+  <input type="checkbox" name="plot" value="pExCal" checked/>Ex (cal.)<br>
+  <input type="checkbox" name="plot" value="pThetaCM" checked/>ThetaCM<br>
+  <input type="checkbox" name="plot" value="pThetaCM_Z" checked/>ThetaCM vs Z<br>
+  <input type="checkbox" name="plot" value="pRecoilXY" checked/>Recoil X vs Y<br>
+  <input type="checkbox" name="plot" value="pRecoilRThetaCM"/>Recoil-R vs ThetaCM<br>
+  <input type="checkbox" name="plot" value="pRecoilRZ"/>Recoil R vs Z<br>
+  <input type="checkbox" name="plot" value="pArrayXY"/>Array X vs Y<br>
+
+  <p></p>
+  
+  <input type="checkbox" name="gate" value="hit==1" checked/>Array Hit<br>
+  <input type="checkbox" name="gate" value="loop<=1" checked/>Loop = 1<br>
+  <input type="checkbox" name="gate" value="thetaCM>10" checked/> ThetaCM > 10 deg<br>
 
-  <input type="checkbox" id="pcExZ" name="pcExZ" value="On"/><label for="pxExZ">e vs z</label>
 
 <p></p>
 
-<input type = "submit" value = "Submit and Run Simulation" />
+<input type = "submit" value = "Submit and Run Simulation" style="width:200px;height:60px;" />
 </form>
 
 
@@ -133,9 +168,9 @@ function addRow() {
   var table = document.getElementById("ExTable"); 
   var nRow = table.rows.length;
   var row = table.insertRow(nRow-1);
-  row.innerHTML = '<td><input type="text" size="5"/></td> \
-                   <td><input type="text" size="5"/></td> \
-                   <td><input type="text" size="5"/></td>';
+  row.innerHTML = '<td><input type="text" name="Ex" size="5"/></td> \
+                   <td><input type="text" name="Jpi" size="5"/></td> \
+                   <td><input type="text" name="Orb" size="5"/></td>';
 }
 
 function deleteRow(){