From 30f192c50747c4fb149d540a14a1a37a9946e2e3 Mon Sep 17 00:00:00 2001 From: Riccardo Manzoni Date: Wed, 10 Jun 2015 18:32:14 +0200 Subject: [PATCH 1/3] modified _cfg to rerun track and PV fitter and BS fitter --- .../BeamFit_LumiBased_NewAlignWorkflow.py | 145 +++++++----------- 1 file changed, 59 insertions(+), 86 deletions(-) diff --git a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py index a1287ff07a913..44142928620d5 100644 --- a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py +++ b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py @@ -1,112 +1,85 @@ import FWCore.ParameterSet.Config as cms process = cms.Process("BSworkflow") -# initialize MessageLogger -process.load("FWCore.MessageLogger.MessageLogger_cfi") - -process.load("RecoVertex.BeamSpotProducer.d0_phi_analyzer_cff") process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring( -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0171/3A22D08D-456A-DF11-961F-001A92811734.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0171/083B6802-236C-DF11-8AC6-0026189437FE.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0170/E6D0589B-136A-DF11-9E90-002618943982.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0170/D098488B-276A-DF11-8069-003048678AF4.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/F466DD94-BF69-DF11-B9B8-00261894390A.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/F2A9245B-BC69-DF11-9C3F-0018F3D096E4.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/DED5E502-B969-DF11-AE8F-002618943964.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/563257AF-C169-DF11-865F-002618943907.root', -'/store/data/Commissioning10/MinimumBias/ALCARECO/TkAlMinBias-May27thReReco_v1/0166/3056B847-C069-DF11-B0B0-003048679010.root' - + "/store/express/Run2015A/StreamExpress/ALCARECO/TkAlMinBias-Express-v1/000/246/959/00000/14174DF2-490A-E511-9862-02163E0143E9.root", ) ) +process.load("FWCore.MessageLogger.MessageLogger_cfi") process.MessageLogger.cerr.FwkReport = cms.untracked.PSet( - reportEvery = cms.untracked.int32(1000000), + reportEvery = cms.untracked.int32(10000), ) - -#process.source = cms.Source('PoolSource', -# debugVerbosity = cms.untracked.uint32(0), -# debugFlag = cms.untracked.bool(False) -# ) +process.MessageLogger.debugModules = ['BeamSpotAnalyzer'] process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(100) #1500 + input = cms.untracked.int32(-1) ) process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) - ) - -# this is for filtering on L1 technical trigger bit -process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff') -process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi') -process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True) -process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND ( 40 OR 41 )') -## -process.load("Configuration.StandardSequences.MagneticField_cff") -process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") -process.GlobalTag.globaltag = 'GR_R_36X_V11::All' #'GR_R_35X_V8::All' -process.load("Configuration.StandardSequences.Geometry_cff") +) +process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi") +process.load("Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff") +process.load('Configuration.Geometry.GeometryRecoDB_cff') +process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff") +process.GlobalTag.globaltag = 'GR_E_V48' -########## RE-FIT TRACKS +## Track refit process.load("RecoTracker.TrackProducer.TrackRefitters_cff") + +# remove the following lines if you run on RECO files process.TrackRefitter.src = 'ALCARECOTkAlMinBias' +process.TrackRefitter.NavigationSchool = '' -## reco PV -process.load("Configuration.StandardSequences.Reconstruction_cff") -process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi") -process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi") -process.offlinePrimaryVertices.TrackLabel = cms.InputTag("TrackRefitter") - -#### remove beam scraping events -process.noScraping= cms.EDFilter("FilterOutScraping", - applyfilter = cms.untracked.bool(True), - debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info - numtrack = cms.untracked.uint32(10), - thresh = cms.untracked.double(0.20) -) -process.p = cms.Path(process.hltLevel1GTSeed + - process.offlineBeamSpot + - process.TrackRefitter + - process.offlinePrimaryVertices+ -# process.noScraping + - process.d0_phi_analyzer) +## PV refit +process.load("TrackingTools.TransientTrack.TransientTrackBuilder_cfi") + +from RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi import offlinePrimaryVertices +process.offlinePrimaryVerticesFromRefittedTrks = offlinePrimaryVertices.clone() +process.offlinePrimaryVerticesFromRefittedTrks.TrackLabel = cms.InputTag("TrackRefitter") +process.offlinePrimaryVerticesFromRefittedTrks.vertexCollections.maxDistanceToBeam = 1 +process.offlinePrimaryVerticesFromRefittedTrks.TkFilterParameters.maxNormalizedChi2 = 20 +process.offlinePrimaryVerticesFromRefittedTrks.TkFilterParameters.minSiliconLayersWithHits = 5 +process.offlinePrimaryVerticesFromRefittedTrks.TkFilterParameters.maxD0Significance = 5.0 +process.offlinePrimaryVerticesFromRefittedTrks.TkFilterParameters.minPixelLayersWithHits = 2 -process.MessageLogger.debugModules = ['BeamSpotAnalyzer'] -################### Primary Vertex -process.offlinePrimaryVertices.PVSelParameters.maxDistanceToBeam = 2 -process.offlinePrimaryVertices.TkFilterParameters.maxNormalizedChi2 = 20 -process.offlinePrimaryVertices.TkFilterParameters.minSiliconLayersWithHits = 5 -process.offlinePrimaryVertices.TkFilterParameters.maxD0Significance = 100 -process.offlinePrimaryVertices.TkFilterParameters.minPixelLayersWithHits = 1 -process.offlinePrimaryVertices.TkClusParameters.TkGapClusParameters.zSeparation = 1 - - -####################### -process.d0_phi_analyzer.BeamFitter.TrackCollection = 'TrackRefitter' -process.d0_phi_analyzer.BeamFitter.MinimumTotalLayers = 6 -process.d0_phi_analyzer.BeamFitter.MinimumPixelLayers = -1 -process.d0_phi_analyzer.BeamFitter.MaximumNormChi2 = 10 -process.d0_phi_analyzer.BeamFitter.MinimumInputTracks = 50 -process.d0_phi_analyzer.BeamFitter.MinimumPt = 1.0 -process.d0_phi_analyzer.BeamFitter.MaximumImpactParameter = 1.0 -process.d0_phi_analyzer.BeamFitter.TrackAlgorithm = cms.untracked.vstring() -#process.d0_phi_analyzer.BeamFitter.Debug = True - -process.d0_phi_analyzer.PVFitter.Apply3DFit = True -process.d0_phi_analyzer.PVFitter.minNrVerticesForFit = 10 -######################### - -process.d0_phi_analyzer.BeamFitter.AsciiFileName = 'BeamFit_LumiBased_NewAlignWorkflow.txt' -process.d0_phi_analyzer.BeamFitter.AppendRunToFileName = False -process.d0_phi_analyzer.BeamFitter.OutputFileName = 'BeamFit_LumiBased_Workflow.root' -#process.d0_phi_analyzer.BeamFitter.SaveNtuple = True -#process.d0_phi_analyzer.BeamFitter.SavePVVertices = True - -# fit as function of lumi sections -process.d0_phi_analyzer.BSAnalyzerParameters.fitEveryNLumi = 1 + +## BeamSpot fit +process.load("RecoVertex.BeamSpotProducer.d0_phi_analyzer_cff") + +process.d0_phi_analyzer.BeamFitter.WriteAscii = True +process.d0_phi_analyzer.BeamFitter.AsciiFileName = 'BeamFit_LumiBased_NewAlignWorkflow_alcareco.txt' +process.d0_phi_analyzer.BeamFitter.AppendRunToFileName = False +process.d0_phi_analyzer.BeamFitter.InputBeamWidth = -1 +process.d0_phi_analyzer.BeamFitter.MaximumImpactParameter = 1.0 +process.d0_phi_analyzer.BeamFitter.MaximumNormChi2 = 10 +process.d0_phi_analyzer.BeamFitter.MinimumInputTracks = 50 +process.d0_phi_analyzer.BeamFitter.MinimumPixelLayers = -1 +process.d0_phi_analyzer.BeamFitter.MinimumPt = 1.0 +process.d0_phi_analyzer.BeamFitter.MinimumTotalLayers = 6 +process.d0_phi_analyzer.BeamFitter.OutputFileName = 'BeamFit_LumiBased_Workflow_alcareco.root' +process.d0_phi_analyzer.BeamFitter.TrackAlgorithm = cms.untracked.vstring() +process.d0_phi_analyzer.BeamFitter.TrackCollection = 'TrackRefitter' +process.d0_phi_analyzer.BeamFitter.SaveFitResults = True +process.d0_phi_analyzer.BeamFitter.SaveNtuple = True +process.d0_phi_analyzer.BeamFitter.SavePVVertices = True + +process.d0_phi_analyzer.PVFitter.Apply3DFit = True +process.d0_phi_analyzer.PVFitter.minNrVerticesForFit = 10 +process.d0_phi_analyzer.PVFitter.nSigmaCut = 50.0 +process.d0_phi_analyzer.PVFitter.VertexCollection = 'offlinePrimaryVerticesFromRefittedTrks' + +process.d0_phi_analyzer.BSAnalyzerParameters.fitEveryNLumi = 1 process.d0_phi_analyzer.BSAnalyzerParameters.resetEveryNLumi = 1 + + +process.p = cms.Path(process.offlineBeamSpot + + process.TrackRefitter + + process.offlinePrimaryVerticesFromRefittedTrks + + process.d0_phi_analyzer) From c8188465500338585afa84310e06766fa35e8479 Mon Sep 17 00:00:00 2001 From: Riccardo Manzoni Date: Thu, 11 Jun 2015 11:09:32 +0200 Subject: [PATCH 2/3] modified magnetic field and global tag conditions, to stick to the 75X defaults --- .../test/BeamFit_LumiBased_NewAlignWorkflow.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py index 44142928620d5..a8570197fb649 100644 --- a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py +++ b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py @@ -23,9 +23,11 @@ ) process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi") -process.load("Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff") +process.load("Configuration.StandardSequences.MagneticField_cff") process.load('Configuration.Geometry.GeometryRecoDB_cff') -process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff") +process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") +# this GT is for the Express, to be consistent with the file above +# in general this GT should be for the ReReco process.GlobalTag.globaltag = 'GR_E_V48' ## Track refit From 81dc5af2837e2ff6d5dccf1859d8cb3eedc7ce44 Mon Sep 17 00:00:00 2001 From: Riccardo Manzoni Date: Thu, 11 Jun 2015 12:40:52 +0200 Subject: [PATCH 3/3] turn off set ntuple saving to default false --- .../BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py index a8570197fb649..43fd668763e79 100644 --- a/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py +++ b/RecoVertex/BeamSpotProducer/test/BeamFit_LumiBased_NewAlignWorkflow.py @@ -69,7 +69,7 @@ process.d0_phi_analyzer.BeamFitter.TrackAlgorithm = cms.untracked.vstring() process.d0_phi_analyzer.BeamFitter.TrackCollection = 'TrackRefitter' process.d0_phi_analyzer.BeamFitter.SaveFitResults = True -process.d0_phi_analyzer.BeamFitter.SaveNtuple = True +process.d0_phi_analyzer.BeamFitter.SaveNtuple = False process.d0_phi_analyzer.BeamFitter.SavePVVertices = True process.d0_phi_analyzer.PVFitter.Apply3DFit = True