def main():
args = sys.argv[1:]
if len(args) < 7:
return usage()
dst_path = args[0]
sample = args[1]
decay = args[2]
mode = args[3]
type = args[4]
ecms = args[5]
dst_num = args[6]
sys.stdout.write('Scanning %s...\n' % dst_path)
dst_list = []
print '***************************************start to search***************************************'
dst_list = search(dst_list, dst_path, '.dst')
print '***************************************searching ending**************************************'
dst_list_groups = group_files_by_num(dst_list, dst_num)
i = 0
for dst_list_fill in dst_list_groups:
file_name = sample + '_' + decay + '_' + mode + '_' + ecms + '_' + str(
i) + '.txt'
f = open(file_name, 'w')
f.write('#include "$ROOTIOROOT/share/jobOptions_ReadRec.txt"\n')
f.write('#include "$MAGNETICFIELDROOT/share/MagneticField.txt"\n')
f.write('#include "$DTAGALGROOT/share/jobOptions_dTag.txt"\n')
f.write('#include "$DDECAYALGROOT/share/jobOptions_DDecay.txt"\n')
f.write('#include "$MEASUREDECMSSVCROOT/share/anaOptions.txt"\n')
f.write('DDecay.Ecms = ' + str(float(ecms) / 1000.) + ';\n')
f.write('DDecay.type = "' + type + '";\n')
f.write('\n')
f.write('// Input REC or DST file name\n')
f.write('EventCnvSvc.digiRootInputFile = {\n')
for dst in dst_list_fill:
if dst != dst_list_fill[-1]:
temp = '"' + dst + '",\n'
f.write(temp)
if dst == dst_list_fill[-1]:
temp = '"' + dst + '"\n'
f.write(temp)
f.write('};\n')
f.write('\n')
f.write(
'// Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )\n'
)
f.write('MessageSvc.OutputLevel = 6;\n')
f.write('\n')
f.write('// Number of events to be processed (default is 10)\n')
f.write('ApplicationMgr.EvtMax = -1;\n')
f.write('\n')
f.write('ApplicationMgr.HistogramPersistency = "ROOT";\n')
f.write(
'NTupleSvc.Output = {\"FILE1 DATAFILE=\'/besfs/groups/cal/dedx/$USER/bes/DstDpi/run/DstDpi/rootfile/'
+ sample + '/' + type + '/' + ecms + '/' + sample + '_' + decay +
'_' + mode + '_' + ecms + '_' + str(i) + '.root' +
'\' OPT=\'NEW\' TYP=\'ROOT\'\"};\n')
f.close()
i = i + 1
print 'All done!'
def main():
args = sys.argv[1:]
if len(args) < 5:
return usage()
dst_path = args[0]
runNo_low = args[1]
runNo_up = args[2]
ecms = args[3]
cms = args[4]
sys.stdout.write('Scanning %s...\n' % dst_path)
for runNo in range(int(runNo_low), int(runNo_up) + 1):
dst_list = []
print '***************************************start to search***************************************'
dst_list = search(dst_list, dst_path, '00' + str(runNo))
print '***************************************searching ending**************************************'
if len(dst_list) > 0 and len(dst_list) < 20:
write_file(ecms, cms, runNo, dst_list)
elif len(dst_list) >= 20:
dst_list_groups = group_files_by_num(dst_list, 20)
i = 0
for dst_list_fill in dst_list_groups:
write_file(ecms, cms, runNo, dst_list_fill, i)
i += 1
else:
print 'runNo: ' + str(runNo) + ' is empty, just ignore it!'
print 'All done!'
def main():
args = sys.argv[1:]
if len(args) < 8:
return usage()
dst_path = args[0]
sample = args[1]
decay = args[2]
mode = args[3]
type = args[4]
ecms = args[5]
cms = args[6]
dst_num = args[7]
sys.stdout.write('Scanning %s...\n' % dst_path)
dst_list = []
print '***************************************start to search***************************************'
dst_list = search(dst_list, dst_path, '.dst')
print '***************************************searching ending**************************************'
dst_list_groups = group_files_by_num(dst_list, dst_num)
i = 0
for dst_list_fill in dst_list_groups:
file_name = sample + '_' + decay + '_' + mode + '_' + ecms + '_' + str(
i) + '.txt'
f = open(file_name, 'w')
f.write('#include "$ROOTIOROOT/share/jobOptions_ReadRec.txt"\n')
f.write('#include "$MAGNETICFIELDROOT/share/MagneticField.txt"\n')
f.write('#include "$DTAGALGROOT/share/jobOptions_dTag.txt"\n')
f.write('#include "$DDECAYALGROOT/share/jobOptions_DDecay.txt"\n')
f.write('#include "$MEASUREDECMSSVCROOT/share/anaOptions.txt"\n')
f.write('DDecay.IsMonteCarlo = true;\n')
f.write('DDecay.Ecms = ' + str(float(cms) / 1000.) + ';\n')
f.write('DDecay.W_m_Kpipi = ' + str(width(ecms[0:4])) + ';\n')
f.write('DDecay.W_rm_Dpipi = ' + str(window(ecms[0:4])) + ';\n')
f.write('\n')
f.write('DTag.NeutralDReconstruction = true;\n')
f.write('DTag.ChargedDReconstruction = true;\n')
f.write('DTag.DsReconstruction = true;\n')
f.write('\n')
f.write('NeutralDSelector.UseMbcCuts = false;\n')
f.write('eutralDSelector.UseDeltaECuts = false;\n')
f.write('NeutralDSelector.UseDeltaMassCuts = true;\n')
f.write('\n')
f.write('ChargedDSelector.UseMbcCuts = false;\n')
f.write('ChargedDSelector.UseDeltaECuts = false;\n')
f.write('ChargedDSelector.UseDeltaMassCuts = true;\n')
f.write('\n')
f.write('DsSelector.UseMbcCuts = false;\n')
f.write('DsSelector.UseDeltaECuts = false;\n')
f.write('DsSelector.UseDeltaMassCuts = true;\n')
f.write('\n')
f.write('// Input REC or DST file name\n')
f.write('EventCnvSvc.digiRootInputFile = {\n')
for dst in dst_list_fill:
if dst != dst_list_fill[-1]:
temp = '"' + dst + '",\n'
f.write(temp)
if dst == dst_list_fill[-1]:
temp = '"' + dst + '"\n'
f.write(temp)
f.write('};\n')
f.write('\n')
f.write(
'// Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )\n'
)
f.write('MessageSvc.OutputLevel =6;\n')
f.write('\n')
f.write('// Number of events to be processed (default is 10)\n')
f.write('ApplicationMgr.EvtMax = -1;\n')
f.write('\n')
f.write('ApplicationMgr.HistogramPersistency = "ROOT";\n')
f.write(
'NTupleSvc.Output = {\"FILE1 DATAFILE=\'/scratchfs/bes/$USER/bes/DDPIPI/v0.2/'
+ sample + '/' + type + '/' + ecms +
'/sys_err/psipp_shape/rootfile/' + sample + '_' + decay + '_' +
mode + '_' + ecms + '_' + str(i) + '.root' +
'\' OPT=\'NEW\' TYP=\'ROOT\'\"};\n')
f.close()
i = i + 1
print 'All done!'
def main():
args = sys.argv[1:]
if len(args) < 7:
return usage()
dst_path = args[0]
sample = args[1]
decay = args[2]
mode = args[3]
type = args[4]
ecms = args[5]
dst_num = args[6]
sys.stdout.write('Scanning %s...\n' % dst_path)
dst_list = []
print '***************************************start to search***************************************'
dst_list = search(dst_list, dst_path, '.dst')
print '***************************************searching ending**************************************'
dst_list_groups = group_files_by_num(dst_list, dst_num)
i = 0
for dst_list_fill in dst_list_groups:
file_name = sample + '_' + decay + '_' + mode + '_' + ecms + '_' + str(
i) + '.txt'
f = open(file_name, 'w')
f.write('#include "$ROOTIOROOT/share/jobOptions_ReadRec.txt"\n')
f.write('#include "$MAGNETICFIELDROOT/share/MagneticField.txt"\n')
f.write('#include "$DTAGALGROOT/share/jobOptions_dTag.txt"\n')
f.write(
'#include "$SIMPLEPIDSVCROOT/share/jobOptions_SimplePIDSvc.txt"\n')
f.write('#include "$VERTEXFITROOT/share/jobOptions_VertexDbSvc.txt"\n')
f.write('#include "$MEASUREDECMSSVCROOT/share/anaOptions.txt"\n')
f.write('\n')
f.write('NeutralDReconstruction.ReadBeamEFromDB = false;\n')
f.write('NeutralDReconstruction.UseCalibBeamE = false;\n')
f.write('NeutralDReconstruction.BeamE = 2.09;\n')
f.write('NeutralDSelector.UseMbcCuts = false;\n')
f.write('NeutralDSelector.UseDeltaECuts = false;\n')
f.write('NeutralDSelector.UseDeltaMassCuts = true;\n')
f.write('NeutralDSelector.DDeltaMassMinCut = -0.12;\n')
f.write('NeutralDSelector.DDeltaMassMaxCut = 0.12;\n')
f.write('ChargedDReconstruction.ReadBeamEFromDB = false;\n')
f.write('ChargedDReconstruction.UseCalibBeamE = false;\n')
f.write('ChargedDSelector.UseMbcCuts = false;\n')
f.write('ChargedDSelector.UseDeltaECuts = false;\n')
f.write('ChargedDSelector.UseDeltaMassCuts = true;\n')
f.write('ChargedDSelector.DDeltaMassMinCut = -0.12;\n')
f.write('ChargedDSelector.DDeltaMassMaxCut = 0.12;\n')
f.write('LocalKaonSelector.useSimplePID = false;\n')
f.write('LocalPionSelector.useSimplePID = false;\n')
f.write('\n')
f.write('ApplicationMgr.DLLs += {"PiD0DmAlg"};\n')
f.write('ApplicationMgr.TopAlg +={ "PiD0Dm" };\n')
f.write('PiD0Dm.AddModesD0 = {0, 1, 3, 4};\n')
f.write('PiD0Dm.AddModesDm = {200, 201, 202, 203, 204};\n')
f.write('PiD0Dm.sample_types = {"signal"};\n')
f.write('PiD0Dm.Ecms = ' + str(float(ecms) / 1000.) + ';\n')
f.write('PiD0Dm.W_D0 = 0.024;\n')
f.write('PiD0Dm.W_Dm = 0.020;\n')
f.write('\n')
f.write('// Input REC or DST file name\n')
f.write('EventCnvSvc.digiRootInputFile = {\n')
for dst in dst_list_fill:
if dst != dst_list_fill[-1]:
temp = '"' + dst + '",\n'
f.write(temp)
if dst == dst_list_fill[-1]:
temp = '"' + dst + '"\n'
f.write(temp)
f.write('};\n')
f.write('\n')
f.write(
'// Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )\n'
)
f.write('MessageSvc.OutputLevel = 6;\n')
f.write('\n')
f.write('// Number of events to be processed (default is 10)\n')
f.write('ApplicationMgr.EvtMax = -1;\n')
f.write('\n')
f.write('ApplicationMgr.HistogramPersistency = "ROOT";\n')
f.write(
'NTupleSvc.Output = {\"FILE1 DATAFILE=\'/besfs/groups/cal/dedx/$USER/bes/DDbarPi-DT/run/DDbarPi/rootfile/'
+ sample + '/' + type + '/' + ecms + '/' + sample + '_' + decay +
'_' + mode + '_' + ecms + '_' + str(i) + '.root' +
'\' OPT=\'NEW\' TYP=\'ROOT\'\"};\n')
f.close()
i = i + 1
print 'All done!'