本文整理汇总了Python中icf.core.PSet.ps方法的典型用法代码示例。如果您正苦于以下问题：Python PSet.ps方法的具体用法？Python PSet.ps怎么用？Python PSet.ps使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类icf.core.PSet的用法示例。
在下文中一共展示了PSet.ps方法的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: AddHistPairWithL1
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def AddHistPairWithL1(cutTree = None,cut = None, RefTrig = None, TestTrig = None, L1ListRef = None,L1ListTest = None):
"""docstring for AddBinedHist"""
out = []
refTrigs = None
if "Mu" in RefTrig:
refTrigs = [TestTrig,RefTrig]
else: refTrigs = [TestTrig]
refPlots = PL_TriggerTurnOns( PSet(DirName = RefTrig+"_For_"+TestTrig,MinObjects =0 ,
MaxObjects = 15,Plots = True, ReWeight = True if "Mu" not in RefTrig else False,
TriggerReWeight = refTrigs, Verbose = False,
ReWeightL1 = False, L1TriggerReWeight = [L1ListRef]).ps())
testTrigPlots = PL_TriggerTurnOns( PSet(DirName = TestTrig+"_From_"+RefTrig,MinObjects = 0,
MaxObjects = 15,Plots = True, ReWeight = True,
TriggerReWeight = [TestTrig], Verbose = False,
ReWeightL1 = True, L1TriggerReWeight = [L1ListTest]).ps())
refTrigPS = PSet(Verbose = False,UsePreScaledTriggers = True,Triggers = [] )
refTrigPS.Triggers = [RefTrig]
refTrigOP = OP_MultiTrigger( refTrigPS.ps() )
testTrigPS = PSet(Verbose = False,UsePreScaledTriggers = True,Triggers = [] )
testTrigPS.Triggers = [TestTrig]
print "RefTrig = %s, testTrig = %s"%(refTrigPS.Triggers[0],testTrigPS.Triggers[0])
testTrigOP = OP_MultiTrigger( testTrigPS.ps() )
cutTree.TAttach(cut,refTrigOP)
cutTree.TAttach(refTrigOP,refPlots)
cutTree.TAttach(refTrigOP,testTrigOP)
cutTree.TAttach(testTrigOP,testTrigPlots)
out.append(refTrigOP)
out.append(refPlots)
out.append(testTrigPlots)
out.append(testTrigOP)
return out
pass示例2: makePlotOp
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def makePlotOp(OP = (), cutTree = None, cut = None, label = ""):
"""docstring for makePlotOp"""
out = []
if OP[1] != None:
plotpset = deepcopy(OP[1])
plotpset.DirName = label
print plotpset.DirName
op = eval(OP[0]+"(plotpset.ps())")
else:
op = eval(OP[0])
out.append(op)
cutTree.TAttach(cut,op)
alpha = OP_CommonAlphaTCut(0.55)
dump = EventDump()
skim_ps=PSet(
SkimName = "myskim",
DropBranches = False,
Branches = [
" keep * "
]
)
skim = SkimOp(skim_ps.ps())
# out.append(skim)
# out.append(skim_ps)
cutTree.TAttach(cut,alpha)
#cutTree.TAttach(cut,dump)
#cutTree.TAttach(cut,skim)
#cutTree.TAttach(alpha,dump)
#cutTree.TAttach(alpha,skim)
out.append(skim)
out.append(alpha)
out.append(dump)
return out
pass示例3: evList_to_pset
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def evList_to_pset(fname, debug=False):
""" Takes an input event list file and converts to a PSet suitable for use in analysis code."""
"""NOTE: made for txt file format 'Run:LS:Event' """
import sys
from icf.core import PSet
with open(fname, "r") as f:
p = PSet()
runList = []
lumiSList = []
evList = []
for line in f:
lineTmp = line.split(":")
#fill each list
runList.append(int(lineTmp[0]))
lumiSList.append(int(lineTmp[1]))
evList.append(int(lineTmp[2]))
#create new entry in the PSet for each list
p._quiet_set("Run", runList)
p._quiet_set("Lumi", lumiSList)
p._quiet_set("Event", evList)
return p.ps()示例4: AddBinedHist
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def AddBinedHist(cutTree = None, OP = (), cut = None, htBins = [],TriggerDict = None,lab = ""):
"""docstring for AddBinedHist"""
out = []
if TriggerDict is not None:
for lower,upper in zip(htBins,htBins[1:]+[None]):
# print "Lower , Upper =", lower , upper
if int(lower) == 325 and upper is None: continue
if int(lower) == 375 and upper is None: continue
if int(lower) == 475 and upper is None: continue
if int(lower) == 675 and upper is None: continue
# print "continue should have happened now"
lowerCut = eval("RECO_CommonHTCut(%d)"%lower)
triggerps = PSet(Verbose = False,
UsePreScaledTriggers = False,
Triggers = [])
triggerps.Triggers = TriggerDict["%d%s"%(lower, "_%d"%upper if upper else "")]
Trigger = OP_MultiTrigger( triggerps.ps() )
out.append(triggerps)
out.append(Trigger)
out.append(lowerCut)
cutTree.TAttach(cut,Trigger)
cutTree.TAttach(Trigger,lowerCut)
if upper:
upperCut = eval("RECO_CommonHTLessThanCut(%d)"%upper)
out.append(upperCut)
cutTree.TAttach(lowerCut,upperCut)
pOps = makePlotOp(cutTree = cutTree, OP = OP, cut = upperCut if upper else lowerCut, label = "%s%d%s"%(lab,lower, "_%d"%upper if upper else ""))
out.append(pOps)
else:
for lower,upper in zip(htBins,htBins[1:]+[None]):
# print "Lower , Upper =", lower , upper
if int(lower) == 325 and upper is None: continue
if int(lower) == 375 and upper is None: continue
if int(lower) == 475 and upper is None: continue
if int(lower) == 675 and upper is None: continue
# print "continue should have happened now"
lowerCut = eval("RECO_CommonHTCut(%d)"%lower)
out.append(lowerCut)
cutTree.TAttach(cut,lowerCut)
if upper:
upperCut = eval("RECO_CommonHTLessThanCut(%d)"%upper)
out.append(upperCut)
cutTree.TAttach(lowerCut,upperCut)
pOps = makePlotOp(cutTree = cutTree, OP = OP, cut = upperCut if upper else lowerCut, label = "%s%d%s"%(lab,lower, "_%d"%upper if upper else ""))
out.append(pOps)
return out
pass示例5: setup_light_tree
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def setup_light_tree(baton, main, data, dont_use_gen, trigger_info):
import wpol.electron_id as electron_id
import libAlex as alex
tree_cfg = PSet(
DoGen = not (data or dont_use_gen),
TriggerInfo = trigger_info,
ElectronCuts80 = electron_id.eff_80,
ElectronCuts70 = electron_id.eff_70,
ElectronCuts60 = electron_id.eff_60,
Isolation = True,
HoverE = True,
DeltaEtaAtVtx = True,
DeltaPhiAtVtx = True,
SigmaIEtaIEta = True,
Conversions = True,
ConversionsExtra = True,
SupressErrors = True
)
baton.tree2 = alex.AlexTree("light_tree_lt4jets", tree_cfg.ps())
main.tree.TAttach(main.GetOp("lt4_jets"), baton.tree2)
return baton示例6: EventDump
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
evDump = EventDump()
# htTriggerEmu = OP_TriggerHT_Emu(250.,40.)
cutTreeData = Tree("Data")
out = []
skim_ps=PSet(
SkimName = "myskim",
DropBranches = False,
Branches = [
" keep * "
]
)
badEvents =runLumiCutter(PSet(Run= [166033,162811, 163233, 163374, 163589, 163759, 163817, 165121, 165364, 165415, 165467, 165472, 165506, 165567, 165570, 165993, 166034, 166049, 166380, 166514, 166565, 166701, 166842, 166864, 166960, 167102, 167282, 167284, 167675, 167807, 167830, 167830, 167898, 167913, 167913, 170876, 170876, 171050, 171050, 171156, 171178, 171178, 171369, 171484, 171484, 171578, 171578, 171812, 171876, 171897, 172033, 172791, 172791, 172799, 172802, 172822, 172822, 172868, 172868, 172949, 172952, 173198, 173241, 173241, 173243, 173380, 173389, 173439, 173657, 173657, 173692, 173692, 173692, 175975, 175990, 176023, 176023, 176201, 176286, 176286, 176304, 176309, 176309, 176547, 176548, 176701, 176702, 176765, 176771, 176771, 176771, 176795, 176796, 176796, 176799, 176844, 176886, 176928, 176933, 176982, 177053, 177074, 177096, 177183, 177183, 177201, 177201, 177222, 177730, 177782, 177782, 177782, 177788, 177875, 178098, 178420, 178421, 178421, 178479, 178479, 178703, 178786, 178786, 178803, 178803, 178920, 178970, 179411, 179411, 179411, 179411, 179434, 179497, 179547, 179889, 180072, 180072, 180076, 180076, 180241, 180241, 180241],Lumi =[0,4, 274, 671, 46, 182, 118, 242, 1132, 812, 113, 493, 57, 540, 696, 778, 203, 57, 637, 12, 12, 79, 7, 318, 60, 6, 38, 891, 517, 1406, 171, 458, 1334, 12, 14, 180, 415, 53, 92, 211, 149, 731, 61, 358, 79, 347, 696, 323, 382, 295, 256, 664, 67, 203, 304, 1922, 616, 1755, 553, 933, 680, 665, 16, 759, 13, 199, 283, 307, 68, 90, 2250, 89, 927, 179, 56, 60, 65, 135, 166, 62, 186, 484, 696, 28, 624, 8, 399, 110, 148, 64, 78, 42, 10, 45, 153, 61, 378, 7, 82, 4, 530, 488, 129, 9, 98, 258, 63, 59, 213, 160, 73, 79, 22, 148, 183, 53, 111, 230, 189, 68, 177, 194, 67, 152, 55, 207, 157, 10, 23, 6, 96, 83, 192, 233, 169, 61, 80, 201, 229, 129, 65, 99]).ps())
skim = SkimOp(skim_ps.ps())
# cutTreeData.Attach(json)
triggers = ["HLT_HT250_AlphaT0p55_v1","HLT_HT250_AlphaT0p55_v2","HLT_HT250_AlphaT0p53_v2","HLT_HT250_AlphaT0p53_v3","HLT_HT250_AlphaT0p53_v4","HLT_HT250_AlphaT0p53_v5","HLT_HT250_AlphaT0p53_v6","HLT_HT250_AlphaT0p55_v2","HLT_HT250_AlphaT0p58_v3","HLT_HT300_AlphaT0p52_v1","HLT_HT300_AlphaT0p52_v2","HLT_HT300_AlphaT0p52_v3","HLT_HT300_AlphaT0p52_v4","HLT_HT300_AlphaT0p52_v5","HLT_HT300_AlphaT0p53_v5","HLT_HT300_AlphaT0p53_v6","HLT_HT300_AlphaT0p53_v6","HLT_HT300_AlphaT0p54_v5","HLT_HT350_AlphaT0p51_v1","HLT_HT350_AlphaT0p51_v2","HLT_HT350_AlphaT0p51_v3","HLT_HT350_AlphaT0p51_v4","HLT_HT350_AlphaT0p51_v5","HLT_HT350_AlphaT0p52_v1","HLT_HT350_AlphaT0p52_v2","HLT_HT350_AlphaT0p52_v2","HLT_HT350_AlphaT0p53_v10","HLT_HT400_AlphaT0p51_v1","HLT_HT400_AlphaT0p51_v2","HLT_HT400_AlphaT0p51_v3","HLT_HT400_AlphaT0p51_v4","HLT_HT400_AlphaT0p51_v5","HLT_HT400_AlphaT0p51_v6","HLT_HT400_AlphaT0p51_v7","HLT_HT400_AlphaT0p51_v7","HLT_HT400_AlphaT0p51_v10"]
triggers = ["HLT_HT250_v5"]
CheckPreOverLapsOp = CheckPreOverLaps()
# for trig in triggers:
# op = JSONOutput("%s"%(trig))
# trigCut = CheckTrigExists( PSet(TrigExistsList = ["%s"%(trig)]).ps() )
# trigCut = OP_MultiTrigger( PSet(Verbose = True,UsePreScaledTriggers = True,Triggers = ["%s"%(trig)] ).ps() )
# cutTreeData.TAttach(json,trigCut)
# cutTreeData.TAttach(trigCut,skim)
# cutTreeData.TAttach(trigCut,op)
# out.append(trigCut)
# out.append(op)
# cutTreeData.TAttach(json,badEvents)
# cutTreeData.FAttach(badEvents,json_lost)示例7: Analysis
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
conf.Common.Electrons.EtaCut = 2.4
conf.Common.Photons.EtCut=30.
conf.Common.Electrons.ApplyID=False
# Create the analysis
a = Analysis("KinCorrection")
tree = Tree("Main")
Trigger=PSet(
McAnal =False,
MSugraScan = 0.,
TriggerBits=PSet(
bit1=""
)
)
ssdlTrigg=SSDLTrigger("Trigger",Trigger.ps())
ZeroMuons = OP_NumComMuons("==", 0)
KinCorrParPlus = PSet(
zMassMin = 71.,
zMassMax = 111.,
BarrelEtaMax = 1.4442,
EndCapEtaMin = 1.56,
EndCapEtaMax = 2.5,
c_ErNu_pt = 0.,
MinElecRescaledPt=0.,
wID=+24., ##24 W+, -24 W-示例8: OP_MultiTrigger
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
Verbose = False,
UsePreScaledTriggers = False,
Triggers = [
"HLT_HT250_AlphaT0p55_v1",
"HLT_HT250_AlphaT0p55_v2",
"HLT_HT250_AlphaT0p53_v2",
"HLT_HT250_AlphaT0p53_v3",
"HLT_HT250_AlphaT0p53_v4",
"HLT_HT250_AlphaT0p53_v5",
"HLT_HT250_AlphaT0p55_v2",
"HLT_HT250_AlphaT0p58_v3",
"HLT_HT250_AlphaT0p60_v3",
]
)
MuTrigger = OP_MultiTrigger(mutriggerps.ps())
#===========================================
#============ Plotting Ops =======================
generic_pset = PSet(
DirName = "MuonAnalysis",
FolderName = "Muon275"
)
skim_ps=PSet(
SkimName = "myskim",
DropBranches = False,
Branches = [
" keep * "
]
)示例9: OP_NumComPhotons
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
numComPhotons = OP_NumComPhotons("<=",0)
# -----------------------------------------------------------------------------
# Test module
from allhadronic.algorithms_cfi import *
ps = PSet(
Verbose = False,
AlphaT = [0.55],
DirName = "Test",
nMax = 6,
Algorithms = AlgorithmsPSet,
)
test = Scaling(ps.ps())
# -----------------------------------------------------------------------------
# Cut flow
cut_flow = Tree("MC")
cut_flow.Attach(selection)
cut_flow.TAttach(selection,numComPhotons)
cut_flow.TAttach(numComPhotons,oddPhoton)
cut_flow.TAttach(oddPhoton,numComElectrons)
cut_flow.TAttach(numComElectrons,oddElectron)
cut_flow.TAttach(oddElectron,numComMuons)
cut_flow.TAttach(numComMuons,oddMuon)
cut_flow.TAttach(oddMuon,test)
# -----------------------------------------------------------------------------示例10: PSet
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
mSuGraScanPlots_before_250_pset = PSet(DirectoryName = "mSuGraScan_before_250")
mSuGraScanPlots_before_350_pset = PSet(DirectoryName = "mSuGraScan_before_350")
mSuGraScanPlots_150_pset = PSet(DirectoryName = "mSuGraScan_150")
mSuGraScanPlots_250_pset = PSet(DirectoryName = "mSuGraScan_250")
mSuGraScanPlots_350_pset = PSet(DirectoryName = "mSuGraScan_350")
#Bkgd Estimation Plots
TanjaBkgdEstPlots_g1jets = PSet(DirectoryName = "BkgdEstPlots_g1jets",
MinObjects=2,
MaxObjects=9,
m0_felder=0,
m12_felder=0
)
plots_bkgdEstPlots_g1jets = OP_BkgdEstPlottingOps( TanjaBkgdEstPlots_g1jets.ps() )
TanjaBkgdEstPlots_2jets = PSet(DirectoryName = "BkgdEstPlots_2jets",
MinObjects=2,
MaxObjects=2,
m0_felder=0,
m12_felder=0
)
plots_bkgdEstPlots_2jets = OP_BkgdEstPlottingOps( TanjaBkgdEstPlots_2jets.ps() )
TanjaBkgdEstPlots_g2jets = PSet(DirectoryName = "BkgdEstPlots_g2jets",
MinObjects=3,
MaxObjects=7,
m0_felder=0,
m12_felder=0
)示例11: __init__
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
def __init__(self,
data,
loose_id,
tight_id,
trigger_bits = (
"HLT_Ele10_LW_L1R", # run < 140041
"HLT_Ele15_SW_L1R", # run <= 143962
"HLT_Ele15_SW_CaloEleId_L1R", # run <= 146427
"HLT_Ele17_SW_CaloEleId_L1R", # run <= 147116
"HLT_Ele17_SW_TightEleId_L1R", # run <= 148818
"HLT_Ele22_SW_TighterEleId_L1R_v2", # run <= 149180
"HLT_Ele22_SW_TighterEleId_L1R_v3" # run >= 149181
),
pf_mht = True,
mht_cuts = [50, 75, 100 ],
met_cuts = [10, 20, 30, 40],
mt_cuts = [30, 50, 70],
lepton_type = "Electron",
ele_jet_dr = 0.3,
lep_ptcut = 25.,
name = "eWPol",
ignore_gen = False,
event_no_dump = True
):
self.name = name
self.jet_dr = wpol.EleJetDRFilter(ele_jet_dr)
self.loose_ele_id = wpol.CustomEleId(loose_id.ps())
self.tight_ele_id = wpol.CustomEleId(tight_id.ps())
self.lep_pt_cut = lep_ptcut
self.pol_plots_cfg = PSet(
DoGen = not (data or ignore_gen),
LeptonType = lepton_type,
METThreshold = 40.,
NLeptons = 1
)
self.ops = []
# Basic event selection
triggers = PSet(Triggers=trigger_bits)
self.ops.append(("basic_plots_pretrig",
wpol.eWPolBasicPlots("eWPolBasicPlots_PreTrigger",
self.pol_plots_cfg.ps())))
self.ops.append(("trigger", fwk.OP_MultiTrigger(
triggers.ps()
)))
self.ops.append(("basic_plots_posttrig",
wpol.eWPolBasicPlots("eWPolBasicPlots_PostTrigger",
self.pol_plots_cfg.ps())))
self.ops.append(("good_events", fwk.OP_GoodEventSelection()))
# Basic object cuts
self.ops.append(("datamc_one_ele", wpol.DataMCPlots("RECO_DataMCPlots_OneEle")))
self.ops.append(("one_ele", fwk.OP_NumComElectrons("==", 1)))
self.ops.append(("basic_plots_gteqoneele",
wpol.eWPolBasicPlots("eWPolBasicPlots_GTEQOneEle",
self.pol_plots_cfg.ps())))
self.ops.append(("good_ele", wpol.ApplyLeptonCut(lepton_type,
lep_ptcut,
self.tight_ele_id,
1,
True)))
self.ops.append(("datamc_one_wp70_ele", wpol.DataMCPlots("RECO_DataMCPlots_OneWP70Ele")))
self.ops.append(("three_charge", wpol.ChargeAgreement()))
self.ops.append(("basic_plots_oneele",
wpol.eWPolBasicPlots("eWPolBasicPlots_OneEle",
self.pol_plots_cfg.ps())))
# self.ops.append(("z_mass", wpol.MassWindowCut(lepton_type,
# PSet( MassLower = 76.,
# MassUpper = 106.,
# MatchCharge = True,
# NetCharge = 0).ps() )))
self.ops.append(("zero_mu", fwk.OP_NumComMuons("==", 0)))
self.ops.append(("lt4_jets", fwk.OP_NumComJets("<", 4)))
self.ops.append(("pol_plots_premht", wpol.RECO_ElPolPlots(
"RECO_ElPolPlots_PreMHT",
self.pol_plots_cfg.ps())))
# Pre MHT Cut
self.ops.append(("basic_plots_premht",
wpol.eWPolBasicPlots("eWPolBasicPlots_PreMHT",
self.pol_plots_cfg.ps())))
# self.ops.append(("datamc_plots_premht", wpol.DataMCPlots(
# "RECO_DataMCPlots_PreMHT2")))
# MHT Cuts and plots
self.mht_cuts = {}
self.pol_plots = {}
self.datamc_plots = {}
for mht in mht_cuts:
if pf_mht:
self.mht_cuts[mht] = wpol.OP_PFMHTCut(mht)
else:
self.mht_cuts[mht] = fwk.RECO_CommonMHTCut(mht)
self.pol_plots["mht%d" % mht] = wpol.RECO_ElPolPlots(
"RECO_ElPolPlots_PostMHT%d" % mht,
self.pol_plots_cfg.ps())
self.datamc_plots["mht%d" % mht] = wpol.DataMCPlots(
"RECO_DataMCPlots_PostMHT%d" % mht)
# MET Cuts and plots
#.........这里部分代码省略.........
示例12: PSet
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
StandardPlots = True,
)
pset5 = PSet(
DirName = "AllCutsTriggerFail",
MinObjects = 0,
MaxObjects = 15,
StandardPlots = True,
)
pset6 = PSet(
DirName = "AllTriggerFail",
MinObjects = 0,
MaxObjects = 15,
StandardPlots = True,
)
AllNoTrigger = WeeklyUpdatePlots(pset1.ps())
AllWithTrigger =WeeklyUpdatePlots(pset2.ps())
AllCutsNoTrigger = WeeklyUpdatePlots(pset3.ps())
AllCutsAfterTrigger =WeeklyUpdatePlots( pset4.ps() )
AllCutsTriggerFail = WeeklyUpdatePlots(pset5.ps())
AllTriggerFail = WeeklyUpdatePlots(pset6.ps())
alphaTnumbers150 = TriggerEffPlots( PSet(DirName = "alphaTnumbers150",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers200 = TriggerEffPlots( PSet(DirName = "alphaTnumbers200",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers250 = TriggerEffPlots( PSet(DirName = "alphaTnumbers250",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers300 = TriggerEffPlots( PSet(DirName = "alphaTnumbers300",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers350 = TriggerEffPlots( PSet(DirName = "alphaTnumbers350",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers400 = TriggerEffPlots( PSet(DirName = "alphaTnumbers400",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers450 = TriggerEffPlots( PSet(DirName = "alphaTnumbers450",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )
alphaTnumbers500 = TriggerEffPlots( PSet(DirName = "alphaTnumbers500",MinObjects =0 ,MaxObjects = 15,EffPlots = True).ps() )示例13: range
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
htBins = [275, 325] + [375+100*i for i in range(6)],
label2 = "btag_%s_%i_NoAlphaT_"%(btags[0],btags[1]),extra = MChiCut))
return out
from ra1objectid.vbtfElectronId_cff import *
from ra1objectid.vbtfMuonId_cff import *
from ra1objectid.ra3PhotonId_cff import *
vbtfElectronIdFilter = Electron_IDFilter( vbtfelectronidWP95ps.ps() )
#vbtfElectronIdFilter = Electron_IDFilter( vbtfelectronidWP90ps.ps() )
ra3PhotonIdFilter = Photon_IDFilter( ra3photonidps.ps() )
muonfilt = CustomVBTFMuID(mu_id.ps()) if switches()["selection"]=="muon" else Muon_IDFilter( vbtfmuonidps.ps() )
# muonfilt = Muon_IDFilter(vbtfmuonidps.ps())
import os
susydir = os.environ['SUSY_WORKING_SW_DIR'] + '/'
if isCmssm(switches()["model"]) :
if switches()["model"] == "tanB10": LOweights = SignalScanLOCrossSectionWeighting(xsToPSet( readLOXS(susydir+"SUSYSignalScan/textfiles/goodModelNames_10_0_1.txt") ).ps())
if switches()["model"] == "tanB40": LOweights = SignalScanLOCrossSectionWeighting(xsToPSet( readLOXS(susydir+"SUSYSignalScan/textfiles/goodModelNames_40_m500_1.txt") ).ps())
def addCutFlowMC(a,cutTreeMC) :
if "tan" in switches()["model"]:a.AddWeightFilter("Weight",LOweights)
if switches()["jes"] != "":
a.AddJetFilter("PreCC",JESUncert)
a.AddMuonFilter("PreCC",muonfilt); print "WARNING: should this be PreCC?"
a.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
a.AddElectronFilter("PreCC",vbtfElectronIdFilter)
a+=cutTreeMC示例14: OP_TriggerCut
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
# ====================================
# trigger = OP_TriggerCut("HLT_Ele15_SW_L1R")
trigger = OP_TriggerCut("HLT_Jet15U")
selection = OP_GoodEventSelection()
HBHE = OP_HadronicHBHEnoiseFilter()
OneElectron = OP_NumComElectrons(">=", 1)
odde = OP_OddElectron()
ZeroMuons = OP_NumComMuons("==", 0)
oddm = OP_OddMuon()
ZeroPhotons = OP_NumComPhotons("==", 0)
LT4Jets = OP_NumComJets(">=", 1)
oddj = OP_OddJet()
BasicPlotsPreMHT = ElectronTree("eTree_PreMHT", pol_plots_cfg.ps())
tree = Tree("Main")
# Main trunk of tree
tree.Attach(trigger)
tree.TAttach(trigger, selection)
tree.TAttach(selection, HBHE)
tree.TAttach(HBHE, OneElectron)
tree.TAttach(OneElectron, ZeroMuons)
tree.TAttach(ZeroMuons, oddj)
tree.TAttach(oddj, BasicPlotsPreMHT)
# tree.TAttach(OneElectron,odde)
# tree.TAttach(odde,ZeroMuons)
# tree.TAttach(ZeroMuons,oddm)示例15: elif
# 需要导入模块: from icf.core import PSet [as 别名]
# 或者: from icf.core.PSet import ps [as 别名]
cut_flow.TAttach(oddMuon,qcdBkgdEst)
elif ( VetoesOption == "NoVetoes" ) :
cut_flow.TAttach(selection_duplicate,qcdBkgdEst)
elif ( VetoesOption == "NoOddVetoes" ) :
cut_flow.TAttach(selection_duplicate,numComPhotons)
cut_flow.TAttach(numComPhotons,numComElectrons)
cut_flow.TAttach(numComElectrons,numComMuons)
cut_flow.TAttach(numComMuons,qcdBkgdEst)
if ( True and FilterOption > -1 ) :
skim_ps = PSet(
SkimName = "Skim",
Branches = ["keep *"],
DropBranches = False,
)
skim = SkimOp( skim_ps.ps() )
cut_flow.TAttach(qcdBkgdEst,skim)
# -----------------------------------------------------------------------------
# dataset
Darren=PSet(
Name="Darren",
Format=("ICF",3),
File=[
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1000_2_5ze.root" ,
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1001_2_QRR.root" ,
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1002_2_56x.root" ,
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1003_2_Z7B.root" ,
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1004_2_jeP.root" ,
"root://xrootd.grid.hep.ph.ic.ac.uk/store/user/henning/ICF/automated/2011_11_11_14_03_49/HT.Run2011A-May10ReReco-v1.AOD/SusyCAF_Tree_1005_2_ObZ.root" ,