if args.noInclusive:
allPhotonRegions = allRegions[args.controlRegion]["regions"]
else:
allPhotonRegions = allRegions[
args.controlRegion]["inclRegion"] + allRegions[
args.controlRegion]["regions"] if photonSelection else allRegions[
args.controlRegion]["regions"]
setup = Setup(year=args.year,
photonSelection=photonSelection
and not "QCD" in args.selectEstimator,
checkOnly=True) #photonselection always false for qcd estimate
# Select estimate
if args.selectEstimator == "Data":
estimate = DataObservation(name="Data", process=setup.processes["Data"])
estimate.isData = True
else:
estimators = EstimatorList(setup, processes=[args.selectEstimator])
estimate = getattr(estimators, args.selectEstimator)
estimate.isData = False
if not estimate:
logger.warning(args.selectEstimator + " not known")
sys.exit(0)
estimate.initCache(setup.defaultCacheDir())
setup3p = setup.sysClone(parameters=parameters3p)
setup3 = setup.sysClone(parameters=parameters3)
setup4p = setup.sysClone(parameters=parameters4p)
allPhotonRegions = allRegions[args.controlRegion]["regions"]
else:
allPhotonRegions = allRegions[
args.controlRegion]["inclRegion"] + allRegions[
args.controlRegion]["regions"] if photonSelection else allRegions[
args.controlRegion]["regions"]
setup = Setup(year=args.year,
photonSelection=photonSelection
and not "QCD" in args.selectEstimator,
checkOnly=True) #photonselection always false for qcd estimate
# Select estimate
if args.selectEstimator == "Data":
estimateFrom = DataObservation(name="Data",
process=setup.processes["Data"])
estimateFrom.isData = True
else:
estimators = EstimatorList(setup, processes=[args.selectEstimator])
estimateFrom = getattr(estimators, args.selectEstimator)
estimateFrom.isData = False
if not estimateFrom:
logger.warning(args.selectEstimator + " not known")
sys.exit(0)
setup = setup.sysClone(parameters=parameters)
def wrapper(arg):
r, channel, setup, addon = arg
logger.info(