pt.AddText('e^{+}e^{-}#rightarrowD^{0}#pi_{0}^{+}D^{*-}')
mbc.SaveAs('./figs/scatter_rm_D_m_Dpi_' + str(ecms) + '_' + D_sample +
'_' + D_type + '.pdf')
# raw_input('Press <Enter> to end...')
if __name__ == '__main__':
args = sys.argv[1:]
if len(args) < 3:
usage()
sys.exit()
ecms = int(args[0])
D_sample = args[1]
D_type = args[2]
path = []
path.append(
'/besfs/groups/cal/dedx/$USER/bes/DstDpi/run/DstDpi/anaroot/data/' +
str(ecms) + '/data_' + str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_signal_before.root')
leg_title = str(ecms) + ' MeV'
min = 3.7
max = 9.0
if D_sample == 'Dplus':
bins = 180
if D_sample == 'D0':
bins = 330
plot(path, leg_title, ecms, min, max, bins, D_sample, D_type)
if __name__ == '__main__':
args = sys.argv[1:]
if len(args) < 3:
usage()
sys.exit()
ecms = int(args[0])
D_sample = args[1]
D_type = args[2]
path = []
path.append(
'/besfs/groups/cal/dedx/$USER/bes/DstDpi/run/DstDpi/anaroot/data/' +
str(ecms) + '/data_' + str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_signal_before.root')
path.append(
'/besfs/groups/cal/dedx/$USER/bes/DstDpi/run/DstDpi/anaroot/data/' +
str(ecms) + '/data_' + str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_background_before.root')
path.append(
'/besfs/groups/cal/dedx/$USER/bes/DstDpi/run/DstDpi/anaroot/mc/D2_2460/'
+ str(ecms) + '/mc_D2_2460_' + str(ecms) + '_DstDpi_' + D_sample +
'_' + convert_name(D_sample, D_type) + '_signal_before.root')
leg_title = str(ecms) + ' MeV'
if D_sample == 'Dplus':
xmin = 2.0
xmax = 2.9
xbins = 180
if D_sample == 'D0':
xmin = 2.0
f_out = TFile(path[1], 'recreate')
h_mDpi.Write()
h_mDbarpi.Write()
h_mDDbar.Write()
h_dalitz.Write()
f_out.Close()
if __name__ == '__main__':
args = sys.argv[1:]
if len(args)<4:
usage()
sys.exit()
ecms = int(args[0])
D_sample = args[1]
D_type = args[2]
sample_type = args[3]
if not os.path.exists('./anaroot/'):
os.makedirs('./anaroot/')
path = []
if sample_type == 'data':
path.append('/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/data/data_'+str(ecms)+'_DstDpi_'+D_sample+'_'+convert_name(D_sample, D_type)+'_signal_angle_before.root')
path.append('./anaroot/PWA_dt.root')
if sample_type == 'background':
path.append('/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/data/data_'+str(ecms)+'_DstDpi_'+D_sample+'_'+convert_name(D_sample, D_type)+'_background_angle_before.root')
path.append('./anaroot/bg1.root')
if sample_type == 'mc':
path.append('/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/mc/mc_'+str(ecms)+'_DstDpi_'+D_sample+'_'+convert_name(D_sample, D_type)+'_angle_before.root')
path.append('./anaroot/PWA_mc.root')
convert(path, sample_type)
if len(args) < 4:
usage()
sys.exit()
ecms = int(args[0])
D_sample = args[1]
D_type = args[2]
sample_type = args[3]
if not os.path.exists('../data/'):
os.makedirs('../data/')
path = []
if sample_type == 'data':
path.append(
'/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/data/data_'
+ str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_signal_angle_before.root')
path.append('../data/PWA_dt.dat')
if sample_type == 'background':
path.append(
'/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/data/data_'
+ str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_background_angle_before.root')
path.append('../data/bg1.dat')
if sample_type == 'mc':
path.append(
'/besfs/groups/cal/dedx/jingmq/bes/DstDpi_momentum/python/anaroot/mc/mc_'
+ str(ecms) + '_DstDpi_' + D_sample + '_' +
convert_name(D_sample, D_type) + '_angle_before.root')
path.append('../data/PWA_mc.dat')
extract(path)