if (args.decond_D is None): decond_D = args.decond else: decond_D = args.decond_D if (args.decond_ecdec is None): decond_ecdec = args.decond else: decond_ecdec = args.decond_ecdec edf, _, eBins = da.get_edf(args.decond)[0:3] DI, _, _, fit = da.get_diffusion(decond_D)[0:4] edD, _, _, eBins_edD = da.get_decD(decond_D, da.DecType.energy)[0:4] edf_edD = da.get_edf(decond_D)[0] sig_I, _, sig_IL, _, eBins_sig = da.get_ec_dec_energy(decond_ecdec, sep_nonlocal=True, threshold=0)[0:5] eBins /= da.const.calorie eBins_edD /= da.const.calorie eBins_sig /= da.const.calorie edf *= da.const.angstrom**3 * da.const.calorie edf_edD *= da.const.angstrom**3 * da.const.calorie DI /= da.const.angstrom**2 / da.const.pico edD /= da.const.angstrom**2 / da.const.pico numPlots = 3 fig, axs = plt.subplots(numPlots, 1, sharex=False, figsize=figsize3) abcPos = (0.03, 0.965)
if (args.decond_D is None): decond_D = args.decond else: decond_D = args.decond_D if (args.decond_ecdec is None): decond_ecdec = args.decond else: decond_ecdec = args.decond_ecdec edf, _, eBins = da.get_edf(args.decond)[0:3] DI, _, _, fit = da.get_diffusion(decond_D)[0:4] edD, _, _, eBins_edD = da.get_decD(decond_D, da.DecType.energy)[0:4] edf_edD = da.get_edf(decond_D)[0] sig_IL, _, eBins_sig = da.get_ec_dec_energy(decond_ecdec)[0:3] eBins /= da.const.calorie eBins_edD /= da.const.calorie eBins_sig /= da.const.calorie edf *= da.const.angstrom**3 * da.const.calorie edf_edD *= da.const.angstrom**3 * da.const.calorie DI /= da.const.angstrom**2 / da.const.pico edD /= da.const.angstrom**2 / da.const.pico numPlots = 3 fig, axs = plt.subplots(numPlots, 1, sharex=False, figsize=figsize3) abcPos = (0.03, 0.965)