Back to home page

EIC code displayed by LXR

 
 

    


Warning, file /include/Geant4/G4INCLCrossSections.hh was not indexed or was modified since last indexation (in which case cross-reference links may be missing, inaccurate or erroneous).

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // INCL++ intra-nuclear cascade model
0027 // Alain Boudard, CEA-Saclay, France
0028 // Joseph Cugnon, University of Liege, Belgium
0029 // Jean-Christophe David, CEA-Saclay, France
0030 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
0031 // Sylvie Leray, CEA-Saclay, France
0032 // Davide Mancusi, CEA-Saclay, France
0033 //
0034 #define INCLXX_IN_GEANT4_MODE 1
0035 
0036 #include "globals.hh"
0037 
0038 #ifndef G4INCLCrossSections_hh
0039 #define G4INCLCrossSections_hh 1
0040 
0041 #include "G4INCLICrossSections.hh"
0042 #include "G4INCLConfig.hh"
0043 
0044 namespace G4INCL {
0045   namespace CrossSections {
0046       G4double elastic(Particle const * const p1, Particle const * const p2);
0047       G4double total(Particle const * const p1, Particle const * const p2);
0048 
0049       G4double NDeltaToNN(Particle const * const p1, Particle const * const p2);
0050       G4double NNToNDelta(Particle const * const p1, Particle const * const p2);
0051       G4double NNToxPiNN(const G4int xpi, Particle const * const p1, Particle const * const p2);
0052       G4double piNToDelta(Particle const * const p1, Particle const * const p2);
0053       G4double piNToxPiN(const G4int xpi, Particle const * const p1, Particle const * const p2);
0054       G4double piNToEtaN(Particle const * const p1, Particle const * const p2);
0055       G4double piNToOmegaN(Particle const * const p1, Particle const * const p2);
0056       G4double piNToEtaPrimeN(Particle const * const p1, Particle const * const p2);
0057         G4double etaNToPiN(Particle const * const p1, Particle const * const p2);
0058         G4double etaNToPiPiN(Particle const * const p1, Particle const * const p2);
0059       G4double etaNToLK(Particle const * const p1, Particle const * const p2);
0060       G4double etaNToSK(Particle const * const p1, Particle const * const p2);
0061       G4double omegaNToPiN(Particle const * const p1, Particle const * const p2);
0062       G4double omegaNToPiPiN(Particle const * const p1, Particle const * const p2);
0063       G4double omegaNToLK(Particle const * const p1, Particle const * const p2);
0064       G4double omegaNToSK(Particle const * const p1, Particle const * const p2);
0065       G4double etaPrimeNToPiN(Particle const * const p1, Particle const * const p2);
0066 
0067       G4double NNToNNEta(Particle const * const p1, Particle const * const p2);
0068       G4double NNToNNEtaExclu(Particle const * const p1, Particle const * const p2);
0069       G4double NNToNNEtaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
0070       G4double NNToNDeltaEta(Particle const * const p1, Particle const * const p2);
0071       G4double NNToNNOmega(Particle const * const p1, Particle const * const p2);
0072       G4double NNToNNOmegaExclu(Particle const * const p1, Particle const * const p2);
0073       G4double NNToNNOmegaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
0074       G4double NNToNDeltaOmega(Particle const * const p1, Particle const * const p2);
0075       
0076       /// \brief Strange cross sections
0077       G4double NNToNLK(Particle const * const p1, Particle const * const p2);
0078       G4double NNToNSK(Particle const * const p1, Particle const * const p2);
0079       G4double NNToNLKpi(Particle const * const p1, Particle const * const p2);
0080       G4double NNToNSKpi(Particle const * const p1, Particle const * const p2);
0081       G4double NNToNLK2pi(Particle const * const p1, Particle const * const p2);
0082       G4double NNToNSK2pi(Particle const * const p1, Particle const * const p2);
0083       G4double NNToNNKKb(Particle const * const p1, Particle const * const p2);
0084       G4double NNToMissingStrangeness(Particle const * const p1, Particle const * const p2);
0085       G4double NDeltaToNLK(Particle const * const p1, Particle const * const p2);
0086       G4double NDeltaToNSK(Particle const * const p1, Particle const * const p2);
0087       G4double NDeltaToDeltaLK(Particle const * const p1, Particle const * const p2);
0088       G4double NDeltaToDeltaSK(Particle const * const p1, Particle const * const p2);
0089       G4double NDeltaToNNKKb(Particle const * const p1, Particle const * const p2);
0090       G4double NpiToLK(Particle const * const p1, Particle const * const p2);
0091       G4double NpiToSK(Particle const * const p1, Particle const * const p2);
0092       G4double p_pimToSzKz(Particle const * const p1, Particle const * const p2);
0093       G4double p_pimToSmKp(Particle const * const p1, Particle const * const p2);
0094       G4double p_pizToSzKp(Particle const * const p1, Particle const * const p2);
0095       G4double NpiToLKpi(Particle const * const p1, Particle const * const p2);
0096       G4double NpiToSKpi(Particle const * const p1, Particle const * const p2);
0097       G4double NpiToLK2pi(Particle const * const p1, Particle const * const p2);
0098       G4double NpiToSK2pi(Particle const * const p1, Particle const * const p2);
0099       G4double NpiToNKKb(Particle const * const p1, Particle const * const p2);
0100       G4double NpiToMissingStrangeness(Particle const * const p1, Particle const * const p2);
0101       G4double NLToNS(Particle const * const p1, Particle const * const p2);
0102       G4double NSToNL(Particle const * const p1, Particle const * const p2);
0103       G4double NSToNS(Particle const * const p1, Particle const * const p2);
0104       G4double NKToNK(Particle const * const p1, Particle const * const p2);
0105       G4double NKToNKpi(Particle const * const p1, Particle const * const p2);
0106       G4double NKToNK2pi(Particle const * const p1, Particle const * const p2);
0107       G4double NKbToNKb(Particle const * const p1, Particle const * const p2);
0108       G4double NKbToSpi(Particle const * const p1, Particle const * const p2);
0109       G4double NKbToLpi(Particle const * const p1, Particle const * const p2);
0110       G4double NKbToS2pi(Particle const * const p1, Particle const * const p2);
0111       G4double NKbToL2pi(Particle const * const p1, Particle const * const p2);
0112       G4double NKbToNKbpi(Particle const * const p1, Particle const * const p2);
0113       G4double NKbToNKb2pi(Particle const * const p1, Particle const * const p2);
0114       G4double NYelastic(Particle const * const p1, Particle const * const p2);
0115       G4double NKbelastic(Particle const * const p1, Particle const * const p2);
0116       G4double NKelastic(Particle const * const p1, Particle const * const p2);
0117       
0118       /// \brief antiparticle cross sections
0119       /// \brief Nucleon-AntiNucleon to Baryon-AntiBaryon cross sections
0120       G4double NNbarElastic(Particle const* const p1, Particle const* const p2);
0121       G4double NNbarCEX(Particle const* const p1, Particle const* const p2);
0122 
0123       G4double NNbarToLLbar(Particle const * const p1, Particle const * const p2);
0124       
0125       /// \brief Nucleon-AntiNucleon to Nucleon-AntiNucleon + pions cross sections
0126       G4double NNbarToNNbarpi(Particle const* const p1, Particle const* const p2);
0127       G4double NNbarToNNbar2pi(Particle const* const p1, Particle const* const p2);
0128       G4double NNbarToNNbar3pi(Particle const* const p1, Particle const* const p2);
0129      
0130       /// \brief Nucleon-AntiNucleon total annihilation cross sections
0131       G4double NNbarToAnnihilation(Particle const* const p1, Particle const* const p2);     
0132       
0133       /** \brief Calculate the slope of the NN DDXS.
0134        *
0135        * \param energyCM energy in the CM frame, in MeV
0136        * \param iso total isospin of the system
0137        *
0138        * \return the slope of the angular distribution, in (GeV/c)^(-2)
0139        */
0140       G4double calculateNNAngularSlope(G4double energyCM, G4int iso);
0141 
0142       /** \brief Compute the "interaction distance".
0143        *
0144        * Defined on the basis of the average value of the N-N cross sections at
0145        * the given kinetic energy.
0146        *
0147        * \return the interaction distance
0148        */
0149       G4double interactionDistanceNN(const ParticleSpecies &aSpecies, const G4double kineticEnergy);
0150       G4double interactionDistanceNbarN(const ParticleSpecies &aSpecies, const G4double kineticEnergy);
0151       G4double interactionDistancenbarN(const ParticleSpecies &aSpecies, const G4double kineticEnergy);
0152       /** \brief Compute the "interaction distance".
0153        *
0154        * Defined on the basis of the average value of the pi-N cross sections at
0155        * the given kinetic energy.
0156        *
0157        * \return the interaction distance
0158        */
0159       G4double interactionDistancePiN(const G4double projectileKineticEnergy);
0160 
0161       /** \brief Compute the "interaction distance".
0162        *
0163        * Defined on the basis of the average value of the K-N cross sections at
0164        * the given kinetic energy.
0165        *
0166        * \return the interaction distance
0167        */
0168       G4double interactionDistanceKN(const G4double projectileKineticEnergy);
0169 
0170       /** \brief Compute the "interaction distance".
0171        *
0172        * Defined on the basis of the average value of the Kbar-N cross sections at
0173        * the given kinetic energy.
0174        *
0175        * \return the interaction distance
0176        */
0177       G4double interactionDistanceKbarN(const G4double projectileKineticEnergy);
0178 
0179       /** \brief Compute the "interaction distance".
0180        *
0181        * Defined on the basis of the average value of the Y-N cross sections at
0182        * the given kinetic energy.
0183        *
0184        * \return the interaction distance
0185        */
0186       G4double interactionDistanceYN(const G4double projectileKineticEnergy);
0187 
0188       void setCrossSections(ICrossSections *c);
0189 
0190       void deleteCrossSections();
0191 
0192       void initialize(Config const * const theConfig);
0193 
0194   }
0195 }
0196 
0197 #endif