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
343 if(
full[particleIndex]) {
346 avv[particleIndex] = (
int) A;
347 zvv[particleIndex] = (
int) Z;
348 enerj[particleIndex] = E;
349 plab[particleIndex] = P;
350 tetlab[particleIndex] = theta;
351 philab[particleIndex] = phi;
352 full[particleIndex] =
true;
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 dm[PACESIZEROWS][PACESIZECOLS]
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
G4double plab[VARNTPSIZE]
G4bool copied[VOLANTSIZE]
G4double getTotalEnergy()
G4double vgsld[ECLDROWS][ECLDCOLS]
G4double ecnz[EC2SUBROWS][EC2SUBCOLS]
G4double enerj[VARNTPSIZE]
G4double getMass(G4int particle)
G4double zpcv[VOLANTSIZE]
G4double alpha[ECLDROWS][ECLDCOLS]
G4double xhe[XHESIZE][EENUCSIZE]
G4double tetlab[VARNTPSIZE]
G4int itypcasc[VARNTPSIZE]
G4double efa[FBCOLS][FBROWS]
G4double total(Particle const *const p1, Particle const *const p2)
G4double ecgnz[ECLDROWS][ECLDCOLS]
G4int getTotalBaryonNumber()
G4double getTotalThreeMomentum()
G4double ecfnz[ECLDROWS][ECLDCOLS]
G4double getMomentumSum()
G4double philab[VARNTPSIZE]
void addParticle(G4double A, G4double Z, G4double E, G4double P, G4double theta, G4double phi)