|
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
|
import FWCore.ParameterSet.Config as cms
#from HLTrigger.HLTfilters.hltHighLevel_cfi import *
#exoticaMuHLT = hltHighLevel
#Define the HLT path to be used.
#exoticaMuHLT.HLTPaths =['HLT_L1MuOpen']
#exoticaMuHLT.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29")
#Define the HLT quality cut
#exoticaHLTMuonFilter = cms.EDFilter("HLTSummaryFilter",
# summary = cms.InputTag("hltTriggerSummaryAOD","","HLT8E29"), # trigger summary
# member = cms.InputTag("hltL3MuonCandidates","","HLT8E29"), # filter or collection
# cut = cms.string("pt>0"), # cut on trigger object
# minN = cms.int32(0) # min. # of passing objects needed
# )
#Define the Reco quality cut
singleRecoMuonPt20Filter = cms.EDFilter("MuonRefSelector",
src = cms.InputTag("muons"),
cut = cms.string('pt > 20.0 && abs(eta) < 2.1' ),
filter = cms.bool(True),
minN = cms.int32(1)
)
singleRecoMuonPt15Filter = cms.EDFilter("MuonRefSelector",
src = cms.InputTag("muons"),
cut = cms.string('pt > 15.0 && abs(eta) < 2.1' ),
filter = cms.bool(True),
minN = cms.int32(1)
)
singleRecoMuonPt10Filter = cms.EDFilter("MuonRefSelector",
src = cms.InputTag("muons"),
cut = cms.string('pt > 10.0 && abs(eta) < 2.1' ),
filter = cms.bool(True),
minN = cms.int32(1)
)
singleRecoMuonPt5Filter = cms.EDFilter("MuonRefSelector",
src = cms.InputTag("muons"),
cut = cms.string('pt > 0.5 && abs(eta) < 2.1 && isGlobalMuon == 1 && abs (track().d0) < 2 ' ),
filter = cms.bool(True),
minN = cms.int32(1)
)
muonSCMerger = cms.EDFilter("EgammaSuperClusterMerger",
src = cms.VInputTag(cms.InputTag('correctedHybridSuperClusters'),
cms.InputTag('correctedMulti5x5SuperClustersWithPreshower'))
)
muonSCCands = cms.EDProducer("ConcreteEcalCandidateProducer",
src = cms.InputTag("muonSCMerger"),
particleType = cms.string('mu+')
)
goodMuonSCs = cms.EDFilter("CandViewRefSelector",
src = cms.InputTag("muonSCCands"),
cut = cms.string('et > 2.0')
)
muonSCPt5Filter = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag("goodMuonSCs"),
minNumber = cms.uint32(2)
)
muonEmClusters = cms.Sequence(
muonSCMerger+muonSCCands+goodMuonSCs+muonSCPt5Filter
)
#Define group sequence, using HLT/Reco quality cut.
#singleMuHLTQualitySeq = cms.Sequence()
singleMuPt20RecoQualitySeq = cms.Sequence(
#singleMuHLT+
singleRecoMuonPt20Filter
)
singleMuPt15RecoQualitySeq = cms.Sequence(
#singleMuHLT+
singleRecoMuonPt15Filter
)
singleMuPt10RecoQualitySeq = cms.Sequence(
#singleMuHLT+
singleRecoMuonPt10Filter
)
singleMuPt5RecoQualitySeq = cms.Sequence(
#singleMuHLT+
# muonEmClusters+
singleRecoMuonPt5Filter
)
|