GenFilters/test/test_d0k4pi_ptdstarcut_gen_13TeV.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('GEN')
0004
0005 # import of standard configurations
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0008 process.load('FWCore.MessageService.MessageLogger_cfi')
0009 process.load('Configuration.EventContent.EventContent_cff')
0010 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0011 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0012 process.load('Configuration.Geometry.GeometrySimDB_cff')
0013 process.load('Configuration.StandardSequences.MagneticField_cff')
0014 process.load('Configuration.StandardSequences.Generator_cff')
0015 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
0016 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0017 process.load('Configuration.StandardSequences.SimIdeal_cff')
0018 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0019 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0020
0021 process.maxEvents = cms.untracked.PSet(
0022     input = cms.untracked.int32(10000)
0023 )
0024
0025 # Input source
0026 process.source = cms.Source("EmptySource")
0027
0028 process.options = cms.untracked.PSet(
0029
0030 )
0031
0032 # Production Info
0033 process.configurationMetadata = cms.untracked.PSet(
0034     version = cms.untracked.string('$Revision: 1.1 $'),
0035     annotation = cms.untracked.string('D*+ -> D0 pi+ -> (K-pi+pi+pi-) pi+ at 13TeV'),
0036     name = cms.untracked.string('$Source: test_d0k4pi_ptdstarcut_gen_13TeV.py,v $')
0037 )
0038
0039 # Output definition
0040
0041 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0042     splitLevel = cms.untracked.int32(0),
0043     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0044     outputCommands = process.RAWSIMEventContent.outputCommands,
0045     fileName = cms.untracked.string('4pi_ptdstarcut_gen.root'),
0046     dataset = cms.untracked.PSet(
0047         filterName = cms.untracked.string(''),
0048         dataTier = cms.untracked.string('GEN')
0049     ),
0050     SelectEvents = cms.untracked.PSet(
0051         SelectEvents = cms.vstring('generation_step')
0052     )
0053 )
0054
0055
0056 # Additional output definition
0057
0058 # Other statements
0059 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0060 from Configuration.AlCa.GlobalTag import GlobalTag
0061 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0062 #process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:startup', '')
0063 #from Configuration.Generator.PythiaUESettings_cfi import * #for 8 TeV
0064 from Configuration.Generator.PythiaUEZ2starSettings_cfi import * #for 13 TeV
0065 #################
0066 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0067     pythiaPylistVerbosity = cms.untracked.int32(1),
0068     filterEfficiency = cms.untracked.double(1.0),
0069     pythiaHepMCVerbosity = cms.untracked.bool(False),
0070     maxEventsToPrint = cms.untracked.int32(2),
0071     comEnergy = cms.double(13000.0),
0072     ExternalDecays = cms.PSet(
0073         EvtGen = cms.untracked.PSet(
0074              operates_on_particles = cms.vint32( 0 ), # 0 (zero) means default list (hardcoded)
0075                                                       # you can put here the list of particles (PDG IDs)
0076                                                       # that you want decayed by EvtGen
0077          use_default_decay = cms.untracked.bool(False),
0078              decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
0079               # decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY.DEC'),
0080               particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
0081               user_decay_file = cms.FileInPath('Dati13/TrackEfficiency/Dstar_D0_K3pi.dec'),
0082               list_forced_decays = cms.vstring('MyD*+','MyD*-')
0083               ),
0084          parameterSets = cms.vstring('EvtGen')
0085     ),
0086     PythiaParameters = cms.PSet(
0087       #process.pythiaUESettingsBlock,
0088       pythiaUESettingsBlock,
0089       ccbarSettings= cms.vstring('MSEL=4     ! ccbar '),
0090       # This is a vector of ParameterSet names to be read, in this order
0091       parameterSets = cms.vstring('pythiaUESettings','ccbarSettings')
0092     )
0093 )
0094
0095 #process.D0DecayFilter = cms.EDFilter("PythiaDauFilter",
0096 #    ParticleID = cms.untracked.int32(421),
0097 #    ChargeConjugation = cms.untracked.bool(True),
0098 #    MinEta = cms.untracked.double(-100.),
0099 #    MaxEta = cms.untracked.double(100.),
0100 #    DaughterIDs = cms.untracked.vint32(-321,211,211,-211),
0101 #    NumberDaughters = cms.untracked.int32(4)
0102 #)
0103
0104 process.DstarFilter = cms.EDFilter("PythiaMomDauFilter",
0105     ParticleID = cms.untracked.int32(413),
0106     DaughterID = cms.untracked.int32(421),
0107     ChargeConjugation = cms.untracked.bool(True),
0108     MinEta = cms.untracked.double(-100.),
0109     MaxEta = cms.untracked.double(100.),
0110     DaughterIDs = cms.untracked.vint32(421,211),
0111     NumberDaughters = cms.untracked.int32(2),
0112     MomMinPt = cms.untracked.double(3.9),
0113     NumberDescendants = cms.untracked.int32(4),
0114     DescendantsIDs = cms.untracked.vint32(-321,211,211,-211)
0115 )
0116
0117 #################
0118
0119 #process.ProductionFilterSequence = cms.Sequence(process.generator*process.DstarFilter*process.D0DecayFilter)
0120 process.ProductionFilterSequence = cms.Sequence(process.generator*process.DstarFilter)
0121
0122 # Path and EndPath definitions
0123 process.generation_step = cms.Path(process.pgen)
0124 #process.simulation_step = cms.Path(process.psim)
0125 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0126 process.endjob_step = cms.EndPath(process.endOfProcess)
0127 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0128
0129 # Schedule definition
0130 #process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.endjob_step,process.RAWSIMoutput_step)
0131 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.endjob_step,process.RAWSIMoutput_step)
0132 # filter all path with the production filter sequence
0133 for path in process.paths:
0134     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
0135
0136 # customisation of the process.
0137
0138 # Automatic addition of the customisation function from SLHCUpgradeSimulations.Configuration.postLS1Customs
0139 from SLHCUpgradeSimulations.Configuration.postLS1Customs import customisePostLS1
0140
0141 #call to customisation function customisePostLS1 imported from SLHCUpgradeSimulations.Configuration.postLS1Customs
0142 process = customisePostLS1(process)
0143
0144 # End of customisation functions
0145
0146
0147