my setting for Drell Yan simulation

 

   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
   }

DY setting:  DY setting: