RecoTrack/test/cosmictrackingPerformanceValidation_cfg.py

0001 #import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TkVal")
0004 process.load("FWCore.MessageService.MessageLogger_cfi")
0005
0006 ### standard includes
0007 process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
0008 process.load("Configuration.StandardSequences.RawToDigi_cff")
0009 process.load("Configuration.EventContent.EventContentCosmics_cff")
0010 process.load("Configuration.StandardSequences.ReconstructionCosmics_cff")
0011 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0012 process.load("SimGeneral.MixingModule.mixNoPU_cfi")
0013 process.load("RecoTracker.TrackProducer.TrackRefitterP5_cfi")
0014
0015
0016 ### conditions
0017 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0018 #process.GlobalTag.globaltag = 'STARTUP_V1::All'
0019 process.GlobalTag.globaltag = 'GLOBALTAG::All'
0020
0021
0022 process.maxEvents = cms.untracked.PSet(
0023     input = cms.untracked.int32(NEVENT)
0024 )
0025 process.source = source
0026
0027 ### validation-specific includes
0028 #process.load("SimTracker.TrackAssociatorProducers.trackAssociatorByChi2_cfi")
0029 #process.load("SimTracker.TrackAssociatorProducers.trackAssociatorByHits_cfi")
0030 process.load("SimTracker.TrackAssociatorProducers.quickTrackAssociatorByHits_cfi")
0031 process.load("Validation.RecoTrack.cosmiccuts_cff")
0032 process.load("Validation.RecoTrack.MultiTrackValidator_cff")
0033 process.load("SimGeneral.TrackingAnalysis.trackingParticles_cfi")
0034 process.load("DQMServices.Components.EDMtoMEConverter_cff")
0035
0036 process.load("Validation.Configuration.postValidation_cff")
0037
0038 process.quickTrackAssociatorByHits.SimToRecoDenominator = cms.string('reco')
0039
0040 ## configuration of the cuts on CTF tracks
0041 process.cutsRecoCTFTracksP5.minHit = cms.int32(NHITRECOTRKMIN)
0042 process.cutsRecoCTFTracksP5.ptMin = cms.double(PTRECOTRKMIN)
0043 process.cutsRecoCTFTracksP5.maxRapidity = cms.double(ETARECOTRKMAX)
0044 process.cutsRecoCTFTracksP5.minRapidity = cms.double(ETARECOTRKMIN)
0045 process.cutsRecoCTFTracksP5.lip = cms.double(LIPRECOTRKMAX)
0046 process.cutsRecoCTFTracksP5.tip = cms.double(TIPRECOTRKMAX)
0047 process.cutsRecoCTFTracksP5.maxChi2 = cms.double(CHISQRECOTRKMAX)
0048
0049 ## configuration of the cuts on CosmicTF tracks
0050 process.cutsRecoCosmicTFTracksP5.minHit = cms.int32(NHITRECOTRKMIN)
0051 process.cutsRecoCosmicTFTracksP5.ptMin = cms.double(PTRECOTRKMIN)
0052 process.cutsRecoCosmicTFTracksP5.maxRapidity = cms.double(ETARECOTRKMAX)
0053 process.cutsRecoCosmicTFTracksP5.minRapidity = cms.double(ETARECOTRKMIN)
0054 process.cutsRecoCosmicTFTracksP5.lip = cms.double(LIPRECOTRKMAX)
0055 process.cutsRecoCosmicTFTracksP5.tip = cms.double(TIPRECOTRKMAX)
0056 process.cutsRecoCosmicTFTracksP5.maxChi2 = cms.double(CHISQRECOTRKMAX)
0057
0058 ## configuration of the cuts on RS tracks
0059 process.cutsRecoRSTracksP5.minHit = cms.int32(NHITRECOTRKMIN)
0060 process.cutsRecoRSTracksP5.ptMin = cms.double(PTRECOTRKMIN)
0061 process.cutsRecoRSTracksP5.maxRapidity = cms.double(ETARECOTRKMAX)
0062 process.cutsRecoRSTracksP5.minRapidity = cms.double(ETARECOTRKMIN)
0063 process.cutsRecoRSTracksP5.lip = cms.double(LIPRECOTRKMAX)
0064 process.cutsRecoRSTracksP5.tip = cms.double(TIPRECOTRKMAX)
0065 process.cutsRecoRSTracksP5.maxChi2 = cms.double(CHISQRECOTRKMAX)
0066
0067 ### configuration MultiTrackValidator ###
0068 process.multiTrackValidator.outputFile = 'val.SAMPLE.root'
0069 process.multiTrackValidator.UseAssociators = cms.bool(True)
0070 process.multiTrackValidator.associators = ['quickTrackAssociatorByHits']
0071 process.multiTrackValidator.parametersDefiner = cms.string('CosmicParametersDefinerForTP')
0072 process.multiTrackValidator.label = ['cutsRecoCTFTracksP5', 'cutsRecoCosmicTFTracksP5', 'cutsRecoRSTracksP5']
0073 process.multiTrackValidator.useLogPt=cms.untracked.bool(False)
0074 process.multiTrackValidator.min = cms.double(MINETA)
0075 process.multiTrackValidator.max = cms.double(MAXETA)
0076 process.multiTrackValidator.nint = cms.int32(NINTETA)
0077 process.multiTrackValidator.minpT = cms.double(MINPT)
0078 process.multiTrackValidator.maxpT = cms.double(MAXPT)
0079 process.multiTrackValidator.nintpT = cms.int32(NINTPT)
0080 process.multiTrackValidator.minDxy = cms.double(MINDXY)
0081 process.multiTrackValidator.maxDxy = cms.double(MAXDXY)
0082 process.multiTrackValidator.nintDxy = cms.int32(NINTDXY)
0083 process.multiTrackValidator.minDz = cms.double(MINDZ)
0084 process.multiTrackValidator.maxDz = cms.double(MAXDZ)
0085 process.multiTrackValidator.nintDz = cms.int32(NINTDZ)
0086 process.multiTrackValidator.minPhi = cms.double(MINPHI)
0087 process.multiTrackValidator.maxPhi = cms.double(MAXPHI)
0088 process.multiTrackValidator.nintPhi = cms.int32(NINTPHI)
0089 process.multiTrackValidator.minVertpos = cms.double(MINVERTPOS)
0090 process.multiTrackValidator.maxVertpos = cms.double(MAXVERTPOS)
0091 process.multiTrackValidator.nintVertpos = cms.int32(NINTVERTPOS)
0092 process.multiTrackValidator.minZpos = cms.double(MINZPOS)
0093 process.multiTrackValidator.maxZpos = cms.double(MAXZPOS)
0094 process.multiTrackValidator.nintZpos = cms.int32(NINTZPOS)
0095 process.multiTrackValidator.dxyRes_rangeMin = cms.double(DXYRESMIN)
0096 process.multiTrackValidator.dxyRes_rangeMax = cms.double(DXYRESMAX)
0097 process.multiTrackValidator.dxyRes_nbin = cms.int32(100)
0098 process.multiTrackValidator.dzRes_rangeMin = cms.double(DZRESMIN)
0099 process.multiTrackValidator.dzRes_rangeMax = cms.double(DZRESMAX)
0100 process.multiTrackValidator.dzRes_nbin = cms.int32(100)
0101 process.multiTrackValidator.phiRes_rangeMin = cms.double(PHIRESMIN)
0102 process.multiTrackValidator.phiRes_rangeMax = cms.double(PHIRESMAX)
0103 process.multiTrackValidator.phiRes_nbin = cms.int32(100)
0104 process.multiTrackValidator.ptRes_rangeMin = cms.double(PTRESMIN)
0105 process.multiTrackValidator.ptRes_rangeMax = cms.double(PTRESMAX)
0106 process.multiTrackValidator.ptRes_nbin = cms.int32(100)
0107 process.multiTrackValidator.cotThetaRes_rangeMin = cms.double(COTTHETARESMIN)
0108 process.multiTrackValidator.cotThetaRes_rangeMax = cms.double(COTTHETARESMAX)
0109 process.multiTrackValidator.cotThetaRes_nbin = cms.int32(100)
0110 process.multiTrackValidator.ptMinTP = cms.double(PTTPMIN)
0111 process.multiTrackValidator.lipTP = cms.double(LIPTPMAX)
0112 process.multiTrackValidator.tipTP = cms.double(TIPTPMAX)
0113 process.multiTrackValidator.minRapidityTP = cms.double(ETATPMIN)
0114 process.multiTrackValidator.maxRapidityTP = cms.double(ETATPMAX)
0115 process.multiTrackValidator.minHitTP = cms.int32(NHITTPMIN)
0116
0117 process.options = cms.untracked.PSet(
0118     wantSummary = cms.untracked.bool(True)
0119 )
0120
0121 process.re_tracking = cms.Sequence(process.siPixelRecHits*process.siStripMatchedRecHits*
0122                                    process.tracksP5*
0123                                    process.cutsRecoCTFTracksP5*process.cutsRecoCosmicTFTracksP5*process.cutsRecoRSTracksP5*
0124                                    process.quickTrackAssociatorByHits*
0125                                    process.multiTrackValidator
0126                                    )
0127
0128 process.re_tracking_and_TP = cms.Sequence(process.trackingParticles*
0129                                    process.siPixelRecHits*process.siStripMatchedRecHits*
0130                                    process.tracksP5*
0131                                    process.cutsRecoCTFTracksP5*process.cutsRecoCosmicTFTracksP5*process.cutsRecoRSTracksP5*
0132                                    process.quickTrackAssociatorByHits*
0133                                    process.multiTrackValidator
0134                                    )
0135
0136 process.only_validation = cms.Sequence(process.cutsRecoCTFTracksP5*process.cutsRecoCosmicTFTracksP5*process.cutsRecoRSTracksP5*
0137                                        process.quickTrackAssociatorByHits*
0138                                        process.multiTrackValidator)
0139
0140 process.only_validation_and_TP = cms.Sequence(process.trackingParticles*
0141                                               process.cutsRecoCTFTracksP5*process.cutsRecoCosmicTFTracksP5*process.cutsRecoRSTracksP5*
0142                                               process.quickTrackAssociatorByHits*
0143                                               process.multiTrackValidator)
0144
0145 ### customized versoin of the OutputModule
0146 ### it save the mininal information which is necessary to perform tracking validation (tracks, tracking particles,
0147 ### digiSimLink,etc..)
0148
0149 process.customEventContent = cms.PSet(
0150      outputCommands = cms.untracked.vstring('drop *')
0151  )
0152
0153 process.customEventContent.outputCommands.extend(process.RecoTrackerRECO.outputCommands)
0154 process.customEventContent.outputCommands.extend(process.BeamSpotRECO.outputCommands)
0155 process.customEventContent.outputCommands.extend(process.SimGeneralFEVTDEBUG.outputCommands)
0156 process.customEventContent.outputCommands.extend(process.RecoLocalTrackerRECO.outputCommands)
0157 process.customEventContent.outputCommands.append('keep *_simSiStripDigis_*_*')
0158 process.customEventContent.outputCommands.append('keep *_simSiPixelDigis_*_*')
0159 process.customEventContent.outputCommands.append('drop SiStripDigiedmDetSetVector_simSiStripDigis_*_*')
0160 process.customEventContent.outputCommands.append('drop PixelDigiedmDetSetVector_simSiPixelDigis_*_*')
0161
0162
0163
0164 process.OUTPUT = cms.OutputModule("PoolOutputModule",
0165                                   process.customEventContent,
0166                                   fileName = cms.untracked.string('output.SAMPLE.root')
0167                                   )
0168
0169
0170 ### final path and endPath
0171 process.p = cms.Path(process.SEQUENCE)
0172 #process.outpath = cms.EndPath(process.OUTPUT)
0173
0174