#!/bin/tcsh
cd .
pythia<<quit 
!show routines
!file go
 include 'pythcom.f'
 include 'bmscom.f'
 include 'pysubs.f'
 include 'pypars.f'
 include 'ludat2.f'
 include 'ludat3.f'
 include 'ludatr.f'
 include 'wwhpcom.f'
c
 i_beam_target_bms=1  ! 1='e- e+'  2='e- gam'  3='gam e+'  4='gam gam'
c
cmrlu(1)=19xxx503 ! ran # seed (0.le.mrlu(1).le.900000000, def=19780503)
 mrlu(1)=19378503 ! ran # seed (0.le.mrlu(1).le.900000000, def=19780503)
c
 mrlu(2)=0
c
 pmas(6,1)=175.
 pmas(23,1)=91.175
 pmas(24,1)=80.22
c
 msel=0         ! process selection:  0=determined by msub 
 msub(1)=1      ! (gamma*/Z)
 msub(22)=0     ! (gamma*/Z) (gamma*/Z)
 msub(25)=0     ! W+ W-
 msub(35)=0     ! e (gamma*/Z)
 msub(36)=0     ! nu W+
 msub(191)=0    ! grace process
 msub(192)=0    ! W+W- ala WWHPZH
c
 mstp(14)=10    ! beam or target photon structure variable
c
 mstp(11)=0     ! pythia initial-state bremsstrahlung:  0=off  1=on
 mstp(42)=1     ! resonance mass:  0=force on mass shell  1=Breit-Wigner
 mstp(43)=3     ! 1= gamma only  2= Z only  3= gamma/Z interference
 mstp(71)=1     ! final state rad master switch:  0=off  1=on
 mstp(111)=1    ! fragmentation & decay:  0=off  1=on
 mstp(123)=1    ! reaction to max wgt viol: 0=stop 1=only_print 2=adjust_wgt
c
 m_beams_bms=1          ! pybms beamstrahlung simulation:  0=off  1=on
 m_brems_bms=1          ! pybms bremsstrahlung simulation:  0=off  1=on
 m_linac_e_bms=0        ! linac energy spread simulation:  0=off 1=on
 wgt_inter_bms=2.d0     ! rejection wgt internal to beams. code
c
c mdme(144,1)=0   ! gamma --> dq dqb
c mdme(145,1)=0   ! gamma --> uq uqb
c mdme(146,1)=0   ! gamma --> sq sqb
c mdme(147,1)=0   ! gamma --> cq cqb
c mdme(148,1)=0   ! gamma --> bq bqb
c mdme(149,1)=1   ! gamma --> tq tqb
c mdme(152,1)=0   ! gamma --> el elb
c mdme(153,1)=0   ! gamma --> mu mub
c mdme(154,1)=0   ! gamma --> ta tab
c
c mdme(156,1)=0   ! Z0 --> dq dqb
c mdme(157,1)=0   ! Z0 --> uq uqb
c mdme(158,1)=0   ! Z0 --> sq sqb
c mdme(159,1)=0   ! Z0 --> cq cqb
c mdme(160,1)=0   ! Z0 --> bq bqb
c mdme(161,1)=1   ! Z0 --> tq tqb
c mdme(164,1)=0   ! Z0 --> el elb
c mdme(165,1)=0   ! Z0 --> ne neb
c mdme(166,1)=0   ! Z0 --> mu mub
c mdme(167,1)=0   ! Z0 --> nm nmb
c mdme(168,1)=0   ! Z0 --> ta tab
c mdme(169,1)=0   ! Z0 --> nt ntb
c
c mdme(172,1)=1   ! W+ --> dqb uq
c mdme(173,1)=1   ! W+ --> dqb cq
c mdme(174,1)=0   ! W+ --> dqb tq
c mdme(176,1)=1   ! W+ --> sqb uq
c mdme(177,1)=1   ! W+ --> sqb cq
c mdme(178,1)=0   ! W+ --> sqb tq
c mdme(180,1)=0   ! W+ --> bqb uq
c mdme(181,1)=0   ! W+ --> bqb cq
c mdme(182,1)=0   ! W+ --> bqb tq
c mdme(188,1)=0   ! W+ --> elb ne
c mdme(189,1)=0   ! W+ --> mub nm
c mdme(190,1)=0   ! W+ --> tab nt
c
c
ci_stdxwrt=0   !  PAW ntuple output format (HEPEVT common contents)
 i_stdxwrt=1   !  stdHEP output format (HEPEVT common contents)
c
c  datfile_bms is the stdHEP (i_stdxwrt=1) or PAW ntuple (i_stdxwrt=0) output file
c
 datfile_bms='pybxxx.dat'
