|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
|
import FWCore.ParameterSet.Config as cms
'''
## 2012 HLT COND.S AS GIVEN BY ARUN IN 53X
from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel as hlt_selector
hlt_selector.throw = cms.bool(False)
TauSkimMuTauMETHLT = hlt_selector.clone()
TauSkimMuTauMETHLT.TriggerResultsTag = cms.InputTag('TriggerResults', '', 'HLT')
TauSkimMuTauMETHLT.andOr=cms.bool(True)
TauSkimMuTauMETHLT.HLTPaths = cms.vstring("HLT_IsoMu15_eta2p1_L1ETM20_v*", "HLT_IsoMu15Rho_eta2p1_L1ETM20_v*")
'''
TauSkimPFTausSelected = cms.EDFilter("PFTauSelector",
src = cms.InputTag("hpsPFTauProducer"),
discriminators = cms.VPSet(
cms.PSet( #discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"), #53X AND 75X
discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFindingNewDMs"), #HTT 2015 TWIKI
selectionCut=cms.double(0.5)
) ## COMMENT OUT IF USING THE PSET BELOW
# ),
#cms.PSet(#discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr"), ## 53X
# #discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseIsolation"), #75X
# discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr3Hits"), #HTT 2015 TWIKI
# selectionCut=cms.double(0.5)
# ),
),
discriminatorContainers = cms.VPSet(),
#cut = cms.string('pt > 22. && abs(eta) < 2.3') #53X
#cut = cms.string('et > 15. && abs(eta) < 2.5') #75X
#cut = cms.string('pt > 18. && abs(eta) < 2.3') #HTT 2015 TWIKI
cut = cms.string('pt > 18. && abs(eta) < 2.3 && hpsPFTauMVA3IsolationChargedIsoPtSum < 2') #MICHAL'S SUGGESTION
)
TauSkimPFTauSkimmedBy1 = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag('TauSkimPFTausSelected'),
minNumber = cms.uint32(1)
)
TauSkimPFTauSkimmedBy2 = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag('TauSkimPFTausSelected'),
minNumber = cms.uint32(2)
)
## MODULE IN 53X ONLY
TauSkimDiTauPairs = cms.EDProducer("CandViewShallowCloneCombiner",
decay = cms.string("TauSkimPFTausSelected TauSkimPFTausSelected"),
checkCharge = cms.bool(False),
cut = cms.string("sqrt((daughter(0).eta-daughter(1).eta)*(daughter(0).eta-daughter(1).eta)+ min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi) ) * min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi) ) )>0.5"),
)
## MODULE IN 53X ONLY
TauSkimDiTauPairFilter = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag("TauSkimDiTauPairs"),
minNumber = cms.uint32(1)
)
## MODULE IN 53X ONLY
TauSkimPFTausSelectedForMuTau = TauSkimPFTausSelected.clone()
TauSkimPFTausSelectedForMuTau.discriminators = cms.VPSet(
cms.PSet( #discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"), # 53X AND 75X
discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFindingNewDMs"), #HTT 2015 TWIKI
selectionCut=cms.double(0.5)
),
# cms.PSet( #discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr"), #53X
# #discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseIsolation"), #75X
# discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr3Hits"), #HTT 2015 TWIKI
# selectionCut=cms.double(0.5)
# ),
#cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByElectronVLooseMVA5"), #HTT 2015 TWIKI (not working!)
# selectionCut=cms.double(0.5)
# ),
)
TauSkimPFTausSelectedForMuTau.discriminatorContainers = cms.VPSet(
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByMuonRejection3"),
workingPoints=cms.vstring("ByLooseMuonRejection3"),
rawValues=cms.vstring(),
selectionCuts=cms.vdouble()
)
)
## MODULE IN 53X ONLY
#TauSkimPFTausSelectedForMuTau.cut = cms.string('pt > 18. && abs(eta) < 2.3') #75X
TauSkimPFTausSelectedForMuTau.cut = cms.string('pt > 18. && abs(eta) < 2.3') #HTT 2015 TWIKI
TauSkimPFTauSkimmedForMuTau = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag('TauSkimPFTausSelectedForMuTau'),
minNumber = cms.uint32(1)
)
## NO MuTauMET MODULES IMPLEMETED HERE FROM 53X
## MODULE IN 53X ONLY
TauSkimMuonSelected = cms.EDFilter("MuonRefSelector",
src = cms.InputTag("muons"),
cut = cms.string("pt > 20 && " + "abs(eta) < 2.1 && " +
"isGlobalMuon && isTrackerMuon" +
#" && globalTrack.isNonnull "+
#" && globalTrack.hitPattern.numberOfValidTrackerHits>=5"+
#" && globalTrack.normalizedChi2<20"+
" && (pfIsolationR03.sumChargedHadronPt/pt) < 0.3"
),
)
## MODULE IN 53X ONLY
TauSkimMuonSkimmedBy1 = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag('TauSkimMuonSelected'),
minNumber = cms.uint32(1)
)
## MODULE IN 53X ONLY
TauSkimMuTauPairs = cms.EDProducer("CandViewShallowCloneCombiner",
decay = cms.string("TauSkimMuonSelected TauSkimPFTausSelectedForMuTau"),
checkCharge = cms.bool(False),
cut = cms.string("sqrt((daughter(0).eta-daughter(1).eta)*(daughter(0).eta-daughter(1).eta)+ min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi) ) * min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi) ) )>0.3"), ## DR CUT LOOSENED FOLLOWING MICHAL'S SUGGESTIONS FROM 0.5 TO 0.3
)
## NO MuTauMET MODULES IMPLEMETED HERE FROM 53X
## MODULE IN 53X ONLY
TauSkimMuTauPairFilter = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag("TauSkimMuTauPairs"),
minNumber = cms.uint32(1)
)
## NO MuTauMET AND DITAU SEQUENCES IMPLEMETED HERE FROM 53X
## MUTAU SEQUENCE IN 53X ONLY
mutauSkimSequence = cms.Sequence(
TauSkimMuonSelected *
TauSkimMuonSkimmedBy1 *
TauSkimPFTausSelectedForMuTau *
TauSkimPFTauSkimmedForMuTau *
TauSkimMuTauPairs *
TauSkimMuTauPairFilter
)
|