###################################################################### 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,