$Copied from VAX on March-31-01. date on VAX was 22-mar-98
$! this was run on the alpha with PYTHIA V5.6 (EXE made in June 1995)
$! PYTHIA LIBRARY MADE IN MARCH 1994. V5.7 was made in Nov. 95 (EXE in Nov. 95).
$!============================================================================
$!  Com file to run PYTHIA top production
$!
$!  Rob Roser/Mark Shaw  25-May-1995
$!
$!============================================================================
$!
$ set noverify
$ setup offline version_7_12 
$ setup cern.v94b
$ setup pythia v5_6 all
$ setup generator development all
$ setup qq v9_0 all
$!
$! define here USR$I2:[ROSER.MC]
$ define here monte$her
$! define here1 USR$I2:[ROSER.MC.random]
$ define here1 csa6$programs
$! define data DATA1$I1:[ROSER.GEN]
$! define data topm$gen
$ define data csa6$topmc4
$ set def here         
$!
 $! -----------------------------------------------------------------
$!  140, 150, 160, 170, 180, 190, 200
$!                       
$!
$ define outfile   data:Pttbar_175_27.genq      !***operator assistance
$!
$! problem with 20!
$!
$! input random #
$     define RNSEED here1:top_rnstat_27.dat           !operator change
$     define pyranin here1:py_top_rnstat_27.dat       !operator change
$!     define RNSEED here1:top_rnstat_22.dat
$!     define pyranin here1:py_top_rnstat_22.dat
$! output random #
$     define RNSTAT here1:top_rnstat_28.dat           !operator change
$     define pyranout here1:py_top_rnstat_28.dat      !operator change
$! -----------------------------------------------------------------
$!
$ set verify
$ run/nodebug here1:pythia              
!
! ---------------------------------------------------------------------------
  in mod pythia
  show input
!  use cleomc partfilt tauola genq   ! as used by Rob
! invert tauola and partfilt
  use cleomc tauola partfilt genq 
  sho path
!
! pythia settings are as described in cdf2920
! Talk to for PYTHIA to configure for b-bbar production
! ttbar generation
! Minimum Pt of t quark approx 7 GeV (mod init,final state radiation)

Talk pythia
    run_number 27                               !*** operator change

!
! random numbers...
    random_number
      write
        file pyranout
        final on
        show
        return
      read
        file pyranin
        final on
        show
        return
      return


    process
        set 6		! ttbar
	list
        return
    cut
	ckin 3 5.0	! min Pt of generated t quark
	ckin 4 750.0	! max Pt of generated t quark
	show
        return
    pythia_par
	mstp 48 1 	! decay t->bW
	mstp 61 1 	! initial state radiation on, 0=off
!	mstp 71 1	! final state radiation on
!
!       Tuned underlying event parameters
!
        mstp 2  2       ! 2nd order alpha_s (D=1)
        mstp 33 2       ! separate K factors (D=0)
        mstp 82 2       ! multiple interactions (D=1)
        parp 82 1.      ! regularization scale of the pt for Mult Int. (D=1.7)
        return
    jetset_par
	mstj 11 3	! use Peterson fragmention for b, Lund for light q's
	mstj 27 1	! include glon rad
	parj 55 -0.0006	! Peterson fragmentation parameter
        particle
                mass 6 175.     !top mass
                return
	decay           ! Turn off long-lived particle decays
	    stable     6 off	!
	    stable   130 on	! Kl
	    stable   310 on	! Ks
	    stable  3122 on	! Lambda
	    stable -3122 on
	    stable  3112 on	! Sigma-
	    stable -3112 on
	    stable  3222 on 	! Sigma+
	    stable -3222 on
	    stable  3312 on 	! Cascade-
	    stable -3312 on
	    stable  3322 on 	! Cascade0
	    stable -3322 on	 
	    stable -3334  on	! Omega-
	    stable -3334 on
            return
        return
return

set ver
talk partfilt
 add
    particle_name (mu-)
    Pt_min 15.0
    eta_max 3.0
    acc
 add
    particle_name (mu+)
    Pt_min 15.0
    eta_max 3.0
    acc
 add
    particle_name (e+)
    Pt_min 15.0
    eta_max 3.0
    acc
 add
    particle_name (e-)
    Pt_min 15.0
    eta_max 3.0
    acc
 show
 return

! ---------------------------------------------------------------------------
!
! Set Cleo parameters
! ===================
 TALK CLEOMC
	INPUT GENP
	OUTPUT GENP
        BC_DECAY ON                      !Bc
        BB_DECAY ON                      !bottomonia
        B_BARYON_DECAY ON                !b baryons
 RETURN
! ---------------------------------------------------------------------------
!
! Set Tauola parameters
! =====================
 TALK TAUOLA/PAR=2 MODIFY BRANCHING_RATIO
	ELECTRON 1.0
	MUON 0.9800
	PION 0.6496
	KAON 0.0372
	KST  0.1388
	3PIPIZ 0.2660
	PI3PIZ 0.0822
	KKZ  0.0178
 RETURN RETURN RETURN
!
! ---------------------------------------------------------------------------
 Output Select A_Bank On
 output select kept_banks gen*
 output select kept_banks hepg
 output select dropped_banks genp
 FILTER/path=1 partfilt ON
 output select event/filter=partfilt
!
!
 output file outfile
! hist open/file=histfile
! hist/mod=topgenp on 
 BEGIN/rep=500/good=20000
! hist write
 SHO TIMING 
!
 EXIT                                         
$ set noverify