32 #define ABLAXX_IN_GEANT4_MODE 1
36 #ifdef ABLAXX_IN_GEANT4_MODE
53 theABLAModel(new
G4Abla(volant, ablaResult)),
76 const G4double eKinRem = (eTotRem - pRem.invariantMass()) /
MeV;
103 result->push_back(product);
116 else if(A == 3 && Z == 2)
return G4He3::He3();
118 else if(A > 0 && Z > 0 && A >= Z) {
121 G4cout <<
"Can't convert particle with A=" << A <<
", Z=" << Z <<
" to G4ParticleDefinition, trouble ahead" <<
G4endl;
143 #endif // ABLAXX_IN_GEANT4_MODE
Class containing ABLA de-excitation code.
CLHEP::Hep3Vector G4ThreeVector
const G4ThreeVector & GetAngularMomentum() const
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
G4double plab[VARNTPSIZE]
Momentum.
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)
G4int avv[VARNTPSIZE]
A (-1 for pions).
G4double enerj[VARNTPSIZE]
Kinetic energy.
std::vector< G4ReactionProduct * > G4ReactionProductVector
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.
G4GLOB_DLL std::ostream G4cout
Evaporation and fission output data.
virtual ~G4AblaInterface()
const G4LorentzVector & GetMomentum() const
static G4Triton * Triton()
G4double tetlab[VARNTPSIZE]
Theta angle.
static G4Proton * Proton()
static G4PionPlus * PionPlus()
G4double GetGroundStateMass() const
static G4Neutron * Neutron()
static const G4double A[nN]
static G4PionZero * PionZero()
static G4Deuteron * Deuteron()
G4ReactionProduct * toG4Particle(G4int A, G4int Z, G4double kinE, G4double px, G4double py, G4double pz) const
Convert an Abla particle to a G4ReactionProduct.
static G4IonTable * GetIonTable()
G4int ntrack
Number of particles.
void initEvapora()
Initialize ABLA evaporation code.
G4double energy(const ThreeVector &p, const G4double m)
static G4PionMinus * PionMinus()
void clear()
Clear and initialize all variables and arrays.
G4double philab[VARNTPSIZE]
Phi angle.
G4ParticleDefinition * toG4ParticleDefinition(G4int A, G4int Z) const
Convert A and Z to a G4ParticleDefinition.
G4double GetExcitationEnergy() const
CLHEP::HepLorentzVector G4LorentzVector