import FWCore.ParameterSet.Config as cms

generator = cms.EDProducer("Pythia6PtGun",
    PGunParameters = cms.PSet(
        ParticleID = cms.vint32(-15),
        AddAntiParticle = cms.bool(True),
        MinPhi = cms.double(-3.14159265359),
        MaxPhi = cms.double(3.14159265359),
        MinPt = cms.double(20.0),
        MaxPt = cms.double(100.0),
        MinEta = cms.double(-2.5),
        MaxEta = cms.double(2.5)
    ),
    PythiaParameters = cms.PSet(
        pythiaTauJets = cms.vstring(
            'MDME(89,1)=0       ! BR=0.178300. tau->nu_ebar,e-,nu_tau',
            'MDME(90,1)=0       ! BR=0.173500. tau->nu_mubar,mu-,nu_tau',
            'MDME(91,1)=0       ! BR=0.113100. tau->nu_tau,pi-',
            'MDME(92,1)=0       ! BR=0.249400. tau->nu_tau,rho-',
            'MDME(93,1)=0       ! BR=0.003000. tau->nu_tau,pi-,pi0',
            'MDME(94,1)=0       ! BR=0.090000. tau->nu_tau,rho-,pi0',
            'MDME(95,1)=0       ! BR=0.002700. tau->nu_tau,pi-,pi0,pi0',
            'MDME(96,1)=0       ! BR=0.010000. tau->nu_tau,rho-,pi0,pi0',
            'MDME(97,1)=0       ! BR=0.001400. tau->nu_tau,pi-,pi0,pi0,pi0',
            'MDME(98,1)=0       ! BR=0.001200. tau->nu_tau,rho-,pi0,pi0,pi0',
            'MDME(99,1)=1       ! BR=0.000250. tau->nu_tau,pi-,K_S0',
            'MDME(100,1)=1      ! BR=0.000250. tau->nu_tau,pi-,K_L0',
            'MDME(101,1)=1      ! BR=0.007100. tau->nu_tau,K-',
            'MDME(102,1)=1      ! BR=0.012000. tau->nu_tau,K*-',
            'MDME(103,1)=1      ! BR=0.000400. tau->nu_tau,K-,pi0',
            'MDME(104,1)=1      ! BR=0.000750. tau->nu_tau,K*-,pi0',
            'MDME(105,1)=1      ! BR=0.000060. tau->nu_tau,K*-,pi0,pi0',
            'MDME(106,1)=1      ! BR=0.000780. tau->nu_tau,K-,K_S0',
            'MDME(107,1)=1      ! BR=0.000780. tau->nu_tau,K-,K_L0',
            'MDME(108,1)=1      ! BR=0.003400. tau->nu_tau,K-,K+,pi-',
            'MDME(109,1)=1      ! BR=0.080000. tau->nu_tau,pi-,rho0',
            'MDME(110,1)=1      ! BR=0.011000. tau->nu_tau,pi-,pi+,pi-',
            'MDME(111,1)=1      ! BR=0.019100. tau->nu_tau,pi-,omega',
            'MDME(112,1)=1      ! BR=0.000060. tau->nu_tau,pi-,eta',
            'MDME(113,1)=1      ! BR=0.005000. tau->nu_tau,rho-,rho0',
            'MDME(114,1)=1      ! BR=0.013300. tau->nu_tau,pi-,rho0,pi0',
            'MDME(115,1)=1      ! BR=0.006700. tau->nu_tau,rho-,pi+,pi-',
            'MDME(116,1)=1      ! BR=0.000500. tau->nu_tau,pi-,pi+,pi-,pi0',
            'MDME(117,1)=1      ! BR=0.003500. tau->nu_tau,rho-,omega',
            'MDME(118,1)=1      ! BR=0.000600. tau->nu_tau,pi-,omega,pi0',
            'MDME(119,1)=1      ! BR=0.001500. tau->nu_tau,rho-,eta',
            'MDME(120,1)=1      ! BR=0.000210. tau->nu_tau,pi-,eta,pi0',
            'MDME(121,1)=1      ! BR=0.000200. tau->nu_tau,rho-,rho0,pi0',
            'MDME(122,1)=1      ! BR=0.000750. tau->nu_tau,pi-,rho0,rho0',
            'MDME(123,1)=1      ! BR=0.000100. tau->nu_tau,pi-,eta,eta',
            'MDME(124,1)=1      ! BR=0.000200. tau->nu_tau,pi-,rho0,pi0,pi0',
            'MDME(125,1)=1      ! BR=0.001100. tau->nu_tau,rho-,rho0,pi0,pi0',
            'MDME(126,1)=1      ! BR=0.000200. tau->nu_tau,pi-,rho+,rho-',
            'MDME(127,1)=1      ! BR=0.000200. tau->nu_tau,pi-,rho+,pi-,pi0',
            'MDME(128,1)=1      ! BR=0.000200. tau->nu_tau,pi-,rho-,pi+,pi0',
            'MDME(129,1)=1      ! BR=0.000220. tau->nu_tau,pi-,rho0,rho0,pi0',
            'MDME(130,1)=1      ! BR=0.000400. tau->nu_tau,K*-,pi0,pi0',
            'MDME(131,1)=1      ! BR=0.000100. tau->nu_tau,K-,pi0,pi0,pi0',
            'MDME(132,1)=1      ! BR=0.002050. tau->nu_tau,pi-,K_S0,pi0',
            'MDME(133,1)=1      ! BR=0.002050. tau->nu_tau,pi-,K_L0,pi0',
            'MDME(134,1)=1      ! BR=0.000690. tau->nu_tau,K-,K_S0,pi0',
            'MDME(135,1)=1      ! BR=0.000690. tau->nu_tau,K-,K_L0,pi0',
            'MDME(136,1)=1      ! BR=0.000250. tau->nu_tau,pi-,K_S0,K_S0',
            'MDME(137,1)=1      ! BR=0.000510. tau->nu_tau,pi-,K_S0,K_L0',
            'MDME(138,1)=1      ! BR=0.000250. tau->nu_tau,pi-,K_L0,K_L0',
            'MDME(139,1)=0      ! BR=0.000000. tau->gamma,tau-',
            'MDME(140,1)=0      ! BR=0.000000. tau->Z0,tau-',
            'MDME(141,1)=0      ! BR=0.000000. tau->W-,nu_tau',
            'MDME(142,1)=0      ! BR=0.000000. tau->h0,tau-'
        ),
        pythiaUESettings = cms.vstring(
            'MSTJ(11)=3     ! Choice of the fragmentation function',
            'MSTJ(22)=2     ! Decay those unstable particles',
            'PARJ(71)=10 .  ! for which ctau  10 mm',
            'MSTP(2)=1      ! which order running alphaS',
            'MSTP(33)=0     ! no K factors in hard cross sections',
            'MSTP(51)=7     ! structure function chosen',
            'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default',
            'MSTP(82)=4     ! Defines the multi-parton model',
            'MSTU(21)=1     ! Check on possible errors during program execution',
            'PARP(82)=1.9409   ! pt cutoff for multiparton interactions',
            'PARP(89)=1960. ! sqrts for which PARP82 is set',
            'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter',
            'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter',
            'PARP(90)=0.16  ! Multiple interactions: rescaling power',
            'PARP(67)=2.5    ! amount of initial-state radiation',
            'PARP(85)=1.0  ! gluon prod. mechanism in MI',
            'PARP(86)=1.0  ! gluon prod. mechanism in MI',
            'PARP(62)=1.25   ! ',
            'PARP(64)=0.2    ! ',
            'MSTP(91)=1     !',
            'PARP(91)=2.1   ! kt distribution',
            'PARP(93)=15.0  ! '
        ),
        parameterSets = cms.vstring(
            'pythiaUESettings',
            'pythiaTauJets'
        )
    )
)


