Setup

Load packages and set their options

Working directory and temporary working space

Utility functions for this application

Finding the strength name for a quadrupole

This function is a bit flaky and depends on the file formatting.  Nevertheless very useful.

Adding quadrupole strength limits

Matching constraints

Better implementation of matching constraints.  This was a bit of a headache.

Inequalities are troublesome in our old approach:

CONSTRAINT, SEQUENCE=LHCB1, RANGE=IP, Y=100., PY=PY1IP;

This is a new version, still a bit tricky.  It allows constraints to be given as lists of 2 or 3 elements.  See examples following.  Nevertheless the main point is that we can give expressions that evaluate to numbers as the last element of each constraint list.  

No checking of s values in the following

CONSTRAINT,SEQUENCE=LHCB1,RANGE=IP2,BETX>100, BETY=4.E-22, ALFX<0.001;

CONSTRAINT,SEQUENCE=LHCB1,RANGE=MQ3,BETX>100, BETY=4.E-22, DX>0.5850972729404622;

CONSTRAINT,SEQUENCE=LHCB1,RANGE=IP2/IP3,BETX>100, BETY=4.E-22, DX>0.5850972729404622;

CONSTRAINT,SEQUENCE=LHCB1,RANGE=IP2/IP3,BETX>100, BETY=4.E-22, DX>0.5850972729404622;

Prototype graphics package

Prototype package for tracing MAD input

Settings and options

Created by Wolfram Mathematica 6.0  (30 August 2007)