def ratios_signal(args):
if args[0] == '537ipb/281ipb' and args[1] == '537ipbv2/281ipb':
tabfilename1 = 'compare_yields_signal_divide_537ipb_9.03_regular12.txt'
title1 = args[0]
headname1 = title1
tabfilename2 = 'compare_yields_signal_divide_537ipbv2_9.03_regular12.txt'
title2 = args[1]
headname2 = title2
else:
raise NameError(args)
tabfile1 = os.path.join(attr.tabpath, tabfilename1)
tabfile2 = os.path.join(attr.tabpath, tabfilename2)
if opts.set and opts.set == 'info':
print 'here'
tab = DHadTable()
tab.column_append_from_tab_file('Mode', tabfile1, 'Mode')
tab.column_append_from_tab_file(title1, tabfile1, headname1)
tab.column_append_from_tab_file(title2, tabfile2, headname2)
tab.column_append_by_diff_sigma_pct('diff(%)', title2, title1)
tab.output(_tabname)
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)
def evt_log(args):
print args
if args[0] not in ['data'] :
raise NameError(args)
dt_type = args[0]
tag = args[1]
label = args[2]
evtpath = attr.evtpath()
if '/' in label:
ver = label.split('/')[0]
label = label.split('/')[1]
evtpath = attr.evtpath(ver)
prefix = 'dir_' + label
evttab = 'entries_evt' + '_'.join(args).replace('/', '_') + '.tab'
tab = DHadTable()
tab.column_append_from_tab_file('Mode', evttab)
tab.output()
sys.exit()
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)
def evt_log(args):
print args
if args[0] not in ['data']:
raise NameError(args)
dt_type = args[0]
tag = args[1]
label = args[2]
evtpath = attr.evtpath()
if '/' in label:
ver = label.split('/')[0]
label = label.split('/')[1]
evtpath = attr.evtpath(ver)
prefix = 'dir_' + label
evttab = 'entries_evt' + '_'.join(args).replace('/', '_') + '.tab'
tab = DHadTable()
tab.column_append_from_tab_file('Mode', evttab)
tab.output()
sys.exit()
def compare_yields_data_signal_divide_537ipb(args):
sys.stdout.write('dhad.table: Creating %s ...\n' % _tabname)
tabpath1 = os.path.join(
attr.base, '7.06', 'tab')
tabfile1 = os.path.join(
tabpath1, 'compare_yields_data_divide_281ipb_537ipb.txt')
tabfile2 = os.path.join(
attr.tabpath, 'compare_yields_signal_divide_537ipb_9.03_regular12.txt')
tab = DHadTable()
tab.column_append_from_tab_file('Mode', tabfile1, 'Mode')
tab.column_append_from_tab_file('Data', tabfile1, '537ipb/281ipb')
tab.column_append_from_tab_file('Signal', tabfile2, '537ipb/281ipb')
tab.column_append_by_diff_sigma_pct('diff(%)', 'Signal', 'Data')
tab.output(_tabname)
def signal_line_shape_syst(args):
'''
--------------------------------------------------
Signal line shape systematics
--------------------------------------------------
'''
tab_name = 'signal_line_shape_syst'
tabbase = os.path.join(attr.base, '7.06', 'tab')
tab = DHadTable()
tab.column_append_from_dict('Mode', 'fname,fnamebar')
tab_file_A = os.path.join(
tabbase, 'yields_M_3.7713_regular_M_3.7723_data_Single.txt')
name_A = 'Mass(+/-0.5 MeV)%'
column = 'max-diff(%)'
tab.column_append_from_tab_file(name_A, tab_file_A, column)
tab_file_B = os.path.join(
tabbase, 'yields_Gamma_0.0261_regular_Gamma_0.0311_data_Single.txt')
name_B = 'Gamma(+/-2.5 MeV)%'
tab.column_append_from_tab_file(name_B, tab_file_B, column)
tab_file_C = os.path.join(
tabbase, 'yields_R_8.3_regular_R_16.3_data_Single.txt')
name_C = 'R(+/- 4)%'
tab.column_append_from_tab_file(name_C, tab_file_C, column)
name_D = 'Total(%)'
tab.column_append_by_add_quadrature3(name_D, name_A,
name_B, name_C, rnd='.01')
tab.rows_join_by_max(name_D)
tab.output(tab_name)
