ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/TreeFiller/python/FullChainExample.py
Revision: 1.1
Committed: Mon Oct 27 19:58:27 2008 UTC (16 years, 6 months ago) by bendavid
Content type: text/x-python
Branch: MAIN
CVS Tags: Mit_006
Log Message:
Added python full chain example

File Contents

# Content
1 # Auto generated configuration file
2 # using:
3 # Revision: 1.77.2.1
4 # Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
5 # with command line options: Configuration/GenProduction/python/PYTHIA6_EWK_Zmumu_10TeV_cff.py --step=GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,RECO --eventcontent=FEVTSIM --datatier=GEN-SIM-RAW-RECO --conditions=FrontierConditions_GlobalTag,IDEAL_V9::All -n 10 --no_exec
6 import FWCore.ParameterSet.Config as cms
7
8 process = cms.Process('RECO')
9
10 # import of standard configurations
11 process.load('Configuration/StandardSequences/Services_cff')
12 process.load('FWCore/MessageService/MessageLogger_cfi')
13 process.load('Configuration/StandardSequences/MixingNoPileUp_cff')
14 process.load('Configuration/StandardSequences/GeometryPilot2_cff')
15 process.load('Configuration/StandardSequences/MagneticField_38T_cff')
16 process.load('Configuration/StandardSequences/Generator_cff')
17 process.load('Configuration/StandardSequences/VtxSmearedEarly10TeVCollision_cff')
18 process.load('Configuration/StandardSequences/Sim_cff')
19 process.load('Configuration/StandardSequences/Digi_cff')
20 process.load('Configuration/StandardSequences/SimL1Emulator_cff')
21 process.load('L1TriggerConfig/L1GtConfigProducers/Luminosity/lumi1030.L1Menu2008_2E30_Unprescaled_cff')
22 process.load('Configuration/StandardSequences/DigiToRaw_cff')
23 process.load('Configuration/StandardSequences/RawToDigi_cff')
24 process.load('Configuration/StandardSequences/Reconstruction_cff')
25 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
26 process.load('Configuration/EventContent/EventContent_cff')
27
28 process.configurationMetadata = cms.untracked.PSet(
29 version = cms.untracked.string('$Revision: 1.3 $'),
30 annotation = cms.untracked.string('Z/gamma* to mumu sample, M(mumu) gt40 GeV -- sqrt(s) = 10TeV'),
31 name = cms.untracked.string('$Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/GenProduction/python/PYTHIA6_EWK_Zmumu_10TeV_cff.py,v $')
32 )
33 process.maxEvents = cms.untracked.PSet(
34 input = cms.untracked.int32(10)
35 )
36 process.options = cms.untracked.PSet(
37 Rethrow = cms.untracked.vstring('ProductNotFound')
38 )
39 # Input source
40 process.source = cms.Source("PythiaSource",
41 pythiaPylistVerbosity = cms.untracked.int32(0),
42 filterEfficiency = cms.untracked.double(0.509),
43 pythiaHepMCVerbosity = cms.untracked.bool(False),
44 comEnergy = cms.untracked.double(10000.0),
45 crossSection = cms.untracked.double(1232.0),
46 maxEventsToPrint = cms.untracked.int32(0),
47 PythiaParameters = cms.PSet(
48 pythiaUESettings = cms.vstring('MSTJ(11)=3 ! Choice of the fragmentation function',
49 'MSTJ(22)=2 ! Decay those unstable particles',
50 'PARJ(71)=10 . ! for which ctau 10 mm',
51 'MSTP(2)=1 ! which order running alphaS',
52 'MSTP(33)=0 ! no K factors in hard cross sections',
53 'MSTP(51)=10042 ! CTEQ6L1 structure function chosen',
54 'MSTP(52)=2 ! work with LHAPDF',
55 'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default',
56 'MSTP(82)=4 ! Defines the multi-parton model',
57 'MSTU(21)=1 ! Check on possible errors during program execution',
58 'PARP(82)=1.8387 ! pt cutoff for multiparton interactions',
59 'PARP(89)=1960. ! sqrts for which PARP82 is set',
60 'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter',
61 'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter',
62 'PARP(90)=0.16 ! Multiple interactions: rescaling power',
63 'PARP(67)=2.5 ! amount of initial-state radiation',
64 'PARP(85)=1.0 ! gluon prod. mechanism in MI',
65 'PARP(86)=1.0 ! gluon prod. mechanism in MI',
66 'PARP(62)=1.25 ! ',
67 'PARP(64)=0.2 ! ',
68 'MSTP(91)=1 !',
69 'PARP(91)=2.1 ! kt distribution',
70 'PARP(93)=15.0 ! '),
71 processParameters = cms.vstring('MSEL=0 !User defined processes',
72 'MSUB(1)=1 !Incl Z0/gamma* production',
73 'MSTP(43)=3 !Both Z0 and gamma*',
74 'MDME(174,1)=0 !Z decay into d dbar',
75 'MDME(175,1)=0 !Z decay into u ubar',
76 'MDME(176,1)=0 !Z decay into s sbar',
77 'MDME(177,1)=0 !Z decay into c cbar',
78 'MDME(178,1)=0 !Z decay into b bbar',
79 'MDME(179,1)=0 !Z decay into t tbar',
80 'MDME(182,1)=0 !Z decay into e- e+',
81 'MDME(183,1)=0 !Z decay into nu_e nu_ebar',
82 'MDME(184,1)=1 !Z decay into mu- mu+',
83 'MDME(185,1)=0 !Z decay into nu_mu nu_mubar',
84 'MDME(186,1)=0 !Z decay into tau- tau+',
85 'MDME(187,1)=0 !Z decay into nu_tau nu_taubar',
86 'CKIN(1)=40. !Minimum sqrt(s_hat) value (=Z mass)'),
87 parameterSets = cms.vstring('pythiaUESettings',
88 'processParameters')
89 )
90 )
91
92 # Output definition
93 process.output = cms.OutputModule("PoolOutputModule",
94 outputCommands = process.FEVTSIMEventContent.outputCommands,
95 fileName = cms.untracked.string('PYTHIA6_EWK_Zmumu_10TeV_cff_py_GEN_SIM_DIGI_L1_DIGI2RAW_RAW2DIGI_RECO.root'),
96 dataset = cms.untracked.PSet(
97 dataTier = cms.untracked.string('GEN-SIM-RAW-RECO'),
98 filterName = cms.untracked.string('')
99 ),
100 SelectEvents = cms.untracked.PSet(
101 SelectEvents = cms.vstring('generation_step')
102 )
103 )
104
105 # Additional output definition
106
107 # Other statements
108 process.GlobalTag.globaltag = 'IDEAL_V9::All'
109
110 #Load MitTreeFiller
111 process.TreeService = cms.Service("TreeService",
112 fileNames = cms.untracked.vstring('mit-full')
113 )
114 process.add_(cms.Service("ObjectService"))
115 process.load("MitProd.TreeFiller.MitTreeFiller_cfi")
116 #disable HLT filler since the full chain example does not (and cannot at the moment) contain HLT info
117 process.MitTreeFiller.MetaInfos.hltActive = False
118
119
120 #Load ElectronID information
121 process.load("MitProd.TreeFiller.ElectronID_cfi")
122
123 #For Jet Corrections
124 process.load("MitProd.TreeFiller.JetCorrections_cfi")
125 process.prefer("L3JetCorrectorIcone5")
126 #enable Jet Corrections for all of our Jet collections
127 process.MitTreeFiller.CaloJets.jetCorrectionsActive = cms.untracked.bool(True)
128 process.MitTreeFiller.ItrCone5Jets.jetCorrectionsActive = cms.untracked.bool(True)
129 process.MitTreeFiller.SisCone5Jets.jetCorrectionsActive = cms.untracked.bool(True)
130 process.MitTreeFiller.SisCone7Jets.jetCorrectionsActive = cms.untracked.bool(True)
131 process.MitTreeFiller.Kt4Jets.jetCorrectionsActive = cms.untracked.bool(True)
132 process.MitTreeFiller.Kt6Jets.jetCorrectionsActive = cms.untracked.bool(True)
133
134 #For JetPlusTracks
135 process.load("JetMETCorrections.Configuration.JetPlusTrackCorrections_cff")
136 process.load("JetMETCorrections.Configuration.ZSPJetCorrections152_cff")
137 process.MitTreeFiller.IC5JetPlusTrack.active = cms.untracked.bool(True)
138
139 #Load Met Corrections
140 process.load("MitProd.TreeFiller.MetCorrections_cfi")
141
142 #Load Flavor Matching Information
143 process.load("MitProd.TreeFiller.JetsMCFlavourMatching_cfi")
144 #Enable Flavor matching for Reco Jets and GenJets
145 process.MitTreeFiller.CaloJets.flavorMatchingActive = cms.untracked.bool(True)
146 process.MitTreeFiller.ItrCone5Jets.flavorMatchingActive = cms.untracked.bool(True)
147 process.MitTreeFiller.SisCone5Jets.flavorMatchingActive = cms.untracked.bool(True)
148 process.MitTreeFiller.SisCone7Jets.flavorMatchingActive = cms.untracked.bool(True)
149 process.MitTreeFiller.Kt4Jets.flavorMatchingActive = cms.untracked.bool(True)
150 process.MitTreeFiller.Kt6Jets.flavorMatchingActive = cms.untracked.bool(True)
151 process.MitTreeFiller.IC5GenJets.flavorMatchingActive = cms.untracked.bool(True)
152 process.MitTreeFiller.SC5GenJets.flavorMatchingActive = cms.untracked.bool(True)
153 process.MitTreeFiller.SC7GenJets.flavorMatchingActive = cms.untracked.bool(True)
154 process.MitTreeFiller.KT4GenJets.flavorMatchingActive = cms.untracked.bool(True)
155 process.MitTreeFiller.KT6GenJets.flavorMatchingActive = cms.untracked.bool(True)
156
157 #produce jet vertex association information
158 process.load("MitProd.TreeFiller.JetVertexAssociation_cfi")
159 #Enable Jet Vertex association for Reco Jet collections
160 process.MitTreeFiller.CaloJets.jetToVertexActive = cms.untracked.bool(True)
161 process.MitTreeFiller.ItrCone5Jets.jetToVertexActive = cms.untracked.bool(True)
162 process.MitTreeFiller.SisCone5Jets.jetToVertexActive = cms.untracked.bool(True)
163 process.MitTreeFiller.SisCone7Jets.jetToVertexActive = cms.untracked.bool(True)
164 process.MitTreeFiller.Kt4Jets.jetToVertexActive = cms.untracked.bool(True)
165 process.MitTreeFiller.Kt6Jets.jetToVertexActive = cms.untracked.bool(True)
166
167 # Path and EndPath definitions
168 process.generation_step = cms.Path(process.pgen)
169 process.simulation_step = cms.Path(process.psim)
170 process.digitisation_step = cms.Path(process.pdigi)
171 process.L1simulation_step = cms.Path(process.SimL1Emulator)
172 process.digi2raw_step = cms.Path(process.DigiToRaw)
173 process.raw2digi_step = cms.Path(process.RawToDigi)
174 process.reconstruction_step = cms.Path(process.reconstruction)
175 process.mitTreeFiller_step = cms.Path(
176 (process.MitEIdSequence
177 + process.MitMetCorrections
178 + process.caloJetMCFlavour
179 + process.jetvertexAssociationSequence
180 + process.ZSPJetCorrections*process.JetPlusTrackCorrections
181 )
182 *process.MitTreeFiller)
183 #process.out_step = cms.EndPath(process.output)
184
185 # Schedule definition
186 process.schedule = cms.Schedule(process.generation_step,process.simulation_step,process.digitisation_step,process.L1simulation_step,process.digi2raw_step,process.raw2digi_step,process.reconstruction_step,process.mitTreeFiller_step)