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)