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)
