def process(dirname, basename, recalc=False, plot=False): inst = basename.split('-')[0] st = basename.split('-')[1] pc = published_constants.PHYSICAL_CONSTANT_RANGES[inst] # find average tension, length T = ( pc[st]['T'][0] + pc[st]['T'][1] ) / 2 L = ( pc['L'][0] + pc['L'][1] ) / 2 #print T, L decayss, f0, B = get_data(dirname, basename) nf, nb, na, stats = mode_decays.handle_decays(decayss, basename, f0, B, T, L, None, plot=plot, recalc=recalc) return nf, nb, na, stats, None
def process(dirname, basename, recalc=False, plot=False): inst = basename.split('-')[0] st = basename.split('-')[1] pc = published_constants.PHYSICAL_CONSTANT_RANGES[inst] # find average tension, length T = (pc[st]['T'][0] + pc[st]['T'][1]) / 2 L = (pc['L'][0] + pc['L'][1]) / 2 #print T, L decayss, f0, B = get_data(dirname, basename) nf, nb, na, stats = mode_decays.handle_decays(decayss, basename, f0, B, T, L, None, plot=plot, recalc=recalc) return nf, nb, na, stats, None
def process(dirname, basename, recalc=False, plot=False): inst = basename.split('-')[0] st = basename.split('-')[1] pc = published_constants.PHYSICAL_CONSTANT_RANGES[inst] # find average tension, length T = (pc[st.upper()]['T'][0] + pc[st.upper()]['T'][1]) / 2 L = (pc['L'][0] + pc['L'][1]) / 2 #print T, L harms = harmonics_data.HarmonicsData(dirname, basename) decays, f0, B, harm_stats = harms.get_data() nf, nb, stats = mode_decays.handle_decays(decays, basename, f0, B, T, L, harm_stats, plot=plot, recalc=recalc) return nf, nb, stats, harm_stats