FullModelReactionDynamics.h

Go to the documentation of this file.

//
// Hadronic Process: Reaction Dynamics
// original by H.P. Wellisch
// Modified by J.L.Chuma 19-Nov-96
// Modified by J.L.Chuma 27-Mar-97
// Modified by J.L.Chuma 30-Apr-97
// Modified by J.L.Chuma 06-Aug-97  to include the original incident particle
//                                  before Fermi motion and evaporation effects

#ifndef SimG4Core_CustomPhysics_FullModelReactionDynamics_h
#define SimG4Core_CustomPhysics_FullModelReactionDynamics_h 1

#include "G4ParticleTypes.hh"
#include "G4DynamicParticle.hh"
#include "G4ReactionProduct.hh"
#include "G4Nucleus.hh"
#include "G4FastVector.hh"
#include "G4HadProjectile.hh"

enum { MYGHADLISTSIZE = 256 };

class FullModelReactionDynamics {
public:
  FullModelReactionDynamics() {}

  ~FullModelReactionDynamics() = default;

  G4bool GenerateXandPt(  // derived from GENXPT
      G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
      G4int &vecLen,
      G4ReactionProduct &modifiedOriginal,  // Fermi motion & evap. effect included
      const G4HadProjectile *originalIncident,
      G4ReactionProduct &currentParticle,
      G4ReactionProduct &targetParticle,
      const G4Nucleus &targetNucleus,
      G4bool &incidentHasChanged,
      G4bool &targetHasChanged,
      G4bool leadFlag,
      G4ReactionProduct &leadingStrangeParticle);

  void SuppressChargedPions(G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
                            G4int &vecLen,
                            const G4ReactionProduct &modifiedOriginal,
                            G4ReactionProduct &currentParticle,
                            G4ReactionProduct &targetParticle,
                            const G4Nucleus &targetNucleus,
                            G4bool &incidentHasChanged,
                            G4bool &targetHasChanged);

  G4bool TwoCluster(  // derived from TWOCLU
      G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
      G4int &vecLen,
      G4ReactionProduct &modifiedOriginal,  // Fermi motion & evap. effect included
      const G4HadProjectile *originalIncident,
      G4ReactionProduct &currentParticle,
      G4ReactionProduct &targetParticle,
      const G4Nucleus &targetNucleus,
      G4bool &incidentHasChanged,
      G4bool &targetHasChanged,
      G4bool leadFlag,
      G4ReactionProduct &leadingStrangeParticle);

  void TwoBody(  // derived from TWOB
      G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
      G4int &vecLen,
      G4ReactionProduct &modifiedOriginal,
      const G4DynamicParticle *originalTarget,
      G4ReactionProduct &currentParticle,
      G4ReactionProduct &targetParticle,
      const G4Nucleus &targetNucleus,
      G4bool &targetHasChanged);

  G4int Factorial(G4int n);

  G4double GenerateNBodyEvent(  // derived from PHASP
      const G4double totalEnergy,
      const G4bool constantCrossSection,
      G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
      G4int &vecLen);

  void ProduceStrangeParticlePairs(G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
                                   G4int &vecLen,
                                   const G4ReactionProduct &modifiedOriginal,
                                   const G4DynamicParticle *originalTarget,
                                   G4ReactionProduct &currentParticle,
                                   G4ReactionProduct &targetParticle,
                                   G4bool &incidentHasChanged,
                                   G4bool &targetHasChanged);

  void NuclearReaction(  // derived from NUCREC
      G4FastVector<G4ReactionProduct, 4> &vec,
      G4int &vecLen,
      const G4HadProjectile *originalIncident,
      const G4Nucleus &aNucleus,
      const G4double theAtomicMass,
      const G4double *massVec);

  FullModelReactionDynamics(const FullModelReactionDynamics &) = delete;
  FullModelReactionDynamics &operator=(const FullModelReactionDynamics &) = delete;

private:
  void Rotate(const G4double numberofFinalStateNucleons,
              const G4ThreeVector &temp,
              const G4ReactionProduct &modifiedOriginal,  // Fermi motion & evap. effect included
              const G4HadProjectile *originalIncident,
              const G4Nucleus &targetNucleus,
              G4ReactionProduct &currentParticle,
              G4ReactionProduct &targetParticle,
              G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
              G4int &vecLen);

  void Defs1(const G4ReactionProduct &modifiedOriginal,
             G4ReactionProduct &currentParticle,
             G4ReactionProduct &targetParticle,
             G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
             G4int &vecLen);

  void AddBlackTrackParticles(const G4double epnb,
                              const G4int npnb,
                              const G4double edta,
                              const G4int ndta,
                              const G4double sprob,
                              const G4double kineticMinimum,
                              const G4double kineticFactor,
                              const G4ReactionProduct &modifiedOriginal,
                              G4double spall,
                              const G4Nucleus &aNucleus,
                              G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
                              G4int &vecLen);

  void MomentumCheck(const G4ReactionProduct &modifiedOriginal,
                     G4ReactionProduct &currentParticle,
                     G4ReactionProduct &targetParticle,
                     G4FastVector<G4ReactionProduct, MYGHADLISTSIZE> &vec,
                     G4int &vecLen);

  G4double normal();
};

#endif