<!-- SPDX-License-Identifier: LGPL-3.0-or-later -->
<!-- Copyright (C) 2023 Sylvester Joosten, Alex Jentsch, Whitney Armstrong, Simon Gardner -->

  <define>

    <constant name="ElectronBeamEnergy" value="9*GeV"  />
    <constant name="IonBeamEnergy"      value="100*GeV" />

    <constant name="FieldScaleFactor" value="1.0"  /> <!-- only use for 275 GeV magnet config ::: ASK FF EXPERTS -->

    <comment>
       Backwards Fields
       Values taken from: https://indico.bnl.gov/event/10974/contributions/51260/
       Relevent values in the tables scale linearly with beam energy
    </comment>


    <!-- Beamline -->
    <constant name="Q1eR_Gradient"  value="-0.739741*tesla/meter/GeV*ElectronBeamEnergy"/>
    <constant name="Q2eR_Gradient"  value="0.66821*tesla/meter/GeV*ElectronBeamEnergy"/>
    <constant name="Q3eR_Gradient"  value="0.21674*tesla/meter/GeV*ElectronBeamEnergy"/>

    <constant name="B2AeR_B"        value="0.0106667*tesla/GeV*ElectronBeamEnergy"/>
    <constant name="B2BeR_B"        value="0.0132222*tesla/GeV*ElectronBeamEnergy"/>

    <!-- Luminosity System -->
    <constant name="LumiSweepMag_B"     value="1/1.2*tesla"/>
    <constant name="LumiAnalyzerMag_B"  value="1/1.2*tesla"/>

    <comment>
       Forward Fields
       These are the values from Scott's tables for the 50cm shifted hadron FF 100 GeV proton lattice. Notice they mostly do not scale with beam energy in the staightforward way from the previous lattice.
    </comment>

    <constant name="B0PF_GradientMax"  value="-8.12238283*tesla/m"/>
    <constant name="B0APF_GradientMax" value="0.0*tesla/m"/>
    <constant name="Q1APF_GradientMax" value="-44.08979*tesla/m"/>
    <constant name="Q1BPF_GradientMax" value="0.0*tesla/m"/>
    <constant name="Q2PF_GradientMax"  value="12.7764668*tesla/m"/>
    <constant name="B1PF_GradientMax"  value="0.0*tesla/m"/>
    <constant name="B1APF_GradientMax" value="0.0*tesla/m"/>
    <constant name="B2PF_GradientMax"  value="0.0*tesla/m"/>
    <constant name="Q0EF_GradientMax" value="-4.774134545*tesla/m"/>
    <constant name="Q1EF_GradientMax"  value="2.127596364*tesla/m"/>

    <constant name="B0PF_Bmax"  value="1.1840539*tesla"/>
    <constant name="B0APF_Bmax" value="-1.1588921*tesla"/>
    <constant name="Q1APF_Bmax" value="0.0*tesla"/>
    <constant name="Q1BPF_Bmax" value="0.0*tesla"/>
    <constant name="Q2PF_Bmax"  value="0.0*tesla"/>
    <constant name="B1PF_Bmax"  value="1.29077532*tesla"/>
    <constant name="B1APF_Bmax" value="0.98178068*tesla"/>
    <constant name="B2PF_Bmax" value="-1.74139319*tesla"/>
    <constant name="Q0EF_BMax" value="0.0*tesla"/>
    <constant name="Q1EF_BMax"  value="0.0*tesla"/>

    <comment>
    These are the ten-sigma cuts for the Roman pots, translated to the physical layout we currently have.
    They are not perfectly ten-sigma for reasons of physical geometry.
    </comment>

    <constant name="offset_central_RP_section" value ="0.71*cm"/>
    <constant name="offset_intermediate_1_RP_section" value ="0.55*cm"/>
    <constant name="offset_intermediate_2_RP_section" value ="0.0*cm"/>
    <constant name="offset_outer_RP_section" value ="0.0*cm"/>

  </define>
