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'