import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter='initialEvents',
pyModuleName=__name__,
pyNameSpace=locals())
#Add smearing
analysisConfigurator.addSmearing('patOverloadedTaus', 'patMuonsForAnalysis',
'convRejElectrons', 'patOverloadedJets',
'METInputTag')
#Create di muon pairs for veto purposes
analysisConfigurator.addDiCandidateModule('diMuons',
'PATMuPairProducer',
'smearedMuons',
'smearedMuons',
'smearedMET',
'smearedJets',
0,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
analysisConfigurator.addSelector(
'osDiMuons', 'PATMuPairSelector',
'charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.userIso(0)))/leg2.pt()<0.2',
'DiMuonCreation', 0, 100)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
analysisConfigurator.addSmearedAll('allMuons','goodMuons','vetoMuons','qcdMuons','allElectrons','goodElectrons','vetoElectrons','qcdElectrons','allJets','cleanJets','goodJets','fwdJets',mpost='')
#analysisConfigurator.addSmearing('patOverloadedTaus','goodMuons','goodElectrons','cleanPatJets','metCorrected',mpost="good")
#analysisConfigurator.addSmearing('patOverloadedTaus','allMuons','allElectrons','cleanPatJets','metCorrected',mpost="all")
#analysisConfigurator.addSmearing('patOverloadedTaus','vetoMuons','vetoElectrons','cleanPatJets','metCorrected',mpost="veto")
#analysisConfigurator.addSmearing('patOverloadedTaus','qcdMuons','qcdElectrons','cleanPatJets','metCorrected',mpost="qcd")
#analysisConfigurator.addSmearing('patOverloadedTaus','selectedPatMuons','selectedPatElectrons','cleanPatJets','metCorrected')
#DiMuonPreSel='leg1.isGlobalMuon && leg1.isTrackerMuon && leg2.isGlobalMuon && leg2.isTrackerMuon && leg1.pt()>20 && leg2.pt()>20 && leg1.charge * leg2.charge<0'
#DiMuonPreSel2='(leg1.isolationR03().sumPt+leg1.isolationR03().emEt+leg1.isolationR03().hadEt)/leg1.pt()<0.15 '
#DiMuonPreSel3='(leg2.isolationR03().sumPt+leg2.isolationR03().emEt+leg2.isolationR03().hadEt)/leg2.pt()<0.15 '
#
##Make Di Muons to VETO DY
##analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuons','smearedMuons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
#analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedGoodMuons','smearedMuonsall','smearedMETall','smearedJetsall',0,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
##analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuonsgood','smearedMuonsgood','smearedMETgood','smearedJetsgood',0,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
#analysisConfigurator.addSelector('diMuonsSel','PATMuPairSelector',DiMuonPreSel+'&&'+DiMuonPreSel2+'&&'+DiMuonPreSel3,'diMuonSel',0,999)
#analysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
#
##Make Muons+MET cleanPatJets
##analysisConfigurator.addCandidateMETModule('MuNuSel','PATMuonNuPairProducer','smearedMuonsgood','smearedMETgood','smearedJetsgood',0,9999,'MuNuSel',genParticles="genDaughters")
##analysisConfigurator.addCandidateMETModule('MuNuVeto','PATMuonNuPairProducer','smearedMuonsveto','smearedMETveto','smearedJetsveto',0,9999,'MuNuVeto',genParticles="genDaughters")
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
zMuMuAnalysisConfigurator = CutSequenceProducer(initialCounter = 'initialZMuMuEvents',
pyModuleName = __name__,
pyNameSpace = locals())
zMuMuAnalysisConfigurator.addDiCandidateModule('zMuMuCandidates','PATMuPairProducer', 'patMuonsForAnalysis','patMuonsForAnalysis','patMETs','patOverloadedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
zMuMuAnalysisConfigurator.addSorter('zMuMuCandidatesSorted','PATMuPairSorter')
zMuMuAnalysisConfigurator.addCandidateMETModule('wCands','PATMuonNuPairProducer','patMuonsForAnalysis','patMETs','selectedPatJets',1,9999,'AtLeastOneWCandidate')
zMuMuAnalysisConfigurator.addDiCandidateModule('zwMuMuMuCandidates','PATMuMuMuNuQuadProducer','zMuMuCandidates','wCands','patMETs','patOverloadedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
zMuMuAnalysisConfigurator.addSelector('zMuMuCandidatesID','PATMuMuMuNuQuadSelector','leg1.leg1.userFloat("isWWMuon")&&leg1.leg2.userFloat("isWWMuon")&&leg1.leg1.pt()>15&&leg1.leg2.pt()>15&&abs(leg1.leg1.eta)<2.4&&abs(leg1.leg2.eta)<2.4','ZMuMuID')
zMuMuAnalysisConfigurator.addSelector('wMuCandidatesID','PATMuMuMuNuQuadSelector','leg2.lepton.userFloat("isWWMuon")&&leg2.lepton.pt()>15&&abs(leg2.lepton.eta)<2.4','wMuMuID')
zMuMuAnalysisConfigurator.addSelector('zMuMuCandidatesIsol1','PATMuMuMuNuQuadSelector','(leg1.leg1.chargedHadronIso()+leg1.leg1.photonIso()+leg1.leg1.neutralHadronIso())/leg1.leg1.pt()<0.25','ZMuMuLeg1Iso')
zMuMuAnalysisConfigurator.addSelector('zMuMuCandidatesIsol2','PATMuMuMuNuQuadSelector','(leg1.leg2.chargedHadronIso()+leg1.leg2.photonIso()+leg1.leg2.neutralHadronIso())/leg1.leg2.pt()<0.25','ZMuMuLeg2Iso')
zMuMuAnalysisConfigurator.addSelector('wMuMuCandidatesIsol2','PATMuMuMuNuQuadSelector','(leg2.lepton.chargedHadronIso()+leg2.lepton.photonIso()+leg2.lepton.neutralHadronIso())/leg2.lepton.pt()<0.25','WLegIso')
zMuMuAnalysisConfigurator.addSelector('zMuMuMuCandidates','PATMuMuMuNuQuadSelector','leg1.charge==0','ZMuMuOS',0,100)
#create the sequence
selectionSequence=zMuMuAnalysisConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
#Add smearing
analysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','patOverloadedJets','METInputTag')
#Create di muon pairs for veto purposes
analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuons','smearedMuons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfigurator.addSelector('osDiMuons','PATMuPairSelector','charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.userIso(0)))/leg2.pt()<0.2','DiMuonCreation',0,100)
analysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
#Make DiTaus
analysisConfigurator.addDiCandidateModule('diTaus','PATMuTauPairProducer','smearedMuons','smearedTaus','smearedMET','smearedJets',1,9999,text='AtLeastOneDiTau',leadingObjectsOnly = False,dR = 0.5,recoMode ="",genParticles='genDaughters')
analysisConfigurator.addSelector('diTausMuonID','PATMuTauPairSelector','leg1.userFloat("isWWMuon")','MuonWW',1)
analysisConfigurator.addSelector('diTausMuonPtEta','PATMuTauPairSelector','leg1.pt()>15&&abs(leg1.eta())<2.1','MuonPtEta',1)
#analysisConfigurator.addSelector('diTausMuonIsolation','PATMuTauPairSelector','(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.1','MuonIsolation',1)
analysisConfigurator.addSelector('diTausTauPtEta','PATMuTauPairSelector','leg2.pt()>20&&abs(leg2.eta())<2.3','TauPtEta',1)
analysisConfigurator.addSelector('diTausDecayFound','PATMuTauPairSelector','leg2.tauID("decayModeFinding")>0','TauDecayFound',1)
analysisConfigurator.addSelector('diTausVLooseTauIsolation','PATMuTauPairSelector','(leg2.tauID("byVLooseCombinedIsolationDeltaBetaCorr")>0)','TauVLooseIsoDB',1)
#analysisConfigurator.addSelector('diTausLooseTauIsolation','PATMuTauPairSelector','(leg2.tauID("byLooseCombinedIsolationDeltaBetaCorr")>0)','TauLooseIsoDB',1)
analysisConfigurator.addSelector('diTausTauElectronVeto','PATMuTauPairSelector','leg2.tauID("againstElectronLoose")','AgainstElectron',1)
analysisConfigurator.addMuTauSVFitSA('muTausSVFit')
analysisConfigurator.addSelector('diTausTauMuonVeto','PATMuTauPairSelector','leg2.tauID("againstMuonTight")','AgainstMuon',1)
analysisConfigurator.addSelector('diTausWRej','PATMuTauPairSelector','mt1MET<40','WRejection',1)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
############################### Ele-Tau ###################################
ETanalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsET',
pyModuleName = __name__,
pyNameSpace = locals())
#ESTausID
ETanalysisConfigurator.addSmearing('patOverloadedTaus','triggeredPatMuons','triggeredPatElectrons','cleanPatJets','patMVAMet','ET')
#create dielectrons
ETanalysisConfigurator.addDiCandidateModule('diElectrons','PATElePairProducer','smearedElectronsET','smearedElectronsET','smearedMETET','','smearedJetsET',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
ETanalysisConfigurator.addSelector('osDiElectrons','PATElePairSelector','abs(leg1.eta())<2.5&&abs(leg2.eta())<2.5&&abs(leg1.userFloat("ipDXY"))<0.045&&abs(leg2.userFloat("ipDXY"))<0.045&&abs(leg1.userFloat("dz"))<0.2&&abs(leg2.userFloat("dz"))<0.2&&leg1.pt()>15&&leg2.pt()>15&&charge==0&&leg2.userFloat("wp95V")>0&&leg1.userFloat("wp95V")>0&&(leg1.userIso(0)+max(leg1.userIso(1)+leg1.neutralHadronIso()-0.5*leg1.userIso(2),0.0))/leg1.pt()<0.3&&(leg2.userIso(0)+max(leg2.userIso(1)+leg2.neutralHadronIso()-0.5*leg2.userIso(2),0.0))/leg2.pt()<0.3','ZEEVEto',0,100)
#Make DiTaus
ETanalysisConfigurator.addDiCandidateModule('eleTaus','PATEleTauPairProducer','smearedElectronsET','smearedTausET','patMVAMet','smearedTausET','smearedJetsET',1,9999,text = 'AtLeastOneEleTau',leadingObjectsOnly = False,dR = 0.5,recoMode = "",genParticles='genDaughters')
ETanalysisConfigurator.addSelector('eleTausEleID',
'PATEleTauPairSelector',
'(abs(leg1.superCluster().eta())<0.8&&leg1.electronID("mvaNonTrigV0")>0.925)||(abs(leg1.superCluster().eta())<1.479&&leg1.electronID("mvaNonTrigV0")>0.975)||(leg1.electronID("mvaNonTrigV0")>0.985)',
'ElectronID',1)
ETanalysisConfigurator.addSelector('eleTausEleConvRej','PATEleTauPairSelector','leg1.userInt("missingHits")==0&&!(leg1.userInt("HasMatchedConversion")>0)','electronConvRej',1)
ETanalysisConfigurator.addSelector('eleTausElePtEta',
'PATEleTauPairSelector',
'leg1.pt()>20&&abs(leg1.eta())<2.1&&abs(leg1.userFloat("dz"))<0.2&&abs(leg1.userFloat("ipDXY"))<0.045','electronPtEta',1)
ETanalysisConfigurator.addSelector('eleTausEleIsolationLoose','PATEleTauPairSelector','(leg1.userIso(0)+max(leg1.userIso(1)+leg1.neutralHadronIso()-0.5*leg1.userIso(2),0.0))/leg1.pt()<0.5','electronIsolationLoose',1)
ETanalysisConfigurator.addSelector('eleTausTauPtEta',
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
llLeadingIso1='(leg1.pt()>20&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0),0.0))/leg1.pt()<0.25)'
llLeadingIso2='(leg2.pt()>10&&(leg2.chargedHadronIso()+max(leg2.photonIso()+leg2.neutralHadronIso()-0.5*leg2.userIso(0),0.0))/leg2.pt()<0.25)'
llLeadingIso=llLeadingIso1+'&&'+llLeadingIso2
mumuLeadingID='leg1.isGlobalMuon()&&leg2.isGlobalMuon()&&leg1.isTrackerMuon()&&leg2.isTrackerMuon()'
mumuLeading=mumuLeadingID+'&& mass>60 && mass<120 && abs(leg1.eta())<2.4 && abs(leg2.eta())<2.4'
eeLeading='charge==0 && leg1.electronID("cicTight")&1==1 && leg2.electronID("cicTight")&1==1 && mass>60 && mass<120 && abs(leg1.eta())<2.5 && abs(leg2.eta())<2.5'
EEanalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsEE',
pyModuleName = __name__,
pyNameSpace = locals())
EEanalysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','selectedPatJets','EE')
EEanalysisConfigurator.addDiCandidateModule('EEdiElectrons','PATElePairProducer', 'smearedElectronsEE','smearedElectronsEE','smearedMETEE','smearedJetsEE',1,genParticles='genDaughters')
EEanalysisConfigurator.addSorter('zEleEleCandidatesSorted','PATElePairSorter')
EEanalysisConfigurator.addSelector('EEosDiElectrons','PATElePairSelector','leg1.gsfTrack().trackerExpectedHitsInner().numberOfHits()<2&&leg2.gsfTrack().trackerExpectedHitsInner().numberOfHits()<2&&(leg1.electronID("cicTight")==1||leg1.electronID("cicTight")==3||leg1.electronID("cicTight")==5||leg1.electronID("cicTight")==7||leg1.electronID("cicTight")==9||leg1.electronID("cicTight")==11||leg1.electronID("cicTight")==13||leg1.electronID("cicTight")==15)&&(leg2.electronID("cicTight")==1||leg2.electronID("cicTight")==3||leg2.electronID("cicTight")==5||leg2.electronID("cicTight")==7||leg2.electronID("cicTight")==9||leg2.electronID("cicTight")==11||leg2.electronID("cicTight")==13||leg2.electronID("cicTight")==15)&&charge==0&&abs(leg1.eta())<2.5&&abs(leg2.eta())<2.5&&mass>40 && leg1.userFloat("SIP3D")<4 && leg2.userFloat("SIP3D")<4','EEDiEleCreation',1)
EEanalysisConfigurator.addSelector('EEeleIso','PATElePairSelector',llLeadingIso,'EELeadingZEleIso')
EEanalysisConfigurator.addSelector('EEaboveThresh','PATElePairSelector','leg1.pt()>20 && leg2.pt()>10','EEAtLeastOneZZCandOverThresholds')
zEleEleSelectionSequence =EEanalysisConfigurator.returnSequence()
MManalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsMM',
pyModuleName = __name__,
pyNameSpace = locals())
#Add smearing
import FWCore.ParameterSet.Config as cms
# Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
############################### Ele-Tau ###################################
ETanalysisConfigurator = CutSequenceProducer(
initialCounter="initialEventsET", pyModuleName=__name__, pyNameSpace=locals()
)
ETanalysisConfigurator.addSmearing(
"triggeredPatTaus", "triggeredPatMuons", "triggeredPatElectrons", "cleanPatJets", "patMVAMetEleTau", "ET"
)
# create dielectrons
ETanalysisConfigurator.addDiCandidateModule(
"diElectrons",
"PATElePairProducer",
"smearedElectronsET",
"smearedElectronsET",
"smearedMETET",
"smearedJetsET",
0,
9999,
text="",
leadingObjectsOnly=False,
dR=0.15,
recoMode="",
)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
tagAndProbeConfigurator = CutSequenceProducer(initialCounter='initialEvents',
pyModuleName=__name__,
pyNameSpace=locals())
##Tag and probe electron+GSF track
tagAndProbeConfigurator.addDiCandidateModule('tagAndProbeEleTracks',
'PATEleTrackPairProducer',
'convRejElectrons',
'gsfTrackCandidates',
'patMETs',
'selectedPatJets',
1,
9999,
text='',
leadingObjectsOnly=True,
dR=0.15,
recoMode="")
tagAndProbeConfigurator.addSorter('eleTracksSorted', 'PATEleTrackPairSorter')
tagAndProbeConfigurator.addSelector(
'selectedTagAndProbeEleTracks', 'PATEleTrackPairSelector',
'abs(leg2.eta())<2.1&&(!(abs(leg2.eta())>1.442&&abs(leg2.eta())<1.566))&&leg1.pt()>15&&charge==0&&leg2.pt()>5&&mass>15&&mass<120&&leg1.userFloat("WWID")>0&&(leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.1',
'tagAndProbeEleTracks', 0, 100)
tagAndProbeSequenceEleTracks = tagAndProbeConfigurator.returnSequence()
##Tag and probe electron+SC
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter='initialEvents',
pyModuleName=__name__,
pyNameSpace=locals())
DiMuonPreSel = 'leg1.isGlobalMuon && leg1.isTrackerMuon && leg2.isGlobalMuon && leg2.isTrackerMuon && leg1.pt()>20 && leg2.pt()>20'
DiMuonPreSel2 = '(leg1.isolationR03().sumPt+leg1.isolationR03().emEt+leg1.isolationR03().hadEt)/leg1.pt()<0.15 '
DiMuonPreSel3 = '(leg2.isolationR03().sumPt+leg2.isolationR03().emEt+leg2.isolationR03().hadEt)/leg2.pt()<0.15 '
#Make Di Muons to VETO DY
analysisConfigurator.addDiCandidateModule('diMuons',
'PATMuPairProducer',
'cleanPatMuons',
'cleanPatMuons',
'systematicsMET',
'cleanPatJets',
0,
9999,
text='',
leadingObjectsOnly=True,
dR=0.15,
recoMode="")
analysisConfigurator.addSelector(
'diMuonsSel', 'PATMuPairSelector',
DiMuonPreSel + '&&' + DiMuonPreSel2 + '&&' + DiMuonPreSel3, 'diMuonSel', 0,
999)
analysisConfigurator.addSorter('diMuonsSorted', 'PATMuPairSorter')
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
tagAndProbeConfigurator = CutSequenceProducer(
initialCounter='initialTagAndProbeEvents',
pyModuleName=__name__,
pyNameSpace=locals())
##Tag and probe
tagAndProbeConfigurator.addDiCandidateModule('tagAndProbeMuTracks',
'PATMuTrackPairProducer',
'patMuonsForAnalysis',
'trackCandidates',
'patMETs',
'selectedPatJets',
1,
9999,
text='',
leadingObjectsOnly=True,
dR=0.15,
recoMode="")
tagAndProbeConfigurator.addSelector(
'selectedTagAndProbeMuTracks', 'PATMuTrackPairSelector',
'calibratedMET.pt()<20 && charge==0&&mass>50&&mass<120&&leg1.userFloat("isWWMuon")&&((leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.2)&&leg1.pt()>25&&leg2.pt()>5&&abs(leg2.eta())<2.4&&abs(leg1.eta())<2.1',
'tagAndProbePairs', 1, 100)
tagAndProbeSequence = tagAndProbeConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
zEleEleAnalysisConfigurator = CutSequenceProducer(initialCounter = 'initialZEleEleEvents',
pyModuleName = __name__,
pyNameSpace = locals())
zEleEleAnalysisConfigurator.addDiCandidateModule('zEleEleCandidates','PATElePairProducer', 'convRejElectrons','convRejElectrons','patMETs','selectedPatJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
zEleEleAnalysisConfigurator.addSorter('zEleEleCandidatesSorted','PATElePairSorter')
zEleEleAnalysisConfigurator.addSelector('zEleEleCandidatesID','PATElePairSelector','leg1.userFloat("WWID")>0&&leg2.userFloat("WWID")>0 &&leg1.pt()>15&&leg2.pt()>15','ZEleEleID')
# zEleEleAnalysisConfigurator.addEleEleSVFit('zEleEleCandidatesSVFit')
zEleEleAnalysisConfigurator.addSelector('zEleEleCandidatesIsol1','PATElePairSelector','leg1.dr03TkSumPt()/leg1.pt()<0.25','ZEleEleLeg1Iso')
zEleEleAnalysisConfigurator.addSelector('zEleEleCandidatesIsol2','PATElePairSelector','leg1.dr03TkSumPt()/leg2.pt()<0.25','ZEleEleLeg2Iso')
zEleEleAnalysisConfigurator.addSelector('zEleEleCandidatesOS','PATElePairSelector','charge==0','ZEleEleOS',0,100)
#create the sequence
zEleEleSelectionSequence=zEleEleAnalysisConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
analysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','patOverloadedJets','METInputTag')
#create dielectrons
analysisConfigurator.addDiCandidateModule('diElectrons','PATElePairProducer', 'smearedElectrons','smearedElectrons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfigurator.addSelector('osDiElectrons','PATElePairSelector','leg1.pt()>15&&leg2.pt()>15&&charge==0&&(leg1.userFloat("WWID")>0||leg1.userFloat("WWMVAID")>0)&&leg2.userFloat("wp95")>0&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0),0.0))/leg1.pt()<0.25&&(leg2.chargedHadronIso()+max(leg2.photonIso()+leg2.neutralHadronIso()-0.5*leg2.userIso(0),0.0))/leg2.pt()<0.25','ZEEVEto',0,100)
#Make DiTaus
analysisConfigurator.addDiCandidateModule('eleTaus','PATEleTauPairProducer', 'smearedElectrons','smearedTaus','smearedMET','smearedJets',1,9999,text = 'AtLeastOneEleTau',leadingObjectsOnly = False,dR = 0.5,recoMode = "",genParticles='genDaughters')
analysisConfigurator.addSelector('eleTausWW','PATEleTauPairSelector','leg1.userFloat("WWID")>0||leg1.userFloat("WWMVAID")>0','WWElectronID',1)
analysisConfigurator.addSelector('eleTausElePtEta','PATEleTauPairSelector','leg1.pt()>20&&abs(leg1.eta())<2.1','electronPtEta',1)
#analysisConfigurator.addSelector('eleTausEleIsolation','PATEleTauPairSelector','(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.particleIso(),0.0))/leg1.pt()<0.1','electronIsolation',1)
analysisConfigurator.addSelector('eleTausTauPtEta','PATEleTauPairSelector','leg2.tauID("decayModeFinding")&&leg2.pt()>20&&abs(leg2.eta())<2.3','TauDecayFound',1)
analysisConfigurator.addSelector('eleTausVLooseIsolation','PATEleTauPairSelector','leg2.tauID("byVLooseCombinedIsolationDeltaBetaCorr")>0','TauLooseIso',1)
#analysisConfigurator.addSelector('eleTausLooseIsolation','PATEleTauPairSelector','leg2.tauID("byLooseCombinedIsolationDeltaBetaCorr")>0','TauLooseIso',1)
analysisConfigurator.addSelector('eleTausTauMuonVeto','PATEleTauPairSelector','leg2.tauID("againstMuonLoose")','AgainstMuon',1)
analysisConfigurator.addSelector('eleTausTauElectronVeto','PATEleTauPairSelector','leg2.tauID("againstElectronMedium")','AgainstElectron',1)
analysisConfigurator.addEleTauSVFitSA('eleTausNSVFit')
import FWCore.ParameterSet.Config as cms
###MU ENRICHED QCD#####################################################################################################
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
zJetsConfigurator = CutSequenceProducer(initialCounter='initialZJetsEvents',
pyModuleName=__name__,
pyNameSpace=locals())
zJetsConfigurator.addDiCandidateModule('diMuonszJetsFR',
'PATMuPairProducer',
'patMuonsForAnalysis',
'patMuonsForAnalysis',
'patMETs',
'selectedPatJets',
1,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
zJetsConfigurator.addSelector(
'osDiMuonsZJetsFR', 'PATMuPairSelector',
'charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&((leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.2)&&((leg2.chargedHadronIso+leg2.photonIso+leg2.neutralHadronIso)/leg2.pt()<0.2)&&mass>35&&mass<120',
'', 0, 100)
zJetsConfigurator.addGeneric('zJetsForFakeRate', 'MuonPairJetFetcher')
zJetsSequence = zJetsConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
###MU ENRICHED QCD#####################################################################################################
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
muEnrichedConfigurator = CutSequenceProducer(initialCounter = 'initialMuEnrichedEvents',
pyModuleName = __name__,
pyNameSpace = locals())
muEnrichedConfigurator.addDiCandidateModule('diMuonsMuEnrichedFR','PATMuPairProducer', 'patMuonsForAnalysis','patMuonsForAnalysis','patMETs','selectedPatJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
muEnrichedConfigurator.addSelector('osDiMuonsMuEnrichedFR','PATMuPairSelector','charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&((leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.26)&&((leg2.chargedHadronIso+leg2.photonIso+leg2.neutralHadronIso)/leg2.pt()<0.3)','',0,0)
#Make Muons+MET
muEnrichedConfigurator.addCandidateMETModule('muEnrichedCands','PATMuonNuPairProducer','patMuonsForAnalysis','patMETs','selectedPatJets',0,9999,'')
muEnrichedConfigurator.addSelector('muEnrichedSelectedCands','PATMuonNuPairSelector','lepton.pt()>15&&abs(lepton.eta())<2.1&&((lepton.chargedHadronIso+lepton.photonIso+lepton.neutralHadronIso)/lepton.pt()>0.5)&&(lepton.triggerObjectMatchesByPath("HLT_Mu9").size>0||(lepton.triggerObjectMatchesByPath("HLT_Mu11").size>0)||(lepton.triggerObjectMatchesByPath("HLT_Mu15_v1").size>0))&&mt<40.','',0,100)
muEnrichedConfigurator.addGeneric('muEnrichedJetsForFakeRate','MuonNuPairJetFetcher')
muEnrichedSequence = muEnrichedConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
tagAndProbeConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
##Tag and probe electron+GSF track
tagAndProbeConfigurator.addDiCandidateModule('tagAndProbeEleTracks','PATEleTrackPairProducer', 'convRejElectrons','gsfTrackCandidates','patMETs','selectedPatJets',1,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
tagAndProbeConfigurator.addSorter('eleTracksSorted','PATEleTrackPairSorter')
tagAndProbeConfigurator.addSelector('selectedTagAndProbeEleTracks','PATEleTrackPairSelector','abs(leg2.eta())<2.1&&(!(abs(leg2.eta())>1.442&&abs(leg2.eta())<1.566))&&leg1.pt()>15&&charge==0&&leg2.pt()>5&&mass>15&&mass<120&&leg1.userFloat("WWID")>0&&(leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.1','tagAndProbeEleTracks',0,100)
tagAndProbeSequenceEleTracks =tagAndProbeConfigurator.returnSequence()
##Tag and probe electron+SC
tagAndProbeConfigurator2 = CutSequenceProducer(initialCounter = 'initialTagAndProbeEvents2',
pyModuleName = __name__,
pyNameSpace = locals())
##Tag and probe electron+SC
tagAndProbeConfigurator2.addDiCandidateModule('tagAndProbeEleSCs','PATEleSCPairProducer', 'convRejElectrons','superClusterCands','patMETs','selectedPatJets',1,9999,text = '',leadingObjectsOnly = True,dR = 0.5,recoMode = "")
tagAndProbeConfigurator.addSorter('eleSCsSorted','PATEleSCPairSorter')
tagAndProbeConfigurator2.addSelector('selectedTagAndProbeEleSCs','PATEleSCPairSelector','calibratedMET.pt()<20 && abs(leg1.eta())<2.5&&abs(leg2.eta())<2.5&&leg1.pt()>20&&leg2.pt()>5&&mass>50&&mass<120&&leg1.userFloat("WWID")>0&&(leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso())/leg1.pt()<0.1','tagAndProbeEleSCs',0,100)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
######################__________________________________MMMT_____________________________________##############################
MMMTanalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsMMMT',
pyModuleName = __name__,
pyNameSpace = locals())
#nSVfitTrackService = cms.Service("NSVfitTrackService")
#Add smearing
MMMTanalysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','selectedPatJets')
MMMTanalysisConfigurator.addDiCandidateModule('MMMTdiMuons','PATMuPairProducer', 'smearedMuons','smearedMuons','smearedMET','smearedJets',1,genParticles='genDaughters')
MMMTanalysisConfigurator.addSelector('MMMTosDiMuons','PATMuPairSelector','charge==0&&leg1.isGlobalMuon()&&leg2.isGlobalMuon()&& mass>40 && abs(leg1.eta())<2.4 && abs(leg2.eta())<2.4','MMMTDiMuonCreation',1)
MMMTanalysisConfigurator.addDiCandidateModule('MMMTmuTau','PATMuTauPairProducer', 'smearedMuons','smearedTaus','smearedMET','smearedJets',1,9999,text = 'MMMTAtLeastOneMuTau',leadingObjectsOnly = False,dR = 0.3,recoMode = "",genParticles='genDaughters')
MMMTanalysisConfigurator.addDiCandidateModule('MMMTzzCands','PATMuMuMuTauQuadProducer','MMMTosDiMuons','MMMTmuTau','smearedMET','smearedJets',1,9999,text='MMMTAtLeastOneZZ',leadingObjectsOnly = False,dR = 0.01,recoMode ="",genParticles='genDaughters')
MMMTanalysisConfigurator.addCrossCleanerModule('MMMTzzCleanedCands','PATMuMuMuTauQuadCrossCleaner',1,9999,text='MMMTAtLeastOneZZCleanedCandidate',dR = 0.1)
MMMTanalysisConfigurator.addSelector('MMMTzzMuID','PATMuMuMuTauQuadSelector','leg1.leg1.isGlobalMuon()&&leg1.leg1.isTrackerMuon()&&leg1.leg2.isTrackerMuon()&&leg1.leg2.isTrackerMuon()','MMMTLeadingZMuID')
MMMTanalysisConfigurator.addSelector('MMMTzzMuIso','PATMuMuMuTauQuadSelector','(leg1.leg1.chargedHadronIso()+max(leg1.leg1.photonIso()+leg1.leg1.neutralHadronIso()-0.5*leg1.leg1.userIso(0),0.0))/leg1.leg1.pt()<0.25&&(leg1.leg2.chargedHadronIso()+max(leg1.leg2.photonIso()+leg1.leg2.neutralHadronIso()-0.5*leg1.leg2.userIso(0),0.0))/leg1.leg2.pt()<0.25','MMMTLeadingZMuIso')
MMMTanalysisConfigurator.addSelector('MMMTzzMuID2','PATMuMuMuTauQuadSelector','leg2.leg1.isGlobalMuon()&&leg2.leg1.isTrackerMuon()','MMMTSecondZMuID')
MMMTanalysisConfigurator.addSelector('MMMTzzTauID','PATMuMuMuTauQuadSelector','leg2.leg2.tauID("decayModeFinding")','MMMTTauDecayModeFinding')
MMMTanalysisConfigurator.addSelector('MMMTzzTauDiscr','PATMuMuMuTauQuadSelector','leg2.leg2.tauID("againstElectronLoose")&&leg2.leg2.tauID("againstMuonTight")','MMMTTauLeptonDiscrimantor')
MMMTanalysisConfigurator.addSelector('MMMTzzMuIso2','PATMuMuMuTauQuadSelector','(leg2.leg1.chargedHadronIso()+max(leg2.leg1.photonIso()+leg2.leg1.neutralHadronIso()-0.5*leg2.leg1.userIso(0),0.0))/leg2.leg1.pt()<0.15','MMMTSecondZMuIso')
MMMTanalysisConfigurator.addSelector('MMMTzzTauIso','PATMuMuMuTauQuadSelector','leg2.leg2.tauID("byLooseCombinedIsolationDeltaBetaCorr")','MMMTTauLooseIsolation')
MMMTanalysisConfigurator.addSorter('MMMTzzCleanedCandsSortedByZMass','PATMuMuMuTauQuadSorterByZMass')
MMMTanalysisConfigurator.addSelector('MMMTzzCleanedCandsAboveThreshold','PATMuMuMuTauQuadSelector','leg1().leg1().pt()>20 && leg1().leg2().pt()>10 && leg2().leg1().pt()>10 &&leg2().leg2().pt()>20 && abs(leg2.leg1.eta())<2.4 && abs(leg2.leg2.eta())<2.3','MMMTAtLeastOneZZCandOverThresholds')
MMMTanalysisConfigurator.addSelector('MMMTzzCleanedCandsMuTauQ','PATMuMuMuTauQuadSelector','leg2.charge()==0','MMMTMuTauCharge')
MMMTanalysisConfigurator.addSelector('MMMTzzCleanedCandsMMMass','PATMuMuMuTauQuadSelector','leg1.mass()>60&&leg1.mass()<120','MMMTMMMass')
MMMTanalysisConfigurator.addSorter('MMMTFinalSel','PATMuMuMuTauQuadSorterByZMass')
MMMTanalysisConfigurator.addSelector('MMMTzzCleanedCandsMTMass','PATMuMuMuTauQuadSelector','leg2.mass()>30&&leg2.mass()<80','MMMTMTMass')
munuKIN='lepton.pt()>25 && abs(lepton.eta())<2.1&&!lepton.pfCandidateRef().isNull()'
munuSel='mt>40&&lepton.userInt("WWID2011")==1&&abs(lepton.userFloat("ipDXY"))<0.02'#(lepton.isolationR03().sumPt+lepton.isolationR03().emEt+lepton.isolationR03().hadEt)/lepton.pt()<0.3&&met.pt()>30&mt>40&&lepton.userFloat("isWWMuon")==1'
elenuKIN='lepton.pt()>30 && abs(lepton.eta())<2.5&&(abs(lepton.eta)<1.4442||abs(lepton.eta)>1.5666)'
elenuSel='mt>40&&(lepton.isEB &&(lepton.sigmaIetaIeta<0.01)&&(-0.8<lepton.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<lepton.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (lepton.isEE &&(lepton.sigmaIetaIeta<0.03)&&(-0.7<lepton.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<lepton.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel='leg1.pt>20&&leg2.pt>20'
DiElePreSel2='(leg1.isEB &&(leg1.sigmaIetaIeta<0.01)&&(-0.8<leg1.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg1.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg1.isEE &&(leg1.sigmaIetaIeta<0.03)&&(-0.7<leg1.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg1.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel2='(leg2.isEB &&(leg2.sigmaIetaIeta<0.01)&&(-0.8<leg2.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg2.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg2.isEE &&(leg2.sigmaIetaIeta<0.03)&&(-0.7<leg2.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg2.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel3='(leg1.chargedHadronIso+leg1.photonIso()+leg1.neutralHadronIso())/leg1.pt()<0.3&&(leg2.chargedHadronIso+leg2.photonIso()+leg2.neutralHadronIso())/leg2.pt()<0.3'
######################__________________________________MNJ_____________________________________##############################
MNJanalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsMNJ',
pyModuleName = __name__,
pyNameSpace = locals())
MNJanalysisConfigurator.addSmearing('cleanPatTaus','cleanPatMuons','cleanPatElectrons','patJetsForAnalysis','systematicsMET','default','MNJ')
#Make Muons+MET
MNJanalysisConfigurator.addCandidateMETModule('wCands','PATMuonNuPairProducer','smearedMuonsMNJ','smearedMETMNJ','smearedJetsMNJ',1,9999,'AtLeastOneWCandidate',genParticles="genDaughters")
#MNJanalysisConfigurator.addSelector('wCandsTight','PATMuonNuPairSelector',munuTight,'tightID',1)
MNJanalysisConfigurator.addSelector('wCandsKIN','PATMuonNuPairSelector',munuKIN,'wCandsKIN',1)
MNJanalysisConfigurator.addSelector('wCandsSel','PATMuonNuPairSelector',munuSel,'wCandsSel',1)
MNJanalysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuonsMNJ','cleanPatMuons','smearedMETMNJ','smearedJetsMNJ',0,9999,text = '',dR=0.5,recoMode = "")
MNJanalysisConfigurator.addSelector('diMuonsSel','PATMuPairSelector',DiMuonPreSel+'&&'+DiMuonPreSel2+'&&'+DiMuonPreSel3,'diMuonSel',0,0)
MNJanalysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
############################### Ele-Tau ###################################
ETanalysisConfigurator = CutSequenceProducer(initialCounter='initialEventsET',
pyModuleName=__name__,
pyNameSpace=locals())
#ESTausID
ETanalysisConfigurator.addSmearing('patOverloadedTaus', 'triggeredPatMuons',
'triggeredPatElectrons', 'cleanPatJets',
'patMVAMet', 'ET')
#create dielectrons
ETanalysisConfigurator.addDiCandidateModule('diElectrons',
'PATElePairProducer',
'smearedElectronsET',
'smearedElectronsET',
'smearedMETET',
'',
'smearedJetsET',
0,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
ETanalysisConfigurator.addSelector(
'osDiElectrons', 'PATElePairSelector',
elenuSelB = "(lepton.chargedHadronIso+max(lepton.photonIso()+lepton.neutralHadronIso()-0.5*lepton.userIso(0),0.0))/lepton.pt()<0.15"
elenuSel = elenuSelA + "&&" + elenuSelB
DiElePreSel = "leg1.pt>20&&leg2.pt>20"
DiElePreSel2 = "(leg1.isEB &&(leg1.sigmaIetaIeta<0.01)&&(-0.8<leg1.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg1.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg1.isEE &&(leg1.sigmaIetaIeta<0.03)&&(-0.7<leg1.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg1.deltaEtaSuperClusterTrackAtVtx<0.01) )"
DiElePreSel2 = "(leg2.isEB &&(leg2.sigmaIetaIeta<0.01)&&(-0.8<leg2.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg2.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg2.isEE &&(leg2.sigmaIetaIeta<0.03)&&(-0.7<leg2.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg2.deltaEtaSuperClusterTrackAtVtx<0.01) )"
DiElePreSel3 = "(leg1.chargedHadronIso+leg1.photonIso()+leg1.neutralHadronIso())/leg1.pt()<0.3&&(leg2.chargedHadronIso+leg2.photonIso()+leg2.neutralHadronIso())/leg2.pt()<0.3"
jetKIN = "leg1.pt>20 && leg2.pt> 20 && abs(leg1.eta)<5 && abs(leg2.eta)<5"
######################__________________________________MNBB_____________________________________##############################
MNBBanalysisConfigurator = CutSequenceProducer(
initialCounter="initialEventsMNBB", pyModuleName=__name__, pyNameSpace=locals()
)
# Make Di Muons to VETO DY
MNBBanalysisConfigurator.addDiCandidateModule(
"diMuons",
"PATMuPairProducer",
"cleanPatMuons",
"cleanPatMuons",
"systematicsMET",
"patJetsForAnalysis",
0,
9999,
text="",
dR=0.5,
#tagAndProbeConfigurator2 = CutSequenceProducer(initialCounter = 'initialTagAndProbeEvents2',
# pyModuleName = __name__,
# pyNameSpace = locals())
##Tag and probe electron+SC
#tagAndProbeConfigurator2.addDiCandidateModule('tagAndProbeEleSCs','PATEleSCPairProducer', 'triggeredPatElectrons','superClusterCands','systematicsMET','cleanPatJets',1,9999,text = '',leadingObjectsOnly = True,dR = 0.5,recoMode = "")
#tagAndProbeConfigurator2.addSorter('eleSCsSorted','PATEleSCPairSorter')
#tagAndProbeConfigurator2.addSelector('selectedTagAndProbeEleSCs','PATEleSCPairSelector','calibratedMET.pt()<20 && abs(leg1.eta())<2.5&&abs(leg2.eta())<2.5&&leg1.pt()>20&&leg2.pt()>5&&mass>50&&mass<120&&leg1.userFloat("WWID")>0&&(leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso())/leg1.pt()<0.1','tagAndProbeEleSCs',0,100)
#tagAndProbeSequenceEleSCs =tagAndProbeConfigurator2.returnSequence()
#Tag and probe electron+electron
tagAndProbeConfigurator3 = CutSequenceProducer(initialCounter = 'initialTagAndProbeEvents3',
pyModuleName = __name__,
pyNameSpace = locals())
tagAndProbeConfigurator3.addDiCandidateModule('tagAndProbeDiElectrons','PATElePairProducer', 'triggeredPatElectrons','triggeredPatElectrons','systematicsMET','cleanPatJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
tagAndProbeConfigurator3.addSorter('dielectronsSortedByZMass','PATElePairSorterByZMass')
tagAndProbeConfigurator3.addSelector('selectedTagAndProbeDiElectrons','PATElePairSelector','calibratedMET.pt()<20&&leg1.pt()>20&&leg2.pt()>5&&charge==0&&(leg1.userIso(0)+max(0.0, leg1.userIso(1)+leg1.neutralHadronIso()-0.5*leg1.userIso(2)))/leg1.pt()<0.10&&leg1.userFloat("WWID")>0&&abs(leg1.eta())<2.5','EEss',0,100)
tagAndProbeSequenceDiElectrons =tagAndProbeConfigurator3.returnSequence()
# tagAndProbeConfigurator4 = CutSequenceProducer(initialCounter = 'initialTagAndProbeEvents4',
# pyModuleName = __name__,
# pyNameSpace = locals())
#
# ##Tag and probe electron+Tau
# tagAndProbeConfigurator4.addDiCandidateModule('tagAndProbeEleTaus','PATEleTauPairProducer', 'triggeredPatElectrons','triggeredPatTaus','systematicsMET','cleanPatJets',1,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
# tagAndProbeConfigurator4.addSelector('selectedTagAndProbeEleTaus','PATEleTauPairSelector','abs(leg2.eta())<2.1&&(!(abs(leg2.eta())>1.442&&abs(leg2.eta())<1.566))&&leg1.pt()>15&&charge==0&&leg2.pt()>5&&mass>15&&mass<120&&leg1.userFloat("WWID")>0&&(leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.1&&(leg2.userFloat("hltOverlapFilterIsoEle15IsoPFTau20")+leg2.userFloat("hltOverlapFilterIsoEle15IsoPFTau15")+leg2.userFloat("hltOverlapFilterIsoEle15TightIsoPFTau20")+leg2.userFloat("hltOverlapFilterIsoEle18MediumIsoPFTau20")+leg2.userFloat("hltOverlapFilterIsoEle20MediumIsoPFTau20"))>0.5','tagAndProbeEleTracks',0,100)
#
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
tagAndProbeConfigurator = CutSequenceProducer(initialCounter = 'initialTagAndProbeEvents',
pyModuleName = __name__,
pyNameSpace = locals())
##Tag and probe
tagAndProbeConfigurator.addDiCandidateModule('tagAndProbeMuTracks','PATMuTrackPairProducer', 'patMuonsForAnalysis','trackCandidates','patMETs','selectedPatJets',1,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
tagAndProbeConfigurator.addSelector('selectedTagAndProbeMuTracks','PATMuTrackPairSelector','calibratedMET.pt()<20 && charge==0&&mass>50&&mass<120&&leg1.userFloat("isWWMuon")&&((leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.2)&&leg1.pt()>25&&leg2.pt()>5&&abs(leg2.eta())<2.4&&abs(leg1.eta())<2.1','tagAndProbePairs',1,100)
tagAndProbeSequence =tagAndProbeConfigurator.returnSequence()
munuKIN='lepton.pt()>25 && abs(lepton.eta())<2.1'
munuSel='lepton.userFloat("isWWMuon")==1'#(lepton.isolationR03().sumPt+lepton.isolationR03().emEt+lepton.isolationR03().hadEt)/lepton.pt()<0.3&&met.pt()>30&mt>40&&lepton.userFloat("isWWMuon")==1'
elenuKIN='lepton.pt()>30 && abs(lepton.eta())<2.5&&(abs(lepton.eta)<1.4442||abs(lepton.eta)>1.5666)'
elenuSel='(lepton.isEB &&(lepton.sigmaIetaIeta<0.01)&&(-0.8<lepton.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<lepton.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (lepton.isEE &&(lepton.sigmaIetaIeta<0.03)&&(-0.7<lepton.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<lepton.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel='( (leg1.pt()>20 && leg2.pt()>10) || (leg1.pt()>10 && leg2.pt()>20))'
DiElePreSel2='(leg1.isEB &&(leg1.sigmaIetaIeta<0.01)&&(-0.8<leg1.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg1.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg1.isEE &&(leg1.sigmaIetaIeta<0.03)&&(-0.7<leg1.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg1.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel2='(leg2.isEB &&(leg2.sigmaIetaIeta<0.01)&&(-0.8<leg2.deltaPhiSuperClusterTrackAtVtx<0.8 ) && ( -0.007<leg2.deltaEtaSuperClusterTrackAtVtx<0.007) ) || (leg2.isEE &&(leg2.sigmaIetaIeta<0.03)&&(-0.7<leg2.deltaPhiSuperClusterTrackAtVtx<0.7 ) && ( -0.01<leg2.deltaEtaSuperClusterTrackAtVtx<0.01) )'
DiElePreSel3='(leg1.chargedHadronIso+leg1.photonIso()+leg1.neutralHadronIso())/leg1.pt()<0.3&&(leg2.chargedHadronIso+leg2.photonIso()+leg2.neutralHadronIso())/leg2.pt()<0.3'
######################__________________________________MNJ_____________________________________##############################
MNJanalysisConfigurator = CutSequenceProducer(initialCounter = 'initialEventsMNJ',
pyModuleName = __name__,
pyNameSpace = locals())
#Make Di Muons to VETO DY
#MNJanalysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'patMuonsForAnalysis','patMuonsForAnalysis','patMETs','patJetsForAnalysis',0,9999,text = '',dR=0.5,recoMode = "")
#MNJanalysisConfigurator.addSelector('diMuonsSel','PATMuPairSelector',DiMuonPreSel+'&&'+DiMuonPreSel2+'&&'+DiMuonPreSel3,'diMuonSel',0,999)
#MNJanalysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
#Make Muons+MET
MNJanalysisConfigurator.addCandidateMETModule('wCands','PATMuonNuPairProducer','patMuonsForAnalysis','patMETs','patJetsForAnalysis',1,9999,'AtLeastOneWCandidate',genParticles="genDaughters")
#MNJanalysisConfigurator.addSelector('wCandsTight','PATMuonNuPairSelector',munuTight,'tightID',1)
MNJanalysisConfigurator.addSelector('wCandsKIN','PATMuonNuPairSelector',munuKIN,'wCandsKIN',1)
MNJanalysisConfigurator.addSelector('wCandsSel','PATMuonNuPairSelector',munuSel,'wCandsSel',1)
import FWCore.ParameterSet.Config as cms
###MU ENRICHED QCD#####################################################################################################
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
wJetsConfigurator = CutSequenceProducer(initialCounter = 'initialWJetsEvents',
pyModuleName = __name__,
pyNameSpace = locals())
wJetsConfigurator.addDiCandidateModule('diMuonswJetsFR','PATMuPairProducer', 'patMuonsForAnalysis','patMuonsForAnalysis','patMETs','selectedPatJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
wJetsConfigurator.addSelector('osDiMuonswJetsFR','PATMuPairSelector','charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&((leg1.chargedHadronIso+leg1.photonIso+leg1.neutralHadronIso)/leg1.pt()<0.26)&&((leg2.chargedHadronIso+leg2.photonIso+leg2.neutralHadronIso)/leg2.pt()<0.3)','',0,0)
#Make Muons+MET
wJetsConfigurator.addCandidateMETModule('wJetsCands','PATMuonNuPairProducer','patMuonsForAnalysis','patMETs','selectedPatJets',0,9999,'')
wJetsConfigurator.addSelector('wJetsSelectedCands','PATMuonNuPairSelector','lepton.pt()>15&&abs(lepton.eta())<2.1&&((lepton.chargedHadronIso+lepton.photonIso+lepton.neutralHadronIso)/lepton.pt()<0.2)&&(lepton.triggerObjectMatchesByPath("HLT_Mu9").size>0||(lepton.triggerObjectMatchesByPath("HLT_Mu11").size>0)||(lepton.triggerObjectMatchesByPath("HLT_Mu11").size>0))&&mt>45.','',0,100)
wJetsConfigurator.addGeneric('wJetsForFakeRate','MuonNuPairJetFetcher')
wJetsSequence = wJetsConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfiguratorEMu = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
analysisConfiguratorEMu.addSmearing('patMuVetoTaus','patMuonsForAnalysis','convRejElectrons','patOverloadedJets','METInputTag')
#MakeDileptonPairs
analysisConfiguratorEMu.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuons','smearedMuons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfiguratorEMu.addSelector('osDiMuons','PATMuPairSelector','charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0),0.0))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.particleIso()))/leg2.pt()<0.2','DiMuonCreation',0,100)
analysisConfiguratorEMu.addSorter('diMuonsSorted','PATMuPairSorter')
analysisConfiguratorEMu.addDiCandidateModule('diElectrons','PATElePairProducer', 'smearedElectrons','smearedElectrons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfiguratorEMu.addSelector('osDiElectrons','PATElePairSelector','leg1.pt()>15&&leg2.pt()>15&&charge==0&&leg1.userFloat("WWID")>0&&leg2.userFloat("wp95")>0&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-leg1.userIso(0),0.0))/leg1.pt()<0.3&&(leg2.chargedHadronIso()+max(leg2.photonIso()+leg2.neutralHadronIso()-leg2.userIso(0),0.0))/leg2.pt()<0.3','ZEEVEto',0,100)
#Make muon.electron pairs
analysisConfiguratorEMu.addDiCandidateModule('eleMuons','PATEleMuPairProducer', 'smearedElectrons','smearedMuons','smearedMET','smearedJets',1,9999,text = 'AtLeastOneEleMu',leadingObjectsOnly = False,dR = 0.5,recoMode = "",genParticles='genDaughters')
analysisConfiguratorEMu.addSelector('eleMuonsElePtEta','PATEleMuPairSelector','leg1.pt()>10&&abs(leg1.eta())<2.5','electronPtEta',1)
analysisConfiguratorEMu.addSelector('eleMuonsMuPtEta','PATEleMuPairSelector','leg2.pt()>10&&abs(leg2.eta())<2.1','muonPtEta',1)
analysisConfiguratorEMu.addSelector('eleMuonsElectronWW','PATEleMuPairSelector','leg1.userFloat("WWID")>0||leg1.userFloat("WWMVAID")>0','WWElectron',1)
analysisConfiguratorEMu.addSelector('eleMuonsMuonWW','PATEleMuPairSelector','leg2.userFloat("isWWMuon")>0','MuonWW',1)
#analysisConfiguratorEMu.addEleMuSVFitSA('eleMuonsNSVFit')
analysisConfiguratorEMu.addSorter('eleMuonsSorted','PATEleMuPairSorter')
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
#Add smearing
analysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','patOverloadedJets','METInputTag')
#Create di muon pairs for veto purposes
analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'smearedMuons','smearedMuons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfigurator.addSelector('osDiMuons','PATMuPairSelector','charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.userIso(0)))/leg2.pt()<0.2','DiMuonCreation',0,0)
analysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
#Make DiTaus
analysisConfigurator.addDiCandidateModule('muJets','PATMuJetPairProducer','smearedMuons','smearedJets','smearedMET','smearedJets',1,9999,text='AtLeastOneDiTau',leadingObjectsOnly = False,dR = 0.5,recoMode ="",genParticles='genDaughters')
analysisConfigurator.addSelector('muJetsMuonID','PATMuJetPairSelector','leg1.userFloat("isWWMuon")','MuonWW',1)
analysisConfigurator.addSelector('muJetsMuonPtEta','PATMuJetPairSelector','leg1.pt()>15&&abs(leg1.eta())<2.1','MuonPtEta',1)
analysisConfigurator.addSelector('muJetsMuonIsolation','PATMuJetPairSelector','(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.15','MuonIsolation',1)
analysisConfigurator.addSelector('muJetsTauPtEta','PATMuJetPairSelector','leg2.pt()>20&&abs(leg2.eta())<2.3','TauPtEta',1)
#create the sequence
selectionSequence =analysisConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
#Make Di Muons to VETO DY
analysisConfigurator.addDiCandidateModule('diMuons','PATMuPairProducer', 'patMuonsForAnalysis','patMuonsForAnalysis','patMETs','selectedPatJets',0,9999,text = '',leadingObjectsOnly = True,dR = 0.15,recoMode = "")
analysisConfigurator.addSorter('diMuonsSorted','PATMuPairSorter')
#Make Muons+MET
analysisConfigurator.addCandidateMETModule('wCands','PATMuonNuPairProducer','patMuonsForAnalysis','patMETs','selectedPatJets',1,9999,'AtLeastOneWCandidate')
analysisConfigurator.addSelector('wCandsGoodMuon','PATMuonNuPairSelector','lepton.pt()>10&&abs(lepton.eta())<2.4','OneLeptonInAcceptance',1)
#analysisConfigurator.addSorter('diTausSorted','PATMuTauPairSorter')
#make Gen Ws
analysisConfigurator.addCandidateMETModule('wGenCands','PATCandNuPairProducer','genMuons','genNeutrinos','selectedPatJets',0,9999,'')
selectionSequence =analysisConfigurator.returnSequence()
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter='initialEvents',
pyModuleName=__name__,
pyNameSpace=locals())
#Add smearing
analysisConfigurator.addSmearing('patOverloadedTaus', 'patMuonsForAnalysis',
'convRejElectrons', 'patOverloadedJets',
'METInputTag')
#Create di muon pairs for veto purposes
analysisConfigurator.addDiCandidateModule('diMuons',
'PATMuPairProducer',
'smearedMuons',
'smearedMuons',
'smearedMET',
'smearedJets',
0,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
analysisConfigurator.addSelector(
'osDiMuons', 'PATMuPairSelector',
'charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(0.0,leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0)))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.userIso(0)))/leg2.pt()<0.2',
'DiMuonCreation', 0, 0)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfiguratorEMu = CutSequenceProducer(initialCounter='initialEvents',
pyModuleName=__name__,
pyNameSpace=locals())
analysisConfiguratorEMu.addSmearing('patMuVetoTaus', 'patMuonsForAnalysis',
'convRejElectrons', 'patOverloadedJets',
'METInputTag')
#MakeDileptonPairs
analysisConfiguratorEMu.addDiCandidateModule('diMuons',
'PATMuPairProducer',
'smearedMuons',
'smearedMuons',
'smearedMET',
'smearedJets',
0,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
analysisConfiguratorEMu.addSelector(
'osDiMuons', 'PATMuPairSelector',
'charge==0&&leg1.isGlobalMuon&&leg2.isGlobalMuon&&leg1.pt()>15&&leg2.pt()>15&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0),0.0))/leg1.pt()<0.2&&(leg2.chargedHadronIso()+max(0.0,leg2.photonIso+leg2.neutralHadronIso()-0.5*leg2.particleIso()))/leg2.pt()<0.2',
'DiMuonCreation', 0, 100)
analysisConfiguratorEMu.addSorter('diMuonsSorted', 'PATMuPairSorter')
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
analysisConfigurator = CutSequenceProducer(initialCounter = 'initialEvents',
pyModuleName = __name__,
pyNameSpace = locals())
nSVfitTrackService = cms.Service("NSVfitTrackService")
analysisConfigurator.addSmearing('patOverloadedTaus','patMuonsForAnalysis','convRejElectrons','patOverloadedJets')
#create dielectrons
analysisConfigurator.addDiCandidateModule('diElectrons','PATElePairProducer', 'smearedElectrons','smearedElectrons','smearedMET','smearedJets',0,9999,text = '',leadingObjectsOnly = False,dR = 0.15,recoMode = "")
analysisConfigurator.addSelector('osDiElectrons','PATElePairSelector','leg1.pt()>15&&leg2.pt()>15&&charge==0&&(leg1.userFloat("WWID")>0||leg1.userFloat("WWMVAID")>0)&&leg2.userFloat("wp95")>0&&(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.userIso(0),0.0))/leg1.pt()<0.25&&(leg2.chargedHadronIso()+max(leg2.photonIso()+leg2.neutralHadronIso()-0.5*leg2.userIso(0),0.0))/leg2.pt()<0.25','ZEEVEto',0,100)
#Make DiTaus
analysisConfigurator.addDiCandidateModule('eleTaus','PATEleTauPairProducer', 'smearedElectrons','smearedTaus','smearedMET','smearedJets',1,9999,text = 'AtLeastOneEleTau',leadingObjectsOnly = False,dR = 0.5,recoMode = "",genParticles='genDaughters')
analysisConfigurator.addSelector('eleTausWW','PATEleTauPairSelector','leg1.userFloat("WWID")>0||leg1.userFloat("WWMVAID")>0','WWElectronID',1)
analysisConfigurator.addSelector('eleTausElePtEta','PATEleTauPairSelector','leg1.pt()>20&&abs(leg1.eta())<2.1','electronPtEta',1)
#analysisConfigurator.addSelector('eleTausEleIsolation','PATEleTauPairSelector','(leg1.chargedHadronIso()+max(leg1.photonIso()+leg1.neutralHadronIso()-0.5*leg1.particleIso(),0.0))/leg1.pt()<0.1','electronIsolation',1)
analysisConfigurator.addSelector('eleTausTauPtEta','PATEleTauPairSelector','leg2.tauID("decayModeFinding")&&leg2.pt()>20&&abs(leg2.eta())<2.3','TauDecayFound',1)
analysisConfigurator.addSelector('eleTausVLooseIsolation','PATEleTauPairSelector','leg2.tauID("byVLooseCombinedIsolationDeltaBetaCorr")>0','TauLooseIso',1)
#analysisConfigurator.addSelector('eleTausLooseIsolation','PATEleTauPairSelector','leg2.tauID("byLooseCombinedIsolationDeltaBetaCorr")>0','TauLooseIso',1)
analysisConfigurator.addSelector('eleTausTauMuonVeto','PATEleTauPairSelector','leg2.tauID("againstMuonLoose")','AgainstMuon',1)
analysisConfigurator.addSelector('eleTausTauElectronVeto','PATEleTauPairSelector','leg2.tauID("againstElectronMedium")','AgainstElectron',1)
import FWCore.ParameterSet.Config as cms
#Import tool that creates the cut sequence
from UWAnalysis.Configuration.tools.CutSequenceProducer import *
zMuMuAnalysisConfigurator = CutSequenceProducer(
initialCounter='initialZMuMuEvents',
pyModuleName=__name__,
pyNameSpace=locals())
zMuMuAnalysisConfigurator.addDiCandidateModule('zMuMuCandidates',
'PATMuPairProducer',
'patMuonsForAnalysis',
'patMuonsForAnalysis',
'patMETs',
'patOverloadedJets',
0,
9999,
text='',
leadingObjectsOnly=False,
dR=0.15,
recoMode="")
zMuMuAnalysisConfigurator.addSorter('zMuMuCandidatesSorted', 'PATMuPairSorter')
zMuMuAnalysisConfigurator.addSelector(
'zMuMuCandidatesID', 'PATMuPairSelector',
'leg1.userFloat("isWWMuon")&&leg2.userFloat("isWWMuon")&&leg1.pt()>15&&leg2.pt()>15&&abs(leg1.eta)<2.4&&abs(leg2.eta)<2.4',
'ZMuMuID')
zMuMuAnalysisConfigurator.addSelector(
'zMuMuCandidatesIsol1', 'PATMuPairSelector',
'(leg1.chargedHadronIso()+leg1.photonIso()+leg1.neutralHadronIso())/leg1.pt()<0.25',
'ZMuMuLeg1Iso')
zMuMuAnalysisConfigurator.addSelector(
