# We can override the configuration we imported
# Here, we lower our Jet Et Cut to 5 GeV
conf.Common.Jets.EtCut = 5.0
# Create the analysis
a = Analysis("Demo")
# With the V2 code you can add linear sets of operations as before
# Here we will use the Tree class to have branching operations
# Create a Tree called Main
tree = Tree("Main")
# Initialise our operations
selection = OP_GoodEventSelection()
json = JSONFilter("Alex's JSON", json_to_pset("../Cert_160404-161216_7TeV_PromptReco_Collisions11_JSON.txt", True))
dijet = OP_NumComJets(">=", 2)
trijet = OP_NumComJets(">=", 3)
plot2 = OP_KinSuiteComPlot("plot2", 10, 2)
plot3 = OP_KinSuiteComPlot("plot3", 10, 2)
# Add the tree to our analysis
a += tree
tree.Attach(selection)
tree.TAttach(selection, json)
tree.TAttach(json, dijet)
tree.TAttach(dijet, plot2)
tree.Attach(trijet)
tree.TAttach(trijet, plot3)
# Define the list of samples to run over
a += tree
if iRecType==0 or iRecType==4 or iRecType==5 : #Z data or W data
trigger_bits = (
"HLT_Ele27_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v1",
"HLT_Ele32_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v1",
"HLT_Ele32_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v2",
"HLT_Ele32_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v3",
"HLT_Ele32_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_v4",
)
triggers = PSet(
Triggers=trigger_bits,
Verbose=False,
UsePreScaledTriggers=False
)
Trigg = OP_MultiTrigger(triggers.ps())
jsonfil=JSONFilter("Cert_160404-167151_7TeV",json_to_pset("/vols/cms03/mjarvis/SUSYv2/WCharge/scripts/json_folder/Cert_160404-166967_7TeV_May10ReRecoPlusPromptReco_840pb.json"))
json_output=JSONOutput("filtered")
tree.Attach(Trigg)
tree.TAttach(Trigg,jsonfil)
tree.TAttach(jsonfil,json_output)
tree.TAttach(json_output,ZeroMuons)
else :
tree.Attach(ZeroMuons)
tree.TAttach(ZeroMuons,Default_Asym) #
if iRecType in [0,4,5]:#ZDATA WDATA
samples =[
SingleElectron_Run2011A_May10ReReco_v1,
SingleElectron_Run2011A_PromptReco_v4_213,
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())
#=========== MC Sample Lists ===============
#MC_2012_Low = Top +WJets+ZJets+DiBoson+DY
MC_2012_Low = TTJets + ZJets
MC_2012_High =TTJets + ZJets
Photon_MC = Photon
#=============================
#===== Data Sample Lists ===========
Had_2012 = [ HT_Run2012A_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob229, HTMHT_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob228,HTMHT_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob238 ]
Muon_2012 = [ SingleMu_Run2012A_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob228, SingleMu_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob229, SingleMu_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob238 ]
Photon_2012 = [ Photon_Run2012A_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob228, SinglePhoton_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob229,SinglePhoton_Run2012B_PromptReco_v1_V17_0_taus_0_doTypeIMetReco_1_zmengJob238 ]
#===================================
def checkSwitches(d) :
assert d["json"] in ["ICHEP_18June.json","Cert_190456-200601_8TeV_PromptReco_Collisions12_JSON_v2.txt"]
def switches():
d = {}
d["reweight_samples"] = [(PU_2012,MC_2012_Low),(PU_2012,MC_2012_High),(PU_2012,Photon_MC),(PU_2012_NW,Top)]
d["data_samples"] = [Had_2012,Muon_2012,Photon_2012]
d["json"] = ["ICHEP_18June.json","Cert_190456-200601_8TeV_PromptReco_Collisions12_JSON_v2.txt"][-1]
checkSwitches(d)
return d
#================= JSON FILTER ========
json = JSONFilter("Json Mask", json_to_pset(switches()["json"]))
# "275_325":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "325_375":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "375_475":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "475_575":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "575_675":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "675_775":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "775_875":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "875":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
#}
#json = JSONFilter("Json Mask", json_to_pset("./Jad_Test.json"))
#json = JSONFilter("Json Mask", json_to_pset("./Golden2011.json"))
#json = JSONFilter("Json Mask", json_to_pset("./Golden29Mar2012.json"))
#json = JSONFilter("Json Mask", json_to_pset("./mergeHTSingleMuDatasets_29032012.json"))
json = JSONFilter("Json Mask", json_to_pset("./json_DCSONLY.txt") )
#json = JSONFilter("Json Mask", json_to_pset(""))
# AlphatTriggerCut(0.52414,50)#
# vertex_reweight = VertexReweighting(
# PSet(
# VertexWeights =[0.20, 0.63, 1.19, 1.57, 1.62, 1.42, 1.09, 0.80 ,0.57, 0.42, 0.30, 0.20]
# # VertexWeights = [0.0, 0.027442995662725636, 0.12682983875287387, 0.28326829632076572, 0.40618954180036759, 0.41605144586432974, 0.33147399297403923, 0.21562021576661147, 0.1140047132529971]
# ).ps())
recHitCut = OP_SumRecHitPtCut(30.)
ZeroMuon = OP_NumComMuons("<=",0)
json_ouput = JSONOutput("filtered")
OneMuon = OP_NumComMuons("==",1)
ZMassCut = RECO_2ndMuonMass(25.0, 91.2, False, "all")
PFMTCut30 = RECO_PFMTCut(30.0)
ScaleSingleJetThreshold=True,
).ps() )
# -----------------------------------------------------------------------------
# Cut flow
cut_flow = Tree(CutFlowOption)
if ( NoCutFlow == True ) :
cut_flow.Attach(qcdBkgdEst)
else :
if ( JsonOption ) : # Use JSON
json = JSONFilter("JSON", json_to_pset(pwd+JsonFileOption))
json_output = JSONOutput("_filtered")
cut_flow.Attach(json)
cut_flow.TAttach(json,json_output)
cut_flow.TAttach(json,selection)
else :
cut_flow.Attach(selection)
if ( DataOption ) : # Data only
cut_flow.TAttach(selection,NoiseFilt)
cut_flow.TAttach(NoiseFilt,selection_duplicate)
else:
cut_flow.TAttach(selection,selection_duplicate)
mht_cut = OP_TriggerMHT_Emu(60.,30.)
"HLT_HT250_AlphaT0p53_v3",
"HLT_HT250_AlphaT0p53_v4",
"HLT_HT250_AlphaT0p53_v5",
"HLT_HT250_AlphaT0p53_v6",
"HLT_HT250_AlphaT0p54_v2",
"HLT_HT250_AlphaT0p54_v3",
"HLT_HT250_AlphaT0p54_v4",
"HLT_HT250_AlphaT0p55_v1",
"HLT_HT250_AlphaT0p62_v1",
"HLT_HT250_AlphaT0p62_v2",
])
DataTrigger = OP_MultiTrigger(datatriggerps.ps())
JetAdd = JetAddition(0.)
json = JSONFilter("Json Mask",
json_to_pset("/home/hep/db1110/public_html/Golden2011.json"))
# AlphatTriggerCut(0.52414,50)#
# vertex_reweight = VertexReweighting(
# PSet(
# VertexWeights =[0.20, 0.63, 1.19, 1.57, 1.62, 1.42, 1.09, 0.80 ,0.57, 0.42, 0.30, 0.20]
# # VertexWeights = [0.0, 0.027442995662725636, 0.12682983875287387, 0.28326829632076572, 0.40618954180036759, 0.41605144586432974, 0.33147399297403923, 0.21562021576661147, 0.1140047132529971]
# ).ps())
# PreScaleWeights = PreScaleReweighting(datatriggerps.ps())
recHitCut = OP_SumRecHitPtCut(30.)
json_ouput = JSONOutput("filtered")
def MakeDataTree(Threshold):
print int(Threshold)
1,1,1,
2,2,2,2,
3,3,3,3,3,
4,4,4,4,4,4,
5,5,5,5,5,5,5,
6,6,6,6,6,6,6,6,
],
)
test = TurnOn(ps.ps())
# -----------------------------------------------------------------------------
# Cut flow
pwd = commands.getoutput('echo $SUSY_WORKING_SW_DIR')
json = JSONFilter("JSON", json_to_pset(pwd+"/allhadronic/python/Cert_160404-163869_7TeV_PromptReco_Collisions11_JSON.txt"))
cut_flow = Tree("CutFlow")
cut_flow.Attach(json)
cut_flow.TAttach(json,selection)
cut_flow.TAttach(selection,NoiseFilt)
cut_flow.TAttach(NoiseFilt,numComPhotons)
cut_flow.TAttach(numComPhotons,oddPhoton)
cut_flow.TAttach(oddPhoton,numComElectrons)
cut_flow.TAttach(numComElectrons,oddElectron)
cut_flow.TAttach(oddElectron,numComMuons)
cut_flow.TAttach(numComMuons,oddMuon)
cut_flow.TAttach(oddMuon,test)
# -----------------------------------------------------------------------------
# Dataset
# "475_575":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "575_675":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "675_775":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "775_875":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
# "875":["HLT_Mu30_v*","HLT_Mu40_eta2p1_v*"],
#}
#json = JSONFilter("Json Mask", json_to_pset("./Jad_Test.json"))
#json = JSONFilter("Json Mask", json_to_pset("./Golden2011.json"))
#json = JSONFilter("Json Mask", json_to_pset("./Golden29Mar2012.json"))
#json = JSONFilter("Json Mask", json_to_pset("./mergeHTSingleMuDatasets_29032012.json"))
#json = JSONFilter("Json Mask", json_to_pset("./Cert_190456-190688_8TeV_PromptReco_Collisions12_JSON.txt") )
#json = JSONFilter("Json Mask", json_to_pset("./Cert_190456-190688_8TeV_PromptReco_Collisions12_JSON_MuonPhys.txt") )
#json = JSONFilter("Json Mask", json_to_pset("json_DCSONLY.txt" ) )
json = JSONFilter("Json Mask", json_to_pset("json_golden2012.txt") )
#json = JSONFilter("Json Mask", json_to_pset(""))
# AlphatTriggerCut(0.52414,50)#
# vertex_reweight = VertexReweighting(
# PSet(
# VertexWeights =[0.20, 0.63, 1.19, 1.57, 1.62, 1.42, 1.09, 0.80 ,0.57, 0.42, 0.30, 0.20]
# # VertexWeights = [0.0, 0.027442995662725636, 0.12682983875287387, 0.28326829632076572, 0.40618954180036759, 0.41605144586432974, 0.33147399297403923, 0.21562021576661147, 0.1140047132529971]
# ).ps())
recHitCut = OP_SumRecHitPtCut(30.)
ZeroMuon = OP_NumComMuons("<=",0)
json_ouput = JSONOutput("filtered")
OneMuon = OP_NumComMuons("==",1)
ZMassCut = RECO_2ndMuonMass(25.0, 91.2, False, "all")
PFMTCut30 = RECO_PFMTCut(30.0)
htTakeMuCut225 = RECO_CommonHTTakeMuCut(225.)
htCut275 = RECO_CommonHTCut(275.)
htCut175 = RECO_CommonHTCut(175.)
htTakeMuCut275 = RECO_CommonHTTakeMuCut(275.)
DeadEcalCutData = OP_DeadECALCut(0.3,0.3,0.5,30.,10,0,"./deadEcal/deadRegionList_GR10_P_V10.txt")
DeadEcalCutMC = OP_DeadECALCut(0.3,0.3,0.5,30.,10,0,"./deadEcal/deadRegionList_START38_V12.txt")
MHT_METCut = OP_MHToverMET(1.25,50.)
MHT_METReverseCut = OP_Reverse_MHToverMET(1.25,50.)
MHT_METCutMC = OP_MHToverMET(1.25,50.)
MET_Filter = OP_METFilters_2012()
json = JSONFilter("Json Mask", json_to_pset("./Jsons/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") )
recHitCut = OP_SumRecHitPtCut(30.)
jet_e2 = OP_NumComJets("==",2)
jet_e3 = OP_NumComJets("==",3)
jet_e4 = OP_NumComJets("==",4)
jet_g3 = OP_NumComJets(">",3)
jet_g4 = OP_NumComJets(">",4)
jet_ge4 = OP_NumComJets(">",3)
jet_le3 = OP_NumComJets("<",4)
jet_ge2 = OP_NumComJets(">",1)
jet_ge4_mu = OP_NumComJets(">",3)
jet_le3_mu = OP_NumComJets("<",4)
do_TP = bTP,
do_RecoilNtuple = False,
rec_Type = 0,
ZDat_File = "./RecoilFit/recoilfits/Recoil_SingleElectron.root",
ZMC_File = "./RecoilFit/recoilfits/Recoil_SingleElectron.root",
W_File = "./RecoilFit/recoilfits/Recoil_SingleElectron.root",
)
Default_Asym=AsymTemplateHistos("Default_Asym",default_Asym.ps())
# Set up CUT flow
a += tree
if bData :
#tree.Attach(Trigg)
#jsonfil=JSONFilter("Cert_160404-167151_7TeV",json_to_pset("/vols/cms02/mjarvis/SUSY2/WCharge/scripts/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON.txt"))
jsonfil=JSONFilter("Cert_160404-167151_7TeV",json_to_pset("/vols/cms02/mjarvis/SUSY2/WCharge/scripts/Cert_160404-167151_7TeV_PromptReco_Collisions11_JSON.txt"))
#tree.TAttach(Trigg,jsonfil)
tree.Attach(jsonfil)
json_output=JSONOutput("filtered")
tree.TAttach(jsonfil,json_output)
tree.TAttach(json_output,ZeroMuons)
else :
tree.Attach(ZeroMuons)
tree.TAttach(ZeroMuons,Default_Asym) #
if bData :
samples =[
SingleElectron_Run2011A_May10ReReco_v1,
SingleElectron_Run2011A_PromptReco_v4_213,
SingleElectron_Run2011A_PromptReco_v4_228,
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())
MinJetPt1=100.,
MinJetPt2=100.,
MinJetPt3=50.,
HtBin=HtBinOption,
Scale=True,
ScaleInSignalRegion=False,
ScaleSingleJetThreshold=True,
).ps())
# -----------------------------------------------------------------------------
# Cut flow
cut_flow = Tree(CutFlowOption)
if (JsonOption): # Use JSON
json = JSONFilter("JSON", json_to_pset(pwd + JsonFileOption))
json_output = JSONOutput("_filtered")
cut_flow.Attach(json)
cut_flow.TAttach(json, json_output)
cut_flow.TAttach(json, selection)
else:
cut_flow.Attach(selection)
if (DataOption): # Data only
cut_flow.TAttach(selection, NoiseFilt)
cut_flow.TAttach(NoiseFilt, selection_duplicate)
else:
cut_flow.TAttach(selection, selection_duplicate)
mht_cut = OP_TriggerMHT_Emu(60., 30.)
oddJet = OP_OddJet()
secondJetET = OP_SecondJetEtCut(100.)
badMuonInJet = OP_BadMuonInJet()
numComLeptons = OP_NumComLeptons("<=",0)
numComPhotons = OP_NumComPhotons("<=",0)
DeadEcalCutData = OP_DeadECALCut(0.3,0.3,0.5,30.,10,0,"./deadRegionList_GR10_P_V10.txt")
# Cross check with the allhadronic analysis
VertexPtOverHT = OP_SumVertexPtOverHT(0.1)
# -----------------------------------------------------------------------------
# Definition of analyses
# Analyses
MHT_METCut = OP_MHToverMET(1.25)
# AK5 Calo
json = JSONFilter("Json Mask", json_to_pset("goldenJSON_13May2011.txt"))
cutTreeData = Tree("Data")
cutTreeData.Attach(json)
cutTreeData.TAttach(json,NoiseFilt)
cutTreeData.TAttach(NoiseFilt,selection)
cutTreeData.TAttach(selection,oddMuon)
cutTreeData.TAttach(oddMuon,oddElectron)
cutTreeData.TAttach(oddElectron,oddPhoton)
cutTreeData.TAttach(oddPhoton,numComLeptons)
cutTreeData.TAttach(numComLeptons,numComPhotons)
cutTreeData.TAttach(numComPhotons,LeadingJetEta)
cutTreeData.TAttach(LeadingJetEta,badMuonInJet)
cutTreeData.TAttach(badMuonInJet,oddJet)
cutTreeData.TAttach(oddJet,LeadingJetCut)
cutTreeData.TAttach(LeadingJetCut,secondJetET)
oneMuon = OP_NumComMuons("==", 1)
diMuon = OP_NumComMuons("==", 2)
zeroMuon = OP_NumComMuons("<=", 0)
ZMassCut = RECO_2ndMuonMass(25.0, 91.2, True, "all")
numComPhotons = OP_NumComPhotons("<=", 0)
muDr = RECO_MuonJetDRCut(0.5)
VertexPtOverHT = OP_SumVertexPtOverHT(0.1)
# -----------------------------------------------------------------------------
# Definition of analyses
# Analyses
MHT_METCut = OP_MHToverMET(1.25, 50.0)
# AK5 Calo
json_ouput = JSONOutput("filtered")
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 * "])
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"]
# 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() )
# cutTreeData.TAttach(json,trigCut)
Bins = [
0
#1,1,1,1,
#1,1,1,1,
#1,1,1,1,
],
)
test = TurnOn(ps.ps())
# -----------------------------------------------------------------------------
# Cut flow
pwd = commands.getoutput('echo $SUSY_WORKING_SW_DIR')
json = JSONFilter("JSON", json_to_pset(pwd+"/hadronic/python/hadronic/ReProcessing_PromptJson_Merged.txt"))
cut_flow = Tree("CutFlow")
cut_flow.Attach(json)
cut_flow.TAttach(json,selection)
cut_flow.TAttach(selection,NoiseFilt)
cut_flow.TAttach(NoiseFilt,numComPhotons)
cut_flow.TAttach(numComPhotons,oddPhoton)
cut_flow.TAttach(oddPhoton,numComElectrons)
cut_flow.TAttach(numComElectrons,oddElectron)
cut_flow.TAttach(oddElectron,numComMuons)
cut_flow.TAttach(numComMuons,oddMuon)
cut_flow.TAttach(oddMuon,test)
# -----------------------------------------------------------------------------
# Dataset
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
"875":["HLT_Mu17_Mu8_v*", ],
}
mu_triggers = {
"275_325":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"325_375":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"375_475":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"475_575":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"575_675":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"675_775":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"775_875":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"875":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
}
json = JSONFilter("Json Mask", json_to_pset("Jsons/json_DCSONLY_23Aug_run2012C.txt") )
recHitCut = OP_SumRecHitPtCut(30.)
ZeroMuon = OP_NumComMuons("<=",0)
json_ouput = JSONOutput("filtered")
OneMuon = OP_NumComMuons("==",1)
ZMassCut = RECO_2ndMuonMass(25.0, 91.2, False, "all")
PFMTCut30 = RECO_PFMTCut(30.0, -1.)
DiMuon = OP_NumComMuons("==",2)
ZMass_2Muons = RECO_2ndMuonMass(25.0, 91.2, True, "OS")
minDRMuonJetCut = RECO_MuonJetDRCut(0.5)
minDRMuonJetCutDiMuon = RECO_MuonJetDRCut(0.5)
Mu45PtCut = OP_LowerMuPtCut(30.0)
Mu50PtCut_HigHT = OP_LowerMuPtCut(50.0)
Mu50PtCut_LowHT275 = OP_LowerMuPtCut(50.0*275./375.)
Mu50PtCut_LowHT325 = OP_LowerMuPtCut(50.0*325./375.)
Mu50PtCut_HigHT_MuTrigPlateau = OP_LowerMuPtCut(50.0)
oneMuon = OP_NumComMuons("==",1)
diMuon = OP_NumComMuons(">=",1)
zeroMuon = OP_NumComMuons("<=",0)
ZMassCut = RECO_2ndMuonMass(25.0, 91.2, True, "all")
numComPhotons = OP_NumComPhotons("<=",0)
muDr = RECO_MuonJetDRCut(0.5)
VertexPtOverHT = OP_SumVertexPtOverHT(0.1)
# -----------------------------------------------------------------------------
# Definition of analyses
# Analyses
MHT_METCut = OP_MHToverMET(1.25,50.)
# AK5 Calo
json_ouput = JSONOutput("filtered")
alphaT = OP_CommonAlphaTCut(0.53)
json = JSONFilter("Json Mask", json_to_pset("./Golden2011.json"))
evDump = EventDump()
json_lost = JSONOutput("lost")
# htTriggerEmu = OP_TriggerHT_Emu(250.,40.)
cutTreeData = Tree("Data")
out = []
badEvents =runLumiCutter(PSet(Run= [162811,162902,162929,163252,163286,163332,163589,163759,163795,163795,163817,165120,165121,165467,165506,165514,165537,165542,165633,165970,166033,166033,166034,166049,166049,166346,166374,166408,166429,166438,166493,166502,166514,166530,166565,166701,166756,166781,166787,166839,166842,166895,166960,167102,167282,167496,167673,167740,167830,167898,167913,167913,170053,170181,170249,170249,170286,170286,175857,175857,175872,175872,175873,175874,175874,175874,175906,175906,175906,175921,175921,175921,175921,175921,175971,175971,175975,175990,175990,175990,176023,176023,176023,176201,176201,176201,176207,176286,176286,176304,176304,176309,176461,176461,176464,176467,176547,176548,176701,176702,176765,176765,176771,176771,176771,176795,176795,176796,176796,176797,176799,176841,176844,176886,176928,176929,176933,176982,177053,177074,177074,177088,177096,177131,177138,177183,177183,177201,177222,177222,177313,177318,177449,177449,177449,177452,177507,177507,177507,177509,177509,177515,177617,177718,177718,177730,177730,177730,177730,177776,177776,177782,177782,177782,177783,177788,177788,177875,177875,178097,178097,178098,178110,178116,178160,178162,178365,178420,178420,178420,178421,178421,178421,178421,178479,178479,178479,178667,178675,178703,178703,178703,178712,178786,178786,178786,178803,178803,178803,178840,178840,178840,178854,178920,178920,178970,178985,179411,179411,179411,179411,179434,179434,179497,179497,179547,179547,179889,179889,179959,180072,180072,180076,180076,180093,180222,180241,180241,180241,180250,],Lumi =[4,25,131,4,108,7,46,182,24,33,118,27,242,113,57,40,13,161,2,9,13,27,203,15,57,7,3,29,4,1,3,2,12,37,12,79,21,2,12,5,7,135,60,6,38,18,2,4,171,4,12,14,9,2,41,45,83,87,25,30,65,67,24,13,33,51,49,51,70,176,182,33,47,56,48,61,179,26,3,56,2,60,65,135,2,40,94,166,62,186,66,223,107,109,9,112,28,121,8,124,110,67,148,64,78,3,42,10,45,197,153,122,61,74,7,171,82,4,42,200,26,2,129,38,71,9,98,63,41,59,62,70,1,235,44,156,42,64,67,110,8,5,3,120,44,213,35,50,55,51,8,160,73,79,242,22,30,148,8,10,7,183,2,180,209,33,209,4,47,53,111,230,34,42,189,3,68,44,11,177,3,31,8,194,20,67,152,2,55,199,2,47,63,207,35,157,6,10,23,6,96,27,83,192,4,233,4,169,18,62,61,80,201,229,98,3,129,65,99,237]).ps())
triggers = ["HLT_DoubleMu8_Mass8_HT200_v4","HLT_DoubleMu8_Mass8_HT200_v5","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu40_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5","HLT_Mu5_HT200_v4","HLT_*","HLT_Mu8_HT200_v3","HLT_Mu8_HT200_v4",]
# cutTreeData.Attach(json)
# for trig in triggers:
# op = JSONOutput("%s"%(trig))
# trigCut = CheckTrigExists( PSet(TrigExistsList = ["%s"%(trig)]).ps() )
# cutTreeData.TAttach(json,trigCut)
# cutTreeData.TAttach(trigCut,op)
# out.append(trigCut)
mu_triggers = {
"275_325":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"325_375":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"375_475":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"475_575":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"575_675":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"675_775":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"775_875":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
"875":["HLT_Mu5_HT200_v3","HLT_Mu5_HT200_v4","HLT_Mu8_HT200_v4","HLT_Mu15_HT200_v2","HLT_Mu15_HT200_v3","HLT_Mu15_HT200_v4","HLT_Mu30_HT200_v1","HLT_Mu30_HT200_v3","HLT_Mu40_HT200_v4","HLT_Mu40_HT300_v4","HLT_Mu40_HT300_v5",],
}
#json = JSONFilter("Json Mask", json_to_pset("./Jad_Test.json"))
json = JSONFilter("Json Mask", json_to_pset("./Golden2011.json"))
#json = JSONFilter("Json Mask", json_to_pset(""))
# AlphatTriggerCut(0.52414,50)#
# vertex_reweight = VertexReweighting(
# PSet(
# VertexWeights =[0.20, 0.63, 1.19, 1.57, 1.62, 1.42, 1.09, 0.80 ,0.57, 0.42, 0.30, 0.20]
# # VertexWeights = [0.0, 0.027442995662725636, 0.12682983875287387, 0.28326829632076572, 0.40618954180036759, 0.41605144586432974, 0.33147399297403923, 0.21562021576661147, 0.1140047132529971]
# ).ps())
zerobtagDiMuon= OP_NumCommonBtagJets("==",0,0.679)
more_than_zero_btagDiMuon= OP_NumCommonBtagJets(">",0,0.679)
onebtagDiMuon= OP_NumCommonBtagJets("==",1,0.679)
more_than_one_btagDiMuon= OP_NumCommonBtagJets(">",1,0.679)
twobtagDiMuon= OP_NumCommonBtagJets("==",2,0.679)
more_than_two_btagDiMuon= OP_NumCommonBtagJets(">",2,0.679)
