* Steering file for * WPHACT v. 1.9 * authors: Elena Accomando, Alessandro Ballestrero, Ezio Maina. ****************************************************************************** * * This program, its steering file, explanations and examples, as * well as other versions of WPHACT can be found with a link to * http://www.to.infn.it/~ballestr/wphact * and its subdirectories * * For questions, comments, remarks and requests, send a mail to * ballestrero@to.infn.it * ****************************************************************************** C This is the first official version in which all fully massive amplitudes C are computed by WPHACT. The user can choose among massless computations, C fully massive ones and processes in which only b's (or taus, or c's) masses C are exactly accounted for. ****************************************************************************** C This example of a steering file contains several comments and explanations C The expert user may use a shorter one in which most comments are dropped. ****************************************************************************** C From this version a new way of giving the input is used. C C From the point of view of the user the sintax is almost identical to C Cern library routine FFREAD. Routines internal to WPHACT are however used, C so one does not need to link the Cern library and, most important, C the real variables can (and must) be given in double precision. C C For the time being, if one changes the value of ireadinput in the data from C 0 to 1, one can still use the old way of answering the read statements. C Syntax: C Lines must not contain more than 80 columns. C * and C characters at the beginning of a line identify it as a comment line C The variable to be read must be specified as the first word of a line C (needs not to be in colun 1). Its value (values) must follow it. C Continuation on following lines of the values is allowed. C Any variable which is not needed for the process under study may remain C without any problem. The program will just ignore it and its value(s). C The advice is to change what is needed for the run at hand, and to leave all C other values used in preceding runs as they are. C The order of the variables to define is irrelevant. The user may change it. C For instance iproc can be defined first and e_cm second, contrary to the C the order below. C In output, however, the order will always be the one of this example. ****************************************************************************** C Input variables and explanations C -------------------------------- C The various final states are described in table proc.ps or can be deduced C from subroutine INITIALIZE. The final particles are enumerated from 3 C to 6 according to the order in which they appear in proc.ps and their C 4-momenta are p3(0:3), p4(0:3), p5(0:3), p6(0:3). C The 4momentum of e+ is p1(0:3), that of e- is p2(0:3) C The first variable to specify is the centre of mass energy: C e_cm ! centre of mass energy (GeV) e_cm 200.d0 C iproc !selects the kind of process. See the table proc.ps or the C subroutine INITIALIZE to find the number of the desired C process. Different final states which have the same diagrams C to compute are grouped in the same iproc, but have a different C ich. iproc 4 C ich ! selects the channel. As explained for iproc ich 1 C imass !imass=0 massless, imass=1 massive. This regards the first C 32 processes. If one choses the remaining ones, which C correspond to Higgs signal and irreducible backgrounds, C the masses of the particles to which the Higgs can decay C are accounted, all others are considered massless. imass 1 ***** The following inputs (until *%%) allow to choose part of the diagrams * The first one (icc3) refers to all processes which contain CC3 diagrams * The second one (iccnc) must always be specified for Mixed Processes. * The remaining ones are instead effective only for fully massive computations * (imass=1) C icc3 ! yes/no CC3 only . Obviously it regards only processes C which contain CC3 diagrams. C The convention for WPHACT is 1=yes, 0=no icc3 0 C IF (iproc.GE.6.AND.iproc.LE.8) !in mixed processes: C iccnc ! 1=CC only, 2=NC only, 3= CC+NC +interf C this flag decides whether all Mixed process is computed C or only its charged current (or neutral current) part. C The choice has to be done both for massive and massless computations iccnc 3 C izz !selects: 1=NC02, 0=all, -1=all-NC02 C NC02 here refers to all diagrams of NC processes with C two Z's decaying to final particles. They are 2 for C no identical particles in the final states, 4 otherwise izz 0 C izg !selects: 1=zeta-gamma, 0=all, -1=all-(zeta-gamma) C with zeta-gamma we refer to all diagrams with one Z and C one gamma* decaying to final particles izg 0 C izg34 !selects: 1=zeta56-gamma34, 0=all,-1=all-(Z56-gam34) C with zeta56-gamma34 we refer to diagrams in which one Z C decays to particles 5 and 6 and one gamma* decaying to C particles 3 and 4 izg34 0 C izg56 !selects: 1=zeta34-gamma56, 0=all,-1=all-(Z34-gam56) C same as before with zeta -> 34 and gamma* -> 56 izg56 0 C inc08 !selects: 1=NC08, 0=all, -1=all-NC08 C NC08 indicates all NC diagrams with two Z's, a Z and C a gamma*, or two gamma* decaying to final particles inc08 0 C igamgam !selects: 1=gamma-gamma, 0=all, -1=all-(gamma-gamma) C with gamma-gamma we indicate NC diagrams for a final C state containing at least a couple e+e-, in which there C are two t-channel photon propagators. igamgam 0 C itch !selects: 1=t-channel, 0=all, -1= all - tchannel C with t-channel we indicate all diagrams where there is at C least one t-channel boson propagator itch 0 *%% C The following flag ismallangle has to be set different from 0 when one wants C to switch from the default phase space mapping which accounts for all C possible s-channel peaking structures, to the one which takes care of t C channel peaks which become important for outgoing electrons near the beam C ismallangle ! 0=no small angle, 1=e- near the beam (180 degrees) ! 2= e- near 180 degrees, e+ near 0 degrees ismallangle 0 C here one chooses wheter to compute Initial State Radiation C The value 0 corresponds to no ISR, 1 to ISR via structure functions, C 2 to ISR generated by QEDPS program (by Munehisa, Kurihara, Fujimoto and C Shimizu) C isr ! yes/no ISR (2=qedps) isr 0 C ibeam !yes/no beamstrahlung ibeam 0 C icoul ! yes/no Coulomb corrections icoul 1 C istrcor ! yes/no 'naive' QCD corrections (also for widths) istrcor 1 C ianc ! yes/no anomalous couplings ianc 0 C IF(ianc.EQ.1) C delz,xf,xz,yf,yz,zz !anomalous couplings parameters delz 0.d0 xf 0.d0 xz 0.d0 yf 0.d0 yz 0.d0 zz 0.d0 C for single W processes C if (iproc.eq.2.or.iproc.eq.4.or.(iproc.eq.7.and.iccnc.ne.2)) C one can choose or not to compute imaginary fermion loop contributions: C ifloop yes/no Im. fermion loop contributions C if ifloop=1, then the following ipr has to be set =2 C for single W processes gauge invariance requires that either ifloop=1 C and ipr=2, or fixed widths are chosen: ipr=0 ifloop 0 C ipr ! widths: 0/1 Z,W fixed/running; 2 Z fixed, W running ipr 1 C iswgcomp ! yes/no s2w and g computed (0= DATA value) C iswgcomp=1 corresponds to G_f scheme. Other schemes can be implemented C by fixing the EW parameters in the data statement. iswgcomp 1 C igwcomp,igzcomp,ighcomp C ! yes/no W,Z,H width computed (0= DATA val.) C If 1 is chosen, the widths are computed by WPHACT in a way consistent C with ME evaluation. 0 corresponds to taking the values from the DATA igwcomp 1 igzcomp 1 ighcomp 1 ***** Cuts are implemented as in WPHACT 1.0. If needed, further details can be * found in the CPC 99(1997)270 where it is described C First choose if you want or not to use the cuts defined by WPHACT: C icut ! yes/no cuts icut 1 ** IF(icut.EQ.1) give all cuts until *%% C e_min ! 4 energy lower cuts (GeV) e_min 0.d0 0.d0 0.d0 0.d0 C e_max ! 4 energy upper cuts (GeV) e_max 5000.d0 5000.d0 5000.d0 5000.d0 C rm_min ! 6 invariant mass lower limits (GeV) C ! (34, 35, 36, 45, 46, 56) rm_min 10.d0 10.d0 10.d0 5.d0 5.d0 5.d0 C rm_max ! 6 invariant mass upper limits (GeV) rm_max 5000.d0 5000.d0 5000.d0 5000.d0 5000.d0 5000.d0 C pt_min ! 4 transverse momenta lower cuts (GeV) pt_min 0.d0 0.d0 0.d0 0.d0 C pt_max ! 4 transverse momenta upper cuts (GeV) pt_max 5000.d0 5000.d0 5000.d0 5000.d0 C icos ! angular cuts may be expressed in degrees (0) or C with the value of the cosine (1) icos 0 C thbeam_min! 4 particle-beam angle lower (in degrees) cuts thbeam_min 5.d0 5.d0 5.d0 5.d0 C thbeam_max! 4 particle-beam angle upper (in degrees) cuts thbeam_max 175.d0 175.d0 175.d0 175.d0 C thsep_min ! 6 particl-particl angle lower(in degrees) cuts thsep_min 0.d0 0.d0 0.d0 0.d0 0.d0 0.d0 C thsep_max ! 6 particl-particl angle upper(in degrees) cuts thsep_max 180.d0 180.d0 180.d0 180.d0 180.d0 180.d0 *%% ** additional cuts may be implemented by the user in WPHACT just after ** the line " * here define additional cuts ", where a commented example ** is reported ***** Distributions are implemented as in WPHACT 1.0. If needed, * further details can be found in the CPC 99(1997)270 * where it is described * First choose if you want or not distributions: C idistr ! yes/no distributions idistr 1 * If distributions have to be automatically computed, they must first be * defined by the user after the lines: " * here define eventual weighted distributions (or use the include file) * include 'abdis.dis' " * where some commented examples of distributions are present. C IF(idistr.EQ.1) , give the number of distributions C ndistr ! number of distributions (limited by WPHACT parameter C ndismax, whose present value is 50) ndistr 2 C DO i=1,ndistr C nsubint(i) !number of sub-intervals with C ! different binning (limited by WPHACT parameter C nintmax, whose present value is 10 ) nsubint(1) 1 C (distr_estrinf(i,j),j=1,nsubint(i)+1) !lower limits of each C ! subint (which coincide with the upper limit of the C ! previous one) + upper limit of the last subint distr_estrinf(1) 0.d0 100.d0 C (nbin_number(i,j),j=1,nsubint(i)) !number of bins in each C subint (limited by WPHACT parameter nbinmax, whose C present value is 500) nbin_number(1) 100 * add other distributions till ndistr total number is reached nsubint(2) 2 distr_estrinf(2) 1.d0 100.d0 150.d0 nbin_number(2) 100 100 nsubint(3) distr_estrinf(3) nbin_number(3) * .... C ENDDO !i C ENDIF ***** Unweighted event generation is implemented as in WPHACT 1.0. If needed, * further details can be found in the CPC 99(1997)270 * where it is described C iflat ! yes/no flat (i.e. unweighted) event generation iflat 0 C IF(iflat.EQ.1)THEN C scalemax ! scale factor for the maximum scalemax 1.1 C istorvegas ! yes/no VEGAS data stored istorvegas 1 C irepeat ! 0=normal;1=repeat 2nd only;2=rep. C ! with nflevts fixed C ! 1=repeat only second iteration using vegas C ! data stored in previous run, C ! 2=repeat, generating an exact number nflevts C ! of flatevents using vegas data stored C ! in previous run irepeat 0 C IF(irepeat.eq.2)THEN C nflevts ! number of events to be generated as C ! described above nflevts C END IF C istormom ! yes/no momenta of flat events written in C ! .dat files istormom 0 C ijetset ! yes/no call to Jetset ijetset 0 C IF(iproc.GE.6.AND.iproc.LE.8)THEN C interf ! 0= mix. interference added to NC; C ! 1= added to CC interf 1 C ENDIF C ENDIF ***** Integration by VEGAS is implemented as in WPHACT 1.0. If needed, * further details can be found in the CPC 99(1997)270 where it * is described C acc ! integration accuracy acc .0001 C iterm ! yes/no thermalization iterm 1 C ncall_term ! thermalization calls per iteration ncall_term 100000 C itmx_term ! thermalization iterations C (limited by WPHACT parameter nitmax, whose present value is 10) itmx_term 3 C ncall ! integration calls per iteration ncall 1000000 C itmx ! integration iterations (2 for flat event generation) C (limited by WPHACT parameter nitmax, whose present value is 10) itmx 3 ********** The following inputs have to be given only for ** Higgs processes and their irreducible backgrounds which correspond ** to processes with iproc greater than 32, for which only the masses ** of the particles to which the higgs can decay are exactly accounted ** for (b's or tau's or c's). C rmb ! quark b mass (GeV) rmb C rmc ! quark c mass (GeV) rmc C rmtau ! tau mass (GeV) rmtau C rmb_run ! quark b mass running(GeV) rmb_run C rmc_run ! quark c mass running(GeV) rmc_run C gamh ! higgs width gamh C alfas_h ! alfas at Higgs mass alfas_h C icch ! 1=Higgs signal, 2=Background, 3=Higgs+Backg+interf icch C IF(icch.NE.2) one has to specify whether one wants SM or MSSM higgs C isusy=0 for SM, isusy=1 for MSSM isusy C IF(isusy.EQ.0) one has to specify the following two inputs: C iha only if (iproc.EQ.39.OR.iproc.EQ.43.OR.iproc.EQ.51).AND.icch.EQ.1) C iha=1 computes all the diagrams with Higgs, C iha=2 computes only the diagram hZ iha 1 C rmh is the higgs mass in GeV C rmh ! Higgs mass (GeV) rmh 100.d0 C IF(isusy.EQ.1) one has to specify all the remaining inputs: C IF((iproc.EQ.39.OR.iproc.EQ.43.OR.iproc.EQ.51).AND. C & icch.EQ.1)THEN C iha ! 1=all, 2=h, 3=A, 4=h-stralhung, 5=A-strahlung iha C ENDIF C irmhcomp ! yes/no rmh computed irmhcomp C IF(irmhcomp.EQ.1)THEN C rma ! Higgs 'A' mass (GeV) rma C tgb ! tan(beta) tgb C iloop ! 1=1-loop, 2=Carena et al. corrections iloop C IF(iloop.EQ.2)THEN C imixing !1=no, 2=maximal, 3=typical, 4=other imixing C IF(imixing.EQ.4)THEN C At At C Ab Ab C rmyou rmyou C ENDIF C ENDIF C ELSE C rmh rmh C rma rma C tgb tgb C ENDIF C ENDIF C ENDIF C ENDIF