Merge pull request cms-sw#29 from amarini/topic_optimization

Topic optimization

script/optimization.py

@@ -289,9 +289,10 @@ def Combine(h1,h2,postfix="",dir="/tmp/%(USER)s"):
 	h1.Clone("sig").Write()
 	h2.Clone("bkg").Write()
 	data= h1.Clone("data") ## no data, but needed to convert datacard into ws, I put an asimov instead
-	data.Add(bkg)
+	data.Add(h2)
 	data.Write()
 	f.Close()
+	print "writing datacard in ",wdir+"datacard%s.txt"%postfix 
 	dat = open(wdir+"datacard%s.txt"%postfix ,"w")
 
 	dat.write("-------------------------------------\n")
@@ -343,8 +344,8 @@ def CleanForCombine(px_s,px_b):
 	RemoveSpikes(px_s,False,3)
 	RemoveSpikes(px_b,False,3)
 
-	px_s.Smooth(20)
-	px_b.Smooth(20)
+	#px_s.Smooth(20)
+	#px_b.Smooth(20)
 	return 
 
 

script/parallel.py

@@ -0,0 +1,37 @@
+#!env python
+
+#import os,sys
+from subprocess import call, check_output
+#import re
+import threading
+import time
+
+from optparse import OptionParser, OptionGroup
+
+usage = "usage: %prog [options] exe1 exe2 ..."
+parser=OptionParser(usage=usage)
+parser.add_option("-n","--ncore" ,dest='ncore',type='int',help="Input Datacard [Default=%default]",default=10)
+
+def call_exe(exe):
+	print " -> Calling "+ exe
+	call(exe,shell=True)
+	print " -> Done "+exe
+
+(opts,args)=parser.parse_args()
+
+threads=[]
+
+for exe in args:
+	while threading.activeCount() >= opts.ncore: 
+		print "sleep ...",
+		time.sleep(3)
+		print "wake up"
+	t = threading.Thread(target=call_exe, args= (exe,)  )
+	t.start()
+	threads.append( t )
+
+for t in threads:
+	t.join()
+
+print "Done"
+