Gabriele Cosmo
481 lines
19 KiB
C++
481 lines
19 KiB
C++
//
|
|
// ********************************************************************
|
|
// * License and Disclaimer *
|
|
// * *
|
|
// * The Geant4 software is copyright of the Copyright Holders of *
|
|
// * the Geant4 Collaboration. It is provided under the terms and *
|
|
// * conditions of the Geant4 Software License, included in the file *
|
|
// * LICENSE and available at http://cern.ch/geant4/license . These *
|
|
// * include a list of copyright holders. *
|
|
// * *
|
|
// * Neither the authors of this software system, nor their employing *
|
|
// * institutes,nor the agencies providing financial support for this *
|
|
// * work make any representation or warranty, express or implied, *
|
|
// * regarding this software system or assume any liability for its *
|
|
// * use. Please see the license in the file LICENSE and URL above *
|
|
// * for the full disclaimer and the limitation of liability. *
|
|
// * *
|
|
// * This code implementation is the result of the scientific and *
|
|
// * technical work of the GEANT4 collaboration. *
|
|
// * By using, copying, modifying or distributing the software (or *
|
|
// * any work based on the software) you agree to acknowledge its *
|
|
// * use in resulting scientific publications, and indicate your *
|
|
// * acceptance of all terms of the Geant4 Software license. *
|
|
// ********************************************************************
|
|
//
|
|
// INCL++ intra-nuclear cascade model
|
|
// Alain Boudard, CEA-Saclay, France
|
|
// Joseph Cugnon, University of Liege, Belgium
|
|
// Jean-Christophe David, CEA-Saclay, France
|
|
// Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
|
|
// Sylvie Leray, CEA-Saclay, France
|
|
// Davide Mancusi, CEA-Saclay, France
|
|
//
|
|
#define INCLXX_IN_GEANT4_MODE 1
|
|
|
|
#include "globals.hh"
|
|
|
|
/** \file G4INCLEventInfo.hh
|
|
* \brief Simple container for output of event results.
|
|
*
|
|
* Contains the results of an INCL cascade.
|
|
*
|
|
* \date 21 January 2011
|
|
* \author Davide Mancusi
|
|
*/
|
|
|
|
#ifndef G4INCLEVENTINFO_HH_HH
|
|
#define G4INCLEVENTINFO_HH_HH 1
|
|
|
|
#include "G4INCLParticleType.hh"
|
|
#ifdef INCL_ROOT_USE
|
|
#include <Rtypes.h>
|
|
#endif
|
|
#include <string>
|
|
#include <vector>
|
|
#include <algorithm>
|
|
|
|
namespace G4INCL {
|
|
#ifndef INCL_ROOT_USE
|
|
typedef G4int Int_t;
|
|
typedef short Short_t;
|
|
typedef G4float Float_t;
|
|
typedef G4double Double_t;
|
|
typedef G4bool Bool_t;
|
|
#endif
|
|
|
|
struct EventInfo {
|
|
EventInfo() :
|
|
nParticles(0),
|
|
event(0),
|
|
eventBias((Float_t)0.0),
|
|
nRemnants(0),
|
|
projectileType(0),
|
|
At(0),
|
|
Zt(0),
|
|
St(0),
|
|
Ap(0),
|
|
Zp(0),
|
|
Sp(0),
|
|
Ep((Float_t)0.0),
|
|
impactParameter((Float_t)0.0),
|
|
nCollisions(0),
|
|
stoppingTime((Float_t)0.0),
|
|
EBalance((Float_t)0.0),
|
|
firstEBalance((Float_t)0.0),
|
|
pLongBalance((Float_t)0.0),
|
|
pTransBalance((Float_t)0.0),
|
|
nCascadeParticles(0),
|
|
transparent(false),
|
|
annihilationP(false),
|
|
annihilationN(false),
|
|
forcedCompoundNucleus(false),
|
|
nucleonAbsorption(false),
|
|
pionAbsorption(false),
|
|
nDecays(0),
|
|
fission(false),
|
|
fissmode(0),
|
|
EStarFis((Float_t)0.0),
|
|
ASad(0),
|
|
ZSad(0),
|
|
nSrcCollisions(0),
|
|
nSrcPairs(0),
|
|
nBlockedCollisions(0),
|
|
nBlockedDecays(0),
|
|
effectiveImpactParameter((Float_t)0.0),
|
|
deltasInside(false),
|
|
sigmasInside(false),
|
|
kaonsInside(false),
|
|
antinucleonsInside(false),
|
|
antikaonsInside(false),
|
|
lambdasInside(false),
|
|
forcedDeltasInside(false),
|
|
forcedDeltasOutside(false),
|
|
forcedPionResonancesOutside(false),
|
|
absorbedStrangeParticle(false),
|
|
forcedSigmaOutside(false),
|
|
forcedStrangeInside(false),
|
|
emitLambda(0),
|
|
emitAntilambda(0),
|
|
emitKaon(false),
|
|
emitAntinucleon(false),
|
|
clusterDecay(false),
|
|
firstCollisionTime((Float_t)0.0),
|
|
firstCollisionXSec((Float_t)0.0),
|
|
firstCollisionSpectatorPosition((Float_t)0.0),
|
|
firstCollisionSpectatorMomentum((Float_t)0.0),
|
|
firstCollisionIsElastic(false),
|
|
nReflectionAvatars(0),
|
|
nCollisionAvatars(0),
|
|
nDecayAvatars(0),
|
|
nUnmergedSpectators(0),
|
|
nEnergyViolationInteraction(0)
|
|
|
|
{
|
|
std::fill_n(A, maxSizeParticles, 0);
|
|
std::fill_n(Z, maxSizeParticles, 0);
|
|
std::fill_n(S, maxSizeParticles, 0);
|
|
std::fill_n(J, maxSizeParticles, 0);
|
|
std::fill_n(PDGCode, maxSizeParticles, 0);
|
|
std::fill_n(ParticleBias, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(EKin, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(px, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(py, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(pz, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(theta, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(phi, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(origin, maxSizeParticles, 0);
|
|
std::fill_n(parentResonancePDGCode, maxSizeParticles, 0);
|
|
std::fill_n(parentResonanceID, maxSizeParticles, 0);
|
|
std::fill_n(emissionTime, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(ARem, maxSizeRemnants, 0);
|
|
std::fill_n(ZRem, maxSizeRemnants, 0);
|
|
std::fill_n(SRem, maxSizeRemnants, 0);
|
|
std::fill_n(EStarRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(JRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(EKinRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(pxRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(pyRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(pzRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(thetaRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(phiRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(jxRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(jyRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(jzRem, maxSizeRemnants, (Float_t)0.0);
|
|
std::fill_n(EKinPrime, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(pzPrime, maxSizeParticles, (Float_t)0.0);
|
|
std::fill_n(thetaPrime, maxSizeParticles, (Float_t)0.0);
|
|
}
|
|
|
|
/** \brief Number of the event */
|
|
static G4ThreadLocal Int_t eventNumber;
|
|
|
|
/** \brief Maximum array size for remnants */
|
|
static const Short_t maxSizeRemnants = 10;
|
|
|
|
/** \brief Maximum array size for emitted particles */
|
|
static const Short_t maxSizeParticles = 1000;
|
|
|
|
/** \brief Number of particles in the final state */
|
|
Short_t nParticles;
|
|
/** \brief Sequential number of the event in the event loop */
|
|
Int_t event;
|
|
/** \brief Particle mass number */
|
|
Short_t A[maxSizeParticles];
|
|
/** \brief Particle charge number */
|
|
Short_t Z[maxSizeParticles];
|
|
/** \brief Particle strangeness number */
|
|
Short_t S[maxSizeParticles];
|
|
/** \brief Particle angular momemtum */
|
|
Short_t J[maxSizeParticles];
|
|
/** \brief PDG numbering of the particles */
|
|
Int_t PDGCode[maxSizeParticles];
|
|
/** \brief Event bias */
|
|
Float_t eventBias;
|
|
/** \brief Particle weight due to the bias */
|
|
Float_t ParticleBias[maxSizeParticles];
|
|
/** \brief Particle kinetic energy [MeV] */
|
|
Float_t EKin[maxSizeParticles];
|
|
/** \brief Particle momentum, x component [MeV/c] */
|
|
Float_t px[maxSizeParticles];
|
|
/** \brief Particle momentum, y component [MeV/c] */
|
|
Float_t py[maxSizeParticles];
|
|
/** \brief Particle momentum, z component [MeV/c] */
|
|
Float_t pz[maxSizeParticles];
|
|
/** \brief Particle momentum polar angle [radians] */
|
|
Float_t theta[maxSizeParticles];
|
|
/** \brief Particle momentum azimuthal angle [radians] */
|
|
Float_t phi[maxSizeParticles];
|
|
/** \brief Origin of the particle
|
|
*
|
|
* Should be -1 for cascade particles, or the number of the remnant for
|
|
* de-excitation particles. */
|
|
Short_t origin[maxSizeParticles];
|
|
/** \brief Particle's parent resonance PDG code */
|
|
Int_t parentResonancePDGCode[maxSizeParticles];
|
|
/** \brief Particle's parent resonance unique ID identifier */
|
|
Int_t parentResonanceID[maxSizeParticles];
|
|
/** \brief History of the particle
|
|
*
|
|
* Condensed information about the de-excitation chain of a particle. For
|
|
* cascade particles, it is just an empty string. For particles arising
|
|
* from the de-excitation of a cascade remnant, it is a string of
|
|
* characters. Each character represents one or more identical steps in
|
|
* the de-excitation process. The currently defined possible character
|
|
* values and their meanings are the following:
|
|
*
|
|
* e: evaporation product
|
|
* E: evaporation residue
|
|
* m: multifragmentation
|
|
* a: light partner in asymmetric fission or IMF emission
|
|
* A: heavy partner in asymmetric fission or IMF emission
|
|
* f: light partner in fission
|
|
* F: heavy partner in fission
|
|
* s: saddle-to-scission emission
|
|
* n: non-statistical emission (decay) */
|
|
std::vector<std::string> history;
|
|
/** \brief Number of remnants */
|
|
Short_t nRemnants;
|
|
/** \brief Projectile particle type */
|
|
Int_t projectileType;
|
|
/** \brief Mass number of the target nucleus */
|
|
Short_t At;
|
|
/** \brief Charge number of the target nucleus */
|
|
Short_t Zt;
|
|
/** \brief Strangeness number of the target nucleus */
|
|
Short_t St;
|
|
/** \brief Mass number of the projectile nucleus */
|
|
Short_t Ap;
|
|
/** \brief Charge number of the projectile nucleus */
|
|
Short_t Zp;
|
|
/** \brief Strangeness number of the projectile nucleus */
|
|
Short_t Sp;
|
|
/** \brief Projectile kinetic energy given as input */
|
|
Float_t Ep;
|
|
/** \brief Impact parameter [fm] */
|
|
Float_t impactParameter;
|
|
/** \brief Number of accepted two-body collisions */
|
|
Int_t nCollisions;
|
|
/** \brief Cascade stopping time [fm/c] */
|
|
Float_t stoppingTime;
|
|
/** \brief Energy-conservation balance [MeV] */
|
|
Float_t EBalance;
|
|
/** \brief First value for the energy-conservation balance [MeV] */
|
|
Float_t firstEBalance;
|
|
/** \brief Longitudinal momentum-conservation balance [MeV/c] */
|
|
Float_t pLongBalance;
|
|
/** \brief Transverse momentum-conservation balance [MeV/c] */
|
|
Float_t pTransBalance;
|
|
/** \brief Number of cascade particles */
|
|
Short_t nCascadeParticles;
|
|
/** \brief True if the event is transparent */
|
|
Bool_t transparent;
|
|
/** \brief True if annihilation at rest on a proton */
|
|
Bool_t annihilationP;
|
|
/** \brief True if annihilation at rest on a neutron */
|
|
Bool_t annihilationN;
|
|
/** \brief True if the event is a forced CN */
|
|
Bool_t forcedCompoundNucleus;
|
|
/** \brief True if the event is a nucleon absorption */
|
|
Bool_t nucleonAbsorption;
|
|
/** \brief True if the event is a pion absorption */
|
|
Bool_t pionAbsorption;
|
|
/** \brief Number of accepted Delta decays */
|
|
Int_t nDecays;
|
|
/** \brief True if the event is fission */
|
|
Bool_t fission;
|
|
/** \brief Fission mode */
|
|
Short_t fissmode;
|
|
/** \brief Excitation energy above fission barrier [MeV] */
|
|
Float_t EStarFis;
|
|
/** \brief Mass number at saddle */
|
|
Short_t ASad;
|
|
/** \brief Charge number at saddle */
|
|
Short_t ZSad;
|
|
/** \brief Mass number at scission */
|
|
std::vector<Int_t> ASci;
|
|
/** \brief Charge number at scission */
|
|
std::vector<Int_t> ZSci;
|
|
/** \brief Number of accepted SRC collisions */
|
|
Int_t nSrcCollisions;
|
|
/** \brief Number of src pairs */
|
|
Int_t nSrcPairs;
|
|
/** \brief Number of two-body collisions blocked by Pauli or CDPP */
|
|
Int_t nBlockedCollisions;
|
|
/** \brief Number of decays blocked by Pauli or CDPP */
|
|
Int_t nBlockedDecays;
|
|
/** \brief Effective (Coulomb-distorted) impact parameter [fm] */
|
|
Float_t effectiveImpactParameter;
|
|
/** \brief Event involved deltas in the nucleus at the end of the cascade */
|
|
Bool_t deltasInside;
|
|
/** \brief Event involved sigmas in the nucleus at the end of the cascade */
|
|
Bool_t sigmasInside;
|
|
/** \brief Event involved kaons in the nucleus at the end of the cascade */
|
|
Bool_t kaonsInside;
|
|
/** \brief Event involved antinucleons in the nucleus at the end of the cascade */
|
|
Bool_t antinucleonsInside;
|
|
/** \brief Event involved antikaons in the nucleus at the end of the cascade */
|
|
Bool_t antikaonsInside;
|
|
/** \brief Event involved lambdas in the nucleus at the end of the cascade */
|
|
Bool_t lambdasInside;
|
|
/** \brief Event involved forced delta decays inside the nucleus */
|
|
Bool_t forcedDeltasInside;
|
|
/** \brief Event involved forced delta decays outside the nucleus */
|
|
Bool_t forcedDeltasOutside;
|
|
/** \brief Event involved forced eta/omega decays outside the nucleus */
|
|
Bool_t forcedPionResonancesOutside;
|
|
/** \brief Event involved forced strange absorption inside the nucleus */
|
|
Bool_t absorbedStrangeParticle;
|
|
/** \brief Event involved forced Sigma Zero decays outside the nucleus */
|
|
Bool_t forcedSigmaOutside;
|
|
/** \brief Event involved forced antiKaon/Sigma absorption inside the nucleus */
|
|
Bool_t forcedStrangeInside;
|
|
/** \brief Number of forced Lambda emit out of the nucleus */
|
|
Int_t emitLambda;
|
|
/** \brief Number of forced Antilambda emit out of the nucleus */
|
|
Int_t emitAntilambda;
|
|
/** \brief Event involved forced Kaon emission */
|
|
Bool_t emitKaon;
|
|
/** \brief Event involved forced Antinucleon emission */
|
|
Bool_t emitAntinucleon;
|
|
/** \brief Event involved cluster decay */
|
|
Bool_t clusterDecay;
|
|
/** \brief Time of the first collision [fm/c] */
|
|
Float_t firstCollisionTime;
|
|
/** \brief Cross section of the first collision (mb) */
|
|
Float_t firstCollisionXSec;
|
|
/** \brief Position of the spectator on the first collision (fm) */
|
|
Float_t firstCollisionSpectatorPosition;
|
|
/** \brief Momentum of the spectator on the first collision (fm) */
|
|
Float_t firstCollisionSpectatorMomentum;
|
|
/** \brief True if the first collision was elastic */
|
|
Bool_t firstCollisionIsElastic;
|
|
/** \brief Number of reflection avatars */
|
|
Int_t nReflectionAvatars;
|
|
/** \brief Number of collision avatars */
|
|
Int_t nCollisionAvatars;
|
|
/** \brief Number of decay avatars */
|
|
Int_t nDecayAvatars;
|
|
/** \brief Number of dynamical spectators that were merged back into the projectile remnant */
|
|
Int_t nUnmergedSpectators;
|
|
/** \brief Number of attempted collisions/decays for which the energy-conservation algorithm failed to find a solution. */
|
|
Int_t nEnergyViolationInteraction;
|
|
/** \brief Emission time [fm/c] */
|
|
Float_t emissionTime[maxSizeParticles];
|
|
/** \brief Remnant mass number */
|
|
Short_t ARem[maxSizeRemnants];
|
|
/** \brief Remnant charge number */
|
|
Short_t ZRem[maxSizeRemnants];
|
|
/** \brief Remnant strangeness number */
|
|
Short_t SRem[maxSizeRemnants];
|
|
/** \brief Remnant excitation energy [MeV] */
|
|
Float_t EStarRem[maxSizeRemnants];
|
|
/** \brief Remnant spin [\f$\hbar\f$] */
|
|
Float_t JRem[maxSizeRemnants];
|
|
/** \brief Remnant kinetic energy [MeV] */
|
|
Float_t EKinRem[maxSizeRemnants];
|
|
/** \brief Remnant momentum, x component [MeV/c] */
|
|
Float_t pxRem[maxSizeRemnants];
|
|
/** \brief Remnant momentum, y component [MeV/c] */
|
|
Float_t pyRem[maxSizeRemnants];
|
|
/** \brief Remnant momentum, z component [MeV/c] */
|
|
Float_t pzRem[maxSizeRemnants];
|
|
/** \brief Remnant momentum polar angle [radians] */
|
|
Float_t thetaRem[maxSizeRemnants];
|
|
/** \brief Remnant momentum azimuthal angle [radians] */
|
|
Float_t phiRem[maxSizeRemnants];
|
|
/** \brief Remnant angular momentum, x component [\f$\hbar\f$] */
|
|
Float_t jxRem[maxSizeRemnants];
|
|
/** \brief Remnant angular momentum, y component [\f$\hbar\f$] */
|
|
Float_t jyRem[maxSizeRemnants];
|
|
/** \brief Remnant angular momentum, z component [\f$\hbar\f$] */
|
|
Float_t jzRem[maxSizeRemnants];
|
|
/** \brief Particle kinetic energy, in inverse kinematics [MeV] */
|
|
Float_t EKinPrime[maxSizeParticles];
|
|
/** \brief Particle momentum, z component, in inverse kinematics [MeV/c] */
|
|
Float_t pzPrime[maxSizeParticles];
|
|
/** \brief Particle momentum polar angle, in inverse kinematics [radians] */
|
|
Float_t thetaPrime[maxSizeParticles];
|
|
|
|
/** \brief Reset the EventInfo members */
|
|
void reset() {
|
|
nParticles = 0;
|
|
event = 0;
|
|
eventBias = (Float_t)0.0;
|
|
history.clear();
|
|
nRemnants = 0;
|
|
projectileType = 0;
|
|
At = 0;
|
|
Zt = 0;
|
|
St = 0;
|
|
Ap = 0;
|
|
Zp = 0;
|
|
Sp = 0;
|
|
Ep = (Float_t)0.0;
|
|
impactParameter = (Float_t)0.0;
|
|
nCollisions = 0;
|
|
stoppingTime = (Float_t)0.0;
|
|
EBalance = (Float_t)0.0;
|
|
firstEBalance = (Float_t)0.0;
|
|
pLongBalance = (Float_t)0.0;
|
|
pTransBalance = (Float_t)0.0;
|
|
nCascadeParticles = 0;
|
|
transparent = false;
|
|
annihilationP = false;
|
|
annihilationN = false;
|
|
forcedCompoundNucleus = false;
|
|
nucleonAbsorption = false;
|
|
pionAbsorption = false;
|
|
nDecays = 0;
|
|
fission = false;
|
|
fissmode = 0;
|
|
EStarFis = (Float_t)0.0;
|
|
ASad = 0;
|
|
ZSad = 0;
|
|
ASci.clear();
|
|
ZSci.clear();
|
|
nSrcCollisions = 0;
|
|
nSrcPairs = 0;
|
|
nBlockedCollisions = 0;
|
|
nBlockedDecays = 0;
|
|
effectiveImpactParameter = (Float_t)0.0;
|
|
deltasInside = false;
|
|
sigmasInside = false;
|
|
kaonsInside = false;
|
|
antinucleonsInside = false;
|
|
antikaonsInside = false;
|
|
lambdasInside = false;
|
|
forcedDeltasInside = false;
|
|
forcedDeltasOutside = false;
|
|
forcedPionResonancesOutside = false;
|
|
absorbedStrangeParticle = false;
|
|
forcedSigmaOutside = false;
|
|
forcedStrangeInside = false;
|
|
emitLambda = 0;
|
|
emitAntilambda = 0;
|
|
emitKaon = false;
|
|
emitAntinucleon = false;
|
|
clusterDecay = false;
|
|
firstCollisionTime = (Float_t)0.0;
|
|
firstCollisionXSec = (Float_t)0.0;
|
|
firstCollisionSpectatorPosition = (Float_t)0.0;
|
|
firstCollisionSpectatorMomentum = (Float_t)0.0;
|
|
firstCollisionIsElastic = false;
|
|
nReflectionAvatars = 0;
|
|
nCollisionAvatars = 0;
|
|
nDecayAvatars = 0;
|
|
nUnmergedSpectators = 0;
|
|
nEnergyViolationInteraction = 0;
|
|
|
|
}
|
|
|
|
/// \brief Move a remnant to the particle array
|
|
void remnantToParticle(const G4int remnantIndex);
|
|
|
|
/// \brief Fill the variables describing the reaction in inverse kinematics
|
|
void fillInverseKinematics(const Double_t gamma);
|
|
};
|
|
}
|
|
|
|
#endif /* G4INCLEVENTINFO_HH_HH */
|