import FWCore.ParameterSet.Config as cms

generator = cms.EDProducer("CosMuoGenProducer",
    RadiusOfTarget = cms.double(8000.0),
    TIFOnly_constant = cms.bool(False),
    TIFOnly_linear = cms.bool(False),
    ElossScaleFactor = cms.double(1.0),
    MaxPhi = cms.double(360.0),
    MaxTheta = cms.double(84.26),
    #MaxTheta = cms.double(89.0),
    MaxT0 = cms.double(12.5),
    ZDistOfTarget = cms.double(15000.0),
    ZCentrOfTarget = cms.double(0.0),
    MinP = cms.double(3.0),
    MinT0 = cms.double(-12.5),
    MTCCHalf = cms.bool(False),
    TrackerOnly = cms.bool(False),
    MinTheta = cms.double(0.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),
    MinPhi = cms.double(0.0),
    PlugVx = cms.double(0.0),
    PlugVz = cms.double(-14000.0), #[mm]

    #material densities in g/cm^3
    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.), #[mm]

   MultiMuon = cms.bool(False),
    # MultiMuon = cms.bool(True),
    MultiMuonFileName = cms.string("CORSIKAmultiMuon.root"),
    MultiMuonFileFirstEvent = cms.int32(1),
    MultiMuonNmin = cms.int32(2),                
                           
    MinEnu = cms.double(10.),                
    MaxEnu = cms.double(10000.),                
    NuProdAlt = cms.double(7.5e6),                       
    AcptAllMu = cms.bool(False),                
    Verbosity = cms.bool(False)
)
