CDC in 80X

Configuration/StandardSequences/python/ReconstructionCosmics_cff.py

@@ -69,6 +69,7 @@
                                              jetsCosmics*
                                              muonsCosmics*
                                              regionalCosmicTracksSeq*
+                                             cosmicDCTracksSeq*
                                              metrecoCosmics*
                                              egammaCosmics*
                                              logErrorHarvester)
@@ -77,12 +78,14 @@
                                              jetsCosmics*
                                              muonsCosmics*
                                              regionalCosmicTracksSeq*
+                                             cosmicDCTracksSeq*
                                              metrecoCosmics*
                                              egammaCosmics*
                                              logErrorHarvester)
 reconstructionCosmics_woTkBHM = cms.Sequence(localReconstructionCosmics*
                                              jetsCosmics*
                                              muonsCosmics*
                                              regionalCosmicTracksSeq*
+                                             cosmicDCTracksSeq*
                                              metrecoCosmics*
                                              egammaCosmics)

Configuration/StandardSequences/python/Reconstruction_cff.py

@@ -43,6 +43,9 @@
 
 from RecoLocalCalo.CastorReco.CastorSimpleReconstructor_cfi import *
 
+# Cosmic During Collisions
+from RecoTracker.SpecialSeedGenerators.cosmicDC_cff import *
+
 localreco = cms.Sequence(trackerlocalreco+muonlocalreco+calolocalreco+castorreco)
 localreco_HcalNZS = cms.Sequence(trackerlocalreco+muonlocalreco+calolocalrecoNZS+castorreco)
 
@@ -91,7 +94,8 @@
                              btagging*
                              recoPFMET*
                              PFTau*
-                             reducedRecHits
+                             reducedRecHits*
+                             cosmicDCTracksSeq
                              )
 
 

FastSimulation/Configuration/python/Reconstruction_AftMix_cff.py

@@ -135,6 +135,7 @@
 ##########################################
 # not commisoned and not relevant in FastSim (?):
 _reco.globalreco.remove(_reco.muoncosmicreco)
+_reco.highlevelreco.remove(_reco.cosmicDCTracksSeq)
 _reco.highlevelreco.remove(_reco.muoncosmichighlevelreco)
 _reco.muons.FillCosmicsIdMap = False
 

RecoMuon/Configuration/python/RecoMuonCosmics_cff.py

@@ -29,6 +29,17 @@
 muons.TrackExtractorPSet.inputTrackCollection = 'ctfWithMaterialTracksP5LHCNavigation'
 muons.CaloExtractorPSet.CenterConeOnCalIntersection = True
 
+#similar to what's in pp configuration
+muonsFromCosmics = muons1stStep.clone()
+muonsFromCosmics.inputCollectionLabels = ['cosmicMuons']
+muonsFromCosmics.inputCollectionTypes = ['outer tracks']
+muonsFromCosmics.TrackExtractorPSet.inputTrackCollection = 'cosmicMuons'
+muonsFromCosmics.TimingFillerParameters.DTTimingParameters.PruneCut = 9999
+muonsFromCosmics.TimingFillerParameters.CSCTimingParameters.PruneCut = 9999
+muonsFromCosmics.fillIsolation = False
+muonsFromCosmics.fillGlobalTrackQuality = False
+muonsFromCosmics.fillGlobalTrackRefits = False
+
 from RecoMuon.MuonIdentification.calomuons_cfi import *
 calomuons.inputTracks = 'ctfWithMaterialTracksP5LHCNavigation'
 calomuons.inputCollection = 'muons'
@@ -67,7 +78,7 @@
 allmuons = cms.Sequence(glbTrackQual*tevMuons*muons*muIsolation*calomuons)
 
 # Final sequence
-muonrecoforcosmics = cms.Sequence(muontrackingforcosmics*allmuons)
+muonrecoforcosmics = cms.Sequence(muontrackingforcosmics*allmuons*muonsFromCosmics)
 muonRecoAllGR = cms.Sequence(muonrecoforcosmics)
 
 # 1 leg mode

RecoTracker/Configuration/python/RecoTrackerP5_EventContent_cff.py

@@ -58,6 +58,9 @@
         'keep *_dedxTruncated40CosmicTF_*_*',
         'keep *_dedxHitInfoCosmicTF_*_*',
         'keep *_dedxHarmonic2CosmicTF_*_*',
+        'keep recoTracks_cosmicDCTracks_*_*',
+        'keep recoTrackExtras_cosmicDCTracks_*_*',
+        'keep TrackingRecHitsOwned_cosmicDCTracks_*_*',
     )
 )
 #RECO content
