| 0795afa3 |
1 | \newpage |
| 2 | \section{Input\label{INPUT}} |
| 3 | |
| 4 | \subsection{Input Format} |
| 5 | |
| 6 | ISAJET is controlled by commands read from the specified input |
| 7 | file by subroutine READIN. (In the interactive version, this file is |
| 8 | first created by subroutine DIALOG.) Syntax errors will generate a |
| 9 | message and stop execution. Based on these commands, subroutine LOGIC |
| 10 | will setup limits for all variables and check for inconsistencies. |
| 11 | Several runs with different parameters can be combined into one job. |
| 12 | The required input format is: |
| 13 | \begin{verbatim} |
| 14 | Title |
| 15 | Ecm,Nevent,Nprint,Njump/ |
| 16 | Reaction |
| 17 | (Optional parameters) |
| 18 | END |
| 19 | (Optional additional runs) |
| 20 | STOP |
| 21 | \end{verbatim} |
| 22 | with all lines starting in column 1 and typed in {\it upper} case. These |
| 23 | lines are explained below. |
| 24 | |
| 25 | Title line: Up to 80 characters long. If the first four letters |
| 26 | are STOP, control is returned to main program. If the first four letters |
| 27 | are SAME, the parameters from previous run are used excepting those |
| 28 | which are explicitly changed. |
| 29 | |
| 30 | Ecm line: This line must always be given even if the title is |
| 31 | SAME. It must give the center of mass energy (Ecm) and the number of |
| 32 | events (Nevent) to be generated. One may also specify the number of |
| 33 | events to be printed (Nprint) and the increment (Njump) for printing. |
| 34 | The first event is always printed if Nprint $>$ 0. For example: |
| 35 | \begin{verbatim} |
| 36 | 800.,1000,10,100/ |
| 37 | \end{verbatim} |
| 38 | generates 1000 events at $E_{\rm cm} = 800\,\GeV$ and prints 10 |
| 39 | events. The events printed are: 1,100,200,\dots. Note that an event |
| 40 | typically takes several pages of output. This line is read with a list |
| 41 | directed format (READ*). |
| 42 | |
| 43 | After Nprint events have been printed, a single line containing the |
| 44 | run number, the event number, and the random number seed is printed |
| 45 | every Njump events (if Njump is nonzero). This seed can be used to start |
| 46 | a new job with the given event if in the new run NSIGMA is set equal to |
| 47 | zero: |
| 48 | \begin{verbatim} |
| 49 | SEED |
| 50 | value/ |
| 51 | NSIGMA |
| 52 | 0/ |
| 53 | \end{verbatim} |
| 54 | In general the same events will only be generated on the same type of |
| 55 | computer. |
| 56 | |
| 57 | Reaction line: This line must be given unless title is SAME, when |
| 58 | it must be omitted. It selects the type of events to be generated. The |
| 59 | present version can generate TWOJET, E+E-, DRELLYAN, MINBIAS, WPAIR, |
| 60 | SUPERSYM, HIGGS, PHOTON, TCOLOR, or WHIGGS events. This line is read |
| 61 | with an A8 format. |
| 62 | |
| 63 | \subsection{Optional Parameters} |
| 64 | |
| 65 | Each optional parameter requires two lines. |
| 66 | The first line is a keyword specifying the parameter and the second |
| 67 | line gives the values for the parameter. The parameters can be given in |
| 68 | any order. Numerical values are read with a list directed format |
| 69 | (READ*), jet and particle types are read with a character format and |
| 70 | must be enclosed in quotes, and logical flags with an L1 format. All |
| 71 | momenta are in GeV and all angles are in radians. |
| 72 | |
| 73 | The parameters can be classified in several groups: |
| 74 | \begin{center} |
| 75 | \begin{tabular}{lllll} |
| 76 | \hline\hline |
| 77 | Jet Limits: & W/H Limits: & Decays: & Constants: & Other: \\ |
| 78 | \hline |
| 79 | JETTYPE1 & HTYPE & FORCE & AMSB & BEAMS \\ |
| 80 | JETTYPE2 & PHIW & FORCE1 & CUTJET & EPOL \\ |
| 81 | JETTYPE3 & QMH & NODECAY & CUTOFF & EEBEAM \\ |
| 82 | MIJLIM & QMW & NOETA & EXTRAD & EEBREM \\ |
| 83 | MTOT & QTW & NOEVOLVE & FRAGMENT & NPOMERON \\ |
| 84 | P & THW & NOFRGMNT & GAUGINO & NSIGMA \\ |
| 85 | PHI & WTYPE & NOGRAV & GMSB & NTRIES \\ |
| 86 | PT & XW & NOPI0 & GMSB2 & PDFLIB \\ |
| 87 | TH & YW & & HMASS & SEED \\ |
| 88 | X & & & HMASSES & STRUC \\ |
| 89 | Y & & & LAMBDA & WFUDGE \\ |
| 90 | WMODE1 & & & MGVTNO & WMMODE \\ |
| 91 | WMODE2 & & & MSSMA & WPMODE \\ |
| 92 | & & & MSSMB & Z0MODE \\ |
| 93 | & & & MSSMC & \\ |
| 94 | & & & MSSMD & \\ |
| 95 | & & & MSSME & \\ |
| 96 | & & & NUSUG1 & \\ |
| 97 | & & & NUSUG2 & \\ |
| 98 | & & & NUSUG3 & \\ |
| 99 | & & & NUSUG4 & \\ |
| 100 | & & & NUSUG5 & \\ |
| 101 | & & & SIGQT & \\ |
| 102 | & & & SIN2W & \\ |
| 103 | & & & SLEPTON & \\ |
| 104 | & & & SQUARK & \\ |
| 105 | & & & SSBCSC & \\ |
| 106 | & & & SUGRA & \\ |
| 107 | & & & SUGRHN & \\ |
| 108 | & & & TCMASS & \\ |
| 109 | & & & TMASS & \\ |
| 110 | & & & WMASS & \\ |
| 111 | \hline\hline |
| 112 | \end{tabular} |
| 113 | \end{center} |
| 114 | |
| 115 | It may be helpful to know that the TWOJET, WPAIR, PHOTON, |
| 116 | SUPERSYM, and WHIGGS processes use the same controlling routines and |
| 117 | so share many of the same variables. In particular, PT limits should |
| 118 | normally be set for these processes, and JETTYPE1 and JETTYPE2 are |
| 119 | used to select the reactions. Similarly, the DRELLYAN, HIGGS, and |
| 120 | TCOLOR processes use the same controlling routines since they all |
| 121 | generate s-channel resonances. The mass limits for these processes are |
| 122 | set by QMW. Normally the QMW limits will surround the $W^\pm$, $Z^0$, |
| 123 | or Higgs mass, but this is not required. (QMH acts like QMW for the |
| 124 | Higgs process.) For historical reasons, JETTYPE1 and JETTYPE2 are used |
| 125 | to select the W decay modes in DRELLYAN, while WMODE1 and WMODE2 select |
| 126 | the W decay modes for WPAIR, HIGGS, and WHIGGS. Also, QTW can be used |
| 127 | to generate DRELLYAN events with non-zero transverse momentum, whereas |
| 128 | HIGGS automatically fixes QTW to be zero. (Of course, non-zero |
| 129 | transverse momentum will be generated by gluon radiation.) |
| 130 | |
| 131 | For example the lines |
| 132 | \begin{verbatim} |
| 133 | P |
| 134 | 40.,50.,10.,100./ |
| 135 | \end{verbatim} |
| 136 | would set limits for the momentum of jet 1 between 40 and 50 GeV, and |
| 137 | for jet 2 between 10 and 100 GeV. As another example the lines |
| 138 | \begin{verbatim} |
| 139 | WTYPE |
| 140 | 'W+'/ |
| 141 | \end{verbatim} |
| 142 | would specify that for DRELLYAN events only W+ events will be generated. |
| 143 | If for a kinematic variable only the lower limit is specified then that |
| 144 | parameter is fixed to the given value. Thus the lines |
| 145 | \begin{verbatim} |
| 146 | P |
| 147 | 40.,,10./ |
| 148 | \end{verbatim} |
| 149 | will fix the momentum for jet 1 to be 40 GeV and for jet 2 to be 10 |
| 150 | GeV. If only the upper limit is specified then the default value is used |
| 151 | for the lower limit. Jet 1 or jet 2 parameters for DRELLYAN events refer |
| 152 | to the W decay products and cannot be fixed. If QTW is fixed to 0, then |
| 153 | standard Drell-Yan events are generated. |
| 154 | |
| 155 | A complete list of keywords and their default values follows. |
| 156 | |
| 157 | \newpage |
| 158 | \begin{center} |
| 159 | \begin{tabular}{lll} |
| 160 | \hline\hline |
| 161 | Keyword & & Explanation \\ |
| 162 | Values & Default values & \\ |
| 163 | \hline |
| 164 | AMSB & & Anomaly-mediated SUSY breaking \\ |
| 165 | $m_0$,$m_{3/2}$,$\tan\beta$,$\sgn\mu$ & none & scalar mass, gravitino mass, \\ |
| 166 | & & VEV ratio, sign \\ |
| 167 | & & \\ |
| 168 | BEAMS & & Initial beams. Allowed are \\ |
| 169 | type$_1$,type$_2$ & 'P','P' & 'P','AP','N','AN'. \\ |
| 170 | & & \\ |
| 171 | CUTJET & & Cutoff mass for QCD jet \\ |
| 172 | $\mu_c$ & 6. & evolution. \\ |
| 173 | & & \\ |
| 174 | CUTOFF & & Cutoff $qt^2=\mu^2Q^\nu$ for \\ |
| 175 | $\mu^2$, $\nu$ & .200,1.0 & DRELLYAN events. \\ |
| 176 | & & \\ |
| 177 | EEBEAM & & impose brem/beamstrahlung \\ |
| 178 | $\sqrt{\hat{s}}_{min}$, $\sqrt{\hat{s}}_{max}$, $\Upsilon$, $\sigma_z$ & |
| 179 | none & min and max subprocess energy, \\ |
| 180 | & & beamstrahlung parameter $\Upsilon$ \\ |
| 181 | & & longitudinal beam size $\sigma_z$ in mm \\ |
| 182 | & & \\ |
| 183 | EEBREM & & impose bremsstrahlung for $e^+e^-$ \\ |
| 184 | $\sqrt{\hat{s}}_{min}$, $\sqrt{\hat{s}}_{max}$ & none & min and max subprocess |
| 185 | energy \\ |
| 186 | & & \\ |
| 187 | EPOL & & Polarization of $e^-$ ($e^+$) beam, \\ |
| 188 | $P_L(e^-),P_L(e^+)$ & 0,0 & $P_L(e)=(n_L-n_R)/(n_L-n_R)$, \\ |
| 189 | & & so that $-1 \le P_L \le 1$ \\ |
| 190 | & & \\ |
| 191 | EXTRAD & & Parameters for EXTRADIM process\\ |
| 192 | $\delta$,$M_D$,UVCUT & None & UVCUT is logical flag \\ |
| 193 | & & \\ |
| 194 | FORCE & & Force decay of particles, \\ |
| 195 | $i,i_1,...,i_5$/ & None & $\pm i \to \pm(i1+...+i5)$. \\ |
| 196 | & & Can call 20 times. \\ |
| 197 | & & See note for $i$ = quark. \\ |
| 198 | & & \\ |
| 199 | FORCE1 & & Force decay $i \to i1+...+i5$. \\ |
| 200 | $i,i_1,...,i_5$/ & None & Can call 40 times. \\ |
| 201 | & & See note for $i$ = quark. \\ |
| 202 | & & \\ |
| 203 | FRAGMENT & & Fragmentation parameters. \\ |
| 204 | $P_{ud}$,\dots & .4,\dots & See also SIGQT, etc. \\ |
| 205 | & & \\ |
| 206 | GAUGINO & & Masses for $\tilde g$, |
| 207 | $\tilde\gamma$, \\ |
| 208 | $m_1$,$m_2$,$m_3,m_4$ & 50,0,100,100 & $\tilde W^+$, and $\tilde Z^0$ \\ |
| 209 | \hline\hline |
| 210 | \end{tabular} |
| 211 | \end{center} |
| 212 | |
| 213 | \newpage |
| 214 | \begin{center} |
| 215 | \begin{tabular}{lll} |
| 216 | \hline\hline |
| 217 | GMSB & & GMSB messenger SUSY breaking, \\ |
| 218 | $\Lambda_m$,$M_m$,$N_5$ & none & mass, number of $5+\bar5$, VEV \\ |
| 219 | $\tan\beta$,$\sgn\mu$,$C_{\rm gr}$ & & ratio, sign, gravitino scale \\ |
| 220 | & & \\ |
| 221 | GMSB2 & & non-minimal GMSB parameters \\ |
| 222 | $\slashchar{R}$,$\delta M_{H_d}^2$,$\delta M_{H_u}^2$,$D_Y(M)$ & 1,0,0,0 & |
| 223 | gaugino mass multiplier \\ |
| 224 | $N_{5_1}$,$N_{5_2}$,$N_{5_3}$ & $N_5$ & Higgs mass shifts, D-term mass$^2$\\ |
| 225 | & & indep. gauge group messengers \\ |
| 226 | & & \\ |
| 227 | HMASS & 0 & Mass for standard Higgs. \\ |
| 228 | $m$ & & \\ |
| 229 | & & \\ |
| 230 | HMASSES & & Higgs meson masses for \\ |
| 231 | $m_1$,\dots,$m_9$ & 0,...,0 & charges 0,0,0,0,0,1,1,2,2. \\ |
| 232 | HTYPE & & One MSSM Higgs type ('HL0', \\ |
| 233 | 'HL0'/ or... & none & 'HH0', or 'HA0') \\ |
| 234 | & & \\ |
| 235 | JETTYPE1 & & )Select types for jets: \\ |
| 236 | 'GL','UP',... & 'ALL' & )'ALL'; 'GL'; 'QUARKS'='UP', \\ |
| 237 | & & )'UB','DN','DB','ST','SB', \\ |
| 238 | JETTYPE2 & & )'CH','CB','BT','BB','TP', \\ |
| 239 | 'GL','UP',... & 'ALL' & )'TB','X','XB','Y','YB'; \\ |
| 240 | & & )'LEPTONS'='E-','E+','MU-', \\ |
| 241 | JETTYPE3 & & )'MU+','TAU-','TAU+'; 'NUS'; \\ |
| 242 | 'GL','UP',... & 'ALL' & )'GM','W+','W-','Z0' \\ |
| 243 | & & ) See note for SUSY types. \\ |
| 244 | & & \\ |
| 245 | LAMBDA & & QCD scale \\ |
| 246 | $\Lambda$ & .2 & \\ |
| 247 | & & \\ |
| 248 | MGVTNO & & Gravitino mass -- ignored for \\ |
| 249 | $M_{\rm gravitino}$ & $10^{20}$~GeV & GMSB model \\ |
| 250 | & & \\ |
| 251 | MIJLIM & & Multimet mass limits \\ |
| 252 | $i$,$j$,$M_{\rm min}$,$M_{\rm max}$ & 0,0,$1\,\GeV$,$1\,\GeV$ & \\ |
| 253 | & & \\ |
| 254 | MSSMA & & MSSM parameters -- \\ |
| 255 | $m(\tilde g)$,$\mu$, & Required & Gluino mass, $\mu$, $A$ mass, \\ |
| 256 | $m(A)$,$\tan\beta$ & & $\tan\beta$ \\ |
| 257 | & & \\ |
| 258 | MSSMB & & MSSM 1st generation -- \\ |
| 259 | $m(q_1)$,$m(d_r)$,$m(u_r)$, & Required & Left and right soft squark and \\ |
| 260 | $m(l_1)$,$m(e_r)$ & & slepton masses \\ |
| 261 | \hline\hline |
| 262 | \end{tabular} |
| 263 | \end{center} |
| 264 | |
| 265 | \newpage |
| 266 | \begin{center} |
| 267 | \begin{tabular}{lll} |
| 268 | \hline\hline |
| 269 | MSSMC & & MSSM 3rd generation -- \\ |
| 270 | $m(q_3)$,$m(b_r)$,$m(t_r)$, & Required & Soft squark masses, slepton \\ |
| 271 | $m(l_3)$,$m(\tau_r)$, & & masses, and squark and slepton \\ |
| 272 | $A_t$,$A_b$,$A_\tau$ & & mixings \\ |
| 273 | & & \\ |
| 274 | MSSMD & & MSSM 2nd generation -- \\ |
| 275 | $m(q_2)$,$m(s_r)$,$m(c_r)$, & from MSSMB & Left and right soft squark and \\ |
| 276 | $m(l_2)$,$m(mu_r)$ & & slepton masses \\ |
| 277 | & & \\ |
| 278 | MSSME & & MSSM gaugino masses -- \\ |
| 279 | $M_1$,$M_2$ & MSSMA + GUT & Default is to scale from gluino\\ |
| 280 | & & \\ |
| 281 | MTOT & & Mass range for multiparton \\ |
| 282 | $M_{\rm min}$,$M_{\rm max}$ & None & processes \\ |
| 283 | & & \\ |
| 284 | NODECAY & & Suppress all decays. \\ |
| 285 | TRUE or FALSE & FALSE & \\ |
| 286 | & & \\ |
| 287 | NOETA & & Suppress eta decays. \\ |
| 288 | TRUE or FALSE & FALSE & \\ |
| 289 | |
| 290 | NOEVOLVE & & Suppress QCD evolution and \\ |
| 291 | TRUE or FALSE & FALSE & hadronization. \\ |
| 292 | & & \\ |
| 293 | NOGRAV & & Suppress gravitino decays in \\ |
| 294 | TRUE or FALSE & FALSE & GMSB model \\ |
| 295 | & & \\ |
| 296 | NOHADRON & & Suppress hadronization of \\ |
| 297 | TRUE or FALSE & FALSE & jets and beam jets. \\ |
| 298 | & & \\ |
| 299 | NONUNU & & Suppress $Z^0$ neutrino decays.\\ |
| 300 | TRUE or FALSE & FALSE & \\ |
| 301 | & & \\ |
| 302 | NOPI0 & &Suppress $\pi^0$ decays. \\ |
| 303 | TRUE or FALSE & FALSE & \\ |
| 304 | & & \\ |
| 305 | NPOMERON & & Allow $n_1<n<n_2$ cut pomerons.\\ |
| 306 | $n_1$,$n_2$ & 1,20 & Controls beam jet mult. \\ |
| 307 | & & \\ |
| 308 | NSIGMA & & Generate n unevolved events \\ |
| 309 | $n$ & 20 & for SIGF calculation. \\ |
| 310 | & & \\ |
| 311 | NTRIES & & Stop if after n tries \\ |
| 312 | $n$ & 1000 & cannot find a good event. \\ |
| 313 | \hline\hline |
| 314 | \end{tabular} |
| 315 | \end{center} |
| 316 | |
| 317 | \newpage |
| 318 | \begin{center} |
| 319 | \begin{tabular}{lll} |
| 320 | \hline\hline |
| 321 | NUSUG1 & & Optional non-universal SUGRA \\ |
| 322 | $M_1$,$M_2$,$M_3$ & none & gaugino masses \\ |
| 323 | & & \\ |
| 324 | NUSUG2 & & Optional non-universal SUGRA \\ |
| 325 | $A_t$,$A_b$,$A_\tau$ & none & $A$ terms \\ |
| 326 | & & \\ |
| 327 | NUSUG3 & & Optional non-universal SUGRA \\ |
| 328 | $M_{H_d}$,$M_{H_u}$ & none & Higgs masses \\ |
| 329 | & & \\ |
| 330 | NUSUG4 & & Optional non-universal SUGRA \\ |
| 331 | $M_{u_L}$,$M_{d_R}$,$M_{u_R}$, & none & 1st/2nd generation masses \\ |
| 332 | $M_{e_L}$,$M_{e_R}$ & & \\ |
| 333 | & & \\ |
| 334 | NUSUG5 & & Optional non-universal SUGRA \\ |
| 335 | $M_{t_L}$,$M_{b_R}$,$M_{t_R}$, & none & 3rd generation masses \\ |
| 336 | $M_{\tau_L}$,$M_{\tau_R}$ & & \\ |
| 337 | & & \\ |
| 338 | P & & Momentum limits for jets. \\ |
| 339 | $p_{\rm min}(1)$,\dots,$p_{\rm max}(3)$ & |
| 340 | 1.,$0.5E_{\rm cm}$ & \\ |
| 341 | & & \\ |
| 342 | PDFLIB & & CERN PDFLIB parton distribution\\ |
| 343 | 'name$_1$',val$_1$,\dots & None & parameters. See PDFLIB manual. \\ |
| 344 | & & \\ |
| 345 | PHI & & Phi limits for jets. \\ |
| 346 | $\phi_{\rm min}(1)$,\dots,$\phi_{\rm max}(3)$ & 0,$2\pi$ & \\ |
| 347 | & & \\ |
| 348 | PHIW & & Phi limits for W. \\ |
| 349 | $\phi_{\rm min}$,$\phi_{\rm max}$ & |
| 350 | 0,$2\pi$ & \\ |
| 351 | & & \\ |
| 352 | PT or PPERP & & $p_t$ limits for jets. \\ |
| 353 | $p_{t,{\rm min}}(1)$,\dots,$p_{t,{\rm max}}(3)$ & |
| 354 | $.05E_{\rm cm}$,$.2E_{\rm cm}$ & Default for TWOJET only. \\ |
| 355 | & & \\ |
| 356 | QMH & & Mass limits for Higgs. \\ |
| 357 | $q_{\rm min}$,$q_{\rm max}$ & |
| 358 | $.05E_{\rm cm}$,$.2E_{\rm cm}$ & Equivalent to QMW. \\ |
| 359 | & & \\ |
| 360 | QMW & & Mass limits for $W$. \\ |
| 361 | $q_{\rm min}$,$q_{\rm max}$ & |
| 362 | $.05E_{\rm cm}$,$.2E_{\rm cm}$ & \\ |
| 363 | & & \\ |
| 364 | QTW & & $q_t$ limits for $W$. Fix |
| 365 | $q_t=0$ \\ |
| 366 | $q_{t,{\rm min}}$,$q_{t,{\rm max}}$ & |
| 367 | .1,$.025E_{\rm cm}$ & for standard Drell-Yan. \\ |
| 368 | & & \\ |
| 369 | SEED & & Random number seed (double \\ |
| 370 | real & 0 & precision if 32 bit). \\ |
| 371 | \hline\hline |
| 372 | \end{tabular} |
| 373 | \end{center} |
| 374 | |
| 375 | \newpage |
| 376 | \begin{center} |
| 377 | \begin{tabular}{lll} |
| 378 | \hline\hline |
| 379 | SIGQT & & Internal $k_t$ parameter for \\ |
| 380 | $\sigma$ & .35 & jet fragmentation. \\ |
| 381 | & & \\ |
| 382 | SIN2W & & Weinberg angle. See WMASS. \\ |
| 383 | $\sin^2(\theta_W)$ & .232 & \\ |
| 384 | & & \\ |
| 385 | SLEPTON & & Masses for $\tilde \nu_e$, |
| 386 | $\tilde e$, $\tilde\nu_\mu$, $\tilde\mu$, $\tilde\nu_\tau$, $\tilde\tau$ \\ |
| 387 | $m_1$,\dots,$m_6$ & 100,\dots,101.8 & \\ |
| 388 | & & \\ |
| 389 | SQUARK & & Masses for $\tilde u$, |
| 390 | $\tilde d$, $\tilde s$, $\tilde c$, $\tilde b$, $\tilde t$ \\ |
| 391 | $m_1$,\dots,$m_6$ & 100.3,...,240. & \\ |
| 392 | & & \\ |
| 393 | SSBCSC & & Alternate mass scale for RGE \\ |
| 394 | $M$ & $M_{GUT}$ & boundary conditions. \\ |
| 395 | & & \\ |
| 396 | STRUC & & Structure functions. CTEQ3L, \\ |
| 397 | name & 'CTEQ3L' & CTEQ2L, EHLQ, OR DO \\ |
| 398 | & & \\ |
| 399 | SUGRA & & Minimal supergravity parameters\\ |
| 400 | $m_0$,$m_{1/2}$,$A_0$, & none & scalar M, gaugino M, trilinear \\ |
| 401 | $\tan\beta$,$\sgn\mu$ & & breaking term, vev ratio, +-1 \\ |
| 402 | TH or THETA & & Theta limits for jets. Do not \\ |
| 403 | $\theta_{\rm min}(1)$,\dots,$\theta_{\rm max}(3)$ & 0,$\pi$ & also set Y. \\ |
| 404 | & & \\ |
| 405 | SUGRHN & & SUGRA see-saw $\nu$-effect \\ |
| 406 | $m_{\nu_\tau}$,$M_N$,$A_n$,$m_{\tilde\nu_R}$ & $0,1E20,0,0$ & nu-mass, |
| 407 | int. scale, \\ |
| 408 | & & GUT scale nu SSB terms \\ |
| 409 | & & \\ |
| 410 | THW & & Theta limits for W. Do not \\ |
| 411 | $\theta_{\rm min}$,$\theta_{\rm max}$ & 0,$\pi$ & also set YW. \\ |
| 412 | & & \\ |
| 413 | TCMASS & & Technicolor mass and width. \\ |
| 414 | $m$,$\Gamma$ & 1000,100 & \\ |
| 415 | & & \\ |
| 416 | TMASS & & t, y, and x quark masses. \\ |
| 417 | $m_t$,$m_y$,$m_x$ & 180.,-1.,-1. & \\ |
| 418 | & & \\ |
| 419 | WFUDGE & & Fudge factor for DRELLYAN \\ |
| 420 | factor & 1.85 & evolution scale. \\ |
| 421 | & & \\ |
| 422 | WMASS & & W and Z masses. See SIN2W. \\ |
| 423 | $M_W$,$M_Z$ & 80.2, 91.19 & \\ |
| 424 | \hline\hline |
| 425 | \end{tabular} |
| 426 | \end{center} |
| 427 | |
| 428 | \newpage |
| 429 | \begin{center} |
| 430 | \begin{tabular}{lll} |
| 431 | \hline\hline |
| 432 | WMMODE & & Decay modes for $W^-$ in parton\\ |
| 433 | 'UP',\dots,'TAU+' & 'ALL' & cascade. See JETTYPE. \\ |
| 434 | & & \\ |
| 435 | WMODE1 & & ) \\ |
| 436 | 'UP','UB',\dots & 'ALL' & )Decay modes for WPAIR. \\ |
| 437 | & & )Same code for quarks and \\ |
| 438 | WMODE2 & & )leptons as JETTYPE. \\ |
| 439 | 'UP','UB',\dots & 'ALL' & ) \\ |
| 440 | & & \\ |
| 441 | WPMODE & & Decay modes for $W^+$ in parton\\ |
| 442 | 'UP',\dots,'TAU+' & 'ALL' & cascade. See JETTYPE. \\ |
| 443 | & & \\ |
| 444 | WTYPE & & Select W type: W+,W-,GM,Z0. \\ |
| 445 | type$_1$,type$_2$ & 'GM','Z0' & Do not mix W+,W- and GM,Z0. \\ |
| 446 | & & \\ |
| 447 | X & & Feynman x limits for jets. \\ |
| 448 | $x_{\rm min}(1)$,\dots,$x_{\rm max}(3)$ & |
| 449 | $-1$,1 & \\ |
| 450 | & & \\ |
| 451 | XGEN & & Jet fragmentation, Peterson \\ |
| 452 | a(1),\dots,a(8) & .96,3,0,.8,.5,... & with $\epsilon=a(n)/m^2$, |
| 453 | $n=4$-8. \\ |
| 454 | & & \\ |
| 455 | XGENSS & & Fragmentation of GLSS, UPSS, \\ |
| 456 | a(1),\dots,a(7) & .5,.5,... & etc. with $\epsilon=a(n)/m**2$ \\ |
| 457 | & & \\ |
| 458 | XW & & Feynman x limits for W. \\ |
| 459 | $x_{\rm min}$,$x_{\rm max}$ & |
| 460 | $-1$,1 & \\ |
| 461 | & & \\ |
| 462 | Y & & Y limits for each jet. \\ |
| 463 | $y_{\rm min}(1)$,\dots,$y_{\rm max}(3)$ & from PT & Do not also set TH. \\ |
| 464 | & & \\ |
| 465 | YW & & Y limits for W. \\ |
| 466 | $y_{\rm min}$,$y_{\rm max}$ & from QTW,QMW & Do not set both YW and THW. \\ |
| 467 | & & \\ |
| 468 | Z0MODE & & Decay modes for $Z^0$ in parton\\ |
| 469 | 'UP',\dots,'TAU+' & 'ALL' & cascade. See JETTYPE. \\ |
| 470 | \hline\hline |
| 471 | \end{tabular} |
| 472 | \end{center} |
| 473 | |
| 474 | \newpage |
| 475 | \subsection{Kinematic and Parton-type Parameters} |
| 476 | |
| 477 | While the TWOJET PT limits and the DRELLYAN QMW limits are |
| 478 | formally optional parameters, they are set by default to be fractions of |
| 479 | $\sqrt{s}$. Thus, for example, the parameter file |
| 480 | \begin{verbatim} |
| 481 | DEFAULT TWOJET JOB |
| 482 | 14000,100,1,100/ |
| 483 | TWOJET |
| 484 | END |
| 485 | STOP |
| 486 | \end{verbatim} |
| 487 | will execute, but it will generate jets between 5\% and 20\% of |
| 488 | $\sqrt{s}$, which is probably not what is wanted. Similarly, the |
| 489 | parameter file |
| 490 | \begin{verbatim} |
| 491 | DEFAULT DRELLYAN JOB |
| 492 | 14000,100,1,100/ |
| 493 | DRELLYAN |
| 494 | END |
| 495 | STOP |
| 496 | \end{verbatim} |
| 497 | will generate $\gamma + Z$ events with masses between 5\% and 20\% of |
| 498 | $\sqrt{s}$, not masses around the $Z$ mass, and transverse momenta |
| 499 | between $1\,{\rm GeV}$ and 2.5\% of $\sqrt{s}$. |
| 500 | |
| 501 | Normally the user should set PT limits for TWOJET, PHOTON, WPAIR, |
| 502 | SUPERSYM, and WHIGGS events and QMW and QTW limits for DRELLYAN, |
| 503 | HIGGS, and TCOLOR events. If these limits are not set, they will be |
| 504 | selected as fractions of $E_{\rm cm}$. This can give nonsense. For |
| 505 | TWOJET the $p_t$ range should usually be less than about a factor of |
| 506 | two except for $b$ and $t$ jets at low $p_t$ to produce uniform |
| 507 | statistics. For $W^+$, $W^-$, or $Z^0$ events or for Higgs events the |
| 508 | QMW (QMH) range should usually include the mass. But one can select |
| 509 | different limits to study, e.g., virtual $W$ production or the effect |
| 510 | of a lighter or heavier Higgs on WW scattering. If only $t$ decays are |
| 511 | selected, then the lower QMW limit must be above the $t$ threshold. |
| 512 | For standard Drell-Yan events QTW should be fixed to zero, |
| 513 | \begin{verbatim} |
| 514 | QTW |
| 515 | 0/ |
| 516 | \end{verbatim} |
| 517 | Transverse momenta will then be generated by initial state gluon |
| 518 | radiation. A range of QTW can also be given. For SUPERSYM either the |
| 519 | masses and decay modes should be specified, or the MSSM, SUGRA, GMSB, or |
| 520 | AMSB parameters should be given. For fourth generation quarks it is |
| 521 | necessary to specify the quark masses. |
| 522 | |
| 523 | Note that if the limits given cover too large a kinematic range, |
| 524 | the program can become very inefficient, since it makes a fit to the |
| 525 | cross section over the specified range. NTRIES has to be increased if |
| 526 | narrow limits are set for X, XW or for jet 1 and jet 2 parameters in |
| 527 | DRELLYAN events. For larger ranges several runs can be combined together |
| 528 | using the integrated cross section per event SIGF/NEVENT as the weight. |
| 529 | This cross section is calculated for each run by Monte Carlo integration |
| 530 | over the specified kinematic limits and is printed at the end of the |
| 531 | run. It is corrected for JETTYPEi, WTYPE, and WMODEi selections; it |
| 532 | cannot be corrected for branching ratios of forced decays or for WPMODE, |
| 533 | WMMODE, or Z0MODE selections, since these can affect an arbitrary number |
| 534 | of particles. |
| 535 | |
| 536 | To generate events over a large range, it is much more efficient |
| 537 | to combine several runs. This is facilitated by using the special job |
| 538 | title SAME as described above. Note that SAME cannot be used to combine |
| 539 | standard DRELLYAN events (QTW fixed equal to 0) and DRELLYAN events with |
| 540 | nonzero QTW. |
| 541 | |
| 542 | The cross sections for multiparton final states in general have |
| 543 | infrared and collinear singularities. To obtain sensible results, it |
| 544 | is in general essential to set limits both on the $p_T$ of each final |
| 545 | parton using PT and on the mass of each pair of partons using MIJLIM. |
| 546 | The default lower limits are all $1\,{\rm GeV}$. Using these default |
| 547 | limits without thought will likely give absurd results. |
| 548 | |
| 549 | For TWOJET, DRELLYAN, and most other processes, the JETTYPEi and |
| 550 | WTYPEi keywords should be used to select the subprocesses to be |
| 551 | included. For $e^+ e^- \to W^+ W^-$, $Z^0 Z^0$, use FORCE and FORCE1 |
| 552 | instead of WMODEi to select the $W$ decay modes. Note that these {\it |
| 553 | do not} change the calculated cross section. (In the E+E- process, the |
| 554 | $W$ and $Z$ decays are currently treated as particle decays, whereas in |
| 555 | the WPAIR and HIGGS processes they are treated as $2 \to 4$ parton |
| 556 | processes.) |
| 557 | |
| 558 | For HIGGS with $W^+W^-$ or $Z^0Z^0$ decays allowed it is |
| 559 | generally necessary to set PT limits for the W's, e.g. |
| 560 | \begin{verbatim} |
| 561 | PT |
| 562 | 50,20000,50,20000/ |
| 563 | \end{verbatim} |
| 564 | If this is not done, then the default lower limit of 1 GeV is used, |
| 565 | and the $t$-channel exchanges will dominate, as they should in the |
| 566 | effective $W$ approximation. Depending on the other parameters, the |
| 567 | program may fail to generate an event in NTRIES tries. |
| 568 | |
| 569 | \subsection{SUSY Parameters} |
| 570 | |
| 571 | SUPERSYM (SUSY) by default generates just gluinos and squarks in |
| 572 | pairs. There are no default masses or decay modes. Masses can be set |
| 573 | using GAUGINO, SQUARK, SLEPTON, and HMASSES. Decay modes can be |
| 574 | specified with FORCE or by modifying the decay table. Left and right |
| 575 | squarks are distinguished but assumed to be degenerate, except for |
| 576 | stops. Since version 7.11, types must be selected with JETTYPEi using |
| 577 | the supersymmetric names, e.g. |
| 578 | \begin{verbatim} |
| 579 | JETTYPE1 |
| 580 | 'GLSS','UPSSL','UPSSR'/ |
| 581 | \end{verbatim} |
| 582 | Use of the corresponding standard model names, e.g. |
| 583 | \begin{verbatim} |
| 584 | JETTYPE1 |
| 585 | 'GL','UP'/ |
| 586 | \end{verbatim} |
| 587 | and generation of pure photinos, winos, and zinos are no longer |
| 588 | supported. |
| 589 | |
| 590 | If MSSMA, MSSMB and MSSMC are given, then the specified parameters |
| 591 | are used to calculate all the masses and decay modes with the ISASUSY |
| 592 | package assuming the minimal supersymmetric extension of the standard |
| 593 | model (MSSM). There are no default values, so you must specify values |
| 594 | for each MSSMi, i=A-C. MSSMD can optionally be used to set the second |
| 595 | generation squark and slepton parameters; if it is omitted, then the |
| 596 | first generation ones are used. MSSME can optionally be used to set the |
| 597 | U(1) and SU(2) gaugino masses; if it is omitted, then the grand |
| 598 | unification values are used. The parameters and the use of the MSSM is |
| 599 | preserved if the title is SAME. FORCE can be used to override the |
| 600 | calculated branching ratios. |
| 601 | |
| 602 | The MSSM option also generates charginos and neutralinos with |
| 603 | cross sections based on the MSSM mixing angles in addition to squarks |
| 604 | and sleptons. These can be selected with JETTYPEi; the complete list of |
| 605 | supersymmetric options is: |
| 606 | \begin{verbatim} |
| 607 | 'GLSS', |
| 608 | 'UPSSL','UBSSL','DNSSL','DBSSL','STSSL','SBSSL','CHSSL','CBSSL', |
| 609 | 'BTSS1','BBSS1','TPSS1','TBSS1', |
| 610 | 'UPSSR','UBSSR','DNSSR','DBSSR','STSSR','SBSSR','CHSSR','CBSSR', |
| 611 | 'BTSS2','BBSS2','TPSS2','TBSS2', |
| 612 | 'W1SS+','W1SS-','W2SS+','W2SS-','Z1SS','Z2SS','Z3SS','Z4SS', |
| 613 | 'NUEL','ANUEL','EL-','EL+','NUML','ANUML',MUL-','MUL+','NUTL', |
| 614 | 'ANUTL','TAU1-','TAU1+','ER-','ER+','MUR-','MUR+','TAU2-','TAU2+', |
| 615 | 'Z0','HL0','HH0','HA0','H+','H-', |
| 616 | 'SQUARKS','GAUGINOS','SLEPTONS','ALL'. |
| 617 | \end{verbatim} |
| 618 | Note that mixing between $L$ and $R$ stop states results in 1 (light) |
| 619 | and 2 (heavy) stop, sbottom and stau eigenstates, which depend on the |
| 620 | input parameters of left- and right- scalar masses, plus $A$ terms, |
| 621 | $\mu$ and $\tan\beta$. The last four JETTYPE's generate respectively |
| 622 | all allowed combinations of squarks and antisquarks, all combinations |
| 623 | of charginos and neutralinos, all combinations of sleptons and |
| 624 | sneutrinos, and all SUSY particles. |
| 625 | |
| 626 | For SUSY Higgs pair production or associated production in E+E-, |
| 627 | select the appropriate JETTYPE's, e.g. |
| 628 | \begin{verbatim} |
| 629 | JETTYPE1 |
| 630 | 'Z0'/ |
| 631 | JETTYPE2 |
| 632 | 'HL0'/ |
| 633 | \end{verbatim} |
| 634 | |
| 635 | As usual, this gives only half the cross section. For single production |
| 636 | of neutral SUSY Higgs in $pp$ and $\bar pp$ reactions, use the HIGGS |
| 637 | process together with the MSSMi, SUGRA, GMSB, or AMSB keywords. You must |
| 638 | specify one and only one Higgs type using |
| 639 | \begin{verbatim} |
| 640 | HTYPE |
| 641 | 'HL0' or 'HH0' or 'HA0'/ <<<<< One only! |
| 642 | \end{verbatim} |
| 643 | If no QMH range is given, one is calculated using $M \pm 5 \Gamma$ for |
| 644 | the selected Higgs. Decays into quarks, leptons, gauge bosons, lighter |
| 645 | Higgs bosons, and SUSY particles are generated using the on-shell |
| 646 | branching ratios from ISASUSY. You can use JETTYPEi to select the |
| 647 | allowed Higgs modes and WMODEi to select the allowed decays of W and Z |
| 648 | bosons. Since heavy SUSY Higgs bosons couple weakly to W pairs, WW |
| 649 | fusion and WW scattering are not included. |
| 650 | |
| 651 | SUGRA can be used instead of MSSMi to generate MSSM decays with |
| 652 | parameters determined from $m_0$, $m_{1/2}$, $A_0$, $\tan\beta$, and |
| 653 | $\sgn\mu=\pm1$ in the minimal supergravity framework. The NUSUGi |
| 654 | keywords can optionally be used to specify additional parameters for |
| 655 | non-universal SUGRA models, while SUGRHN is used to specify the |
| 656 | parameterf of an optional right-handed neutrino. Similarly, the GMSB |
| 657 | keyword is used to specify the $\Lambda$, $M_m$, $N_5$, $\tan\beta$, |
| 658 | $\sgn\mu=\pm1$, and $C_{\rm grav}$ parameters of the minimal Gauge |
| 659 | Mediated SUSY Breaking model. GMSB2 can optionally be used to specify |
| 660 | additional parameters of non-minimal GMSB models. The AMSB keyword is |
| 661 | used to specify $m_0$, $m_{3/2}$, $\tan\beta$, and $\sgn\mu$ for the |
| 662 | minimal Anomaly Mediated SUSY Breaking model. Note that $m_{3/2}$ is |
| 663 | much larger than the weak scale, typically 50~TeV. |
| 664 | |
| 665 | WHIGGS is used to generate $W$ plus neutral Higgs events. For the |
| 666 | Standard Model the JETTYPE is \verb|HIGGS|. If any of the SUSY models |
| 667 | is specified, then the appropriate SUSY Higgs type should be used, |
| 668 | most likely \verb|HL0|. In either case WMODEi is used to specify the |
| 669 | $W$ decay modes. The Higgs is treated as a particle; its decay modes |
| 670 | can be set using FORCE. |
| 671 | |
| 672 | \subsection{Forced Decay Modes} |
| 673 | |
| 674 | The FORCE keyword requires special care. Its list must contain the |
| 675 | numerical particle IDENT codes, e.g. |
| 676 | \begin{verbatim} |
| 677 | FORCE |
| 678 | 140,130,-120/ |
| 679 | \end{verbatim} |
| 680 | The charge-conjugate mode is also forced for its antiparticle. Thus the |
| 681 | above example forces both $\bar D^0 \to K^+ \pi^-$ and $D^0 \to K^- |
| 682 | \pi^+$. If only a specific decay is wanted one should use the FORCE1 |
| 683 | command; e.g. |
| 684 | \begin{verbatim} |
| 685 | FORCE1 |
| 686 | 140,130,-120/ |
| 687 | \end{verbatim} |
| 688 | only forces $\bar D^0 \to K^+ \pi^-$. |
| 689 | |
| 690 | To force a heavy quark decay one must generally separately force |
| 691 | each hadron containing it. If the decay is into three leptons or quarks, |
| 692 | then the real or virtual W propagator is inserted automatically. Since |
| 693 | Version 7.30, top and fourth generation quarks are treated as |
| 694 | particles and decayed directly rather than first being made into |
| 695 | hadrons. Thus for example |
| 696 | \begin{verbatim} |
| 697 | FORCE1 |
| 698 | 6,-12,11,5/ |
| 699 | \end{verbatim} |
| 700 | forces all top quarks to decay into an positron, neutrino and a |
| 701 | b-quark (which will be hadronized). For the physical top mass, the |
| 702 | positron and neutrino will come from a real W. Note that forcing $t |
| 703 | \to W^+ b$ and $W^+ \to e^+ \nu_e$ does {\it not} give the same |
| 704 | result; the first uses the correct $V-A$ matrix element, while the |
| 705 | second decays the $W$ according to phase space. |
| 706 | |
| 707 | Forced modes included in the decay table or generated by ISASUSY |
| 708 | will automatically be put into the correct order and will use the |
| 709 | correct matrix element. Modes not listed in the decay table are |
| 710 | allowed, but caution is advised because a wrong decay mode can cause |
| 711 | an infinite loop or other unexpected effects. |
| 712 | |
| 713 | FORCE (FORCE1) can be called at most 20 (40) times in any run plus |
| 714 | all subsequent 'SAME' runs. If it is called more than once for a given |
| 715 | parent, all calls are listed, and the last call is used. Note that FORCE |
| 716 | applies to particles only, but that for gamma, W+, W-, Z0 and |
| 717 | supersymmetric particles the same IDENT codes are used both as jet types |
| 718 | and as particles. |
| 719 | |
| 720 | \subsection{Parton Distributions} |
| 721 | |
| 722 | The default parton distributions are fit CTEQ3L from the CTEQ |
| 723 | Collaboration using lowest order QCD. The CTEQ and the older EHLQ and |
| 724 | Duke-Owens distributions can be selected using the STRUC keyword. |
| 725 | |
| 726 | If PDFLIB support is enabled (see Section 4), then any of the |
| 727 | distributions in the PDFLIB compilation by H. Plothow-Besch can be |
| 728 | selected using the PDFLIB keyword and giving the proper parameters, |
| 729 | which are identical to those described in the PDFLIB manual and are |
| 730 | simply passed to the routine PDFSET. For example, to select fit 29 |
| 731 | (CTEQ3L) by the CTEQ group, leaving all other parameters with their |
| 732 | default values, use |
| 733 | \begin{verbatim} |
| 734 | PDFLIB |
| 735 | 'CTEQ',29D0/ |
| 736 | \end{verbatim} |
| 737 | Note that the fit-number and the other parameters are of type DOUBLE |
| 738 | PRECISION (REAL on 64-bit machines). There is no internal passing of |
| 739 | parameters except for those which control the printing of messages. |
| 740 | |
| 741 | \subsection{Multiparton Processes} |
| 742 | |
| 743 | For multiparton final states one should in general set limits |
| 744 | on the total mass \verb|MTOT| of the final state, on the minimum |
| 745 | \verb|PT| of each light parton, and on the minimum mass \verb|MIMLIM| |
| 746 | of each pair of light partons. Limits for \verb|PT| are set in the |
| 747 | ususal way. Limits for the mass $M_{ij}$ of partons $i,j$ are set using |
| 748 | \begin{verbatim} |
| 749 | MIJLIM |
| 750 | i,j,Mmin,Mmax |
| 751 | \end{verbatim} |
| 752 | If $i=j=0$, the limit is applied to all jet pairs. For example the |
| 753 | following parameter file generates \verb|ZJJ| events at the LHC with a |
| 754 | mimimum $p_T$ of $20\,\GeV$ and a minimum mass of $20\,\GeV$ for all |
| 755 | jet pairs: |
| 756 | \begin{verbatim} |
| 757 | GENERATE ZJJ with PTMIN = 20 GEV AND MMIN = 20 GEV |
| 758 | 14000,100,1,100/ |
| 759 | ZJJ |
| 760 | PT |
| 761 | 20,7000,20,7000,20,7000/ |
| 762 | MIJLIM |
| 763 | 0,0,20,7000/ |
| 764 | MTOT |
| 765 | 100,500/ |
| 766 | NSIGMA |
| 767 | 200/ |
| 768 | NTRIES |
| 769 | 10000/ |
| 770 | END |
| 771 | STOP |
| 772 | \end{verbatim} |
| 773 | The default lower limits for \verb|PT| and \verb|MIJLIM| are |
| 774 | $1\,\GeV$. While these limits are sufficient to make the cross |
| 775 | sections finite, they will in general not give physically sensible |
| 776 | results. Thus, {\it the user must think carefully about what limits |
| 777 | should be set.} |
git://git.uio.no / u / mrichter / AliRoot.git / blame