Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
import FWCore.ParameterSet.Config as cms
from GeneratorInterface.GenFilters.CosmicGenFilterHelix_cff import * 

generator = cms.EDProducer("CosMuoGenProducer",
    ZCentrOfTarget = cms.double(0.0),
    MinP = cms.double(10.0),
    MinP_CMS = cms.double(-1.0), ##negative means MinP_CMS = MinP. Only change this if you know what you are doing!
    MaxP = cms.double(3000.0),
    MinTheta = cms.double(0.0),
    MaxTheta = cms.double(84.),
    MinPhi = cms.double(0.0),
    MaxPhi = cms.double(360.0),
    MinT0 = cms.double(-12.5),
    MaxT0 = cms.double(12.5),
    PlugVx = cms.double(0.0),
    PlugVz = cms.double(-14000.0),                
    MinEnu = cms.double(10.),                
    MaxEnu = cms.double(10000.),                
    NuProdAlt = cms.double(7.5e6),                       
    AcptAllMu = cms.bool(False), 
    ElossScaleFactor = cms.double(1.0),
    RadiusOfTarget = cms.double(8000.0),
    ZDistOfTarget = cms.double(15000.0),
    TrackerOnly = cms.bool(False),
    TIFOnly_constant = cms.bool(False),
    TIFOnly_linear = cms.bool(False),
    MTCCHalf = cms.bool(False),
    Verbosity = cms.bool(False),
    RhoAir = cms.double(0.001214),
    RhoWall = cms.double(2.5),
    RhoRock = cms.double(2.5),
    RhoClay = cms.double(2.3),
    RhoPlug = cms.double(2.5),
    ClayWidth = cms.double(50000.),
                           
    MultiMuon = cms.bool(False),
    # MultiMuon = cms.bool(True),
    MultiMuonFileName = cms.string("CORSIKAmultiMuon.root"),
    MultiMuonFileFirstEvent = cms.int32(1),
    MultiMuonNmin = cms.int32(2),              
)


#Filter
ProductionFilterSequence = cms.Sequence(generator*cosmicInPixelLoose)