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 omegaNToPiN(Particle const * const p1, Particle const * const p2);
0060       G4double omegaNToPiPiN(Particle const * const p1, Particle const * const p2);
0061       G4double etaPrimeNToPiN(Particle const * const p1, Particle const * const p2);
0062 
0063       G4double NNToNNEta(Particle const * const p1, Particle const * const p2);
0064       G4double NNToNNEtaExclu(Particle const * const p1, Particle const * const p2);
0065       G4double NNToNNEtaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
0066       G4double NNToNDeltaEta(Particle const * const p1, Particle const * const p2);
0067       G4double NNToNNOmega(Particle const * const p1, Particle const * const p2);
0068       G4double NNToNNOmegaExclu(Particle const * const p1, Particle const * const p2);
0069       G4double NNToNNOmegaxPi(const G4int xpi, Particle const * const p1, Particle const * const p2);
0070       G4double NNToNDeltaOmega(Particle const * const p1, Particle const * const p2);
0071       
0072       /// \brief Strange cross sections
0073       G4double NNToNLK(Particle const * const p1, Particle const * const p2);
0074       G4double NNToNSK(Particle const * const p1, Particle const * const p2);
0075       G4double NNToNLKpi(Particle const * const p1, Particle const * const p2);
0076       G4double NNToNSKpi(Particle const * const p1, Particle const * const p2);
0077       G4double NNToNLK2pi(Particle const * const p1, Particle const * const p2);
0078       G4double NNToNSK2pi(Particle const * const p1, Particle const * const p2);
0079       G4double NNToNNKKb(Particle const * const p1, Particle const * const p2);
0080       G4double NNToMissingStrangeness(Particle const * const p1, Particle const * const p2);
0081       G4double NDeltaToNLK(Particle const * const p1, Particle const * const p2);
0082       G4double NDeltaToNSK(Particle const * const p1, Particle const * const p2);
0083       G4double NDeltaToDeltaLK(Particle const * const p1, Particle const * const p2);
0084       G4double NDeltaToDeltaSK(Particle const * const p1, Particle const * const p2);
0085       G4double NDeltaToNNKKb(Particle const * const p1, Particle const * const p2);
0086       G4double NpiToLK(Particle const * const p1, Particle const * const p2);
0087       G4double NpiToSK(Particle const * const p1, Particle const * const p2);
0088       G4double p_pimToSzKz(Particle const * const p1, Particle const * const p2);
0089       G4double p_pimToSmKp(Particle const * const p1, Particle const * const p2);
0090       G4double p_pizToSzKp(Particle const * const p1, Particle const * const p2);
0091       G4double NpiToLKpi(Particle const * const p1, Particle const * const p2);
0092       G4double NpiToSKpi(Particle const * const p1, Particle const * const p2);
0093       G4double NpiToLK2pi(Particle const * const p1, Particle const * const p2);
0094       G4double NpiToSK2pi(Particle const * const p1, Particle const * const p2);
0095       G4double NpiToNKKb(Particle const * const p1, Particle const * const p2);
0096       G4double NpiToMissingStrangeness(Particle const * const p1, Particle const * const p2);
0097       G4double NLToNS(Particle const * const p1, Particle const * const p2);
0098       G4double NSToNL(Particle const * const p1, Particle const * const p2);
0099       G4double NSToNS(Particle const * const p1, Particle const * const p2);
0100       G4double NKToNK(Particle const * const p1, Particle const * const p2);
0101       G4double NKToNKpi(Particle const * const p1, Particle const * const p2);
0102       G4double NKToNK2pi(Particle const * const p1, Particle const * const p2);
0103       G4double NKbToNKb(Particle const * const p1, Particle const * const p2);
0104       G4double NKbToSpi(Particle const * const p1, Particle const * const p2);
0105       G4double NKbToLpi(Particle const * const p1, Particle const * const p2);
0106       G4double NKbToS2pi(Particle const * const p1, Particle const * const p2);
0107       G4double NKbToL2pi(Particle const * const p1, Particle const * const p2);
0108       G4double NKbToNKbpi(Particle const * const p1, Particle const * const p2);
0109       G4double NKbToNKb2pi(Particle const * const p1, Particle const * const p2);
0110       G4double NYelastic(Particle const * const p1, Particle const * const p2);
0111       G4double NKbelastic(Particle const * const p1, Particle const * const p2);
0112       G4double NKelastic(Particle const * const p1, Particle const * const p2);
0113       
0114       /// \brief antiparticle cross sections
0115       /// \brief Nucleon-AntiNucleon to Baryon-AntiBaryon cross sections
0116       G4double NNbarElastic(Particle const* const p1, Particle const* const p2);
0117       G4double NNbarCEX(Particle const* const p1, Particle const* const p2);
0118 
0119       G4double NNbarToLLbar(Particle const * const p1, Particle const * const p2);
0120       
0121       /// \brief Nucleon-AntiNucleon to Nucleon-AntiNucleon + pions cross sections
0122       G4double NNbarToNNbarpi(Particle const* const p1, Particle const* const p2);
0123       G4double NNbarToNNbar2pi(Particle const* const p1, Particle const* const p2);
0124       G4double NNbarToNNbar3pi(Particle const* const p1, Particle const* const p2);
0125      
0126       /// \brief Nucleon-AntiNucleon total annihilation cross sections
0127       G4double NNbarToAnnihilation(Particle const* const p1, Particle const* const p2);     
0128       
0129       /** \brief Calculate the slope of the NN DDXS.
0130        *
0131        * \param energyCM energy in the CM frame, in MeV
0132        * \param iso total isospin of the system
0133        *
0134        * \return the slope of the angular distribution, in (GeV/c)^(-2)
0135        */
0136       G4double calculateNNAngularSlope(G4double energyCM, G4int iso);
0137 
0138       /** \brief Compute the "interaction distance".
0139        *
0140        * Defined on the basis of the average value of the N-N cross sections at
0141        * the given kinetic energy.
0142        *
0143        * \return the interaction distance
0144        */
0145       G4double interactionDistanceNN(const ParticleSpecies &aSpecies, const G4double kineticEnergy);
0146 
0147       /** \brief Compute the "interaction distance".
0148        *
0149        * Defined on the basis of the average value of the pi-N cross sections at
0150        * the given kinetic energy.
0151        *
0152        * \return the interaction distance
0153        */
0154       G4double interactionDistancePiN(const G4double projectileKineticEnergy);
0155 
0156       /** \brief Compute the "interaction distance".
0157        *
0158        * Defined on the basis of the average value of the K-N cross sections at
0159        * the given kinetic energy.
0160        *
0161        * \return the interaction distance
0162        */
0163       G4double interactionDistanceKN(const G4double projectileKineticEnergy);
0164 
0165       /** \brief Compute the "interaction distance".
0166        *
0167        * Defined on the basis of the average value of the Kbar-N cross sections at
0168        * the given kinetic energy.
0169        *
0170        * \return the interaction distance
0171        */
0172       G4double interactionDistanceKbarN(const G4double projectileKineticEnergy);
0173 
0174       /** \brief Compute the "interaction distance".
0175        *
0176        * Defined on the basis of the average value of the Y-N cross sections at
0177        * the given kinetic energy.
0178        *
0179        * \return the interaction distance
0180        */
0181       G4double interactionDistanceYN(const G4double projectileKineticEnergy);
0182 
0183       void setCrossSections(ICrossSections *c);
0184 
0185       void deleteCrossSections();
0186 
0187       void initialize(Config const * const theConfig);
0188 
0189   }
0190 }
0191 
0192 #endif