def make_it_easier(*args):
furman_pivi_surface_LHC = {
'use_modified_sigmaE': False,
'use_ECLOUD_theta0_dependence': False,
'use_ECLOUD_energy': False,
'conserve_energy': False,
'exclude_rediffused': True,
'choice': 'poisson',
'M_cut': M_cut,
'p_n': np.array(args[0:10]),
# 'eps_n': np.array(args[10:]),
'eps_n': 1.634185 / (np.array(args[0:10]) - 1),
# Parameters for backscattered electrons
'p1EInf': 0.02,
'p1Ehat': 0.496,
'eEHat': 0.,
'w': 60.86,
'p': 1.,
'e1': 0.26,
'e2': 2.,
'sigmaE': 2.,
# Parameters for rediffused electrons
'p1RInf': 0.2,
'eR': 0.041,
'r': 0.104,
'q': 0.5,
'r1': 0.26,
'r2': 2.,
# Parameters for true secondaries
'deltaTSHat': 1.8848,
'eHat0': 332.,
's': 1.35,
't1': 0.5, # t1 and t2 based on taylor expansion
't2': 1., # of PyECLOUD formula for E_max(theta)
't3': 0.7,
't4': 1.
}
sey_mod_FP = fp.SEY_model_furman_pivi(
furman_pivi_surface_LHC,
E_th=None,
sigmafit=None,
mufit=None,
switch_no_increase_energy=0,
thresh_low_energy=None,
secondary_angle_distribution='cosine_3D')
chamb = ellip_cham_geom_object(1., 1., flag_verbose_file=False)
impact_management_object = impact_management(chamb=chamb,
sey_mod=sey_mod_FP,
Dx_hist=.1,
scrub_en_th=25.,
Nbin_En_hist=100,
En_hist_max=3000,
flag_seg=False,
cos_angle_width=0.05,
flag_cos_angle_hist=True)
return sey_mod_FP, impact_management_object, furman_pivi_surface_LHC
'q': 0.5,
'r1': 0.26,
'r2': 2.,
'deltaTSHat': 1.8848,
'eHat0': 276.8,
's': 1.54,
't1': 0.66,
't2': 0.8,
't3': 0.7,
't4': 1.
}
sey_mod = fp.SEY_model_furman_pivi(E_th=35.,
sigmafit=1.0828,
mufit=1.6636,
secondary_angle_distribution='cosine_3D',
switch_no_increase_energy=0,
thresh_low_energy=-1,
furman_pivi_surface=furman_pivi_surface_LHC)
def extract_sey_curves(n_rep, E_impact_eV_test, cos_theta_test, charge, mass):
deltas = {}
for etype in list(sey_mod.event_types.keys()):
etype_name = sey_mod.event_types[etype]
deltas[etype_name] = np.zeros(
(len(cos_theta_test), len(E_impact_eV_test)))
print('Extracting SEY curves...')
for i_ct, ct in enumerate(cos_theta_test):
print(('%d/%d' % (i_ct + 1, len(cos_theta_test))))
'eHat0': 276.8,
's': 1.54,
't1': 0.66,
't2': 0.8,
't3': 0.7,
't4': 1.
}
flag_costheta_Emax_shift = True
flag_costheta_delta_scale = True
test_obj = fp.SEY_model_furman_pivi(
E_th=35.,
sigmafit=1.0828,
mufit=1.6636,
secondary_angle_distribution='cosine_3D',
switch_no_increase_energy=0,
thresh_low_energy=-1,
furman_pivi_surface=furman_pivi_surface_scaled,
flag_costheta_delta_scale=flag_costheta_delta_scale,
flag_costheta_Emax_shift=flag_costheta_Emax_shift)
qq = 0.5 # From FP paper
sigma_e = 2.
E_0_single = 1600.
E_0 = np.array([E_0_single] * int(1e5))
energy = np.linspace(0.001, E_0_single, num=int(1e5))
alpha = 0.3
round_to_digits = 4
plt.close('all')
'deltaTSHat': 1.8848,
'eHat0': 332.,
's': 1.35,
't1': 0.5, # t1 and t2 based on taylor expansion
't2': 1., # of PyECLOUD formula for E_max(theta)
't3': 0.7,
't4': 1.
}
flag_costheta_Emax_shift = True
flag_costheta_delta_scale = True
sey_mod_initial = fp.SEY_model_furman_pivi(
furman_pivi_surface_LHC,
E_th=None,
sigmafit=None,
mufit=None,
switch_no_increase_energy=0,
thresh_low_energy=None,
secondary_angle_distribution='cosine_3D',
flag_costheta_Emax_shift=flag_costheta_Emax_shift,
flag_costheta_delta_scale=flag_costheta_delta_scale)
def make_it_easier(*args):
furman_pivi_surface_LHC = {
'use_modified_sigmaE': False,
'use_ECLOUD_theta0_dependence': False,
'use_ECLOUD_energy': False,
'conserve_energy': False,
'exclude_rediffused': True,
'choice': 'poisson',
'M_cut': M_cut,