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)))
