######################################################################
print "***** stau1 mass =", str(m1_fixed), "GeV, tan(beta) =", str(
tb), "*****\n"
print "***** scan initialization *****"
m2_min = max(100, m1_fixed)
m2_max = 1000.
theta_min = 0.
theta_max = pi / 2.
# Prepare output
fresults = open(output, 'w')
# Initialize a Lilith object
lilithcalc = lilith.Lilith(verbose=False, timer=False)
# Read experimental data
lilithcalc.readexpinput(exp_input)
######################################################################
# Scan routine
######################################################################
m2logLmin = 10000
max = -1
print "***** running scan *****"
for theta in np.linspace(theta_min, theta_max, grid_subdivisions):
fresults.write('\n')
for m2 in np.linspace(m2_min, m2_max, grid_subdivisions):
CGa=CGa,
Cg=Cg,
BRinv=BRinv,
precision=myprecision)
lilithcalc.computelikelihood(userinput=myXML_user_input)
return lilithcalc.l
######################################################################
# Calculations
######################################################################
print("***** initializating *****")
# Initialize a Lilith object
lilithcalc = lilith.Lilith(verbose, timer)
# Read experimental data
lilithcalc.readexpinput(myexpinput)
# Initialize the fit; parameter starting values and limits
if iminuit_version_f < 2.0:
m = Minuit(getL,
CGa=1,
limit_CGa=(0, 3),
Cg=1,
limit_Cg=(0, 3),
BRinv=0.2,
limit_BRinv=(0, 0.9),
errordef=1,
error_CGa=0.1,
error_Cg=0.1,
