00001 # AlCaReco for muon based alignment using beam-halo muons in the CSC overlap regions 00002 00003 import FWCore.ParameterSet.Config as cms 00004 00005 # import HLTrigger.HLTfilters.hltHighLevel_cfi 00006 # ALCARECOMuAlBeamHaloOverlapsHLT = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone( 00007 # HLTPaths = ['HLT_CSCBeamHaloOverlapRing1', 'HLT_CSCBeamHaloOverlapRing2'], 00008 # throw = False # tolerate triggers stated above, but not available 00009 # ) 00010 00011 # from RecoMuon.Configuration.RecoMuonCosmics_cff import * # required as of 2_1_X 00012 # from Geometry.CommonDetUnit.bareGlobalTrackingGeometry_cfi import * 00013 # from RecoMuon.DetLayers.muonDetLayerGeometry_cfi import * 00014 00015 # import RecoLocalMuon.CSCSegment.cscSegments_cfi 00016 # cscSegmentsALCARECOBHO = RecoLocalMuon.CSCSegment.cscSegments_cfi.cscSegments.clone() 00017 # cscSegmentsALCARECOBHO.algo_type = 4 # Choice of the building algo: 1 SK, 2 TC, 3 DF, 4 ST, ... 00018 # cscSegmentsALCARECOBHO.inputObjects = 'csc2DRecHits' 00019 00020 # import RecoMuon.MuonSeedGenerator.CosmicMuonSeedProducer_cfi 00021 # CosmicMuonSeedALCARECOBHO = RecoMuon.MuonSeedGenerator.CosmicMuonSeedProducer_cfi.CosmicMuonSeed.clone() 00022 # CosmicMuonSeedALCARECOBHO.EnableDTMeasurement = False 00023 # CosmicMuonSeedALCARECOBHO.CSCRecSegmentLabel = 'cscSegmentsALCARECOBHO' 00024 00025 # import RecoMuon.CosmicMuonProducer.cosmicMuons_cfi 00026 # cosmicMuonsALCARECOBHO = RecoMuon.CosmicMuonProducer.cosmicMuons_cfi.cosmicMuons.clone() 00027 # cosmicMuonsALCARECOBHO.MuonSeedCollectionLabel = 'CosmicMuonSeedALCARECOBHO' 00028 # cosmicMuonsALCARECOBHO.TrajectoryBuilderParameters.EnableDTMeasurement = False 00029 # cosmicMuonsALCARECOBHO.TrajectoryBuilderParameters.CSCRecSegmentLabel = 'cscSegmentsALCARECOBHO' 00030 00031 # reconstructAsCosmicMuonsALCARECOBHO = cms.Sequence(cscSegmentsALCARECOBHO * CosmicMuonSeedALCARECOBHO * cosmicMuonsALCARECOBHO) 00032 00033 ALCARECOMuAlBeamHaloOverlapsEnergyCut = cms.EDFilter("AlignmentCSCBeamHaloSelectorModule", 00034 filter = cms.bool(True), 00035 # src = cms.InputTag("cosmicMuonsALCARECOBHO"), 00036 src = cms.InputTag("cosmicMuons"), # get cosmicMuons from global-run reconstruction 00037 minStations = cms.uint32(0), # no "energy cut" yet 00038 minHitsPerStation = cms.uint32(4) 00039 ) 00040 00041 aLCARECOMuAlBeamHaloOverlaps = cms.EDFilter("AlignmentCSCOverlapSelectorModule", 00042 filter = cms.bool(True), 00043 src = cms.InputTag("ALCARECOMuAlBeamHaloOverlapsEnergyCut"), 00044 minHitsPerChamber = cms.uint32(4), 00045 station = cms.int32(0) # all stations: I'll need to split it by station (8 subsamples) offline 00046 ) 00047 00048 # seqALCARECOMuAlBeamHaloOverlaps = cms.Sequence(ALCARECOMuAlBeamHaloOverlapsHLT + reconstructAsCosmicMuonsALCARECOBHO * ALCARECOMuAlBeamHaloOverlapsEnergyCut * aLCARECOMuAlBeamHaloOverlaps) 00049 # seqALCARECOMuAlBeamHaloOverlaps = cms.Sequence(reconstructAsCosmicMuonsALCARECOBHO * ALCARECOMuAlBeamHaloOverlapsEnergyCut * aLCARECOMuAlBeamHaloOverlaps) 00050 seqALCARECOMuAlBeamHaloOverlaps = cms.Sequence(ALCARECOMuAlBeamHaloOverlapsEnergyCut * aLCARECOMuAlBeamHaloOverlaps) 00051