forked from PyCOMPLETE/PyPARIS_sim_class
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Simulation_parameters.py.example
147 lines (114 loc) · 3.26 KB
/
Simulation_parameters.py.example
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
from scipy.constants import c
####################
# 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 = 16
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
n_non_parallelizable = 2 #rf and aperture
# Transverse Damper Settings
enable_transverse_damper = False
dampingrate_x = 100.
dampingrate_y = 100.
if enable_transverse_damper: n_non_parallelizable += 1
###################
# Beam Parameters #
###################
bunch_from_file = None
intensity = 1.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
# Numerical Parameters
n_slices = 500
z_cut = 2.5e-9/2*c # For slicing
macroparticles_per_slice = 5000
n_macroparticles = macroparticles_per_slice*n_slices
#################
# Stop Criteria #
#################
# 1. Turns
N_turns = 128 # Per job
N_turns_target = 20000
# 2. Losses
sim_stop_frac = 0.9
# 3. Emittance Growth
flag_check_emittance_growth = True
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}
# 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 = True
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