#!/bin/tcsh
cd .
ln -fs dat_filename dat_softlink
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
mrlu(1)=19xxx503 ! ran # seed (0.le.mrlu(1).le.900000000, def=19780503)
cmrlu(1)=19342503 ! 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
mdme(144,1)=0 ! gamma --> dq dqb
mdme(145,1)=0 ! gamma --> uq uqb
mdme(146,1)=0 ! gamma --> sq sqb
mdme(147,1)=0 ! gamma --> cq cqb
mdme(148,1)=0 ! gamma --> bq bqb
mdme(149,1)=1 ! gamma --> tq tqb
mdme(152,1)=0 ! gamma --> el elb
mdme(153,1)=0 ! gamma --> mu mub
mdme(154,1)=0 ! gamma --> ta tab
c
mdme(156,1)=0 ! Z0 --> dq dqb
mdme(157,1)=0 ! Z0 --> uq uqb
mdme(158,1)=0 ! Z0 --> sq sqb
mdme(159,1)=0 ! Z0 --> cq cqb
mdme(160,1)=0 ! Z0 --> bq bqb
mdme(161,1)=1 ! Z0 --> tq tqb
mdme(164,1)=0 ! Z0 --> el elb
mdme(165,1)=0 ! Z0 --> ne neb
mdme(166,1)=0 ! Z0 --> mu mub
mdme(167,1)=0 ! Z0 --> nm nmb
mdme(168,1)=0 ! Z0 --> ta tab
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=1
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
pe_wwhp=-1.d0 ! initial state e- polarization
c 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=15. ! min. e+e- ecm following brems/beams (GeV)
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)
xmhat_min_bms=2.*pmas(6,1)+0.5
cxmhat_min_bms=2.*pmas(24,1)+0.5
cxmhat_min_bms=2.* pmas(23,1)+0.5
cxmhat_min_bms=pmas(24,1)
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
xnb_bms=100.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=1000.
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 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=5540
Call INIPYTH
Call go
Call HEP2NT(' SUM')
Call pybms_exit
Stop
End
quit
echo "end of runpyb"