| 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 |
|
|
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)
|
