def cutFlow(cutTreeMC, model):
    out = []

    cutTreeMC.Attach(count_total)
    if isCmssm(model):
        cutTreeMC.TAttach(count_total, selection)
        out.append(tripleScale(model=switches()["model"], cutTree=cutTreeMC, cut=count_total, label="before"))
    elif isSms(model):
        cutTreeMC.TAttach(count_total, smsFilter)
        out.append(smsOps(model=switches()["model"], cutTree=cutTreeMC, cut=sms_filter, label="before"))
        cutTreeMC.TAttach(smsFilter, selection)

    if switches()["selection"] != "muon":
        cutTreeMC.TAttach(selection, oddMuon)
        cutTreeMC.TAttach(oddMuon, oddElectron)
        cutTreeMC.TAttach(oddElectron, oddPhoton)
        cutTreeMC.TAttach(oddPhoton, numComLeptons)
        cutTreeMC.TAttach(numComLeptons, numComPhotons)
        cutTreeMC.TAttach(numComPhotons, LeadingJetEta)
        cutTreeMC.TAttach(LeadingJetEta, badMuonInJet)
        cutTreeMC.TAttach(badMuonInJet, oddJet)
        cutTreeMC.TAttach(oddJet, numComJetsGeq2)
        cutTreeMC.TAttach(numComJetsGeq2, secondJetET)
        cutTreeMC.TAttach(secondJetET, alphaT)
        cutTreeMC.TAttach(alphaT, deadECAL_MC)
        cutTreeMC.TAttach(deadECAL_MC, MHToverMET)
    else:
        cutTreeMC.TAttach(selection, oddMuon)
        cutTreeMC.TAttach(oddMuon, oddElectron)
        cutTreeMC.TAttach(oddElectron, oddPhoton)
        cutTreeMC.TAttach(oddPhoton, numComElectrons)
        cutTreeMC.TAttach(numComElectrons, numComMuons)
        cutTreeMC.TAttach(numComMuons, minDRMuonJetCut)
        cutTreeMC.TAttach(minDRMuonJetCut, ZMassCut)
        cutTreeMC.TAttach(ZMassCut, CommonMHTovHT)
        cutTreeMC.TAttach(CommonMHTovHT, PFMTCut30)
        cutTreeMC.TAttach(PFMTCut30, numComPhotons)
        cutTreeMC.TAttach(numComPhotons, LeadingJetEta)
        cutTreeMC.TAttach(LeadingJetEta, badMuonInJet)
        cutTreeMC.TAttach(badMuonInJet, oddJet)
        cutTreeMC.TAttach(oddJet, numComJetsGeq2)
        cutTreeMC.TAttach(numComJetsGeq2, secondJetET)
        cutTreeMC.TAttach(secondJetET, alphaT)
        cutTreeMC.TAttach(alphaT, deadECAL_MC)
        cutTreeMC.TAttach(deadECAL_MC, MHToverMET)

    out.append(
        addBinnedStuff(model=switches()["model"], cutTree=cutTreeMC, cut=MHToverMET, htBins=[250, 300, 350, 450])
    )
    out.append(
        addBinnedStuff(
            model=switches()["model"],
            cutTree=cutTreeMC,
            cut=MHToverMET,
            htBins=[275, 325] + [375 + 100 * i for i in range(6)],
        )
    )
    return out
def sample():
    if isCmssm(switches()["model"]):
        return eval("PhysicsProcesses_mSUGRA_%sFall10v1" % switches()["model"].replace("tanB", "tanbeta"))
    elif isSms(switches()["model"]):
        return eval(
            "PhysicsProcesses_Topology%s_38xFall10_spadhi%s"
            % (switches()["model"], "_new" if switches()["model"] == "T2" else "")
        )
    else:
        return None
Exemplo n.º 3
0
def sample():
    if isCmssm(switches()["model"]):
        return eval("PhysicsProcesses_mSUGRA_%sFall10v1" %
                    switches()["model"].replace("tanB", "tanbeta"))
    elif isSms(switches()["model"]):
        return eval("PhysicsProcesses_Topology%s_38xFall10_spadhi%s" %
                    (switches()["model"],
                     "_new" if switches()["model"] == "T2" else ""))
    else:
        return None
