def getCouplingsFromSample(sample):
couplings = ["DC1A", "DC1V", "DC2A", "DC2V", "ctZ", "ctZI", "cpt", "cpQM"]
l = []
for c in couplings:
val = getCouplingFromName(sample.name, c)
if val:
l.append((c, val))
else:
l.append((c, 0))
return l
y_par = ''
allowedIntervals = []
indirectConstraints = []
else:
raise NotImplementedError
else:
raise NotImplementedError
print "Number of results", len(signals)
var1_values = []
for s in signals:
s.var1 = getCouplingFromName(s.name, x_var)
if s.var1 != 0.0 and (s.var1 + x_shift) not in var1_values:
var1_values.append(s.var1 + x_shift)
var1_values = sorted(var1_values)
x = []
x_forRange = []
z = []
res_dic = {}
print "Searching for bestfit point"
bestNLL = 9999.