-
Notifications
You must be signed in to change notification settings - Fork 8
/
submitAllLimits.py
83 lines (73 loc) · 3.51 KB
/
submitAllLimits.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
import os,sys
import time
masses =[m*100 for m in range(10,45+1)]
models = ["altqW"]
#channels = ["qVHPnew","qVLPnew"]
channels =["qVnew","qWHP","qWLP","qZHP","qZLP"]
#workspaceBKG = "CMS_jj_bkg_qV_13TeV"
workspaceBKG = "CMS_jj_bkg_qValt_13TeV"
for channel in channels:
for model in models:
if model.find("q") != -1:
masses =[m*100 for m in range(12,60+1)]
#masses = [4000]
for mass in masses:
#workspaceSignal = "CMS_jj_"+model+"_"+str(int(mass))+"_13TeV"
workspaceSignal = "CMS_jj_qW_"+str(int(mass))+"_13TeV"
command = 'qsub -q all.q submitLimits.sh CMS_jj_'+model+'_'+str(int(mass))+'_13TeV_CMS_jj_'+channel+' '+str(int(mass))+' CMS_jj_'+str(int(mass))+'_'+model+'_13TeV_CMS_jj_'+channel+'_asymptoticCLs_new.root '+workspaceBKG+' '+workspaceSignal
print command
os.system(command)
#time.sleep(1)
#masses =[m*100 for m in range(12,43+1)]
##models = ["ZprimeWW","BulkWW","BulkZZ","WZ"]
#models = ["BulkWW"]
##channels = ["WWHP","WWLP","ZZHP","ZZLP","WZHP","WZLP","VVnew"]
##channels = ["ZZHP","ZZLP","VVnew"]
##models =["BulkWW"]
##masses =[4000]
##channels = ["VVHPnew","VVLPnew","VVnew"]
##channels = ["VVnew"]
#channels = ["WWHP","WWLP","ZZHP","ZZLP","WZHP","WZLP"]
#workspaceBKG = " CMS_jj_bkg_VV_13TeV"
#for channel in channels:
#for model in models:
#if model.find("q") != -1:
#masses =[m*100 for m in range(12,60+1)]
#for mass in masses:
#workspaceSignal = "CMS_jj_"+model+"_"+str(int(mass))+"_13TeV"
#command = 'qsub -q all.q submitLimits.sh CMS_jj_'+model+'_'+str(int(mass))+'_13TeV_CMS_jj_'+channel+' '+str(int(mass))+' CMS_jj_'+str(int(mass))+'_'+model+'_13TeV_CMS_jj_'+channel+'_asymptoticCLs_new.root '+workspaceBKG+' '+workspaceSignal
#print command
#os.system(command)
##time.sleep(1)
#masses =[m*100 for m in range(12,43+1)]
#models = ["qW"]
##channels = ["qVnew"]
#channels =["qWHP","qWLP","qZHP","qZLP","qVnew"]
#workspaceBKG = "CMS_jj_bkg_qV_13TeV"
#for channel in channels:
#for model in models:
#if model.find("q") != -1:
#masses =[m*100 for m in range(12,60+1)]
#for mass in masses:
#workspaceSignal = "CMS_jj_"+model+"_"+str(int(mass))+"_13TeV"
#command = 'qsub -q all.q submitLimits.sh CMS_jj_'+model+'_'+str(int(mass))+'_13TeV_CMS_jj_'+channel+' '+str(int(mass))+' CMS_jj_'+str(int(mass))+'_'+model+'_13TeV_CMS_jj_'+channel+'_asymptoticCLs_new.root '+workspaceBKG+' '+workspaceSignal
#print command
#os.system(command)
##time.sleep(1)
## alternaive fit function !
#masses =[m*100 for m in range(12,43+1)]
##masses =[4100]
#models = ["altBulkWW"]
##channels = ["WWHP","WWLP","ZZHP","ZZLP","WZHP","WZLP"]
#channels = ["VVnew"]
#workspaceBKG = " CMS_jj_bkg_VValt_13TeV"
#for channel in channels:
#for model in models:
#if model.find("q") != -1:
#masses =[m*100 for m in range(12,41+1)]
#for mass in masses:
#workspaceSignal = "CMS_jj_"+model+"_"+str(int(mass))+"_13TeV"
#command = 'qsub -q all.q submitLimits.sh CMS_jj_'+model+'_'+str(int(mass))+'_13TeV_CMS_jj_'+channel+' '+str(int(mass))+' CMS_jj_'+str(int(mass))+'_'+model+'_13TeV_CMS_jj_'+channel+'_asymptoticCLs_new.root '+workspaceBKG+' '+workspaceSignal
#print command
#os.system(command)
##time.sleep(1)