TIsajet/Icommon.h

   1 #ifndef ROOT_ICommon
   2 #define ROOT_ICommon
   3
   4 #ifndef __CFORTRAN_LOADED
   5 //*KEEP,cfortran.
   6 #include "cfortran.h"
   7 //*KEND.
   8 #endif
   9
  10 extern "C" {
  11
  12 /*========================================================================*/
  13 /* COMMON/DYLIM/QMIN,QMAX,QTMIN,QTMAX,YWMIN,YWMAX,XWMIN,XWMAX,THWMIN,     */
  14 /*     &        THWMAX,PHWMIN,PHWMAX,SETLMQ(12)                           */
  15 /*------------------------------------------------------------------------*/
  16 typedef struct {
  17     Float_t qmin;
  18     Float_t qmax;
  19     Float_t qtmin;
  20     Float_t qtmax;
  21     Float_t ywmin;
  22     Float_t ywmax;
  23     Float_t xwmin;
  24     Float_t xwmax;
  25     Float_t thwmin;
  26     Float_t thwmax;
  27     Float_t phwmin;
  28     Float_t phwmax;
  29     Bool_t setlmq[12];
  30     Bool_t ywset;
  31     Bool_t thwset;
  32 } DylimCommon;
  33
  34 #define DYLIM COMMON_BLOCK(DYLIM, dylim)
  35 COMMON_BLOCK_DEF(DylimCommon,DYLIM);
  36
  37 /**************************************************************************/
  38 /*           D E S C R I P T I O N :                                      */
  39 /*------------------------------------------------------------------------*/
  40 /* qmin, qmax : Mass limits for W                                         */
  41 /* qtmin, qtmax : q-subscript-t limits for W                              */
  42 /* ywmin, ywmax : Y limits for W. Cannot be set simultaneously with       */
  43 /*                thwmin, thwmax.                                         */
  44 /* xwmin, xwmax : Feynman x limit for W                                   */
  45 /* thwmin, thwmax : Theta limits for W. Cannot be set simultaneously      */
  46 /*                  with ywmin, ywmax.                                    */
  47 /* phwmin, phwmax : Phi limits for W                                      */
  48 /* setlmq :                                                               */
  49 /* ywset, thwset : Logical flags internal to the interface.               */
  50 /**************************************************************************/
  51
  52
  53
  54 /*========================================================================*/
  55 /* COMMON/EEPAR/SGMXEE,PLEP,PLEM,RSHMIN,RSHMAX,UPSLON,SIGZ,IBREM,IBEAM    */
  56 /*------------------------------------------------------------------------*/
  57 typedef struct {
  58     Float_t sgmxee;
  59     Float_t plep;
  60     Float_t plem;
  61     Float_t rshmin;
  62     Float_t rshmax;
  63     Float_t upslon;
  64     Float_t sigz;
  65     Bool_t ibrem;
  66     Bool_t ibeam;
  67 } EeparCommon;
  68
  69 #define EEPAR COMMON_BLOCK(EEPAR, eepar)
  70 COMMON_BLOCK_DEF(EeparCommon,EEPAR);
  71
  72 /**************************************************************************/
  73 /*           D E S C R I P T I O N :                                      */
  74 /*------------------------------------------------------------------------*/
  75 /* sgmxee                                                                 */
  76 /* plep : Polarisation of positron beam (P-subscript-L of e-superscript-+ */
  77 /* plem : Polarisation of electron beam (P-subscript-L of e-superscript-- */
  78 /*                                                                        */
  79 /* Brem- and beam-strahlung parameters                                    */
  80 /* rshmin : minimum subprocess energy (root-s-hyphen-subscript-min)       */
  81 /* rshmax : maximum subprocess energy (root-s-hyphen-subscript-max)       */
  82 /* upslon : beamstrahlung parameter (UPSILON)                             */
  83 /* sigz : Longitudinal beam size in mm (sigma-subscript-z)                */
  84 /* ibrem : True if EEBREM used                                            */
  85 /* ibeam : True if EEBEAM used                                            */
  86 /**************************************************************************/
  87
  88
  89
  90 /*========================================================================*/
  91 /* COMMON/FORCE/NFORCE,IFORCE(MXFORC),MFORCE(5,MXFORC),LOOK2(2,MXFORC),   */
  92 /*     &        LOOKST(MXFORC),MEFORC(MXFORC)                             */
  93 /*------------------------------------------------------------------------*/
  94 typedef struct {
  95     static const Int_t mxforc = 40;
  96     Int_t nforce;
  97     Int_t iforce[mxforc];
  98     Int_t mforce[mxforc][5];
  99     Int_t look2[mxforc][2];
 100     Int_t lookst[mxforc];
 101     Int_t meforc[mxforc];
 102 } ForceCommon;
 103
 104 #define FORCE COMMON_BLOCK(FORCE, force)
 105 COMMON_BLOCK_DEF(ForceCommon,FORCE);
 106
 107 /**************************************************************************/
 108 /*           D E S C R I P T I O N :                                      */
 109 /*------------------------------------------------------------------------*/
 110 /* mxforce : Size of forced-decay array                                   */
 111 /* nforce : Number of forced decay paths                                  */
 112 /* iforce : ID code for particle of forced decay                          */
 113 /* mforce : Mode(s) of forced decay                                       */
 114 /* look2                                                                  */
 115 /* lookst                                                                 */
 116 /* meforc : Matrix element of the forced decay                            */
 117 /**************************************************************************/
 118
 119
 120
 121 /*========================================================================*/
 122 /* COMMON/FRGPAR/PUD,PBARY,SIGQT,PEND,XGEN(8),PSPIN1(8),PMIX1(3,2),       */
 123 /*     &         PMIX2(3,2),XGENSS(9)                                     */
 124 /*------------------------------------------------------------------------*/
 125 typedef struct {
 126     Float_t pud;
 127     Float_t pbary;
 128     Float_t sigqt;
 129     Float_t pend;
 130     Float_t xgen[8];
 131     Float_t pspin1[8];
 132     Float_t pmix1[2][3];
 133     Float_t pmix2[2][3];
 134     Float_t xgenss[9];
 135     Float_t *pmixx1[6];
 136     Float_t *pmixx2[6];
 137     Float_t *frpar[32];
 138 } FrgparCommon;
 139
 140 #define FRGPAR COMMON_BLOCK(FRGPAR, frgpar)
 141 COMMON_BLOCK_DEF(FrgparCommon,FRGPAR);
 142
 143 /**************************************************************************/
 144 /*           D E S C R I P T I O N :                                      */
 145 /*------------------------------------------------------------------------*/
 146 /* Fragmentation parameters :                                             */
 147 /* pud :                                                                  */
 148 /* pbary :                                                                */
 149 /* sigqt : Internal k-subscript-t parameter for jet fragmentation (sigma) */
 150 /* pend :                                                                 */
 151 /* xgen : Jet fragmentation, Peterson with epsilon = a(n) / m, n = 4-8    */
 152 /* pspin1 :                                                               */
 153 /* pmix1 :                                                                */
 154 /* pmix2 :                                                                */
 155 /* xgenss : Fragmentation of GLSS, UPSS etc with epsilon = a(n)/m-squared */
 156 /**************************************************************************/
 157
 158
 159
 160 /*========================================================================*/
 161 /* COMMON/HCON/ANWWWW(4,4,4),ADWWWW(2,4),AIWWWW(4),HMASS,HGAM,HGAMS(29),  */
 162 /*     &       ETAHGG,MATCHH(29),ZSTARS(4,2),IHTYPE,HGAMSS(85,85)         */
 163 /*------------------------------------------------------------------------*/
 164 typedef struct {
 165     Float_t anwwww[4][4][4];
 166     Float_t adwwww[4][2];
 167     Float_t aiwwww[4];
 168     Float_t hmass;
 169     Float_t hgam;
 170     Float_t hgams[29];
 171     Float_t etahgg;
 172     Int_t matchh[29];
 173     Float_t zstars[2][4];
 174     Int_t ihtype;
 175     Float_t hgamss[85][85];
 176 } HconCommon;
 177
 178 #define HCON COMMON_BLOCK(HCON, hcon)
 179 COMMON_BLOCK_DEF(HconCommon,HCON);
 180
 181 /**************************************************************************/
 182 /*           D E S C R I P T I O N :                                      */
 183 /*------------------------------------------------------------------------*/
 184 /* anwwww :                                                               */
 185 /* adwwww :                                                               */
 186 /* aiwwww :                                                               */
 187 /* hmass :                                                                */
 188 /* hgam :                                                                 */
 189 /* hgams :                                                                */
 190 /* etahgg:                                                                */
 191 /* matchh :                                                               */
 192 /* zstars :                                                               */
 193 /* ihtype : MSSM Higgs type, either HL0, HH0 or HA0                       */
 194 /* hgamss :                                                               */
 195 /**************************************************************************/
 196
 197
 198
 199 /*========================================================================*/
 200 /* COMMON/JETLIM/PMIN(MXLIM),PMAX(MXLIM),PTMIN(MXLIM),PTMAX(MXLIM),       */
 201 /*     &         YJMIN(MXLIM),YJMAX(MXLIM),PHIMIN(MXLIM),PHIMAX(MXLIM),   */
 202 /*     &         XJMIN(MXLIM),XJMAX(MXLIM),THMIN(MXLIM),THMAX(MXLIM),     */
 203 /*     &         SETLMJ(12*MXLIM)                                         */
 204 /*------------------------------------------------------------------------*/
 205 typedef struct {
 206     static const Int_t mxlim = 8;
 207     Float_t pmin[mxlim];
 208     Float_t pmax[mxlim];
 209     Float_t ptmin[mxlim];
 210     Float_t ptmax[mxlim];
 211     Float_t yjmin[mxlim];
 212     Float_t yjmax[mxlim];
 213     Float_t phimin[mxlim];
 214     Float_t phimax[mxlim];
 215     Float_t xjmin[mxlim];
 216     Float_t xjmax[mxlim];
 217     Float_t thmin[mxlim];
 218     Float_t thmax[mxlim];
 219     Bool_t setlmj[12*mxlim];
 220 } JetlimCommon;
 221
 222 #define JETLIM COMMON_BLOCK(JETLIM, jetlim)
 223 COMMON_BLOCK_DEF(JetlimCommon,JETLIM);
 224
 225 /**************************************************************************/
 226 /*           D E S C R I P T I O N :                                      */
 227 /*------------------------------------------------------------------------*/
 228 /* Changes to mxlim should be echoed in MGLIMS and some of its access     */
 229 /* routines.                                                              */
 230 /* pmin, pmax : Momentum limits for jets.                                 */
 231 /* ptmin, ptmax : p-subscript-t limits for jets.                          */
 232 /* yjmin, yjmax : Y limits for each jet; not simultaneous to thmin, thmax.*/
 233 /* phimin, phimax : Phi limits for jets.                                  */
 234 /* xjmin, xjmax : Feynman x limits for jets.                              */
 235 /* thmin, thmax : Theta limits for jets, not simultaneous with yjmin,     */
 236 /*                yjmax.                                                  */
 237 /* setlmj :                                                               */
 238 /* yjset, thset : Logical flags internal to the interface.                 */
 239 /**************************************************************************/
 240
 241
 242
 243 /*========================================================================*/
 244 /* COMMON/JETPAR/P(3),PT(3),YJ(3),PHI(3),XJ(3),TH(3),CTH(3),STH(3),       */
 245 /*     &         JETTYP(3),SHAT,THAT,QSQ,X1,X2,PBEAM(2),                  */
 246 /*     &         QMW,QW,QTW,YW,XW,THW,QTMW,PHIW,SHAT1,THAT1,UHAT1,JWTYP,  */
 247 /*     &         INITYP(2),ISIGS,PBEAMS(5)                                */
 248 /*------------------------------------------------------------------------*/
 249 typedef struct {
 250     Float_t p[3];
 251     Float_t pt[3];
 252     Float_t yj[3];
 253     Float_t phi[3];
 254     Float_t xj[3];
 255     Float_t th[3];
 256     Float_t cth[3];
 257     Float_t sth[3];
 258     Int_t jettyp[3];
 259     Float_t shat;
 260     Float_t that;
 261     Float_t uhat;
 262     Float_t qsq;
 263     Float_t x1;
 264     Float_t x2;
 265     Float_t pbeam[2];
 266     Float_t qmw;
 267     Float_t qw;
 268     Float_t qtw;
 269     Float_t yw;
 270     Float_t xw;
 271     Float_t thw;
 272     Float_t qtmw;
 273     Float_t phiw;
 274     Float_t shat1;
 275     Float_t that1;
 276     Float_t uhat1;
 277     Int_t jwtyp;
 278     Float_t alfqsq;
 279     Float_t cthw;
 280     Float_t sthw;
 281     Float_t q0w;
 282     Int_t inityp[2];
 283     Int_t isigs;
 284     Float_t pbeams[5];
 285 } JetparCommon;
 286
 287 #define JETPAR COMMON_BLOCK(JETPAR, jetpar)
 288 COMMON_BLOCK_DEF(JetparCommon,JETPAR);
 289
 290
 291 /**************************************************************************/
 292 /*           D E S C R I P T I O N :                                      */
 293 /*------------------------------------------------------------------------*/
 294 /* jwtyp : W type, either W+-, Z0 or GM.                                  */
 295 /**************************************************************************/
 296
 297
 298
 299 /*========================================================================*/
 300 /* COMMON/KEYS/IKEY, KEYON, KEYS(MXKEYS)                                 */
 301 /*------------------------------------------------------------------------*/
 302
 303 typedef struct {
 304     static const Int_t mxkeys = 20;
 305     Bool_t keys[mxkeys];
 306     Bool_t keyon;
 307     Int_t ikey;
 308     Char_t *reac;
 309 } KeysCommon;
 310
 311 #define KEYS COMMON_BLOCK(KEYS, keys)
 312 COMMON_BLOCK_DEF(KeysCommon,KEYS);
 313
 314 /**************************************************************************/
 315 /*           D E S C R I P T I O N :                                      */
 316 /*------------------------------------------------------------------------*/
 317 /* keys : Flag for process type.                                          */
 318 /**************************************************************************/
 319
 320
 321
 322 /*========================================================================*/
 323 /* COMMON/KKGRAV/NEXTRAD,MASSD,KKGSD,SURFD,UVCUT                          */
 324 /*------------------------------------------------------------------------*/
 325
 326 typedef struct {
 327     Int_t nextrad;
 328     Float_t massd;
 329     Float_t kkgsd;
 330     Float_t surfd;
 331     Bool_t uvcut;
 332 } KkgravCommon;
 333
 334 #define KKGRAV COMMON_BLOCK(KKGRAV, kkgrav)
 335 COMMON_BLOCK_DEF(KkgravCommon,KKGRAV);
 336
 337 /**************************************************************************/
 338 /*           D E S C R I P T I O N :                                      */
 339 /*------------------------------------------------------------------------*/
 340 /* EXTRADIM process parameters :                                          */
 341 /* nextrad : delta                                                        */
 342 /* massd : M-subscript-D                                                  */
 343 /* uvcut : Logical flag                                                   */
 344 /*                                                                        */
 345 /* kkgsd :                                                                */
 346 /* surfd :                                                                */
 347 /**************************************************************************/
 348
 349
 350
 351 /*========================================================================*/
 352 /* COMMON/MBGEN/POMWT(LIMPOM),POMGEN(LIMPOM),MNPOM,MXPOM,PDIFFR,NPOM,     */
 353 /*     &        XBARY(2),DXBARY(2),XPOM(LIMPOM,2)                         */
 354 /*------------------------------------------------------------------------*/
 355 typedef struct {
 356     static const Int_t limpom = 20;
 357     Float_t pomwt[limpom];
 358     Float_t pomgen[limpom];
 359     Int_t mnpom;
 360     Int_t mxpom;
 361     Float_t pdiffr;
 362     Int_t npom;
 363     Float_t xbary[2];
 364     Float_t dxbary[2];
 365     Float_t xpom[2][limpom];
 366 } MbgenCommon;
 367
 368 #define MBGEN COMMON_BLOCK(MBGEN, mbgen)
 369 COMMON_BLOCK_DEF(MbgenCommon,MBGEN);
 370
 371 /**************************************************************************/
 372 /*           D E S C R I P T I O N :                                      */
 373 /*------------------------------------------------------------------------*/
 374 /* pomwt                                                                  */
 375 /* pomgen                                                                 */
 376 /* mnpom, mxpom : Min and max number of cut pomerons.                     */
 377 /* pdiffr                                                                 */
 378 /* npom                                                                   */
 379 /* xbary                                                                  */
 380 /* dxbary                                                                 */
 381 /* xpom                                                                   */
 382 /**************************************************************************/
 383
 384
 385
 386 /*========================================================================*/
 387 /* COMMON/MGLIMS/EHMGMN,EHMGMX,YHMGMN,YHMGMX,AMIJMN(MXLIM,MXLIM),         */
 388 /*     &         AMIJMX(MXLIM,MXLIM),FIXMIJ(MXLIM,MXLIM)                  */
 389 /*------------------------------------------------------------------------*/
 390 typedef struct {
 391     static const Int_t mxlim = 8;
 392     Float_t ehmgmn;
 393     Float_t ehmgmx;
 394     Float_t yhmgmn;
 395     Float_t yhmgmx;
 396     Float_t amijmn[mxlim][mxlim];
 397     Float_t amijmx[mxlim][mxlim];
 398     Bool_t fixmij[mxlim][mxlim];
 399 } MglimsCommon;
 400
 401 #define MGLIMS COMMON_BLOCK(MGLIMS, mglims)
 402 COMMON_BLOCK_DEF(MglimsCommon,MGLIMS);
 403
 404 /**************************************************************************/
 405 /*           D E S C R I P T I O N :                                      */
 406 /*------------------------------------------------------------------------*/
 407 /* Changes in mxlim should be echoed in JETLIM and in access routines for */
 408 /* amijmn and amijmx.                                                     */
 409 /* Limits for MadGraph multiparton processes                              */
 410 /* ehmgmn, ehmgmx : Mass range                                            */
 411 /* yhmgmn, yhmgmx :                                                       */
 412 /* amijmn, amijmx : Multimet mass limits                                  */
 413 /* fixmij :                                                               */
 414 /**************************************************************************/
 415
 416
 417
 418 /*========================================================================*/
 419 /* COMMON/NODCAY/NODCAY,NOETA,NOPI0,NONUNU,NOEVOL,NOHADR,NOGRAV           */
 420 /*------------------------------------------------------------------------*/
 421 typedef struct {
 422     Bool_t nodcay;
 423     Bool_t noeta;
 424     Bool_t nopi0;
 425     Bool_t nonunu;
 426     Bool_t noevol;
 427     Bool_t nohadr;
 428     Bool_t nograv;
 429 } NodcayCommon;
 430
 431 #define NODCAY COMMON_BLOCK(NODCAY, nodcay)
 432 COMMON_BLOCK_DEF(NodcayCommon,NODCAY);
 433
 434 /**************************************************************************/
 435 /*           D E S C R I P T I O N :                                      */
 436 /*------------------------------------------------------------------------*/
 437 /* Booleans used to forbid decay channels.                                */
 438 /* nodcay : Suppress all decays                                           */
 439 /* noeta  : Suppress eta decays                                           */
 440 /* nopi0  : Suppress pi-zero decays                                       */
 441 /* nonunu : Suppess Z-zero neutrino decays                                */
 442 /* noevol : Suppress QCD evolution and hadronisation                      */
 443 /* nohadr : Suppress hadronisation of jets and beam jets                  */
 444 /* nograv : Suppress gravitino decays in GMSB model                       */
 445 /**************************************************************************/
 446
 447
 448
 449 /*========================================================================*/
 450 /* COMMON/PARTCL/NPTCL,PPTCL(5, MXPTCL),IORIG(MXPTCL),IDENT(MXPTCL),      */
 451 /*     &         IDCAY(MXPTCL)                                            */
 452 /*------------------------------------------------------------------------*/
 453 typedef struct {
 454     static const Int_t mxptcl = 4000;
 455     static const Int_t ipack = 10000;
 456     Int_t nptcl;
 457     Float_t pptcl[mxptcl][5];
 458     Int_t iorig[mxptcl];
 459     Int_t ident[mxptcl];
 460     Int_t idcay[mxptcl];
 461 } PartclCommon;
 462
 463 #define PARTCL COMMON_BLOCK(PARTCL, partcl)
 464 COMMON_BLOCK_DEF(PartclCommon,PARTCL);
 465
 466 /**************************************************************************/
 467 /*           D E S C R I P T I O N :                                      */
 468 /*------------------------------------------------------------------------*/
 469 /* Stores HEPEVT data.                                                    */
 470 /* nptcl : Number of particles.                                           */
 471 /* pptcl : x, y, z, and 0 momentum of particles, and mass.                */
 472 /* iorig : origin of particles.                                           */
 473 /* ident : ID code.                                                       */
 474 /* idcay : ID of decay products.                                          */
 475 /* mxptcl : Max number of particles.                                      */
 476 /* ipack : Packing integer, used in decoding idcay.                       */
 477 /**************************************************************************/
 478
 479
 480
 481 /*========================================================================*/
 482 /* COMMON/PRIMAR/NJET,SCM,HALFE,ECM,IDIN(2),NEVENT,NTRIES,NSIGMA          */
 483 /*------------------------------------------------------------------------*/
 484 typedef struct {
 485     Int_t njet;
 486     Float_t scm;
 487     Float_t halfe;
 488     Float_t ecm;
 489     Int_t idin[2];
 490     Int_t nevent;
 491     Int_t ntries;
 492     Int_t nsigma;
 493 } PrimarCommon;
 494
 495 #define PRIMAR COMMON_BLOCK(PRIMAR, primar)
 496 COMMON_BLOCK_DEF(PrimarCommon,PRIMAR);
 497
 498 /**************************************************************************/
 499 /*           D E S C R I P T I O N :                                      */
 500 /*------------------------------------------------------------------------*/
 501 /* njet :                                                                 */
 502 /* scm : Square of the center-of-mass energy.                             */
 503 /* halfe : Half the center-of-mass energy.                                */
 504 /* ecm : Center-of-mass energy.                                           */
 505 /* idin : Beam types 1 and 2.                                             */
 506 /* nevent : Number of events to build.                                    */
 507 /* ntries : Maximum number of tries to find a good event.                 */
 508 /* nsigma : Number of unevolved events for SIGF calculation.              */
 509 /**************************************************************************/
 510
 511
 512
 513 /*========================================================================*/
 514 /* COMMON/PRTOUT/NEVPRT,NJUMP                                             */
 515 /*------------------------------------------------------------------------*/
 516 typedef struct {
 517     Int_t nevprt;
 518     Int_t njump;
 519 } PrtoutCommon;
 520
 521 #define PRTOUT COMMON_BLOCK(PRTOUT, prtout)
 522 COMMON_BLOCK_DEF(PrtoutCommon,PRTOUT);
 523
 524 /**************************************************************************/
 525 /*           D E S C R I P T I O N :                                      */
 526 /*------------------------------------------------------------------------*/
 527 /* nevprt : Number of events to print                                     */
 528 /* njump : Print every njump events.                                      */
 529 /**************************************************************************/
 530
 531
 532
 533 /*========================================================================*/
 534 /* COMMON/QCDPAR/ALAM,ALAM2,CUTJET,ISTRUC                                 */
 535 /*------------------------------------------------------------------------*/
 536 typedef struct {
 537     Float_t alam;
 538     Float_t alam2;
 539     Float_t cutjet;
 540     Int_t istruc;
 541 } QcdparCommon;
 542
 543 #define QCDPAR COMMON_BLOCK(QCDPAR, qcdpar)
 544 COMMON_BLOCK_DEF(QcdparCommon,QCDPAR);
 545
 546 /**************************************************************************/
 547 /*           D E S C R I P T I O N :                                      */
 548 /*------------------------------------------------------------------------*/
 549 /* alam : QCD scale (LAMBDA)                                              */
 550 /* alam2 : Square of alam                                                 */
 551 /* cutjet : Cutoff mass for QCD jet evolution (mu-subscript-c)            */
 552 /* istruc : Structure functions CTEQ3L, CTEQ2L, EHLQ or DO.               */
 553 /**************************************************************************/
 554
 555
 556
 557 /*========================================================================*/
 558 /* COMMON/QLMASS/AMLEP(100),NQLEP,NMES,NBARY                              */
 559 /*------------------------------------------------------------------------*/
 560 typedef struct {
 561     static const Int_t nqlep = 61;
 562     static const Int_t nmes = 2;
 563     static const Int_t nbary = 2;
 564     Float_t amlep[100];
 565 } QlmassCommon;
 566
 567 #define QLMASS COMMON_BLOCK(QLMASS, qlmass)
 568 COMMON_BLOCK_DEF(QlmassCommon,QLMASS);
 569
 570 /**************************************************************************/
 571 /*           D E S C R I P T I O N :                                      */
 572 /*------------------------------------------------------------------------*/
 573 /* amlep : (C++ indices, Fortran is one greater)                          */
 574 /* [5-7] : t, y and x quark mass (m-subscript-t, -y, -x)                  */
 575 /* [21-26] : Masses of u, d, s, c, b and t, all -tilde                    */
 576 /* [29] : g-tilde mass                                                    */
 577 /* [30] : gamma-tilde mass                                                */
 578 /* [31-36] : Masses for nu-e, e, nu-mu, mu, nu-tau and tau, all -tilde    */
 579 /* [39] : W-plus-tilde mass                                               */
 580 /* [40] : Z-zero-tilde mass                                               */
 581 /* [63-71] : Higgs meson masses, charges 0,0,0,0,0,1,1,2,2                */
 582 /**************************************************************************/
 583
 584
 585
 586 /*========================================================================*/
 587 /* COMMON/SEED/XSEED                                                      */
 588 /*------------------------------------------------------------------------*/
 589 typedef struct {
 590     Char_t xseed[24];
 591 } SeedCommon;
 592
 593 #define SEED COMMON_BLOCK(SEED, seed)
 594 COMMON_BLOCK_DEF(SeedCommon,SEED);
 595
 596 /**************************************************************************/
 597 /*           D E S C R I P T I O N :                                      */
 598 /*------------------------------------------------------------------------*/
 599 /* Random number seed                                                     */
 600 /**************************************************************************/
 601
 602
 603
 604 /*========================================================================*/
 605 /* COMMON/SUGNU/XNUSUG(18)                                                */
 606 /*------------------------------------------------------------------------*/
 607 typedef struct {
 608     Float_t xnusug[18];
 609 } SugnuCommon;
 610
 611 #define SUGNU COMMON_BLOCK(SUGNU, sugnu)
 612 COMMON_BLOCK_DEF(SugnuCommon,SUGNU);
 613
 614 /**************************************************************************/
 615 /*           D E S C R I P T I O N :                                      */
 616 /*------------------------------------------------------------------------*/
 617 /* Non-universal SUGRA terms : (C++ indexes again)                        */
 618 /* xnusug[0-2] : Gaugino masses                                           */
 619 /* [3-5] : A terms (A-subscript-tau, -b, -t)                              */
 620 /* [6-7] : Higgs masses (M-subscript-H-subscript-d, -u)                   */
 621 /* [8-12] : 1st / 2nd generation masses (M-subscript-eR, -eL, -dR,        */
 622 /*                                                  -uR, -uL)             */
 623 /* [13-17] : 3rd generation masses (M-subscript-tau R, -tau L, -bR,       */
 624 /*                                             -tR, -tL)                  */
 625 /**************************************************************************/
 626
 627
 628
 629 /*========================================================================*/
 630 /* COMMON/TCPAR/TCMRHO,TCGRHO                                             */
 631 /*------------------------------------------------------------------------*/
 632 typedef struct {
 633     Float_t tcmrho;
 634     Float_t tcgrho;
 635 } TcparCommon;
 636
 637 #define TCPAR COMMON_BLOCK(TCPAR, tcpar)
 638 COMMON_BLOCK_DEF(TcparCommon,TCPAR);
 639
 640 /**************************************************************************/
 641 /*           D E S C R I P T I O N :                                      */
 642 /*------------------------------------------------------------------------*/
 643 /* tcmrho : Technicolor mass                                              */
 644 /* tcgrho : Technicolor width                                             */
 645 /**************************************************************************/
 646
 647
 648
 649 /*========================================================================*/
 650 /* COMMON/TYPES/LOC(100),NTYP,NJTTYP,NWWTYP(2),NWMODE(3)                  */
 651 /*------------------------------------------------------------------------*/
 652 typedef struct {
 653     static const Int_t mxtype = 8;
 654     Int_t loc[100];
 655     Int_t ntyp;
 656     Int_t njttyp[mxtype];
 657     Int_t nwwtyp[2];
 658     Int_t nwmode[3];
 659 } TypesCommon;
 660
 661 #define TYPES COMMON_BLOCK(TYPES, types)
 662 COMMON_BLOCK_DEF(TypesCommon,TYPES);
 663
 664 /**************************************************************************/
 665 /*           D E S C R I P T I O N :                                      */
 666 /*------------------------------------------------------------------------*/
 667 /* Note : Any change in mxtype should be copied in xtypes.                */
 668 /* loc :                                                                  */
 669 /* ntyp :                                                                 */
 670 /* njttyp : Number of jet types in jet [index].                           */
 671 /* nwwtyp :                                                               */
 672 /* nwmode :                                                               */
 673 /**************************************************************************/
 674
 675
 676
 677 /*========================================================================*/
 678 /* COMMON/XMSSM/GOMSSM,GOSUG,GOGMSB,GOAMSB,AL3UNI,                        */
 679 /*     &        XGLSS,XMUSS,XHASS,XTBSS,                                  */
 680 /*     &        XQ1SS,XDRSS,XURSS,XL1SS,XERSS,                            */
 681 /*     &        XQ2SS,XSRSS,XCRSS,XL2SS,XMRSS,                            */
 682 /*     &        XQ3SS,XBRSS,XTRSS,XL3SS,XTARSS,XATSS,XABSS,XATASS,        */
 683 /*     &        XM1SS,XM2SS,XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU,              */
 684 /*     &        XLAMGM,XMESGM,XN5GM,XCMGV,XMGVTO,                         */
 685 /*     &        XRSLGM,XDHDGM,XDHUGM,XDYGM,XN51GM,XN52GM,XN53GM,          */
 686 /*     &        XMN3NR,XMAJNR,XANSS,XNRSS,XSBCS                           */
 687 /*------------------------------------------------------------------------*/
 688
 689 typedef struct {
 690     Bool_t gomssm;
 691     Bool_t gosug;
 692     Bool_t gogmsb;
 693     Bool_t goamsb;
 694     Bool_t al3uni;
 695     Float_t xglss;
 696     Float_t xmuss;
 697     Float_t xhass;
 698     Float_t xtbss;
 699     Float_t xq1ss;
 700     Float_t xdrss;
 701     Float_t xurss;
 702     Float_t xl1ss;
 703     Float_t xerss;
 704     Float_t xq2ss;
 705     Float_t xsrss;
 706     Float_t xcrss;
 707     Float_t xl2ss;
 708     Float_t xmrss;
 709     Float_t xq3ss;
 710     Float_t xbrss;
 711     Float_t xtrss;
 712     Float_t xl3ss;
 713     Float_t xtarss;
 714     Float_t xatss;
 715     Float_t xabss;
 716     Float_t xatass;
 717     Float_t xm1ss;
 718     Float_t xm2ss;
 719     Float_t xm0su;
 720     Float_t xmhsu;
 721     Float_t xa0su;
 722     Float_t xtgbsu;
 723     Float_t xsmusu;
 724     Float_t xlamgm;
 725     Float_t xmesgm;
 726     Float_t xn5gm;
 727     Float_t xcmgv;
 728     Float_t mgvto;
 729     Float_t xrslgm;
 730     Float_t xdhdgm;
 731     Float_t xdhugm;
 732     Float_t xdygm;
 733     Float_t xn51gm;
 734     Float_t xn52gm;
 735     Float_t xn53gm;
 736     Float_t xmn3nr;
 737     Float_t xmajnr;
 738     Float_t xanss;
 739     Float_t xnrss;
 740     Float_t xsbcs;
 741
 742 } XmssmCommon;
 743
 744 #define XMSSM COMMON_BLOCK(XMSSM, xmssm)
 745 COMMON_BLOCK_DEF(XmssmCommon,XMSSM);
 746
 747 /**************************************************************************/
 748 /*           D E S C R I P T I O N :                                      */
 749 /*------------------------------------------------------------------------*/
 750 /* Keyword flags :                                                        */
 751 /* gomssm : True if any of the MSSM* keywords or SUGRA, GMSB or AMSB have */
 752 /*          been used.                                                    */
 753 /* gosug  : True if SUGRA has been used.                                  */
 754 /* gogmsb : True if GMSB has been used.                                   */
 755 /* goamsb : True if AMSB has been used.                                   */
 756 /* al3uni : True if AL3UNI has been used.                                 */
 757 /*                                                                        */
 758 /* MSSM parameters :                                                      */
 759 /* xglss : Gluino mass                                                    */
 760 /* xmuss : mu                                                             */
 761 /* xhass : A mass                                                         */
 762 /* xtbss : tan beta                                                       */
 763 /*                                                                        */
 764 /* MSSM first generation :                                                */
 765 /* xq1ss : q-subscript-1 mass                                             */
 766 /* xdrss : d-subscript-r mass                                             */
 767 /* xurss : u-subscript-r mass                                             */
 768 /* xl1ss : l-subscript-1 mass                                             */
 769 /* xerss : e-subscript-1 mass                                             */
 770 /*                                                                        */
 771 /* MSSM second generation :                                               */
 772 /* xq2ss : q-subscript-2 mass                                             */
 773 /* xsrss : s-subscript-r mass                                             */
 774 /* xcrss : c-subscript-r mass                                             */
 775 /* xl2ss : l-subscript-2 mass                                             */
 776 /* xmrss : mu-subscript-r mass                                            */
 777 /*                                                                        */
 778 /* MSSM third generation :                                                */
 779 /* xq3ss : q-subscript-3 mass                                             */
 780 /* xbrss : b-subscript-r mass                                             */
 781 /* xtrss : t-subscript-r mass                                             */
 782 /* xl3ss : l-subscript-3 mass                                             */
 783 /* xtarss : tau-subscript-r mass                                          */
 784 /* xatss, xabss, xatass : A-subscript-t, -b, -tau mixings.                */
 785 /*                                                                        */
 786 /* MSSM gaugino masses                                                    */
 787 /* xm1ss, xm2ss : M-subscript-1, -2                                       */
 788 /*                                                                        */
 789 /* Anomaly-mediated SUSY breaking / Minimal supergravity parameters :     */
 790 /* xm0su : scalar mass (m-subscript-0)                                    */
 791 /* xmhsu : gravitino mass (m-subscript-three halves)                      */
 792 /* xa0su : trilinear breaking term (A-subscript-0)                        */
 793 /* xtgbsu : VEV ratio (tan beta)                                          */
 794 /* xsmusu : sign (sgn mu)                                                 */
 795 /*                                                                        */
 796 /* GMSB messinger SUSY breaking                                           */
 797 /* xlamgm : mass (LAMBDA-subscript-m)                                     */
 798 /* xmesgm :                                                               */
 799 /* xn5gm : number of 5 + 5-bar (M-subscript-m)                            */
 800 /* xcmgv : gravitino scale (C-subscript-gr)                               */
 801 /*                                                                        */
 802 /* Non-minimal GMSB parameters                                            */
 803 /* xrslgm : gaugino mass multiplier (R-slash)                             */
 804 /* xdhdgm : Higgs mass shift (delta M-squared-subscript-Hd)               */
 805 /* xdhugm : Higgs mass shift (delta M-squared-subscript-Hu)               */
 806 /* xdygm : D-term mass-squared (D-subscript-Y of M)                       */
 807 /* xn51gm, xn52gm, xn53gm : independent gauge group messengers            */
 808 /*                          (N-subscript-5-subscript-1, -2, -3)           */
 809 /* SUGRA see-saw nu-effect                                                */
 810 /* xmn3nr : nu-mass (m-subscript-nu-subscript-tau)                        */
 811 /* xmajnr : int. scale (M-subscript-N)                                    */
 812 /* xanss : GUT scale (A-subscript-n)                                      */
 813 /* xnrss : nu SSB terms (m-subscript-nu-tilde-subscript-R)                */
 814 /*                                                                        */
 815 /* xsbcs                                                                  */
 816 /*                                                                        */
 817 /* Gravitino mass                                                         */
 818 /* xmgvto : gravitino mass (M-subscript-gravitino)                        */
 819 /**************************************************************************/
 820
 821
 822
 823 /*========================================================================*/
 824 /* COMMON/XTYPES/PARTYP(40),TITLE(10),JETYP(30,MXTYPE),WWTYP(30,2),       */
 825 /*     &         WMODES(30,3)                                             */
 826 /*------------------------------------------------------------------------*/
 827 typedef struct {
 828     static const Int_t mxtype = 8;
 829     Char_t* partyp[40];
 830     Char_t* title;
 831     Char_t* jetyp[mxtype][30];
 832     Char_t* wwtyp[2][30];
 833     Char_t* wmodes[3][30];
 834 } XtypesCommon;
 835
 836 #define XTYPES COMMON_BLOCK(XTYPES, xtypes)
 837 COMMON_BLOCK_DEF(XtypesCommon,XTYPES);
 838
 839 /**************************************************************************/
 840 /*           D E S C R I P T I O N :                                      */
 841 /*------------------------------------------------------------------------*/
 842 /* Note : Any change in mxtype should be copied in types.                 */
 843 /* partyp                                                                 */
 844 /* title                                                                  */
 845 /* jetyp                                                                  */
 846 /* wwtyp : Decay modes for WPAIR                                          */
 847 /* wmodes : Decay modes for weak force carriers                           */
 848 /**************************************************************************/
 849
 850
 851
 852 /*========================================================================*/
 853 /* COMMON/WCON/SIN2W,WMASS(4),WGAM(4),AQ(12,4),BQ(12,4),COUT(4),          */
 854 /*     &       MATCH(25,4),WCBR(25,4),CUTOFF,CUTPOW,TBRWW(4,2),           */
 855 /*     &       RBRWW(12,4,2),EZ,AQDP(12,4),BQDP(12,4),EZDP,WFUDGE         */
 856 /*------------------------------------------------------------------------*/
 857 typedef struct {
 858
 859     Float_t sin2w;
 860     Float_t wmass[4];
 861     Float_t wgam[4];
 862     Float_t aq[4][12];
 863     Float_t bq[4][12];
 864     Float_t cout[4];
 865     Int_t match;
 866     Float_t wcbr[4][25];
 867     Float_t cutoff;
 868     Float_t cutpow;
 869     Float_t tbrww[2][4];
 870     Float_t rbrww[2][4][12];
 871     Float_t ez;
 872     Float_t aqdp[4][12];
 873     Float_t bqdp[4][12];
 874     Float_t ezdp;
 875     Float_t wfudge;
 876 } WconCommon;
 877
 878 #define WCON COMMON_BLOCK(WCON, wcon)
 879 COMMON_BLOCK_DEF(WconCommon, WCON);
 880
 881 /**************************************************************************/
 882 /*           D E S C R I P T I O N :                                      */
 883 /*------------------------------------------------------------------------*/
 884 /* sin2w : Weinberg angle (sin-squared(theta-subscript-W))                */
 885 /* wmass : W and Z masses                                                 */
 886 /* wgam  :                                                                */
 887 /* aq    :                                                                */
 888 /* bq    :                                                                */
 889 /* cout  :                                                                */
 890 /* match                                                                  */
 891 /* wcbr                                                                   */
 892 /* cutoff, cutpow : mu-square and nu, respectively, in formula            */
 893 /*                  cutoff q*t-square = mu-square*Q-power-nu              */
 894 /*                  for DRELLYAN events.                                  */
 895 /* tbrww                                                                  */
 896 /* rbrww                                                                  */
 897 /* ez                                                                     */
 898 /* aqdp                                                                   */
 899 /* bqdp                                                                   */
 900 /* ezdp                                                                   */
 901 /* wfudge : Fudge factor for DRELLYAN evolution scale                     */
 902 /**************************************************************************/
 903 }
 904
 905 #endif
 906
 907 // Endfile.
 908
 909
 910
 911
 912
 913
 914
 915