Hadr02/include/G4UrQMD1_3Interface.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // *                                                                  *
0021 // * Parts of this code which have been  developed by Abdel-Waged     *
0022 // * et al under contract (31-465) to the King Abdul-Aziz City for    *
0023 // * Science and Technology (KACST), the National Centre of           *
0024 // * Mathematics and Physics (NCMP), Saudi Arabia.                    *
0025 // *                                                                  *
0026 // * By using,  copying,  modifying or  distributing the software (or *
0027 // * any work based  on the software)  you  agree  to acknowledge its *
0028 // * use  in  resulting  scientific  publications,  and indicate your *
0029 // * acceptance of all terms of the Geant4 Software license.          *
0030 // ********************************************************************
0031 //
0032 /// \file hadronic/Hadr02/include/G4UrQMD1_3Interface.hh
0033 /// \brief Definition of the G4UrQMD1_3Interface class
0034 //
0035 //
0036
0037 #ifndef G4UrQMD1_3Interface_hh
0038 #define G4UrQMD1_3Interface_hh
0039
0040 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0041 //
0042 // MODULE:           G4UrQMD1_3Model.hh
0043 //
0044 // Version:          0.B
0045 // Date:             20/12/12
0046 // Author:           Kh. Abdel-Waged and Nuha Felemban
0047 // Revised by:       V.V. Uzhinskii
0048 //                   SPONSERED BY
0049 // Customer:         KAUST/NCMP
0050 // Contract:         31-465
0051 //
0052 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0053 //
0054 //
0055 // Class Description
0056 //
0057 //
0058 // Class Description - End
0059 //
0060 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0061 ///////////////////////////////////////////////////////////////////////////////
0062
0063 #include "G4SystemOfUnits.hh"
0064 #include "globals.hh"
0065
0066 //  coms
0067 //
0068 const G4int nmax = 500;
0069 const G4int nspl = 500;
0070 const G4int smax = 500;
0071 //  comres
0072 const G4int minnuc = 1;
0073 const G4int minmes = 100;
0074 const G4int maxmes = 132;
0075 const G4int numnuc = 16;
0076 const G4int numdel = 10;
0077 const G4int maxnuc = minnuc + numnuc - 1;
0078 const G4int mindel = minnuc + maxnuc;
0079 const G4int maxdel = mindel + numdel - 1;
0080 const G4int minres = minnuc + 1;
0081 const G4int maxres = maxdel;
0082 const G4int numlam = 13;
0083 const G4int numsig = 9;
0084 const G4int numcas = 6;
0085 const G4int numome = 1;
0086 const G4int minlam = mindel + numdel;
0087 const G4int maxlam = minlam + numlam - 1;
0088 const G4int minsig = minlam + numlam;
0089 const G4int maxsig = minsig + numsig - 1;
0090 const G4int mincas = minsig + numsig;
0091 const G4int maxcas = mincas + numcas - 1;
0092 const G4int minome = mincas + numcas;
0093 const G4int maxome = minome + numome - 1;
0094 const G4int minbar = minnuc;
0095 const G4int maxbar = maxome;
0096 const G4int offmeson = minmes;
0097 const G4int maxmeson = maxmes;
0098 const G4int maxbra = 11;
0099 const G4int maxbrm = 25;
0100 const G4int maxbrs1 = 10;
0101 const G4int maxbrs2 = 3;
0102 const G4int nsigs = 10;
0103 const G4int itblsz = 100;
0104 const G4int maxreac = 13;
0105 const G4int maxpsig = 12;
0106 //
0107 // comwid
0108 //
0109 const G4int widnsp = 120;
0110 const G4double mintab = 0.10;
0111 const G4double maxtab1 = 5.0;
0112 const G4double maxtab2 = 50.0;
0113 const G4int tabver = 9;
0114 //
0115 // options
0116 //
0117 const G4int numcto = 400;
0118 const G4int numctp = 400;
0119 const G4int maxstables = 20;
0120 //
0121 // colltab (collision tables)
0122 //
0123 const G4int ncollmax = 100;
0124 //
0125 // inputs
0126 //
0127 const G4int aamax = 300;
0128 //
0129 // newpart (new created particles)
0130 //
0131 const G4int mprt = 200;
0132 const G4int oprt = 2;
0133 //
0134 // boxinc
0135 //
0136 const G4int bptmax = 20;
0137 //
0138
0139 // This next line is required as the default version of FORTRAN LOGICAL is
0140 // four bytes long, whereas storage for G4bool is one byte.
0141 //
0142 // comnorm
0143 const G4int n = 400;
0144 //
0145 // comstr
0146 const G4int njspin = 8;
0147 //
0148 // iso
0149 const G4int jmax = 7;
0150
0151 // This next line is required as the default version of FORTRAN LOGICAL is
0152 // four bytes long, whereas storage for G4bool is one byte.
0153 //
0154
0155 typedef G4int ftnlogical;
0156
0157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0158 //
0159 // Standard common block for UrQMD
0160 // Common options for coms.f
0161 //  20 commons
0162 //
0163 //
0164 struct ccurqmd13urqmdparams
0165 {
0166     G4int u_at, u_zt, u_ap, u_zp;
0167     G4double u_elab, u_imp;
0168     G4int u_sptar, u_spproj;
0169 };
0170
0171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0172
0173 struct ccurqmd13sys
0174 {
0175     G4int npart, nbar, nmes, ctag, nsteps, uid_cnt, ranseed, event, ap, at, zp, zt, eos, dectag,
0176       nhardres, nsoftres, ndecres, nelcoll, nblcoll;
0177 };
0178
0179 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0180
0181 struct ccurqmd13rsys
0182 {
0183     G4double time, acttime, bdist, bimp, bmin, ebeam, ecm;
0184 };
0185
0186 struct ccurqmd13comseed
0187 {
0188     ftnlogical firstseed;
0189 };
0190
0191 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0192
0193 struct ccurqmd13logic
0194 {
0195     ftnlogical lsct[nmax], logSky, logYuk, logCb, logPau;
0196 };
0197
0198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0199
0200 struct ccurqmd13mdprop
0201 {
0202     G4double r0_t[nmax], rx_t[nmax], ry_t[nmax], rz_t[nmax];
0203 };
0204
0205 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0206
0207 struct ccurqmd13cuts
0208 {
0209     G4double cutmax, cutPau, cutCb, cutYuk, cutSky, cutdww;
0210 };
0211
0212 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0213
0214 struct ccurqmd13spdata
0215 {
0216     G4double spx[nspl], spPauy[nspl], outPau[nspl], spCby[nspl], outCb[nspl], spYuky[nspl],
0217       outYuk[nspl], spSkyy[nspl], outSky[nspl], spdwwy[nspl], outdww[nspl];
0218 };
0219
0220 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0221
0222 struct ccurqmd13isys
0223 {
0224     G4int spin[nmax], ncoll[nmax], charge[nmax], ityp[nmax], lstcoll[nmax], iso3[nmax],
0225       origin[nmax], strid[nmax], uid[nmax];
0226 };
0227
0228 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0229
0230 struct ccurqmd13coor
0231 {
0232     G4double r0[nmax], rx[nmax], ry[nmax], rz[nmax], p0[nmax], px[nmax], py[nmax], pz[nmax],
0233       fmass[nmax], rww[nmax], dectime[nmax];
0234 };
0235
0236 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0237
0238 struct ccurqmd13frag
0239 {
0240     G4double tform[nmax], xtotfac[nmax];
0241 };
0242
0243 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0244
0245 struct ccurqmd13aios
0246 {
0247     G4double airx[nmax], airy[nmax], airz[nmax], aipx[nmax], aipy[nmax], aipz[nmax], aorx[4][nmax],
0248       aory[4][nmax], aorz[4][nmax], aopx[4][nmax], aopy[4][nmax], aopz[4][nmax];
0249 };
0250
0251 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0252
0253 struct ccurqmd13pots
0254 {
0255     G4double Cb0, Yuk0, Pau0, Sky20, Sky30, gamSky, gamYuk, drPau, dpPau, gw, sgw, delr, fdel, dt,
0256       da, db, dtimestep;
0257 };
0258
0259 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0260
0261 struct ccurqmd13scoor
0262 {
0263     G4double r0s[smax], rxs[smax], rys[smax], rzs[smax], p0s[smax], pxs[smax], pys[smax], pzs[smax],
0264       sfmass[smax];
0265 };
0266
0267 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0268
0269 struct ccurqmd13sisys
0270 {
0271     G4int sspin[smax], scharge[smax], sityp[smax], siso3[smax], suid[smax];
0272 };
0273
0274 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0275
0276 struct ccurqmd13ssys
0277 {
0278     G4int nspec;
0279 };
0280
0281 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0282
0283 struct ccurqmd13rtdelay
0284 {
0285     G4double p0td[nmax][2], pxtd[nmax][2], pytd[nmax][2], pztd[nmax][2], fmasstd[nmax][2];
0286 };
0287
0288 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0289
0290 struct ccurqmd13itdelay
0291 {
0292     G4int ityptd[nmax][2], iso3td[nmax][2];
0293 };
0294
0295 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0296
0297 struct ccurqmd13svinfo
0298 {
0299     G4int itypt[2], uidt[2], origint[2], iso3t[2];
0300 };
0301
0302 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0303
0304 struct ccurqmd13ffermi
0305 {
0306     G4double ffermpx[nmax], ffermpy[nmax], ffermpz[nmax];
0307 };
0308
0309 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0310
0311 struct ccurqmd13peq
0312 {
0313     G4double peq1, peq2;
0314 };
0315
0316 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0317 //
0318 // Definition for Collision Term
0319 // Commons  comres
0320 // 4 commons
0321 //
0322
0323 struct ccurqmd13versioning
0324 {
0325     char versiontag[45];
0326 };
0327
0328 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0329
0330 struct ccurqmd13resonances
0331 {
0332     G4double massres[maxbar - minbar + 1], widres[maxbar - minbar + 1];
0333     G4double massmes[maxmes - minmes + 1];
0334     G4double widmes[maxmes - minmes + 1];
0335     G4double mmesmn[maxmes - minmes + 1];
0336     G4double branres[maxdel - minnuc][maxbra + 1];
0337     G4double branmes[maxmes - minmes][maxbrm + 1];
0338
0339     G4double branbs1[maxsig - minlam][maxbrs1 + 1];
0340     G4double branbs2[maxcas - mincas][maxbrs2 + 1];
0341
0342     G4int bs1type[maxbrs1 + 1][4], bs2type[maxbrs2 + 1][4];
0343     G4int lbs1[maxsig - minlam][maxbrs1 + 1];
0344     G4int lbs2[maxcas - mincas][maxbrs2 + 1];
0345     G4int lbm[maxmes - minmes][maxbrm + 1];
0346
0347     G4int jres[maxbar - minbar + 1];
0348     G4int jmes[maxmes - minmes + 1];
0349     G4int lbr[maxdel - minnuc][maxbra + 1];
0350     G4int brtype[maxbra + 1][4];
0351     G4int pares[maxbar - minbar + 1], pames[maxmes - minmes + 1];
0352     G4int bmtype[maxbrm + 1][4];
0353     G4int isores[maxbar - minbar + 1], isomes[maxmes - minmes + 1];
0354     G4int strres[maxbar - minbar + 1], strmes[maxmes - minmes + 1];
0355     G4int mlt2it[maxmes - minmes];
0356 };
0357
0358 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0359
0360 struct ccurqmd13sigtabi
0361 {
0362     G4int sigmaln[maxreac][2][maxpsig];
0363     G4int sigmainf[20][nsigs];
0364 };
0365
0366 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0367
0368 struct ccurqmd13sigtabr
0369 {
0370     G4double sigmas[itblsz][nsigs], sigmascal[5][nsigs];
0371 };
0372
0373 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0374
0375 // comwid
0376 struct ccurqmd13decaywidth
0377 {
0378     G4double tabx[widnsp];
0379     G4double fbtaby[2][maxbar - minbar + 1][widnsp];
0380     G4double pbtaby[maxbra + 1][maxbar - minbar + 1][2][widnsp];
0381     G4double fmtaby[2][maxmes - minmes + 1][widnsp];
0382     G4double pmtaby[maxbrm + 1][maxmes - minmes + 1][2][widnsp];
0383     G4int wtabflg;
0384 };
0385
0386 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0387
0388 struct ccurqmd13brwignorm
0389 {
0390     G4double bwbarnorm[maxbar - minbar + 1];
0391     G4double bwmesnorm[maxmes - minmes + 1];
0392 };
0393
0394 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0395
0396 struct ccurqmd13xsections
0397 {
0398     G4double tabxnd[widnsp];
0399     G4double frrtaby[maxdel - 1][2][2][widnsp];
0400 };
0401
0402 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0403
0404 struct ccurqmd13tabnames
0405 {
0406     char tabname[77];
0407 };
0408
0409 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0410 //
0411 // options
0412 //
0413 struct ccurqmd13options
0414 {
0415     G4int CTOption[numcto];
0416     G4double CTParam[numctp];
0417 };
0418
0419 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0420
0421 struct ccurqmd13optstrings
0422 {
0423     char ctodc[numcto][2];
0424     char ctpdc[numctp][2];
0425 };
0426
0427 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0428
0429 struct ccurqmd13loptions
0430 {
0431     ftnlogical fixedseed, bf13, bf14, bf15, bf16, bf17, bf18, bf19, bf20;
0432 };
0433
0434 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0435
0436 struct ccurqmd13stables
0437 {
0438     G4int nstable;
0439     G4int stabvec[maxstables];
0440 };
0441
0442 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0443 //
0444 // colltab
0445 //
0446 struct ccurqmd13colltab
0447 {
0448     G4double cttime[ncollmax + 1], ctsqrts[ncollmax], ctsigtot[ncollmax], tmin;
0449     G4int cti1[ncollmax], cti2[ncollmax];
0450     G4int nct, actcol;
0451     ftnlogical ctvalid[ncollmax];
0452     G4int ctsav[ncollmax];
0453     G4int nsav, apt;
0454     G4double ctcolfluc[ncollmax];
0455 };
0456
0457 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0458 //
0459 // inputs
0460 //
0461 struct ccurqmd13inputs
0462 {
0463     G4int nevents, spityp[2], prspflg;
0464     G4int trspflg, spiso3[2], outsteps, bflag, srtflag, efuncflag;
0465     G4int nsrt, firstev, npb;
0466 };
0467
0468 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0469
0470 struct ccurqmd13input2
0471 {
0472     G4double srtmin, srtmax, pbeam, betann, betatar, betapro;
0473     G4double pbmin, pbmax;
0474 };
0475
0476 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0477
0478 struct ccurqmd13protarints
0479 {
0480     G4int pt_iso3[2][aamax], pt_ityp[2][aamax], pt_spin[2][aamax];
0481     G4int pt_charge[2][aamax], pt_aa[2], pt_uid[2][aamax];
0482 };
0483
0484 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0485
0486 struct ccurqmd13protarreals
0487 {
0488     G4double pt_r0[2][aamax], pt_rx[2][aamax], pt_ry[2][aamax], pt_rz[2][aamax], pt_fmass[2][aamax],
0489       pt_dectime[2][aamax];
0490     G4double pt_p0[2][aamax], pt_px[2][aamax], pt_py[2][aamax], pt_pz[2][aamax];
0491     G4double pt_rho[2][aamax];
0492     G4double pt_pmax[2][aamax];
0493 };
0494
0495 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0496
0497 // newpart
0498 struct ccurqmd13inewpart
0499 {
0500     G4int itypnew[mprt], i3new[mprt], itot[mprt], inew[mprt], nexit;
0501     G4int iline, strcount, pslot[oprt], nstring1, nstring2, sidnew[mprt], itypold[oprt],
0502       iso3old[oprt];
0503 };
0504
0505 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0506
0507 struct ccurqmd13rnewpart
0508 {
0509     G4double pnew[mprt][5], xnew[mprt][4], betax, betay, betaz, pold[oprt][5], p0nn, pxnn, pynn,
0510       pznn, pnn, mstring[2], pnnout, xtotfacold[oprt];
0511 };
0512
0513 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0514
0515 struct ccurqmd13fnewpart
0516 {
0517     G4double leadfac[mprt];
0518 };
0519
0520 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0521 //
0522 // boxinc
0523 //
0524 struct ccurqmd13boxic
0525 {
0526     G4int cbox;
0527     G4int boxflag;
0528     G4int mbox;
0529     G4int bptityp[bptmax], bptiso3[bptmax], bptpart[bptmax];
0530     G4int edensflag, para, solid, mbflag, mtest;
0531 };
0532
0533 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0534
0535 struct ccurqmd13boxrc
0536 {
0537     G4double bptpmax[bptmax];
0538     G4double edens;
0539     G4double lbox;
0540     G4double lboxhalbe;
0541     G4double lboxd;
0542     G4double mbp0, mbpx, mbpy, mbpz;
0543 };
0544
0545 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0546 // comnorm
0547 struct ccurqmd13normsplin
0548 {
0549     G4double x_norm[n][4], y_norm[n][4];
0550     G4double y2a[n][4], y2b[n][4], dx;
0551 };
0552
0553 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0554 // comstr
0555 struct ccurqmd13FRGSPA
0556 {
0557     G4double pjspns, pmix1s[njspin][3], pmix2s[njspin][3], pbars, parqls, parrs;
0558 };
0559
0560 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0561
0562 struct ccurqmd13FRGCPA
0563 {
0564     G4double pjspnc, pmix1c[njspin][3], pmix2c[njspin][3], pbarc;
0565 };
0566
0567 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0568
0569 struct ccurqmd13coparm
0570 {
0571     G4double parm[njspin];
0572 };
0573
0574 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0575
0576 struct ccurqmd13const
0577 {
0578     G4double pi;
0579 };
0580
0581 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0582 //// freezeout
0583 //
0584 struct ccurqmd13frcoor
0585 {
0586     G4double frr0[nmax], frrx[nmax], frry[nmax], frrz[nmax], frp0[nmax], frpx[nmax], frpy[nmax],
0587       frpz[nmax];
0588 };
0589
0590 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0591 //  input
0592 struct ccurqmd13values
0593 {
0594     G4double valint[1];
0595 };
0596
0597 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0598 // cascinit
0599 struct ccurqmd13ini
0600 {
0601     ftnlogical bcorr;
0602 };
0603
0604 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0605 // iso
0606 struct ccurqmd13factorials
0607 {
0608     G4double logfak[101];
0609 };
0610
0611 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0612 struct ccurqmd13cgks
0613 {
0614     G4double cgktab[jmax + 1][2 * jmax + 1][2 * jmax + 1][jmax + 1][jmax + 1];
0615 };
0616
0617 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0618 // UrQMD
0619 //
0620 struct ccurqmd13energies
0621 {
0622     G4double ekinbar, ekinmes, esky2, esky3, eyuk, ecb, epau;
0623 };
0624
0625 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0626 // urqmd
0627 extern "C"
0628 {
0629   extern int time_();
0630   extern void loginit_();
0631   extern void sseed_(int*);
0632   extern void uinit_(int*);
0633   extern void urqmd_();
0634   extern int pdgid_(int*, int*);  // ityp
0635
0636   extern void g4urqmdblockdata_();
0637
0638   // urqmdparams
0639   extern struct ccurqmd13urqmdparams urqmdparams_;
0640   // coms
0641   extern struct ccurqmd13sys sys_;
0642   extern struct ccurqmd13rsys rsys_;
0643   extern struct ccurqmd13comseed comseed_;
0644   extern struct ccurqmd13logic logic_;
0645   extern struct ccurqmd13mdprop mdprop_;
0646   extern struct ccurqmd13cuts cuts_;
0647   extern struct ccurqmd13spdata spdata_;
0648   extern struct ccurqmd13isys isys_;
0649   extern struct ccurqmd13coor coor_;
0650   extern struct ccurqmd13frag frag_;
0651   extern struct ccurqmd13aios aios_;
0652   extern struct ccurqmd13pots pots_;
0653   extern struct ccurqmd13scoor scoor_;
0654   extern struct ccurqmd13sisys sisys_;
0655   extern struct ccurqmd13ssys ssys_;
0656   extern struct ccurqmd13rtdelay rtdelay_;
0657   extern struct ccurqmd13itdelay itdelay_;
0658   extern struct ccurqmd13svinfo svinfo_;
0659   extern struct ccurqmd13ffermi ffermi_;
0660   extern struct ccurqmd13peq peq_;
0661   // comres
0662   extern struct ccurqmd13versioning versioning_;
0663   extern struct ccurqmd13resonances resonances_;
0664   extern struct ccurqmd13sigtabi sigtabi_;
0665   extern struct ccurqmd13sigtabr sigtabr_;
0666
0667   // comwid
0668   extern struct ccurqmd13decaywidth decaywidth_;
0669   extern struct ccurqmd13brwignorm brwignorm_;
0670   extern struct ccurqmd13xsections xsections_;
0671   extern struct ccurqmd13tabnames tabnames_;
0672   // options
0673   extern struct ccurqmd13options options_;
0674   extern struct ccurqmd13optstrings optstrings_;
0675   extern struct ccurqmd13loptions loptions_;
0676   extern struct ccurqmd13stables stables_;
0677   // colltab
0678   extern struct ccurqmd13colltab colltab_;
0679   // inputs
0680   extern struct ccurqmd13inputs inputs_;
0681   extern struct ccurqmd13input2 input2_;
0682   extern struct ccurqmd13protarints protarints_;
0683   extern struct ccurqmd13protarreals protarreals_;
0684   // newpart
0685   extern struct ccurqmd13inewpart inewpart_;
0686   extern struct ccurqmd13rnewpart rnewpart_;
0687   extern struct ccurqmd13fnewpart fnewpart_;
0688   // bocinc
0689   extern struct ccurqmd13boxic boxic_;
0690   extern struct ccurqmd13boxrc boxrc_;
0691   // comnorm
0692   struct ccurqmd13normsplin normsplin_;
0693   // comstr
0694   struct ccurqmd13FRGSPA FRGSPA_;
0695   struct ccurqmd13FRGCPA FRGCPA_;
0696   struct ccurqmd13coparm coparm_;
0697   struct ccurqmd13const const_;
0698   // freezeout
0699   struct ccurqmd13frcoor frcoor_;
0700   // urqmd
0701   extern struct ccurqmd13energies energies_;
0702   // input
0703   extern struct ccurqmd13values values_;
0704   // cascinit
0705   extern struct ccurqmd13ini ini_;
0706   // iso
0707   extern struct ccurqmd13factorials factorials_;
0708   extern struct ccurqmd13cgks cgks_;
0709 }
0710
0711 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0712 #endif