optimized muon validation

Validation/RecoMuon/python/associators_cff.py

@@ -41,20 +41,26 @@
 seedsOfSTAmuons.L2seedsCollection = cms.InputTag("ancientMuonSeed")
 seedsOfSTAmuons_seq = cms.Sequence( seedsOfSTAmuons )
 
+# select probe tracks
+import PhysicsTools.RecoAlgos.recoTrackSelector_cfi
+probeTracks = PhysicsTools.RecoAlgos.recoTrackSelector_cfi.recoTrackSelector.clone()
+probeTracks.quality = cms.vstring('highPurity')
+probeTracks.tip = cms.double(3.5)
+probeTracks.lip = cms.double(30.)
+probeTracks.ptMin = cms.double(4.0)
+probeTracks.minRapidity = cms.double(-2.4)
+probeTracks.maxRapidity = cms.double(2.4)
+probeTracks_seq = cms.Sequence( probeTracks )
+
 #
 # Associators for Full Sim + Reco:
 #
 
 tpToTkmuTrackAssociation = cms.EDProducer('TrackAssociatorEDProducer',
     associator = cms.string('TrackAssociatorByHits'),
     label_tp = cms.InputTag('mix', 'MergedTrackTruth'),
-    label_tr = cms.InputTag('generalTracks')
-)
-
-tpToTkMuonTrackAssociation = cms.EDProducer('TrackAssociatorEDProducer',
-    associator = cms.string('MuonAssociatorByHits'),
-    label_tp = cms.InputTag('mix', 'MergedTrackTruth'),
-    label_tr = cms.InputTag('generalTracks','')
+#    label_tr = cms.InputTag('generalTracks')
+    label_tr = cms.InputTag('probeTracks')
 )
 
 tpToStaTrackAssociation = cms.EDProducer('TrackAssociatorEDProducer',
@@ -169,7 +175,8 @@
 tpToL3MuonAssociation = SimMuon.MCTruth.MuonAssociatorByHits_cfi.muonAssociatorByHits.clone()
 
 tpToTkMuonAssociation.tpTag = 'mix:MergedTrackTruth'
-tpToTkMuonAssociation.tracksTag = 'generalTracks'
+#tpToTkMuonAssociation.tracksTag = 'generalTracks'
+tpToTkMuonAssociation.tracksTag = 'probeTracks'
 tpToTkMuonAssociation.UseTracker = True
 tpToTkMuonAssociation.UseMuon = False
 
@@ -314,11 +321,10 @@
 
 muonAssociation_seq = cms.Sequence(
     extractedMuonTracks_seq
-    + seedsOfSTAmuons_seq
-    +(tpToTkMuonAssociation)
+    +seedsOfSTAmuons_seq
+    +probeTracks_seq+(tpToTkMuonAssociation+tpToTkmuTrackAssociation)
     +(tpToStaSeedAssociation+tpToStaMuonAssociation+tpToStaUpdMuonAssociation+tpToGlbMuonAssociation)
-    +(tpToTkmuTrackAssociation)
-#   +(tpToTkmuTrackAssociation+tpToStaTrackAssociation+tpToStaUpdTrackAssociation+tpToGlbTrackAssociation)
+#   +(tpToStaTrackAssociation+tpToStaUpdTrackAssociation+tpToGlbTrackAssociation)
 )
 muonAssociationTEV_seq = cms.Sequence(
     (tpToTevFirstMuonAssociation+tpToTevPickyMuonAssociation+tpToTevDytMuonAssociation)
@@ -348,7 +354,8 @@
 tpToTkmuTrackAssociationFS = cms.EDProducer('TrackAssociatorEDProducer',
     associator = cms.string('TrackAssociatorByHits'),
     label_tp = cms.InputTag('mix', 'MergedTrackTruth'),
-    label_tr = cms.InputTag('generalTracks')
+#    label_tr = cms.InputTag('generalTracks')
+    label_tr = cms.InputTag('probeTracks')
 )
 
 tpToStaTrackAssociationFS = cms.EDProducer('TrackAssociatorEDProducer',
@@ -441,7 +448,8 @@
 baseMuonAssociatorFS.ROUList = ['famosSimHitsTrackerHits']
 
 
-tpToTkMuonAssociationFS   = baseMuonAssociatorFS.clone()
+tpToTkMuonAssociationFS  = baseMuonAssociatorFS.clone()
+tpToStaSeedAssociationFS = baseMuonAssociatorFS.clone()
 tpToStaMuonAssociationFS  = baseMuonAssociatorFS.clone()
 tpToStaUpdMuonAssociationFS  = baseMuonAssociatorFS.clone()
 tpToStaRefitMuonAssociationFS  = baseMuonAssociatorFS.clone()
@@ -455,11 +463,18 @@
 tpToL2UpdMuonAssociationFS   = baseMuonAssociatorFS.clone()
 tpToL3MuonAssociationFS   = baseMuonAssociatorFS.clone()
 
-tpToTkMuonAssociationFS.tracksTag = 'generalTracks'
+#tpToTkMuonAssociationFS.tracksTag = 'generalTracks'
+tpToTkMuonAssociationFS.tracksTag = 'probeTracks'
 tpToTkMuonAssociationFS.tpTag = 'mix:MergedTrackTruth'
 tpToTkMuonAssociationFS.UseTracker = True
 tpToTkMuonAssociationFS.UseMuon = False
 
+tpToStaSeedAssociationFS.tpTag = 'mix:MergedTrackTruth'
+tpToStaSeedAssociationFS.tracksTag = 'seedsOfSTAmuons'
+tpToStaSeedAssociationFS.UseTracker = False
+tpToStaSeedAssociationFS.UseMuon = True
+
+tpToStaMuonAssociationFS.tpTag = 'mix:MergedTrackTruth'
 tpToStaMuonAssociationFS.tracksTag = 'standAloneMuons'
 tpToStaMuonAssociationFS.UseTracker = False
 tpToStaMuonAssociationFS.UseMuon = True
@@ -526,10 +541,11 @@
 
 muonAssociationFastSim_seq = cms.Sequence(
         extractedMuonTracks_seq
-        +(tpToTkMuonAssociationFS+tpToStaMuonAssociationFS+tpToStaUpdMuonAssociationFS+tpToGlbMuonAssociationFS)
+        +seedsOfSTAmuons_seq
+        +probeTracks+(tpToTkMuonAssociationFS+tpToTkmuTrackAssociationFS) 
+        +(tpToStaSeedAssociationFS+tpToStaMuonAssociationFS+tpToStaUpdMuonAssociationFS+tpToGlbMuonAssociationFS)
         +(tpToStaRefitMuonAssociationFS+tpToStaRefitUpdMuonAssociationFS)
         +(tpToTevFirstMuonAssociationFS+tpToTevPickyMuonAssociationFS+tpToTevDytMuonAssociationFS)
-        +tpToTkmuTrackAssociationFS
 #        +tpToStaTrackAssociationFS+tpToStaUpdTrackAssociationFS+tpToGlbTrackAssociationFS
 #        +tpToTevFirstTrackAssociationFS+tpToTevPickyTrackAssociationFS
         )

Validation/RecoMuon/python/muonValidationFastSim_cff.py

@@ -51,6 +51,18 @@
 tevMuonDytTrackVTrackAssocFS.label_tp_effic = 'mix:MergedTrackTruth'
 tevMuonDytTrackVTrackAssocFS.label_tp_fake = 'mix:MergedTrackTruth'
 
+from Validation.RecoMuon.muonValidation_cff import trkProbeTrackVMuonAssoc 
+trkProbeTrackVMuonAssocFS = trkProbeTrackVMuonAssoc.clone()
+trkProbeTrackVMuonAssocFS.associatormap = 'tpToTkMuonAssociationFS'
+trkProbeTrackVMuonAssocFS.label_tp_effic = 'mix:MergedTrackTruth'
+trkProbeTrackVMuonAssocFS.label_tp_fake = 'mix:MergedTrackTruth'
+
+from Validation.RecoMuon.muonValidation_cff import staSeedTrackVMuonAssoc
+staSeedTrackVMuonAssocFS = staSeedTrackVMuonAssoc.clone() 
+staSeedTrackVMuonAssocFS.associatormap = 'tpToStaSeedAssociationFS' 
+staSeedTrackVMuonAssocFS.label_tp_effic = 'mix:MergedTrackTruth'
+staSeedTrackVMuonAssocFS.label_tp_fake = 'mix:MergedTrackTruth'
+
 from Validation.RecoMuon.muonValidation_cff import staMuonTrackVMuonAssoc
 
 staMuonTrackVMuonAssocFS = staMuonTrackVMuonAssoc.clone()
@@ -112,7 +124,8 @@
 
 
 # Muon validation sequence
-muonValidationFastSim_seq = cms.Sequence(trkMuonTrackVTrackAssocFS
+muonValidationFastSim_seq = cms.Sequence(trkProbeTrackVMuonAssocFS+trkMuonTrackVTrackAssocFS
+                                         +staSeedTrackVMuonAssocFS
                                          +staMuonTrackVMuonAssocFS+staUpdMuonTrackVMuonAssocFS+glbMuonTrackVMuonAssocFS
                                          +staRefitMuonTrackVMuonAssocFS+staRefitUpdMuonTrackVMuonAssocFS
                                          +tevMuonFirstTrackVMuonAssocFS+tevMuonPickyTrackVMuonAssocFS+tevMuonDytTrackVMuonAssocFS

Validation/RecoMuon/python/muonValidation_cff.py

@@ -8,18 +8,11 @@
 trkMuonTrackVTrackAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
 trkMuonTrackVTrackAssoc.associatormap = 'tpToTkmuTrackAssociation'
 trkMuonTrackVTrackAssoc.associators = ('TrackAssociatorByHits',)
-trkMuonTrackVTrackAssoc.label = ('generalTracks',)
+#trkMuonTrackVTrackAssoc.label = ('generalTracks',)
+trkMuonTrackVTrackAssoc.label = ('probeTracks',)
 trkMuonTrackVTrackAssoc.usetracker = True
 trkMuonTrackVTrackAssoc.usemuon = False
 
-
-trkMuonTrackVMuonTrackAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
-trkMuonTrackVMuonTrackAssoc.associatormap = 'tpToTkMuonAssociation'
-trkMuonTrackVMuonTrackAssoc.associators = ('MuonAssociatorByHits',)
-trkMuonTrackVMuonTrackAssoc.label = ('generalTracks',)
-trkMuonTrackVMuonTrackAssoc.usetracker = True
-trkMuonTrackVMuonTrackAssoc.usemuon = False
-
 trkCosmicMuonTrackVTrackAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
 trkCosmicMuonTrackVTrackAssoc.associatormap = 'tpToTkCosmicTrackAssociation'
 trkCosmicMuonTrackVTrackAssoc.associators = ('TrackAssociatorByHits',)
@@ -104,12 +97,14 @@
 glbCosmicMuonTrackVTrackAssoc.usetracker = True
 glbCosmicMuonTrackVTrackAssoc.usemuon = True
 
-staMuonTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
-staMuonTrackVMuonAssoc.associatormap = 'tpToStaMuonAssociation'
-staMuonTrackVMuonAssoc.associators = ('MuonAssociationByHits',)
-staMuonTrackVMuonAssoc.label = ('standAloneMuons',)
-staMuonTrackVMuonAssoc.usetracker = False
-staMuonTrackVMuonAssoc.usemuon = True
+trkProbeTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
+#trkMuonTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
+trkProbeTrackVMuonAssoc.associatormap = 'tpToTkMuonAssociation' 
+trkProbeTrackVMuonAssoc.associators = ('MuonAssociationByHits',)
+##trkMuonTrackVMuonAssoc.label = ('generalTracks',)
+trkProbeTrackVMuonAssoc.label = ('probeTracks',)
+trkProbeTrackVMuonAssoc.usetracker = True
+trkProbeTrackVMuonAssoc.usemuon = False
 
 staSeedTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
 staSeedTrackVMuonAssoc.associatormap = 'tpToStaSeedAssociation'
@@ -118,6 +113,12 @@
 staSeedTrackVMuonAssoc.usetracker = False
 staSeedTrackVMuonAssoc.usemuon = True
 
+staMuonTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
+staMuonTrackVMuonAssoc.associatormap = 'tpToStaMuonAssociation'
+staMuonTrackVMuonAssoc.associators = ('MuonAssociationByHits',)
+staMuonTrackVMuonAssoc.label = ('standAloneMuons',)
+staMuonTrackVMuonAssoc.usetracker = False
+staMuonTrackVMuonAssoc.usemuon = True
 
 staUpdMuonTrackVMuonAssoc = Validation.RecoMuon.MuonTrackValidator_cfi.muonTrackValidator.clone()
 staUpdMuonTrackVMuonAssoc.associatormap = 'tpToStaUpdMuonAssociation'
@@ -315,13 +316,11 @@
 recoMuonVMuAssoc_tgt.primaryVertex = 'offlinePrimaryVertices'
 
 # Muon validation sequence
-muonValidation_seq = cms.Sequence(trkMuonTrackVTrackAssoc
-                                 +staMuonTrackVMuonAssoc+staUpdMuonTrackVMuonAssoc+glbMuonTrackVMuonAssoc
-                                  +recoMuonVMuAssoc_trk+recoMuonVMuAssoc_sta+recoMuonVMuAssoc_glb+recoMuonVMuAssoc_tgt)
 
-muonValidationTest_seq = cms.Sequence(trkMuonTrackVMuonTrackAssoc
-                                 +staMuonTrackVMuonAssoc+staSeedTrackVMuonAssoc+staUpdMuonTrackVMuonAssoc+glbMuonTrackVMuonAssoc
-                                  +recoMuonVMuAssoc_trk+recoMuonVMuAssoc_sta+recoMuonVMuAssoc_glb+recoMuonVMuAssoc_tgt)
+muonValidation_seq = cms.Sequence(trkProbeTrackVMuonAssoc+trkMuonTrackVTrackAssoc
+                                 +staSeedTrackVMuonAssoc
+                                 +staMuonTrackVMuonAssoc+staUpdMuonTrackVMuonAssoc+glbMuonTrackVMuonAssoc
+                                 +recoMuonVMuAssoc_trk+recoMuonVMuAssoc_sta+recoMuonVMuAssoc_glb+recoMuonVMuAssoc_tgt)
 
 muonValidationTEV_seq = cms.Sequence(tevMuonFirstTrackVMuonAssoc+tevMuonPickyTrackVMuonAssoc+tevMuonDytTrackVMuonAssoc)
 
@@ -334,8 +333,7 @@
 
 # The muon association and validation sequence
 
-recoMuonValidation = cms.Sequence((muonAssociation_seq*muonValidationTest_seq)
-#                                 +(muonAssociation_seq*muonValidationTest_seq)
+recoMuonValidation = cms.Sequence((muonAssociation_seq*muonValidation_seq)
                                  +(muonAssociationTEV_seq*muonValidationTEV_seq)
                                  +(muonAssociationSET_seq*muonValidationSET_seq)
                                  +(muonAssociationRefit_seq*muonValidationRefit_seq)