ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/TreeFiller/python/FullChainExample.py
Revision: 1.3
Committed: Sun Mar 22 08:48:59 2009 UTC (16 years, 1 month ago) by loizides
Content type: text/x-python
Branch: MAIN
CVS Tags: HEAD
Changes since 1.2: +1 -1 lines
State: FILE REMOVED
Log Message:
Cleanup. Moved MitProd dependent python files in MitEdm to MitProd. Removed old set of fillers. Intention is to keep three top level MitTreeFiller_template files for RECO, AOD and CRAFT data.

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.2 $'),
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/UserCode/MitProd/TreeFiller/python/FullChainExample.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 #Load Mit vProducer
120 process.load("MitEdm.Producers.vProducer_cff")
121
122 #Load Mit Mvf Conversion producer
123 process.load("MitEdm.Producers.conversionProducer_cff")
124
125 # compute ECAL shower shape variables
126 process.load("Configuration.StandardSequences.Geometry_cff")
127 process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
128 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
129 process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
130
131 process.load("FWCore.MessageLogger.MessageLogger_cfi")
132
133 #Load ElectronID information
134 process.load("MitProd.TreeFiller.ElectronID_cfi")
135 #Load Electron Isolation. For Jurassic isolation veto fix.
136 process.load("RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff")
137
138
139 #For Jet Corrections (Summer08 Jet corrections)
140 process.load("JetMETCorrections.Configuration.L2L3Corrections_Summer08_cff")
141 process.prefer("L3JetCorrectorIC5Calo")
142
143 #enable Jet Corrections for all of our Jet collections
144 process.MitTreeFiller.CaloJets.jetCorrectionsActive = cms.untracked.bool(True)
145 process.MitTreeFiller.ItrCone5Jets.jetCorrectionsActive = cms.untracked.bool(True)
146 process.MitTreeFiller.SisCone5Jets.jetCorrectionsActive = cms.untracked.bool(True)
147 process.MitTreeFiller.SisCone7Jets.jetCorrectionsActive = cms.untracked.bool(True)
148 process.MitTreeFiller.Kt4Jets.jetCorrectionsActive = cms.untracked.bool(True)
149 process.MitTreeFiller.Kt6Jets.jetCorrectionsActive = cms.untracked.bool(True)
150
151 #For JetPlusTracks
152 process.load("JetMETCorrections.Configuration.JetPlusTrackCorrections_cff")
153 process.load("JetMETCorrections.Configuration.ZSPJetCorrections152_cff")
154 process.MitTreeFiller.IC5JetPlusTrack.active = cms.untracked.bool(True)
155
156 #Load Met Corrections
157 process.load("MitProd.TreeFiller.MetCorrections_cfi")
158
159 #Load Flavor Matching Information
160 process.load("MitProd.TreeFiller.JetsMCFlavourMatching_cfi")
161 #Enable Flavor matching for Reco Jets and GenJets
162 process.MitTreeFiller.CaloJets.flavorMatchingActive = cms.untracked.bool(True)
163 process.MitTreeFiller.ItrCone5Jets.flavorMatchingActive = cms.untracked.bool(True)
164 process.MitTreeFiller.SisCone5Jets.flavorMatchingActive = cms.untracked.bool(True)
165 process.MitTreeFiller.SisCone7Jets.flavorMatchingActive = cms.untracked.bool(True)
166 process.MitTreeFiller.Kt4Jets.flavorMatchingActive = cms.untracked.bool(True)
167 process.MitTreeFiller.Kt6Jets.flavorMatchingActive = cms.untracked.bool(True)
168 process.MitTreeFiller.IC5GenJets.flavorMatchingActive = cms.untracked.bool(True)
169 process.MitTreeFiller.SC5GenJets.flavorMatchingActive = cms.untracked.bool(True)
170 process.MitTreeFiller.SC7GenJets.flavorMatchingActive = cms.untracked.bool(True)
171 process.MitTreeFiller.KT4GenJets.flavorMatchingActive = cms.untracked.bool(True)
172 process.MitTreeFiller.KT6GenJets.flavorMatchingActive = cms.untracked.bool(True)
173
174 #produce jet vertex association information
175 process.load("MitProd.TreeFiller.JetVertexAssociation_cfi")
176 #Enable Jet Vertex association for Reco Jet collections
177 process.MitTreeFiller.CaloJets.jetToVertexActive = cms.untracked.bool(True)
178 process.MitTreeFiller.ItrCone5Jets.jetToVertexActive = cms.untracked.bool(True)
179 process.MitTreeFiller.SisCone5Jets.jetToVertexActive = cms.untracked.bool(True)
180 process.MitTreeFiller.SisCone7Jets.jetToVertexActive = cms.untracked.bool(True)
181 process.MitTreeFiller.Kt4Jets.jetToVertexActive = cms.untracked.bool(True)
182 process.MitTreeFiller.Kt6Jets.jetToVertexActive = cms.untracked.bool(True)
183
184 #Load track detector associator for Track-ECal association
185 process.load("MitProd.TreeFiller.TrackEcalAssociation_cfi")
186 process.MitTreeFiller.TrackAssociatorParameters = cms.untracked.PSet(process.TrackAssociatorParameters)
187 #enable Track-Ecal assocation in fillers
188 process.MitTreeFiller.GeneralTracks.ecalAssocActive = True
189 process.MitTreeFiller.StandaloneMuonTracks.ecalAssocActive = True
190 process.MitTreeFiller.StandaloneMuonTracksWVtxConstraint.ecalAssocActive = True
191 process.MitTreeFiller.GlobalMuonTracks.ecalAssocActive = True
192 process.MitTreeFiller.ConversionInOutTracks.ecalAssocActive = True
193 process.MitTreeFiller.ConversionOutInTracks.ecalAssocActive = True
194 process.MitTreeFiller.GsfTracks.ecalAssocActive = True
195
196 #load gsf track to general track associator
197 process.load("MitEdm.Producers.gsfTrackAssociator_cff")
198 process.MitTreeFiller.Electrons.gsfTrackAssocName = 'gsfTrackAssociator'
199
200 # Path and EndPath definitions
201 process.generation_step = cms.Path(process.pgen)
202 process.simulation_step = cms.Path(process.psim)
203 process.digitisation_step = cms.Path(process.pdigi)
204 process.L1simulation_step = cms.Path(process.SimL1Emulator)
205 process.digi2raw_step = cms.Path(process.DigiToRaw)
206 process.raw2digi_step = cms.Path(process.RawToDigi)
207 process.reconstruction_step = cms.Path(process.reconstruction)
208 process.mitTreeFiller_step = cms.Path(
209 process.eleIsolationSequence *
210 process.gsfTrackAssociator *
211 process.vProducer *
212 process.conversionProducer *
213 ( process.MitEIdSequence
214 + process.MitMetCorrections
215 + process.caloJetMCFlavour
216 + process.jetvertexAssociationSequence
217 + process.ZSPJetCorrections*process.JetPlusTrackCorrections
218 )
219 *process.MitTreeFiller
220 *process.vFiller
221 *process.conversionFiller
222 )
223 #process.out_step = cms.EndPath(process.output)
224
225 # Schedule definition
226 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)