Configuration/python/BAMBUExample_FastChain.py

# $Id: BAMBUExample_FastChain.py,v 1.2 2009/07/22 11:26:19 loizides Exp $
#
# Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v 
# with command line options:
#  TTbar_cfi -s GEN,FASTSIM --conditions FrontierConditions_GlobalTag,MC_31X_V3::All
#  -n 100 --no_exec --no_output --pileup NoPileUp

import FWCore.ParameterSet.Config as cms

process = cms.Process('HLT')

# import of standard configurations
process.load('FastSimulation/Configuration/RandomServiceInitialization_cff')
process.load('FastSimulation.PileUpProducer.PileUpSimulator10TeV_cfi')
process.load('Configuration/StandardSequences/MagneticField_38T_cff')
process.load('Configuration/StandardSequences/Generator_cff')
process.load('FastSimulation/Configuration/FamosSequences_cff')
process.load('FastSimulation/Configuration/HLT_1E31_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedParameters_cfi')
process.load('FastSimulation/Configuration/CommonInputs_cff')
process.load('FastSimulation/Configuration/EventContent_cff')

process.configurationMetadata = cms.untracked.PSet(
    version = cms.untracked.string('CMSSW_3_2_1 + Mit_010'),
    annotation = cms.untracked.string('TTbar_cfi nevts:100'),
    name = cms.untracked.string('BambuFastChainExample')
)

process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(100)
)
process.options = cms.untracked.PSet(
    Rethrow = cms.untracked.vstring('ProductNotFound')
)

# input source
process.source = cms.Source("EmptySource")

# output definition
#process.output = cms.OutputModule("PoolOutputModule",
#    splitLevel = cms.untracked.int32(0),
#    outputCommands = process.RECOSIMEventContent.outputCommands,
#    fileName = cms.untracked.string('TTbar_cfi_GEN_FASTSIM.root'),
#    dataset = cms.untracked.PSet(
#       dataTier = cms.untracked.string(''),
#       filterName = cms.untracked.string('')
#    ),
#    SelectEvents = cms.untracked.PSet(
#       SelectEvents = cms.vstring('generation_step')
#    )
#)

# Additional output definition

# other statements
process.famosPileUp.PileUpSimulator = process.PileUpSimulatorBlock.PileUpSimulator
process.famosPileUp.PileUpSimulator.averageNumber = 0
process.famosSimHits.SimulateCalorimetry = True
process.famosSimHits.SimulateTracking = True
process.famosSimHits.ActivateDecays.comEnergy = 10000
process.simulation = cms.Sequence(process.simulationWithFamos)
process.HLTEndSequence = cms.Sequence(process.reconstructionWithFamos)

# set correct vertex smearing
process.Early10TeVCollisionVtxSmearingParameters.type = cms.string("BetaFunc")
process.famosSimHits.VertexGenerator = process.Early10TeVCollisionVtxSmearingParameters
process.famosPileUp.VertexGenerator = process.Early10TeVCollisionVtxSmearingParameters

# apply tracker misalignment
process.famosSimHits.ApplyAlignment = True
process.misalignedTrackerGeometry.applyAlignment = True

process.GlobalTag.globaltag = 'MC_31X_V3::All'
process.generator = cms.EDFilter("Pythia6GeneratorFilter",
    pythiaPylistVerbosity = cms.untracked.int32(0),
    filterEfficiency = cms.untracked.double(1.0),
    pythiaHepMCVerbosity = cms.untracked.bool(False),
    comEnergy = cms.double(10000.0),
    maxEventsToPrint = cms.untracked.int32(0),
    PythiaParameters = cms.PSet(
        pythiaUESettings = cms.vstring('MSTJ(11)=3     ! Choice of the fragmentation function', 
            'MSTJ(22)=2     ! Decay those unstable particles', 
            'PARJ(71)=10 .  ! for which ctau  10 mm', 
            'MSTP(2)=1      ! which order running alphaS', 
            'MSTP(33)=0     ! no K factors in hard cross sections', 
            'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)', 
            'MSTP(52)=2     ! work with LHAPDF', 
            'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default', 
            'MSTP(82)=4     ! Defines the multi-parton model', 
            'MSTU(21)=1     ! Check on possible errors during program execution', 
            'PARP(82)=1.8387   ! pt cutoff for multiparton interactions', 
            'PARP(89)=1960. ! sqrts for which PARP82 is set', 
            'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter', 
            'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter', 
            'PARP(90)=0.16  ! Multiple interactions: rescaling power', 
            'PARP(67)=2.5    ! amount of initial-state radiation', 
            'PARP(85)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(86)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(62)=1.25   ! ', 
            'PARP(64)=0.2    ! ', 
            'MSTP(91)=1      !', 
            'PARP(91)=2.1   ! kt distribution', 
            'PARP(93)=15.0  ! '),
        processParameters = cms.vstring('MSEL      = 0     ! User defined processes', 
            'MSUB(81)  = 1     ! qqbar to QQbar', 
            'MSUB(82)  = 1     ! gg to QQbar', 
            'MSTP(7)   = 6     ! flavour = top', 
            'PMAS(6,1) = 175.  ! top quark mass'),
        parameterSets = cms.vstring('pythiaUESettings', 
            'processParameters')
    )
)
#process.ProductionFilterSequence = cms.Sequence(process.generator)

# path and end path definitions
process.generation_step = cms.Path(process.generator +
                                   cms.SequencePlaceholder("randomEngineStateProducer") +
                                   process.GeneInfo+process.genJetMET)
process.reconstruction = cms.Path(process.reconstructionWithFamos)
#process.out_step = cms.EndPath(process.output)

# load MitTreeFiller 
process.TreeService = cms.Service("TreeService",
    fileNames = cms.untracked.vstring('mit-full'),
)
process.add_(cms.Service("ObjectService"))

