Geant4  10.01.p02
G4Abla.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 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // ABLAXX statistical de-excitation model
27 // Pekka Kaitaniemi, HIP (translation)
28 // Christelle Schmidt, IPNL (fission code)
29 // Davide Mancusi, CEA (contact person INCL/ABLA)
30 // Aatos Heikkinen, HIP (project coordination)
31 //
32 #define ABLAXX_IN_GEANT4_MODE 1
33 
34 #include "globals.hh"
35 
36 #ifndef G4Abla_hh
37 #define G4Abla_hh 1
38 
39 #ifdef ABLAXX_IN_GEANT4_MODE
40 #include "globals.hh"
41 #else
42 #include "G4INCLGeant4Compat.hh"
43 #include "G4INCLConfig.hh"
44 #endif
45 
46 #include "G4AblaRandom.hh"
47 #include "G4AblaDataDefs.hh"
48 #include "G4AblaFissionBase.hh"
49 
54 class G4Abla {
55 
56 public:
64 #ifdef ABLAXX_IN_GEANT4_MODE
65  G4Abla(G4Volant *aVolant, G4VarNtp *aVarntp);
66 #else
67  G4Abla(G4INCL::Config *config, G4Volant *aVolant, G4VarNtp *aVarntp);
68 #endif
69 
73  ~G4Abla();
74 
76  G4Abla(G4Abla const &other);
77 
79  G4Abla &operator=(G4Abla const &other);
80 
84  void setVerboseLevel(G4int level);
85 
90  return volant;
91  }
92 
107  void breakItUp(G4int nucleusA, G4int nucleusZ, G4double nucleusMass, G4double excitationEnergy,
108  G4double angularMomentum, G4double recoilEnergy, G4double momX, G4double momY, G4double momZ,
109  G4int eventnumber);
110 
111  // Evaporation
112 public:
117  void initEvapora();
118 
131  void qrot(G4double z, G4double a, G4double bet, G4double sig, G4double u, G4double *qr);
132 
137  void mglw(G4double a, G4double z, G4double *el);
138 
142  void mglms(G4double a, G4double z, G4int refopt4, G4double *el);
143 
147  G4double spdef(G4int a, G4int z, G4int optxfis);
148 
152  G4double fissility(int a,int z, int optxfis);
153 
157  void evapora(G4double zprf, G4double aprf, G4double *ee_par, G4double jprf,
158  G4double *zf_par, G4double *af_par, G4double *mtota_par,
159  G4double *pleva_par, G4double *pxeva_par, G4double *pyeva_par,
160  G4int *ff_par, G4int *inttype_par, G4int *inum_par);
161 
165  void direct(G4double zprf,G4double a, G4double ee, G4double jprf,
166  G4double *probp_par, G4double *probn_par, G4double *proba_par,
167  G4double *probf_par, G4double *ptotl_par, G4double *sn_par, G4double *sbp_par, G4double *sba_par, G4double *ecn_par,
168  G4double *ecp_par,G4double *eca_par, G4double *bp_par, G4double *ba_par, G4int, G4int inum, G4int itest);
169 
173  void densniv(G4double a, G4double z, G4double ee, G4double esous, G4double *dens, G4double bshell, G4double bs, G4double bk,
174  G4double *temp, G4int optshp, G4int optcol, G4double defbet);
175 
182  G4double bfms67(G4double zms, G4double ams);
183 
191  void lpoly(G4double x, G4int n, G4double pl[]);
192 
199  G4double eflmac(G4int ia, G4int iz, G4int flag, G4int optshp);
200 
205  void appariem(G4double a, G4double z, G4double *del);
206 
211  void parite(G4double n, G4double *par);
212 
217  G4double tau(G4double bet, G4double homega, G4double ef, G4double t);
218 
223  G4double cram(G4double bet, G4double homega);
224 
230  G4double bipol(int iflag, G4double y);
231 
241  void barfit(G4int iz, G4int ia, G4int il, G4double *sbfis, G4double *segs, G4double *selmax);
242 
246  G4double haz(G4int k);
247  void standardRandom(G4double *rndm, G4long *seed);
248 
253 
257  G4double fd(G4double E);
258 
262  G4double f(G4double E);
263 
268 
273 
277  void guet(G4double *x_par, G4double *z_par, G4double *find_par);
278 
279 public:
280  // Coordinate system transformations:
281  void lorab(G4double gam, G4double eta, G4double ein, G4double pin[],
282  G4double *eout, G4double pout[]);
283 
284  void translab(G4double gamrem, G4double etrem, G4double csrem[4], G4int nopart, G4int ndec);
285  void translabpf(G4double masse1, G4double t1, G4double p1, G4double ctet1,
286  G4double phi1, G4double gamrem, G4double etrem, G4double R[][4],
287  G4double *plab1, G4double *gam1, G4double *eta1, G4double csdir[]);
288 
289  void rotab(G4double R[4][4], G4double pin[4], G4double pout[4]);
290 
291  // Utils
292  G4int min(G4int a, G4int b);
294  G4int max(G4int a, G4int b);
296 
297  G4int nint(G4double number);
298  G4int secnds(G4int x);
299  G4int mod(G4int a, G4int b);
303  G4int idnint(G4double value);
307  return ec2sub;
308  }
309 
310 private:
313 
325 #ifndef ABLAXX_IN_GEANT4_MODE
326  G4INCL::Config *theConfig;
327 #endif
328 };
329 
330 #endif
G4Ec2sub * ec2sub
Definition: G4Abla.hh:318
G4Ecld * ecld
Definition: G4Abla.hh:319
Class containing ABLA de-excitation code.
Definition: G4Abla.hh:54
void parite(G4double n, G4double *par)
PROCEDURE FOR CALCULATING THE PARITY OF THE NUMBER N.
Definition: G4Abla.cc:2561
G4double eflmac(G4int ia, G4int iz, G4int flag, G4int optshp)
This function will calculate the liquid-drop nuclear mass for spheri configuration according to the p...
Definition: G4Abla.cc:2385
void setVerboseLevel(G4int level)
Set verbosity level.
Definition: G4Abla.cc:76
The INCL configuration object.
Definition: G4INCLConfig.hh:60
G4Abla(G4Volant *aVolant, G4VarNtp *aVarntp)
This constructor is used by standalone test driver and the Geant4 interface.
Definition: G4Abla.cc:45
G4Volant * getVolant()
Get the internal output data structure pointer.
Definition: G4Abla.hh:89
void lorab(G4double gam, G4double eta, G4double ein, G4double pin[], G4double *eout, G4double pout[])
Definition: G4Abla.cc:3489
G4int nint(G4double number)
Definition: G4Abla.cc:3631
G4double utilabs(G4double a)
Definition: G4Abla.cc:3733
G4double z
Definition: TRTMaterials.hh:39
G4double haz(G4int k)
Random numbers.
Definition: G4Abla.cc:3520
G4Ec2sub * getFrldmTable()
Definition: G4Abla.hh:306
G4Ald * ald
Definition: G4Abla.hh:316
void translabpf(G4double masse1, G4double t1, G4double p1, G4double ctet1, G4double phi1, G4double gamrem, G4double etrem, G4double R[][4], G4double *plab1, G4double *gam1, G4double *eta1, G4double csdir[])
Definition: G4Abla.cc:3428
void mglms(G4double a, G4double z, G4int refopt4, G4double *el)
Mglms.
Definition: G4Abla.cc:971
long G4long
Definition: G4Types.hh:80
G4VarNtp * varntp
Definition: G4Abla.hh:324
G4double fmaxhaz(G4double T)
tirage aleatoire dans une maxwellienne
Definition: G4Abla.cc:3032
void lpoly(G4double x, G4int n, G4double pl[])
This subroutine calculates the ordinary legendre polynomials of order 0 to n-1 of argument x and stor...
Definition: G4Abla.cc:2367
G4double a
Definition: TRTMaterials.hh:39
G4double fissility(int a, int z, int optxfis)
Calculation of fissility parameter.
Definition: G4Abla.cc:1064
G4double cram(G4double bet, G4double homega)
KRAMERS FAKTOR - REDUCTION OF THE FISSION PROBABILITY INDEPENDENT OF EXCITATION ENERGY.
Definition: G4Abla.cc:2615
void direct(G4double zprf, G4double a, G4double ee, G4double jprf, G4double *probp_par, G4double *probn_par, G4double *proba_par, G4double *probf_par, G4double *ptotl_par, G4double *sn_par, G4double *sbp_par, G4double *sba_par, G4double *ecn_par, G4double *ecp_par, G4double *eca_par, G4double *bp_par, G4double *ba_par, G4int, G4int inum, G4int itest)
Calculation of particle emission probabilities.
Definition: G4Abla.cc:1457
G4Opt * opt
Definition: G4Abla.hh:322
int G4int
Definition: G4Types.hh:78
void mglw(G4double a, G4double z, G4double *el)
Model de la goutte liquide de c.
Definition: G4Abla.cc:944
G4Fb * fb
Definition: G4Abla.hh:320
G4int mod(G4int a, G4int b)
Definition: G4Abla.cc:3675
void evapora(G4double zprf, G4double aprf, G4double *ee_par, G4double jprf, G4double *zf_par, G4double *af_par, G4double *mtota_par, G4double *pleva_par, G4double *pxeva_par, G4double *pyeva_par, G4int *ff_par, G4int *inttype_par, G4int *inum_par)
Main evaporation routine.
Definition: G4Abla.cc:1098
void densniv(G4double a, G4double z, G4double ee, G4double esous, G4double *dens, G4double bshell, G4double bs, G4double bk, G4double *temp, G4int optshp, G4int optcol, G4double defbet)
Level density parameters.
Definition: G4Abla.cc:2127
G4double spdef(G4int a, G4int z, G4int optxfis)
Definition: G4Abla.cc:1014
G4int verboseLevel
Definition: G4Abla.hh:311
void breakItUp(G4int nucleusA, G4int nucleusZ, G4double nucleusMass, G4double excitationEnergy, G4double angularMomentum, G4double recoilEnergy, G4double momX, G4double momY, G4double momZ, G4int eventnumber)
Main interface to the de-excitation code.
Definition: G4Abla.cc:102
Options.
G4Eenuc * eenuc
Definition: G4Abla.hh:317
Evaporation and fission output data.
Shell corrections and deformations.
G4Fiss * fiss
Definition: G4Abla.hh:321
G4double iz
Definition: TRTMaterials.hh:39
const G4double p1
G4int ilast
Definition: G4Abla.hh:312
G4int idnint(G4double value)
Definition: G4Abla.cc:3713
G4int idint(G4double a)
Definition: G4Abla.cc:3699
const G4int n
~G4Abla()
Basic destructor.
Definition: G4Abla.cc:82
void standardRandom(G4double *rndm, G4long *seed)
G4double bipol(int iflag, G4double y)
CALCULATION OF THE SURFACE BS OR CURVATURE BK OF A NUCLEUS RELATIVE TO THE SPHERICAL CONFIGURATION BA...
Definition: G4Abla.cc:2632
void translab(G4double gamrem, G4double etrem, G4double csrem[4], G4int nopart, G4int ndec)
Definition: G4Abla.cc:3202
G4Abla & operator=(G4Abla const &other)
Dummy assignment operator.
Masses.
G4double dmin1(G4double a, G4double b, G4double c)
Definition: G4Abla.cc:3719
void initEvapora()
Initialize ABLA evaporation code.
Definition: G4Abla.cc:682
G4double pace2(G4double a, G4double z)
Definition: G4Abla.cc:3093
Fission barriers.
G4Volant * volant
Definition: G4Abla.hh:323
G4int max(G4int a, G4int b)
Definition: G4Abla.cc:3621
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
RISE TIME IN WHICH THE FISSION WIDTH HAS REACHED 90 PERCENT OF ITS FINAL VALUE.
Definition: G4Abla.cc:2586
G4double bfms67(G4double zms, G4double ams)
This subroutine calculates the fission barriers of the liquid-drop model of Myers and Swiatecki (1967...
Definition: G4Abla.cc:2350
G4double f(G4double E)
FONCTION INTEGRALE DE FD(E)
Definition: G4Abla.cc:3026
void guet(G4double *x_par, G4double *z_par, G4double *find_par)
Definition: G4Abla.cc:3131
G4AblaFissionBase * fissionModel
Definition: G4Abla.hh:314
double G4double
Definition: G4Types.hh:76
static unsigned bk
Definition: csz_inflate.cc:360
G4double expohaz(G4int k, G4double T)
TIRAGE ALEATOIRE DANS UNE EXPONENTIELLLE : Y=EXP(-X/T)
Definition: G4Abla.cc:3012
G4Pace * pace
Definition: G4Abla.hh:315
void barfit(G4int iz, G4int ia, G4int il, G4double *sbfis, G4double *segs, G4double *selmax)
THIS SUBROUTINE RETURNS THE BARRIER HEIGHT BFIS, THE GROUND-STATE ENERGY SEGS, IN MEV...
Definition: G4Abla.cc:2687
void appariem(G4double a, G4double z, G4double *del)
Procedure for calculating the pairing correction to the binding energy of a specific nucleus...
Definition: G4Abla.cc:2531
void qrot(G4double z, G4double a, G4double bet, G4double sig, G4double u, G4double *qr)
Coefficient of collective enhancement including damping Input: z,a,bet,sig,u Output: qr - collective ...
Definition: G4Abla.cc:898
G4int secnds(G4int x)
Definition: G4Abla.cc:3659
G4int min(G4int a, G4int b)
Definition: G4Abla.cc:3601
G4double dmod(G4double a, G4double b)
void rotab(G4double R[4][4], G4double pin[4], G4double pout[4])
Definition: G4Abla.cc:3505
G4double dint(G4double a)
Definition: G4Abla.cc:3685
G4double fd(G4double E)
DISTRIBUTION DE MAXWELL.
Definition: G4Abla.cc:3019