q = schema.newQuery()
l1_info = lumiQueryAPI.trgBybitnameByrun(q, run, l1_seed)
del q
l1_prescales_seen = set()
for ls, (count, deadtime, bit_number, prescale) in l1_info.iteritems():
if ls in lumis:
l1_prescales_seen.add(prescale)
l1_prescales[(run, ls)] = prescale
print 'Run %i: HLT path, prescales: %s %s L1 seed, prescales: %s %s' % (
run, hlt_path, repr(sorted(hlt_prescales_seen)), l1_seed,
repr(sorted(l1_prescales_seen)))
session.transaction().commit()
del session
del svc
from SUSYBSMAnalysis.Zprime2muAnalysis.tools import to_pickle
to_pickle((l1_prescales, hlt_prescales), path + '.gzpickle')
'''
from SUSYBSMAnalysis.Zprime2muAnalysis.tools import from_pickle
l1, hlt = from_pickle('Mu15.gzpickle')
rls = sorted(set(l1.keys()) & set(hlt.keys()))
max_seen = 0
for rl in rls:
prescale = l1[rl]*hlt[rl]
if prescale > max_seen:
print 'new max prescale seen for', rl, ':', prescale
max_seen = prescale
'''
# Get the L1 prescale too.
q = schema.newQuery()
l1_info = lumiQueryAPI.trgBybitnameByrun(q, run, l1_seed)
del q
l1_prescales_seen = set()
for ls, (count, deadtime, bit_number, prescale) in l1_info.iteritems():
if ls in lumis:
l1_prescales_seen.add(prescale)
l1_prescales[(run,ls)] = prescale
print 'Run %i: HLT path, prescales: %s %s L1 seed, prescales: %s %s' % (run, hlt_path, repr(sorted(hlt_prescales_seen)), l1_seed, repr(sorted(l1_prescales_seen)))
session.transaction().commit()
del session
del svc
from SUSYBSMAnalysis.Zprime2muAnalysis.tools import to_pickle
to_pickle((l1_prescales, hlt_prescales), path + '.gzpickle')
'''
from SUSYBSMAnalysis.Zprime2muAnalysis.tools import from_pickle
l1, hlt = from_pickle('Mu15.gzpickle')
rls = sorted(set(l1.keys()) & set(hlt.keys()))
max_seen = 0
for rl in rls:
prescale = l1[rl]*hlt[rl]
if prescale > max_seen:
print 'new max prescale seen for', rl, ':', prescale
max_seen = prescale
'''
tot_lumi_in_slot += lumi
# if we've accumulated enough lumi, save this slot and make a new one
if tot_lumi_in_slot > lumi_per_slot:
slots.append((tot_lumi_in_slot, curr_slot))
curr_slot = defaultdict(list)
tot_lumi_in_slot = 0
slots.append((tot_lumi_in_slot, curr_slot))
print 'last slot has lumi', tot_lumi_in_slot
return slots, lookup
lumi_per_slot = 491.7
slots, lookup = build_slots(lumi_per_slot)
to_pickle((slots, lookup),
'Run2011MuonsOnly.slots.%iipb.gzpickle' % lumi_per_slot)
#slots, lookup = from_pickle('Run2011MuonsOnly.slots.%iipb.gzpickle' % lumi_per_slot)
num_slots = len(slots)
mass_ranges = [
(60, 120),
(120, 200),
(200, 400),
(120, 600),
]
histos = [
ROOT.TH1F('%s_%s' % m, 'event count for mass range %s-%s GeV' % m,
num_slots, 0, num_slots) for m in mass_ranges
]
tot_lumi_in_slot += lumi
# if we've accumulated enough lumi, save this slot and make a new one
if tot_lumi_in_slot > lumi_per_slot:
slots.append((tot_lumi_in_slot, curr_slot))
curr_slot = defaultdict(list)
tot_lumi_in_slot = 0
slots.append((tot_lumi_in_slot, curr_slot))
print 'last slot has lumi', tot_lumi_in_slot
return slots, lookup
lumi_per_slot = 491.7
slots, lookup = build_slots(lumi_per_slot)
to_pickle((slots, lookup), 'Run2011MuonsOnly.slots.%iipb.gzpickle' % lumi_per_slot)
#slots, lookup = from_pickle('Run2011MuonsOnly.slots.%iipb.gzpickle' % lumi_per_slot)
num_slots = len(slots)
mass_ranges = [
( 60,120),
(120,200),
(200,400),
(120,600),
]
histos = [ROOT.TH1F('%s_%s' % m, 'event count for mass range %s-%s GeV' % m, num_slots, 0, num_slots) for m in mass_ranges]
f = ROOT.TFile('data/Run2011MuonsOnly/ana_datamc_data.root')
t = f.SimpleNtupler.Get('t')
#cut = ROOT.TTreeFormula('cut', 'OurSelNewNoSign && SameSign', t)
