Forgot file

Author: giadarol

examples/006_impedance_instability/Simulation_parameters.py

@@ -0,0 +1,161 @@
+from scipy.constants import c
+
+
+write_buffer_every = 20
+
+
+####################
+# Machine Settings #
+####################
+
+machine_configuration = 'LHC-collision'
+
+# # Use this part for optics from file
+# # n_segments needs to be None if optics_pickle_file is specified
+# optics_pickle_file = 'lhc2018_25cm_only_triplets_IR15_b1_optics.pkl'
+# n_segments = None
+# beta_x =  None
+# beta_y =  None
+# Q_x = None
+# Q_y = None
+
+# # Use this part for smooth machine
+optics_pickle_file = None
+n_segments = 1
+beta_x =  400.0
+beta_y =  400.0
+Q_x = 62.27
+Q_y = 60.295
+
+Qp_x = 0.
+Qp_y = 0.
+
+octupole_knob = 0.
+
+V_RF = 12e6
+
+
+# Transverse Damper Settings
+enable_transverse_damper = False
+dampingrate_x = 50.
+dampingrate_y = 100.
+
+
+###################
+# Beam Parameters #
+###################
+
+bunch_from_file = None
+
+intensity = 2e+11
+
+epsn_x = 2.5e-6
+epsn_y = 2.5e-6
+
+sigma_z = 1.2e-9/4*c
+
+x_kick_in_sigmas = 0.1
+y_kick_in_sigmas = 0.1
+recenter_all_slices = True
+
+# Numerical Parameters
+n_slices = 500
+z_cut = 2.5e-9/2*c # For slicing
+macroparticles_per_slice = 100
+n_macroparticles = macroparticles_per_slice*n_slices
+
+
+#################
+# Stop Criteria #  
+#################
+
+# 1. Turns
+N_turns = 1000 # Per job
+N_turns_target = 1000
+# 2. Losses
+sim_stop_frac = 0.9
+# 3. Emittance Growth
+flag_check_emittance_growth = False
+epsn_x_max_growth_fraction = 0.5
+epsn_y_max_growth_fraction = epsn_x_max_growth_fraction
+
+
+######################
+# Footprint Settings #
+######################
+
+footprint_mode = False
+n_macroparticles_for_footprint_map = 500000
+n_macroparticles_for_footprint_track = 5000
+
+
+####################
+# E-Cloud Settings #
+####################
+
+# General E-Cloud Settings
+chamb_type = 'polyg'
+x_aper = 2.300000e-02
+y_aper = 1.800000e-02
+filename_chm = 'LHC_chm_ver.mat'
+Dt_ref = 5e-12
+pyecl_input_folder = './pyecloud_config'
+sey = 1.30
+
+# Transverse Multigrid Parameters
+PyPICmode = 'ShortleyWeller_WithTelescopicGrids'
+N_min_Dh_main = 10.
+Dh_sc_ext = .8e-3
+f_telescope = 0.3
+N_nodes_discard = 5.
+target_size_internal_grid_sigma = 10.
+target_Dh_internal_grid_sigma = 0.2
+custom_target_grid_arcs = None
+
+# # Uncomment for custom grid
+# custom_target_grid_arcs = {
+#     'x_min_target': -3e-3,
+#     'x_max_target': 3e-3,
+#     'y_min_target': -3.1e-3,
+#     'y_max_target': 3.1e-3,
+#     'Dh_target': 7e-5}
+
+force_interp_at_substeps_interacting_slices = True
+
+# Enable Kicks Different Planes
+enable_kick_x = True
+enable_kick_y = True
+
+# Dedicated Dipole E-Cloud Settings
+enable_arc_dip = False
+fraction_device_dip = 0.65
+init_unif_edens_flag_dip = 1
+init_unif_edens_dip = 1.000000e+12
+N_MP_ele_init_dip = 500000
+N_mp_max_dip = N_MP_ele_init_dip*4
+B_multip_dip = [8.33] #T
+
+# Dedicated Quadrupole E-Cloud Settings
+enable_arc_quad = False
+fraction_device_quad = 26.000000e-02 #7.000000e-02
+N_mp_max_quad = 2000000
+B_multip_quad = [0., 188.2] #T
+folder_path = '../../LHC_ecloud_distrib_quads/'
+filename_state = 'combined_distribution_sey%.2f_%.1fe11ppb_7tev.mat'%(sey,intensity/1e11)
+filename_init_MP_state_quad = folder_path + filename_state
+
+# Dedicated Kick Element Settings
+enable_eclouds_at_kick_elements = False
+path_buildup_simulations_kick_elements = '/home/kparasch/workspace/Triplets/ec_headtail_triplets/simulations_PyECLOUD/!!!NAME!!!_sey1.35'
+name_MP_state_file_kick_elements = 'MP_state_9.mat'
+orbit_factor = 6.250000e-01
+
+######################
+# Impedance settings #
+######################
+
+enable_impedance = True
+resonator_R_shunt = 25e6
+resonator_frequency = 2e9
+resonator_Q = 1.
+n_slices_wake = 500