source = PythiaSource 
{ 
   // untracked int32 pythiaPylistVerbosity = 1
   // untracked bool pythiaHepMCVerbosity = false
   // untracked int32 maxEventsToPrint = 2
  PSet PythiaParameters = {
	

    # This is a vector of ParameterSet names to be read, in this order
    vstring parameterSets = {   
      "processParameters"
    }

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

    vstring processParameters = {
        "MSEL=0                  !User defined processes",
        "MSUB(1)=1               !Incl Z0/gamma* production",
        "MSTP(43)=2              !Only Z0", 
        "MDME(174,1)=0           !Z decay into d dbar",        
        "MDME(175,1)=0           !Z decay into u ubar",
        "MDME(176,1)=0           !Z decay into s sbar",
        "MDME(177,1)=0           !Z decay into c cbar",
        "MDME(178,1)=0           !Z decay into b bbar",
        "MDME(179,1)=0           !Z decay into t tbar",
        "MDME(182,1)=0           !Z decay into e- e+",
        "MDME(183,1)=0           !Z decay into nu_e nu_ebar",
        "MDME(184,1)=1           !Z decay into mu- mu+",
        "MDME(185,1)=0           !Z decay into nu_mu nu_mubar",
        "MDME(186,1)=0           !Z decay into tau- tau+",
        "MDME(187,1)=0           !Z decay into nu_tau nu_taubar"
	
 	   }
	}	
  
	}