TkAlCaRecoProducers/test/Alca_BeamFit_Workflow.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("alcaBeamSpotWorkflow")
0004 # initialize MessageLogger
0005 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0006
0007 process.MessageLogger.cerr = cms.untracked.PSet(enable = cms.untracked.bool(False))
0008 process.MessageLogger.cout = cms.untracked.PSet(
0009     threshold = cms.untracked.string('INFO'),
0010     default = cms.untracked.PSet(
0011        limit = cms.untracked.int32(0)
0012     ),
0013     AlcaBeamSpotProducer = cms.untracked.PSet(
0014         #reportEvery = cms.untracked.int32(100) # every 1000th only
0015     limit = cms.untracked.int32(0)
0016     )
0017 )
0018 #process.MessageLogger.cout.enableStatistics = cms.untracked.bool(True)
0019
0020 process.load("Calibration.TkAlCaRecoProducers.AlcaBeamSpotProducer_cff")
0021
0022 process.source = cms.Source("PoolSource",
0023     fileNames = cms.untracked.vstring(
0024 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0171/3A22D08D-456A-DF11-961F-001A92811734.root',
0025 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0171/083B6802-236C-DF11-8AC6-0026189437FE.root',
0026 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0170/E6D0589B-136A-DF11-9E90-002618943982.root',
0027 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0170/D098488B-276A-DF11-8069-003048678AF4.root',
0028 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/F466DD94-BF69-DF11-B9B8-00261894390A.root',
0029 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/F2A9245B-BC69-DF11-9C3F-0018F3D096E4.root',
0030 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/DED5E502-B969-DF11-AE8F-002618943964.root',
0031 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/563257AF-C169-DF11-865F-002618943907.root',
0032 '/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/3056B847-C069-DF11-B0B0-003048679010.root'
0033
0034     )
0035 )
0036
0037 process.maxEvents = cms.untracked.PSet(
0038     input = cms.untracked.int32(10000) #1500
0039 )
0040
0041 process.options = cms.untracked.PSet(
0042     wantSummary = cms.untracked.bool(True)
0043 )
0044
0045 # this is for filtering on L1 technical trigger bit
0046 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0047 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0048 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0049 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND ( 40 OR 41 ) AND NOT (36 OR 37 OR 38 OR 39)')
0050
0051 ##
0052 process.load("Configuration.StandardSequences.MagneticField_cff")
0053 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0054 process.GlobalTag.globaltag = 'GR_R_38X_V9::All' #'GR_R_35X_V8::All'
0055 process.load("Configuration.StandardSequences.GeometryDB_cff")
0056
0057
0058 ## reco PV
0059 process.load("Configuration.StandardSequences.Reconstruction_cff")
0060 process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
0061 process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
0062 process.offlinePrimaryVertices.TrackLabel = cms.InputTag("ALCARECOTkAlMinBias")
0063
0064 #### remove beam scraping events
0065 process.noScraping= cms.EDFilter("FilterOutScraping",
0066                                  applyfilter = cms.untracked.bool(True),
0067                                  debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0068                                  numtrack = cms.untracked.uint32(10),
0069                                  thresh = cms.untracked.double(0.20)
0070 )
0071
0072
0073
0074 ################### Primary Vertex
0075 process.offlinePrimaryVertices.PVSelParameters.maxDistanceToBeam = 2
0076 process.offlinePrimaryVertices.TkFilterParameters.maxNormalizedChi2 = 20
0077 process.offlinePrimaryVertices.TkFilterParameters.minSiliconLayersWithHits = 5
0078 process.offlinePrimaryVertices.TkFilterParameters.maxD0Significance = 100
0079 process.offlinePrimaryVertices.TkFilterParameters.minPixelLayersWithHits = 1
0080 process.offlinePrimaryVertices.TkClusParameters.TkGapClusParameters.zSeparation = 1
0081
0082
0083 #######################
0084 process.alcaBeamSpotProducer.BeamFitter.TrackCollection = 'ALCARECOTkAlMinBias'
0085 process.alcaBeamSpotProducer.BeamFitter.MinimumTotalLayers = 6
0086 process.alcaBeamSpotProducer.BeamFitter.MinimumPixelLayers = -1
0087 process.alcaBeamSpotProducer.BeamFitter.MaximumNormChi2 = 10
0088 process.alcaBeamSpotProducer.BeamFitter.MinimumInputTracks = 50
0089 process.alcaBeamSpotProducer.BeamFitter.MinimumPt = 1.0
0090 process.alcaBeamSpotProducer.BeamFitter.MaximumImpactParameter = 1.0
0091 process.alcaBeamSpotProducer.BeamFitter.TrackAlgorithm =  cms.untracked.vstring()
0092 #process.alcaBeamSpotProducer.BeamFitter.Debug = True
0093
0094 process.alcaBeamSpotProducer.PVFitter.Apply3DFit = True
0095 process.alcaBeamSpotProducer.PVFitter.minNrVerticesForFit = 10
0096 #########################
0097
0098
0099 # fit as function of lumi sections
0100 process.alcaBeamSpotProducer.AlcaBeamSpotProducerParameters.fitEveryNLumi = 1
0101 process.alcaBeamSpotProducer.AlcaBeamSpotProducerParameters.resetEveryNLumi = 1
0102
0103 process.out = cms.OutputModule( "PoolOutputModule",
0104                                 fileName = cms.untracked.string( 'AlcaBeamSpot.root' ),
0105                                 outputCommands = cms.untracked.vstring("keep *")
0106                               )
0107
0108
0109 process.e = cms.EndPath( process.out )
0110
0111 process.p = cms.Path(process.hltLevel1GTSeed +
0112                      process.offlineBeamSpot +
0113 #                     process.TrackRefitter +
0114                      process.offlinePrimaryVertices+
0115 #                     process.noScraping +
0116                      process.alcaBeamSpotProducer)