os.system('ln -s ../../Source/pdf.inc') os.chdir('..') os.chdir("..") return #####################################################################################" ## ## LAUNCH PROGRAM ## ###################################################################################### if (__name__ == '__main__'): go_to_main_dir() check_for_help(sys.argv) MWparam = MW_info('MadWeight_card.dat') MWparam.set_run_opt(sys.argv) for dir in MWparam.MW_listdir: if not os.path.exists("SubProcesses/" + dir + "/call_TF.f"): print("Currently no transfer function loaded ...") listdir = os.listdir('./Source/MadWeight/transfer_function/data') print('Available transfer functions:\n ', end=' ') print('\n '.join([ content[3:-4] for content in listdir if (content.startswith('TF') and content.endswith('dat')) ])) name = input('Choose your transfer Function\n') P_dir, MW_dir = detect_SubProcess(P_mode=1) os.chdir('./Source/MadWeight/transfer_function') change_tf.create_TF_main(name, 0, MW_dir)
f(x)="""+self.function+""" fit f(x) '"""+graph.tag+""".dat' u 1:2:3 via a,b,c """ return text ################################################################################################ ## MAIN ## ################################################################################################ if (__name__=="__main__"): from MW_param import go_to_main_dir,MW_info go_to_main_dir() MWparam=MW_info('MadWeight_card.dat') MWparam.name='fermi' if 0: obj=Likelihood(MWparam,auto=1) elif 1: obj=Likelihood(MWparam,auto=0) file_list=['./Events/fermi2/fermi2_weights.out','./Events/fermi3/fermi3_weights.out'] obj.charge_weight(filename=file_list) #self.info['mw_run']['4'] = normalize with cross obj.correct_with_cross('./Events/fermi2/fermi2_cross_weights.out') obj.correct_with_acc() obj.define_plots(MWparam.info['mw_parameter']['1']) #option -> gives how param card are generated if obj.run_type=='1': obj.define_section_plots() obj.plot() elif 1:
os.system('ln -s ../../Source/pdf.inc') os.chdir('..') os.chdir("..") return #####################################################################################" ## ## LAUNCH PROGRAM ## ###################################################################################### if(__name__=='__main__'): go_to_main_dir() check_for_help(sys.argv) MWparam=MW_info('MadWeight_card.dat') MWparam.set_run_opt(sys.argv) for dir in MWparam.MW_listdir: if not os.path.exists("SubProcesses/"+dir+"/call_TF.f"): print "Currently no transfer function loaded ..." listdir=os.listdir('./Source/MadWeight/transfer_function/data') print 'Available transfer functions:\n ', print '\n '.join([content[3:-4] for content in listdir if (content.startswith('TF') and content.endswith('dat'))]) name=raw_input('Choose your transfer Function\n') P_dir,MW_dir=detect_SubProcess(P_mode=1) os.chdir('./Source/MadWeight/transfer_function') change_tf.create_TF_main(name,0,MW_dir) print os.getcwd() os.chdir('../../..') break
f(x)=""" + self.function + """ fit f(x) '""" + graph.tag + """.dat' u 1:2:3 via a,b,c """ return text ################################################################################################ ## MAIN ## ################################################################################################ if (__name__ == "__main__"): from MW_param import go_to_main_dir, MW_info go_to_main_dir() MWparam = MW_info('MadWeight_card.dat') MWparam.name = 'fermi' if 0: obj = Likelihood(MWparam, auto=1) elif 1: obj = Likelihood(MWparam, auto=0) file_list = [ './Events/fermi2/fermi2_weights.out', './Events/fermi3/fermi3_weights.out' ] obj.charge_weight(filename=file_list ) #self.info['mw_run']['4'] = normalize with cross obj.correct_with_cross('./Events/fermi2/fermi2_cross_weights.out') obj.correct_with_acc() obj.define_plots(MWparam.info['mw_parameter']
gg.writelines(param + "." + event + "\t" + str(0) + "\t" + str(0) + "\n") ff.close() gg.close() #c) write the cross section file hh = open( './Events/' + run_name + '/' + run_name + '_cross_weights.out', 'w') for param in self.info.actif_param: if total_cross[str(param)][0]: hh.writelines( str(param) + "\t" + str(total_cross[str(param)][0]) + "\t" + str(sqrt(total_cross[str(param)][1])) + "\n") hh.close() return ########################################################################### ################### TEST ################## ############################################################################ if __name__ == '__main__': print os.getcwd() go_to_main_dir() MW_content = MW_info('MadWeight_card.dat') #collect_schedular(MW_content) Collect_All_dir(MW_content)
print "collecting data" collect_schedular(MWparam) return list_pass_event #####################################################################################" ## ## LAUNCH PROGRAM ## ###################################################################################### if (__name__ == '__main__'): go_to_main_dir() check_for_help(sys.argv) MWparam = MW_info('MadWeight_card.dat') MWparam.set_run_opt(sys.argv) MWparam.old_name = MWparam['mw_acc']['old_run'] if MWparam.run_opt['param']: create_param.Param_card(run_name=MWparam) MWparam.update_nb_card() Cards.create_include_file(MWparam) update_cuts_status(MWparam) pass_event = Launch_all_SubProcess(MWparam) if MWparam.run_opt['plot']: passing_event('./Events/'+MWparam.old_name+'/'+MWparam.old_name+'_weights.out',\ './Events/'+MWparam.name+'/'+MWparam.name+'_weights.out' , pass_event)
if MWparam.run_opt['collect']: print "collecting data" collect_schedular(MWparam) return list_pass_event #####################################################################################" ## ## LAUNCH PROGRAM ## ###################################################################################### if(__name__=='__main__'): go_to_main_dir() check_for_help(sys.argv) MWparam=MW_info('MadWeight_card.dat') MWparam.set_run_opt(sys.argv) MWparam.old_name=MWparam['mw_acc']['old_run'] if MWparam.run_opt['param']: create_param.Param_card(run_name=MWparam) MWparam.update_nb_card() Cards.create_include_file(MWparam) update_cuts_status(MWparam) pass_event=Launch_all_SubProcess(MWparam) if MWparam.run_opt['plot']: passing_event('./Events/'+MWparam.old_name+'/'+MWparam.old_name+'_weights.out',\ './Events/'+MWparam.name+'/'+MWparam.name+'_weights.out'