Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 23:31:52

0001 import FWCore.ParameterSet.Config as cms
0002 import FWCore.ParameterSet.VarParsing as VarParsing
0003 
0004 process = cms.Process("MCvertices")
0005 
0006 #prepare options
0007 
0008 options = VarParsing.VarParsing("analysis")
0009 
0010 options.register ('globalTag',
0011                   "DONOTEXIST",
0012                   VarParsing.VarParsing.multiplicity.singleton, # singleton or list
0013                   VarParsing.VarParsing.varType.string,          # string, int, or float
0014                   "GlobalTag")
0015 #options.globalTag = "DONOTEXIST"
0016 
0017 options.parseArguments()
0018 
0019 #
0020 process.options = cms.untracked.PSet(
0021     wantSummary = cms.untracked.bool(True),
0022     fileMode = cms.untracked.string("FULLMERGE")
0023     )
0024 
0025 process.load("FWCore.MessageService.MessageLogger_cfi")
0026 
0027 process.MessageLogger.cout.enable = cms.untracked.bool(True)
0028 process.MessageLogger.cout.threshold = cms.untracked.string("INFO")
0029 process.MessageLogger.cout.default = cms.untracked.PSet(
0030     limit = cms.untracked.int32(10000000)
0031     )
0032 process.MessageLogger.cout.FwkReport = cms.untracked.PSet(
0033     reportEvery = cms.untracked.int32(10000)
0034     )
0035 
0036 process.MessageLogger.cerr.enable = cms.untracked.bool(True)
0037 process.MessageLogger.cerr.threshold = cms.untracked.string("WARNING")
0038 process.MessageLogger.cerr.default = cms.untracked.PSet(
0039     limit = cms.untracked.int32(10000000)
0040     )
0041 process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
0042     reportEvery = cms.untracked.int32(100000)
0043     )
0044 
0045 #----Remove too verbose PrimaryVertexProducer
0046 
0047 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptive")
0048 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptiveNoBS")
0049 
0050 #----Remove too verbose BeamSpotOnlineProducer
0051 
0052 process.MessageLogger.suppressInfo.append("testBeamSpot")
0053 process.MessageLogger.suppressInfo.append("onlineBeamSpot")
0054 process.MessageLogger.suppressWarning.append("testBeamSpot")
0055 process.MessageLogger.suppressWarning.append("onlineBeamSpot")
0056 
0057 #----Remove too verbose TrackRefitter
0058 
0059 process.MessageLogger.suppressInfo.append("newTracksFromV0")
0060 process.MessageLogger.suppressInfo.append("newTracksFromOtobV0")
0061 
0062 
0063 #------------------------------------------------------------------
0064 
0065 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0066 
0067 process.source = cms.Source("PoolSource",
0068                     fileNames = cms.untracked.vstring(),
0069 #                    skipBadFiles = cms.untracked.bool(True),
0070                     inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0071                     )
0072 
0073 
0074 process.source.fileNames = cms.untracked.vstring(options.inputFiles)
0075 
0076 #process.load("Validation.RecoVertex.pvRecoSequence_cff")
0077 process.load("Validation.RecoVertex.pvSelectionSequence_cff")
0078 
0079 #process.load("Validation.RecoVertex.mcvertexweight_cfi")
0080 process.load("Validation.RecoVertex.mcverticesanalyzer_cfi")
0081 process.load("Validation.RecoVertex.mcvsrecoverticesanalyzer_cfi")
0082 
0083 #process.mcvertexweight45mm = process.mcvertexweight.clone()
0084 #process.mcvertexweight45mm.weighterConfig.finalSigma = cms.double(4.5)
0085 
0086 process.mcvsreconofakeverticesanalyzer = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("noFakeVertices"))
0087 process.mcvsrecogoodverticesanalyzer = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVertices"))
0088 
0089 process.mcvsrecoverticesanalyzervisible = process.mcvsrecoverticesanalyzer.clone(useVisibleVertices = cms.bool(True))
0090 process.mcvsreconofakeverticesanalyzervisible = process.mcvsreconofakeverticesanalyzer.clone(useVisibleVertices = cms.bool(True))
0091 process.mcvsrecogoodverticesanalyzervisible = process.mcvsrecogoodverticesanalyzer.clone(useVisibleVertices = cms.bool(True))
0092 #process.mcvsrecogoodverticesD0s5 = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVerticesD0s5"))
0093 #process.mcvsrecogoodverticesD0s51mm = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVerticesD0s51mm"))
0094 #process.mcvsrecogoodverticesDA100um = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVerticesDA100um"))
0095 #process.mcvsrecogoodverticesDA100umV7 = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVerticesDA100umV7"))
0096 #process.mcvsrecogoodverticesDA100umV8 = process.mcvsrecoverticesanalyzer.clone(pvCollection = cms.InputTag("goodVerticesDA100umV8"))
0097 
0098 
0099 #process.weightedmcverticesanalyzer = process.mcverticesanalyzer.clone(useWeight = cms.bool(True))
0100 #process.weightedmcvsrecoverticesanalyzer = process.mcvsrecoverticesanalyzer.clone(useWeight = cms.bool(True))
0101 #process.weightedmcvsreconofakeverticesanalyzer = process.mcvsreconofakeverticesanalyzer.clone(useWeight = cms.bool(True))
0102 #process.weightedmcvsrecogoodverticesanalyzer = process.mcvsrecogoodverticesanalyzer.clone(useWeight = cms.bool(True))
0103 #process.weightedmcvsrecogoodverticesD0s5 = process.mcvsrecogoodverticesD0s5.clone(useWeight = cms.bool(True))
0104 #process.weightedmcvsrecogoodverticesD0s51mm = process.mcvsrecogoodverticesD0s51mm.clone(useWeight = cms.bool(True))
0105 #process.weightedmcvsrecogoodverticesDA100um = process.mcvsrecogoodverticesDA100um.clone(useWeight = cms.bool(True))
0106 #process.weightedmcvsrecogoodverticesDA100umV7 = process.mcvsrecogoodverticesDA100umV7.clone(useWeight = cms.bool(True))
0107 #process.weightedmcvsrecogoodverticesDA100umV8 = process.mcvsrecogoodverticesDA100umV8.clone(useWeight = cms.bool(True))
0108 
0109 
0110 #process.weighted45mmmcverticesanalyzer = process.mcverticesanalyzer.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0111 #process.weighted45mmmcvsrecoverticesanalyzer = process.mcvsrecoverticesanalyzer.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0112 #process.weighted45mmmcvsreconofakeverticesanalyzer = process.mcvsreconofakeverticesanalyzer.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0113 #process.weighted45mmmcvsrecogoodverticesanalyzer = process.mcvsrecogoodverticesanalyzer.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0114 #process.weighted45mmmcvsrecogoodverticesD0s5 = process.mcvsrecogoodverticesD0s5.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0115 #process.weighted45mmmcvsrecogoodverticesD0s51mm = process.mcvsrecogoodverticesD0s51mm.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0116 #process.weighted45mmmcvsrecogoodverticesDA100um = process.mcvsrecogoodverticesDA100um.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0117 #process.weighted45mmmcvsrecogoodverticesDA100umV7 = process.mcvsrecogoodverticesDA100umV7.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0118 #process.weighted45mmmcvsrecogoodverticesDA100umV8 = process.mcvsrecogoodverticesDA100umV8.clone(useWeight = cms.bool(True), weightProduct=cms.InputTag("mcvertexweight45mm"))
0119 
0120 process.load("Validation.RecoVertex.anotherprimaryvertexanalyzer_cfi")
0121 
0122 process.primaryvertexanalyzer.vHistogramMakerPSet.runHisto=cms.untracked.bool(False)
0123 process.primaryvertexanalyzer.vHistogramMakerPSet.histoParameters = cms.untracked.PSet(
0124     nBinX = cms.untracked.uint32(2000), xMin=cms.untracked.double(-0.5), xMax=cms.untracked.double(0.5),
0125     nBinY = cms.untracked.uint32(2000), yMin=cms.untracked.double(-0.5), yMax=cms.untracked.double(0.5),
0126     nBinZ = cms.untracked.uint32(300), zMin=cms.untracked.double(-30.), zMax=cms.untracked.double(30.)
0127     )
0128 
0129 process.pvnofakevertices = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("noFakeVertices"))
0130 process.pvgoodvertices = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVertices"))
0131 #process.pvgoodverticesD0s5 = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVerticesD0s5"))
0132 #process.pvgoodverticesD0s51mm = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVerticesD0s51mm"))
0133 #process.pvgoodverticesDA100um = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVerticesDA100um"))
0134 #process.pvgoodverticesDA100umV7 = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVerticesDA100umV7"))
0135 #process.pvgoodverticesDA100umV8 = process.primaryvertexanalyzer.clone(pvCollection=cms.InputTag("goodVerticesDA100umV8"))
0136 
0137 
0138 process.p626 = cms.Path(
0139 #    process.seqPVReco +
0140     process.seqSimplePVSelection 
0141     + process.mcverticesanalyzer
0142     + process.mcvsrecoverticesanalyzer + process.mcvsreconofakeverticesanalyzer + process.mcvsrecogoodverticesanalyzer
0143     + process.mcvsrecoverticesanalyzervisible + process.mcvsreconofakeverticesanalyzervisible + process.mcvsrecogoodverticesanalyzervisible
0144 #    +  process.mcvsrecogoodverticesD0s5 +process.mcvsrecogoodverticesD0s51mm +
0145 #    process.mcvsrecogoodverticesDA100um +process.mcvsrecogoodverticesDA100umV7 +process.mcvsrecogoodverticesDA100umV8 +
0146     + process.pvnofakevertices + process.pvgoodvertices 
0147 #    + process.pvgoodverticesD0s5 + process.pvgoodverticesD0s51mm 
0148 #    + process.pvgoodverticesDA100um + process.pvgoodverticesDA100umV7 + process.pvgoodverticesDA100umV8
0149                         ) 
0150 
0151 #process.p520 = cms.Path(process.mcvertexweight+
0152 #                        process.seqPVReco + process.seqPVSelection +
0153 #                        process.weightedmcverticesanalyzer + process.weightedmcvsrecoverticesanalyzer +
0154 #                        process.weightedmcvsreconofakeverticesanalyzer + process.weightedmcvsrecogoodverticesanalyzer +
0155 #                        process.weightedmcvsrecogoodverticesD0s5 +process.weightedmcvsrecogoodverticesD0s51mm +
0156 #                        process.weightedmcvsrecogoodverticesDA100um + process.weightedmcvsrecogoodverticesDA100umV7 + process.weightedmcvsrecogoodverticesDA100umV8
0157 #                        )
0158 
0159 #process.p450 = cms.Path(process.mcvertexweight45mm+
0160 #                        process.seqPVReco + process.seqPVSelection +
0161 #                        process.weighted45mmmcverticesanalyzer + process.weighted45mmmcvsrecoverticesanalyzer +
0162 #                        process.weighted45mmmcvsreconofakeverticesanalyzer + process.weighted45mmmcvsrecogoodverticesanalyzer +
0163 #                        process.weighted45mmmcvsrecogoodverticesD0s5 +process.weighted45mmmcvsrecogoodverticesD0s51mm +
0164 #                        process.weighted45mmmcvsrecogoodverticesDA100um + process.weighted45mmmcvsrecogoodverticesDA100umV7 + process.weighted45mmmcvsrecogoodverticesDA100umV8
0165 #                        )
0166 
0167 #----GlobalTag ------------------------
0168 
0169 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0170 from Configuration.AlCa.GlobalTag import GlobalTag
0171 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0172 
0173 
0174 process.TFileService = cms.Service('TFileService',
0175                                    fileName = cms.string('mcvertices.root')
0176                                    )
0177