Geant4_10
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);
301  return ec2sub;
302  }
303 
304 private:
305  G4int verboseLevel;
306  G4int ilast;
307 
308  G4AblaFissionBase *fissionModel;
309  G4Pace *pace;
310  G4Ald *ald;
311  G4Eenuc *eenuc;
312  G4Ec2sub *ec2sub;
313  G4Ecld *ecld;
314  G4Fb *fb;
315  G4Fiss *fiss;
316  G4Opt *opt;
317  G4Volant *volant;
318  G4VarNtp *varntp;
319 #ifndef ABLAXX_IN_GEANT4_MODE
320  G4INCL::Config *theConfig;
321 #endif
322 };
323 
324 #endif
Definition: G4Abla.hh:54
void parite(G4double n, G4double *par)
Definition: G4Abla.cc:2561
G4double eflmac(G4int ia, G4int iz, G4int flag, G4int optshp)
Definition: G4Abla.cc:2385
void setVerboseLevel(G4int level)
Definition: G4Abla.cc:76
G4Abla(G4Volant *aVolant, G4VarNtp *aVarntp)
Definition: G4Abla.cc:45
TTree * t1
Definition: plottest35.C:26
G4Volant * getVolant()
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
tuple a
Definition: test.py:11
G4double utilabs(G4double a)
Definition: G4Abla.cc:3733
DOUBLE PRECISION function rndm(RDUMMY)
Definition: dpm25nulib.f:1460
G4double haz(G4int k)
Definition: G4Abla.cc:3520
G4Ec2sub * getFrldmTable()
Definition: G4Abla.hh:300
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)
Definition: G4Abla.cc:971
long G4long
Definition: G4Types.hh:80
G4double fmaxhaz(G4double T)
Definition: G4Abla.cc:3032
void lpoly(G4double x, G4int n, G4double pl[])
Definition: G4Abla.cc:2367
tuple x
Definition: test.py:50
function bk(X)
Definition: hijing1.383.f:5001
G4double fissility(int a, int z, int optxfis)
Definition: G4Abla.cc:1064
G4double cram(G4double bet, G4double homega)
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)
Definition: G4Abla.cc:1457
int G4int
Definition: G4Types.hh:78
void mglw(G4double a, G4double z, G4double *el)
Definition: G4Abla.cc:944
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)
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)
Definition: G4Abla.cc:2127
G4double spdef(G4int a, G4int z, G4int optxfis)
Definition: G4Abla.cc:1014
Double_t y
Definition: plot.C:279
void breakItUp(G4int nucleusA, G4int nucleusZ, G4double nucleusMass, G4double excitationEnergy, G4double angularMomentum, G4double recoilEnergy, G4double momX, G4double momY, G4double momZ, G4int eventnumber)
Definition: G4Abla.cc:102
tuple b
Definition: test.py:12
Char_t n[5]
tuple pl
Definition: readPY.py:5
G4double iz
Definition: TRTMaterials.hh:39
G4int idnint(G4double value)
Definition: G4Abla.cc:3713
G4int idint(G4double a)
Definition: G4Abla.cc:3699
~G4Abla()
Definition: G4Abla.cc:82
void standardRandom(G4double *rndm, G4long *seed)
G4double bipol(int iflag, G4double y)
Definition: G4Abla.cc:2632
void translab(G4double gamrem, G4double etrem, G4double csrem[4], G4int nopart, G4int ndec)
Definition: G4Abla.cc:3202
G4double dmin1(G4double a, G4double b, G4double c)
Definition: G4Abla.cc:3719
void initEvapora()
Definition: G4Abla.cc:682
G4double pace2(G4double a, G4double z)
Definition: G4Abla.cc:3093
G4int max(G4int a, G4int b)
Definition: G4Abla.cc:3621
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
Definition: G4Abla.cc:2586
tuple z
Definition: test.py:28
G4double bfms67(G4double zms, G4double ams)
Definition: G4Abla.cc:2350
const XML_Char int const XML_Char * value
Definition: expat.h:331
G4double f(G4double E)
Definition: G4Abla.cc:3026
void guet(G4double *x_par, G4double *z_par, G4double *find_par)
Definition: G4Abla.cc:3131
double G4double
Definition: G4Types.hh:76
tuple c
Definition: test.py:13
G4double expohaz(G4int k, G4double T)
Definition: G4Abla.cc:3012
void barfit(G4int iz, G4int ia, G4int il, G4double *sbfis, G4double *segs, G4double *selmax)
Definition: G4Abla.cc:2687
void appariem(G4double a, G4double z, G4double *del)
Definition: G4Abla.cc:2531
void qrot(G4double z, G4double a, G4double bet, G4double sig, G4double u, G4double *qr)
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)
Definition: G4Abla.cc:3019