PFTracking/test/fastSimWithDisplacedVertexFinder_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003
0004 process = cms.Process("TEST")
0005
0006
0007 process.maxEvents = cms.untracked.PSet(
0008     input = cms.untracked.int32(20000)
0009 )
0010
0011 #generation
0012 #process.load("RecoParticleFlow.Configuration.source_singleTau_cfi")
0013 #process.load("RecoParticleFlow.PFTracking.source_jetGun_NuclearTest_cfi")
0014 process.load("RecoParticleFlow.PFTracking.source_particleGun_NuclearTest_cfi")
0015 #process.generator.PGunParameters.ParticleID = cms.vint32(22)
0016 # process.load("FastSimulation.Configuration.SimpleJet_cfi")
0017
0018 #fastsim
0019 process.load("FastSimulation.Configuration.RandomServiceInitialization_cff")
0020 process.load("FastSimulation.Configuration.CommonInputs_cff")
0021 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0022 from Configuration.AlCa.autoCond import autoCond
0023 process.GlobalTag.globaltag = autoCond['mc']
0024
0025
0026 process.ecalRecHit.doMiscalib = True
0027
0028 process.fastSimProducer.SimulateCalorimetry = True
0029 for layer in process.fastSimProducer.detectorDefinition.BarrelLayers:
0030     layer.interactionModels = cms.untracked.vstring("pairProduction", "nuclearInteraction", "bremsstrahlung", "energyLoss", "multipleScattering", "trackerSimHits")
0031 for layer in process.fastSimProducer.detectorDefinition.ForwardLayers:
0032     layer.interactionModels = cms.untracked.vstring("pairProduction", "nuclearInteraction", "bremsstrahlung", "energyLoss", "multipleScattering", "trackerSimHits")
0033 process.famosPileUp.PileUpSimulator.averageNumber = 0.0
0034
0035 process.fastSimProducer.VertexGenerator.BetaStar = 0.00001
0036 process.fastSimProducer.VertexGenerator.SigmaZ = 0.00001
0037
0038 # Parametrized magnetic field (new mapping, 4.0 and 3.8T)
0039 process.load("Configuration.StandardSequences.MagneticField_40T_cff")
0040 #process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0041 process.VolumeBasedMagneticFieldESProducer.useParametrizedTrackerField = True
0042
0043 process.load("RecoParticleFlow.PFProducer.particleFlowSimParticle_cff")
0044 process.load("RecoParticleFlow.PFTracking.particleFlowDisplacedVertexCandidate_cff")
0045 process.load("RecoParticleFlow.PFTracking.particleFlowDisplacedVertex_cff")
0046
0047 process.dump = cms.EDAnalyzer("EventContentAnalyzer")
0048
0049 process.displacedVertexSelector = cms.EDFilter(
0050     "PFDisplacedVertexSelector",
0051     src = cms.InputTag("particleFlowDisplacedVertex"),
0052 #    src = cms.InputTag("toto"),
0053 #    cut = cms.string("position.x>1000.0")
0054     cut = cms.string("isNucl"),
0055     filter = cms.bool(True)
0056     )
0057
0058 process.simVertexSelector = cms.EDFilter(
0059     "SimVertexSelector",
0060     src = cms.InputTag("fastSimProducer"),
0061 #    cut = cms.string("position.x>1000.0")
0062     cut = cms.string("position.rho>2.5"),
0063 #    cut = cms.string(""),
0064     filter = cms.bool(True)
0065     )
0066
0067 #process.selector = cms.Path(process.displacedVertexSelector*process.simVertexSelector)
0068 process.selector = cms.Path(process.displacedVertexSelector)
0069
0070 process.particleFlow.rejectTracks_Bad =  cms.bool(True)
0071 process.particleFlow.rejectTracks_Step45 = cms.bool(True)
0072 process.particleFlowBlock.useConversions = cms.bool(False)
0073 process.particleFlowBlock.useV0 = cms.bool(False)
0074
0075 #process.particleFlow.blocks = cms.InputTag("particleFlowBlock", "WithoutNI", "TEST")
0076
0077
0078
0079
0080 #process.ak7PFJets.src = cms.InputTag("particleFlow", "WithoutNI", "TEST")
0081 #process.ak7PFJets.inputEtMin = 0.0
0082
0083 #process.ak5PFJets.src = cms.InputTag("particleFlow", "WithoutNI", "TEST")
0084 #process.ak5PFJets.inputEtMin = 0.0
0085
0086 #process.PFJetMet = cms.Sequence(
0087 #    process.ak7PFJets +
0088 #    process.ak5PFJets
0089 #    )
0090
0091 #process.famosTauTaggingSequence = cms.Sequence()
0092 #process.famosBTaggingSequence = cms.Sequence()
0093 #process.famosPFTauTaggingSequence = cms.Sequence()
0094
0095
0096
0097 process.particleFlow.usePFNuclearInteractions = cms.bool(True)
0098 process.particleFlowBlock.useNuclear = cms.bool(True)
0099 process.particleFlow.iCfgCandConnector.bCalibPrimary = cms.bool(True)
0100 #process.particleFlow.correctSecondary = cms.bool(True)
0101
0102 process.printList = cms.EDAnalyzer("ParticleListDrawer",
0103                                    src = cms.InputTag("genParticles"),
0104                                    printOnlyHardInteraction  = cms.untracked.bool(False),
0105                                    maxEventsToPrint = cms.untracked.int32(10)
0106                                    )
0107
0108
0109 process.highPtJets = cms.EDFilter(
0110     "CandViewSelector",
0111     src = cms.InputTag("ak7PFJets"),
0112     cut = cms.string( "pt()>60" )
0113     )
0114
0115 process.filterHighPtJets = cms.EDFilter(
0116     "CandCountFilter",
0117     src = cms.InputTag("highPtJets"),
0118     minNumber = cms.uint32(1),
0119     )
0120
0121 process.p1 = cms.Path(
0122 #    process.famosWithCaloTowersAndParticleFlow +
0123     process.generator +
0124     process.famosWithEverything +# +
0125     process.displacedVertexSelector +
0126     process.caloJetMetGen +
0127     process.particleFlowSimParticle+
0128     process.printList#+
0129 #    process.highPtJets+
0130 #    process.filterHighPtJets
0131     )
0132
0133
0134
0135
0136 process.load("FastSimulation.Configuration.EventContent_cff")
0137 process.aod = cms.OutputModule("PoolOutputModule",
0138     process.AODSIMEventContent,
0139     fileName = cms.untracked.string('aod.root')
0140 )
0141
0142 process.load("FastSimulation.Configuration.EventContent_cff")
0143 process.reco = cms.OutputModule("PoolOutputModule",
0144     process.RECOSIMEventContent,
0145     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring("p1")),
0146     fileName = cms.untracked.string('reco_tauGunWithNI.root')
0147 )
0148
0149 process.load("RecoParticleFlow.Configuration.Display_EventContent_cff")
0150 process.display = cms.OutputModule("PoolOutputModule",
0151     process.DisplayEventContent,
0152     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring("p1")),
0153     fileName = cms.untracked.string('display_withoutNI.root')
0154 )
0155
0156 process.outpath = cms.EndPath(process.reco)
0157
0158 process.schedule = cms.Schedule(
0159     process.p1,
0160     process.outpath
0161 )
0162
0163
0164
0165 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0166 process.options = cms.untracked.PSet(
0167     makeTriggerResults = cms.untracked.bool(False),
0168     wantSummary = cms.untracked.bool(False),
0169     Rethrow = cms.untracked.vstring('Unknown',
0170         'ProductNotFound',
0171         'DictionaryNotFound',
0172         'InsertFailure',
0173         'Configuration',
0174         'LogicError',
0175         'UnimplementedFeature',
0176         'InvalidReference',
0177         'NullPointerError',
0178         'NoProductSpecified',
0179         'EventTimeout',
0180         'EventCorruption',
0181         'ModuleFailure',
0182         'ScheduleExecutionFailure',
0183         'EventProcessorFailure',
0184         'FileInPathError',
0185         'FatalRootError',
0186         'NotFound')
0187 )
0188 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0189 #