BeamSpotProducer/test/BeamFit_LumiBased_Workflow.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("BSworkflow")
0004 # initialize MessageLogger
0005 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0006
0007 process.load("RecoVertex.BeamSpotProducer.d0_phi_analyzer_cff")
0008
0009 process.source = cms.Source("PoolSource",
0010     fileNames = cms.untracked.vstring(
0011 '/store/express/Run2010B/StreamExpress/ALCARECO/TkAlMinBias-v2/000/147/984/00B7AE46-58D8-DF11-9A23-001D09F292D1.root'
0012 )
0013 )
0014
0015 process.MessageLogger.cerr.FwkReport  = cms.untracked.PSet(
0016     reportEvery = cms.untracked.int32(1000000),
0017 )
0018
0019 #process.source = cms.Source('PoolSource',
0020 #                            debugVerbosity = cms.untracked.uint32(0),
0021 #                            debugFlag = cms.untracked.bool(False)
0022 #                            )
0023
0024 process.maxEvents = cms.untracked.PSet(
0025     input = cms.untracked.int32(-1) #1500
0026 )
0027
0028 process.options = cms.untracked.PSet(
0029     wantSummary = cms.untracked.bool(True)
0030     )
0031
0032 # this is for filtering on L1 technical trigger bit
0033 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0034 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0035 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0036 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND ( 40 OR 41 ) AND NOT (36 OR 37 OR 38 OR 39)')
0037
0038 ## reco PV
0039 process.load("Configuration.StandardSequences.MagneticField_cff")
0040 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0041 process.GlobalTag.globaltag = 'GR_R_38X_V11::All'
0042
0043 process.load("Configuration.StandardSequences.Reconstruction_cff")
0044 process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
0045 process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
0046 process.offlinePrimaryVertices.TrackLabel = cms.InputTag("ALCARECOTkAlMinBias")
0047
0048 #### remove beam scraping events
0049 process.noScraping= cms.EDFilter("FilterOutScraping",
0050                                  applyfilter = cms.untracked.bool(True),
0051                                  debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0052                                  numtrack = cms.untracked.uint32(10),
0053                                  thresh = cms.untracked.double(0.20)
0054 )
0055
0056 process.p = cms.Path(
0057 #    process.hltLevel1GTSeed +
0058 #    process.offlineBeamSpot +
0059 #    process.offlinePrimaryVertices+
0060 #    process.noScraping +
0061     process.d0_phi_analyzer)
0062
0063 process.MessageLogger.debugModules = ['BeamSpotAnalyzer']
0064
0065 ################### Primary Vertex
0066 process.offlinePrimaryVertices.PVSelParameters.maxDistanceToBeam = 2
0067 process.offlinePrimaryVertices.TkFilterParameters.maxNormalizedChi2 = 20
0068 process.offlinePrimaryVertices.TkFilterParameters.minSiliconLayersWithHits = 5
0069 process.offlinePrimaryVertices.TkFilterParameters.maxD0Significance = 100
0070 process.offlinePrimaryVertices.TkFilterParameters.minPixelLayersWithHits = 1
0071 process.offlinePrimaryVertices.TkClusParameters.TkGapClusParameters.zSeparation = 1
0072
0073
0074 #######################
0075 process.d0_phi_analyzer.BeamFitter.TrackCollection = 'ALCARECOTkAlMinBias'
0076 process.d0_phi_analyzer.BeamFitter.MinimumTotalLayers = 6
0077 process.d0_phi_analyzer.BeamFitter.MinimumPixelLayers = -1
0078 process.d0_phi_analyzer.BeamFitter.MaximumNormChi2 = 10
0079 process.d0_phi_analyzer.BeamFitter.MinimumInputTracks = 50
0080 process.d0_phi_analyzer.BeamFitter.MinimumPt = 1.0
0081 process.d0_phi_analyzer.BeamFitter.MaximumImpactParameter = 1.0
0082 process.d0_phi_analyzer.BeamFitter.TrackAlgorithm =  cms.untracked.vstring()
0083 #process.d0_phi_analyzer.BeamFitter.Debug = True
0084
0085 process.d0_phi_analyzer.PVFitter.Apply3DFit = True
0086 process.d0_phi_analyzer.PVFitter.minNrVerticesForFit = 10
0087 #########################
0088
0089 process.d0_phi_analyzer.BeamFitter.AsciiFileName = 'BeamFit_LumiBased_Workflow.txt'
0090 process.d0_phi_analyzer.BeamFitter.AppendRunToFileName = False
0091 process.d0_phi_analyzer.BeamFitter.OutputFileName = 'BeamFit_LumiBased_Workflow.root'
0092 #process.d0_phi_analyzer.BeamFitter.SaveNtuple = True
0093 process.d0_phi_analyzer.BeamFitter.SavePVVertices = True
0094
0095 # fit as function of lumi sections
0096 process.d0_phi_analyzer.BSAnalyzerParameters.fitEveryNLumi = 1
0097 process.d0_phi_analyzer.BSAnalyzerParameters.resetEveryNLumi = 1