BeamSpotProducer/test/BeamFit_RunBased_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/data/BeamCommissioning09/MinimumBias/RECO/Dec19thReReco_341_v1/0002/FE0A36FF-11ED-DE11-A8EE-002618943849.root',
0012 #    '/store/data/BeamCommissioning09/MinimumBias/RECO/Dec19thReReco_341_v1/0002/F89D50B7-07ED-DE11-B079-0026189437ED.root',
0013 #    '/store/data/BeamCommissioning09/MinimumBias/RECO/Dec19thReReco_341_v1/0002/F0E44CBD-07ED-DE11-A0F7-002618943914.root',
0014 #    '/store/data/BeamCommissioning09/MinimumBias/RECO/Dec19thReReco_341_v1/0002/DAB587B8-07ED-DE11-8369-00304867C1BC.root',
0015 #    '/store/data/BeamCommissioning09/MinimumBias/RECO/Dec19thReReco_341_v1/0002/C8E03DBC-07ED-DE11-88ED-00248C0BE005.root'
0016     '/store/data/BeamCommissioning09/MinimumBias/ALCARECO/StreamTkAlMinBias-Dec19thReReco_341_v1/0006/DEC0A7AE-2CEE-DE11-A801-002618943945.root',
0017     '/store/data/BeamCommissioning09/MinimumBias/ALCARECO/StreamTkAlMinBias-Dec19thReReco_341_v1/0006/D49332EC-F5ED-DE11-8F55-002618943950.root'
0018     )
0019 )
0020
0021 process.maxEvents = cms.untracked.PSet(
0022     input = cms.untracked.int32(-1) #1500
0023 )
0024
0025 process.options = cms.untracked.PSet(
0026     wantSummary = cms.untracked.bool(True)
0027     )
0028
0029 # this is for filtering on L1 technical trigger bit
0030 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0031 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0032 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0033 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND ( 40 OR 41 )')
0034 ##
0035
0036 ## reco PV
0037 process.load("Configuration.StandardSequences.MagneticField_cff")
0038 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0039 process.GlobalTag.globaltag = 'GR10_P_V2::All'
0040
0041 process.load("Configuration.StandardSequences.Reconstruction_cff")
0042 process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
0043 process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
0044 process.offlinePrimaryVertices.TrackLabel = cms.InputTag("ALCARECOTkAlMinBias")
0045
0046 #### remove beam scraping events
0047 process.noScraping= cms.EDFilter("FilterOutScraping",
0048                                  applyfilter = cms.untracked.bool(True),
0049                                  debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0050                                  numtrack = cms.untracked.uint32(10),
0051                                  thresh = cms.untracked.double(0.20)
0052 )
0053
0054 process.p = cms.Path(process.hltLevel1GTSeed +
0055                      process.offlineBeamSpot +
0056                      process.offlinePrimaryVertices+
0057 #                     process.noScraping +
0058                      process.d0_phi_analyzer)
0059
0060 process.MessageLogger.debugModules = ['BeamSpotAnalyzer']
0061
0062 #######################
0063 process.d0_phi_analyzer.BeamFitter.TrackCollection = 'ALCARECOTkAlMinBias'
0064 process.d0_phi_analyzer.BeamFitter.MinimumTotalLayers = 6
0065 process.d0_phi_analyzer.BeamFitter.MinimumPixelLayers = -1
0066 process.d0_phi_analyzer.BeamFitter.MaximumNormChi2 = 10
0067 process.d0_phi_analyzer.BeamFitter.MinimumInputTracks = 2
0068 process.d0_phi_analyzer.BeamFitter.MinimumPt = 1.0
0069 process.d0_phi_analyzer.BeamFitter.MaximumImpactParameter = 1.0
0070 process.d0_phi_analyzer.BeamFitter.TrackAlgorithm =  cms.untracked.vstring()
0071 process.d0_phi_analyzer.BeamFitter.InputBeamWidth = -1 # 0.0400
0072 process.d0_phi_analyzer.BeamFitter.Debug = True
0073
0074 process.d0_phi_analyzer.PVFitter.Apply3DFit = True
0075 #########################
0076
0077 process.d0_phi_analyzer.BeamFitter.AsciiFileName = 'BeamFit_RunBased_Workflow.txt'
0078 process.d0_phi_analyzer.BeamFitter.AppendRunToFileName = False
0079 process.d0_phi_analyzer.BeamFitter.OutputFileName = 'BeamFit_RunBased_Workflow.root'
0080 #process.d0_phi_analyzer.BeamFitter.SaveNtuple = True
0081
0082 # fit as function of lumi sections
0083 process.d0_phi_analyzer.BSAnalyzerParameters.fitEveryNLumi = -1
0084 process.d0_phi_analyzer.BSAnalyzerParameters.resetEveryNLumi = -1