python/test/template_pion_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("CastorTest")
0004 process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0005
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007
0008
0009 #process.load("Configuration.StandardSequences.GeometryExtended_cff")
0010 process.load("SimG4CMS.Forward.castorGeometryXML_cfi")
0011 process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
0012 process.load("Configuration.StandardSequences.MagneticField_cff")
0013 process.load("SimG4Core.Application.g4SimHits_cfi")
0014 process.load("Configuration.StandardSequences.Generator_cff")
0015
0016
0017 process.RandomNumberGeneratorService.theSource.initialSeed = RAND_sourceSeed
0018 process.RandomNumberGeneratorService.generator.initialSeed = RAND_generator
0019 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = RAND_VtxSmeared
0020 process.RandomNumberGeneratorService.g4SimHits.initialSeed = RAND_g4SimHits
0021 process.RandomNumberGeneratorService.mix.initialSeed = RAND_mix
0022
0023
0024 process.maxEvents = cms.untracked.PSet(
0025     input = cms.untracked.int32(MAXEVENTS)
0026 )
0027
0028
0029 process.source = cms.Source("EmptySource")
0030
0031 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0032     PGunParameters = cms.PSet(
0033         PartID = cms.vint32(211),
0034         MaxEta = cms.double(-5.9),
0035         MaxPhi = cms.double(3.14),
0036         MinEta = cms.double(-6.0),
0037         MinE = cms.double(MINIMUM_ENERGY),
0038         MinPhi = cms.double(-3.14),
0039         MaxE = cms.double(MAXIMUM_ENERGY)
0040     ),
0041     AddAntiParticle = cms.bool(True),
0042     Verbosity = cms.untracked.int32(0)
0043 )
0044
0045 process.CaloSD = cms.PSet(
0046     NeutronThreshold = cms.double(30.0),
0047     ProtonThreshold  = cms.double(30.0),
0048     IonThreshold     = cms.double(30.0),
0049     SuppressHeavy    = cms.bool(False),
0050     EminTrack        = cms.double(1.0),
0051     TmaxHit          = cms.double(1000.0),
0052     HCNames          = cms.vstring('ZDCHITS'),
0053     EminHits         = cms.vdouble(0.0),
0054     TmaxHits         = cms.vdouble(2000.0),
0055     UseResponseTables= cms.vint32(0),
0056     BeamPosition     = cms.double(0.0),
0057     CorrectTOFBeam   = cms.bool(False),
0058     DetailedTiming   = cms.untracked.bool(False),
0059     UseMap           = cms.untracked.bool(True),
0060     Verbosity        = cms.untracked.int32(0),
0061     CheckHits        = cms.untracked.int32(25)
0062 )
0063 process.o1 = cms.OutputModule("PoolOutputModule",
0064     fileName = cms.untracked.string('POOLOUTPUTFILE')
0065 )
0066
0067
0068 process.common_maximum_timex = cms.PSet( # need to be localy redefined
0069    MaxTrackTime  = cms.double(500.0),  # need to be localy redefined
0070    MaxTimeNames  = cms.vstring(), # need to be localy redefined
0071    MaxTrackTimes = cms.vdouble()  # need to be localy redefined
0072 )
0073 process.p1 = cms.Path(process.generator*process.pgen*process.VtxSmeared*process.g4SimHits)
0074 process.outpath = cms.EndPath(process.o1)
0075 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP'
0076 process.g4SimHits.Generator.ApplyEtaCuts = False
0077
0078
0079
0080 process.g4SimHits.StackingAction = cms.PSet(
0081    process.common_heavy_suppression,
0082    process.common_maximum_timex,        # need to be localy redefined
0083    TrackNeutrino = cms.bool(False),
0084    KillDeltaRay  = cms.bool(False),
0085    KillHeavy     = cms.bool(False),
0086    SaveFirstLevelSecondary = cms.untracked.bool(True),
0087    SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0088    SavePrimaryDecayProductsAndConversionsInCalo    = cms.untracked.bool(True),
0089   SavePrimaryDecayProductsAndConversionsInMuon    = cms.untracked.bool(True)
0090 )
0091
0092 process.g4SimHits.SteppingAction = cms.PSet(
0093    process.common_maximum_timex, # need to be localy redefined
0094    KillBeamPipe            = cms.bool(True),
0095    CriticalEnergyForVacuum = cms.double(2.0),
0096    CriticalDensity         = cms.double(1e-15),
0097    EkinNames               = cms.vstring(),
0098    EkinThresholds          = cms.vdouble(),
0099    EkinParticles           = cms.vstring(),
0100    Verbosity               = cms.untracked.int32(0)
0101 )
0102
0103 process.g4SimHits.G4Commands = ['/tracking/verbose 1']
0104 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0105     type = cms.string('CastorTestAnalysis'),
0106     CastorTestAnalysis = cms.PSet(
0107         EventNtupleFileName = cms.string('NTUPLEFILE'),
0108         Verbosity = cms.int32(0),
0109         StepNtupleFileName = cms.string('stepNtuple.root'),
0110         StepNtupleFlag = cms.int32(0),
0111         EventNtupleFlag = cms.int32(1)
0112     )
0113 ))
0114