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)