示例#1
0
def analyze(event, jet_trigger_filter, **kw):
    for jet in event.jets():
        if accept_jet(event, jet) and jet_trigger_filter(event, jet):
            for track in filter(event.charge_filter, event.matchedPairs()):
                if accept_track(track) and abs(track.DeltaPhi(jet))>2.0:
                    z = track.ptPr()/shifted(event.year, jet.Pt())
                    yield (z, mcasym.num('ZERO', event))
示例#2
0
def analyze(event, jet_trigger_filter, **kw):
    for jet in event.jets():
        if accept_jet(event, jet) and jet_trigger_filter(event, jet):
            for track in filter(event.charge_filter, event.matchedPairs()):
                if accept_track(track) and abs(track.DeltaPhi(jet)) > 2.0:
                    true_pt = shifted(event.year, jet.Pt())
                    z = track.ptPr() / true_pt
                    weight = minbias_jet_weight_true(true_pt)
                    yield (z, weight * mcasym.num('MIN', event))
示例#3
0
文件: ZERO.py 项目: kocolosk/analysis
def analyze(event, jet_trigger_filter, **kw):
    for jet in event.jets():
        if accept_jet(event, jet) and jet_trigger_filter(event, jet):
            for track in filter(event.charge_filter, event.matchedPairs()):
                if accept_track(track) and abs(track.DeltaPhi(jet))>2.0:
                    true_pt = shifted(event.year, jet.Pt())
                    z = track.ptPr()/true_pt
                    weight = minbias_jet_weight_true(true_pt)
                    yield (z, weight*mcasym.num('ZERO', event))
示例#4
0
def analyze(event, **kw):
    for track in event.matchedPairs():
        if event.charge_filter(track) and accept_track(track):
            yield (track.pt(), mcasym.num('P070', event))
示例#5
0
def analyze(event, **kw):
    for track in event.pythiaRecord():
        if event.charge_filter(track) and accept_track(track):
            yield (track.vec.Pt(), mcasym.num('P045', event))
示例#6
0
def analyze(event, **kw):
    for track in event.matchedPairs():
        if event.charge_filter(track) and accept_track(track):
            yield (track.pt(), mcasym.num("P045", event))
示例#7
0
文件: M045.py 项目: kocolosk/analysis
def analyze(event, **kw):
    for track in event.pythiaRecord():
        if event.charge_filter(track) and accept_track(track):
            yield (track.vec.Pt(), mcasym.num('M045', event))