Beispiel #1
0
def para_md(args):

    para = 'md'
    namestyle = 'fname'
    rnd = '0.0000001'
    tag = 's'

    tab_A = args[0]
    tab_B = args[1]

    latel_A = tab_A
    latel_B = tab_B

    label = tools.parseopts_set(opts.set, 'label')

    if label:
        labels = label.split(',')
        label_A = labels[0]
        label_B = labels[1]
       
    tab = DHadTable()
    tab.column_append_from_dict('Mode', namestyle)
    tab.column_append_from_fit_files(label_A, tab_A, para, tag, rnd)
    tab.column_append_from_fit_files(label_B, tab_B, para, tag, rnd)
    tab.column_append_by_diff('diff[MeV]', label_B,label_A, rnd=rnd, factor=1000)

    tab.output(_tabname)
Beispiel #2
0
def diff(args):

    label= None
    if opts.set:
        for li in opts.set.split(':'):
            name = li.split('=')[0]
            value = li.split('=')[1]
            sys.stdout.write('Set %s = %s \n' % (name, value))
            if name == 'label':
                label = value

    if label:
        labels = label.split(',')
        label_A = labels[0]
        label_B = labels[1]

    tab_A = 'compare_%s.txt' % args[0]
    tab_B = 'compare_%s.txt' % args[1]

    tabfile_A = os.path.join(attr.tabpath, tab_A)
    tabfile_B = os.path.join(attr.tabpath, tab_B)

    tab = DHadTable()
    tab.column_append_from_tab_file('Mode', tabfile_A, 'Mode')
    tab.column_append_from_tab_file(label_A, tabfile_A, 'diff(%)')
    tab.column_append_from_tab_file(label_B, tabfile_B, 'diff(%)')
    tab.column_append_by_diff('diff', label_B,label_A)
       
    tab.output(_tabname)
Beispiel #3
0
def fsr_syst(args):
    '''
    --------------------
       FSR systematics
    --------------------
    '''
    tab_name = 'FSR_syst'

    tabbase = os.path.join(attr.base, '7.06', 'tab')
    
    tab = DHadTable()
    tab.column_append_from_dict('Mode', 'fname,fnamebar')
    tab_file = os.path.join(tabbase, 'noFSR_single_signal_eff.txt')
    columnA    = 'Efficiency without FSR(%)'
    tab.column_append_from_tab_file(columnA, tab_file)
    columnB     = 'Efficiency with FSR(%)'
    tab.column_append_from_tab_file(columnB, tab_file)
    tab.column_append_by_diff('diff(%)', columnA, columnB, factor=100)
    tab.column_append_by_times_number('diff X 30%', 'diff(%)', '0.3',
                                      rnd='.01')
    tab.rows_join_by_max('diff X 30%')
    tab.output(tab_name)