source = PythiaSource 
{ 
  untracked int32 pythiaPylistVerbosity = 0
  untracked bool pythiaHepMCVerbosity = false
  untracked int32 maxEventsToPrint = 0

  PSet PythiaParameters = {

    # This is a vector of ParameterSet names to be read, in this order
    # The first two are in the include files below
    # The last one are simply my additional parameters
    vstring parameterSets = {
      "pythiaUESettings", 
      "processParameters" 
    }

        include 
"Configuration/Spring08Production/data/PythiaUESettings.cfi"

        vstring processParameters = {
      'MSTP(1) = 4',
      'MSEL = 7',
      'PMAS(5,1) = 4.4D0    ! b quarks mass',
      'PMAS(6,1) = 175.0D0  ! t quarks mass',
      'PMAS(7,1) = 300.0D0  ! bprime quarks mass',
      'PMAS(7,2) = 1.0D0', 
      'PMAS(7,3) = 10.0D0',
      'VCKM(1,4) = 0.00078961D0', 
      'VCKM(2,4) = 0.01354896D0',
      'VCKM(3,4) = 0.04700224D0',
      'VCKM(4,4) = 0.93857344D0',
      'VCKM(4,1) = 0.00001936D0',
      'VCKM(4,2) = 0.01290496D0',
      'VCKM(4,3) = 0.04840000D0',
      'MDME(56,1)=0     ! g b4',
      'MDME(57,1)=0     ! gamma b4',
      'MDME(58,1)=0     ! Z0 b4',
      'MDME(59,1)=0     ! W u',
      'MDME(60,1)=0     ! W c',
      'MDME(61,1)=1     ! W t',
      'MDME(62,1)=0     ! W t4',
      'MDME(63,1)=0     ! h0 b4',
      'MDME(64,1)=-1    ! H- c',
      'MDME(65,1)=-1    ! H- t',
      'BRAT(56)  = 0.0D0',
      'BRAT(57)  = 0.0D0',
      'BRAT(58)  = 0.0D0',
      'BRAT(59)  = 0.0D0',
      'BRAT(60)  = 0.0D0',
      'BRAT(61)  = 1.0D0',
      'BRAT(62)  = 0.0D0',
      'BRAT(63)  = 0.0D0',
      'BRAT(64)  = 0.0D0',
      'BRAT(65)  = 0.0D0',
      'MDME(174,1)=1           !Z decay into d dbar',
      'MDME(175,1)=1           !Z decay into u ubar',
      'MDME(176,1)=1           !Z decay into s sbar',
      'MDME(177,1)=1           !Z decay into c cbar',
      'MDME(178,1)=1           !Z decay into b bbar',
      'MDME(179,1)=1           !Z decay into t tbar',
      'MDME(180,1)=-1          !Z decay into b4 b4bar',
      'MDME(181,1)=-1          !Z decay into t4 t4bar',
      'MDME(182,1)=1           !Z decay into e- e+',
      'MDME(183,1)=1           !Z decay into nu_e nu_ebar',
      'MDME(184,1)=1           !Z decay into mu- mu+',
      'MDME(185,1)=1           !Z decay into nu_mu nu_mubar',
      'MDME(186,1)=1           !Z decay into tau- tau+',
      'MDME(187,1)=1           !Z decay into nu_tau nu_taubar',
      'MDME(188,1)=-1          !Z decay into tau4 tau4bar',
      'MDME(189,1)=-1          !Z decay into nu_tau4 nu_tau4bar',
      'MDME(190,1)=1           !W decay into u dbar',
      'MDME(191,1)=1           !W decay into c dbar',
      'MDME(192,1)=1           !W decay into t dbar',
      'MDME(193,1)=-1          !W decay into t4 dbar',
      'MDME(194,1)=1           !W decay into u sbar',
      'MDME(195,1)=1           !W decay into c sbar',
      'MDME(196,1)=1           !W decay into t sbar',
      'MDME(197,1)=-1          !W decay into t4 sbar',
      'MDME(198,1)=1           !W decay into u bbar',
      'MDME(199,1)=1           !W decay into c bbar',
      'MDME(200,1)=1           !W decay into t bbar',
      'MDME(201,1)=-1          !W decay into t4 bbar',
      'MDME(202,1)=-1          !W decay into u b4bar',
      'MDME(203,1)=-1          !W decay into c b4bar',
      'MDME(204,1)=-1          !W decay into t b4bar',
      'MDME(205,1)=-1          !W decay into t4 b4bar',
      'MDME(206,1)=1           !W decay into e- nu_e',
      'MDME(207,1)=1           !W decay into mu nu_mu',
      'MDME(208,1)=1           !W decay into tau nu_tau',
      'MDME(209,1)=-1          !W decay into tau4 nu_tau4'   
        }

  }
}