Geant4 Cross Reference |
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 27 #ifndef G4IonsShenCrossSection_h 28 #define G4IonsShenCrossSection_h 29 // 30 // Class Description 31 // Implementation of formulas 32 // Shen et al. Nuc. Phys. A 491 130 (1989); 33 // Total Reaction Cross Section for Heavy-Ion Collisions 34 // 35 // Class Description - End 36 // 18-Sep-2003 First version is written by T. Koi 37 // 12-Nov-2003 Set upper limit at 10 GeV/n 38 // 12-Nov-2003 Insted of the lower limit, 39 // 0 is returned to a partilce with energy lowae than 10 MeV/n 40 // 15-Nov-2006 Change upper limit to 1 TeV/n 41 // However above 10GeV/n XS become constant. 42 // 23-Dec-2006 Isotope dependence added by D. Wright 43 // 19-Aug-2011 V.Ivanchenko move to new design and make x-section per element 44 // 45 46 #include "globals.hh" 47 #include "G4Proton.hh" 48 49 #include "G4VCrossSectionDataSet.hh" 50 51 class G4IonsShenCrossSection : public G4VCrossSectionDataSet 52 { 53 public: 54 55 G4IonsShenCrossSection(); 56 57 virtual ~G4IonsShenCrossSection(); 58 59 virtual 60 G4bool IsElementApplicable(const G4DynamicParticle* aDP, 61 G4int Z, const G4Material*); 62 63 virtual 64 G4double GetElementCrossSection(const G4DynamicParticle*, 65 G4int Z, const G4Material*); 66 67 virtual 68 G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A, 69 const G4Isotope* iso = 0, 70 const G4Element* elm = 0, 71 const G4Material* mat = 0); 72 73 virtual void CrossSectionDescription(std::ostream&) const; 74 75 private: 76 const G4double upperLimit; 77 // const G4double lowerLimit; 78 const G4double r0; 79 80 G4double calEcmValue(const G4double, const G4double, const G4double); 81 G4double calCeValue(const G4double); 82 }; 83 84 #endif 85