Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/inclxx/incl_physics/include/G4INCLCrossSectionsINCL46.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 // INCL++ intra-nuclear cascade model
 27 // Alain Boudard, CEA-Saclay, France
 28 // Joseph Cugnon, University of Liege, Belgium
 29 // Jean-Christophe David, CEA-Saclay, France
 30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
 31 // Sylvie Leray, CEA-Saclay, France
 32 // Davide Mancusi, CEA-Saclay, France
 33 //
 34 #define INCLXX_IN_GEANT4_MODE 1
 35 
 36 #include "globals.hh"
 37 
 38 /** \file G4INCLCrossSectionsINCL46.hh
 39  * \brief Cross sections used in INCL4.6
 40  *
 41  * \date 25nd October 2013
 42  * \author Davide Mancusi
 43  */
 44 
 45 #ifndef G4INCLCROSSSECTIONSINCL46_HH
 46 #define G4INCLCROSSSECTIONSINCL46_HH
 47 
 48 #include "G4INCLICrossSections.hh"
 49 
 50 namespace G4INCL {
 51   /// \brief Cross sections used in INCL4.6
 52 
 53   class CrossSectionsINCL46 : public ICrossSections{
 54     public:
 55 
 56       /// \brief Elastic particle-particle cross section
 57       virtual G4double elastic(Particle const * const p1, Particle const * const p2);
 58 
 59       /// \brief Total (elastic+inelastic) particle-particle cross section
 60       virtual G4double total(Particle const * const p1, Particle const * const p2);
 61 
 62       /// \brief Cross section for NDelta->NN
 63       virtual G4double NDeltaToNN(Particle const * const p1, Particle const * const p2);
 64 
 65       /// \brief Cross section for NN->NDelta
 66       virtual G4double NNToNDelta(Particle const * const p1, Particle const * const p2);
 67 
 68       /// \brief Cross section for piN->NDelta
 69       virtual G4double piNToDelta(Particle const * const p1, Particle const * const p2);
 70 
 71       /// \brief Cross section for X pion production - NN Channel
 72       virtual G4double NNToxPiNN(const G4int xpi, Particle const * const p1, Particle const * const p2);
 73 
 74       /// \brief Cross section for X pion production - piN Channel
 75       virtual G4double piNToxPiN(const G4int xpi, Particle const * const p1, Particle const * const p2);
 76     
 77       /// \brief Cross sections for mesonic resonance production - piN Channel
 78       virtual G4double piNToEtaN(Particle const * const p1, Particle const * const p2);
 79       virtual G4double piNToOmegaN(Particle const * const p1, Particle const * const p2);
 80       virtual G4double piNToEtaPrimeN(Particle const * const p1, Particle const * const p2);
 81     
 82       /// \brief Cross sections for mesonic resonance absorption on nucleons - piN Channel
 83       virtual G4double etaNToPiN(Particle const * const p1, Particle const * const p2);
 84       virtual G4double omegaNToPiN(Particle const * const p1, Particle const * const p2);
 85       virtual G4double etaPrimeNToPiN(Particle const * const p1, Particle const * const p2);
 86       
 87       /// \brief Cross sections for mesonic resonance absorption on nucleon - pipiN Channel
 88       virtual G4double etaNToPiPiN(Particle const * const p1, Particle const * const p2);     
 89       virtual G4double omegaNToPiPiN(Particle const * const p1, Particle const * const p2);     
 90       
 91       /// \brief Cross section for Eta production - NN entrance channel
 92       virtual G4double NNToNNEta(Particle const * const p1, Particle const * const p2);
 93       
 94       /// \brief Cross section for Eta production  (exclusive) - NN entrance channel
 95       virtual G4double NNToNNEtaExclu(Particle const * const p1, Particle const * const p2);
 96     
 97       /// \brief Cross section for X pion production - NNEta Channel
 98       virtual G4double NNToNNEtaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
 99 
100       /// \brief Cross section for N-Delta-Eta production - NNEta Channel
101       virtual G4double NNToNDeltaEta(Particle const * const p1, Particle const * const p2);
102       
103       
104       /// \brief elastic scattering for Nucleon-Strange Particles cross sections
105       virtual G4double NYelastic(Particle const * const p1, Particle const * const p2);
106       virtual G4double NKbelastic(Particle const * const p1, Particle const * const p2);
107       virtual G4double NKelastic(Particle const * const p1, Particle const * const p2);
108       
109       /// \brief Nucleon-Nucleon to Stange particles cross sections
110       virtual G4double NNToNLK(Particle const * const p1, Particle const * const p2);
111       virtual G4double NNToNSK(Particle const * const p1, Particle const * const p2);
112       virtual G4double NNToNLKpi(Particle const * const p1, Particle const * const p2);
113       virtual G4double NNToNSKpi(Particle const * const p1, Particle const * const p2);
114       virtual G4double NNToNLK2pi(Particle const * const p1, Particle const * const p2);
115       virtual G4double NNToNSK2pi(Particle const * const p1, Particle const * const p2);
116       virtual G4double NNToNNKKb(Particle const * const p1, Particle const * const p2);
117       virtual G4double NNToMissingStrangeness(Particle const * const p1, Particle const * const p2);
118       
119       /// \brief Nucleon-Delta to Stange particles cross sections
120       virtual G4double NDeltaToNLK(Particle const * const p1, Particle const * const p2);
121       virtual G4double NDeltaToNSK(Particle const * const p1, Particle const * const p2);
122       virtual G4double NDeltaToDeltaLK(Particle const * const p1, Particle const * const p2);
123       virtual G4double NDeltaToDeltaSK(Particle const * const p1, Particle const * const p2);
124       
125       virtual G4double NDeltaToNNKKb(Particle const * const p1, Particle const * const p2);
126       
127       /// \brief Nucleon-Pion to Stange particles cross sections
128       virtual G4double NpiToLK(Particle const * const p1, Particle const * const p2);
129       virtual G4double NpiToSK(Particle const * const p1, Particle const * const p2);
130       virtual G4double p_pimToSzKz(Particle const * const p1, Particle const * const p2);
131       virtual G4double p_pimToSmKp(Particle const * const p1, Particle const * const p2);
132       virtual G4double p_pizToSzKp(Particle const * const p1, Particle const * const p2);
133       virtual G4double NpiToLKpi(Particle const * const p1, Particle const * const p2);
134       virtual G4double NpiToSKpi(Particle const * const p1, Particle const * const p2);
135       virtual G4double NpiToLK2pi(Particle const * const p1, Particle const * const p2);
136       virtual G4double NpiToSK2pi(Particle const * const p1, Particle const * const p2);
137       virtual G4double NpiToNKKb(Particle const * const p1, Particle const * const p2);
138       virtual G4double NpiToMissingStrangeness(Particle const * const p1, Particle const * const p2);
139       
140       /// \brief Nucleon-Hyperon cross sections
141       virtual G4double NLToNS(Particle const * const p1, Particle const * const p2);
142       virtual G4double NSToNL(Particle const * const p1, Particle const * const p2);
143       virtual G4double NSToNS(Particle const * const p1, Particle const * const p2);
144       
145       /// \brief Nucleon-Kaon inelastic cross sections
146       virtual G4double NKToNK(Particle const * const p1, Particle const * const p2);
147       virtual G4double NKToNKpi(Particle const * const p1, Particle const * const p2);
148       virtual G4double NKToNK2pi(Particle const * const p1, Particle const * const p2);
149       
150       /// \brief Nucleon-antiKaon inelastic cross sections
151       virtual G4double NKbToNKb(Particle const * const p1, Particle const * const p2);
152       virtual G4double NKbToSpi(Particle const * const p1, Particle const * const p2);
153       virtual G4double NKbToLpi(Particle const * const p1, Particle const * const p2);
154       virtual G4double NKbToS2pi(Particle const * const p1, Particle const * const p2);
155       virtual G4double NKbToL2pi(Particle const * const p1, Particle const * const p2);
156       virtual G4double NKbToNKbpi(Particle const * const p1, Particle const * const p2);
157       virtual G4double NKbToNKb2pi(Particle const * const p1, Particle const * const p2);
158     
159       /// \brief Cross section for Eta production - NN entrance channel
160       virtual G4double NNToNNOmega(Particle const * const particle1, Particle const * const particle2);
161       
162       /// \brief Cross section for Eta production  (exclusive) - NN entrance channel
163       virtual G4double NNToNNOmegaExclu(Particle const * const particle1, Particle const * const particle2);
164     
165       /// \brief Cross section for X pion production - NNEta Channel
166       virtual G4double NNToNNOmegaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
167     
168       /// \brief Cross section for N-Delta-Eta production - NNEta Channel
169       virtual G4double NNToNDeltaOmega(Particle const * const p1, Particle const * const p2);
170       
171       /// \brief antiparticle cross sections
172       /// \brief Nucleon-AntiNucleon to Baryon-AntiBaryon cross sections
173       virtual G4double NNbarElastic(Particle const* const p1, Particle const* const p2);
174       virtual G4double NNbarCEX(Particle const* const p1, Particle const* const p2);
175 
176       virtual G4double NNbarToLLbar(Particle const * const p1, Particle const * const p2);
177       
178       /// \brief Nucleon-AntiNucleon to Nucleon-AntiNucleon + pions cross sections
179       virtual G4double NNbarToNNbarpi(Particle const* const p1, Particle const* const p2);
180       virtual G4double NNbarToNNbar2pi(Particle const* const p1, Particle const* const p2);
181       virtual G4double NNbarToNNbar3pi(Particle const* const p1, Particle const* const p2);
182      
183       /// \brief Nucleon-AntiNucleon total annihilation cross sections
184       virtual G4double NNbarToAnnihilation(Particle const* const p1, Particle const* const p2);     
185     
186       /** \brief Calculate the slope of the NN DDXS.
187        *
188        * \param energyCM energy in the CM frame, in MeV
189        * \param iso total isospin of the system
190        *
191        * \return the slope of the angular distribution
192        */
193       virtual G4double calculateNNAngularSlope(G4double energyCM, G4int iso);
194 
195     protected:
196 
197       /// \brief Internal implementation of the elastic cross section
198       G4double elasticNNLegacy(Particle const * const part1, Particle const * const part2);
199 
200       /// \brief Internal function for the delta-production cross section
201       G4double deltaProduction(const G4int isospin, const G4double pLab);
202 
203       G4double spnPiPlusPHE(const G4double x);
204       G4double spnPiMinusPHE(const G4double x);
205     
206   };
207 }
208 
209 #endif
210