Skims/python/muonTagProbeFilters_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #from HLTrigger.HLTfilters.hltHighLevel_cfi import *

0004 #exoticaMuHLT = hltHighLevel

0005 #Define the HLT path to be used.

0006 #exoticaMuHLT.HLTPaths =['HLT_L1MuOpen']

0007 #exoticaMuHLT.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29")

0008
0009 #Define the HLT quality cut

0010 #exoticaHLTMuonFilter = cms.EDFilter("HLTSummaryFilter",

0011 #    summary = cms.InputTag("hltTriggerSummaryAOD","","HLT8E29"), # trigger summary

0012 #    member  = cms.InputTag("hltL3MuonCandidates","","HLT8E29"),      # filter or collection

0013 #    cut     = cms.string("pt>0"),                     # cut on trigger object

0014 #    minN    = cms.int32(0)                  # min. # of passing objects needed

0015 # )

0016
0017
0018 #Define the Reco quality cut

0019 from SimGeneral.HepPDTESSource.pythiapdt_cfi import *
0020
0021 # Make the charged candidate collections from tracks

0022 allTracks = cms.EDProducer("TrackViewCandidateProducer",
0023                                src = cms.InputTag("generalTracks"),
0024                                particleType = cms.string('mu+'),
0025                                cut = cms.string('pt > 0'),
0026                            filter = cms.bool(True)
0027                            )
0028
0029 staTracks = cms.EDProducer("TrackViewCandidateProducer",
0030                                src = cms.InputTag("standAloneMuons","UpdatedAtVtx"),
0031                                particleType = cms.string('mu+'),
0032                                cut = cms.string('pt > 0.5 && abs(d0) < 2.0 && abs(vz) < 25.0'),
0033                            filter = cms.bool(True)
0034                            )
0035
0036 # Make the input candidate collections

0037 tagCands = cms.EDFilter("MuonRefSelector",
0038                             src = cms.InputTag("muons"),
0039                             cut = cms.string('isGlobalMuon > 0 && pt > 1.0 && abs(eta) < 2.1'),
0040                         filter = cms.bool(True)
0041                         )
0042
0043 # Standalone muon tracks (probes)

0044 staCands = cms.EDFilter("RecoChargedCandidateRefSelector",
0045                             src = cms.InputTag("staTracks"),
0046                             cut = cms.string('pt > 0.5 && abs(eta) < 2.1'),
0047                         filter = cms.bool(True)
0048                         )
0049
0050 # Tracker muons (to be matched)

0051 tkProbeCands = cms.EDFilter("RecoChargedCandidateRefSelector",
0052                                 src = cms.InputTag("allTracks"),
0053                                 cut = cms.string('pt > 0.5'),
0054                             filter = cms.bool(True)
0055                             )
0056
0057 # Match track and stand alone candidates

0058 # to get the passing probe candidates

0059 TkStaMap = cms.EDProducer("TrivialDeltaRViewMatcher",
0060                             src = cms.InputTag("tkProbeCands"),
0061                             distMin = cms.double(0.15),
0062                             matched = cms.InputTag("staCands"),
0063                         filter = cms.bool(True)
0064                         )
0065
0066 # Use the producer to get a list of matched candidates

0067 TkStaMatched = cms.EDProducer("RecoChargedCandidateMatchedProbeMaker",
0068                                 Matched = cms.untracked.bool(True),
0069                                 ReferenceSource = cms.untracked.InputTag("staCands"),
0070                                 ResMatchMapSource = cms.untracked.InputTag("TkStaMap"),
0071                                 CandidateSource = cms.untracked.InputTag("tkProbeCands"),
0072                             filter = cms.bool(True)
0073                             )
0074
0075 TkStaUnmatched = cms.EDProducer("RecoChargedCandidateMatchedProbeMaker",
0076                                   Matched = cms.untracked.bool(False),
0077                                   ReferenceSource = cms.untracked.InputTag("staCands"),
0078                                   ResMatchMapSource = cms.untracked.InputTag("TkStaMap"),
0079                                   CandidateSource = cms.untracked.InputTag("tkProbeCands"),
0080                               filter = cms.bool(True)
0081                               )
0082
0083 # Make the tag probe association map

0084 JPsiMMTagProbeMap = cms.EDProducer("TagProbeMassProducer",
0085                                      MassMaxCut = cms.untracked.double(4.5),
0086                                      TagCollection = cms.InputTag("tagCands"),
0087                                      MassMinCut = cms.untracked.double(1.5),
0088                                      ProbeCollection = cms.InputTag("tkProbeCands"),
0089                                      PassingProbeCollection = cms.InputTag("TkStaMatched")
0090                                  )
0091
0092 JPsiMMTPFilter = cms.EDFilter("TagProbeMassEDMFilter",
0093                         tpMapName = cms.string('JPsiMMTagProbeMap')
0094                         )
0095
0096 ZMMTagProbeMap = cms.EDProducer("TagProbeMassProducer",
0097                                  MassMaxCut = cms.untracked.double(120.0),
0098                                  TagCollection = cms.InputTag("tagCands"),
0099                                  MassMinCut = cms.untracked.double(50.0),
0100                                  ProbeCollection = cms.InputTag("tkProbeCands"),
0101                                  PassingProbeCollection = cms.InputTag("TkStaMatched")
0102                                  )
0103
0104 ZMMTPFilter = cms.EDFilter("TagProbeMassEDMFilter",
0105                         tpMapName = cms.string('ZMMTagProbeMap')
0106                         )
0107
0108
0109 #Define group sequence, using HLT/Reco quality cut.

0110 #exoticaMuHLTQualitySeq = cms.Sequence()

0111 tagProbeSeq = cms.Sequence(allTracks+staTracks*tagCands+tkProbeCands+staCands*TkStaMap*TkStaMatched)
0112
0113 muonJPsiMMRecoQualitySeq = cms.Sequence(
0114     #exoticaMuHLT+

0115     tagProbeSeq+JPsiMMTagProbeMap+JPsiMMTPFilter
0116     )
0117
0118 muonZMMRecoQualitySeq = cms.Sequence(
0119     #exoticaMuHLT+

0120     tagProbeSeq+ZMMTagProbeMap+ZMMTPFilter
0121     )
0122