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()
#NC
is_nc = True
################################################################################
my_proc4 = fmwk.ana_processor()
for i in xrange(1,len(sys.argv)):
my_proc4.add_input_file(sys.argv[i])
my_proc4.set_io_mode(fmwk.storage_manager.kBOTH)
out_name4 = "not_nc_delta_rad_1delta_1gamma.root"
my_proc4.set_output_file(out_name4)
my_proc4.enable_filter()
my_proc4.add_process(fmwk.ccnc(is_nc))
idrs4 = fmwk.is_delta_rad_spec(is_not)
idrs4.add_delta_parents(2214)
idrs4.add_delta_parents(2114)
my_proc4.add_process(idrs4)
my_proc4.run()
################################################################################
my_proc5 = fmwk.ana_processor()
for i in xrange(1,len(sys.argv)):
my_proc5.add_input_file(sys.argv[i])
msg = '\n'
msg += "Usage 1: %s $INPUT_ROOT_FILE\n" % sys.argv[0]
msg += '\n'
sys.stderr.write(msg)
sys.exit(1)
from ROOT import larlite as fmwk
from ROOT import gRandom
my_proc = fmwk.ana_processor()
#for i in xrange(1,len(sys.argv)):
#my_proc.add_input_file(sys.argv[i])
with open(sys.argv[2]) as file:
for line in file:
line = line.strip('\n')
my_proc.add_input_file(line)
break
my_proc.set_io_mode(fmwk.storage_manager.kBOTH)
my_proc.set_output_file(sys.argv[1])
my_proc.enable_filter()
my_proc.add_process(fmwk.ccnc(True))
my_proc.add_process(fmwk.delta_rad_filter())
my_proc.run()
from ROOT import larlite as fmwk
from ROOT import ertool
from seltool.trackpidDef import GetTrackPidInstance
gROOT.SetBatch(True)
# Create ana_processor instance
my_proc = fmwk.ana_processor()
my_proc.enable_filter(True)
# Set input root file
for x in xrange(1, len(sys.argv)):
print sys.argv[x]
my_proc.add_input_file(sys.argv[x])
# Specify IO mode
my_proc.set_io_mode(fmwk.storage_manager.kREAD)
my_proc.add_process(fmwk.ccnc(1));
my_proc.add_process(fmwk.GetPOT());
my_proc.run()
# done!
print
print "Finished running ana_processor event loop!"
print
sys.exit(0)
my_proc = fmwk.ana_processor()
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")
msg += '\n'
sys.stderr.write(msg)
sys.exit(1)
from ROOT import larlite as fmwk
from ROOT import gRandom
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)
file_name = "ccnc.root"
my_proc.set_output_file(file_name)
my_proc.enable_filter()
is_nc = False
my_proc.add_process(fmwk.ccnc(is_nc))
my_proc.run()
print
print "Output in",
print file_name
print
