knob_settings['par_crab1'] = 0. knob_settings['par_crab5'] = 0. ######################## # Build MAD-X instance # ######################## # Start mad Madx = pm.Madxp mad = Madx(command_log="mad_collider.log") # Set verbose flag mad.globals.par_verbose = int(configuration['verbose_mad_parts']) # Build sequence (alse creates link to optics_toolkit and calls it) ost.build_sequence(mad, beam=beam_to_configure) # Set twiss formats for MAD-X parts (macro from opt. toolkit) mad.input('exec, twiss_opt;') # Apply optics ost.apply_optics(mad, optics_file=optics_file) # Modify weak or strong beam crabbing mad.globals.ahcrab_l1b1 -= knob_settings[ 'par_crab_L1B1'] / 100. * mad.globals.ahcrab_l1b1 mad.globals.ahcrab_l1b2 -= knob_settings[ 'par_crab_L1B2'] / 100. * mad.globals.ahcrab_l1b2 mad.globals.ahcrab_r1b1 -= knob_settings[ 'par_crab_R1B1'] / 100. * mad.globals.ahcrab_r1b1 mad.globals.ahcrab_r1b2 -= knob_settings[
knob_settings['par_crab1'] = 0. knob_settings['par_crab5'] = 0. ######################## # Build MAD-X instance # ######################## # Start mad Madx = pm.Madxp mad = Madx(command_log="mad_collider.log") # Set verbose flag mad.globals.par_verbose = int(configuration['verbose_mad_parts']) # Build sequence (alse creates link to optics_toolkit and calls it) ost.build_sequence(mad, beam=beam_to_configure, configuration=configuration) # Set twiss formats for MAD-X parts (macro from opt. toolkit) mad.input('exec, twiss_opt;') # Apply optics ost.apply_optics(mad, optics_file=optics_file) # Attach beam to sequences mad.globals.nrj = configuration['beam_energy_tot'] particle_type = 'proton' if 'particle_mass' in configuration.keys(): particle_mass = configuration['particle_mass'] particle_type = 'ion' else: