git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Obsolete - removed
[u/mrichter/AliRoot.git] / TIsajet / Icommon.h
CommitLineData
7dae8538 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
10extern "C" {
11
12/*========================================================================*/
13/* COMMON/DYLIM/QMIN,QMAX,QTMIN,QTMAX,YWMIN,YWMAX,XWMIN,XWMAX,THWMIN, */
14/* & THWMAX,PHWMIN,PHWMAX,SETLMQ(12) */
15/*------------------------------------------------------------------------*/
16typedef 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)
35COMMON_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/*------------------------------------------------------------------------*/
57typedef 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)
70COMMON_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/*------------------------------------------------------------------------*/
94typedef 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)
105COMMON_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/*------------------------------------------------------------------------*/
125typedef 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)
141COMMON_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/*------------------------------------------------------------------------*/
164typedef 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)
179COMMON_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/*------------------------------------------------------------------------*/
205typedef 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)
223COMMON_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/*------------------------------------------------------------------------*/
249typedef 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)
288COMMON_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
303typedef 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)
312COMMON_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
326typedef 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)
335COMMON_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/*------------------------------------------------------------------------*/
355typedef 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)
369COMMON_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/*------------------------------------------------------------------------*/
390typedef 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)
402COMMON_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/*------------------------------------------------------------------------*/
421typedef 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)
432COMMON_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/*------------------------------------------------------------------------*/
453typedef 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)
464COMMON_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/*------------------------------------------------------------------------*/
484typedef 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)
496COMMON_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/*------------------------------------------------------------------------*/
516typedef struct {
517 Int_t nevprt;
518 Int_t njump;
519} PrtoutCommon;
520
521#define PRTOUT COMMON_BLOCK(PRTOUT, prtout)
522COMMON_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/*------------------------------------------------------------------------*/
536typedef 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)
544COMMON_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/*------------------------------------------------------------------------*/
560typedef 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)
568COMMON_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/*------------------------------------------------------------------------*/
589typedef struct {
590 Char_t xseed[24];
591} SeedCommon;
592
593#define SEED COMMON_BLOCK(SEED, seed)
594COMMON_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/*------------------------------------------------------------------------*/
607typedef struct {
608 Float_t xnusug[18];
609} SugnuCommon;
610
611#define SUGNU COMMON_BLOCK(SUGNU, sugnu)
612COMMON_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/*------------------------------------------------------------------------*/
632typedef struct {
633 Float_t tcmrho;
634 Float_t tcgrho;
635} TcparCommon;
636
637#define TCPAR COMMON_BLOCK(TCPAR, tcpar)
638COMMON_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/*------------------------------------------------------------------------*/
652typedef 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)
662COMMON_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
689typedef 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)
745COMMON_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/*------------------------------------------------------------------------*/
827typedef 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)
837COMMON_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/*------------------------------------------------------------------------*/
857typedef 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)
879COMMON_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