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