gtar -xvzf koralw-1.03-export.tar.gz
you do the following
cd koralw-1.03-export/demo
make demo
The above will compile and run demonstarion program on any unix machine
and compare the pre-prepared output (hp) with the actual output.
There are also three more complicated demos (make demo1, make demo2, make demo2a) with external data files, some more tests and histograms. In particular demo2a runs with anomalous couplings set to wild anomalous values.
koralw-1.03-export/korww low-level MC generator, technics! koralw-1.03-export/model matrix element, physics! koralw-1.03-export/glib histograming package and utilities koralw-1.03-export/kwlund JETSET and PHOTOS koralw-1.03-export/demo demonstration program koralw-1.03-export/tauola TAUOLA koralw-1.03-export/interfaces TAUOLA, JETSET and PHOTOS interfaces koralw-1.03-export/semian semianalytical routine KORWAN koralw-1.03-export/ampli4f dummy external matrix element
KeyRad = 1000*KeyCul+100*KeyNLL+10*KeyFSR+KeyISR
KeyCul=0 Coulomb corr. OFF
KeyCul=0 Coulomb corr. ON
KeyNLL=0 sets Next-to Leading alpha/pi terms to zero
KeyNLL=1 alpha/pi in YFS formfactor is kept
KeyFSR Fnitial State Radiation switch, INACTIVE
KeyISR=0 Initial State Radiation OFF,
KeyISR=1 Initial State Radiation ON.
KeyPhy = 100000*KeyWu +10000*KeyRed +1000*KeySpn
+100*KeyZet +10*KeyMas +KeyBra
KeyBra= 0 Born branching ratios, no mixing
KeyBra= 1 Branching ratios with mixing and QCD
KeyMas= 0 Massles kinematics for W decay products
KeyMas= 1 Massive kinematics for W decay products
KeyZet= 0 Z width in Z propagator: s/M_Z *GAMM_Z
KeyZet= 1 Z width in Z propagator: M_Z *GAMM_Z
KeyZet= 2 no (0) Z width in Z propagator.
KeySpn= 0-off, 1-on for spin effects in W decays
KeyRed= Reduction of massive FS to massles Matr.El.
KeyRed= 0 fine
KeyRed= 1 crude, 4-mom. non conserving
KeyWu= 0 W width in W propagator: s/M_W *GAMM_W
KeyWu= 1 W width in W propagator: M_W *GAMM_W
KeyWu= 2 no (0) W width in W propagator.
KeyTek= 10*KeyRnd +KeyWgt
KeyWgt =0, WTMOD=1 useful for apparatus Monte Carlo.
KeyWgt =1, WTMOD varying, option faster and safer
KeyWgt =2, WTMOD=1 for internal matrix el, AND varying for external
matrix el.
KeyRnd = not implemented ???????????
KeyMis = 100*KeyACC +10*Key4f +KeyMix
miscelaneus auxiliarym key for tests
actualy KeyMix=0 causes sinw2 to be according LEP2 presription
while KeyMix=1 prowides sinw2 as in 1.00 version.
Key4f = 0, external matrix el. OFF
Key4f = 1, external matrix el. ON
KeyACC = 0, anomalous WWV couplings in internal matr. el. OFF
KeyACC = 1, anomalous WWV couplings in internal matr. el. ON
KeyDwm decay channel of W-
KeyDwp decay channel of W+
0-all chann. according to br. ratios
1-ud, 2-cd, 3-us, 4-cs, 5-ub, 6-cb, 7-e, 8-mu, 9-tau
negative number is reset to some default value
nout output unit number of the generator
jak1 TAUOLA, tau of W+ decay channel
jak2 TAUOLA, tau of W- decay channel
itdkrc TAUOLA, radiative corr. in leponic tau decays switch
ifphot PHOTOS, activation switch
ifhadm JETSET, W- hadronisation activation switch
ifhadp JETSET, W+ hadronisation activation switch
!!!! IFHADM and IFHADP are glued together (for now) !!!!
!-------
CMSENE : CMS energy [GeV]
GFERMI : Fermi constant
ALFWIN : alpha QED at WW treshold scale (inverse)
AMAZ : Z mass [GeV]
GAMMZ : Z width [GeV]
AMAW : W mass [GeV]
GAMMW : W width [GeV], For GAMMW < 0 RE-CALCULATED inside program
VVMIN : Photon spectrum parameter (dummy IR cut)
VVMAX : Photon spectrum parameter (max generated v, v = 1-sprim/s)
WTMAX : max wt used for rejection in non weighted mode,
!-------
g1(2) : g_1^Z, anomalous coupling for WWZ vertex
kap(2) : kappa_Z, anomalous coupling for WWZ vertex
lam(2) : lambda_Z, anomalous coupling for WWZ vertex
g4(2) : g_4^Z, anomalous coupling for WWZ vertex
g5(2) : g_5^Z, anomalous coupling for WWZ vertex
kapt(2): kappa-tilde_Z, anomalous coupling for WWZ vertex
lamt(2): lambda-tilde_Z, anomalous coupling for WWZ vertex
!-------
g1(1) : g_1^g, anomalous coupling for WWg vertex
kap(1) : kappa_g, anomalous coupling for WWg vertex
lam(1) : lambda_g, anomalous coupling for WWg vertex
g4(1) : g_4^g, anomalous coupling for WWg vertex
g5(1) : g_5^g, anomalous coupling for WWg vertex
kapt(1): kappa-tilde_g, anomalous coupling for WWg vertex
lamt(1): lambda-tilde_g, anomalous coupling for WWg vertex
These inputs are sent to KORALW via arguments NPAR and XPAR:
NPAR(1)=KeyRad
NPAR(2)=KeyPhy
NPAR(3)=KeyTek
NPAR(4)=KeyMis
NPAR(5)=KEYDWM
NPAR(6)=KEYDWP
NPAR(7)=NOUT
NPAR(21)=JAK1
NPAR(22)=JAK2
NPAR(23)=ITDKRC
NPAR(24)=IFPHOT
NPAR(25)=IFHADM
NPAR(26)=IFHADP
XPAR(1)=CMSENE
XPAR(2)=GMU
XPAR(3)=ALFWIN
XPAR(4)=AMAZ
XPAR(5)=GAMMZ
XPAR(6)=AMAW
XPAR(7)=GAMMW
XPAR(8)=VVMIN
XPAR(9)=VVMAX
XPAR(10)=WTMAX
XPAR(21) - XPAR(57):
g1(1) = DCMPLX(xpar(21),xpar(31))
kap(1) = DCMPLX(xpar(22),xpar(32))
lam(1) = DCMPLX(xpar(23),xpar(33))
g4(1) = DCMPLX(xpar(24),xpar(34))
g5(1) = DCMPLX(xpar(25),xpar(35))
kapt(1) = DCMPLX(xpar(26),xpar(36))
lamt(1) = DCMPLX(xpar(27),xpar(37))
!----
g1(2) = DCMPLX(xpar(41),xpar(51))
kap(2) = DCMPLX(xpar(42),xpar(52))
lam(2) = DCMPLX(xpar(43),xpar(53))
g4(2) = DCMPLX(xpar(44),xpar(54))
g5(2) = DCMPLX(xpar(45),xpar(55))
kapt(2) = DCMPLX(xpar(46),xpar(56))
lamt(2) = DCMPLX(xpar(47),xpar(57))
Other inputs in subroutine FILINP:
AMEL : electron mass for QED part [GeV]
ALFINV : inverse coupling const for QED part
GPICOB : conversion to picobarns
BR(1) - BR(9) : branching ratios
Note! It is assumed that \sum BR(I) = 1
AMAFIN(1) - AMAFIN(16) : masses of final fermions [GeV]
production: e**4 = 16*pi**2*alphaw**2
decays: g**4 = e**4/sinw2**2
this is normalized to the value for (ev)(ev) channel
for i,j channels the above decay c.c. is replaced by:
g**4* br(i)/br(el) *br(j)/br(el)
and for total x-sect over various channels by:
g**4* 1/br(el) * 1/br(el) = g**4/br(el)**2
INPUT
svar: CMS energy squared [GeV].
vvmin: minimal v variable, in most cases should be set to 0.
vvmax: maximal v variable.
keymod: Defines type of Structure Functions used for IS Radiation
keymod= 0 No ISR, Born,
keymod=300 Zero Order, YFS style,
keymod=301 First Order, YFS style,
keymod=302 Second Order, YFS style,
keymod=303 Third Order, YFS style,
keymod=502 Second Order, Gribov-Kuraev-Fadin style,
keymod=310 First Order YFS Beta0 only,
keymod=311 First Order YFS Beta1 only,
keymod=320 Second Order YFS Beta0 only,
keymod=321 Second Order YFS Beta1 only,
keymod=322 Second Order YFS Beta2 only,
keymod<0 as (-keymod) but multiplied by v,
(d\sigma/d\log v), for tests only,
keymod>10000 Initialisation only.
keypre: Defines precision level of the computation.
For keymod=0 (No ISR, Born), in (e\nu_e) channel:
keypre=1 absolute error 1d-5 pb,
keypre=2 absolute error 1d-6 pb,
keypre=3 absolute error 1d-7 pb,
For keymod>0 (ISR), in (e\nu_e) channel:
keypre=1 absolute error 3d-5 pb,
keypre=2 absolute error 1d-5 pb,
keypre=3 absolute error 1d-6 pb,
keypre=4 absolute error 1d-7 pb.
OUTPUT
xsect: cross-section in picobarns.
errabs: absolute error in picobarns.