def run(delta_pdgs, is_not, output_name, ccnc = 0):
if os.path.exists(output_name):
shutil.rmtree(output_name)
os.mkdir(output_name)
file_increment = 10
increment_counter = 0
file_counter = 1
number_of_arguments = len(sys.argv) - 1
loop_range = 1
if(number_of_arguments > file_increment):
loop_range = number_of_arguments / file_increment
if(number_of_arguments % file_increment):
loop_range += 1
my_proc = fmwk.ana_processor()
for i in xrange(0, loop_range):
my_proc.reset()
increment_counter += file_increment
while(file_counter <= increment_counter):
if file_counter >= len(sys.argv):
break
my_proc.add_input_file(sys.argv[file_counter])
file_counter += 1
my_proc.set_io_mode(fmwk.storage_manager.kBOTH)
my_proc.set_output_file(output_name + "/" + output_name + "_" + str(i) + ".root")
my_proc.enable_filter()
if(ccnc == 1):
my_proc.add_process(fmwk.ccnc(True))
elif(ccnc == 2):
my_proc.add_process(fmwk.ccnc(False))
idrs = fmwk.is_delta_rad_spec(is_not)
for i in xrange(0, len(delta_pdgs)):
idrs.add_delta_parents(delta_pdgs[i])
my_proc.add_process(idrs)
my_proc.run()
is_not = True
#ALL
################################################################################
my_proc = fmwk.ana_processor()
for i in xrange(1,len(sys.argv)):
my_proc.add_input_file(sys.argv[i])
my_proc.set_io_mode(fmwk.storage_manager.kBOTH)
out_name = "not_delta_rad_1delta_1gamma.root"
my_proc.set_output_file(out_name)
my_proc.enable_filter()
idrs = fmwk.is_delta_rad_spec(is_not)
idrs.add_delta_parents(2214)
idrs.add_delta_parents(2114)
my_proc.add_process(idrs)
my_proc.run()
################################################################################
my_proc2 = fmwk.ana_processor()
for i in xrange(1,len(sys.argv)):
my_proc2.add_input_file(sys.argv[i])
my_proc2.set_io_mode(fmwk.storage_manager.kBOTH)
out_name2 = "not_delta_rad_1deltaplus_1gamma.root"
while(file_counter <= increment_counter):
if(file_counter >= len(sys.argv)):
break
my_proc.add_input_file(sys.argv[file_counter])
file_counter += 1
my_proc.set_io_mode(fmwk.storage_manager.kBOTH)
my_proc.set_output_file(output_name + "/" + output_name + "_" + str(i) + ".root")
my_proc.enable_filter()
if(ccnc == 1):
my_proc.add_process(fmwk.ccnc(True))
elif(ccnc == 2):
my_proc.add_process(fmwk.ccnc(False))
idrs = fmwk.is_delta_rad_spec(is_not)
for(i in xrange(0, len(delta_pdgs))):
idrs.add_delta_parents(delta_pdgs[i])
my_proc.add_process(idrs)
my_proc.run()
run([2214, 2114], False, "delta_rad_1delta_1gamma")