Difference between revisions of "AP-CEBAF-Optics-2015-02-28"

From CASA Wiki
Jump to: navigation, search
 
(57 intermediate revisions by 2 users not shown)
Line 1: Line 1:
This is a wiki page for references, information, and procedures about Feb 27-Mar 1 work to transversely match CEBAF Twiss parameters at the injector (~0L07), 1E03, and xE01 locations.
+
This is a wiki page for references, information, and procedures about Feb 27-Mar 4 2015 work to transversely match CEBAF Twiss parameters at the injector (~0L07), 1E03, and xE01 locations.
 +
 
 +
[http://casa.jlab.org/bteam/a.pdf example PDF link for Audrey]
  
 
== Philosophy ==
 
== Philosophy ==
The philosophy was developed at a 15:00 meeting in MCC, including Arne Freyberger, Mike Spata, Yves Roblin, Mike TIefenback, Jay Benesch, Todd Satogata, and Operations:
+
The philosophy was developed at a Feb 28 15:00 meeting in MCC, including Arne Freyberger, Mike Spata, Yves Roblin, Mike Tiefenback, Jay Benesch, Todd Satogata, and Operations:
 
* use zigzag harp scans to identify quadrupole ranges for good optics/emittance measurements at each harp.
 
* use zigzag harp scans to identify quadrupole ranges for good optics/emittance measurements at each harp.
 
* modify qsUtility configuration files in the golden area and acquire optics/emittance measurements with traditional qsUtility scans.
 
* modify qsUtility configuration files in the golden area and acquire optics/emittance measurements with traditional qsUtility scans.
Line 20: Line 22:
 
* confirm harp is working and set harp parameters
 
* confirm harp is working and set harp parameters
 
* perform zigzag measurements to establish scan quadrupoles and ranges
 
* perform zigzag measurements to establish scan quadrupoles and ranges
* edit qsUtility configuration files
+
* edit qsUtility configuration files if necessary
 +
** configuration files are in /a/opsdata/fiefdata/fileio/q/qsUtility/config/F3-0/GoldenTemplates
 
* perform qsUtility scans
 
* perform qsUtility scans
* evaluate qsUtility data to get measured Twiss parameters and remittances
+
* evaluate qsUtility data to get measured Twiss parameters and emittances
 
* back-propagate measured Twiss parameters to start of upstream xS04 quadrupole with elegant
 
* back-propagate measured Twiss parameters to start of upstream xS04 quadrupole with elegant
 
** A script to do this is in ~satogata/backpropagate.pl
 
** A script to do this is in ~satogata/backpropagate.pl
 
* use qsUtility betatron match module to calculate xS04-xS10 match
 
* use qsUtility betatron match module to calculate xS04-xS10 match
 +
** configuration files are in /a/opsdata/fiefdata/fileio/q/qsUtility/match_config/F3-0
 
* evaluate and install match
 
* evaluate and install match
 
* re-scan harp and evaluate match effectiveness
 
* re-scan harp and evaluate match effectiveness
 +
* After match is installed, one has to set the recombiner quads back to design in order to proceed to the next arc.
 
* perform confirmation fopt
 
* perform confirmation fopt
  
 
== Data and References ==
 
== Data and References ==
* [https://docs.google.com/spreadsheets/d/1AMVnjOjYZAqxLSQI_bzgt4cbmZlRY0DIh6uUaYCadPk  Spring 2015 Optics Spreadsheet ]
+
* [https://docs.google.com/spreadsheets/d/1AMVnjOjYZAqxLSQI_bzgt4cbmZlRY0DIh6uUaYCadPk  Spring 2015 Optics Spreadsheet ] (Google docs)
 
* 0L injector match (E=107.2 MeV, BDL=536.365*K1 for L=0.15m magnets)
 
* 0L injector match (E=107.2 MeV, BDL=536.365*K1 for L=0.15m magnets)
** 0L07: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQD0L06 horizontal scan centered near 1800, range +/-270 (BDL 3.35+/-0.5)
+
** 0L07 H: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQD0L06 horizontal scan centered near 1800G, range +/-270G (BDL 3.35+/-0.5)
** 0L08: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB0L07 vertical scan centered near -3600, range +/-500 (BDL -6.71+/-1.0)
+
** 0L08 V: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB0L07 vertical scan centered near -3600G, range +/-500G (BDL -6.71+/-1.0)
* 1E03 match (E=1057.2 MeV, BDL=5289.659*K1 for L=0.15m magnets):
+
** Match calculation results and loading: [https://logbooks.jlab.org/entry/3322599 e3322599] using C-laser, IHA0L07_2015-02-28_23:26 and IHA0L07_2015-03-01_00:01 datasets.
** 1E03: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB1E02 horizontal scan centered near 3100, range +/-700 (BDL 0.586+/-0.133)
+
** 2015-03-10 C laser after injector work, before 750 MHz separator work: [https://logbooks.jlab.org/entry/3324993 horizontal analysis], [https://logbooks.jlab.org/entry/3324999 horizontal analysis 2], [https://logbooks.jlab.org/entry/3325015 vertical analysis], [https://logbooks.jlab.org/entry/3325021 NL-inj match], [https://logbooks.jlab.org/entry/3325024 horizontal analysis after match], and [https://logbooks.jlab.org/entry/3325028 vertical analysis after match].
** 1E03: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB1E01 vertical scan centered near -2900, range +/-1200 (BDL -0.548+/-0.227)
+
* 1E03 match (E=1057.2 MeV, BDL=5289.659*K1 for L=0.15m magnets, be careful MQN1S04 is L=0.30m, others are L=0.15m):
* 2E01 match (E=2007.2 MeV, BDL=10042.95*K1 for L=0.15m magnets):
+
** 1E03 H: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB1E02 horizontal scan centered near 3100G, range +/-700G (BDL 0.586+/-0.133)
* 3E01 match (E=2957.2 MeV, BDL=14796.24*K1 for L=0.15m magnets):
+
** 1E03 V: [https://logbooks.jlab.org/entry/3322505 zigzags] -> MQB1E01 vertical scan centered near -2900G, range +/-1200G (BDL -0.548+/-0.227)
* 4E01 match (E=3907.2 MeV, BDL=19549.52*K1 for L=0.15m magnets):
+
** Analysis and match calculation: [https://logbooks.jlab.org/entry/3322630 e3322630]
* 5E01 match (E=4857.2 MeV, BDL=24302.81*K1 for L=0.15m magnets):
+
** Installed and confirmed match with fopt: [https://logbooks.jlab.org/entry/3322693 e3322693] [https://logbooks.jlab.org/entry/3322695 e3322695]
* 6E01 match (E=5807.2 MeV, BDL=29056.10*K1 for L=0.15m magnets):
+
* 2E01 match (E=2007.2 MeV, BDL=20085.89*K1 for L=0.30m magnets):
* 8E01 match (E=7707.2 MeV, BDL=38562.68*K1 for L=0.15m magnets):
+
** 2E01 H: [https://logbooks.jlab.org/entry/3322505 zigzags] [https://logbooks.jlab.org/entry/3322697 zigzags]-> MQC2S10 horizontal scan centered near 4600G, range +/-1000G (BDL 0.229+/-0.05)
* 9E01 match (E=8657.2 MeV, BDL=43315.97*K1 for L=0.15m magnets):
+
** 2E01 V: No zigzag solution found. Todd explored space with qsUtilityConfig and found reasonable scan of MQC2S08 K=-0.6 to 0.6 with MQC2S09 K=-0.1 m^-2 and MQC2S10 K=0.2 m^-2 [http://casa.jlab.org/bteam/2015-02-28-IHA2E01-y-scan.png (image)]
* AE01 match (E=9607.2 MeV, BDL=48069.25*K1 for L=0.15m magnets):
+
** In order to test that with the zigzag do the following (yves):
* 5C00 match (E=10557.2 MeV, BDL=52822.54*K1 for L=0.15m magnets):
+
# caput MQC2S09.BDL -2009
* 1C match (E=9607.2 MeV, fifth pass, BDL=48069.25*K1 for L=0.15m magnets):
+
# caput MQC2S10.BDL 4017
* 2C match (E=2007.2 MeV, first pass, BDL=10042.95*K1 for L=0.15m magnets):
+
# scan_script.csh run IHA2E01 MQC2S08 Y 12051.000 -12051.000 9
 +
then restore the 2S09 and 2S10 quads .
 +
** 2E01 [https://logbooks.jlab.org/entry/3322698 horizontal scan] and [https://logbooks.jlab.org/entry/3322700 analysis]
 +
** 2E01 [https://logbooks.jlab.org/entry/3322703 vertical scan] and [https://logbooks.jlab.org/entry/3322705 analysis]
 +
** 2E01 [https://logbooks.jlab.org/entry/3322706 backpropagation] and fitting
 +
* 3E01 match (E=2957.2 MeV, BDL=29592.47*K1 for L=0.30m magnets):
 +
** 3E01 H: broad scan IHA3E01_2015-02-23_22:27.sdds of MQA3S10 with MQB3E01=-0.7991 shows minimum near K1=1.1 m^-2, decent range K1=1.6 to 0.6. Put into IHA3E01-x-GOLD.xml.
 +
** 3E01 V: broad scan IHA3E01_2015-02-23_22:27.sdds of MQA3S10 with MQB3E01=-0.7991 shows minimum near K1=1.35 m^-2, decent range K1=2.1 to 0.65. Put into IHA3E01-y-GOLD.xml.
 +
** 3E01 [https://logbooks.jlab.org/entry/3322957 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/3322956 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3322958 match calculation]
 +
** 3E01 [https://logbooks.jlab.org/entry/3322968 fopt after 3E match] and [https://logbooks.jlab.org/entry/3322961 viewer walk], all looking good
 +
* 4E01 match (E=3907.2 MeV, BDL=39099.05*K1 for L=0.30m magnets):
 +
** The 02242015 4E01 scans indicate a need for vertical quad centering at MQA4S09. (done Mar 3 2015 13:30)
 +
** 4E01 H: broad scan IHA4E01_2015-02-24_00:38.sdds of MQA4S09 with MQA4S10=0.2538 and MQB4E01=-0.5713 shows minimum near K1=0.0 m^-2, decent range K1=0.5 to -0.5. Profiles look ratty.
 +
** 4E01 V: broad scan IHA4E01_2015-02-24_00:38.sdds of MQA4S09 with MQA4S10=0.2538 and MQB4E01=-0.5713 shows minimum near K1=-0.45 m^-2, decent range K1=0.9 to -1.8.
 +
** 4E01 [https://logbooks.jlab.org/entry/3322982 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/3322985 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3322994 corrected match calculation]
 +
** 4E01 [https://logbooks.jlab.org/entry/3322999 fopt after 4E match] and [https://logbooks.jlab.org/entry/3323000 viewer walk], all looking good
 +
* 5E01 match (E=4857.2 MeV, BDL=48605.63*K1 for L=0.30m magnets):
 +
** Initial zigzag scans with [https://logbooks.jlab.org/entry/3323003 design ranges] indicated a need for quad centering.
 +
** IHA5E01 harp appears broken: [https://logbooks.jlab.org/entry/3323035 e3323035]
 +
* 6E01 match (E=5807.2 MeV, BDL=58112.20*K1 for L=0.30m magnets):
 +
** 6E01 [https://logbooks.jlab.org/entry/3323043 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/333045 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3323046 match download]
 +
** 6E01 [https://logbooks.jlab.org/entry/3323051 fopt after 6E match] and [https://logbooks.jlab.org/entry/3323050 viewer walk].
 +
** fopt was sufficiently problematic (particularly in vertical), and data quality was sufficiently concerning (particularly in horizontal), that Yves and Todd decided to re-do the 6E01 zigzag measurements and match. 6S04-10 quads were restored to design before additional measurements.
 +
** 6E01 [https://logbooks.jlab.org/entry/3323057 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/3323058 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3323061 match download] (removing quads S07/8/9 from fit to improve convergence).
 +
** 6E01 [https://logbooks.jlab.org/entry/3323068 fopt after 6E match] and restoring 6A15-19 quads to design.
 +
* 7E01 has no harp (it was relocated to 1E03 several runs ago for dogleg edge focusing measurements)
 +
* 8E01 match (E=7707.2 MeV, BDL=77125.36*K1 for L=0.30m magnets):
 +
** 8E01 harp has encoder problem: [https://logbooks.jlab.org/entry/3323074 e3323074] but it could be worked around with manual intervention.
 +
** 8E01 [https://logbooks.jlab.org/entry/3323077 horizontal zigzag analysis] after manual intervention. Vertical data had a [https://logbooks.jlab.org/entry/3323084 bad range], then we were down for a 1L18 vacuum valve problem.
 +
** 8E01 [https://logbooks.jlab.org/entry/3323286 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/3323291 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3323292 propagation] and [https://logbooks.jlab.org/entry/3323293 match download]
 +
** 8E01 [https://logbooks.jlab.org/entry/3323304 fopt after 8E01 match, 8R restored to design]
 +
* 9E01 match (E=8657.2 MeV, BDL=86631.93*K1 for L=0.30m magnets):
 +
** 9E01 [https://logbooks.jlab.org/entry/3323335 horizontal zigzag analysis] and [https://logbooks.jlab.org/entry/3323339 vertical zigzag analysis] and [https://logbooks.jlab.org/entry/3323341 propagation] and [https://logbooks.jlab.org/entry/3323343 match download]
 +
** 9E01 [https://logbooks.jlab.org/entry/3323353 fopt after 9E01 match, 9R restored to design]
 +
** 9E01 [https://logbooks.jlab.org/entry/3323357 remeasure at 9E01 after match, '''beautifully on design''']
 +
* AE01 match (E=9607.2 MeV, BDL=96138.51*K1 for L=0.30m magnets):
 +
** AE01 [https://logbooks.jlab.org/entry/3323363 horizontal zigzag analysis]
 +
** AE01 vertical scans required IPMAS09 quad centering, then [https://logbooks.jlab.org/entry/3323370 the IHAAE01 IOC crashed]
 +
* 5C00 match (E=10557.2 MeV, BDL=105645.09*K1 for L=0.30m magnets):
 +
* 1C match (E=9607.2 MeV, fifth pass, BDL=96138.51*K1 for L=0.30m magnets):
 +
** data taken and match performed [https://logbooks.jlab.org/entry/3323781 horizontal and vertical data analysis].  Ranges had to be adjusted. Ideal ranges seems to be
 +
** 0.25 <K1<0.50 for MQNAT06 taking the X data.
 +
** -0.7 <K1<-1.2 for MQK1C04 taking the Y data.
 +
* 2C match (E=2007.2 MeV, first pass, BDL=20085.89*K1 for L=0.30m magnets):
 
* rayTrace data: [https://logbooks.jlab.org/entry/3322485 3322485 H] [https://logbooks.jlab.org/entry/3322488 3322488 H] [https://logbooks.jlab.org/entry/3322490 3322490 V]
 
* rayTrace data: [https://logbooks.jlab.org/entry/3322485 3322485 H] [https://logbooks.jlab.org/entry/3322488 3322488 H] [https://logbooks.jlab.org/entry/3322490 3322490 V]
 +
 +
== Useful Tools ==
 +
'''backpropagate.pl''': Todd has written a tool to backpropagate Twiss parameters from the front of a given quadrupole (where they were measured using harp scans) to the front of an earlier quadrupole (which is the start of the matching section, almost always the nS04 quadrupole here). There is no coupling, so you can do just a single plane by just using easy numbers (e.g. beta=1 and alpha=0) for the plane you are not interested in. Usage us
 +
  # ~satogata/backpropagate/backpropagate.pl <measQuad> <backQuad> <betax> <alphax> <betay> <alphay>
 +
It will generate a back-propagation elegant lattice and elegant file in your current directory, and will produce errors if you try to give it a quadrupole name that does not exist. An example of backpropagation command line and output for the 2E region is [https://logbooks.jlab.org/entry/3322706 here].
 +
 +
'''setquadstomachine''': During the backpropagation or the matching calculation, if one wants to take into account quadrupoles that are off design, one can use the script ''~roblin/tools/setquadstomachine'' which together with a snap file of the region in question produces a sdds file suitable for inclusion in the ELEGANT command file via load_parameters. usage is
 +
  # ~roblin/tools/setquadstomachine <commandfile> <snapfile>
 +
and it will make QUADparams.sdds ready for load_parameters. Read the comments in the script for further details (or page Yves).
 +
 +
'''analyzeZigZagRange.pl''': Todd has written a tool to help analyze zigzag scans for quad scan ranges. Usage is, for example,
 +
# ~satogata/tools/analyzeZigZagRange.pl IHA1E03.SA.02282015_00:28:42
 +
This tool prints out the quadrupole that is being scanned, the list of scans and corresponding BDL and K values for that quad, and uses gnuplot to plot the zigzag scan data. One can then count the peaks and look up the corresponding BDL and K values for use in qsUtility configuration files.
 +
 +
'''zigzag2qsUtility''': Dennis has a tool that converts zigzag scans to input files for qsUtility, though at present the qsUtility file that is generated needs to have the quadrupole K values added manually. An example for IHA3E01 is
 +
# zigzag2qsUtility -scan_file=IHA3E01.SA.03032015_10:48:30 -show_plots -exclusion_zone=0.15 -plane=X -lattice_name=test -momentum=2200 -start_element=MQA3S10 -pass=2

Latest revision as of 14:56, 19 July 2016

This is a wiki page for references, information, and procedures about Feb 27-Mar 4 2015 work to transversely match CEBAF Twiss parameters at the injector (~0L07), 1E03, and xE01 locations.

example PDF link for Audrey

Philosophy

The philosophy was developed at a Feb 28 15:00 meeting in MCC, including Arne Freyberger, Mike Spata, Yves Roblin, Mike Tiefenback, Jay Benesch, Todd Satogata, and Operations:

  • use zigzag harp scans to identify quadrupole ranges for good optics/emittance measurements at each harp.
  • modify qsUtility configuration files in the golden area and acquire optics/emittance measurements with traditional qsUtility scans.
  • use qsUtility to analyze the scans and get measured horizontal/vertical Twiss parameters at the start of the respective scan quads.
  • back-propagate measured Twiss parameters to the start of the upstream xS04 quadrupole.
  • use the qsUtility betatron match module (which uses ced2elegant and elegant) to calculate the match, with match point at the harp, using xS04-xS10 quadrupoles for matching to distribute matching effort.
  • evaluate and install the match.
  • re-scan the harp and compare beam sizes to sqrt(emit_meas * beta_design)

Prerequisites at Each Measurement Location

  • 4-8 uA tune beam to first dumplette downstream of match point.
  • quad centering of BPM/quadrupole pairs for qsUtility scan quadrupoles.
  • clean steering in region of measurement and downstream to maximize aperture, minimize beam losses.
  • clean ORFP dispersion in region of measurement and correction. In particular, dispersion in both planes must be consistent with zero at scan quadrupoles and harp.

Procedure

  • confirm harp is working and set harp parameters
  • perform zigzag measurements to establish scan quadrupoles and ranges
  • edit qsUtility configuration files if necessary
    • configuration files are in /a/opsdata/fiefdata/fileio/q/qsUtility/config/F3-0/GoldenTemplates
  • perform qsUtility scans
  • evaluate qsUtility data to get measured Twiss parameters and emittances
  • back-propagate measured Twiss parameters to start of upstream xS04 quadrupole with elegant
    • A script to do this is in ~satogata/backpropagate.pl
  • use qsUtility betatron match module to calculate xS04-xS10 match
    • configuration files are in /a/opsdata/fiefdata/fileio/q/qsUtility/match_config/F3-0
  • evaluate and install match
  • re-scan harp and evaluate match effectiveness
  • After match is installed, one has to set the recombiner quads back to design in order to proceed to the next arc.
  • perform confirmation fopt

Data and References

  • Spring 2015 Optics Spreadsheet (Google docs)
  • 0L injector match (E=107.2 MeV, BDL=536.365*K1 for L=0.15m magnets)
  • 1E03 match (E=1057.2 MeV, BDL=5289.659*K1 for L=0.15m magnets, be careful MQN1S04 is L=0.30m, others are L=0.15m):
    • 1E03 H: zigzags -> MQB1E02 horizontal scan centered near 3100G, range +/-700G (BDL 0.586+/-0.133)
    • 1E03 V: zigzags -> MQB1E01 vertical scan centered near -2900G, range +/-1200G (BDL -0.548+/-0.227)
    • Analysis and match calculation: e3322630
    • Installed and confirmed match with fopt: e3322693 e3322695
  • 2E01 match (E=2007.2 MeV, BDL=20085.89*K1 for L=0.30m magnets):
    • 2E01 H: zigzags zigzags-> MQC2S10 horizontal scan centered near 4600G, range +/-1000G (BDL 0.229+/-0.05)
    • 2E01 V: No zigzag solution found. Todd explored space with qsUtilityConfig and found reasonable scan of MQC2S08 K=-0.6 to 0.6 with MQC2S09 K=-0.1 m^-2 and MQC2S10 K=0.2 m^-2 (image)
    • In order to test that with the zigzag do the following (yves):
# caput MQC2S09.BDL -2009
# caput MQC2S10.BDL 4017
# scan_script.csh run IHA2E01 MQC2S08 Y 12051.000 -12051.000 9

then restore the 2S09 and 2S10 quads .

Useful Tools

backpropagate.pl: Todd has written a tool to backpropagate Twiss parameters from the front of a given quadrupole (where they were measured using harp scans) to the front of an earlier quadrupole (which is the start of the matching section, almost always the nS04 quadrupole here). There is no coupling, so you can do just a single plane by just using easy numbers (e.g. beta=1 and alpha=0) for the plane you are not interested in. Usage us

 # ~satogata/backpropagate/backpropagate.pl <measQuad> <backQuad> <betax> <alphax> <betay> <alphay>

It will generate a back-propagation elegant lattice and elegant file in your current directory, and will produce errors if you try to give it a quadrupole name that does not exist. An example of backpropagation command line and output for the 2E region is here.

setquadstomachine: During the backpropagation or the matching calculation, if one wants to take into account quadrupoles that are off design, one can use the script ~roblin/tools/setquadstomachine which together with a snap file of the region in question produces a sdds file suitable for inclusion in the ELEGANT command file via load_parameters. usage is

 # ~roblin/tools/setquadstomachine <commandfile> <snapfile>

and it will make QUADparams.sdds ready for load_parameters. Read the comments in the script for further details (or page Yves).

analyzeZigZagRange.pl: Todd has written a tool to help analyze zigzag scans for quad scan ranges. Usage is, for example,

# ~satogata/tools/analyzeZigZagRange.pl IHA1E03.SA.02282015_00:28:42

This tool prints out the quadrupole that is being scanned, the list of scans and corresponding BDL and K values for that quad, and uses gnuplot to plot the zigzag scan data. One can then count the peaks and look up the corresponding BDL and K values for use in qsUtility configuration files.

zigzag2qsUtility: Dennis has a tool that converts zigzag scans to input files for qsUtility, though at present the qsUtility file that is generated needs to have the quadrupole K values added manually. An example for IHA3E01 is

# zigzag2qsUtility -scan_file=IHA3E01.SA.03032015_10:48:30 -show_plots -exclusion_zone=0.15 -plane=X -lattice_name=test -momentum=2200 -start_element=MQA3S10 -pass=2