$ Copied from VAX on 31-mar-01. dat on VAX was 14-mar-98
$! This Rob's com file modified to run at LBL. I am running it on PYTHIA V5.6
$! not V5.7 as the com file says. Also the EXE is called V5_7, but it is
$! V5_6, as I understand from Rob Roser
$!
$!============================================================================
$!  Com file to run PYTHIA top production
$!
$!  Rob Roser/Mark Shaw  25-May-1995
$!
$!============================================================================
$!
$ set noverify
$ setup offline version_7_12 
$!  don't worry about this setup, it just defines logicals.  5.7 does not exist
$!  on this alpha cluster
$ setup pythia v5_6 all
$ setup generator development all
$ setup qq v9_0 all
$!
$ define here monte$her
$! define here1 USR$I2:[ROSER.MC.random]
$ define here1 csa6$programs
$ define data topm$gen
$ set def here         
$! define CDFPYTC USR$I2:[ROSER.MC]CVTABLE.DBT;1 this is for V5_7
$!
$! -----------------------------------------------------------------
$!  140, 150, 160, 170, 180, 190, 200
$!   x    x    x    x    x    x    x
$!
$! define output file, and input/output random number files.
$! ======================================================================
$! ======================================================================
$  define outfil1     data:Ptt_56N_175_21.genq      
$  define outfil2     data:Ptt_56N_175_22.genq      
$  define outfil3     data:Ptt_56N_175_23.genq      
$!  define outfil4     data:Ptt_57_175_04.genq      
$! -----------------------------------------------------------------
$     define RNSEED   here1:top_rnstat_16.dat 
$     define RNSTAT   here1:top_rnstat_116.dat
$! -----------------------------------------------------------------
$! -----------------------------------------------------------------
$!
$ set verify
$! Rob run with V5_7 by Rob, 
$ run/nodebug csa6$programs:PYTHIA.EXE              
!
  in mod pythia
  show input
  use cleomc tauola partfilt genq 
  sho path
!

Talk pythia
    run_number 21                               !*** operator change

    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 2 	! decay t->bW
	mstp 61 1 	! initial state radiation on, 0=off
	mstp 71 1	! final state radiation on
	mstp 52 2 	! MRSD0' S.F.
	mstp 51 330 	! MRSD0' S.F.
!
!       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.006	! 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 outfil1
 BEGIN/rep=500/good=10

!  ======================================================================
!  ======================================================================

Talk pythia
    run_number 22                               !*** operator change
return
!
 output file outfil2
 BEGIN/rep=500/good=10
!!  ======================================================================
!!  ======================================================================
!
Talk pythia
    run_number 23                               !*** operator change
!return
!
 output file outfil3
 BEGIN/rep=500/good=10
!!  ======================================================================
!!  ======================================================================
!
!Talk pythia
!    run_number 4                               !*** operator change
!return
!
! output file outfil4
! BEGIN/rep=500/good=5000
!
 SHO TIMING 
!
 EXIT                                         
$ set noverify