from StopsDilepton.analysis.regions import defaultRegions, reducedRegionsA, reducedRegionsB, reducedRegionsAB, reducedRegionsNew, reducedRegionsC from StopsDilepton.analysis.Cache import Cache setup.verbose = False setup.parameters['metMin'] = args.metMin setup.parameters['metSigMin'] = args.metSigMin if args.regions == "defaultRegions": regions = defaultRegions elif args.regions == "reducedRegionsA": regions = reducedRegionsA elif args.regions == "reducedRegionsB": regions = reducedRegionsB elif args.regions == "reducedRegionsAB": regions = reducedRegionsAB elif args.regions == "reducedRegionsNew": regions = reducedRegionsNew elif args.regions == "reducedRegionsC": regions = reducedRegionsC else: raise Exception("Unknown regions setup") estimators = constructEstimatorList( ['TTJets', 'other', args.estimateDY, args.estimateTTZ]) if args.multiIsoWP != "": multiIsoWPs = ['VL', 'L', 'M', 'T', 'VT'] wpMu, wpEle = args.multiIsoWP.split(',') from StopsDilepton.tools.objectSelection import multiIsoLepString setup.externalCuts.append( multiIsoLepString(wpMu, wpEle, ('l1_index', 'l2_index'))) setup.prefixes.append('multiIso' + args.multiIsoWP.replace(',', '')) if args.relIso04 > 0: setup.externalCuts.append("&&".join([ "LepGood_relIso04[" + ist + "]<" + str(args.relIso04) for ist in ('l1_index', 'l2_index') ])) setup.prefixes.append('relIso04sm' + str(int(100 * args.relIso04)))
argParser.add_argument('--logLevel', action='store', default='INFO', nargs='?', choices=['CRITICAL', 'ERROR', 'WARNING', 'INFO', 'DEBUG', 'TRACE', 'NOTSET'], help="Log level for logging") argParser.add_argument("--regions", action='store', default='defaultRegions', nargs='?', choices=["defaultRegions","reducedRegionsA","reducedRegionsB","reducedRegionsAB","reducedRegionsNew"], help="which regions setup?") argParser.add_argument("--signal", action='store', default='T2tt', nargs='?', choices=["T2tt","DM"], help="which signal to plot?") argParser.add_argument("--estimateDY", action='store', default='DY', nargs='?', choices=["DY","DY-DD"], help="which DY estimate?") argParser.add_argument("--estimateTTZ", action='store', default='TTZ', nargs='?', choices=["TTZ","TTZ-DD","TTZ-DD-Top16009"], help="which DY estimate?") args = argParser.parse_args() if args.regions == "defaultRegions": regions = defaultRegions elif args.regions == "reducedRegionsA": regions = reducedRegionsA elif args.regions == "reducedRegionsB": regions = reducedRegionsB elif args.regions == "reducedRegionsAB": regions = reducedRegionsAB elif args.regions == "reducedRegionsNew": regions = reducedRegionsNew else: raise Exception("Unknown regions setup") detailedEstimators = constructEstimatorList(['TTJets','other-detailed', args.estimateDY, args.estimateTTZ]) signalSetup = setup.sysClone(parameters={'useTriggers':False}) for estimator in detailedEstimators: estimator.style = styles.fillStyle( getattr( color, estimator.name.split('-')[0] ) ) from StopsDilepton.samples.cmgTuples_FastSimT2tt_mAODv2_25ns_postProcessed import * from StopsDilepton.samples.cmgTuples_FullSimTTbarDM_mAODv2_25ns_postProcessed import * signalEstimators = [ MCBasedEstimate(name=s.name, sample={channel:s for channel in allChannels}, cacheDir=setup.defaultCacheDir() ) for s in ([T2tt_450_0] if args.signal == "T2TT" else [TTbarDMJets_scalar_Mchi1_Mphi100])] for estimator in signalEstimators: estimator.style = styles.lineStyle( getattr(color, estimator.name ), width=2 ) estimators = detailedEstimators + signalEstimators for e in estimators: e.initCache(setup.defaultCacheDir())
argParser.add_argument("--estimateTTZ", action='store', default='TTZ', nargs='?', choices=["TTZ", "TTZ-DD", "TTZ-DD-Top16009"], help="which DY estimate?") args = argParser.parse_args() if args.regions == "defaultRegions": regions = defaultRegions elif args.regions == "reducedRegionsA": regions = reducedRegionsA elif args.regions == "reducedRegionsB": regions = reducedRegionsB elif args.regions == "reducedRegionsAB": regions = reducedRegionsAB elif args.regions == "reducedRegionsNew": regions = reducedRegionsNew else: raise Exception("Unknown regions setup") detailedEstimators = constructEstimatorList( ['TTJets', 'other-detailed', args.estimateDY, args.estimateTTZ]) signalSetup = setup.sysClone(parameters={'useTriggers': False}) for estimator in detailedEstimators: estimator.style = styles.fillStyle( getattr(color, estimator.name.split('-')[0])) from StopsDilepton.samples.cmgTuples_FastSimT2tt_mAODv2_25ns_postProcessed import * from StopsDilepton.samples.cmgTuples_FullSimTTbarDM_mAODv2_25ns_postProcessed import * signalEstimators = [ MCBasedEstimate(name=s.name, sample={channel: s for channel in allChannels}, cacheDir=setup.defaultCacheDir()) for s in ([T2tt_450_0] if args.signal ==
if args.regions == "defaultRegions": regions = defaultRegions elif args.regions == "reducedRegionsA": regions = reducedRegionsA elif args.regions == "reducedRegionsB": regions = reducedRegionsB elif args.regions == "reducedRegionsAB": regions = reducedRegionsAB elif args.regions == "reducedRegionsNew": regions = reducedRegionsNew elif args.regions == "reducedRegionsC": regions = reducedRegionsC else: raise Exception("Unknown regions setup") estimators = constructEstimatorList(["TTJets", "other", args.estimateDY, args.estimateTTZ]) if args.multiIsoWP != "": multiIsoWPs = ["VL", "L", "M", "T", "VT"] wpMu, wpEle = args.multiIsoWP.split(",") from StopsDilepton.tools.objectSelection import multiIsoLepString setup.externalCuts.append(multiIsoLepString(wpMu, wpEle, ("l1_index", "l2_index"))) setup.prefixes.append("multiIso" + args.multiIsoWP.replace(",", "")) if args.relIso04 > 0: setup.externalCuts.append( "&&".join(["LepGood_relIso04[" + ist + "]<" + str(args.relIso04) for ist in ("l1_index", "l2_index")]) ) setup.prefixes.append("relIso04sm" + str(int(100 * args.relIso04)))