def scan_top_scalingttH_floatingBRs_constrainedbbZZ(args):
config = basic_config(args)
config.datacard = approval.ktcg.floatingBRs[config.decay_channel]
config.tags.append('scalingttH')
config.tags.append('floatingBRs')
config.tags.append('constrainedbbZZ')
# config.nPoints = 30*30
# Around 4.5 hours (good max of 6 hours) for 8 points (combWithHbb, observed)
# Reasonably constant, can go up to 12 points
# Keep same ranges for hgg/hzz for now
# config.nPointsPerJob = 12
# config.set_parameter_range('ct', -4.2, 4.2)
# config.set_parameter_range('cg', -0.30, 0.30)
config.nPoints = 100 * 100
config.nPointsPerJob = 12
config.set_parameter_range('ct', -6.0, 6.0)
config.set_parameter_range('cg', -0.45, 0.45)
if args.hzz:
config.nPointsPerJob = 300
config.queue = 'short.q'
if args.asimov:
config.nPointsPerJob = 10
config.set_parameter_range('ct', -8.2, 8.2)
config.set_parameter_range('cg', -0.50, 0.50)
differentialutils.run_postfit_scan(config)
def scan_top_scalingttH_cgonly(args):
args = differentialutils.set_one_decay_channel(args, 'combWithHbb')
config = basic_config(args)
# config.datacard = 'out/workspaces_May17/combWithHbb_Top_reweighted_scalingttH.root'
config.datacard = 'out/workspaces_May22/combWithHbb_Top_reweighted_scalingttH_floatingBRs.root'
config.tags.append('scalingttH')
config.tags.append('cgonly')
config.suppress_output = True
config.queue = 'short.q'
config.FLOAT_OTHER_POIS = False
config.del_parameter_range('ct')
config.del_poi('ct')
config.set_parameter('ct', 1.0, hard=True)
config.freeze_parameter('ct')
if args.asimov:
config.set_parameter_range('cg', -1.6, 0.2)
else:
config.set_parameter_range('cg', -0.1, 0.4)
config.nPoints = 25
config.nPointsPerJob = 1
# differentialutils.run_postfit_fastscan_scan(config)
differentialutils.run_postfit_scan(
config,
postfit_file=
'out/Scan_May25_Top_combWithHbb_scalingttH_cgonly_asimov/postfit_and_fastscan/higgsCombine_POSTFIT_ASIMOV_combWithHbb_Top_reweighted_scalingttH_floatingBRs.MultiDimFit.mH125.root'
)
def scan_yukawa_G2A(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.nPoints = 50*50
config.set_parameter_range('kappab', -40., 40.)
config.set_parameter_range('kappac', -70., 70.)
config.tags.append('G2A')
config.datacard = yukawa_G.G2A
differentialutils.run_postfit_scan(config)
def scan_topctcb_scalingbbHttH_floatingBRs(args):
args = differentialutils.set_one_decay_channel(args, 'combWithHbb')
config = basic_config_ctcb(args)
config.datacard = 'out/workspaces_May29/combWithHbb_TopCtCb_reweighted_scalingbbHttH_floatingBRs.root'
config.tags.append('scalingbbHttH')
config.tags.append('floatingBRs')
config.set_parameter_range('ct', -10.0, 10.0)
config.set_parameter_range('cb', -20., 20.)
config.nPoints = 40 * 40
config.nPointsPerJob = 9
differentialutils.run_postfit_scan(config)
def scan_top_scalingttH_floatingBRs_NONconstrainedbbZZ(args):
args = differentialutils.set_one_decay_channel(args, 'combWithHbb')
config = basic_config(args)
config.datacard = 'out/workspaces_May30/combWithHbb_Top_reweighted_scalingttH_floatingBRs.root'
config.tags.append('scalingttH')
config.tags.append('floatingBRs')
config.tags.append('NONconstrainedbbZZ')
config.set_parameter_range('ct', -1.7, 1.7)
config.set_parameter_range('cg', -0.10, 0.10)
# config.nPoints = 10*10
# config.nPointsPerJob = 1
# config.queue = 'short.q'
config.nPoints = 30 * 30
config.nPointsPerJob = 7
config.queue = 'all.q'
differentialutils.run_postfit_scan(config)
def scan_top_BRdependent(args):
args = set_combWithHbb_and_asimov(args)
config = basic_config(args)
config.datacard = LatestPaths.ws.top.BRcouplingDependency
config.nPoints = 80 * 80
config.nPointsPerJob = 10
config.PhysicsModelParameterRanges = [
'ct={0},{1}'.format(-8.5, 8.5), 'cg={0},{1}'.format(-0.65, 0.65)
]
config.tags.append('couplingDependentBR')
config.PhysicsModelParameters.append('kappa_V=0.999')
# config.freezeNuisances.append('kappa_V')
config.floatNuisances.append('kappa_V')
config.PhysicsModelParameterRanges.append('kappa_V=-1000.0,1.0')
config.tags.append('floatKappaV')
differentialutils.run_postfit_scan(config)
def scan_topctcb_scalingbbHttH_floatingBRs_constrainedbbZZ(args):
config = basic_config_ctcb(args)
config.datacard = approval.ktkb.floatingBRs[config.decay_channel]
config.tags.append('scalingbbHttH')
config.tags.append('floatingBRs')
config.tags.append('constrainedbbZZ')
# # around max 5, occasionally 5.5 hours for 7 points; can go up to 12
# config.set_parameter_range('ct', -7.0, 7.0)
# config.set_parameter_range('cb', -17., 17.)
# config.nPoints = 25*25
# config.nPointsPerJob = 10
config.set_parameter_range('ct', -3.5, 3.5)
config.set_parameter_range('cb', -15.5, 15.5)
config.nPoints = 60 * 60
config.nPointsPerJob = 10
if args.hzz:
config.nPointsPerJob = 300
config.queue = 'short.q'
differentialutils.run_postfit_scan(config)
def scan_yukawa_BRdependent(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
config.datacard = LatestPaths.ws.yukawa.BRcouplingDependency
config.nPoints = 60*60
config.nPointsPerJob = 10
config.PhysicsModelParameterRanges = [
'kappab={0},{1}'.format(-25., 25),
'kappac={0},{1}'.format(-45., 45)
]
config.tags.append('couplingDependentBR')
config.PhysicsModelParameters.append('kappa_V=0.999')
# config.freezeNuisances.append('kappa_V')
config.floatNuisances.append('kappa_V')
config.PhysicsModelParameterRanges.append('kappa_V=-1000.0,1.0')
config.tags.append('floatKappaV')
differentialutils.run_postfit_scan(config)
def scan_yukawa_NONscalingbbH_floatingBRs(args):
# args = differentialutils.set_one_decay_channel(args, 'combination')
config = basic_config(args)
config.tags.append('NONscalingbbH')
config.tags.append('floatingBRs')
config.set_parameter_range('kappab', -50., 50.)
config.set_parameter_range('kappac', -90., 90.)
config.nPoints = 35*35
# For asimov, around 15 min per point for the combination
config.nPointsPerJob = 5
config.queue = 'short.q'
if args.hzz:
config.set_parameter_range('kappab', -10., 20.)
config.set_parameter_range('kappac', -50., 50.)
config.nPointsPerJob = 300
config.queue = 'short.q'
config.datacard = approval.floatingBRs[config.decay_channel]
differentialutils.run_postfit_scan(config)
def scan_top_scalingttH_floatingBRs(args):
raise RuntimeError('Probably not what I want!')
args = differentialutils.set_one_decay_channel(args, 'combWithHbb')
config = basic_config(args)
# config.datacard = 'out/workspaces_May17/combWithHbb_Top_reweighted_scalingttH.root'
config.datacard = 'out/workspaces_May22/combWithHbb_Top_reweighted_scalingttH_floatingBRs.root'
config.tags.append('scalingttH')
config.tags.append('floatingBRs')
# Runs out of CPU
# config.set_parameter_range('ct', -2.5, 2.5)
# config.set_parameter_range('cg', -0.15, 0.15)
# config.nPoints = 40*40
# config.nPointsPerJob = 12
# Even faster
config.set_parameter_range('ct', -1.7, 1.7)
config.set_parameter_range('cg', -0.10, 0.10)
config.nPoints = 10 * 10
config.nPointsPerJob = 4
# differentialutils.run_postfit_fastscan_scan(config)
differentialutils.run_postfit_scan(config)
def projection_ktcg_scan_floatingBRs(args):
Config = CombineConfigKTCG_floatingBRs_scenario2 if args.scenario2 else CombineConfigKTCG_floatingBRs
config = Config(args)
differentialutils.run_postfit_scan(config)
def projection_kbkc_scan_floatingBRs(args):
args = differentialutils.force_asimov(args)
Config = CombineConfigKBKC_floatingBRs_scenario2 if args.scenario2 else CombineConfigKBKC_floatingBRs
config = Config(args)
differentialutils.run_postfit_scan(config)
