예제 #1
0
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
예제 #2
0
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
예제 #3
0
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