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,