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: