nCores = 8 outputDir = "/raid/raid7/lucien/" + out_path nEvents = -1 disableProgressBar = False justEndSequence = True componentList = sigSamples for dataset in componentList: if dataset.isMC: dataset.lumi = 77.30 for component in dataset.componentList: component.maxEvents = nEvents anaSkimmer = AnalysisSkimmer("AnalysisSkimmer") genWeightCounter = GenWeightCounter("GenWeightCounter") sequence = Sequence() sequence.add(anaSkimmer) sequence.add(genWeightCounter) efficiencyEndModule = EfficiencyEndModule("/home/lucien/public_html/" + out_path) endSequence = EndSequence(skipHadd=justEndSequence) endSequence.add(efficiencyEndModule) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "EventWeight.root"
from HZZ4l.Weighter.FakeRateWeighter import FakeRateWeighter from NanoAOD.Weighter.XSWeighter import XSWeighter # Stealing module from NanoAOD framework WrongFCSkimmer_SR = WrongFCSkimmer("WrongFC-SR",cutflow="SR") WrongFCSkimmer_CR = WrongFCSkimmer("WrongFC-CR",cutflow="CR") WrongFCSkimmer_3P1F = WrongFCSkimmer("WrongFC-CR",cutflow="3P1F") zxSkimmer = ZXCRSkimmer("ZXCRSkimmer") dataMCWeighter = DataMCWeighter("DataMCWeighter") nloWeighter = NLOWeighter("NLOWeighter") xsWeighter = XSWeighter("XSWeighter") fakeRateWeighter = FakeRateWeighter("FakeRateWeighter") wrongFC_sr_sequence = Sequence() wrongFC_sr_sequence.add(WrongFCSkimmer_SR) wrongFC_sr_sequence.add(xsWeighter) wrongFC_sr_sequence.add(dataMCWeighter) wrongFC_sr_sequence.add(fakeRateWeighter) wrongFC_cr_sequence = Sequence() wrongFC_cr_sequence.add(WrongFCSkimmer_CR) wrongFC_cr_sequence.add(xsWeighter) wrongFC_cr_sequence.add(dataMCWeighter) wrongFC_cr_sequence.add(fakeRateWeighter) wrongFC_3p1f_sequence = Sequence() wrongFC_3p1f_sequence.add(WrongFCSkimmer_3P1F) wrongFC_3p1f_sequence.add(xsWeighter) wrongFC_3p1f_sequence.add(dataMCWeighter)
leptonJetProducer = LeptonJetProducer("LeptonJetProducer", "Run2016") #leptonCatProducer = LeptonCatProducer("CategoryProducer") leptonCatProducer = SRProducer("CategoryProducer") nJet40Producer = NJet40Producer("NJet40Producer") leptonProducer = LeptonProducer("LeptonProducer") jetProducer = JetProducer("JetProducer") phoProducer = PhotonProducer("PhotonProducer") variableProducer = VariableProducer("VariableProducer") hltWeighter = HLTWeighter("HLTWeighter") tightLooseHLTWeighter = HLTWeighter("HLTWeighter", cutflow="TightLoose") tightLooseBTagWeighter = BTagSFWeighter("BTagSFWeighter", cutflow="TightLoose", fileTag="_CSVv2_Moriond17_B_H") sr_sequence = Sequence() sr_sequence.add(metSkimmer) sr_sequence.add(leptonProducer) sr_sequence.add(jetProducer) sr_sequence.add(baselineSkimmer) sr_sequence.add(leptonCatProducer) sr_sequence.add(llHtSkimmer) sr_sequence.add(hltSkimmer) sr_sequence.add(xsWeighter) rpv_sequence = Sequence() rpv_sequence.add(leptonProducer) rpv_sequence.add(jetProducer) rpv_sequence.add(metFilter) rpv_sequence.add(baselineSkimmer) rpv_sequence.add(leptonCatProducer)
#from DarkZ.Skimmer.AnalysisSkimmer import AnalysisSkimmer #from DarkZ.Skimmer.BlindSkimmer import BlindSkimmer from DarkZ.Weighter.DataMCWeighter import DataMCWeighter from DarkZ.Weighter.NLOWeighter import NLOWeighter from Wto3l.Skimmer.FinalstateSkimmer import FinalstateSkimmer from Wto3l.Skimmer.EventSkimmer import EventSkimmer from Wto3l.Weighter.FakerateWeighter import FakerateWeighter from Wto3l.Producer.DeltaR_Producer import DeltaR_Producer from NanoAOD.Weighter.XSWeighter import XSWeighter # Stealing module from NanoAOD framework dataMCWeighter = DataMCWeighter("DataMCWeighter") nloWeighter = NLOWeighter("NLOWeighter") xsWeighter = XSWeighter("XSWeighter") finalstateSkimmer = FinalstateSkimmer("FinalstateSkimmer") eventSkimmer = EventSkimmer("EventSkimmer") fakerateWeighter = FakerateWeighter("FakerateWeighter") deltar_producer = DeltaR_Producer("DeltaR_Producer") #bliSkimmer = BlindSkimmer("BlindSkimmer") Wto3l_sequence = Sequence() Wto3l_sequence.add(xsWeighter) Wto3l_sequence.add(dataMCWeighter) #Wto3l_sequence.add(nloWeighter) Wto3l_sequence.add(deltar_producer) Wto3l_sequence.add(finalstateSkimmer) Wto3l_sequence.add(eventSkimmer) Wto3l_sequence.add(fakerateWeighter) #Wto3l_sequence.add(bliSkimmer)
from Core.Sequence import Sequence from Zprime.Skimmer.AnalysisSkimmer import AnalysisSkimmer from Zprime.Producer.VariableProducer import VariableProducer from DarkZ.Weighter.DataMCWeighter import DataMCWeighter from DarkZ.Weighter.NLOWeighter import NLOWeighter from NanoAOD.Weighter.XSWeighter import XSWeighter # Stealing module from NanoAOD framework SRSkimmer = AnalysisSkimmer("SRSkimmer") m4lcrSkimmer = AnalysisSkimmer("m4lcrSkimmer", "m4lCR") dataMCWeighter = DataMCWeighter("DataMCWeighter") nloWeighter = NLOWeighter("NLOWeighter") xsWeighter = XSWeighter("XSWeighter") varProducer = VariableProducer("VariableProducer") signal_sequence = Sequence() signal_sequence.add(SRSkimmer) signal_sequence.add(xsWeighter) signal_sequence.add(dataMCWeighter) signal_sequence.add(nloWeighter) signal_sequence.add(varProducer) m4lcr_sequence = Sequence() m4lcr_sequence.add(m4lcrSkimmer) m4lcr_sequence.add(xsWeighter) m4lcr_sequence.add(dataMCWeighter) m4lcr_sequence.add(nloWeighter) m4lcr_sequence.add(varProducer)
from LJMet.Weighter.XSWeighter import XSWeighter from LJMet.Skimmer.AnalysisSkimmer import AnalysisSkimmer from LJMet.Weighter.DataMCWeighter import DataMCWeighter from LJMet.Producer.CategoryProducer import CategoryProducer from Core.Sequence import Sequence xsWeighter = XSWeighter("XSWeighter") dataMCWeighter = DataMCWeighter("DataMCWeighter") srSkimmer = AnalysisSkimmer("AnalysisSkimmer", cutflow="SR") crSkimmer = AnalysisSkimmer("AnalysisSkimmer", cutflow="CR") preSkimmer = AnalysisSkimmer("AnalysisSkimmer", cutflow="Preselection") catProducer = CategoryProducer("CategoryProducer") sr_sequence = Sequence() sr_sequence.add(srSkimmer) sr_sequence.add(xsWeighter) sr_sequence.add(dataMCWeighter) sr_sequence.add(catProducer) cr_sequence = Sequence() cr_sequence.add(crSkimmer) cr_sequence.add(xsWeighter) cr_sequence.add(dataMCWeighter) cr_sequence.add(catProducer) pre_sequence = Sequence() pre_sequence.add(preSkimmer) pre_sequence.add(xsWeighter) pre_sequence.add(dataMCWeighter) pre_sequence.add(catProducer)
#LambdaFunc("x: (abs(x.eta) < 1.4 and x.mvaID_WP90 > -0.02) or (abs(x.eta) > 1.4 and abs(x.eta) < 2.5 and x.mvaID_WP90 > -0.26)"), LambdaFunc("x: x.mvaID_WP90"), #LambdaFunc("x: x.pfRelIso03_all < 0.35"), LambdaFunc("x: True"), LambdaFunc("x: x.pt > 0. and abs(x.eta) < 2.5"), LambdaFunc( "ev,x: all([x.electronIdx != e.getIndex() for e in ev.selElectrons])"), #LambdaFunc("x: x.cutBasedBitmap == 1"), #LambdaFunc("x: x.pfRelIso03_all < 0.35"), ) variableProducer = VariableProducer("VariableProducer", ) reolvedSRSkimmer = ResolvedSRSkimmer("SRSkimmer", ) mergedSRSkimmer = MergedSRSkimmer("SRSkimmer", ) # ____________________________________________________________________________________________________________________________________________ || sr_resolved_sequence = Sequence() sr_resolved_sequence.add(muProducer) sr_resolved_sequence.add(elProducer) sr_resolved_sequence.add(phoProducer) sr_resolved_sequence.add(reolvedSRSkimmer) sr_resolved_sequence.add(xsWeighter) sr_merged_sequence = Sequence() sr_merged_sequence.add(muProducer) sr_merged_sequence.add(elProducer) sr_merged_sequence.add(phoProducer) sr_merged_sequence.add(mergedSRSkimmer) sr_merged_sequence.add(xsWeighter) # ____________________________________________________________________________________________________________________________________________ ||
), ], ) resonaceZ1Skimmer = NarrowResonanceSkimmer( "NarrowResonanceSkimmer", [ BaseObject( "NarrowResonanceSelection", selFunc=LambdaFunc("x: (abs(x.idL3[0]) == 11 and abs(x.idL4[0]) == 11 and x.massZ1[0] > 8.50 and x.massZ1[0] < 11.0) or (abs(x.idL3[0]) == 13 and abs(x.idL4[0]) == 13 and x.massZ1[0] > 8.50 and x.massZ1[0] < 11.0)"), ), ], ) #____________________________________________________________________________________________________ # Make Sequences darkphoton_signal_sequence = Sequence() darkphoton_signal_sequence.add(blindSkimmer) darkphoton_signal_sequence.add(xsWeighter) darkphoton_signal_sequence.add(nloWeighter) darkphoton_signal_sequence.add(dataMCWeighter) darkphoton_signal_sequence.add(variableProducer) darkphoton_signal_sequence.add(fakeRateWeighter) darkphoton_signal_sequence.add(darkPhotonSRSkimmer) ## Signal-Region Skimmer darkphoton_signal_sequence.add(resonaceZ1Skimmer) darkphoton_signal_sequence.add(resonaceZ2Skimmer) darkphoton_signal_unblind_sequence = Sequence() darkphoton_signal_unblind_sequence.add(xsWeighter) darkphoton_signal_unblind_sequence.add(nloWeighter) darkphoton_signal_unblind_sequence.add(dataMCWeighter) darkphoton_signal_unblind_sequence.add(variableProducer)
), ), ], ) dataMCWeighter = DataMCWeighter("DataMCWeighter") nloWeighter = NLOWeighter("NLOWeighter") xsWeighter = XSWeighter("XSWeighter") bliSkimmer = BlindSkimmer("BlindSkimmer") fakeRateWeighter = FakeRateWeighter("FakeRateWeighter") fakeRateWeighter_map = FakeRateWeighter("FakeRateWeighter", "make_map") fakeRateWeighter_syst = FakeRateWeighter("FakeRateWeighter", "syst") variableProducer = VariableProducer("VariableProducer") darkphoton_signal_sequence = Sequence() darkphoton_signal_sequence.add(darkPhotonSRSkimmer) darkphoton_signal_sequence.add(variableProducer) darkphoton_signal_sequence.add(resonaceSkimmer) darkphoton_fullm4l_sequence = Sequence() darkphoton_fullm4l_sequence.add(darkPhotonFullm4lSkimmer) darkphoton_fullm4l_sequence.add(variableProducer) darkphoton_fullm4l_sequence.add(resonaceSkimmer) zx_map_sequence = Sequence() zx_map_sequence.add(darkPhotonSRSkimmer) zx_map_sequence.add(variableProducer) darkphoton_cr_sequence = Sequence() darkphoton_cr_sequence.add(darkPhotonCRSkimmer)
from Core.Sequence import Sequence from CSC.Producer.StandAloneMuonProducer import StandAloneMuonProducer from CSC.Producer.SegmentProducer import SegmentProducer from CSC.Weighter.EventWeighter import EventWeighter eventWeighter = EventWeighter("EventWeighter") standAloneMuonProducer = StandAloneMuonProducer("StandAloneMuonProducer") segmentProducer = SegmentProducer("SegmentProducer") MuonSequence = Sequence() MuonSequence.add(eventWeighter) MuonSequence.add(standAloneMuonProducer) #MuonSequence.add(segmentProducer)