| 1 |
# $Id: BAMBUExample_FastChain.py,v 1.6 2009/11/05 21:38:09 bendavid Exp $
|
| 2 |
#
|
| 3 |
# Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
|
| 4 |
# with command line options:
|
| 5 |
# TTbar_cfi -s GEN,FASTSIM --conditions FrontierConditions_GlobalTag,MC_31X_V3::All
|
| 6 |
# -n 100 --no_exec --no_output --pileup NoPileUp
|
| 7 |
|
| 8 |
import FWCore.ParameterSet.Config as cms
|
| 9 |
|
| 10 |
process = cms.Process('HLT')
|
| 11 |
|
| 12 |
# import of standard configurations
|
| 13 |
process.load('FastSimulation/Configuration/RandomServiceInitialization_cff')
|
| 14 |
process.load('FastSimulation.PileUpProducer.PileUpSimulator10TeV_cfi')
|
| 15 |
process.load('Configuration/StandardSequences/MagneticField_38T_cff')
|
| 16 |
process.load('Configuration/StandardSequences/Generator_cff')
|
| 17 |
process.load('FastSimulation/Configuration/FamosSequences_cff')
|
| 18 |
process.load('FastSimulation/Configuration/HLT_1E31_cff')
|
| 19 |
process.load('IOMC.EventVertexGenerators.VtxSmearedParameters_cfi')
|
| 20 |
process.load('FastSimulation/Configuration/CommonInputs_cff')
|
| 21 |
process.load('FastSimulation/Configuration/EventContent_cff')
|
| 22 |
|
| 23 |
process.configurationMetadata = cms.untracked.PSet(
|
| 24 |
version = cms.untracked.string('CMSSW_3_2_1 + Mit_012d'),
|
| 25 |
annotation = cms.untracked.string('TTbar'),
|
| 26 |
name = cms.untracked.string('BambuFastChainExample')
|
| 27 |
)
|
| 28 |
|
| 29 |
process.maxEvents = cms.untracked.PSet(
|
| 30 |
input = cms.untracked.int32(100)
|
| 31 |
)
|
| 32 |
process.options = cms.untracked.PSet(
|
| 33 |
Rethrow = cms.untracked.vstring('ProductNotFound')
|
| 34 |
)
|
| 35 |
|
| 36 |
# input source
|
| 37 |
process.source = cms.Source("EmptySource")
|
| 38 |
|
| 39 |
# output definition
|
| 40 |
#process.output = cms.OutputModule("PoolOutputModule",
|
| 41 |
# splitLevel = cms.untracked.int32(0),
|
| 42 |
# outputCommands = process.RECOSIMEventContent.outputCommands,
|
| 43 |
# fileName = cms.untracked.string('TTbar_cfi_GEN_FASTSIM.root'),
|
| 44 |
# dataset = cms.untracked.PSet(
|
| 45 |
# dataTier = cms.untracked.string(''),
|
| 46 |
# filterName = cms.untracked.string('')
|
| 47 |
# ),
|
| 48 |
# SelectEvents = cms.untracked.PSet(
|
| 49 |
# SelectEvents = cms.vstring('generation_step')
|
| 50 |
# )
|
| 51 |
#)
|
| 52 |
|
| 53 |
# Additional output definition
|
| 54 |
|
| 55 |
# other statements
|
| 56 |
process.famosPileUp.PileUpSimulator = process.PileUpSimulatorBlock.PileUpSimulator
|
| 57 |
process.famosPileUp.PileUpSimulator.averageNumber = 0
|
| 58 |
process.famosSimHits.SimulateCalorimetry = True
|
| 59 |
process.famosSimHits.SimulateTracking = True
|
| 60 |
process.famosSimHits.ActivateDecays.comEnergy = 10000
|
| 61 |
process.simulation = cms.Sequence(process.simulationWithFamos)
|
| 62 |
process.HLTEndSequence = cms.Sequence(process.reconstructionWithFamos)
|
| 63 |
|
| 64 |
# set correct vertex smearing
|
| 65 |
process.Early10TeVCollisionVtxSmearingParameters.type = cms.string("BetaFunc")
|
| 66 |
process.famosSimHits.VertexGenerator = process.Early10TeVCollisionVtxSmearingParameters
|
| 67 |
process.famosPileUp.VertexGenerator = process.Early10TeVCollisionVtxSmearingParameters
|
| 68 |
|
| 69 |
# apply tracker misalignment
|
| 70 |
process.famosSimHits.ApplyAlignment = True
|
| 71 |
process.misalignedTrackerGeometry.applyAlignment = True
|
| 72 |
|
| 73 |
process.GlobalTag.globaltag = 'MC_31X_V9::All'
|
| 74 |
process.generator = cms.EDFilter("Pythia6GeneratorFilter",
|
| 75 |
pythiaPylistVerbosity = cms.untracked.int32(0),
|
| 76 |
filterEfficiency = cms.untracked.double(1.0),
|
| 77 |
pythiaHepMCVerbosity = cms.untracked.bool(False),
|
| 78 |
comEnergy = cms.double(10000.0),
|
| 79 |
maxEventsToPrint = cms.untracked.int32(0),
|
| 80 |
PythiaParameters = cms.PSet(
|
| 81 |
pythiaUESettings = cms.vstring('MSTJ(11)=3 ! Choice of the fragmentation function',
|
| 82 |
'MSTJ(22)=2 ! Decay those unstable particles',
|
| 83 |
'PARJ(71)=10 . ! for which ctau 10 mm',
|
| 84 |
'MSTP(2)=1 ! which order running alphaS',
|
| 85 |
'MSTP(33)=0 ! no K factors in hard cross sections',
|
| 86 |
'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
|
| 87 |
'MSTP(52)=2 ! work with LHAPDF',
|
| 88 |
'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default',
|
| 89 |
'MSTP(82)=4 ! Defines the multi-parton model',
|
| 90 |
'MSTU(21)=1 ! Check on possible errors during program execution',
|
| 91 |
'PARP(82)=1.8387 ! pt cutoff for multiparton interactions',
|
| 92 |
'PARP(89)=1960. ! sqrts for which PARP82 is set',
|
| 93 |
'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter',
|
| 94 |
'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter',
|
| 95 |
'PARP(90)=0.16 ! Multiple interactions: rescaling power',
|
| 96 |
'PARP(67)=2.5 ! amount of initial-state radiation',
|
| 97 |
'PARP(85)=1.0 ! gluon prod. mechanism in MI',
|
| 98 |
'PARP(86)=1.0 ! gluon prod. mechanism in MI',
|
| 99 |
'PARP(62)=1.25 ! ',
|
| 100 |
'PARP(64)=0.2 ! ',
|
| 101 |
'MSTP(91)=1 !',
|
| 102 |
'PARP(91)=2.1 ! kt distribution',
|
| 103 |
'PARP(93)=15.0 ! '),
|
| 104 |
processParameters = cms.vstring('MSEL = 0 ! User defined processes',
|
| 105 |
'MSUB(81) = 1 ! qqbar to QQbar',
|
| 106 |
'MSUB(82) = 1 ! gg to QQbar',
|
| 107 |
'MSTP(7) = 6 ! flavour = top',
|
| 108 |
'PMAS(6,1) = 175. ! top quark mass'),
|
| 109 |
parameterSets = cms.vstring('pythiaUESettings',
|
| 110 |
'processParameters')
|
| 111 |
)
|
| 112 |
)
|
| 113 |
#process.ProductionFilterSequence = cms.Sequence(process.generator)
|
| 114 |
|
| 115 |
# path and end path definitions
|
| 116 |
process.generation_step = cms.Path(process.generator +
|
| 117 |
cms.SequencePlaceholder("randomEngineStateProducer") +
|
| 118 |
process.GeneInfo+process.genJetMET)
|
| 119 |
process.reconstruction = cms.Path(process.reconstructionWithFamos)
|
| 120 |
#process.out_step = cms.EndPath(process.output)
|
| 121 |
|
| 122 |
process.add_(cms.Service("ObjectService"))
|
| 123 |
|
| 124 |
process.load("MitProd.TreeFiller.MitTreeFiller_cfi")
|
| 125 |
|
| 126 |
# enable fillers of MC Truth information
|
| 127 |
process.MitTreeFiller.MCParticles.active = True
|
| 128 |
process.MitTreeFiller.MCEventInfo.active = True
|
| 129 |
process.MitTreeFiller.IC5GenJets.active = True
|
| 130 |
process.MitTreeFiller.SC5GenJets.active = True
|
| 131 |
process.MitTreeFiller.SC7GenJets.active = True
|
| 132 |
process.MitTreeFiller.KT4GenJets.active = True
|
| 133 |
process.MitTreeFiller.KT6GenJets.active = True
|
| 134 |
process.MitTreeFiller.GenMet.active = True
|
| 135 |
|
| 136 |
# hack out unavailable stuff, pending proper fastsim and aod sequences
|
| 137 |
process.MitTreeFiller.MetaInfos.hltActive = False
|
| 138 |
process.MitTreeFiller.MCParticles.simActive = False
|
| 139 |
process.MitTreeFiller.ConversionInOutTracks.active = False
|
| 140 |
process.MitTreeFiller.ConversionOutInTracks.active = False
|
| 141 |
process.MitTreeFiller.Conversions.active = False
|
| 142 |
process.MitTreeFiller.ConversionInOutElectronsStable.active = False
|
| 143 |
process.MitTreeFiller.ConversionOutInElectronsStable.active = False
|
| 144 |
|
| 145 |
process.bambu_step = cms.Path(process.MitTreeFiller)
|
| 146 |
|
| 147 |
# schedule definition
|
| 148 |
process.schedule = cms.Schedule(process.generation_step)
|
| 149 |
process.schedule.extend(process.HLTSchedule)
|
| 150 |
process.schedule.extend([process.reconstruction,process.bambu_step])
|
| 151 |
|
| 152 |
# special treatment in case of production filter sequence
|
| 153 |
#for path in process.paths:
|
| 154 |
#getattr(process,path)._seq = process.ProductionFilterSequence*getattr(process,path)._seq
|