Exemplo n.º 4
0
def checkSwitches(d) :
  assert d["model"] in ["T1","T2","tanB3","tanB10","tanB50"]
  assert d["jes"] in ["","+ve","-ve"]
  assert isCmssm(d["model"]) or isSms(d["model"])
Exemplo n.º 5
0
def cutFlow(cutTreeMC, model):
    out = []

    cutTreeMC.Attach(count_total)
    if isCmssm(model):
        cutTreeMC.TAttach(count_total, selection)
        out.append(
            tripleScale(model=switches()["model"],
                        cutTree=cutTreeMC,
                        cut=count_total,
                        label="before"))
    elif isSms(model):
        cutTreeMC.TAttach(count_total, smsFilter)
        out.append(
            smsOps(model=switches()["model"],
                   cutTree=cutTreeMC,
                   cut=sms_filter,
                   label="before"))
        cutTreeMC.TAttach(smsFilter, selection)

    if switches()["selection"] != "muon":
        cutTreeMC.TAttach(selection, oddMuon)
        cutTreeMC.TAttach(oddMuon, oddElectron)
        cutTreeMC.TAttach(oddElectron, oddPhoton)
        cutTreeMC.TAttach(oddPhoton, numComLeptons)
        cutTreeMC.TAttach(numComLeptons, numComPhotons)
        cutTreeMC.TAttach(numComPhotons, LeadingJetEta)
        cutTreeMC.TAttach(LeadingJetEta, badMuonInJet)
        cutTreeMC.TAttach(badMuonInJet, oddJet)
        cutTreeMC.TAttach(oddJet, numComJetsGeq2)
        cutTreeMC.TAttach(numComJetsGeq2, secondJetET)
        cutTreeMC.TAttach(secondJetET, alphaT)
        cutTreeMC.TAttach(alphaT, deadECAL_MC)
        cutTreeMC.TAttach(deadECAL_MC, MHToverMET)
    else:
        cutTreeMC.TAttach(selection, oddMuon)
        cutTreeMC.TAttach(oddMuon, oddElectron)
        cutTreeMC.TAttach(oddElectron, oddPhoton)
        cutTreeMC.TAttach(oddPhoton, numComElectrons)
        cutTreeMC.TAttach(numComElectrons, numComMuons)
        cutTreeMC.TAttach(numComMuons, minDRMuonJetCut)
        cutTreeMC.TAttach(minDRMuonJetCut, ZMassCut)
        cutTreeMC.TAttach(ZMassCut, CommonMHTovHT)
        cutTreeMC.TAttach(CommonMHTovHT, PFMTCut30)
        cutTreeMC.TAttach(PFMTCut30, numComPhotons)
        cutTreeMC.TAttach(numComPhotons, LeadingJetEta)
        cutTreeMC.TAttach(LeadingJetEta, badMuonInJet)
        cutTreeMC.TAttach(badMuonInJet, oddJet)
        cutTreeMC.TAttach(oddJet, numComJetsGeq2)
        cutTreeMC.TAttach(numComJetsGeq2, secondJetET)
        cutTreeMC.TAttach(secondJetET, alphaT)
        cutTreeMC.TAttach(alphaT, deadECAL_MC)
        cutTreeMC.TAttach(deadECAL_MC, MHToverMET)

    out.append(
        addBinnedStuff(model=switches()["model"],
                       cutTree=cutTreeMC,
                       cut=MHToverMET,
                       htBins=[250, 300, 350, 450]))
    out.append(
        addBinnedStuff(model=switches()["model"],
                       cutTree=cutTreeMC,
                       cut=MHToverMET,
                       htBins=[275, 325] + [375 + 100 * i for i in range(6)]))
    return out
Exemplo n.º 6
0
def checkSwitches(d):
    assert d["model"] in ["T1", "T2", "tanB3", "tanB10", "tanB50"]
    assert d["jes"] in ["", "+ve", "-ve"]
    assert isCmssm(d["model"]) or isSms(d["model"])