Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4UrQMD1_3Interface.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * *
21 // * Parts of this code which have been developed by Abdel-Waged *
22 // * et al under contract (31-465) to the King Abdul-Aziz City for *
23 // * Science and Technology (KACST), the National Centre of *
24 // * Mathematics and Physics (NCMP), Saudi Arabia. *
25 // * *
26 // * By using, copying, modifying or distributing the software (or *
27 // * any work based on the software) you agree to acknowledge its *
28 // * use in resulting scientific publications, and indicate your *
29 // * acceptance of all terms of the Geant4 Software license. *
30 // ********************************************************************
31 //
34 //
35 #ifndef G4UrQMD1_3Interface_hh
36 #define G4UrQMD1_3Interface_hh
37 
38 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
39 //
40 // MODULE: G4UrQMD1_3Model.hh
41 //
42 // Version: 0.B
43 // Date: 20/12/12
44 // Author: Kh. Abdel-Waged and Nuha Felemban
45 // Revised by: V.V. Uzhinskii
46 // SPONSERED BY
47 // Customer: KAUST/NCMP
48 // Contract: 31-465
49 //
50 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
51 //
52 //
53 // Class Description
54 //
55 //
56 // Class Description - End
57 //
58 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
60 #include "globals.hh"
61 // coms
62 //
63 const G4int nmax = 500;
64 const G4int nspl = 500;
65 const G4int smax = 500;
66 // comres
67 const G4int minnuc=1;
68 const G4int minmes=100;
69 const G4int maxmes=132;
70 const G4int numnuc=16;
71 const G4int numdel=10;
72 const G4int maxnuc=minnuc+numnuc-1;
73 const G4int mindel=minnuc+maxnuc;
74 const G4int maxdel=mindel+numdel-1;
75 const G4int minres=minnuc+1;
76 const G4int maxres=maxdel;
77 const G4int numlam=13;
78 const G4int numsig=9;
79 const G4int numcas=6;
80 const G4int numome=1;
81 const G4int minlam=mindel+numdel;
82 const G4int maxlam=minlam+numlam-1;
83 const G4int minsig=minlam+numlam;
84 const G4int maxsig=minsig+numsig-1;
85 const G4int mincas=minsig+numsig;
86 const G4int maxcas=mincas+numcas-1;
87 const G4int minome=mincas+numcas;
88 const G4int maxome=minome+numome-1;
89 const G4int minbar=minnuc;
90 const G4int maxbar=maxome;
91 const G4int offmeson=minmes;
92 const G4int maxmeson=maxmes;
93 const G4int maxbra=11;
94 const G4int maxbrm=25;
95 const G4int maxbrs1=10;
96 const G4int maxbrs2=3;
97 const G4int nsigs = 10;
98 const G4int itblsz= 100;
99 const G4int maxreac = 13;
100 const G4int maxpsig = 12;
101 //
102 //comwid
103 //
104 const G4int widnsp=120;
105 const G4double mintab=0.10;
106 const G4double maxtab1=5.0;
107 const G4double maxtab2=50.0;
108 const G4int tabver=9;
109 //
110 // options
111 //
112 const G4int numcto=400;
113 const G4int numctp=400;
114 const G4int maxstables=20;
115 //
116 // colltab (collision tables)
117 //
118 const G4int ncollmax = 100;
119 //
120 // inputs
121 //
122 const G4int aamax=300;
123 //
124 // newpart (new created particles)
125 //
126 const G4int mprt=200;
127 const G4int oprt=2;
128 //
129 // boxinc
130 //
131 const G4int bptmax=20;
132 //
133 
134 // This next line is required as the default version of FORTRAN LOGICAL is
135 // four bytes long, whereas storage for G4bool is one byte.
136 //
137 // comnorm
138 const G4int n = 400;
139 //
140 // comstr
141 const G4int njspin=8;
142 //
143 //iso
144 const G4int jmax=7;
145 
146 // This next line is required as the default version of FORTRAN LOGICAL is
147 // four bytes long, whereas storage for G4bool is one byte.
148 //
149 
150 typedef G4int ftnlogical;
151 
152 //
153 // Standard common block for UrQMD
154 // Common options for coms.f
155 // 20 commons
156 //
157 //
158 struct ccurqmd13urqmdparams
159 {
160 G4int u_at,u_zt,u_ap,u_zp;
161 G4double u_elab,u_imp;
162 G4int u_sptar,u_spproj;
163 };
164 
165 struct ccurqmd13sys
166 {
167 
168 G4int npart, nbar, nmes, ctag,nsteps,uid_cnt,
169  ranseed,event,ap,at,zp,zt,eos,dectag,
170  nhardres, nsoftres, ndecres, nelcoll, nblcoll;
171 };
172 
173 
174 struct ccurqmd13rsys
175 {
176 G4double time,acttime,bdist,bimp,bmin,ebeam,ecm;
177 };
178 
179 struct ccurqmd13comseed
180 {
181  ftnlogical
182  firstseed;
183 };
184 
185 
186 
187 struct ccurqmd13logic
188 {
189  ftnlogical
190  lsct[nmax], logSky, logYuk, logCb, logPau;
191 
192 };
193 
194 struct ccurqmd13mdprop
195 {
196  G4double
197  r0_t[nmax], rx_t[nmax], ry_t[nmax], rz_t[nmax];
198 };
199 
200 
201 struct ccurqmd13cuts
202 {
203  G4double
204  cutmax, cutPau, cutCb, cutYuk, cutSky, cutdww;
205 };
206 
207 
208 struct ccurqmd13spdata
209 {
210  G4double
211  spx[nspl], spPauy[nspl], outPau[nspl],
212  spCby[nspl], outCb[nspl],
213  spYuky[nspl], outYuk[nspl],
214  spSkyy[nspl], outSky[nspl],
215  spdwwy[nspl], outdww[nspl];
216 };
217 
218 
219 struct ccurqmd13isys
220 {
221 
222 G4int spin[nmax],ncoll[nmax],charge[nmax],ityp[nmax],
223  lstcoll[nmax],
224  iso3[nmax],origin[nmax],strid[nmax],uid[nmax];
225 };
226 
227 
228 struct ccurqmd13coor
229 {
230 G4double
231  r0[nmax], rx[nmax], ry[nmax], rz[nmax],
232  p0[nmax], px[nmax], py[nmax], pz[nmax],
233  fmass[nmax], rww[nmax],dectime[nmax];
234 };
235 
236 struct ccurqmd13frag
237 {
238 G4double
239  tform[nmax], xtotfac[nmax];
240 };
241 
242 struct ccurqmd13aios
243 {
244 G4double
245  airx[nmax], airy[nmax], airz[nmax],
246  aipx[nmax], aipy[nmax], aipz[nmax],
247  aorx [4][nmax], aory[4][nmax], aorz[4][nmax],
248  aopx[4][nmax], aopy[4][nmax], aopz[4][nmax];
249 };
250 
251 
252 struct ccurqmd13pots
253 {
254 
255 G4double
256  Cb0, Yuk0, Pau0, Sky20, Sky30, gamSky,
257  gamYuk, drPau, dpPau, gw, sgw, delr, fdel,
258  dt,da, db,dtimestep;
259 };
260 
261 struct ccurqmd13scoor
262 {
263 
264 G4double
265  r0s[smax], rxs[smax], rys[smax], rzs[smax],
266  p0s[smax], pxs[smax], pys[smax], pzs[smax],
267  sfmass[smax];
268 };
269 
270 struct ccurqmd13sisys
271 {
272  G4int
273  sspin[smax], scharge[smax], sityp[smax], siso3[smax],
274  suid[smax];
275 };
276 
277 
278 struct ccurqmd13ssys
279 {
280  G4int nspec;
281 };
282 
283 struct ccurqmd13rtdelay
284 {
285 G4double
286  p0td[nmax][2],pxtd[nmax][2],pytd[nmax][2],pztd[nmax][2],
287  fmasstd[nmax][2];
288 };
289 
290 
291 struct ccurqmd13itdelay
292 {
293 G4int
294  ityptd[nmax][2],iso3td[nmax][2];
295 };
296 
297 struct ccurqmd13svinfo
298 {
299 G4int
300  itypt[2],uidt[2],origint[2],iso3t[2];
301 };
302 
303 
304 struct ccurqmd13ffermi
305 {
306 G4double
307  ffermpx[nmax], ffermpy[nmax], ffermpz[nmax];
308 };
309 
310 struct ccurqmd13peq
311 {
312 G4double peq1, peq2;
313 };
314 
315 //
316 // Definition for Collision Term
317 // Commons comres
318 // 4 commons
319 //
320 
321 struct ccurqmd13versioning
322 {
323 char versiontag[45];
324 };
325 
326 
327 struct ccurqmd13resonances
328 {
329 
330  G4double massres[maxbar-minbar+1],widres[maxbar-minbar+1];
331  G4double massmes[maxmes-minmes+1];
332  G4double widmes[maxmes-minmes+1];
333  G4double mmesmn[maxmes-minmes+1];
334  G4double branres[maxdel-minnuc][maxbra+1];
335  G4double branmes[maxmes-minmes][maxbrm+1];
336 
337  G4double branbs1[maxsig-minlam][maxbrs1+1];
338  G4double branbs2[maxcas-mincas][maxbrs2+1];
339 
340  G4int bs1type[maxbrs1+1][4],bs2type[maxbrs2+1][4];
341  G4int lbs1[maxsig-minlam][maxbrs1+1];
342  G4int lbs2[maxcas-mincas][maxbrs2+1];
343  G4int lbm[maxmes-minmes][maxbrm+1];
344 
345  G4int jres[maxbar-minbar+1];
346  G4int jmes[maxmes-minmes+1];
347  G4int lbr[maxdel-minnuc][maxbra+1];
348  G4int brtype[maxbra+1][4];
349  G4int pares[maxbar-minbar+1],pames[maxmes-minmes+1];
350  G4int bmtype[maxbrm+1][4];
351  G4int isores[maxbar-minbar+1], isomes[maxmes-minmes+1];
352  G4int strres[maxbar-minbar+1],strmes[maxmes-minmes+1];
353  G4int mlt2it[maxmes-minmes];
354 };
355 
356 struct ccurqmd13sigtabi
357 {
358 G4int sigmaln[maxreac][2][maxpsig];
359 G4int sigmainf[20][nsigs];
360 };
361 
362 
363 struct ccurqmd13sigtabr
364 {
365 
366 G4double sigmas[itblsz][nsigs],sigmascal[5][nsigs];
367 };
368 
369 //comwid
370 struct ccurqmd13decaywidth
371 {
372 G4double tabx [widnsp];
373 G4double fbtaby [2][maxbar-minbar+1][widnsp];
374 G4double pbtaby[maxbra+1][maxbar-minbar+1][2][widnsp];
375 G4double fmtaby [2][maxmes-minmes+1][widnsp];
376 G4double pmtaby [maxbrm+1][maxmes-minmes+1][2][widnsp];
377 G4int wtabflg;
378 
379 };
380 
381 
382 struct ccurqmd13brwignorm
383 {
384 G4double bwbarnorm[maxbar-minbar+1];
385 G4double bwmesnorm[maxmes-minmes+1];
386 };
387 
388 
389 struct ccurqmd13xsections
390 {
391 G4double tabxnd [widnsp];
392 G4double frrtaby[maxdel-1][2][2][widnsp];
393 };
394 
395 
396 struct ccurqmd13tabnames
397 {
398 char tabname[77];
399 };
400 //-----------------
401 //
402 // options
403 //
404 struct ccurqmd13options
405 {
406 G4int CTOption[numcto];
407 
408 G4double CTParam[numctp];
409 };
410 
411 
412 struct ccurqmd13optstrings
413 {
414 char ctodc[numcto][2];
415 char ctpdc[numctp][2];
416 };
417 
418 struct ccurqmd13loptions
419 {
420 ftnlogical
421  fixedseed,bf13,bf14,bf15,bf16,bf17,bf18,bf19,
422  bf20;
423 };
424 
425 struct ccurqmd13stables
426 {
427 G4int nstable;
428 G4int stabvec[maxstables];
429 };
430 
431 //
432 //colltab
433 //
434 struct ccurqmd13colltab
435 {
436 G4double
437  cttime[ncollmax+1],ctsqrts[ncollmax],
438  ctsigtot[ncollmax],tmin;
439 G4int
440  cti1[ncollmax],cti2[ncollmax];
441 G4int
442  nct,actcol;
443 ftnlogical
444  ctvalid[ncollmax];
445 G4int
446  ctsav[ncollmax];
447 G4int
448  nsav,apt;
449 G4double
450  ctcolfluc[ncollmax];
451 };
452 
453 //
454 // inputs
455 //
456 struct ccurqmd13inputs
457 {
458 G4int nevents,spityp[2],prspflg;
459 G4int trspflg,spiso3[2],outsteps,bflag,srtflag,efuncflag;
460 G4int nsrt,firstev,npb;
461 };
462 
463 struct ccurqmd13input2
464 {
465 G4double srtmin,srtmax,pbeam,betann,betatar,betapro;
466 G4double pbmin,pbmax;
467 };
468 
469 struct ccurqmd13protarints
470 {
471 
472 G4int pt_iso3[2][aamax],pt_ityp[2][aamax],pt_spin[2][aamax];
473 
474 G4int pt_charge[2][aamax],pt_aa[2],pt_uid[2][aamax];
475 };
476 
477 struct ccurqmd13protarreals
478 {
479 G4double pt_r0[2][aamax],pt_rx[2][aamax],pt_ry[2][aamax],
480  pt_rz[2][aamax],pt_fmass[2][aamax],pt_dectime[2][aamax];
481 G4double pt_p0[2][aamax],pt_px[2][aamax],pt_py[2][aamax],
482  pt_pz[2][aamax];
483 G4double pt_rho[2][aamax];
484 G4double pt_pmax[2][aamax];
485 };
486 // newpart
487 struct ccurqmd13inewpart
488 {
489 G4int itypnew[mprt],i3new[mprt],itot[mprt],inew[mprt],nexit;
490 G4int iline,strcount,pslot[oprt],nstring1, nstring2,
491  sidnew[mprt],itypold[oprt],iso3old[oprt];
492 };
493 
494 struct ccurqmd13rnewpart
495 {
496 G4double pnew[mprt][5],xnew[mprt][4],betax,betay,betaz,
497  pold[oprt][5],p0nn,pxnn,pynn,pznn,pnn, mstring[2],
498  pnnout,xtotfacold[oprt];
499 
500 };
501 
502 struct ccurqmd13fnewpart
503 {
504 G4double leadfac[mprt];
505 };
506 //
507 // boxinc
508 //
509 struct ccurqmd13boxic
510 {
511 G4int cbox;
512 G4int boxflag;
513 G4int mbox;
514 G4int bptityp[bptmax],bptiso3[bptmax],bptpart[bptmax];
515 G4int edensflag,para,solid, mbflag,mtest;
516 };
517 
518 struct ccurqmd13boxrc
519 {
520 G4double bptpmax[bptmax];
521 G4double edens;
522 G4double lbox;
523 G4double lboxhalbe;
524 G4double lboxd;
525 G4double mbp0, mbpx, mbpy, mbpz;
526 };
527 // comnorm
528 struct ccurqmd13normsplin
529 {
530 G4double x_norm[n][4],y_norm[n][4];
531 G4double y2a[n][4],y2b[n][4], dx;
532 };
533 // comstr
534 struct ccurqmd13FRGSPA
535 {
536 G4double
537  pjspns, pmix1s[njspin][3], pmix2s[njspin][3]
538  , pbars, parqls, parrs;
539 };
540 struct ccurqmd13FRGCPA
541 {
542 G4double
543  pjspnc, pmix1c[njspin][3], pmix2c[njspin][3], pbarc;
544 };
545 
546 struct ccurqmd13coparm
547 {
548 G4double parm[njspin];
549 };
550 
551 struct ccurqmd13const
552 {
553 G4double pi;
554 };
556 //
557 struct ccurqmd13frcoor
558 {
559 G4double frr0[nmax], frrx[nmax], frry[nmax], frrz[nmax],
560  frp0[nmax], frpx[nmax], frpy[nmax], frpz[nmax];
561 };
562 
563 // input
564 struct ccurqmd13values
565 {
566 G4double valint[1];
567 };
568 
569 // cascinit
570 struct ccurqmd13ini
571 {
572 ftnlogical bcorr;
573 };
574 
575 
576 // iso
577 struct ccurqmd13factorials
578 {
579 G4double logfak[101];
580 };
581 struct ccurqmd13cgks
582 {
583 
584 G4double cgktab[jmax+1][2*jmax+1][2*jmax+1][jmax+1][jmax+1];
585 };
586 
587 // UrQMD
588 //
589 struct ccurqmd13energies
590 {
591 G4double ekinbar, ekinmes, esky2, esky3,
592  eyuk, ecb, epau;
593 };
594 
595 
596 
597 // urqmd
598 extern "C"
599 {
600 extern int time_ ();
601 extern void loginit_();
602 extern void sseed_ (int*);
603 extern void uinit_ (int*);
604 extern void urqmd_ ();
605 extern int pdgid_ (int*, int*); //ityp
606 
607 extern void g4urqmdblockdata_ ();
608 
609 // urqmdparams
610 extern struct ccurqmd13urqmdparams urqmdparams_;
611 //coms
612 extern struct ccurqmd13sys sys_;
613 extern struct ccurqmd13rsys rsys_;
614 extern struct ccurqmd13comseed comseed_;
615 extern struct ccurqmd13logic logic_;
616 extern struct ccurqmd13mdprop mdprop_;
617 extern struct ccurqmd13cuts cuts_;
618 extern struct ccurqmd13spdata spdata_;
619 extern struct ccurqmd13isys isys_;
620 extern struct ccurqmd13coor coor_;
621 extern struct ccurqmd13frag frag_;
622 extern struct ccurqmd13aios aios_;
623 extern struct ccurqmd13pots pots_;
624 extern struct ccurqmd13scoor scoor_;
625 extern struct ccurqmd13sisys sisys_;
626 extern struct ccurqmd13ssys ssys_;
627 extern struct ccurqmd13rtdelay rtdelay_;
628 extern struct ccurqmd13itdelay itdelay_;
629 extern struct ccurqmd13svinfo svinfo_;
630 extern struct ccurqmd13ffermi ffermi_;
631 extern struct ccurqmd13peq peq_;
632 //comres
633 extern struct ccurqmd13versioning versioning_;
634 extern struct ccurqmd13resonances resonances_;
635 extern struct ccurqmd13sigtabi sigtabi_;
636 extern struct ccurqmd13sigtabr sigtabr_;
637 
638 //comwid
639 extern struct ccurqmd13decaywidth decaywidth_;
640 extern struct ccurqmd13brwignorm brwignorm_;
641 extern struct ccurqmd13xsections xsections_;
642 extern struct ccurqmd13tabnames tabnames_;
643 //options
644 extern struct ccurqmd13options options_;
645 extern struct ccurqmd13optstrings optstrings_;
646 extern struct ccurqmd13loptions loptions_;
647 extern struct ccurqmd13stables stables_;
648 //colltab
649 extern struct ccurqmd13colltab colltab_;
650 //inputs
651 extern struct ccurqmd13inputs inputs_;
652 extern struct ccurqmd13input2 input2_;
653 extern struct ccurqmd13protarints protarints_;
654 extern struct ccurqmd13protarreals protarreals_;
655 //newpart
656 extern struct ccurqmd13inewpart inewpart_;
657 extern struct ccurqmd13rnewpart rnewpart_;
658 extern struct ccurqmd13fnewpart fnewpart_;
659 //bocinc
660 extern struct ccurqmd13boxic boxic_;
661 extern struct ccurqmd13boxrc boxrc_;
662 // comnorm
663 struct ccurqmd13normsplin normsplin_;
664 //comstr
665 struct ccurqmd13FRGSPA FRGSPA_;
666 struct ccurqmd13FRGCPA FRGCPA_;
667 struct ccurqmd13coparm coparm_;
668 struct ccurqmd13const const_;
669 // freezeout
670 struct ccurqmd13frcoor frcoor_;
671 //urqmd
672 extern struct ccurqmd13energies energies_;
673 //input
674 extern struct ccurqmd13values values_;
675 // cascinit
676 extern struct ccurqmd13ini ini_;
677 //iso
678 extern struct ccurqmd13factorials factorials_;
679 extern struct ccurqmd13cgks cgks_;
680 
681 }
682 
683 #endif
Generated on Sat May 25 2013 14:32:56 for Geant4 by   doxygen 1.8.4