ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4INCLCrossSectionsINCL46.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4INCLCrossSectionsINCL46.hh
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 
45 #ifndef G4INCLCROSSSECTIONSINCL46_HH
46 #define G4INCLCROSSSECTIONSINCL46_HH
47 
48 #include "G4INCLICrossSections.hh"
49 
50 namespace G4INCL {
52 
53  class CrossSectionsINCL46 : public ICrossSections{
54  public:
55 
57  virtual G4double elastic(Particle const * const p1, Particle const * const p2);
58 
60  virtual G4double total(Particle const * const p1, Particle const * const p2);
61 
63  virtual G4double NDeltaToNN(Particle const * const p1, Particle const * const p2);
64 
66  virtual G4double NNToNDelta(Particle const * const p1, Particle const * const p2);
67 
69  virtual G4double piNToDelta(Particle const * const p1, Particle const * const p2);
70 
72  virtual G4double NNToxPiNN(const G4int xpi, Particle const * const p1, Particle const * const p2);
73 
75  virtual G4double piNToxPiN(const G4int xpi, Particle const * const p1, Particle const * const p2);
76 
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 
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 
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 
92  virtual G4double NNToNNEta(Particle const * const p1, Particle const * const p2);
93 
95  virtual G4double NNToNNEtaExclu(Particle const * const p1, Particle const * const p2);
96 
98  virtual G4double NNToNNEtaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
99 
101  virtual G4double NNToNDeltaEta(Particle const * const p1, Particle const * const p2);
102 
103 
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 
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 
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 
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 
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 
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 
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 
160  virtual G4double NNToNNOmega(Particle const * const particle1, Particle const * const particle2);
161 
163  virtual G4double NNToNNOmegaExclu(Particle const * const particle1, Particle const * const particle2);
164 
166  virtual G4double NNToNNOmegaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
167 
169  virtual G4double NNToNDeltaOmega(Particle const * const p1, Particle const * const p2);
170 
178  virtual G4double calculateNNAngularSlope(G4double energyCM, G4int iso);
179 
180  protected:
181 
183  G4double elasticNNLegacy(Particle const * const part1, Particle const * const part2);
184 
186  G4double deltaProduction(const G4int isospin, const G4double pLab);
187 
188  G4double spnPiPlusPHE(const G4double x);
189  G4double spnPiMinusPHE(const G4double x);
190 
191  };
192 }
193 
194 #endif