CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
phase2TkCustomsBE5D.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 #GEN-SIM so far...
3 def customise(process):
4  if hasattr(process,'DigiToRaw'):
5  process=customise_DigiToRaw(process)
6  if hasattr(process,'RawToDigi'):
7  process=customise_RawToDigi(process)
8  if hasattr(process,'reconstruction'):
9  process=customise_Reco(process)
10  if hasattr(process,'digitisation_step'):
11  process=customise_Digi(process)
12  if hasattr(process,'dqmoffline_step'):
13  process=customise_DQM(process)
14  if hasattr(process,'dqmHarvesting'):
15  process=customise_harvesting(process)
16  if hasattr(process,'validation_step'):
17  process=customise_Validation(process)
18  process=customise_condOverRides(process)
19 
20  return process
21 
22 def customise_Digi(process):
23  process.mix.digitizers.pixel.MissCalibrate = False
24  process.mix.digitizers.pixel.LorentzAngle_DB = False
25  process.mix.digitizers.pixel.killModules = False
26  process.mix.digitizers.pixel.useDB = False
27  process.mix.digitizers.pixel.DeadModules_DB = False
28  process.mix.digitizers.pixel.NumPixelBarrel = cms.int32(10)
29  process.mix.digitizers.pixel.NumPixelEndcap = cms.int32(7)
30  process.mix.digitizers.pixel.ThresholdInElectrons_FPix = cms.double(2000.0)
31  process.mix.digitizers.pixel.ThresholdInElectrons_BPix = cms.double(2000.0)
32  process.mix.digitizers.pixel.ThresholdInElectrons_BPix_L1 = cms.double(2000.0)
33  process.mix.digitizers.pixel.thePixelColEfficiency_BPix4 = cms.double(0.999)
34  process.mix.digitizers.pixel.thePixelEfficiency_BPix4 = cms.double(0.999)
35  process.mix.digitizers.pixel.thePixelChipEfficiency_BPix4 = cms.double(0.999)
36  process.mix.digitizers.pixel.thePixelColEfficiency_FPix3 = cms.double(0.999)
37  process.mix.digitizers.pixel.thePixelEfficiency_FPix3 = cms.double(0.999)
38  process.mix.digitizers.pixel.thePixelChipEfficiency_FPix3 = cms.double(0.999)
39  process.mix.digitizers.pixel.AddPixelInefficiencyFromPython = cms.bool(False)
40  process.mix.digitizers.strip.ROUList = cms.vstring("g4SimHitsTrackerHitsPixelBarrelLowTof",
41  'g4SimHitsTrackerHitsPixelEndcapLowTof')
42  process.digitisation_step.remove(process.simSiStripDigiSimLink)
43  process.mergedtruth.simHitCollections.tracker = []
44  return process
45 
46 
47 def customise_DigiToRaw(process):
48  process.digi2raw_step.remove(process.siPixelRawData)
49  process.digi2raw_step.remove(process.rpcpacker)
50  return process
51 
52 def customise_RawToDigi(process):
53  process.raw2digi_step.remove(process.siPixelDigis)
54  return process
55 
56 def customise_Reco(process):
57  #this may be a trimmed out process with only local reco
58  #if so, don't use the customize stuff
59  ## need changes to mixedtriplets step to use for imcreasing high eta efficiency
60  process.reconstruction.remove(process.pixelLessStepSeedClusterMask)
61  process.reconstruction.remove(process.castorreco)
62  process.reconstruction.remove(process.CastorTowerReco)
63  process.reconstruction.remove(process.ak7BasicJets)
64  process.reconstruction.remove(process.ak7CastorJetID)
65  #process.iterTracking.remove(process.PixelLessStep)
66  #process.iterTracking.remove(process.TobTecStep)
67  process.MixedTripletStep.remove(process.mixedTripletStepSeedsB)
68  process.mixedTripletStepSeeds = cms.EDProducer("SeedCombiner",
69  seedCollections = cms.VInputTag(cms.InputTag("mixedTripletStepSeedsA"))
70  )
71  process.load("RecoTracker.IterativeTracking.HighPtTripletStep_cff")
72  from RecoTracker.IterativeTracking.HighPtTripletStep_cff import HighPtTripletStep
73  process.iterTracking = cms.Sequence(process.InitialStep*
74  process.HighPtTripletStep*
75  process.LowPtTripletStep*
76  process.PixelPairStep*
77  process.DetachedTripletStep*
78  process.MixedTripletStep*
79  #process.PixelLessStep*
80  #process.TobTecStep*
81  process.earlyGeneralTracks*
82  process.muonSeededStep*
83  process.preDuplicateMergingGeneralTracks*
84  process.generalTracksSequence*
85  process.ConvStep*
86  process.conversionStepTracks
87  )
88 
89 
90 
91 
92  process.convClusters.oldClusterRemovalInfo=cms.InputTag("mixedTripletStepClusters")
93  process.convClusters.trajectories=cms.InputTag("mixedTripletStepTracks")
94  process.convClusters.overrideTrkQuals= cms.InputTag("mixedTripletStep")
95  process.mixedTripletStepSeedLayersA.layerList = cms.vstring('BPix1+BPix2+BPix3',
96  'BPix1+BPix2+FPix1_pos',
97  'BPix1+BPix2+FPix1_neg',
98  'BPix1+FPix1_pos+FPix2_pos',
99  'BPix1+FPix1_neg+FPix2_neg',
100  'BPix2+FPix1_pos+FPix2_pos',
101  'BPix2+FPix1_neg+FPix2_neg')
102 
103  process.earlyGeneralTracks.setsToMerge = cms.VPSet( cms.PSet( tLists=cms.vint32(0,1,2,3,4), pQual=cms.bool(True) ))
104  process.earlyGeneralTracks.hasSelector=cms.vint32(1,1,1,1,1)
105  process.earlyGeneralTracks.selectedTrackQuals = cms.VInputTag(
106  cms.InputTag("initialStepSelector","initialStep"),
107  cms.InputTag("highPtTripletStepSelector","highPtTripletStep"),
108  cms.InputTag("lowPtTripletStepSelector","lowPtTripletStep"),
109  cms.InputTag("pixelPairStepSelector","pixelPairStep"),
110  cms.InputTag("mixedTripletStep")
111  )
112  process.earlyGeneralTracks.TrackProducers = cms.VInputTag(
113  cms.InputTag("initialStepTracks"),
114  cms.InputTag("highPtTripletStepTracks"),
115  cms.InputTag("lowPtTripletStepTracks"),
116  cms.InputTag("pixelPairStepTracks"),
117  cms.InputTag("mixedTripletStepTracks")
118  )
119 
120  process.tripletElectronSeedLayers.BPix.skipClusters=cms.InputTag('mixedTripletStepSeedClusterMask')
121  process.tripletElectronSeedLayers.FPix.skipClusters=cms.InputTag('mixedTripletStepSeedClusterMask')
122  process.tripletElectronClusterMask.oldClusterRemovalInfo=cms.InputTag('mixedTripletStepSeedClusterMask')
123 
124  process.initialStepSeedClusterMask.oldClusterRemovalInfo=cms.InputTag("mixedTripletStepClusters")
125  process.newCombinedSeeds.seedCollections = cms.VInputTag(cms.InputTag('initialStepSeeds'),
126  cms.InputTag("highPtTripletStepSeeds"),
127  cms.InputTag('pixelPairStepSeeds'),
128  cms.InputTag('mixedTripletStepSeeds'),
129  cms.InputTag('tripletElectronSeeds'),
130  cms.InputTag('pixelPairElectronSeeds'),
131  cms.InputTag('stripPairElectronSeeds') )
132  process.stripPairElectronSeedLayers.layerList = cms.vstring('BPix4+BPix5') # Optimize later
133  process.stripPairElectronSeedLayers.BPix = cms.PSet(
134  HitProducer = cms.string('siPixelRecHits'),
135  TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4PixelPairs'),
136  skipClusters = cms.InputTag("pixelPairStepClusters"),
137  )
138  process.regionalCosmicTrackerSeeds.OrderedHitsFactoryPSet.LayerPSet.layerList = cms.vstring('BPix10+BPix9') # Optimize later
139  process.regionalCosmicTrackerSeeds.OrderedHitsFactoryPSet.LayerPSet.BPix = cms.PSet(
140  HitProducer = cms.string('siPixelRecHits'),
141  TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4PixelPairs'),
142  skipClusters = cms.InputTag("pixelPairStepClusters"),
143  )
144  from RecoPixelVertexing.PixelTriplets.quadrupletseedmerging_cff import PixelSeedMergerQuadruplets
145  process.pixelTracks.SeedMergerPSet = cms.PSet(
146  layerList = PixelSeedMergerQuadruplets,
147  addRemainingTriplets = cms.bool(False),
148  mergeTriplets = cms.bool(True),
149  ttrhBuilderLabel = cms.string('PixelTTRHBuilderWithoutAngle')
150  )
151  process.initialStepSeedClusterMask.oldClusterRemovalInfo=cms.InputTag("mixedTripletStepClusters")
152 
153  # Need this line to stop error about missing siPixelDigis.
154  process.MeasurementTracker.inactivePixelDetectorLabels = cms.VInputTag()
155  process.load("SLHCUpgradeSimulations.Geometry.recoFromSimDigis_cff")
156  # Use with latest pixel geometry. Only used for seeds, so we can use the Phase1Tk file.
157  # We will need to turn it off for any steps that use the outer pixels as seeds.
158  process.ClusterShapeHitFilterESProducer.PixelShapeFile = cms.string('RecoPixelVertexing/PixelLowPtUtilities/data/pixelShape_Phase1Tk.par')
159  # Now make sure we us CPE Generic
160  process.mergedDuplicateTracks.TTRHBuilder = 'WithTrackAngle'
161  process.ctfWithMaterialTracks.TTRHBuilder = 'WithTrackAngle'
162  process.PixelCPEGenericESProducer.UseErrorsFromTemplates = cms.bool(False)
163  process.PixelCPEGenericESProducer.TruncatePixelCharge = cms.bool(False)
164  process.PixelCPEGenericESProducer.LoadTemplatesFromDB = cms.bool(False)
165  process.PixelCPEGenericESProducer.Upgrade = cms.bool(True)
166  process.PixelCPEGenericESProducer.IrradiationBiasCorrection = False
167  process.PixelCPEGenericESProducer.DoCosmics = False
168  process.siPixelRecHits.CPE = cms.string('PixelCPEGeneric')
169  #the quadruplet merger configuration
170  # from this PSet the quadruplet merger uses only the layer list so these could probably be removed
171  PixelSeedMergerQuadruplets.BPix.TTRHBuilder = cms.string("PixelTTRHBuilderWithoutAngle" )
172  PixelSeedMergerQuadruplets.BPix.HitProducer = cms.string("siPixelRecHits" )
173  PixelSeedMergerQuadruplets.FPix.TTRHBuilder = cms.string("PixelTTRHBuilderWithoutAngle" )
174  PixelSeedMergerQuadruplets.FPix.HitProducer = cms.string("siPixelRecHits" )
175 
176  process.highPtTripletStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
177  process.detachedTripletStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
178  process.initialStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
179  process.pixelPairStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
180  process.lowPtTripletStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
181  process.convStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
182  process.mixedTripletStepTracks.TTRHBuilder=cms.string('WithTrackAngle')
183  process.muonSeededSeedsInOut.TrackerRecHitBuilder = cms.string('WithTrackAngle')
184  process.muonSeededTracksInOut.TTRHBuilder = cms.string('WithTrackAngle')
185  process.muons1stStep.TrackerKinkFinderParameters.TrackerRecHitBuilder=cms.string('WithTrackAngle')
186  process.regionalCosmicTracks.TTRHBuilder=cms.string('WithTrackAngle')
187  process.cosmicsVetoTracksRaw.TTRHBuilder=cms.string('WithTrackAngle')
188 
189  return process
190 
191 def customise_condOverRides(process):
192  process.load('SLHCUpgradeSimulations.Geometry.fakeConditions_BarrelEndcap5D_cff')
193  process.trackerNumberingSLHCGeometry.layerNumberPXB = cms.uint32(20)
194  process.trackerTopologyConstants.pxb_layerStartBit = cms.uint32(20)
195  process.trackerTopologyConstants.pxb_ladderStartBit = cms.uint32(12)
196  process.trackerTopologyConstants.pxb_moduleStartBit = cms.uint32(2)
197  process.trackerTopologyConstants.pxb_layerMask = cms.uint32(15)
198  process.trackerTopologyConstants.pxb_ladderMask = cms.uint32(255)
199  process.trackerTopologyConstants.pxb_moduleMask = cms.uint32(1023)
200  process.trackerTopologyConstants.pxf_diskStartBit = cms.uint32(18)
201  process.trackerTopologyConstants.pxf_bladeStartBit = cms.uint32(12)
202  process.trackerTopologyConstants.pxf_panelStartBit = cms.uint32(10)
203  process.trackerTopologyConstants.pxf_moduleMask = cms.uint32(255)
204  return process
205 
206 
207 def l1EventContent(process):
208  #extend the event content
209 
210  alist=['RAWSIM','FEVTDEBUG','FEVTDEBUGHLT','GENRAW','RAWSIMHLT','FEVT']
211  for a in alist:
212  b=a+'output'
213  if hasattr(process,b):
214  getattr(process,b).outputCommands.append('keep PSimHits_g4SimHits_*_*')
215  getattr(process,b).outputCommands.append('keep SimTracks_g4SimHits_*_*')
216  getattr(process,b).outputCommands.append('keep SimVertexs_g4SimHits_*_*')
217  getattr(process,b).outputCommands.append('keep *_simSiPixelDigis_*_*')
218  getattr(process,b).outputCommands.append('keep *_genParticles_*_*')
219  getattr(process,b).outputCommands.append('keep *_L1TkBeams_*_*')
220  getattr(process,b).outputCommands.append('keep *_L1TkClustersFromPixelDigis_*_*')
221  getattr(process,b).outputCommands.append('keep *_L1TkClustersFromSimHits_*_*')
222  getattr(process,b).outputCommands.append('keep *_L1TkStubsFromPixelDigis_*_*')
223  getattr(process,b).outputCommands.append('keep *_L1TkStubsFromSimHits_*_*')
224  getattr(process,b).outputCommands.append('keep *_siPixelRecHits_*_*')
225  #drop some bigger collections we don't think we need
226  getattr(process,b).outputCommands.append('drop PSimHits_g4SimHits_EcalHitsEB_*')
227  getattr(process,b).outputCommands.append('drop PSimHits_g4SimHits_EcalHitsEE_*')
228  getattr(process,b).outputCommands.append('drop *_L1TkStubsFromSimHits_StubsFail_*')
229  return process
230 
HCalCustoms.customise_Validation
def customise_Validation
Definition: HCalCustoms.py:237
HCalCustoms.customise_harvesting
def customise_harvesting
Definition: HCalCustoms.py:234
HCalCustoms.customise_DQM
def customise_DQM
process.load("RecoLocalCalo.HcalRecProducers.HcalSimpleReconstructor_ho_cfi")
Definition: HCalCustoms.py:218
quadrupletseedmerging_cff
Generated for CMSSW Reference Manual by doxygen 1.8.5