BeamSpotProducer/test/BeamFit_LumiBased_template.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     {DATAFILE}
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 )')
0037 ##
0038
0039 ## reco PV
0040 process.load("Configuration.StandardSequences.MagneticField_cff")
0041 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0042 process.GlobalTag.globaltag = 'GR10_P_V5::All'
0043
0044 process.load("Configuration.StandardSequences.Reconstruction_cff")
0045 process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
0046 process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
0047 process.offlinePrimaryVertices.TrackLabel = cms.InputTag("ALCARECOTkAlMinBias")
0048
0049 #### remove beam scraping events
0050 process.noScraping= cms.EDFilter("FilterOutScraping",
0051                                  applyfilter = cms.untracked.bool(True),
0052                                  debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0053                                  numtrack = cms.untracked.uint32(10),
0054                                  thresh = cms.untracked.double(0.20)
0055 )
0056
0057 process.p = cms.Path(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.minSiliconHits = 6
0069 process.offlinePrimaryVertices.TkFilterParameters.maxD0Significance = 100
0070 process.offlinePrimaryVertices.TkFilterParameters.minPixelHits = 1
0071 process.offlinePrimaryVertices.TkClusParameters.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 = 2
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.InputBeamWidth = -1 # 0.0400
0084 #process.d0_phi_analyzer.BeamFitter.Debug = True
0085
0086 process.d0_phi_analyzer.PVFitter.Apply3DFit = True
0087 process.d0_phi_analyzer.PVFitter.minNrVerticesForFit = 10
0088 #########################
0089
0090 process.d0_phi_analyzer.BeamFitter.AsciiFileName = '{OUTPUT_FILENAME}'
0091 #process.d0_phi_analyzer.BeamFitter.AppendRunToFileName = False
0092 process.d0_phi_analyzer.BeamFitter.OutputFileName = 'BeamFit_RunBased_Workflow.root'
0093 #process.d0_phi_analyzer.BeamFitter.SaveNtuple = 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