folder = "tmp/output_++--_FC=40.0_p=200.0_10"
emittance_trans = 9
emittance_long = 0.0
n_spills = 1000
n_per_spill = 10
momentum = 200.0
physics = "standard"
energy = (momentum**2.+105.658**2)**0.5

if emittance_trans < 1e-9:
    transverse = {"transverse_mode":"pencil"}
else:
    transverse = {
          "transverse_mode":"constant_solenoid",
          "emittance_4d":emittance_trans,
          "normalised_angular_momentum":0.0,
          "bz":4.e-3,
    }

if emittance_long < 1e-6:
    longitudinal = {
            "longitudinal_mode":"pencil",
            "momentum_variable":"p",
    }
else:
    longitudinal = {
            "longitudinal_mode":"twiss",
            "momentum_variable":"p",
            "beta_l":10.,
            "alpha_l":-1.,
            "emittance_l":emittance_long,
    }

simulation_reference_particle = {
    "random_seed": 0,
    "energy":226.,
    "particle_id":-13,
    "time": 0.0,
    "position":{"x":0.0, "y":0.0, "z":-4050.001},
    "momentum":{"x":0.0, "y":0.0, "z":1.0}
}
simulation_geometry_filename = "StepPi.dat"
verbose_level = 1

spill_generator_number_of_spills = n_spills

beam = {
    "particle_generator":"counter",
    "random_seed":0,
    "definitions":[{
       "reference":simulation_reference_particle,
       "random_seed_algorithm":"incrementing_random",
       "n_particles_per_spill":n_per_spill,
       "transverse":transverse,
       "longitudinal":longitudinal,
       "coupling":{"coupling_mode":"none"}
       #"spin":{"x":0.0, "y":0.0, "z":1.0}
    },]
}
physics_processes = physics
particle_decay = False

output_root_file_name = "maus_output_et="+str(emittance_trans)+"_el="+str(emittance_long)+"_"+str(n_spills*n_per_spill)+"_"+physics+".root"