if ( mode == "pp:DY" )

   {

     // to filter DY(e-/e+), for some reason, this method doesn't work, use starfiltermaker instead

     //command("gfilt dyFilter");

     pythia6->SetFrame("CMS", 500.0 );

     pythia6->SetBlue("proton");

     pythia6->SetYell("proton");

     if ( tune )

     {

         pythia6->PyTune( tune );

         PyPars_t &pypars = pythia6->pypars();

         printf("PARP(90) was %f, replacing it with 0.213",pypars.parp(90));

         pypars.parp(90)=0.213;

      }

     // Setup pythia process

     PySubs_t &pysubs = pythia6->pysubs();

     PyDat3_t &pydat3 = pythia6->pydat3();

     pysubs.msel = 0;

     pysubs.msub(1)=1; //qq->ll & qq->Z

     // set limits of allowed pseudorapidity, parton level

     //pysubs.ckin(13)=2.0;

     //pysubs.ckin(14)=5.0;

     // set limits of allowed pseudorapidity, parton level

     //pysubs.ckin(15)=2.0;

     //pysubs.ckin(16)=5.0;

 

     // Turn off other γ*/z0 decay channels to speed up DY->e-/e+ pair,

     pydat3.mdme(174,1)=0; // d      + dbar,        off

     pydat3.mdme(175,1)=0; // u      + ubar,        off

     pydat3.mdme(176,1)=0; // s      + sbar,        off

     pydat3.mdme(177,1)=0; // c      + cbar,        off

     pydat3.mdme(178,1)=0; // b      + bbar,        off

     pydat3.mdme(179,1)=0; // t      + tbar,        off

     pydat3.mdme(182,1)=1; // e-     + e+,          on, allow decays to e-/e+

     pydat3.mdme(183,1)=0; // nu_e   + nu_ebar,     off

     pydat3.mdme(184,1)=0; // mu-    + mu+,         off

     pydat3.mdme(185,1)=0; // nu_mu  + nu_mubar,    off

     pydat3.mdme(186,1)=0; // tau-   + tau+,        off

     pydat3.mdme(187,1)=0; // nu_tau + nu_taubar,   off

     pydat3.mdme(180,1)=0; // not developed channel, cross section is 0

     pydat3.mdme(181,1)=0; // not developed channel, cross section is 0

     pydat3.mdme(188,1)=0; // not developed channel, cross section is 0

     pydat3.mdme(189,1)=0; // not developed channel, cross section is 0

     //pydat3.mdcy(102,1)=1; // set pi0(id=102) decay

     //pydat3.mdcy(pycomp(111),1)=1; // set pi0(id=102) decay

     //pythia6 -> SetDecayFlag( 111, 1 ); // pi0

     //pysubs.ckin(3)=1.0; // parton level, lower limit momentum, GeV

   }