#!/bin/tcsh
cd .
bbh<<quit 
!show routines
!file go
 include 'ludatr.f'
 include 'bbhcom.f'
c
c
cmrlu(1)=19xxx503 ! ran # seed (0.le.mrlu(1).le.900000000, def=19780503)
 mrlu(1)=19011503 ! ran # seed (0.le.mrlu(1).le.900000000, def=19780503)
c
 mrlu(2)=0
c
 i_lepton_bbh=11 ! =0 hadronc prod;  =11,13,15 for e,mu,tau pairs resp.
c
ci_stdxwrt=0     !  PAW ntuple output format (HEPEVT common contents)
 i_stdxwrt=1     !  stdHEP output format (HEPEVT common contents)
c
c  datfile_bbh is the stdHEP (i_stdxwrt=1) or PAW ntuple (i_stdxwrt=0) output file
c
 datfile_bbh='bbhxxx_stage'
c
c
c
c
cc c_pzr.... are parton distribution files
c
 c_pzrn_inp_bbh='/nfs/juno/work/timb/bbh009_pzrn'
c
 c_pzrdat_inp_bbh='/nfs/juno/work/timb/bbh009_pzrdat'
c
 c_pzrdat2_inp_bbh='/nfs/juno/work/timb/bbh009_pzrdat2'
c
c
c
c
 c_pzrn_out_bbh='bbhxxx_pzrn'
c
 c_pzrdat_out_bbh='bbhxxx_pzrdat'
c
 c_pzrdat2_out_bbh='bbhxxx_pzrdat2'
c
c
c  NOTE: For gam gam --> hadrons you must make an initial run with f_read_bbh=F
c        in order to produce the parton distribution files bbhxxx_pzrn, bbhxxx_pzrdat, 
c        and bbhxxx_pzrdat2.   Then set f_read_bbh=T and use the generated files as input.  The 
c        files bbhxxx_pzrn, bbhxxx_pzrdat, and bbhxxx_pzrdat2 only depend on pt_min_bbh and
c        ecm_bbh.    Thus, you only need to produce a new set of files if you change
c        pt_min_bbh and/or ecm_bbh.
c
c        For lepton pair production (i_lepton_bbh=11,13,15) use f_read_bbh=.false.  .
c        Lepton pair-production does not include mass effects, so you should run with 
c        pt_min_bbh and cm_isa_bbh well above pair-production threshold.
c
cf_read_bbh=.true.       ! read quark/gluon distrib.
 f_read_bbh=.false.      ! generate quark/gluon distrib.
c
cpt_min_bbh=1.6d0        ! minimum mini-jet or lepton pair pt in GeV
cpt_min_bbh=3.2d0        ! minimum mini-jet or lepton pair pt in GeV
 pt_min_bbh=0.2d0        ! minimum mini-jet or lepton pair pt in GeV
c
c
c
c
c  NOTE: For now there is no means to turn off minimum bias events. 
c        Thus, if you set f_dir_bbh=.true, f_1res_bbh=.false., 
c        and f_2res_bbh=.false., you will still get mostly 
c        minimum bias events.
c
c
 f_dir_bbh=.true.        ! include direct processes
cf_dir_bbh=.false.       ! don't include direct processes
c
 f_1res_bbh=.true.       ! include 1-resolved processes
cf_1res_bbh=.false.      ! don't include 1-resolved processes
c
 f_2res_bbh=.true.       ! include 2-resolved processes
cf_2res_bbh=.false.      ! don't include 2-resolved processes
c
c
c
c
 f_frag_bbh=.true.       ! fragment quarks and gluons
cf_frag_bbh=.false.      ! don't fragment quarks and gluons
c
cf_dltf_bbh=.true.       ! gam-gam ecm spec is delta func at max ecm
 f_dltf_bbh=.false.      ! gam-gam ecm spec is NOT delta func at max ecm
c
cf_star_bbh=.true.       ! use h2res for xsec when 1 jet pair pt>pt_min
 f_star_bbh=.false.      ! don't use h2res for xsec
c
 fbeam_bbh=.true.        ! include beamsstr. photons in photon lumi
cfbeam_bbh=.false.       ! don't include beamsstr. photons in ph lumi
c
 fbrem_bbh=.true.        ! include bremsstr. photons in photon lumi
cfbrem_bbh=.false.       ! don't include bremsstr. photons in ph lumi
c
 f_full_bbh=.true.       ! use full expression for w-w spectrum 
cf_full_bbh=.false.      ! only the leading term is kept
c
 f_pskn_bbh=.true.       !  q2 is taken from eqn. 3.12 of slac-pub-5873
cf_pskn_bbh=.false.      !  q2=q2_cut
c
c
ccm_isa_bbh=2.0d0        ! min gam-gam ecm in GeV for event gen.
ccm_isa_bbh=5.0d0        ! min gam-gam ecm in GeV for event gen.
 cm_isa_bbh=0.4d0        ! min gam-gam ecm in GeV for event gen.
 xk_bea_bbh=0.001d0      ! minimum beamstr photon energy in GeV
 cm_dlt_bbh=1000.d0      ! gam-gam ecm in GeV for f_dltf_bbh=.true. 
c
 i_cuts_bbh=0            ! =0 no y cuts,   =1 implement y cuts
 yy_cut_bbh=1.5d0        ! y cut for dir & 1-res if i_cut=1
c                          (no y cut for 2-res processes)
c
c
c
c
c JLC/NLC0500 8/8/98
c   ecm_bbh=515.d0
c   xn_bbh=0.95d10
c   sx_bbh=330.d-9
c   sy_bbh=4.9d-9
c   sz_bbh=120.d-6
c   betx_bbh=12.d-3
c   bety_bbh=0.12d-3
c
c JLC/NLC1000 8/8/98
    ecm_bbh=1008.d0
    xn_bbh=0.95d10
    sx_bbh=235.d-9
    sy_bbh=3.9d-9
    sz_bbh=120.d-6
    betx_bbh=12.d-3
    bety_bbh=0.15d-3
c
c
c
c
c NLC1500 (Modified RRuth Snowmass 96)
c   ecm_bbh=1500.d0
c   xn_bbh=1.10d10
c   sx_bbh=202.d-9
c   sy_bbh=5.10d-9
c   sz_bbh=150.d-6
c   betx_bbh=12.d-3
c   bety_bbh=0.20d-3
c
 NEVENT_STD=100
 Call go
 Call HEP2NT(' SUM')
 Call bbh_exit 
 Stop
 End
quit
echo "end of runbbh"