# 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)