本文整理汇总了Python中icf.core.PSet类的典型用法代码示例。如果您正苦于以下问题：Python PSet类的具体用法？Python PSet怎么用？Python PSet使用的例子？那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了PSet类的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: makePlotOp
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示例2: evList_to_pset
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()示例3: AddBinedHist
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示例4: setup_light_tree
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示例5: AddHistPairWithL1
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示例6: PSet
# Configuration
# ====================================
tag = "DeltaInv"
mht_cut = 100.0
gen_match = False
lepton_type = "Electron"
my_ntuple_path = "/home/alex/Research/Samples/ICF"
id_working_point = eid.eff_70
id_working_point_noIso = eid.eff_80
sig_id = PSet(
Cuts=id_working_point_noIso,
useIso=False,
Isolation=True,
HoverE=True,
DeltaAtVtx=True,
SigmaIEtaIEta=True,
Conversions=True,
)
bg_id = PSet(
Cuts=id_working_point_noIso,
useIso=False,
Isolation=True,
HoverE=True,
DeltaAtVtx=False,
SigmaIEtaIEta=True,
Conversions=True,
)
示例7: OP_MHToverMET
# Analyses
MHT_METCut = OP_MHToverMET(1.25,50.)
# AK5 Calo
json_ouput = JSONOutput("filtered")
json_lost = JSONOutput("lost")
alphaT = OP_CommonAlphaTCut(0.53)
# json = JSONFilter("Json Mask", json_to_pset("/home/hep/db1110/public_html/Golden2011.json"))
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() )示例8: PSet
if exc.errno == errno.EEXIST:
pass
else: raise
#============= Triggers ==============
mutriggerps = PSet(
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"示例9: OP_UnCorLeadJetCut
unCorLeadJetCut = OP_UnCorLeadJetCut(40.)
# -----------------------------------------------------------------------------
# Common plots
plots_comjet = OP_ComJetPlots("CommonJetPlots",6)
plots_objkin = OP_ObjKinPlots("ObjectKinePlots",100,6)
plots_common = OP_CommonPlots("CommonPlots")
plots_kinsuite = OP_KinSuiteComPlot("KinSuitePlots",6,2)
pset1 = PSet(
DirName = "MT2",
MinObjects = 2,
MaxObjects = 10,
Mt2DiJet = False,#MT2 plots of Dijet system NOT USED HERE
Mt2Multi_Jet = False,#As above for multiJet
Mt2Compare = False,#Compares the above methods
MHT = True, #Plots MHT Variables
DalitzMT2 = False, # A testing/idea for putting MT2 in with Dalitz plots, also False
Development = True, #Development area, here is where the JetMET PAS plots have been developed, Needs tidying, will be put in to a File of its own, for the moment everything you want will be in here
)
pset2 = PSet(
DirName = "MT2after",
MinObjects = 2,
MaxObjects = 10,
Mt2DiJet = False,#MT2 plots of Dijet system NOT USED HERE
Mt2Multi_Jet = False,#As above for multiJet
Mt2Compare = False,#Compares the above methods
MHT = True, #Plots MHT Variables
DalitzMT2 = False, # A testing/idea for putting MT2 in with Dalitz plots, also False示例10: __init__
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
#.........这里部分代码省略.........
示例11: PSet
# -----------------------------------------------------------------------------
# Necessary includes
import setupSUSY
from libFrameworkSUSY import *
from libHadronic import *
from icf.core import PSet,Analysis
from icf.config import defaultConfig
from copy import deepcopy
from icf.JetCorrections import *
mSuGraScanPlots_beforeAll_pset = PSet(DirectoryName = "mSuGraScan_beforeAll")
#mSuGra scan plots
mSuGraScanPlots_before_150_pset = PSet(DirectoryName = "mSuGraScan_before_150")
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,示例12: Analysis
conf.Common.Muons.PtCut = 15.0
conf.Common.Muons.TrkIsoCut = -1.
conf.Common.Muons.CombIsoCut = 0.1
conf.Common.Electrons.PtCut = 15.0
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,示例13: PSet
Pythia8=[QCD_Pt_0to15_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_15to20_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_20to30_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_30to50_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_50to80_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_80to120_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_120to170_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_170to230_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_230to300_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_300to380_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_380to470_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_470to600_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_600to800_7TeV_pythia8_Summer10_START36_V10_S09_v1,QCD_Pt_800to1000_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_1000to1400_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_1400to1800_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_1800to2200_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_2200to2600_7TeV_pythia8_Summer10_START36_V10_S09_v2,QCD_Pt_3000to3500_7TeV_pythia8_Summer10_START36_V10_S09_v2]
# CMSSW_3_8_4_patch3 V14-00-02 samples
from montecarlo.QCD_Pythia6_384patch3_V14_00_02_ALL import *
from montecarlo.QCD_Pythia8_384patch3_V14_00_02_ALL import *
from data.Jet_15pb_WithTP_json221010 import *
AllMC = QCD_Pythia6_384patch3_V14_00_02_ALL+QCD_Pythia8_384patch3_V14_00_02_ALL+MC+Pythia8
# -----------------------------------------------------------------------------
#Plot the common plots!
pset1 = PSet(
DirName = "AllNoTrigger",
MinObjects = 0,
MaxObjects = 15,
StandardPlots = True,
)
pset2 = PSet(
DirName = "AllWithTrigger",
MinObjects = 0,
MaxObjects = 15,
StandardPlots = True,
)
pset3 = PSet(
DirName = "AllCutsNoTrigger",
MinObjects = 0,
MaxObjects = 15,
StandardPlots = True,
)
pset4 = PSet(示例14: OP_OddPhoton
oddPhoton = OP_OddPhoton()
numComElectrons = OP_NumComElectrons("<=",0)
numComMuons = OP_NumComMuons("<=",0)
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)示例15: elif
elif ( VetoesOption == "NoEGamma" ) :
cut_flow.TAttach(selection_duplicate,numComMuons)
cut_flow.TAttach(numComMuons,oddMuon)
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" ,