HAIntranuke2025.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::Intranuke
 
\brief    The INTRANUKE intranuclear hadron transport MC.
          Is a concrete implementation of the EventRecordVisitorI interface.
 
\ref      R.Merenyi et al., Phys.Rev.D45 (1992)
          R.D.Ransome, Nucl.Phys.B 139 (2005)
 
          Current INTRANUKE development is led by S.Dytman and H.Gallagher.
          The original INTRANUKE cascade MC was developed (in fortran) for the
          NeuGEN MC by R.Edgecock, G.F.Pearce, W.A.Mann, R.Merenyi and others.
 
\author   Steve Dytman <dytman+@pitt.edu>, Pittsburgh University
          Aaron Meyer <asm58@pitt.edu>, Pittsburgh University
          Alex Bell, Pittsburgh University
          Hugh Gallagher <gallag@minos.phy.tufts.edu>, Tufts University
          Costas Andreopoulos <c.andreopoulos \at cern.ch> STFC, Rutherford Lab
 
\created  September 20, 2005
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
 
\created on  september,2025 by Mohamed Ismail <msi10@pitt.edu>
no changes from the 2018 class version, major changes made in INukeHadroData2025 for new hA pion splines
Add data, use hN for high pion KE, use INCL for low energy.  Use splines for channel and total reac xs
to improve accuracy.  Also, smooth results to avoid discontinuities.
*/
//____________________________________________________________________________
 
#ifndef _HA_INTRANUKE_2025_H_
#define _HA_INTRANUKE_2025_H_
 
#include <TGenPhaseSpace.h>
 
#include "Physics/NuclearState/NuclearModelI.h"
#include "Framework/Algorithm/AlgFactory.h"
#include "Framework/EventGen/EventRecordVisitorI.h"
#include "Physics/HadronTransport/INukeMode.h"
#include "Physics/HadronTransport/INukeHadroFates2025.h"
#include "Physics/HadronTransport/Intranuke2025.h"
 
class TLorentzVector;
class TVector3;
 
namespace genie {
 
class GHepParticle;
class INukeHadroData2025;
class PDGCodeList;
 
class HAIntranuke2025 : public Intranuke2025 {
 
friend class IntranukeTester;
 
public :
  HAIntranuke2025();
  HAIntranuke2025(string config);
 ~HAIntranuke2025();
 
  void ProcessEventRecord(GHepRecord * event_rec) const;
 
  virtual string GetINukeMode() const {return "hA2025";};
  virtual string GetGenINukeMode() const {return "hA";};
 
private:
 
  void LoadConfig (void);
 
  void  SimulateHadronicFinalState           (GHepRecord* ev, GHepParticle* p) const;
  void  SimulateHadronicFinalStateKinematics (GHepRecord* ev, GHepParticle* p) const;
 
  INukeFateHA_t HadronFateHA     (const GHepParticle* p) const;
  //INukeFateHA_t HadronFateOset   (void) const;
  void          Inelastic        (GHepRecord* ev, GHepParticle* p, INukeFateHA_t fate) const;
  void          ElasHA           (GHepRecord* ev, GHepParticle* p, INukeFateHA_t fate) const;
  void          InelasticHA      (GHepRecord* ev, GHepParticle* p, INukeFateHA_t fate) const;
  double        PiBounce         (void) const;
  double        PnBounce         (void) const;
  int           HandleCompoundNucleus(GHepRecord* ev, GHepParticle* p, int mom) const;           
 
  mutable int nuclA;     ///< value of A for the target nucleus in hA mode
  mutable unsigned int fNumIterations;
};
 
}      // genie namespace
 
#endif // _HA_INTRANUKE_2025_H_