Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:03:38

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 def customiseCommon(process):
0004 
0005     
0006     #####################################################################################################
0007     ####
0008     ####  Top level replaces for handling strange scenarios of early collisions
0009     ####
0010 
0011     ## TRACKING:
0012     ## Skip events with HV off
0013     process.newSeedFromTriplets.ClusterCheckPSet.MaxNumberOfPixelClusters=2000
0014     process.newSeedFromPairs.ClusterCheckPSet.MaxNumberOfStripClusters=10000
0015     process.secTriplets.ClusterCheckPSet.MaxNumberOfPixelClusters=2000
0016     process.fifthSeeds.ClusterCheckPSet.MaxNumberOfStripClusters = 10000
0017     process.fourthPLSeeds.ClusterCheckPSet.MaxNumberOfStripClusters=10000
0018 
0019     ###### FIXES TRIPLETS FOR LARGE BS DISPLACEMENT ######
0020 
0021     ### prevent bias in pixel vertex
0022     process.pixelVertices.useBeamConstraint = False
0023     
0024     ### pixelTracks
0025     #---- new parameters ----
0026     process.pixelTracks.RegionFactoryPSet.RegionPSet.nSigmaZ  = cms.double(4.06) # was originHalfLength = 15.9; translated assuming sigmaZ ~ 3.8
0027     
0028     ### 0th step of iterative tracking
0029     #---- replaces ----
0030     process.newSeedFromTriplets.RegionFactoryPSet.ComponentName = 'GlobalRegionProducerFromBeamSpot' # was GlobalRegionProducer
0031     #---- new parameters ----
0032     process.newSeedFromTriplets.RegionFactoryPSet.RegionPSet.nSigmaZ   = cms.double(4.06)  # was originHalfLength = 15.9; translated assuming sigmaZ ~ 3.8
0033     process.newSeedFromTriplets.RegionFactoryPSet.RegionPSet.beamSpot = cms.InputTag("offlineBeamSpot")
0034 
0035     ### 2nd step of iterative tracking
0036     #---- replaces ----
0037     process.secTriplets.RegionFactoryPSet.ComponentName = 'GlobalRegionProducerFromBeamSpot' # was GlobalRegionProducer
0038     #---- new parameters ----
0039     process.secTriplets.RegionFactoryPSet.RegionPSet.nSigmaZ  = cms.double(4.47)  # was originHalfLength = 17.5; translated assuming sigmaZ ~ 3.8
0040     process.secTriplets.RegionFactoryPSet.RegionPSet.beamSpot = cms.InputTag("offlineBeamSpot")
0041 
0042     ## Primary Vertex
0043     process.offlinePrimaryVerticesWithBS.PVSelParameters.maxDistanceToBeam = 2
0044     process.offlinePrimaryVerticesWithBS.TkFilterParameters.maxNormalizedChi2 = 20
0045     process.offlinePrimaryVerticesWithBS.TkFilterParameters.minSiliconHits = 6
0046     process.offlinePrimaryVerticesWithBS.TkFilterParameters.maxD0Significance = 100
0047     process.offlinePrimaryVerticesWithBS.TkFilterParameters.minPixelHits = 1
0048     process.offlinePrimaryVerticesWithBS.TkClusParameters.zSeparation = 1
0049     process.offlinePrimaryVertices.PVSelParameters.maxDistanceToBeam = 2
0050     process.offlinePrimaryVertices.TkFilterParameters.maxNormalizedChi2 = 20
0051     process.offlinePrimaryVertices.TkFilterParameters.minSiliconHits = 6
0052     process.offlinePrimaryVertices.TkFilterParameters.maxD0Significance = 100
0053     process.offlinePrimaryVertices.TkFilterParameters.minPixelHits = 1
0054     process.offlinePrimaryVertices.TkClusParameters.zSeparation = 1
0055 
0056     ## ECAL 
0057     process.ecalRecHit.ChannelStatusToBeExcluded = [ 1, 2, 3, 4, 8, 9, 10, 11, 12, 13, 14, 78, 142 ]
0058 
0059     ## HCAL temporary fixes
0060     process.hfreco.samplesToAdd = 4
0061     
0062     ## EGAMMA
0063     process.gsfElectrons.applyPreselection = cms.bool(False)
0064     process.photons.minSCEtBarrel = 2.
0065     process.photons.minSCEtEndcap =2.
0066     process.photonCore.minSCEt = 2.
0067     process.conversionTrackCandidates.minSCEt =1.
0068     process.conversions.minSCEt =1.
0069     process.trackerOnlyConversions.AllowTrackBC = cms.bool(False)
0070     process.trackerOnlyConversions.AllowRightBC = cms.bool(False)
0071     process.trackerOnlyConversions.MinApproach = cms.double(-.25)
0072     process.trackerOnlyConversions.DeltaCotTheta = cms.double(.07)
0073     process.trackerOnlyConversions.DeltaPhi = cms.double(.2)
0074     
0075     ###
0076     ###  end of top level replacements
0077     ###
0078     ###############################################################################################
0079 
0080     return (process)
0081 
0082 
0083 ##############################################################################
0084 def customisePPData(process):
0085     process= customiseCommon(process)
0086     process.hfreco.firstSample=3
0087     ##Preshower
0088     process.ecalPreshowerRecHit.ESBaseline = 0
0089 
0090     ##Preshower algo for data is different than for MC
0091     process.ecalPreshowerRecHit.ESRecoAlgo = cms.untracked.int32(1)
0092 
0093     return process
0094 
0095 
0096 ##############################################################################
0097 def customisePPMC(process):
0098     process=customiseCommon(process)
0099     process.hfreco.firstSample=1
0100     
0101     return process
0102 
0103 ##############################################################################
0104 def customiseCosmicData(process):
0105     process.ecalPreshowerRecHit.ESBaseline = cms.int32(0)
0106     process.ecalPreshowerRecHit.ESRecoAlgo = cms.untracked.int32(1)
0107     
0108     return process
0109 
0110 ##############################################################################
0111 def customiseCosmicMC(process):
0112     
0113     return process
0114         
0115 
0116 ##############################################################################
0117 def customiseExpress(process):
0118     process= customisePPData(process)
0119 
0120     return process