ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/pgunnell/rivet_cfg.py
Revision: 1.1
Committed: Wed Jan 11 10:14:23 2012 UTC (13 years, 3 months ago) by pgunnell
Content type: text/x-python
Branch: MAIN
CVS Tags: HEAD
Log Message:
MPI_on

File Contents

# User Rev Content
1 pgunnell 1.1 # Auto generated configuration file
2     # using:
3     # Revision: 1.341
4     # Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
5     # with command line options: Configuration/GenProduction/python/MinBias_TuneZ2_7TeV_pythia6_cff.py -s GEN --datatier=GEN-SIM-RAW --conditions auto:mc --eventcontent RAWSIM --no_exec -n 10000 --python_filename=rivet_cfg.py --customise=Configuration/GenProduction/rivet_customize.py
6     import FWCore.ParameterSet.Config as cms
7    
8     process = cms.Process('GEN')
9    
10     # import of standard configurations
11     process.load('Configuration.StandardSequences.Services_cff')
12     process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
13     process.load('FWCore.MessageService.MessageLogger_cfi')
14     process.load('Configuration.EventContent.EventContent_cff')
15     process.load('SimGeneral.MixingModule.mixNoPU_cfi')
16     process.load('Configuration.StandardSequences.GeometryDB_cff')
17     process.load('Configuration.StandardSequences.MagneticField_38T_cff')
18     process.load('Configuration.StandardSequences.Generator_cff')
19     process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic7TeV2011Collision_cfi')
20     process.load('GeneratorInterface.Core.genFilterSummary_cff')
21     process.load('Configuration.StandardSequences.EndOfProcess_cff')
22     process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
23    
24     process.maxEvents = cms.untracked.PSet(
25     input = cms.untracked.int32(600000)
26     # output = cms.untracked.int32(1)
27     )
28    
29     # Input source
30     process.source = cms.Source("EmptySource")
31    
32     process.options = cms.untracked.PSet(
33    
34     )
35    
36     # Production Info
37     process.configurationMetadata = cms.untracked.PSet(
38     version = cms.untracked.string('$Revision: 1.3 $'),
39     annotation = cms.untracked.string('PYTHIA6-MinBias TuneZ2 at 7TeV'),
40     name = cms.untracked.string('$Source: /cvs/CMSSW/CMSSW/Configuration/GenProduction/python/MinBias_TuneZ2_7TeV_pythia6_cff.py,v $')
41     )
42    
43     # Output definition
44    
45     process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
46     splitLevel = cms.untracked.int32(0),
47     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
48     outputCommands = process.RAWSIMEventContent.outputCommands,
49     fileName = cms.untracked.string('MinBias_TuneZ2_7TeV_pythia6_cff_py_GEN.root'),
50     dataset = cms.untracked.PSet(
51     filterName = cms.untracked.string(''),
52     dataTier = cms.untracked.string('GEN-SIM-RAW')
53     ),
54     SelectEvents = cms.untracked.PSet(
55     SelectEvents = cms.vstring('generation_step')
56     )
57     )
58    
59     # Additional output definition
60    
61     # Other statements
62     process.GlobalTag.globaltag = 'MC_44_V5::All'
63    
64     process.generator = cms.EDFilter("Pythia6GeneratorFilter",
65     pythiaPylistVerbosity = cms.untracked.int32(1),
66     filterEfficiency = cms.untracked.double(1.0),
67     pythiaHepMCVerbosity = cms.untracked.bool(False),
68     comEnergy = cms.double(7000.0),
69     crossSection = cms.untracked.double(71260000000.0),
70     maxEventsToPrint = cms.untracked.int32(0),
71     PythiaParameters = cms.PSet(
72     pythiaUESettings = cms.vstring('MSTU(21)=1 ! Check on possible errors during program execution',
73     'MSTJ(22)=2 ! Decay those unstable particles',
74     'PARJ(71)=10 . ! for which ctau 10 mm',
75     'MSTP(33)=0 ! no K factors in hard cross sections',
76     'MSTP(2)=1 ! which order running alphaS',
77     'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
78     'MSTP(52)=2 ! work with LHAPDF',
79     'PARP(82)=1.832 ! pt cutoff for multiparton interactions',
80     'PARP(89)=1800. ! sqrts for which PARP82 is set',
81     'PARP(90)=0.275 ! Multiple interactions: rescaling power',
82     'MSTP(95)=6 ! CR (color reconnection parameters)',
83     'PARP(77)=1.016 ! CR',
84     'PARP(78)=0.538 ! CR',
85     'PARP(80)=0.1 ! Prob. colored parton from BBR',
86     'PARP(83)=0.356 ! Multiple interactions: matter distribution parameter',
87     'PARP(84)=0.651 ! Multiple interactions: matter distribution parameter',
88     'PARP(62)=1.025 ! ISR cutoff',
89     'MSTP(91)=1 ! Gaussian primordial kT',
90     'PARP(93)=10.0 ! primordial kT-max',
91     'MSTP(81)=21 ! multiple parton interactions 1 is Pythia default',
92     'MSTP(82)=4 ! Defines the multi-parton model'),
93     processParameters = cms.vstring('MSEL=5 ! User defined processes',
94     'CKIN(3)=20 ! Minimum value of pthat',
95     #'MSUB(11)=1 ! Min bias process',
96     #'MSUB(12)=1 ! Min bias process',
97     #'MSUB(13)=1 ! Min bias process',
98     #'MSUB(28)=1 ! Min bias process',
99     #'MSUB(53)=1 ! Min bias process',
100     #'MSUB(68)=1 ! Min bias process',
101     #'MSUB(92)=1 ! Min bias process, single diffractive',
102     #'MSUB(93)=1 ! Min bias process, single diffractive',
103     #'MSUB(94)=1 ! Min bias process, double diffractive',
104     #'MSUB(95)=1 ! Min bias process'
105     ),
106     parameterSets = cms.vstring('pythiaUESettings',
107     'processParameters')
108     )
109     )
110    
111     #process.filter12 = cms.EDFilter("MCParticlePairFilter",
112     # status = cms.untracked.vint32(1,1),
113     # MinPt = cms.untracked.vdouble(2.5,2.5),
114     # MaxEta = cms.untracked.vdouble(2.5,2.5),
115     # MinEta = cms.untracked.vdouble(-2.5,-2.5),
116     # ParticleCharge = cms.untracked.int32(-1),
117     # ParticleID1 = cms.untracked.vint32(-5,5),
118     # ParticleID2 = cms.untracked.vint32(5,-5)
119     #)
120    
121     process.ProductionFilterSequence = cms.Sequence(process.generator)
122    
123    
124     #process.ProductionFilterSequence = cms.Sequence(process.generator)
125    
126     # Path and EndPath definitions
127     process.generation_step = cms.Path(process.pgen)
128     process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
129     process.endjob_step = cms.EndPath(process.endOfProcess)
130     process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
131    
132     # Schedule definition
133     process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.endjob_step,process.RAWSIMoutput_step)
134     # filter all path with the production filter sequence
135     for path in process.paths:
136     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
137    
138     # customisation of the process.
139    
140     # Automatic addition of the customisation function from Configuration.GenProduction.rivet_customize
141     process.load('GeneratorInterface.RivetInterface.rivetAnalyzer_cfi')
142    
143     from Configuration.GenProduction.rivet_customize import customise
144    
145     # My analyzer
146     process.rivetAnalyzer.AnalysisNames = cms.vstring('CMS_2011_SIM')
147     process.rivetAnalyzer.OutputFile = cms.string('mcfile.aida')
148     process.rivetAnalyzer.CrossSection = cms.double(1133000)
149    
150     #call to customisation function customise imported from Configuration.GenProduction.rivet_customize
151     process = customise(process)
152    
153     # End of customisation functions