c
c
 chfile_bms='pybxxx_hbook'
 grcfile_bms='grc538.dat'
 i_hbook_bms=0
 i_just_plot_bms=0
 i_allms_plot_bms=0
 i_beams_plot_bms=0
 i_brems_plot_bms=0
 xlow_hist_bms=0.0d0
 ncall_xsec_bms=500
 itmx_xsec_bms=3
 ncall_beams_bms=100000
 itmx_beams_bms=5
 nprn_bms=0
 i_just_calculate_bms=0
c
   nexternal=0
c
   nexternal=nexternal+1
   isub_ext(nexternal)=191
   title_ext(nexternal)=' e- + e+ --> X (grace)'
c
   nexternal=nexternal+1
   isub_ext(nexternal)=192
   title_ext(nexternal)=' e- + e+ --> W+W- (wwhpzh)'
c
c
ccccccccc    BEGIN parameters for isub=192  e+e- --> W+W-  wwhpzh version
c
c  pe_wwhp=1.d0       ! initial state e- polarization
c  pe_wwhp=-1.d0      ! initial state e- polarization
   pe_wwhp=0.d0       ! initial state e- polarization
   ctmn_wwhp=-1.d0   ! minimum cos(W- production angle)
   ctmx_wwhp=1.d0    ! maximum cos(W- production angle)
c
ccccccccc    END   parameters for isub=192  e+e- --> W+W-  wwhpzh version
c
   roots_p(4)=1050.
   roots_p(5)=1200.
   roots_p(6)=1350.
   roots_p(7)=1425.
   roots_p(8)=1500.
c
   cross_p(1)=0.
c
c N_100  0% polarization
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.2528
c   cross_p(3)=0.5748
c   cross_p(4)=0.8685
c   cross_p(5)=1.1231
c   cross_p(6)=1.5336
c   cross_p(7)=1.5547
c   cross_p(8)=1.7596
c
c N_100  100% e-(L)
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.3984
c   cross_p(3)=1.0153
c   cross_p(4)=1.6227
c   cross_p(5)=2.0861
c   cross_p(6)=2.7570
c   cross_p(7)=3.0805
c   cross_p(8)=3.3887
c
c N_1000  0% polarization
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.2626
c   cross_p(3)=0.7287
c   cross_p(4)=1.2998
c   cross_p(5)=1.9502
c   cross_p(6)=2.7990
c   cross_p(7)=3.2016
c   cross_p(8)=3.7832
c
c N_1000  100% e-(L)
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.4141
c   cross_p(3)=1.2626
c   cross_p(4)=2.4796
c   cross_p(5)=3.7514
c   cross_p(6)=5.4967
c   cross_p(7)=6.4013
c   cross_p(8)=7.5084
c
c N_let  0% polarization
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.2557
c   cross_p(3)=0.6481
c   cross_p(4)=1.0883
c   cross_p(5)=1.5308
c   cross_p(6)=2.0495
c   cross_p(7)=2.3227
c   cross_p(8)=2.6654
c
c N_let  100% e-(L)
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.4103
c   cross_p(3)=1.1799
c   cross_p(4)=2.0666
c   cross_p(5)=2.9778
c   cross_p(6)=4.0422
c   cross_p(7)=4.6665
c   cross_p(8)=5.2724
c
c E_100
c   roots_p(1)=2.*pmas(6,1)
c   roots_p(2)=600.
c   roots_p(3)=850.
c   cross_p(2)=0.3118
c   cross_p(3)=1.3620
c   cross_p(4)=2.6161
c   cross_p(5)=3.7245
c   cross_p(6)=4.9386
c   cross_p(7)=5.5749
c   cross_p(8)=6.2311
c
c W_100
c   roots_p(1)=2.*pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=3.6865
c   cross_p(3)=11.5321
c   cross_p(4)=23.5431
c   cross_p(5)=31.2959
c   cross_p(6)=39.0539
c   cross_p(7)=43.1393
c   cross_p(8)=47.0980
c
c EW_100  0<pt34<3000
c   roots_p(1)=2.*pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=136.299
c   cross_p(3)=469.633
c   cross_p(4)=899.819
c   cross_p(5)=1124.895
c   cross_p(6)=1360.391
c   cross_p(7)=1432.977
c   cross_p(8)=1577.932
c
c EW_100  0<pt34<30
c   roots_p(1)=2.*pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=135.425
c   cross_p(3)=451.259
c   cross_p(4)=834.096
c   cross_p(5)=1019.134
c   cross_p(6)=1201.060
c   cross_p(7)=1280.035
c   cross_p(8)=1377.132
c
c EW_100  30<pt34<30000
c   roots_p(1)=2.*pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=0.
c   cross_p(3)=14.827
c   cross_p(4)=69.364
c   cross_p(5)=112.728
c   cross_p(6)=157.312
c   cross_p(7)=182.023
c   cross_p(8)=208.978
c
c WZ_100
c   roots_p(1)=pmas(23,1)+pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=2.950
c   cross_p(3)=13.091
c   cross_p(4)=32.090
c   cross_p(5)=42.439
c   cross_p(6)=53.784
c   cross_p(7)=59.183
c   cross_p(8)=65.239
c
c WBZ_100
c   roots_p(1)=pmas(23,1)+pmas(24,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=2.104
c   cross_p(3)=12.510
c   cross_p(4)=31.976
c   cross_p(5)=42.084
c   cross_p(6)=54.297
c   cross_p(7)=60.808
c   cross_p(8)=67.711
c
c WWZ_100
c   roots_p(1)=pmas(23,1)+2.*pmas(24,1)
c   roots_p(2)=550.
c   roots_p(3)=800.
c   cross_p(2)=48.511
c   cross_p(3)=61.769
c   cross_p(4)=62.956
c   cross_p(5)=60.212
c   cross_p(6)=57.525
c   cross_p(7)=56.353
c   cross_p(8)=55.703
c
c ZZZ_100
c   roots_p(1)=3.*pmas(23,1)
c   roots_p(2)=550.
c   roots_p(3)=800.
c   cross_p(2)=1.1640
c   cross_p(3)=1.0858
c   cross_p(4)=0.9437
c   cross_p(5)=0.8602
c   cross_p(6)=0.7668
c   cross_p(7)=0.7613
c   cross_p(8)=0.7248
c
c Z_100
c   roots_p(1)=2.*pmas(23,1)
c   roots_p(2)=460.
c   roots_p(3)=750.
c   cross_p(2)=0.593
c   cross_p(3)=3.112
c   cross_p(4)=8.221
c   cross_p(5)=11.471
c   cross_p(6)=14.984
c   cross_p(7)=16.891
c   cross_p(8)=18.732
c
c TT  100% e-(L)
    roots_p(1)=2.*pmas(6,1)
    roots_p(2)=600.
    roots_p(3)=850.
    cross_p(2)=609.823
    cross_p(3)=326.842
    cross_p(4)=218.310
    cross_p(5)=168.429
    cross_p(6)=133.745
    cross_p(7)=120.270
    cross_p(8)=108.704
c
c TTZ_100
c   roots_p(1)=2.*pmas(6,1)+pmas(23,1)
c   roots_p(2)=650.
c   roots_p(3)=850.
c   cross_p(2)=4.061
c   cross_p(3)=5.041
c   cross_p(4)=4.845
c   cross_p(5)=4.516
c   cross_p(6)=4.123
c   cross_p(7)=3.952
c   cross_p(8)=3.772
c
cxmhat_min_bms=2.            ! min. e+e- ecm following brems/beams (GeV)
cxmhat_min_bms=roots_p(1)    ! min. e+e- ecm following brems/beams (GeV)
cxmhat_min_bms=2.*pmas(6,1)  ! min. e+e- ecm following brems/beams (GeV)
 xmhat_min_bms=2.*pmas(24,1)+0.5
c
 i_pyup_sigh=3                ! cross-section calculation by pyup_sigh:
c                               1 for constant cross-section; 2 for 1/s;
c                               3 if specified by roots_p & cross_p
c
    freq_bms=120.d0
c
c   ecm=192.
c
c
c
c
c JLC/NLC0500 8/8/98
   ecm=500.
   xn_bms=0.95d10
   sx_bms=330.d-9
   sy_bms=4.9d-9
   sz_bms=120.d-6
   betx_bms=12.d-3
   bety_bms=0.12d-3
c
c
c
c JLC/NLC1000 8/8/98
c   ecm=1008.
c   xn_bms=0.95d10
c   sx_bms=235.d-9
c   sy_bms=3.9d-9
c   sz_bms=120.d-6
c   betx_bms=12.d-3
c   bety_bms=0.15d-3
c
c
c
c NLC1500 (Modified RRuth Snowmass 96)
c    ecm=1500.
c    xnb_bms=100.d0
c    xn_bms=1.10d10
c    sx_bms=202.d-9
c    sy_bms=5.10d-9
c    sz_bms=150.d-6
c    betx_bms=12.d-3
c    bety_bms=0.20d-3
c
cNEVENT_STD=1000
 NEVENT_STD=10
 Call INIPYTH
 Call go
 Call HEP2NT(' SUM')
 Call pybms_exit 
 Stop
 End
quit
echo "end of runpyb"