@@ -110,6 +113,9 @@
         'keep *_dedxTruncated40CosmicTF_*_*',
         'keep *_dedxHitInfoCosmicTF_*_*',
         'keep *_dedxHarmonic2CosmicTF_*_*',
+        'keep recoTracks_cosmicDCTracks_*_*',
+        'keep recoTrackExtras_cosmicDCTracks_*_*',
+        'keep TrackingRecHitsOwned_cosmicDCTracks_*_*',
     )
 )
 #AOD content

RecoTracker/Configuration/python/RecoTrackerP5_cff.py

@@ -77,3 +77,6 @@
 # REGIONAL RECONSTRUCTION
 from RecoTracker.Configuration.RecoTrackerNotStandard_cff import *
 regionalCosmicTrackerSeeds.RegionInJetsCheckPSet = cms.PSet( doJetsExclusionCheck   = cms.bool( False ) )
+
+# CDC Reconstruction
+from RecoTracker.SpecialSeedGenerators.cosmicDC_cff import *

RecoTracker/Configuration/python/RecoTracker_EventContent_cff.py

@@ -22,6 +22,9 @@
         'keep *_dedxHitInfo_*_*',
         'keep *_dedxHarmonic2_*_*',
         'keep *_trackExtrapolator_*_*',
+        'keep recoTracks_cosmicDCTracks_*_*',
+        'keep recoTrackExtras_cosmicDCTracks_*_*',
+        'keep TrackingRecHitsOwned_cosmicDCTracks_*_*',
      )
 )
 #RECO content
@@ -44,7 +47,10 @@
         'keep *_dedxTruncated40_*_*',
         'keep *_dedxHitInfo_*_*',
         'keep *_dedxHarmonic2_*_*',
-        'keep *_trackExtrapolator_*_*'
+        'keep *_trackExtrapolator_*_*',
+        'keep recoTracks_cosmicDCTracks_*_*',
+        'keep recoTrackExtras_cosmicDCTracks_*_*',
+        'keep TrackingRecHitsOwned_cosmicDCTracks_*_*',
     )
 )
 #AOD content

RecoTracker/SpecialSeedGenerators/python/cosmicDC_cff.py

@@ -0,0 +1,37 @@
+import FWCore.ParameterSet.Config as cms
+
+# seeding
+import RecoTracker.SpecialSeedGenerators.outInSeedsFromStandaloneMuons_cfi
+import TrackingTools.KalmanUpdators.Chi2MeasurementEstimator_cfi 
+hitCollectorForCosmicDCSeeds = TrackingTools.KalmanUpdators.Chi2MeasurementEstimator_cfi.Chi2MeasurementEstimator.clone(
+    ComponentName = cms.string('hitCollectorForCosmicDCSeeds'),
+    MaxChi2 = cms.double(100.0), ## was 30 ## TO BE TUNED
+    nSigma  = cms.double(4.),    ## was 3  ## TO BE TUNED 
+    MaxDisplacement = cms.double(100),
+    MaxSagitta = cms.double(-1.0),
+    MinimalTolerance = cms.double(0.5),
+    appendToDataLabel = cms.string(''),
+)
+cosmicDCSeeds = RecoTracker.SpecialSeedGenerators.outInSeedsFromStandaloneMuons_cfi.outInSeedsFromStandaloneMuons.clone(
+    src = cms.InputTag("muonsFromCosmics"),
+    cut = cms.string("pt > 2 && abs(eta)<1.2 && phi<0"),
+    hitCollector = cms.string('hitCollectorForCosmicDCSeeds'),
+    fromVertex = cms.bool(False),
+    maxEtaForTOB = cms.double(1.5),
+    minEtaForTEC = cms.double(0.7),
+)
+
+# Ckf pattern
+import RecoTracker.CkfPattern.CkfTrackCandidatesP5_cff
+cosmicDCCkfTrackCandidates = RecoTracker.CkfPattern.CkfTrackCandidatesP5_cff.ckfTrackCandidatesP5.clone(
+    src = cms.InputTag( "cosmicDCSeeds" ),
+)
+
+# Track producer
+import RecoTracker.TrackProducer.CTFFinalFitWithMaterialP5_cff
+cosmicDCTracks = RecoTracker.TrackProducer.CTFFinalFitWithMaterialP5_cff.ctfWithMaterialTracksCosmics.clone(
+    src = cms.InputTag( "cosmicDCCkfTrackCandidates" ),
+)
+
+# Final Sequence
+cosmicDCTracksSeq = cms.Sequence( cosmicDCSeeds * cosmicDCCkfTrackCandidates * cosmicDCTracks )