Olivier ShelbayaThe package xml2optr (Paul Jung) can be used to generate segments of the /acc database into useable transoptr sy.f stacks. To generate a transoptr sequence from /acc-dB values with xml2optr (see /acc-utilities/xml2optr), go to the desired output folder on your system and type:
xml2optr -tf tf/hebt-dragon.xml ios-mws-hebt2-dragon
where you've also created a tf/ subfolder which contains the tunefile hebt-dragon.xml. This will automatically generate sy.f, data.dat as specified in the xml tunefile. All requisite field maps for RF elements will also be placed in the folder. Assuming proper transoptr configuration on your system, you can then simply type optr, and TRANSOPTR should execute successfully.
Most common issue preventing optr: mode-4 while trying to run RF elements. Change optr to mode-5 (3D space charge), by editing character 2, line 2 in data.dat to read 5 instead of 4.
File hebt-dragon.xml:
<root xmlns:xi='http://www.w3.org/2001/XInclude' path="ios-mws-hebt2-dragon">
<notes>These starting parameters correspond to what's in TRI-BN-19-06, which originally was coded in as a twissmatch.
I've run it with vective(1) and extracted the sigma matrix elements to define the numbers below.
see conf. prof. below for output DTL long. emit
http://accelconf.web.cern.ch/AccelConf/pac97/papers/pdf/5W036.PDF
LANA simulations predict ez = 1.7 pi keV/u*ns (awful units)
November 24, 2019
Starting s-parameters are from TRI-BN-19-06 for HEBT-Q1 sequence.
</notes>
<optr s11="0.21435*cm" s22="6.44178*mrad" s33="0.543046*cm" s44="9.23282*mrad" s55="1.0*cm" s66="0.01*rad" r12="0.794" r34="0.986" r56="-0.500" start="startOf_t3d_tune_hebt" end="HEBT2:IV8"/>
<tune mass="21.9914*u" chargestate="4.0" energy="4.3983*MeV"/>
<xi:include href='tune/ios-mws-hebt2-dragon_design200.xml' xpointer="xpointer(//root/tune)" parse="xml"/>
</root>