data_datasets = [["/BTau/Run2010B-Dec22ReReco_v1/AOD", "RUN2010B"], ["/BTau/Run2010A-Dec22ReReco_v1/AOD", "RUN2010A"]]


#
# Run s8_monitor_input over QCD without skipping events
#
dijet10QCDMiter = Miter()
for dataset in exports["qcd"].iterMany():
    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(),
Example #2
0
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..",
                            data = "1",
                            input_files = exports['data'].getValues( dataset = sample,trigger=trigger[1] ))
                singleMonitorMiter.add( monitor_node,dataset = sample,
                                                trigger = trigger[1],
                                                bin     = bin[0],
                                                opoint  = opoint,
                                                njetpt  = True,
                                                njeteta = True,
                                                privert = True,
                                                type    = "data" )
                skiplessDataMiter.add( monitor_node, dataset = sample,
                                                     trigger = trigger[1],
                                                     bin     = bin[0],
Example #3
0
    for bin in jet_bins:
        for sample in qcd_datasets:
            for trigger in trigger_list:
                step_postfix = "-%s-%s-%s-%s-noskip" % (sample[1],trigger[1],bin[0],opoint)
                if not opoint in skiplessS8Monitor:
                    skiplessS8Monitor[ opoint ] = {}
                if not bin[0] in skiplessS8Monitor[ opoint ]:
                    skiplessS8Monitor[ opoint ][ bin[0] ] = {}
                if not trigger[1] in skiplessS8Monitor[ opoint ][ bin[0] ]:
                    skiplessS8Monitor[ opoint ][ bin[0] ][ trigger[1] ] = []

                currNode = run_monitor_input_helper( g,
                            jet_pt = bin[1],
                            tag = opoint,
                            fileKey = sample[1],
                            trigger_name = trigger[4],
                            step_postfix = step_postfix,
                            muon_pt = "6..",
                            input_files = qcdTreeMiter.getOneValue( dataset =  sample[1] ),
                            simulate_trigger = triggers_to_simulate[ trigger[1] ])
                skiplessS8Monitor[ opoint ][ bin[0] ][ trigger[1] ].append( currNode )
                singleMonitorMiter.add( currNode,dataset = sample[1],
                                                trigger = trigger[1],
                                                bin     = bin[0],
                                                opoint  = opoint,
                                                njetpt  = True,
                                                njeteta = True,
                                                privert = True,
                                                type    = "skipless" )