def tripleScale(model="tanB10", cutTree=None, cut=None, label="", SM=False):
from libSUSYSignalScan import OP_mSuGraPlottingOps, OP_SMPlottingOps
out = []
for scale in q2_scales: # see onelepton_settings.py
scale2 = scale.replace(".", "")
if scale2[-1] == "0":
scale2 = scale2[:-1]
pset = deepcopy(msugra_pset)
mSUGRAFile = "%s/SUSYSignalScan/textfiles/scale_xsection_nlo%s_m0_m12_%s_1v1.txt" % (
susyDir(),
scale,
model.replace("tanB", ""),
)
ps = getMSUGRAPSet(mSUGRAFile)
pset._quiet_set("NLO", ps)
pset.DirectoryName = "Counter_BSMGrid_%s_scale%s" % (label, scale2)
if SM:
op = OP_SMPlottingOps(pset.ps())
else:
op = OP_mSuGraPlottingOps(pset.ps())
out.append(op)
out.append(ps)
if cut is not None:
cutTree.TAttach(cut, op)
else:
cutTree.Attach(op)
return (op, out)
from copy import deepcopy
import libWPol as wpol
import libHadronic as had
from icf.core import PSet, Analysis, susyDir
from icf.utils import json_to_pset
#from onelepton_settings import *
from onelepton_dataRA4_sync_settings import *
from onelepton.filters import reweightVertices
from onelepton.utils import getDeadECALList
json = JSONFilter("Cert_160404-165542",
json_to_pset("%s/onelepton/json/Json_Latest.txt" % susyDir()))
# lepton specific cuts
TwoMu = OP_NumComMuons("==",2)
selection = OP_GoodEventSelection()
skim_ps=PSet(
SkimName = "myskim",
DropBranches = False,
Branches = [
" keep * "
]
)
skim = SkimOp(skim_ps.ps())
ZPlots = OP_Zplots("jdt")
# Create analysis cut tree
Triggers=[
"HLT_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_HT250_v*",
"HLT_Ele10_CaloIdL_CaloIsoVL_TrkIdVL_TrkIsoVL_HT200_v*",
"HLT_HT200_Ele5_CaloIdVL_TrkIdVL_CaloIsoVL_TrkIsoVL_PFMHT35_v*",
"HLT_HT250_Ele5_CaloIdVL_TrkIdVL_CaloIsoVL_TrkIsoVL_PFMHT35_v*",
"HLT_HT300_Ele5_CaloIdVL_TrkIdVL_CaloIsoVL_TrkIsoVL_PFMHT40_v*",
"HLT_HT350_Ele5_CaloIdVL_TrkIdVL_CaloIsoVL_TrkIsoVL_PFMHT45_v*",
"HLT_HT300_Ele5_CaloIdVL_CaloIsoVL_TrkIdVL_TrkIsoVL_PFMHT40_v*",
"HLT_HT350_Ele5_CaloIdVL_CaloIsoVL_TrkIdVL_TrkIsoVL_PFMHT45_v*",
],
Verbose=False,
UsePreScaledTriggers=False,
)
json = JSONFilter("Cert_160404-165542", json_to_pset("%s/onelepton/json/RA4muonSync.txt" % susyDir()))
triggerData4X = OP_MultiTrigger(trg_set2.ps())
triggerData4XCtrl = OP_MultiTrigger(trg_set4.ps())
triggerData42X_PromptReco = OP_MultiTrigger(trg_set5.ps())
triggerData42XCtrl_PromptReco = OP_MultiTrigger(trg_set6.ps())
triggerData42X_ReReco = OP_MultiTrigger(trg_set2.ps())
triggerData42XCtrl_ReReco = OP_MultiTrigger(trg_set4.ps())
etriggerData42X_ReReco = OP_MultiTrigger(etrg_set1.ps())
etriggerData42X_PromptReco = OP_MultiTrigger(etrg_set2.ps())
etriggerData42XCtrl_ReReco = OP_TriggerCut("HLT_Ele27_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v3")
etriggerData42XCtrl_PromptReco = OP_MultiTrigger(etrg_set3.ps())
def setupSUSYWeighting(a):
xs = susy_utils.readLOXS("%s/SUSYSignalScan/textfiles/goodModelNames_10_0_1.txt" % susyDir())
xsps = susy_utils.xsToPSet(xs)
from libSUSYSignalScan import SignalScanLOCrossSectionWeighting
weight = SignalScanLOCrossSectionWeighting(xsps.ps())
a.AddWeightFilter("Weight", weight)
return [weight]
UsePreScaledTriggers=False,
)
etrg_RA4sync_May10 = PSet(
Triggers=[
"HLT_Ele10_CaloIdL_CaloIsoVL_TrkIdVL_TrkIsoVL_HT200_v*",
"HLT_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_HT200_v*",
"HLT_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_HT250_v*",
"HLT_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_HT250_PFMHT25_v*",
],
Verbose=False,
UsePreScaledTriggers=False,
)
json = JSONFilter("Cert_160404-165542", json_to_pset("%s/onelepton/json/Json_Latest.txt" % susyDir()))
# json_to_pset("%s/onelepton/json/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON_v3.txt" % susyDir())) ## for EleHad RA4 Sync exercise
triggerData4X = OP_MultiTrigger(trg_set2.ps())
triggerData4XCtrl = OP_MultiTrigger(trg_set4.ps())
triggerData42X_PromptReco = OP_MultiTrigger(trg_set5.ps())
triggerData42XCtrl_PromptReco = OP_MultiTrigger(trg_set6.ps())
triggerData42X_ReReco = OP_MultiTrigger(trg_set2.ps())
triggerData42XCtrl_ReReco = OP_MultiTrigger(trg_set4.ps())
etriggerData42X_ReReco = OP_MultiTrigger(etrg_set1.ps())
etriggerData42X_PromptReco = OP_MultiTrigger(etrg_set2.ps())
