def crossSections(opts, tabname, label):
bffilename = 'bf_stat_sys'
bffile = tools.set_file(
extbase=attr.brfpath, prefix='dir_'+label, comname=bffilename)
brf = BrfFile(bffile)
tab = THSTable()
tab.row_append(['Quantity', 'Value'])
names = ['sigma(D0D0bar)', 'sigma(D+D-)',
'sigma(DDbar)', 'chg/neu']
for name in names:
row = [name, brf.parsed[name]]
if name == 'chg/neu':
tab.row_append(row, rnd='.001')
else:
name = row[0]
value = row[1]
value = tab.cell_trim(
cell=value, factor=0.001, rnd='.001',
opt='(cell){~\\rm nb}')
row = [name, value]
tab.row_append(row)
tab.output(tabname, trans_dict=attr.cross_sections_dict)
sys.stdout.write('%s\n' % brf.parsed['coeff_ddbar'])
def mode_dependent_syst(opts, tabname):
tab = THSTable()
tabprefix = 'dir_818ipbv12'
modenames = [modes[mode]['sname'] for mode in modes]
modenames.insert(0, 'Source')
tab.row_append(modenames)
tab.row_append(cbx.mode_syst_signal_shape(tabprefix))
tab.row_append(cbx.mode_syst_tracking(tabprefix))
tab.row_append(cbx.mode_syst_ks(tabprefix), 'Eff($\KS$)')
tab.row_append(cbx.mode_syst_pi0(tabprefix), 'Eff($\pi^0$)')
tab.row_append(cbx.mode_syst_pid_pion(tabprefix), 'PID - $\pipm$')
tab.row_append(cbx.mode_syst_pid_kaon(tabprefix), 'PID - $\Kpm$ ')
tab.row_append(cbx.mode_syst_lepton_veto(tabprefix))
tab.row_append(cbx.mode_syst_de(tabprefix), '$|\DeltaE|$ (*)')
tab.row_append(cbx.mode_syst_bkg_shape(tabprefix), 'Bkgd shape')
tab.row_append(cbx.mode_syst_fsr(tabprefix))
tab.row_append(cbx.mode_syst_substructure(tabprefix))
tab.row_append(cbx.mode_syst_multcand(tabprefix))
texhead = r'''Source & $K\pi$ & $K\pi\pi^0$ & $K\pi\pi\pi$ & $K\pi\pi$ & $K\pi\pi\pi^0$ & $K_S^0\pi$ & $K_S^0\pi\pi^0$ & $K_S^0\pi\pi\pi$ & $KK\pi$ '''
tab.output(tabname, texhead)