def scan_top_scalingttH_couplingdependentBRs(args):
config = basic_config(args)
config.datacard = approval.ktcg.couplingdependentBRs[config.decay_channel]
config.tags.append('scalingttH')
config.tags.append('couplingdependentBRs')
config.nPoints = 30 * 30
# On combWithHbb, observed: around 4 hours for 10 points, but lots of variance
# Good maximum per job is 10 (tiny bit of loss)
config.nPointsPerJob = 10
config.set_parameter_range('ct', -0.4, 4.4)
config.set_parameter_range('cg', -0.28, 0.1)
if args.hgg:
config.set_parameter_range('ct', -0.3, 2.6)
config.set_parameter_range('cg', -0.08, 0.12)
elif args.hzz:
config.set_parameter_range('ct', -0.8, 3.2)
config.set_parameter_range('cg', -0.15, 0.16)
config.nPointsPerJob = 300
config.queue = 'short.q'
if args.asimov and args.hgg:
# Special scan for Vittorio
config.nPoints = 20 * 20
config.nPointsPerJob = 5
config.set_parameter_range('ct', -0.3, 5.3)
config.set_parameter_range('cg', -0.12, 0.2)
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_BRdependent_and_profiledTotalXS(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
config.datacard = LatestPaths.ws_combined_Yukawa_couplingDependentBR_profiledTotalXS
config.subDirectory += '_couplingDependentBR_profiledTotalXS'
config.fix_parameter_at_value('kappa_V', 0.999)
differentialutils.run_postfit_fastscan_scan(config)
def scan_top(args):
config = basic_config(args)
config.datacard = LatestPaths.ws.top.nominal[
differentialutils.get_decay_channel_tag(args)]
config.nPoints = 110 * 110
config.nPointsPerJob = int(0.5 * config.nPointsPerJob)
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_profiledTotalXS(args):
args = set_combWithHbb_and_asimov(args)
config = basic_config(args)
config.datacard = 'out/workspaces_Mar13/combWithHbb_Top_reweighted_profiledTotalXS.root'
config.tags.append('profiledTotalXS')
config.nPoints = 150 * 150
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G1A(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.set_parameter_range('kappab', -7., 7.)
config.set_parameter_range('kappac', -14., 14.)
config.tags.append('G1A')
config.datacard = yukawa_G.G1A
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G0B(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.queue = 'short.q'
config.nPointsPerJob = 600
config.tags.append('G0B')
config.datacard = yukawa_G.G0B
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_fitOnlyNormalization(args):
differentialutils.assert_asimov(args)
config.datacard = (
LatestPaths.ws_combined_Top_profiledTotalXS_fitOnlyNormalization
if not args.highpt else
LatestPaths.ws_combined_TopHighPt_profiledTotalXS_fitOnlyNormalization)
config.tags.append('fitOnlyNormalization')
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_unreweighted(args):
config = basic_config(args)
datacard_dict = LatestPaths.ws.yukawa.unreweighted
config.datacard = datacard_dict[differentialutils.get_decay_channel_tag(args)]
# if args.combination:
# config.datacard = 'out/workspaces_Mar09/combination_Yukawa_nominal.root'
# elif args.hgg:
# config.datacard = 'out/workspaces_Mar09/hgg_Yukawa_nominal.root'
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_scalingttH(args):
config = basic_config(args)
config.datacard = approval.ktcg.fixedBRs[config.decay_channel]
config.tags.append('scalingttH')
config.nPoints = 30 * 30
config.set_parameter_range('ct', -0.4, 4.4)
config.set_parameter_range('cg', -0.28, 0.1)
config.nPointsPerJob = 6
differentialutils.run_postfit_fastscan_scan(config)
def scan_topctcb_scalingbbHttH(args):
config = basic_config_ctcb(args)
config.datacard = approval.ktkb.fixedBRs[config.decay_channel]
config.tags.append('scalingbbHttH')
config.set_parameter_range('ct', -1.55, 1.55)
config.set_parameter_range('cb', -20., 20.)
config.nPoints = 50 * 50
config.nPointsPerJob = 12
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_hzz_testrun(real_args):
args = deepcopy(real_args)
differentials.core.set_one_decay_channel(args, 'hzz')
args.asimov = True
config = basic_config(args)
config.datacard = LatestPaths.ws.top.nominal[
differentialutils.get_decay_channel_tag(args)]
config.nPoints = 10 * 10
config.nPointsPerJob = 50
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G1B(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.queue = 'short.q'
config.nPointsPerJob = 600
config.set_parameter_range('kappab', -3., 3.)
config.set_parameter_range('kappac', -10., 10.)
config.tags.append('G1B')
config.datacard = yukawa_G.G1B
differentialutils.run_postfit_fastscan_scan(config)
def test_dopoints_option_implementation(args):
args = differentialutils.set_one_decay_channel(args, 'hzz', asimov=True)
decay_channel = differentialutils.get_decay_channel_tag(args)
config = CombineConfigKBKC(args)
config.datacard = workspaces.couplingdependentBRs[decay_channel]
config.tags.append('couplingdependentBRs')
config.range_couplingdependentBRs(args)
config.nPoints = 10 * 10
config.nPointsPerJob = 50
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_lumiStudy(args):
args = set_combWithHbb_and_asimov(args)
config = basic_config(args)
config.datacard = 'out/workspaces_Mar13/combWithHbb_Top_reweighted_lumiScale.root'
config.freezeNuisances.append('lumiScale')
config.hardPhysicsModelParameters.append('lumiScale=8.356546')
config.subDirectory += '_lumiStudy'
# config.hardPhysicsModelParameters.append( 'lumiScale=83.56546' )
# config.subDirectory += '_lumiStudy3000fb'
differentialutils.run_postfit_fastscan_scan(config)
def scan_topctcb_lumiscale(args):
args = differentialutils.set_one_decay_channel(args,
'combWithHbb',
asimov=True)
config = basic_config_ctcb(args)
config.datacard = LatestPaths.ws.topctcb.lumiScale
config.freezeNuisances.append('lumiScale')
config.hardPhysicsModelParameters.append('lumiScale=8.356546')
config.tags.append('lumiStudy')
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_lumiScale(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
config.datacard = LatestPaths.ws.yukawa.lumiScale
config.freezeNuisances.append('lumiScale')
config.hardPhysicsModelParameters.append( 'lumiScale=8.356546' )
config.tags.append('lumi300fb')
# config.hardPhysicsModelParameters.append( 'lumiScale=83.56546' )
# config.tags.append('lumi3000fb')
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_profiledTotalXS(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
# config.datacard = LatestPaths.ws_combined_Yukawa_profiledTotalXS
config.datacard = LatestPaths.ws.yukawa.profiledTotalXS
config.tags.append('profiledTotalXS')
config.nPoints = 60*60
config.nPointsPerJob = 15
config.deltaNLLCutOff = 120.
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G0B_reweighted(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.nPoints = 50*50
config.queue = 'short.q'
config.nPointsPerJob = 600
config.set_parameter_range('kappab', -5., 8.)
config.set_parameter_range('kappac', -20., 20.)
config.tags.append('G0B_reweighted')
config.datacard = yukawa_G.G0B_reweighted
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_lastBinDropped(real_args):
args = deepcopy(real_args)
args.asimov = True
config = basic_config(args)
if args.combination:
config.datacard = 'out/workspaces_Mar07/combination_Top_reweighted_lastBinDroppedHgg.root'
elif args.combWithHbb:
config.datacard = 'out/workspaces_Mar07/combWithHbb_Top_reweighted_lastBinDroppedHgg.root'
else:
raise NotImplementedError('Run with --combination or --combWithHbb')
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_last2BinsDropped(real_args):
args = deepcopy(real_args)
args.asimov = True
config = basic_config(args)
if args.combWithHbb:
config.datacard = LatestPaths.ws.top.last2BinsDropped.combWithHbb
elif args.combination:
config.datacard = LatestPaths.ws.top.last2BinsDropped.combination
else:
raise NotImplementedError('Run with --combination or --combWithHbb')
config.tags.append('last2BinsDropped')
config.nPoints = 150 * 150
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G1BKV(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.queue = 'short.q'
config.nPointsPerJob = 350
config.set_parameter_range('kappab', -7., 7.)
config.set_parameter_range('kappac', -14., 14.)
config.set_parameter_range('kappa_V', -10000., 1.)
config.hardPhysicsModelParameters.append('kappa_V=0.999')
config.floatNuisances.append('kappa_V')
config.tags.append('G1BKV')
config.datacard = yukawa_G.G1B
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_G1A_reweighted_noTheoryUnc_scaledByMuTotalXS(args):
args = differentialutils.set_one_decay_channel(args, 'combination', asimov=True)
config = basic_config(args)
config.nPoints = 50*50
config.set_parameter_range('kappab', -2., 2.)
config.set_parameter_range('kappac', -8., 8.)
config.nPointsPerJob = 16
config.queue = 'short.q'
config.tags.append('G1A_reweighted_noTheoryUnc_scaledByMuTotalXS')
config.datacard = yukawa_G.G1A_reweighted_noTheoryUnc_scaledByMuTotalXS
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_NONscalingbbH_couplingdependentBRs(args):
# args = differentialutils.set_one_decay_channel(args, 'combination')
config = basic_config(args)
config.tags.append('NONscalingbbH')
config.tags.append('couplingdependentBRs')
config.set_parameter_range('kappab', -2., 2.)
config.set_parameter_range('kappac', -8., 8.)
config.nPoints = 30*30
config.nPointsPerJob = 5
config.queue = 'short.q'
if args.hzz:
config.set_parameter_range('kappab', -5., 5.)
config.set_parameter_range('kappac', -15., 15.)
config.nPointsPerJob = 300
config.queue = 'short.q'
config.datacard = approval.couplingdependentBRs[config.decay_channel]
differentialutils.run_postfit_fastscan_scan(config)
def scan_topctcb_scalingbbHttH_couplingdependentBRs(args):
config = basic_config_ctcb(args)
config.datacard = approval.ktkb.couplingdependentBRs[config.decay_channel]
config.tags.append('scalingbbHttH')
config.tags.append('couplingdependentBRs')
config.nPoints = 35 * 35
# Up to max ~7 hours for 10 points (combWithHbb, observed), can go up to 12 points
config.set_parameter_range('ct', -0.3, 2.5)
config.set_parameter_range('cb', -2.8, 2.8)
config.nPointsPerJob = 12
if args.hgg:
# kt-plane asym means hgg goes up to reasonably high kt
config.nPoints = 25 * 25
config.set_parameter_range('ct', 0.2, 3.3)
config.set_parameter_range('cb', -1.8, 1.8)
config.nPointsPerJob = 10
if args.hzz:
config.set_parameter_range('ct', -3.3, 3.3)
config.set_parameter_range('cb', -3.3, 3.3)
config.nPointsPerJob = 300
config.queue = 'short.q'
differentialutils.run_postfit_fastscan_scan(config)
def projection_ktcg_scan_couplingdependentBRs(args):
Config = CombineConfigKTCG_couplingdependentBRs_scenario2 if args.scenario2 else CombineConfigKTCG_couplingdependentBRs
config = Config(args)
# if args.combWithHbb:
# # config.minimizer_settings_reasonable()
# # config.freeze_mcstat_uncertainties()
# # if args.scenario2:
# # reusable_postfit = 'out/Scan_projection_ktcg_Oct04_combWithHbb_scenario2_couplingdependentBRs_asimov/postfit_and_fastscan/higgsCombine_POSTFIT_ASIMOV_projection_ktcgkb_combWithHbb_s2groups_Jul18_couplingModel_reweighted_scenario2_couplingdependentBRs.MultiDimFit.mH125.root'
# # else:
# # reusable_postfit = 'out/Scan_projection_ktcg_Oct04_combWithHbb_couplingdependentBRs_asimov/postfit_and_fastscan/higgsCombine_POSTFIT_ASIMOV_projection_ktcgkb_combWithHbb_Jul10_couplingModel_reweighted_couplingdependentBRs.MultiDimFit.mH125.root'
# # differentialutils.run_fastscan_scan_reused_postfit(
# # config,
# # reusable_postfit
# # )
# # differentialutils.run_postfit_fastscan_scan(config)
# config.minimizer_settings_andrew()
# differentialutils.scan_directly(config)
# else:
# differentialutils.run_postfit_fastscan_scan(config)
differentialutils.run_postfit_fastscan_scan(config)
def scan_top_BRdependent_and_profiledTotalXS(args):
differentialutils.assert_asimov(args)
config.datacard = LatestPaths.ws_combined_Top_couplingDependentBR_profiledTotalXS
config.subDirectory += '_couplingDependentBR_profiledTotalXS'
config.fix_parameter_at_value('kappa_V', 0.999)
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_noTheoryUnc(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
config.datacard = LatestPaths.ws.yukawa.noTheoryUnc
config.tags.append('noTheoryUnc')
differentialutils.run_postfit_fastscan_scan(config)
def scan_yukawa_fitOnlyNormalization(real_args):
args = set_combination_and_asimov(real_args)
config = basic_config(args)
config.datacard = LatestPaths.ws_combined_Yukawa_profiledTotalXS_fitOnlyNormalization
config.tags.append('fitOnlyNormalization')
differentialutils.run_postfit_fastscan_scan(config)
def projection_kbkc_scan_couplingdependentBRs(args):
args = differentialutils.force_asimov(args)
Config = CombineConfigKBKC_couplingdependentBRs_scenario2 if args.scenario2 else CombineConfigKBKC_couplingdependentBRs
config = Config(args)
differentialutils.run_postfit_fastscan_scan(config)
def scan_topctcb(args):
config = basic_config_ctcb(args)
config.datacard = LatestPaths.ws.topctcb.nominal[
differentialutils.get_decay_channel_tag(args)]
differentialutils.run_postfit_fastscan_scan(config)
