ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitProd/TreeFiller/python/FullChainExample.py
Revision: 1.2
Committed: Mon Jan 12 11:57:17 2009 UTC (16 years, 4 months ago) by bendavid
Content type: text/x-python
Branch: MAIN
CVS Tags: Mit_008pre2, Mit_008pre1, Mit_006b, Mit_006a
Changes since 1.1: +47 -7 lines
Log Message:
brought full example up to date with latest framework changes

File Contents

# User Rev Content
1 bendavid 1.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 bendavid 1.2 version = cms.untracked.string('$Revision: 1.1 $'),
30 bendavid 1.1 annotation = cms.untracked.string('Z/gamma* to mumu sample, M(mumu) gt40 GeV -- sqrt(s) = 10TeV'),
31 bendavid 1.2 name = cms.untracked.string('$Source: /cvs_server/repositories/CMSSW/UserCode/MitProd/TreeFiller/python/FullChainExample.py,v $')
32 bendavid 1.1 )
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 bendavid 1.2 #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 bendavid 1.1
133     #Load ElectronID information
134     process.load("MitProd.TreeFiller.ElectronID_cfi")
135 bendavid 1.2 #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 bendavid 1.1
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 bendavid 1.2 #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 bendavid 1.1 # 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 bendavid 1.2 process.eleIsolationSequence *
210     process.gsfTrackAssociator *
211     process.vProducer *
212     process.conversionProducer *
213     ( process.MitEIdSequence
214 bendavid 1.1 + process.MitMetCorrections
215     + process.caloJetMCFlavour
216     + process.jetvertexAssociationSequence
217     + process.ZSPJetCorrections*process.JetPlusTrackCorrections
218     )
219 bendavid 1.2 *process.MitTreeFiller
220     *process.vFiller
221     *process.conversionFiller
222     )
223 bendavid 1.1 #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)