if sys.argv[i].find("--signif")!=-1: significance = True
if sys.argv[i].find("--work")!=-1: workspaceFlag = True
if sys.argv[i].find("--no-work")!=-1: workspaceFlag = False
if sys.argv[i].find("--expected-a-priori")!=-1: expected_a_priori = True
if refXsecFile is not None:
print "INFO: Input ref xsec file!"
gluinoFile = rt.TFile.Open(refXsecFile,"READ")
gluinoHistName = refXsecFile.split("/")[-1].split(".")[0]
gluinoHist = gluinoFile.Get(gluinoHistName)
nJobs = 1 # do 1 toy each job => 1 toy
print box, model, queue
gchipairs = getGChiPairs(model)
gchipairs = reversed(gchipairs)
pwd = os.environ['PWD']
submitDir = "submit"+model+fitRegion+box
outputDir = "output"+model+fitRegion+box
os.system("mkdir -p %s"%(submitDir))
os.system("mkdir -p %s"%(outputDir))
# for compting what jobs are left:
doneFile = open(done)
if significance and nToys<0:
outFileList = [outFile.replace("higgsCombine","").replace(".ProfileLikelihood.mH120.root\n","") for outFile in doneFile.readlines() if outFile.find("ProfileLikelihood")!=-1]
nToysPerJob = int(nToys/nJobs)
nXsec = 5 # do 5 xsec points + 0 lower values + 1 higher value
if asymptoticFile is not None:
print "INFO: Input ref xsec file!"
asymptoticRootFile = rt.TFile.Open(asymptoticFile,"READ")
expMinus2 = asymptoticRootFile.Get("xsecUL_ExpMinus2_%s_%s"%(model,box))
expPlus2 = asymptoticRootFile.Get("xsecUL_ExpPlus2_%s_%s"%(model,box))
expMinus = asymptoticRootFile.Get("xsecUL_ExpMinus_%s_%s"%(model,box))
expPlus = asymptoticRootFile.Get("xsecUL_ExpPlus_%s_%s"%(model,box))
exp = asymptoticRootFile.Get("xsecUL_Exp_%s_%s"%(model,box))
print box, model, queue
gchipairs = getGChiPairs(model,fine=False)
gchipairs = reversed(gchipairs)
pwd = os.environ['PWD']
submitDir = "submit"+model+fitRegion+box
outputDir = "output"+model+fitRegion+box
os.system("mkdir -p %s"%(submitDir))
os.system("mkdir -p %s"%(outputDir))
# for compting what jobs are left:
doneFile = open(done)
outFileList = []
if step2:
