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
0
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(alpha, dump)
    # cutTree.TAttach(alpha,skim)
    out.append(alpha)
    out.append(dump)
    return out
    pass
Пример #3
0
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
            # 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
            # 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
0
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()
Пример #5
0
def createRatioPlots(ratio,
                     truth=False,
                     pt_default=50.,
                     ht_default=350.,
                     pt1_default=125.,
                     pt2_default=100.):

    # Define pT bins to be used
    pt_min = 20
    pt_max = 50
    pt_step = 10

    # Define x3 from pT and HT defaults for analysis
    x3_default = (2. * pt_default) / (ht_default + pt_default)
    factor = (2. - x3_default) / x3_default

    # Determine maximum values for HT and Meff (from max pT value)
    ht_max = pt_max * 1. * factor
    meff_max = pt_max * 1. + ht_max

    # Iterate through AlphaT values
    for ii in range(0, 11):
        alphat = 500 + ii * 5

        # Iterate through pT values
        for pt in range(pt_min, pt_max + pt_step, pt_step):

            # Define HT and Meff from pT value
            ht = pt * 1. * factor
            meff = pt * 1. + ht

            # Define scaled pT values appropriate for given Meff
            pt1 = pt1_default * meff / meff_max
            pt2 = pt2_default * meff / meff_max

            # Iterate through scaled pT values
            for jj in range(0.8, 1.2, 0.1):

                # Define scaled pT values appropriate for given Meff
                pt1 *= jj
                pt2 *= jj

                # Define histograms
                dir = "Ratio" + str(alphat) + "Pt" + str(pt) + "Scale" + str(
                    jj)
                ratio.append(
                    RobPlottingOps(
                        PSet(
                            DirName=dir,
                            MinObjects=2,
                            MaxObjects=8,
                            Ratio=True,
                            MinJetPt=pt * 1.,
                            MinJetPt1=pt1,
                            MinJetPt2=pt2,
                            UseGen=truth,
                            MinGenPt=pt_min * 1.,
                            MaxGenMET=10.,
                            AlphaTcut=(alphat / 1000.),
                        ).ps()))
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
Пример #7
0
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
Пример #8
0
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
Пример #9
0
def metScaleSystematic(a, shiftdir, lepton):
    """ Apply the pfmet scale systematic as a MET filter.
    a : analysis object
    shiftdir : True (False) for upwards (downwards) shift
    lepton : 'Muon' or 'Electron'
    """
    pfmet_jecunc = PSet(
        #the location of the JECUncertainty txt file
        jecuncfile = "GR_R_42_V19_AK5PF_Uncertainty.txt",
        jecuncfileresidual = "GR_R_42_V19_AK5PF_L2L3Residual.txt",
        #either Muon or Electron
        lepton = lepton,
        #the jet threshold to consider for the JEC Uncertainties
        pfjetthresh = 10.0,
        #the percentage/100 to shift the unclustered energy by
        unclusteredshift = 0.1,
        #whether or not to apply the shift upwards (True) or downwards (False)
        shiftup = shiftdir
        )
    pfmetjecunc = wpol.pfMETJECUnc(pfmet_jecunc.ps())
    a.AddMETFilter("pfMET", pfmetjecunc)
    return pfmetjecunc
def CorrectionPset(susydir, corrfile):
    CorFile = open(susydir + corrfile, "r")
    Eta1 = []
    Eta2 = []
    Corr = []
    for aRow in CorFile:
        rows = aRow.split()
        for coloum in rows:
            Eta1.append(float(rows[0]))
            Eta2.append(float(rows[1]))
            Corr.append(float(rows[5]))
    ps = PSet(EtaLow=Eta1, EtaHi=Eta2, Correction=Corr)
    return ps
Пример #11
0
from data.Multijet_json291010 import *
from data.MultiJet_Run2010B_PromptReco_v2_408 import *
from data.JetMET_Run2010A_PromptReco_v4_11092010 import *
from data.JetMET_Run2010A_PromptReco_v4_24092010 import *
from data.edward_data_skim_allJetTo_Json151010 import *
from data.Jet_MultiJet_JetMet_38_JSonUpTo151010 import *
from data.MultiJet_json221010 import *
#from icf.utils import local_ntuple

import icf.utils as utils

id_working_point = eid.eff_80
pf_ele_id = PSet(
    IsoPF = 0.2,
    EleMva = 0.4,
    LostHits = 1.,
    ElecDxy = 0.04,
    ElecDz = 1.
    )

pf_mu_id = PSet(
    IsoPF = 0.2,
    ValidHits = 11.,
    MuonID = True,
    MuonDxy = 0.02,
    MuonDz = 1.
    )
sig_id = PSet(
    Cuts = id_working_point,
    Isolation = True,
    HoverE = True,
Пример #12
0
# 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() )
Пример #13
0
# conf_ak5_jMCpt.Common.print_out()
anal_ak5_jptMC = Analysis("AK5JPT")
addCutFlowMC(anal_ak5_jptMC)

# AK7 Calo

conf_ak7_caloMC = deepcopy(defaultConfig)
conf_ak7_caloMC.Ntuple = deepcopy(ak7_calo)
conf_ak7_caloMC.XCleaning = deepcopy(default_cc)
conf_ak7_caloMC.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak7_caloMC = Analysis("AK7Calo")
addCutFlowMC(anal_ak7_caloMC)
outDir = "../results_" + strftime("%d_%b_%H") + "//NoSmear37/"
ensure_dir(outDir)

testFile = PSet(
    Name="SingleEvent",
    Format=("ICF", 3),
    File="/home/hep/elaird1/84_darrens_event/event.root",
    Weight=1.0,
)

anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, MC)
# anal_ak5_pfMC.Run(outDir,conf_ak5_pfMC,MC)
# anal_ak5_pfMC.Run("../results_"+strftime("%d_%b_%H")+"//NoSmear",conf_ak5_pfMC,[QCD_AllPtBins_7TeV_Pythia])
# anal_ak5_jptMC.Run("../results_"+strftime("%d_%b_%H")+"//NoSmear",conf_ak5_jptMC,MC)
# anal_ak5_jptMC.Run("../results_"+strftime("%d_%b_%H")+"//NoSmear",conf_ak5_jptMC,[QCD_AllPtBins_7TeV_Pythia])
# anal_ak7_caloMC.Run("../results_"+strftime("%d_%b_%H")+"//NoSmear",conf_ak7_caloMC,MC)
# anal_ak7_caloMC.Run("../results_"+strftime("%d_%b_%H")+"//NoSmear",conf_ak7_caloMC,[QCD_AllPtBins_7TeV_Pythia])
Пример #14
0
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(
Пример #15
0
conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

from data.Run2011.HT_Run2011_promptReco_DCS import *
from data.Run2011.RA1ToBurn import *
from data.Run2011.HT42_incomplete import *
from data.Run2011.HT_Run2011A import *
from data.Run2011.HTRun2011AB import *
test = PSet(Name="test",
            Format=("ICF", 3),
            Weight=1.0,
            File="../../Ntuples/375_newFormat.root")
outdir = "../results_" + strftime("%d_%b_%H") + "/Data/"
ensure_dir(outdir)
anal_ak5_caloData.Run(outdir, conf_ak5_caloData, [HTRun2011AB])  #HT_Run2011A])
Пример #16
0
from libHadronic import *
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from copy import deepcopy
# from montecarlo.LMx import *
# from bryn_samples import *
# from montecarlo.WJets_Madgraph import *

# -----------------------------------------------------------------------------
# Samples
#from local_samples import *
# from met_samples import *
# from JetMetTau139375_Uncleaned import *

f = PSet(Name="Testing",
         File=["../SusyCAF_Tree_9_1.root"],
         Format=("ICF", 2),
         Weight=1.0)

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
Пример #17
0
conf.Common.Jets.EtCut = 20.0
conf.Common.Jets.EtaCut = 5.0
conf.Common.Muons.EtaCut = 2.4
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,
    EndCapEtaMax=2.5,
    MinElecRescaledPt=0.,
    c_ErNu_pt=0.,
    wID=+24.,  ##24 W+, -24 W-
    mW=80.398,
    mZ=90.188,
Пример #18
0
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from copy import deepcopy
from icf.JetCorrections import *

# -----------------------------------------------------------------------------
# Samples
#import yours in your running script

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #19
0
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
Пример #20
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[LM4])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[LM5])

#QCD
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[QCD_AllPtBins_7TeV_Pythia])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,QCD_Pythia8_ALL)
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,QCD_MadGraph_ALL)
#
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,QCD_Pythia6_384patch3_V14_00_02_ALL)
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,QCD_Pythia8_384patch3_V14_00_02_ALL)

#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[PhotonJet_AllPtBins_7TeV_Pythia])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,PhotonJet_MadGraph_ALL)

#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[ttbarTauola])
#wjets_madgraph_vols.FirstEntry = 0
#wjets_madgraph_vols.LastEntry = 199999
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[wjets_madgraph])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[wjets_madgraph_vols])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[zjets_madgraph])
#anal_ak5_caloMC.Run("../results_100_100_50_eff",conf_ak5_calo,[Zinvisible_jets])

myMagical14events = PSet(  #from 15pb-1
    Name="14magicalEvents",
    Weight=1.,
    Format=("ICF", 2),
    File=[
        "/vols/cms02/whyntie/public/ntuple_14eventsPassingAlphaT_15invpb.root",
    ])
#anal_ak5_caloData.Run("../results_100_100_50_eff",conf_ak5_calo,[myMagical14events]) #15.0865/pb
Пример #22
0
#!/usr/bin/env python
"""
Plot settings for Bryn's plotting operations.
"""

from libbryn import *

from icf.core import PSet

from copy import deepcopy

#plot configuration
pset_2om = PSet(
    DirName="250_infGev_2om",
    MinObjects=2,
    MaxObjects=15,
    StandardPlots=True,
)

pset_2eq = PSet(
    DirName="250_infGev_2eq",
    MinObjects=2,
    MaxObjects=2,
    StandardPlots=True,
)

pset_3om = PSet(
    DirName="250_infGev_3om",
    MinObjects=3,
    MaxObjects=15,
    StandardPlots=True,
Пример #23
0
# Notes:
# * A few variables need the 38 ntuples (pixel hits, number of chambers).

from icf.core import PSet

vbtfmuonidps = PSet(
    Verbose = False,
    LooseCuts = PSet(
        #ID variables
        RequireTrackerMuon = False,
        RequireGlobalMuon  = True,
        ChiSquaredCut = 999999.,
        NumChambersCut = 0,
        NumGlobalHits = 0,
        NumTrackerHits = 0,
        NumPixelHits = 0,
        DxyCut = 9999999.,
        ),#end of LooseCuts
    TightCuts = PSet(
        #ID variables
        RequireTrackerMuon = True,
        RequireGlobalMuon  = True,
        ChiSquaredCut = 10.0,
        NumChambersCut = 2,
        NumGlobalHits = 1,
        NumTrackerHits = 10,
        NumPixelHits = 1,
        DxyCut = 0.2, #it's in cm, right? This should be 2mm.
        ),#end of TightCuts
    )#end of VBTF electron ID PSet
Пример #24
0
        "HLT_HT550_v6",
        "HLT_HT550_v7",
        "",
        ]
    
NoQcdOption = False
if ( dataset == 2 or dataset == 3 ) :
    NoQcdOption = True
    
# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons=PSet(
Prefix="electron",
Suffix="Pat",
LooseID="EIDLoose",
TightID="EIDTight",
)

default_ntuple.Muons=PSet(
Prefix="muon",
Suffix="Pat",
LooseID="IsGlobalMuon",
TightID="IDGlobalMuonPromptTight",
)

default_ntuple.Taus=PSet(
Prefix="tau",
Suffix="Pat",
LooseID="TauIdbyTaNCfrOnePercent",
TightID="TauIdbyTaNCfrTenthPercent"
Пример #25
0
from samples.LMx import *
from samples.tt import *
from samples.zjets_madgraph_filepset import *
#import samples.JetMetTau as data
from samples.wjets import *
#from icf.utils import local_ntuple

import icf.utils as utils

id_working_point = eid.eff_80

sig_id = PSet(
    Cuts = id_working_point,
    Isolation = True,
    HoverE = True,
    DeltaEtaAtVtx = True,
    DeltaPhiAtVtx = True,
    SigmaIEtaIEta = True,
    Conversions = True,
    Ctrl = False
    )

bg_id = PSet(
    Cuts = id_working_point,
    Isolation = True,
    HoverE = True,
    DeltaEtaAtVtx = False,
    DeltaPhiAtVtx = False,
    SigmaIEtaIEta = True,
    Conversions = True,
    Ctrl = True
)
Пример #26
0
class eWPol:
    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

        self.met_cuts = {}
        for mht in mht_cuts:
            for met in met_cuts:
                self.met_cuts[(mht, met)] = fwk.OP_PFMETCut(met)
                self.pol_plots["mht%d_met%d" % (mht, met)] = wpol.RECO_ElPolPlots(
                    "RECO_ElPolPlots_PostMHT%dMET%d" % (mht, met),
                    self.pol_plots_cfg.ps())

        # MT Cuts and plots
        self.mt_cuts = {}
        for mht in mht_cuts:
            for mt in mt_cuts:
                self.mt_cuts[(mht, mt)] = wpol.ePFMTCut(lepton_type, mt)
                self.pol_plots["mht%d_mt%d" % (mht, mt)] = wpol.RECO_ElPolPlots(
                    "RECO_ElPolPlots_PostMHT%dMT%d" % (mht, mt),
                    self.pol_plots_cfg.ps())
		self.datamc_plots["mht%d_mt%d" % (mht, mt)] = wpol.DataMCPlots(
		    "RECO_DataMCPlots_PostMHT%dMT%d" % (mht, mt))
        self.cuts = {}
        self.tree = fwk.Tree(self.name)
        if event_no_dump:
            self.event_no_dump = wpol.EventNoDump("EventNoDump_MT50", "event_no")

    def buildTree(self,
                  basic_plots=True,
                  pol_plots=True,
                  vetos = [],
                  tree = None,
                  root_op = None):
        if tree is None:
            tree = self.tree

        veto_list = vetos
        if not basic_plots:
            veto_list.append("basic_plots")
        if not pol_plots:
            veto_list.append("pol_plots")
        ops_ = []
        for name, op in self.ops:
            veto = False
            print name
            for v in vetos:
                if name.startswith(v):
                    print "Removing operation '%s' from tree" % name
                    veto = True
            if not veto:
                ops_.append((name, op))
        if root_op is None:
            tree.Attach(ops_[0][1])
        else:
            tree.TAttach(root_op, ops_[0][1])

        idx = 1
        for name, op in ops_[1:]:
            tree.TAttach(ops_[idx-1][1], op)
            idx += 1
        for mht, mht_op in self.mht_cuts.iteritems():
             tree.TAttach(ops_[idx-1][1], mht_op)
             if pol_plots:
                 tree.TAttach(mht_op, self.pol_plots["mht%d" % mht])
                 tree.TAttach(mht_op, self.datamc_plots["mht%d" % mht])
        for (mht, met), met_op in self.met_cuts.iteritems():
             tree.TAttach(self.mht_cuts[mht], met_op)
             if pol_plots:
                 tree.TAttach(met_op, self.pol_plots["mht%d_met%d" % (mht, met)])
        for(mht, mt), mt_op in self.mt_cuts.iteritems():
             tree.TAttach(self.mht_cuts[mht], mt_op)
             if pol_plots:
                 tree.TAttach(mt_op, self.pol_plots["mht%d_mt%d" % (mht, mt)])
                 tree.TAttach(mt_op, self.datamc_plots["mht%d_mt%d" % (mht, mt)])
	     # Fringe hack
             if mht == 50 and self.event_no_dump is not None:
	         tree.TAttach(mt_op, self.event_no_dump)

    def installFilters(self, a):
        a.AddElectronFilter("PostCC", self.loose_ele_id)
        a.AddJetFilter("PostCC", self.jet_dr)

    def GetOp(self, name):
        for n, o in self.ops:
            if n == name:
                return o
        return None
Пример #27
0
import setupSUSY

from libFrameworkSUSY import *
from libHadronic import *
from liballhadronic import *
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from copy import deepcopy

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)

default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)

default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
Пример #28
0
# -----------------------------------------------------------------------------
# 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,
Пример #29
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #30
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #31
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #32
0
    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_2600to3000_7TeV_pythia8_Summer10_START36_V10_S09_v2,
    QCD_Pt_3000to3500_7TeV_pythia8_Summer10_START36_V10_S09_v2
]

QCD_PYTHIA6 = [
    QCD_Pt80_7TeV_Pythia, QCD_Pt170_7TeV_Pythia, QCD_Pt300_7TeV_Pythia,
    QCD_Pt470_7TeV_Pythia, QCD_Pt800_7TeV_Pythia, QCD_Pt1400_7TeV_Pythia
]

pset_standard_hadronic_smallPFMET_noLowerLimit = PSet(
    DirName="HadronicCommon_smallPFMET_noLowerLimit",
    MinObjects=1,
    MaxObjects=10,
    StandardPlots=True,
    CountPlots=False,
    DeadECALPlots=True,
    BabyPlots=True)

pset_standard_hadronic_beforeAlphaT = deepcopy(
    pset_standard_hadronic_smallPFMET_noLowerLimit)
pset_standard_hadronic_beforeAlphaT.DirName = "HadronicCommon_beforeAlphaT"

pset_standard_hadronic_PFMET = deepcopy(
    pset_standard_hadronic_smallPFMET_noLowerLimit)
pset_standard_hadronic_PFMET.DirName = "HadronicCommon_PFMET"

pset_standard_hadronic_PFMET_babyMHT = deepcopy(
    pset_standard_hadronic_smallPFMET_noLowerLimit)
pset_standard_hadronic_PFMET_babyMHT.DirName = "HadronicCommon_PFMET_babyMHT"
Пример #33
0
)
ElHad_Prompt_deadecal_ps._quiet_set("Accept", False)
ElHad_Prompt_Cut_BE = EventSelector(ElHad_May10ReReco_deadecal_ps.ps())

ElHad_Prompt_TP_deadecal_ps = getDeadECALList(
    "../python/deadecal/ElectronHad_Run2011A_PromptReco_v4_missed1Job_filteredEventsTP_cfi.py"
)
ElHad_Prompt_TP_deadecal_ps._quiet_set("Accept", False)
ElHad_Prompt_Cut_TP = EventSelector(ElHad_May10ReReco_TP_deadecal_ps.ps())


# Trigger
trg_set1 = PSet(Triggers=["HLT_Mu5_HT200_v*"], Verbose=False, UsePreScaledTriggers=False)
trg_set2 = PSet(
    #   Triggers = ["HLT_Mu8_HT200_v*","HLT_Mu15_HT200_v*","HLT_Mu20_HT200_v*","HLT_Mu30_HT200_v*","HLT_HT200_Mu5_PFMHT35_v*""HLT_HT*_Mu*_PFMHT*_v*"],
    Triggers=["HLT_Mu*_HT*_v*", "HLT_HT*_Mu*_PFMHT*_v*"],
    Verbose=False,
    UsePreScaledTriggers=False,
)
trg_set4 = PSet(
    #    Triggers = ["HLT_Mu15_v*","HLT_Mu20_v*","HLT_IsoMu17_v*"],
    Triggers=["HLT_Mu15_HT200_v*", "HLT_Mu8_HT200_v*"],
    Verbose=False,
    UsePreScaledTriggers=False,
)

# Electron triggers
etrg_set1 = PSet(
    Triggers=[
        "HLT_Ele10_CaloIdL_CaloIsoVL_TrkIdVL_TrkIsoVL_HT200_v*",
        "HLT_Ele10_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_HT200_v*",
    ],
Пример #34
0
QCD_Pt_300to380_7TeV_pythia8_Summer10_START36_V10_S09_v1,
QCD_Pt_30to50_7TeV_pythia8_Summer10_START36_V10_S09_v2,
QCD_Pt_380to470_7TeV_pythia8_Summer10_START36_V10_S09_v1,
QCD_Pt_470to600_7TeV_pythia8_Summer10_START36_V10_S09_v1,
QCD_Pt_50to80_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_80to120_7TeV_pythia8_Summer10_START36_V10_S09_v1]

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons=PSet(
Prefix="electron",
Suffix="Pat",
LooseID="EIDLoose",
TightID="EIDTight",
)
default_ntuple.Muons=PSet(
Prefix="muon",
Suffix="Pat",
LooseID="IsGlobalMuon",
TightID="IDGlobalMuonPromptTight",
)
default_ntuple.Taus=PSet(
Prefix="tau",
Suffix="Pat",
LooseID="TauIdbyTaNCfrOnePercent",
TightID="TauIdbyTaNCfrTenthPercent"
)
default_ntuple.Jets=PSet(
Пример #35
0
    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

        self.met_cuts = {}
        for mht in mht_cuts:
            for met in met_cuts:
                self.met_cuts[(mht, met)] = fwk.OP_PFMETCut(met)
                self.pol_plots["mht%d_met%d" % (mht, met)] = wpol.RECO_ElPolPlots(
                    "RECO_ElPolPlots_PostMHT%dMET%d" % (mht, met),
                    self.pol_plots_cfg.ps())

        # MT Cuts and plots
        self.mt_cuts = {}
        for mht in mht_cuts:
            for mt in mt_cuts:
                self.mt_cuts[(mht, mt)] = wpol.ePFMTCut(lepton_type, mt)
                self.pol_plots["mht%d_mt%d" % (mht, mt)] = wpol.RECO_ElPolPlots(
                    "RECO_ElPolPlots_PostMHT%dMT%d" % (mht, mt),
                    self.pol_plots_cfg.ps())
		self.datamc_plots["mht%d_mt%d" % (mht, mt)] = wpol.DataMCPlots(
		    "RECO_DataMCPlots_PostMHT%dMT%d" % (mht, mt))
        self.cuts = {}
        self.tree = fwk.Tree(self.name)
        if event_no_dump:
            self.event_no_dump = wpol.EventNoDump("EventNoDump_MT50", "event_no")
Пример #36
0
from data.JetMET_Run2010A_PromptReco_v4_24092010 import *
from data.edward_data_skim_allJetTo_Json151010 import *
from data.Jet_MultiJet_JetMet_38_JSonUpTo151010 import *
from data.MultiJet_json221010 import *
#from icf.utils import local_ntuple

import icf.utils as utils

id_working_point = eid.eff_80

sig_id = PSet(
    Cuts = id_working_point,
    Isolation = True,
    HoverE = True,
    DeltaEtaAtVtx = True,
    DeltaPhiAtVtx = True,
    SigmaIEtaIEta = True,
    Conversions = True,
    CorrEEMisalig = True,
    Ctrl = False,
    Mu = True
    )

bg_id = PSet(
    Cuts = id_working_point,
    Isolation = True,
    HoverE = True,
    DeltaEtaAtVtx = False,
    DeltaPhiAtVtx = False,
    SigmaIEtaIEta = True,
    Conversions = True,
    CorrEEMisalig = True,
Пример #37
0
from libHadronic import *
from icf.core import PSet,Analysis
from icf.config import defaultConfig
from copy import deepcopy

# -----------------------------------------------------------------------------
# Samples
  # put yout samples in here or import a file containing your data samples
from met_samples import *
# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons=PSet(
Prefix="electron",
Suffix="Pat",
LooseID="EIDLoose",
TightID="EIDTight",
)
default_ntuple.Muons=PSet(
Prefix="muon",
Suffix="Pat",
LooseID="IsGlobalMuon",
TightID="IDGlobalMuonPromptTight",
)
default_ntuple.Taus=PSet(
Prefix="tau",
Suffix="Pat",
)
default_ntuple.Jets=PSet(
Prefix="ic5Jet",
Suffix="Pat",
    ps = PSet(EtaLow=Eta1, EtaHi=Eta2, Correction=Corr)
    return ps


YourSusyDir = commands.getoutput('echo $SUSY_WORKING_SW_DIR')
caloResCorrections = CorrectionPset(
    YourSusyDir, "/framework/python/icf/ResidualJetEnergyCorrections.txt")
ak5pfResCorrections = CorrectionPset(
    YourSusyDir, "/framework/python/icf/Spring10DataV2_L2L3Residual_AK5PF.txt")

evDiagPs = PSet(
    LatexFilename="diagSlides.tex",
    Triggers=[
        "HLT_HT100U",
        "HLT_HT120U",
        "HLT_HT140U",
        "HLT_HT150U_v3",
        "HLT_HT160U_v1",
        "HLT_HT160U_v3",
    ],
    #    CaloJetCorrections = caloResCorrections.ps(),
)
evDiag = OP_PFDiagnosticSlideMaker(evDiagPs.ps(), caloResCorrections.ps(),
                                   ak5pfResCorrections.ps())
#cutTreeData.TAttach(MHToverMET2om_350plus,evDisplay)
cutTreeData.TAttach(MHToverMET2om_350plus, evId)
cutTreeData.TAttach(MHToverMET2om_350plus, evDiag)
#
#
# number of common jets = 2
cutTreeData.TAttach(count_total, numComJets2Please)
cutTreeData.TAttach(numComJets2Please, count_numComJets2eq)
import setupSUSY
from libFrameworkSUSY import *
from libHadronic import *
from icf.core import PSet,Analysis
from time import strftime
from time import strftime
from Binned_DiMuon_Config import *

mu_id = PSet(
    doJetLoop = False,
    MuID = "Tight",
    MinPt = 10.,
    MaxEta = 2.5,
    MaxIsolation = 0.1,
    DRMuJet = 0.3,
    MaxGlbTrkDxy = 0.02,
    MinGlbTrkNumOfValidHits = 11,
    SegMatch2GlbMu = 1,
    PixelHitsOnInrTrk = 1,
    MaxInrTrkDz = 1.
        )

"""mu_id = PSet(
    doJetLoop = False,
    MuID = "Tight",
    MinPt = 10.,
    MaxEta = 2.5,
    MaxIsolation = 5.0,
    DRMuJet = 0.3,
    MaxGlbTrkDxy = 0.02,
Пример #40
0
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,
Пример #41
0
#conf.Ntuple.Taus.Prefix="tauHPS"
conf.Ntuple.Taus.Suffix="PF"
conf.Common.Taus.TightID = False
conf.Common.Taus.ApplyID = False
conf.Common.Taus.PtCut =  15.0
conf.Common.Taus.EtaCut = 2.4



IsoTightTau=2
IsoLooseTau=0
UseHPS=True
tau_id = PSet(
    TauSource=conf.Ntuple.Taus.Prefix,
    tauIDbyIso = IsoTightTau,
    AgainstElectron = True,
    AgainstMuon = True,
    UseHps = UseHPS
)


ele_id = PSet(
    OldNtu = False,
    Cuts = lepconf.ele_ptThres_vbtf80,
    IsolationSR = False,      # high pt
    IsoBarEC = True,
    HoverE = True,
    DeltaAtVtx = True,
    SigmaIEtaIEta = True,
    Conversions = True,
    ConversionsExtra = True,    # not in ntuple V00-08
Пример #42
0
from validation.cutCounters_cff import * #TW for convenience...
from copy import deepcopy
from SUSYSignalScan.plottingStuff_cfg import tripleScale,smsOps,addBinnedStuff,isCmssm,isSms
from hadronic.samples_cff import *
from montecarlo.PhysicsProcesses_mSUGRA_tanbeta10Fall10v1 import *
from montecarlo.PhysicsProcesses_mSUGRA_tanbeta3Fall10v1 import *
from SUSYSignalScan.utils import *
import os

mu_id = PSet(
   MuID = "Tight",
   MinPt = 10.,
   MaxEta = 2.5,
   MaxIsolation = 0.1,
   DRMuJet = 0.3,
   MaxGlbTrkDxy = 0.02,
   MinGlbTrkNumOfValidHits = 11,
   SegMatch2GlbMu = 1,
   PixelHitsOnInrTrk = 1,
   MaxInrTrkDz = 1.,
   doJetLoop = False,
   )


def mkdir(path) :
  try: os.makedirs(path)
  except OSError as e :
    if e.errno!=17 : raise e

def checkSwitches(d) :
  assert d["model"] in ["T1","T2","T2bb","T2tt","tanB3","tanB10","tanB40","tanB50"]
Пример #43
0
#!/usr/bin/env python

# standard imports
import setupSUSY
from libFrameworkSUSY import *
import libOneLepton as one
import libWPol as wpol
from copy import deepcopy
from icf.core import PSet, Analysis

# The leptonic specicif settings&files
from onelepton_settings import *
from onelepton_MC_samples import *
from onelepton.filters import metScaleSystematic, vertexReweighting, metResolutionSystematic

syst_ps = PSet(Polarisation=True, Lepton="Muon", METResolution=True)
syst_plots = one.SystematicsPlots("Systematics", syst_ps.ps())

syst = Analysis("Systematics")
tree = Tree("Main")
syst += tree

tree.Attach(syst_plots)

syst.Run("./systematics", conf, samplesMC)
Пример #44
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #45
0
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

tedSkim = PSet(Name="TedSkim",
               Format=("ICF", 3),
               Weight=1.0,
               File="../resultsWithSingleTop//Data/AK5Calo_Jets.root")
#"~elaird1/public_html/73_candidates/v3/350_bin/calo.root")

from data.Run2011.HT_Run2011_promptReco_DCS import *
from data.Run2011.HT42_incomplete import *
from data.Run2011.RA1ToBurn import *
outDir = "../results/Data43"
ensure_dir(outDir)
anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HT42_incomplete])

# from data.MultiJet_Run2010B_Nov4ReReco_v1 import *

#
# anal_ak5_pfData.Run("../results/Data",conf_ak5_pfData,[HT_Run2011_promptReco_DCS])
# anal_ak5_jptData.Run("../results/",conf_ak5_jptData,data)
Пример #46
0
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)
Пример #47
0
    QCD_Pt_2200to2600_7TeV_pythia8_Summer10_START36_V10_S09_v2,
    QCD_Pt_3000to3500_7TeV_pythia8_Summer10_START36_V10_S09_v2
]

MCS = [QCD_Pt_230to300_7TeV_pythia8_Summer10_START36_V10_S09_v2]

ttWPlottingOpsMHT120 = RECO_ttWPlotting("ttW_AfterMHT120")
ttWPlottingOps100 = RECO_ttWPlotting("ttW_After250")
ttWPlottingOps350 = RECO_ttWPlotting("ttW_After350")
ttWPlottingOpsAT = RECO_ttWPlotting("ttW_AfterAT")

psetS = PSet(DirName="After_350Gev",
             MinObjects=2,
             MaxObjects=9,
             StandardPlots=True,
             DeadECALPlots=False,
             CleaningControlPlots=False,
             MECPlots=False,
             DeadECALFile="./deadRegionList_START36_V9.txt",
             DeadECAL_MinJetPtCut=10.,
             DeadECAL_MinBiasCut=0.5)

HadStandard350 = HadronicCommonPlots(psetS.ps())
psetMHT = PSet(
    #  DirName      = "HadronicCommon_1",
    #  MinObjects   = 2,
    #  MaxObjects   = 6,
    #  StandardPlots     = False,
    #  DeadECALPlots = True,
    #  CleaningControlPlots = True,
    #  MECPlots = True,
    #  DeadECALFile = "./deadRegionList_START36_V9.txt",
      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"
Пример #49
0
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

tedSkim = PSet(Name="TedSkim",
               Format=("ICF", 3),
               Weight=1.0,
               File="../results_" + strftime("%d_%b_%H") +
               "/WithSingleTop//Data/AK5Calo_Jets.root")
#"~elaird1/public_html/73_candidates/v3/350_bin/calo.root")

from data.Run2011.HT_Run2011_promptReco_DCS import *
from data.Run2011.HT42_incomplete import *
from data.Run2011.RA1ToBurn import *
from data.Run2011.HTRun2011AB import *

outDir = "../results_" + strftime("%d_%b_%H") + "//Data37"
ensure_dir(outDir)
# HT42_incomplete.LastEntry = 100
# HT42_incomplete.File = HT42_incomplete.File[0:1]
anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HTRun2011AB])
Пример #50
0
#conf.Ntuple.Taus.Prefix="tauHPS"
#conf.Ntuple.Taus.Suffix="PF"
conf.Ntuple.Taus.Suffix="PF"
conf.Common.Taus.TightID = False
conf.Common.Taus.ApplyID = False
conf.Common.Taus.PtCut =  15.0
conf.Common.Taus.EtaCut = 2.3


IsoTightTau=1
IsoLooseTau=0

tau_id = PSet(
    TauSource=conf.Ntuple.Taus.Prefix,
    tauIDbyIso = IsoTightTau,
#     tauIDbyDecay = True,
     AgainstElectron = True,
     AgainstMuon = True
)


ele_id = PSet(
    Cuts = lepconf.ele_ptThres_vbtf80,
)

muon_id = PSet(#we're missing the ecalvetoDep and hcalvetoDep in the ntuples (for MIP requirement) - AGB 22/10/10
   Cuts = lepconf.muonCuts,
)
elecSelFilter=CustomEleId(ele_id.ps())
muonSelFilter=CustomMuonId(muon_id.ps())
tauSelFilter=CustomTauId(tau_id.ps())
Пример #51
0
# 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,
)
Пример #52
0
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from icf.utils import json_to_pset
from allhadronic.golden_cff import *
from copy import deepcopy

pwd = commands.getoutput('echo $SUSY_WORKING_SW_DIR')
output_dir = pwd + "/results/v60"

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)

default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)

default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
Пример #53
0
    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" ,
Пример #54
0
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from copy import deepcopy
from icf.JetCorrections import *

# -----------------------------------------------------------------------------
# Samples
#import yours in your running script

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
    Suffix="Pat",
Пример #55
0
#from icf.core import PSet

# This is the default configuration
from icf.config import defaultConfig


# import default configuration
conf=defaultConfig

ele_id = PSet(
    Cuts = lepconf.ele_ptThres_vbtf80,
    IsolationSR = True,      # high pt
    HoverE = True,
    DeltaAtVtx = True,
    SigmaIEtaIEta = True,
    Conversions = True,
    ConversionsExtra = False,    # not in ntuple V08
    IsolationPF = False,
    LowPtSel = True,
    LowPtID = True,
    LowPtIso = True    
)



muon_id = PSet(
   Cuts = lepconf.muonCuts,
   IsolationPF = False,
   IsolationSR = True,
   LowPtIso = True,
   MuonID = True,
Пример #56
0
import setupSUSY
from libFrameworkSUSY import *
from libHadronic import *
from liballhadronic import *
from icf.core import PSet, Analysis
from icf.config import defaultConfig
from copy import deepcopy

# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)

default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)

default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
Пример #57
0
# Common cut definitions
#Avaiable criteria for MC and for Data are at current slightly different Hence the making of two trees
#DataOnly!

from icf.JetCorrections import *
corPset =  CorrectionPset("ResidualJetEnergyCorrections.txt")
corPset =  CorrectionPset("Spring10DataV2_L2L3Residual_AK5PF.txt")
JetCorrections = JESCorrections( corPset.ps() ,True)
NoiseFilt= OP_HadronicHBHEnoiseFilter()
selection = OP_GoodEventSelection()



HT_Trigger_PS = PSet(
    Verbose = False,
    UsePreScaledTriggers = False,
    Triggers = [
        ]
    )


Cross_Trigger_PS = PSet(
    Verbose = False,
    UsePreScaledTriggers = False,
    Triggers =[]
    )


def PreScaledPair(cutTree = None, cut = None, NumeratorTrig = None, DenominatorTrig = None,Label = "", Debug = False):
  if str(type(NumeratorTrig)) != "<type 'list'>":
    NumeratorTrig = [NumeratorTrig]
  if str(type(DenominatorTrig)) != "<type 'list'>":
Пример #58
0
    try:
        os.makedirs(path)
    except OSError as exc:  # Python >2.5
        if exc.errno == errno.EEXIST:
            pass
        else:
            raise


# -----------------------------------------------------------------------------
# Reading the collections from the ntuple

default_ntuple = deepcopy(defaultConfig.Ntuple)
default_ntuple.Electrons = PSet(
    Prefix="electron",
    Suffix="Pat",
    LooseID="EIDLoose",
    TightID="EIDTight",
)
default_ntuple.Muons = PSet(
    Prefix="muon",
    Suffix="Pat",
    LooseID="IsGlobalMuon",
    TightID="IDGlobalMuonPromptTight",
)
default_ntuple.SecMuons = PSet(Prefix="muon", Suffix="PF")
default_ntuple.Taus = PSet(Prefix="tau",
                           Suffix="Pat",
                           LooseID="TauIdbyTaNCfrOnePercent",
                           TightID="TauIdbyTaNCfrTenthPercent")
default_ntuple.Jets = PSet(
    Prefix="ic5Jet",
Пример #59
0
# This is the default configuration
from icf.config import defaultConfig


# import default configuration
conf=defaultConfig

ele_id = PSet(
     ParticleFlow=True,
    OldNtu = False,
    Cuts = lepconf.ele_ptThres_vbtf80,
    IsolationSR = True,      # high pt
    HoverE = True,
    DeltaAtVtx = True,
    SigmaIEtaIEta = True,
    Conversions = True,
    ConversionsExtra = True,    # not in ntuple V08
    IsolationPF = False,
    LowPtSel = True,
    LowPtID = True,
    LowPtIso = True,
    ChargeComp = True,
    D0sel = True
)



muon_id = PSet(
    ParticleFlow=True,
   Cuts = lepconf.muonCuts,
   IsolationPF = False,
Пример #60
0
#config.XCleaning.ResolveConflicts=True
#config.XCleaning.Jets.PtCut=40.
#config.XCleaning.Jets.EtaCut=3.
#config.XCleaning.Photons.EtCut=40.
#config.XCleaning.Photons.IDReq=2
config.Common.ApplyXCleaning=False
config.Common.Jets.PtCut=40.
config.Common.Jets.EtaCut=2.5
config.Common.Photons.EtCut=80.
config.Common.Photons.IDReq=2

# -----------------------------------------------------------------------------
# User Variables

SelectedJetCuts = PSet(
  Et = 80.,
  DeltaR = 0.4 
)
SelectedPhotonCuts = PSet(
  Et = 80.,
  EtaEBMax = 1.4442,
  EtaEEMin = 1.56,
  EtaEEMax = 2.5,
  IDReq = 3
)
selected_jets = SelectedJets(SelectedJetCuts.ps())
selected_photons = SelectedPhotons(SelectedPhotonCuts.ps())

# -----------------------------------------------------------------------------
# Cut/Filter Operations

Trigger_PSet = PSet(