Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 11:21:23

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 #