Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:09:15

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