ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/Configuration/python/BAMBUExample_FastChain.py
Revision: 1.1
Committed: Tue Jul 21 23:39:19 2009 UTC (15 years, 9 months ago) by bendavid
Content type: text/x-python
Branch: MAIN
Log Message:
Added working example

File Contents

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