~max-mcginley/third-order-lie-map/generator

 

static const double stepSizePlot = 10.0;

 

static const char* lattice_locale = "lattices";

 

static std::string lattice_filepath(const char* name) {return std::string(lattice_locale) + std::string("/") + std::string(name) + std::string("/StepPi.") + std::string(name) + std::string(".dat");}

static const int numberOfParticles = 96;

 

static const double stepSize = 2.0;

static const double posScale = 0.001;

static const double momScale = 0.001;

 

static const double lengthScale = 1.0/299.792; //One millimeter in terms of base units ns ie 1/c in ns mm^-1

static const double lightSpeed = 299792458.;

static const double kiloTeslaScale = 89856.7; // One kilotesla. Energy is in MeV. Base units are 1ns and 1Mev. 1kT = 10e3 kg A^-1 s^-2 = kiloTeslaScale * MeV ns^-1

static const double sigma = 2000.; //a convenient length scale, in mm

static const double muonMass = 105.65837;

static const double refMomentum = 200.;

static const double refEnergy = std::sqrt(refMomentum*refMomentum + muonMass*muonMass);

static const double stepSizeFieldPlot = 0.03125;

static const double zRoundingValueB = 64.0; //Rounds z to the nearest 1/<value>. Ensure <val> > 1.0/stepSizeFieldPlot   

static const double reducedStep = 0.25;

static const double zRoundingValue2 = 8.0;

static const double peakB = 0.002;

static const double bumpLength = 3000.0;

static const int numberOfEmittanceParticles = 10000;

static const double defaultPosMomRation = std::sqrt(sigma/refMomentum);

static const int inOutValues = 40;

static const double maxEmittanceIn = 1200.0;

static const int hamiltonianCoefficientCodes[19] = {256, 1152, 2048, 641, 200, 1537, 1096, 130, 1040, 1026, 144, 585, 515, 74, 529, 88, 4, 18, 32};

 

static bool streaming = true;

 

static void turnOffStream() {streaming = false;}

static void turnOnStream() {streaming = true;}

 

static std::ostream& getCurrentStream() {

        //if (streaming) return std::cout;

        return std::cerr;

}