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.particleFlow.iCfgCandConnector.bCalibPrimary = cms.bool(True)
0099 #process.particleFlow.correctSecondary = cms.bool(True)
0100
0101 process.printList = cms.EDAnalyzer("ParticleListDrawer",
0102                                    src = cms.InputTag("genParticles"),
0103                                    printOnlyHardInteraction  = cms.untracked.bool(False),
0104                                    maxEventsToPrint = cms.untracked.int32(10)
0105                                    )
0106
0107
0108 process.highPtJets = cms.EDFilter(
0109     "CandViewSelector",
0110     src = cms.InputTag("ak7PFJets"),
0111     cut = cms.string( "pt()>60" )
0112     )
0113
0114 process.filterHighPtJets = cms.EDFilter(
0115     "CandCountFilter",
0116     src = cms.InputTag("highPtJets"),
0117     minNumber = cms.uint32(1),
0118     )
0119
0120 process.p1 = cms.Path(
0121 #    process.famosWithCaloTowersAndParticleFlow +
0122     process.generator +
0123     process.famosWithEverything +# +
0124     process.displacedVertexSelector +
0125     process.caloJetMetGen +
0126     process.particleFlowSimParticle+
0127     process.printList#+
0128 #    process.highPtJets+
0129 #    process.filterHighPtJets
0130     )
0131
0132
0133
0134
0135 process.load("FastSimulation.Configuration.EventContent_cff")
0136 process.aod = cms.OutputModule("PoolOutputModule",
0137     process.AODSIMEventContent,
0138     fileName = cms.untracked.string('aod.root')
0139 )
0140
0141 process.load("FastSimulation.Configuration.EventContent_cff")
0142 process.reco = cms.OutputModule("PoolOutputModule",
0143     process.RECOSIMEventContent,
0144     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring("p1")),
0145     fileName = cms.untracked.string('reco_tauGunWithNI.root')
0146 )
0147
0148 process.load("RecoParticleFlow.Configuration.Display_EventContent_cff")
0149 process.display = cms.OutputModule("PoolOutputModule",
0150     process.DisplayEventContent,
0151     SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring("p1")),
0152     fileName = cms.untracked.string('display_withoutNI.root')
0153 )
0154
0155 process.outpath = cms.EndPath(process.reco)
0156
0157 process.schedule = cms.Schedule(
0158     process.p1,
0159     process.outpath
0160 )
0161
0162
0163
0164 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0165 process.options = cms.untracked.PSet(
0166     makeTriggerResults = cms.untracked.bool(False),
0167     wantSummary = cms.untracked.bool(False),
0168     Rethrow = cms.untracked.vstring('Unknown',
0169         'ProductNotFound',
0170         'DictionaryNotFound',
0171         'InsertFailure',
0172         'Configuration',
0173         'LogicError',
0174         'UnimplementedFeature',
0175         'InvalidReference',
0176         'NullPointerError',
0177         'NoProductSpecified',
0178         'EventTimeout',
0179         'EventCorruption',
0180         'ModuleFailure',
0181         'ScheduleExecutionFailure',
0182         'EventProcessorFailure',
0183         'FileInPathError',
0184         'FatalRootError',
0185         'NotFound')
0186 )
0187 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0188 #