ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4INCLCrossSectionsMultiPionsAndResonances.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4INCLCrossSectionsMultiPionsAndResonances.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 G4INCLCROSSSECTIONSMULTIPIONSANDRESONANCES_HH
46 #define G4INCLCROSSSECTIONSMULTIPIONSANDRESONANCES_HH
47 
48 #include "G4INCLCrossSectionsMultiPions.hh"
49 //#include <limits>
50 
51 namespace G4INCL {
53 
54  class CrossSectionsMultiPionsAndResonances : public CrossSectionsMultiPions {
55  public:
56  CrossSectionsMultiPionsAndResonances();
57 
59  virtual G4double elastic(Particle const * const p1, Particle const * const p2);
60 
62  virtual G4double total(Particle const * const p1, Particle const * const p2);
63 
65  virtual G4double piNToxPiN(const G4int xpi, Particle const * const p1, Particle const * const p2);
66 
68  virtual G4double piNToEtaN(Particle const * const p1, Particle const * const p2);
69  virtual G4double piNToOmegaN(Particle const * const p1, Particle const * const p2);
70  virtual G4double piNToEtaPrimeN(Particle const * const p1, Particle const * const p2);
71 
73  virtual G4double etaNToPiN(Particle const * const p1, Particle const * const p2);
74  virtual G4double omegaNToPiN(Particle const * const p1, Particle const * const p2);
75  virtual G4double etaPrimeNToPiN(Particle const * const p1, Particle const * const p2);
76 
78  virtual G4double etaNToPiPiN(Particle const * const p1, Particle const * const p2);
79 
81  virtual G4double NNToNNEta(Particle const * const particle1, Particle const * const particle2);
82 
84  virtual G4double NNToNNEtaExclu(Particle const * const particle1, Particle const * const particle2);
85 
87  virtual G4double NNToNNOmega(Particle const * const particle1, Particle const * const particle2);
88 
90  virtual G4double NNToNNOmegaExclu(Particle const * const particle1, Particle const * const particle2);
91 
93  virtual G4double NNToxPiNN(const G4int xpi, Particle const * const p1, Particle const * const p2);
94 
96  virtual G4double NNToNNEtaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
97 
99  virtual G4double NNToNDeltaEta(Particle const * const p1, Particle const * const p2);
100 
102  virtual G4double NNToNNOmegaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
103 
105  virtual G4double NNToNDeltaOmega(Particle const * const p1, Particle const * const p2);
106 
107 
108  protected:
110  static const G4int nMaxPiNN;
111 
113  static const G4int nMaxPiPiN;
114 
116  const HornerC7 s11pzHC;
118  const HornerC8 s01ppHC;
120  const HornerC4 s01pzHC;
122  const HornerC4 s11pmHC;
124  const HornerC5 s12pmHC;
126  const HornerC3 s12ppHC;
128  const HornerC4 s12zzHC;
130  const HornerC4 s02pzHC;
132  const HornerC6 s02pmHC;
134  const HornerC4 s12mzHC;
135 
137  static const G4double s11pzOOT;
139  static const G4double s01ppOOT;
141  static const G4double s01pzOOT;
143  static const G4double s11pmOOT;
145  static const G4double s12pmOOT;
147  static const G4double s12ppOOT;
149  static const G4double s12zzOOT;
151  static const G4double s02pzOOT;
153  static const G4double s02pmOOT;
155  static const G4double s12mzOOT;
156 
157 
159  G4double piMinuspToEtaN(Particle const * const p1, Particle const * const p2);
160  G4double piMinuspToEtaN(const G4double ECM);
161  G4double piMinuspToOmegaN(Particle const * const p1, Particle const * const p2);
162  G4double piMinuspToOmegaN(const G4double ECM);
163 // G4double piPluspOnePi(Particle const * const p1, Particle const * const p2);
164 // G4double piMinuspOnePi(Particle const * const p1, Particle const * const p2);
165 // G4double piPluspTwoPi(Particle const * const p1, Particle const * const p2);
166 // G4double piMinuspTwoPi(Particle const * const p1, Particle const * const p2);
167 
169 // virtual G4double piNOnePi(Particle const * const p1, Particle const * const p2);
170 
172 // virtual G4double piNTwoPi(Particle const * const p1, Particle const * const p2);
173 
176 
178  virtual G4double NNToNNEtaIso(const G4double ener, const G4int iso);
179 
181  virtual G4double NNToNNEtaExcluIso(const G4double ener, const G4int iso);
182 
184  virtual G4double NNToNNEtaOnePi(Particle const * const part1, Particle const * const part2);
186  virtual G4double NNToNNEtaOnePiOrDelta(Particle const * const part1, Particle const * const part2);
188  virtual G4double NNToNNEtaTwoPi(Particle const * const part1, Particle const * const part2);
190  virtual G4double NNToNNEtaThreePi(Particle const * const part1, Particle const * const part2);
192  virtual G4double NNToNNEtaFourPi(Particle const * const part1, Particle const * const part2);
193 
194 
196  virtual G4double NNToNNOmegaIso(const G4double ener, const G4int iso);
197 
199  virtual G4double NNToNNOmegaExcluIso(const G4double ener, const G4int iso);
200 
202  virtual G4double NNToNNOmegaOnePi(Particle const * const part1, Particle const * const part2);
204  virtual G4double NNToNNOmegaOnePiOrDelta(Particle const * const part1, Particle const * const part2);
206  virtual G4double NNToNNOmegaTwoPi(Particle const * const part1, Particle const * const part2);
208  virtual G4double NNToNNOmegaThreePi(Particle const * const part1, Particle const * const part2);
210  virtual G4double NNToNNOmegaFourPi(Particle const * const part1, Particle const * const part2);
211 
212 
214  virtual G4double etaNElastic(Particle const * const p1, Particle const * const p2);
215  virtual G4double omegaNElastic(Particle const * const p1, Particle const * const p2);
216 
217 
219  virtual G4double omegaNInelastic(Particle const * const p1, Particle const * const p2);
220 
222  virtual G4double omegaNToPiPiN(Particle const * const p1, Particle const * const p2);
223 
224  };
225 }
226 
227 #endif