from StopsDilepton.analysis.Region import Region from StopsDilepton.analysis.defaultAnalysis import setup #Always taking data lumi setup.lumi = { channel: setup.sample['Data'][channel]['lumi'] for channel in setup.sample['Data'].keys() } setup.loadChains() from StopsDilepton.analysis.DataDrivenDYEstimate import DataDrivenDYEstimate estimateDY = DataDrivenDYEstimate(name='DY-DD', cacheDir=None) regionDY = Region('dl_mt2ll', (140, -1)) for channel in ['MuMu']: res = estimateDY.cachedEstimate(regionDY, channel, setup) print "\n Result in ", channel, " for estimate ", estimateDY.name, regionDY, ":", res #, 'jer',jer, 'jec', jec
def getRegionsFromThresholds(var, vals): return [Region(var, (vals[i], vals[i + 1])) for i in range(len(vals) - 1)] + [Region(var, (vals[-1], -1))]
#!/usr/bin/env python from StopsDilepton.analysis.Region import Region from StopsDilepton.analysis.estimators import setup from StopsDilepton.samples.cmgTuples_Data25ns_mAODv2_postProcessed import * from StopsDilepton.analysis.DataDrivenTTZEstimate import DataDrivenTTZEstimate setup.dataSamples = { "EE": DoubleEG_Run2015D, "MuMu": DoubleMuon_Run2015D, "EMu": MuonEG_Run2015D } setup.verbose = True estimateTTZ = DataDrivenTTZEstimate(name='TTZ-DD', cacheDir=None) regionTTZ = Region('dl_mt2ll', (0, -1)) for channel, sample in setup.dataSamples.iteritems(): setup.lumi[channel] = sample.lumi res = estimateTTZ.cachedEstimate(regionTTZ, channel, setup) print "\n Result in ", channel, " for estimate ", estimateTTZ.name, regionTTZ, ":", res #, 'jer',jer, 'jec', jec
defaultRegions = getRegions3D(mt2llThresholds=[0, 140, 240], mt2blblThresholds=[0, 100, 200], mt2bbThresholds=[70, 170, 270]) reducedRegionsA = getRegions3D(mt2llThresholds=[0, 140, 240], mt2blblThresholds=[0, 100, 200], mt2bbThresholds=[70, 170]) reducedRegionsB = getRegions3D(mt2llThresholds=[0, 140, 240], mt2blblThresholds=[0, 100], mt2bbThresholds=[70, 170, 270]) reducedRegionsAB = getRegions3D(mt2llThresholds=[0, 140, 240], mt2blblThresholds=[0, 100], mt2bbThresholds=[70, 170]) reducedRegionsNew = [ Region("dl_mt2ll", (0, 100)) + Region("dl_mt2bb", (0, -1)) + Region("dl_mt2blbl", (0, -1)) ] # whole mt2ll < 140 as control region reducedRegionsNew += getRegions3D(mt2llThresholds=[100, 140, 240], mt2blblThresholds=[0, 100, 200], mt2bbThresholds=[70, 170]) reducedRegionsC = [ Region("dl_mt2ll", (0, 100)) + Region("dl_mt2bb", (0, -1)) + Region("dl_mt2blbl", (0, -1)) ] # whole mt2ll < 140 as control region reducedRegionsC += getRegions3D(mt2llThresholds=[100, 140, 240], mt2blblThresholds=[0, 150], mt2bbThresholds=[70, 170])