process.load("MitProd.TreeFiller.MitTreeFiller_cfi")

#enable fillers of MC Truth information
process.MitTreeFiller.MCParticles.active = True
process.MitTreeFiller.MCEventInfo.active = True
process.MitTreeFiller.IC5GenJets.active  = True
process.MitTreeFiller.SC5GenJets.active  = True
process.MitTreeFiller.SC7GenJets.active  = True
process.MitTreeFiller.KT4GenJets.active  = True
process.MitTreeFiller.KT6GenJets.active  = True
process.MitTreeFiller.GenMet.active      = True

# hack out unavailable stuff, pending proper fastsim and aod sequences
process.MitTreeFiller.MetaInfos.hltActive                   = False
process.MitTreeFiller.MCParticles.simActive                 = False
process.MitTreeFiller.ConversionInOutTracks.active          = False
process.MitTreeFiller.ConversionOutInTracks.active          = False
process.MitTreeFiller.Conversions.active                    = False
process.MitTreeFiller.ConversionInOutElectronsStable.active = False
process.MitTreeFiller.ConversionOutInElectronsStable.active = False

process.bambu_step  = cms.Path(process.MitTreeFiller)

# schedule definition
process.schedule = cms.Schedule(process.generation_step)
process.schedule.extend(process.HLTSchedule)
process.schedule.extend([process.reconstruction,process.bambu_step])

# special treatment in case of production filter sequence  
#for path in process.paths: 
#getattr(process,path)._seq = process.ProductionFilterSequence*getattr(process,path)._seq
Revision:	1.3
Committed:	Wed Jul 22 15:16:24 2009 UTC (15 years, 9 months ago) by bendavid
Content type:	text/x-python
Branch:	MAIN
Changes since 1.2:	+13 -3 lines
Log Message:	Fixed example
#	User	Rev	Content
1	bendavid	1.3	# $Id: BAMBUExample_FastChain.py,v 1.2 2009/07/22 11:26:19 loizides Exp $
2	loizides	1.2	#
3	bendavid	1.1	# Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
4	loizides	1.2	# 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	bendavid	1.1	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	loizides	1.2	version = cms.untracked.string('CMSSW_3_2_1 + Mit_010'),
25	bendavid	1.1	annotation = cms.untracked.string('TTbar_cfi nevts:100'),
26	loizides	1.2	name = cms.untracked.string('BambuFastChainExample')
27	bendavid	1.1	)
28	loizides	1.2
29	bendavid	1.1	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	loizides	1.2
36			# input source
37	bendavid	1.1	process.source = cms.Source("EmptySource")
38
39	loizides	1.2	# output definition
40	bendavid	1.1	#process.output = cms.OutputModule("PoolOutputModule",
41	loizides	1.2	# 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	bendavid	1.1	#)
52
53			# Additional output definition
54
55	loizides	1.2	# other statements
56	bendavid	1.1	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	loizides	1.2
69			# apply tracker misalignment
70	bendavid	1.1	process.famosSimHits.ApplyAlignment = True
71			process.misalignedTrackerGeometry.applyAlignment = True
72
73			process.GlobalTag.globaltag = 'MC_31X_V3::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	loizides	1.2	# path and end path definitions
116			process.generation_step = cms.Path(process.generator +
117			cms.SequencePlaceholder("randomEngineStateProducer") +
118			process.GeneInfo+process.genJetMET)
119	bendavid	1.1	process.reconstruction = cms.Path(process.reconstructionWithFamos)
120			#process.out_step = cms.EndPath(process.output)
121
122	loizides	1.2	# load MitTreeFiller
123	bendavid	1.1	process.TreeService = cms.Service("TreeService",
124			fileNames = cms.untracked.vstring('mit-full'),
125			)
126			process.add_(cms.Service("ObjectService"))
127
128	bendavid	1.3	process.load("MitProd.TreeFiller.MitTreeFiller_cfi")
129
130			#enable fillers of MC Truth information
131			process.MitTreeFiller.MCParticles.active = True
132			process.MitTreeFiller.MCEventInfo.active = True
133			process.MitTreeFiller.IC5GenJets.active = True
134			process.MitTreeFiller.SC5GenJets.active = True
135			process.MitTreeFiller.SC7GenJets.active = True
136			process.MitTreeFiller.KT4GenJets.active = True
137			process.MitTreeFiller.KT6GenJets.active = True
138			process.MitTreeFiller.GenMet.active = True
139	bendavid	1.1
140	loizides	1.2	# hack out unavailable stuff, pending proper fastsim and aod sequences
141	bendavid	1.1	process.MitTreeFiller.MetaInfos.hltActive = False
142			process.MitTreeFiller.MCParticles.simActive = False
143			process.MitTreeFiller.ConversionInOutTracks.active = False
144			process.MitTreeFiller.ConversionOutInTracks.active = False
145			process.MitTreeFiller.Conversions.active = False
146			process.MitTreeFiller.ConversionInOutElectronsStable.active = False
147			process.MitTreeFiller.ConversionOutInElectronsStable.active = False
148
149	bendavid	1.3	process.bambu_step = cms.Path(process.MitTreeFiller)
150	bendavid	1.1
151	loizides	1.2	# schedule definition
152	bendavid	1.1	process.schedule = cms.Schedule(process.generation_step)
153			process.schedule.extend(process.HLTSchedule)
154			process.schedule.extend([process.reconstruction,process.bambu_step])
155	loizides	1.2
156	bendavid	1.1	# special treatment in case of production filter sequence
157			#for path in process.paths:
158	loizides	1.2	#getattr(process,path)._seq = process.ProductionFilterSequence*getattr(process,path)._seq