32 #define ABLAXX_IN_GEANT4_MODE 1
38 #ifndef G4AblaDataDefs_hh
39 #define G4AblaDataDefs_hh 1
41 #ifdef ABLAXX_IN_GEANT4_MODE
60 #define PACESIZEROWS 500
61 #define PACESIZECOLS 500
76 #define EC2SUBROWS 154
198 #define EENUCSIZE 2002
218 #define VOLANTSIZE 301
250 for(
G4int i = 0; i <=
iv; i++) {
258 G4double totA = 0.0, totZ = 0.0, totP = 0.0;
260 for(
G4int i = 0; i <=
iv; i++) {
261 if(i == 0 &&
acv[i] != 0) {
279 #define VARNTPSIZE 301
353 ntrack = particleIndex + 1;
362 G4int baryonNumber = 0;
365 baryonNumber +=
avv[i];
377 energy += std::sqrt(std::pow(
plab[i], 2) + std::pow(
getMass(i), 2));
403 const G4double protonMass = 938.272;
404 const G4double neutronMass = 939.565;
408 if(
avv[particle] == 1 &&
zvv[particle] == 1) mass = protonMass;
409 if(
avv[particle] == 1 &&
zvv[particle] == 0) mass = neutronMass;
410 if(
avv[particle] == -1) mass = pionMass;
411 if(
avv[particle] > 1)
412 mass =
avv[particle] * protonMass +
zvv[particle] * neutronMass;
420 G4int nProton = 0, nNeutron = 0;
421 G4int nPiPlus = 0, nPiZero = 0, nPiMinus = 0;
422 G4int nH2 = 0, nHe3 = 0, nAlpha = 0;
423 G4int nFragments = 0;
424 G4int nParticles = 0;
427 if(
avv[i] == 1 &&
zvv[i] == 1) nProton++;
428 if(
avv[i] == 1 &&
zvv[i] == 0) nNeutron++;
429 if(
avv[i] == -1 &&
zvv[i] == 1) nPiPlus++;
430 if(
avv[i] == -1 &&
zvv[i] == 0) nPiZero++;
431 if(
avv[i] == -1 &&
zvv[i] == -1) nPiMinus++;
432 if(
avv[i] == 2 &&
zvv[i] == 1) nH2++;
433 if(
avv[i] == 3 &&
zvv[i] == 2) nHe3++;
434 if(
avv[i] == 4 &&
zvv[i] == 2) nAlpha++;
435 if(
zvv[i] > 2) nFragments++;
G4double spectatorEx
Spectator nucleus excitation energy for light ion projectile support.
G4int targetZ
Target charge number.
G4double dm[PACESIZEROWS][PACESIZECOLS]
G4int spectatorZ
Spectator nucleus charge number for light ion projectile support.
G4int projType
Projectile type.
G4double plab[VARNTPSIZE]
Momentum.
G4int spectatorA
Spectator nucleus mass number for light ion projectile support.
G4bool copied[VOLANTSIZE]
G4bool full[VARNTPSIZE]
The state of the index: true = reserved false = free.
G4double getTotalEnergy()
Return total energy.
G4int avv[VARNTPSIZE]
A (-1 for pions).
void dump()
Dump debugging output.
G4double vgsld[ECLDROWS][ECLDCOLS]
Difference between deformed ground state and ldm value.
G4double ecnz[EC2SUBROWS][EC2SUBCOLS]
G4double enerj[VARNTPSIZE]
Kinetic energy.
G4bool needsFermiBreakup
Does this nucleus require Fermi break-up treatment? Only applicable when used together with Geant4...
G4int izfis
Z of fiss nucleus.
G4double getMass(G4int particle)
G4int targetA
Target mass number.
G4double spectatorP2
Spectator nucleus momentum y-component.
double A(double temperature)
Evaporation and fission output data.
Shell corrections and deformations.
G4double spectatorT
Spectator nucleus kinetic energy.
G4double zpcv[VOLANTSIZE]
G4double alpha[ECLDROWS][ECLDCOLS]
Alpha ground state deformation (this is not beta2!) beta2 = std::sqrt(5/(4pi)) * alpha.
G4int jremn
Remnant Intrinsic Spin.
G4double massini
A of the remnant.
G4double xhe[XHESIZE][EENUCSIZE]
G4double tetlab[VARNTPSIZE]
Theta angle.
G4int itypcasc[VARNTPSIZE]
emitted in cascade (0) or evaporation (1).
G4Fiss()
Options and parameters for fission channel.
G4double masp
Projectile spectator A, Z, Eex;.
void dump()
Dump the contents of the ecnz data table.
G4double spectatorP1
Spectator nucleus momentum x-component.
G4double mzini
Z of the remnant.
G4int kfis
Fission 1/0=Y/N.
G4int ntrack
Number of particles.
G4double efa[FBCOLS][FBROWS]
G4double total(Particle const *const p1, Particle const *const p2)
G4double ecgnz[ECLDROWS][ECLDCOLS]
Ground state shell correction frldm for a spherical ground state.
G4int getTotalBaryonNumber()
Baryon number conservation check.
G4double energy(const ThreeVector &p, const G4double m)
G4int mulnevap
Evaporation n multip.
G4double getTotalThreeMomentum()
Return total three momentum.
G4double ecfnz[ECLDROWS][ECLDCOLS]
Shell correction for the saddle point (now: == 0).
G4double getMomentumSum()
void clear()
Clear and initialize all variables and arrays.
G4double philab[VARNTPSIZE]
Phi angle.
G4double spectatorM
Spectator nucleus mass.
G4double bimpact
Impact parameter.
G4int iafis
A of fiss nucleus.
void addParticle(G4double A, G4double Z, G4double E, G4double P, G4double theta, G4double phi)
Add a particle to the INCL/ABLA final output.
G4double spectatorP3
Spectator nucleus momentum z-component.
G4double estfis
Excit energy at fis.
G4int mulncasc
Cascade n multip.
G4double projEnergy
Projectile energy.
G4int mulntot
Total n multip.
G4double exini
Excitation energy.