dataTriggerMiter.add( trigger_map['hltdijet20'], bin="100..110" )
dataTriggerMiter.add( trigger_map['hltdijet40'], bin="100..110" )
dataTriggerMiter.add( trigger_map['hltdijet70'], bin="100..110" )
dataTriggerMiter.add( trigger_map['hltdijet20'], bin="110..120" )
dataTriggerMiter.add( trigger_map['hltdijet40'], bin="110..120" )
dataTriggerMiter.add( trigger_map['hltdijet70'], bin="110..120" )
dataTriggerMiter.add( trigger_map['hltdijet110'], bin="120.." )
skiplessDataMiter = Miter()
for opoint in operating_points:
for bin in jet_bins:
for tmap in dataTriggerMiter.iterGrouped( 'bin' ):
if tmap.vals[0][1]['bin'] != bin[1]:
continue
for onetmap in tmap.iterMany():
sample = onetmap[0][5]
trigger = onetmap[0]
step_postfix = "-%s-%s-%s-%s" % (sample,trigger[1],bin[0],opoint)
monitor_node =\
run_monitor_input_helper( g,
jet_pt = bin[1],
tag = opoint,
fileKey = sample,
# trigger_name = trigger[4],
step_postfix = step_postfix,
muon_pt = "6..",
step_postfix = "-%s-%s-%s-djet10mc" % (dataset[1]["dataset"], bin[0], opoint)
currNode = run_monitor_input_helper(
g,
jet_pt=bin[1],
tag=opoint,
fileKey=dataset[1]["dataset"],
trigger_name="HLT_BTagMu_DiJet10U",
step_postfix=step_postfix,
muon_pt="6..",
input_files=dataset[0],
)
dijet10QCDMiter.add(currNode, bin=bin[0], dataset=dataset[1]["dataset"], opoint=opoint)
singleMonitorMiter = Miter()
rootQCDMiter = Miter()
# merge the dijet10 QCD into one guy
for currmerge in dijet10QCDMiter.iterGrouped("bin", "opoint"):
step_postfix = "-%s-%s" % (currmerge.vals[0][1]["bin"], currmerge.vals[0][1]["opoint"])
mergeNode = merge_with_root_qcd_helper(
g,
"root-qcd" + step_postfix,
step_postfix,
inputNodes=currmerge.getValues(),
triggerName="hltjet10u",
lumiMiter=exports["luminosityMiter"],
)
rootQCDMiter.add(mergeNode, bin=currmerge.vals[0][1]["bin"], opoint=currmerge.vals[0][1]["opoint"])
singleMonitorMiter.add(
mergeNode,
trigger="hltdijet10u",
bin=currmerge.vals[0][1]["bin"],
opoint=currmerge.vals[0][1]["opoint"],
addFileNamesToCommandLine = True,
noEmptyFiles = True)
g.addEdge( luminosityListNode, luminositySumNode, sumLuminosity )
#
# Now, compute the lumi per-trigger. If there's just one dataset that contains
# the trigger we're looking for, we can just pass through from the previous
# luminosity node. Otherwise, we need to add them together (usually don't add
# them)
#
#luminositySumByTrigger = {}
luminositySumMiter = Miter()
#for trigger in luminosityByTrigger:
for trigger in luminosityMiter.iterGrouped( 'trigger' ):
# should make a wrapper for this magic
triggerName = trigger.vals[0][1][ 'trigger' ]
lumiSum = Node( "lumisum-%s" % triggerName )
g.addNode( lumiSum )
luminositySumMiter.add( lumiSum, trigger = triggerName ) #ByTrigger[ trigger ] = lumiSum
if len( trigger ) == 1:
g.addEdge( trigger.getOneValue(), lumiSum, NullEdge() )
elif len( trigger ) > 1:
lumiCollect = Node( "lumicollect-%s" % triggerName )
g.addNode( lumiCollect )
for dataset in trigger:
g.addEdge( dataset.getOneValue(), lumiCollect, NullEdge() )
noEmptyFiles=True,
)
g.addEdge(luminosityListNode, luminositySumNode, sumLuminosity)
#
# Now, compute the lumi per-trigger. If there's just one dataset that contains
# the trigger we're looking for, we can just pass through from the previous
# luminosity node. Otherwise, we need to add them together (usually don't add
# them)
#
# luminositySumByTrigger = {}
luminositySumMiter = Miter()
# for trigger in luminosityByTrigger:
for trigger in luminosityMiter.iterGrouped("trigger"):
# should make a wrapper for this magic
triggerName = trigger.vals[0][1]["trigger"]
lumiSum = Node("lumisum-%s" % triggerName)
g.addNode(lumiSum)
luminositySumMiter.add(lumiSum, trigger=triggerName) # ByTrigger[ trigger ] = lumiSum
if len(trigger) == 1:
g.addEdge(trigger.getOneValue(), lumiSum, NullEdge())
elif len(trigger) > 1:
lumiCollect = Node("lumicollect-%s" % triggerName)
g.addNode(lumiCollect)
for dataset in trigger:
g.addEdge(dataset.getOneValue(), lumiCollect, NullEdge())
dataTriggerMiter.add(trigger_map["hltdijet20"], bin="100..110")
dataTriggerMiter.add(trigger_map["hltdijet40"], bin="100..110")
dataTriggerMiter.add(trigger_map["hltdijet70"], bin="100..110")
dataTriggerMiter.add(trigger_map["hltdijet20"], bin="110..120")
dataTriggerMiter.add(trigger_map["hltdijet40"], bin="110..120")
dataTriggerMiter.add(trigger_map["hltdijet70"], bin="110..120")
dataTriggerMiter.add(trigger_map["hltdijet110"], bin="120..")
skiplessDataMiter = Miter()
for opoint in operating_points:
for bin in jet_bins:
for tmap in dataTriggerMiter.iterGrouped("bin"):
if tmap.vals[0][1]["bin"] != bin[1]:
continue
for onetmap in tmap.iterMany():
sample = onetmap[0][5]
trigger = onetmap[0]
step_postfix = "-%s-%s-%s-%s" % (sample, trigger[1], bin[0], opoint)
monitor_node = run_monitor_input_helper(
g,
jet_pt=bin[1],
tag=opoint,
fileKey=sample,
# trigger_name = trigger[4],
step_postfix=step_postfix,
muon_pt="6..",
for bin in jet_bins:
for opoint in operating_points:
step_postfix = "-%s-%s-%s-djet10mc" % (dataset[1]['dataset'],bin[0],opoint)
currNode = run_monitor_input_helper( g,
jet_pt = bin[1],
tag = opoint,
fileKey = dataset[1]['dataset'],
trigger_name = "HLT_BTagMu_DiJet10U",
step_postfix = step_postfix,
muon_pt = "6..",
input_files = dataset[0] )
dijet10QCDMiter.add( currNode, bin = bin[0], dataset=dataset[1]['dataset'], opoint=opoint )
singleMonitorMiter = Miter()
rootQCDMiter = Miter()
# merge the dijet10 QCD into one guy
for currmerge in dijet10QCDMiter.iterGrouped( 'bin', 'opoint' ):
step_postfix = "-%s-%s" % ( currmerge.vals[0][1]['bin'], currmerge.vals[0][1]['opoint'] )
mergeNode = merge_with_root_qcd_helper( g, "root-qcd" + step_postfix, step_postfix,
inputNodes = currmerge.getValues(),
triggerName = "hltjet10u",
lumiMiter = exports['luminosityMiter']
)
rootQCDMiter.add( mergeNode,
bin = currmerge.vals[0][1]['bin'],
opoint = currmerge.vals[0][1]['opoint'] )
singleMonitorMiter.add( mergeNode,\
trigger = 'hltdijet10u',
bin = currmerge.vals[0][1]['bin'],
opoint = currmerge.vals[0][1]['opoint'],
njetpt = True,
njeteta = True,
