Geant4  10.00.p03
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 
78  void setVerboseLevel(G4int level);
79 
84  return volant;
85  }
86 
101  void breakItUp(G4int nucleusA, G4int nucleusZ, G4double nucleusMass, G4double excitationEnergy,
102  G4double angularMomentum, G4double recoilEnergy, G4double momX, G4double momY, G4double momZ,
103  G4int eventnumber);
104 
105  // Evaporation
106 public:
111  void initEvapora();
112 
125  void qrot(G4double z, G4double a, G4double bet, G4double sig, G4double u, G4double *qr);
126 
131  void mglw(G4double a, G4double z, G4double *el);
132 
136  void mglms(G4double a, G4double z, G4int refopt4, G4double *el);
137 
141  G4double spdef(G4int a, G4int z, G4int optxfis);
142 
146  G4double fissility(int a,int z, int optxfis);
147 
151  void evapora(G4double zprf, G4double aprf, G4double *ee_par, G4double jprf,
152  G4double *zf_par, G4double *af_par, G4double *mtota_par,
153  G4double *pleva_par, G4double *pxeva_par, G4double *pyeva_par,
154  G4int *ff_par, G4int *inttype_par, G4int *inum_par);
155 
159  void direct(G4double zprf,G4double a, G4double ee, G4double jprf,
160  G4double *probp_par, G4double *probn_par, G4double *proba_par,
161  G4double *probf_par, G4double *ptotl_par, G4double *sn_par, G4double *sbp_par, G4double *sba_par, G4double *ecn_par,
162  G4double *ecp_par,G4double *eca_par, G4double *bp_par, G4double *ba_par, G4int, G4int inum, G4int itest);
163 
167  void densniv(G4double a, G4double z, G4double ee, G4double esous, G4double *dens, G4double bshell, G4double bs, G4double bk,
168  G4double *temp, G4int optshp, G4int optcol, G4double defbet);
169 
176  G4double bfms67(G4double zms, G4double ams);
177 
185  void lpoly(G4double x, G4int n, G4double pl[]);
186 
193  G4double eflmac(G4int ia, G4int iz, G4int flag, G4int optshp);
194 
199  void appariem(G4double a, G4double z, G4double *del);
200 
205  void parite(G4double n, G4double *par);
206 
211  G4double tau(G4double bet, G4double homega, G4double ef, G4double t);
212 
217  G4double cram(G4double bet, G4double homega);
218 
224  G4double bipol(int iflag, G4double y);
225 
235  void barfit(G4int iz, G4int ia, G4int il, G4double *sbfis, G4double *segs, G4double *selmax);
236 
240  G4double haz(G4int k);
241  void standardRandom(G4double *rndm, G4long *seed);
242 
247 
251  G4double fd(G4double E);
252 
256  G4double f(G4double E);
257 
262 
267 
271  void guet(G4double *x_par, G4double *z_par, G4double *find_par);
272 
273 public:
274  // Coordinate system transformations:
275  void lorab(G4double gam, G4double eta, G4double ein, G4double pin[],
276  G4double *eout, G4double pout[]);
277 
278  void translab(G4double gamrem, G4double etrem, G4double csrem[4], G4int nopart, G4int ndec);
279  void translabpf(G4double masse1, G4double t1, G4double p1, G4double ctet1,
280  G4double phi1, G4double gamrem, G4double etrem, G4double R[][4],
281  G4double *plab1, G4double *gam1, G4double *eta1, G4double csdir[]);
282 
283  void rotab(G4double R[4][4], G4double pin[4], G4double pout[4]);
284 
285  // Utils
286  G4int min(G4int a, G4int b);
288  G4int max(G4int a, G4int b);
290 
291  G4int nint(G4double number);
292  G4int secnds(G4int x);
293  G4int mod(G4int a, G4int b);
297  G4int idnint(G4double value);
301  return ec2sub;
302  }
303 
304 private:
307 
319 #ifndef ABLAXX_IN_GEANT4_MODE
320  G4INCL::Config *theConfig;
321 #endif
322 };
323 
324 #endif
G4Ec2sub * ec2sub
Definition: G4Abla.hh:312
G4Ecld * ecld
Definition: G4Abla.hh:313
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:67
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:83
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:300
G4Ald * ald
Definition: G4Abla.hh:310
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:318
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:316
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:314
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:305
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:311
Evaporation and fission output data.
Shell corrections and deformations.
G4Fiss * fiss
Definition: G4Abla.hh:315
G4double iz
Definition: TRTMaterials.hh:39
G4int ilast
Definition: G4Abla.hh:306
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
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:317
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:308
double G4double
Definition: G4Types.hh:76
G4double expohaz(G4int k, G4double T)
TIRAGE ALEATOIRE DANS UNE EXPONENTIELLLE : Y=EXP(-X/T)
Definition: G4Abla.cc:3012
G4Pace * pace
Definition: G4Abla.hh:309
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