Alignment/CommonAlignmentProducer/test/TrackerAPForPCL/alignment_mille.py - master - githubmirror/cmssw

Alignment/CommonAlignmentProducer/test/TrackerAPForPCL/alignment_mille.py
229 строк9 KB
Christopher Jones
Update MessageLogger_cfi to new syntax
07 дек 2020, 21:37
07 дек 2020, 21:37eb0725a
#-- Common selection based on CRUZET 2015 Setup mp1553

import FWCore.ParameterSet.Config as cms
#from Configuration.AlCa.autoCond import *

process = cms.Process("Alignment")



process.options = cms.untracked.PSet(
    Rethrow = cms.untracked.vstring("ProductNotFound"), # do not accept this exception
    wantSummary = cms.untracked.bool(True)
    )

# initialize  MessageLogger
process.load("FWCore.MessageService.MessageLogger_cfi")

process.MessageLogger.files.alignment = cms.untracked.PSet(
    DEBUG = cms.untracked.PSet(
        limit = cms.untracked.int32(-1)
        ),
    INFO = cms.untracked.PSet(
        limit = cms.untracked.int32(-1)
        ),
    WARNING = cms.untracked.PSet(
        limit = cms.untracked.int32(10)
        ),
    ERROR = cms.untracked.PSet(
        limit = cms.untracked.int32(-1)
        ),
    Alignment = cms.untracked.PSet(
        limit = cms.untracked.int32(-1),
        ),
    enableStatistics = cms.untracked.bool(True)
    )

   
process.MessageLogger.cerr.enable = cms.untracked.bool(False)
process.options   = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )


#-- Magnetic field
process.load('Configuration.StandardSequences.MagneticField_cff')
#process.load("Configuration/StandardSequences/MagneticField_38T_cff") ## FOR 3.8T
#process.load("Configuration.StandardSequences.MagneticField_0T_cff")  ## FOR 0T

#-- Load geometry
process.load("Configuration.Geometry.GeometryRecoDB_cff")

#-- Global Tag
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.connect = "frontier://FrontierProd/CMS_CONDITIONS"
process.GlobalTag.globaltag = "GR_P_V56"

#-- initialize beam spot
process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")

#-- Set APEs to ZERO
import CalibTracker.Configuration.Common.PoolDBESSource_cfi
process.conditionsInTrackerAlignmentErrorRcd = CalibTracker.Configuration.Common.PoolDBESSource_cfi.poolDBESSource.clone(
    connect = cms.string('frontier://FrontierProd/CMS_CONDITIONS'),
    toGet = cms.VPSet(
        cms.PSet(
            record = cms.string('TrackerAlignmentErrorExtendedRcd'),
            tag = cms.string('TrackerIdealGeometryErrorsExtended210_mc')
            )
        )
    )
process.prefer_conditionsInTrackerAlignmentErrorRcd = cms.ESPrefer("PoolDBESSource", "conditionsInTrackerAlignmentErrorRcd")

#-- Load Local Alignment Object
#process.GlobalTag.toGet = cms.VPSet(
#    cms.PSet(record = cms.string("TrackerAlignmentRcd"),
#             tag = cms.string("testTag"),
#             connect = cms.untracked.string("sqlite_file:/afs/cern.ch/work/c/chmartin/private/Alignment/pp3.8T_2015_Alignment/Local_DB/TkAlignment.db")
#             )
#    )
        
#-- AlignmentTrackSelector
import Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi
process.HighPuritySelector = Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi.AlignmentTrackSelector.clone(
    applyBasicCuts = True,
    #filter = True,
    src = 'ALCARECOTkAlMinBias',
    trackQualities = ["highPurity"],
    pMin = 4.,
    )

process.load("Alignment.CommonAlignmentProducer.AlignmentTrackSelector_cfi")
process.AlignmentTrackSelector.src = 'HitFilteredTracks' # adjust to input file
process.AlignmentTrackSelector.applyBasicCuts = True
process.AlignmentTrackSelector.pMin = 8.
process.AlignmentTrackSelector.ptMin = 1.0
process.AlignmentTrackSelector.etaMin = -999.
process.AlignmentTrackSelector.etaMax = 999.
process.AlignmentTrackSelector.nHitMin = 8
process.AlignmentTrackSelector.nHitMin2D = 2
process.AlignmentTrackSelector.chi2nMax = 9999.
process.AlignmentTrackSelector.applyMultiplicityFilter = False 
process.AlignmentTrackSelector.maxMultiplicity = 1
process.AlignmentTrackSelector.applyNHighestPt = False
process.AlignmentTrackSelector.nHighestPt = 1
process.AlignmentTrackSelector.seedOnlyFrom = 0
process.AlignmentTrackSelector.applyIsolationCut = False
process.AlignmentTrackSelector.minHitIsolation = 0.8
process.AlignmentTrackSelector.applyChargeCheck = False
process.AlignmentTrackSelector.minHitChargeStrip = 30.
#Special option for PCL
process.AlignmentTrackSelector.minHitsPerSubDet.inPIXEL = 2


#-- new track hit filter
# TrackerTrackHitFilter takes as input the tracks/trajectories coming out from TrackRefitter1
process.load("RecoTracker.FinalTrackSelectors.TrackerTrackHitFilter_cff")
process.TrackerTrackHitFilter.src = 'TrackRefitter1'
process.TrackerTrackHitFilter.useTrajectories= True  # this is needed only if you require some selections; but it will work even if you don't ask for them
process.TrackerTrackHitFilter.minimumHits = 8
process.TrackerTrackHitFilter.replaceWithInactiveHits = True
process.TrackerTrackHitFilter.rejectBadStoNHits = True
process.TrackerTrackHitFilter.commands = cms.vstring("keep PXB","keep PXE","keep TIB","keep TID","keep TOB","keep TEC")#,"drop TID stereo","drop TEC stereo")
process.TrackerTrackHitFilter.stripAllInvalidHits = False
process.TrackerTrackHitFilter.StoNcommands = cms.vstring("ALL 18.0")
process.TrackerTrackHitFilter.rejectLowAngleHits = True
process.TrackerTrackHitFilter.TrackAngleCut = 0.17# in rads, starting from the module surface; .35 for cosmcics ok, .17 for collision tracks
process.TrackerTrackHitFilter.usePixelQualityFlag = True #False

#-- TrackFitter
import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cff
process.HitFilteredTracks = RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cff.ctfWithMaterialTracks.clone(
       	src = 'TrackerTrackHitFilter',
    	TrajectoryInEvent = True,
    	TTRHBuilder = 'WithAngleAndTemplate', #should already be default
        NavigationSchool = cms.string('')
)
    
#-- Alignment producer
process.load("Alignment.CommonAlignmentProducer.AlignmentProducerAsAnalyzer_cff")
process.AlignmentProducer.ParameterBuilder.Selector = cms.PSet(
    alignParams = cms.vstring(
        "PixelHalfBarrels,111111",
        "PXECHalfCylinders,111111",
        )
    )

process.AlignmentProducer.doMisalignmentScenario = False #True


process.AlignmentProducer.checkDbAlignmentValidity = False
process.AlignmentProducer.applyDbAlignment = True
process.AlignmentProducer.tjTkAssociationMapTag = 'TrackRefitter2'

process.AlignmentProducer.algoConfig = process.MillePedeAlignmentAlgorithm
process.AlignmentProducer.algoConfig.mode = 'full'
process.AlignmentProducer.algoConfig.mergeBinaryFiles = cms.vstring()
process.AlignmentProducer.algoConfig.binaryFile = 'milleBinary_40.dat'
process.AlignmentProducer.algoConfig.TrajectoryFactory = cms.PSet(
      #process.BrokenLinesBzeroTrajectoryFactory
      process.BrokenLinesTrajectoryFactory 
      )
process.AlignmentProducer.algoConfig.pedeSteerer.pedeCommand = 'pede'
process.AlignmentProducer.algoConfig.pedeSteerer.method = 'inversion  5  0.8'
process.AlignmentProducer.algoConfig.pedeSteerer.options = cms.vstring(
    #'regularisation 1.0 0.05', # non-stated pre-sigma 50 mrad or 500 mum 
     'entries 500',
     'chisqcut  30.0  4.5', #,
     'threads 1 1' #, 
     #'outlierdownweighting 3','dwfractioncut 0.1'
     #'outlierdownweighting 5','dwfractioncut 0.2'
    )

#-- TrackRefitter
process.load("RecoTracker.TrackProducer.TrackRefitters_cff")
process.TrackRefitter1 = RecoTracker.TrackProducer.TrackRefitter_cfi.TrackRefitter.clone(
    src ='HighPuritySelector',#'ALCARECOTkAlMinBias',
    NavigationSchool = cms.string(''),
    TrajectoryInEvent = True,
    TTRHBuilder = "WithAngleAndTemplate" #default
    )

process.TrackRefitter2 = process.TrackRefitter1.clone(
    src = 'AlignmentTrackSelector',
#    TTRHBuilder = 'WithTrackAngle'
    )
    
process.source = cms.Source("PoolSource",
		skipEvents = cms.untracked.uint32(0),
		fileNames = cms.untracked.vstring(
			'/store/express/Run2015B/StreamExpress/ALCARECO/TkAlMinBias-Express-v1/000/251/718/00000/426E4511-292A-E511-9FB8-02163E0133E0.root'
    ),
#    lumisToSkip = cms.untracked.VLuminosityBlockRange('251718:1-251718:20'),
)


#process.maxEvents = cms.untracked.PSet(
#    input = cms.untracked.int32(20000)
#) 

#process.dump = cms.EDAnalyzer("EventContentAnalyzer")                                                         
#process.p  = cms.Path(process.dump) 


process.p = cms.Path(process.HighPuritySelector
		     *process.offlineBeamSpot
                     *process.TrackRefitter1
                     *process.TrackerTrackHitFilter
                     *process.HitFilteredTracks
                     *process.AlignmentTrackSelector
                     *process.TrackRefitter2
                     *process.AlignmentProducer
                     )
process.AlignmentProducer.saveToDB = True
from CondCore.DBCommon.CondDBSetup_cfi import *
process.PoolDBOutputService = cms.Service(
        "PoolDBOutputService",
            CondDBSetup,
            connect = cms.string('sqlite_file:TkAlignment.db'),
            toPut = cms.VPSet(cms.PSet(
                record = cms.string('TrackerAlignmentRcd'),
                tag = cms.string('testTag')
            )#,
                              #     #                  cms.PSet(
                              #     #  record = cms.string('TrackerAlignmentErrorRcd'),
                              #     #  tag = cms.string('testTagAPE') # needed is saveApeToDB = True
                              #     #)
                                                     )
            )                     
                     
# MPS needs next line as placeholder for pede _cfg.py:
#MILLEPEDEBLOCK