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
singlePhotonPt20Filter = cms.EDFilter("PhotonSelector",
src = cms.InputTag("photons"),
cut = cms.string('pt > 20 && abs(eta) < 2.7' ),
filter = cms.bool(True)
)
singlePhotonPt15Filter = cms.EDFilter("PhotonSelector",
src = cms.InputTag("photons"),
cut = cms.string('pt > 15 && abs(eta) < 2.7' ),
filter = cms.bool(True)
)
singlePhotonPt10Filter = cms.EDFilter("PhotonSelector",
src = cms.InputTag("photons"),
cut = cms.string('pt > 10 && abs(eta) < 2.7' ),
filter = cms.bool(True)
)
singlePhotonPt5Filter = cms.EDFilter("PhotonSelector",
src = cms.InputTag("photons"),
cut = cms.string('pt > 5 && abs(eta) < 2.7' ),
filter = cms.bool(True)
)
egSuperClusterMerger = cms.EDProducer("EgammaSuperClusterMerger",
src = cms.VInputTag(cms.InputTag('correctedHybridSuperClusters'),
cms.InputTag('correctedMulti5x5SuperClustersWithPreshower'))
)
egSuperClusterCands = cms.EDProducer("ConcreteEcalCandidateProducer",
src = cms.InputTag("egSuperClusterMerger"),
particleType = cms.string('gamma')
)
eggoodSuperClusters = cms.EDFilter("CandViewRefSelector",
src = cms.InputTag("egSuperClusterCands"),
cut = cms.string('et > 3.0')
)
egSuperClusterPt5Filter = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag("eggoodSuperClusters"),
minNumber = cms.uint32(2)
)
oneEmCluster = cms.Sequence(
egSuperClusterMerger+egSuperClusterCands+eggoodSuperClusters+egSuperClusterPt5Filter
)
gammaJet = cms.EDFilter("EtaPtMinCandViewSelector",
src = cms.InputTag("iterativeCone5CaloJets"),
ptMin = cms.double(5),
etaMin = cms.double(-2),
etaMax = cms.double(2)
)
gammaJetFilter = cms.EDFilter("CandViewCountFilter",
src = cms.InputTag("gammaJet"),
minNumber = cms.uint32(2)
)
#Define group sequence, using HLT/Reco quality cut.
#exoticaMuHLTQualitySeq = cms.Sequence()
singlePhotonPt20QualitySeq = cms.Sequence(
#exoticaMuHLT+
singlePhotonPt20Filter
)
singlePhotonPt15QualitySeq = cms.Sequence(
#exoticaMuHLT+
singlePhotonPt15Filter
)
singlePhotonPt10QualitySeq = cms.Sequence(
#exoticaMuHLT+
singlePhotonPt10Filter
)
singlePhotonPt5QualitySeq = cms.Sequence(
#exoticaMuHLT+
oneEmCluster+gammaJet+gammaJetFilter+singlePhotonPt5Filter
)
|