def GenerateProfLHoodInterval(workspaceFileName): print "*"*80 print "Generating Profile Likelihood Interval" print "*"*80 print "Looking at workspace in file:",workspaceFileName outDir = "results_proflhood" print outDir mkdirs(outDir) command = "code/bin/profileLHoodInterval "+workspaceFileName +" "+outDir system(command) return outDir
def ExtractSignificance(workspaceFileName): print "*"*80 print "Extract Significance" print "*"*80 print "Looking at workspace in file:",workspaceFileName outDir = "results_significance" print outDir mkdirs(outDir) command = "code/bin/extractSignificance "+workspaceFileName +" "+outDir print command system(command) return outDir
minusTwoSigma= intervalHist.GetBinContent(2)
minusOneSigma= intervalHist.GetBinContent(3)
plusOneSigma= intervalHist.GetBinContent(4)
plusTwoSigma= intervalHist.GetBinContent(5)
print "Final Measurement:"
plusOneSigmaDelta = (plusOneSigma-measuredValue)
minusOneSigmaDelta = -(minusOneSigma-measuredValue)
measuredPrintOut = str(measuredValue)+" + "+ str(plusOneSigmaDelta) +" - "+ str(minusOneSigmaDelta)
print measuredPrintOut
if outputDir!="":
mkdirs(outputDir)
ofile = open(outputDir+"/measuredValue.txt",'w')
ofile.write(measuredPrintOut)
ofile.close()
canvas = TCanvas("canvas","canvas",600,600)
object =None
index = 0
rooCurveTotName=""
rooCurveTot = None
#don't remove immediately because it change the size of the list
objectIndicesToRemove = []
while True:
object = rooPlotTot.getObject(index)
