Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-01-18 09:58:34

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 #ifndef G4LENDCrossSection_h
0027 #define G4LENDCrossSection_h 1
0028 
0029 // Class Description
0030 // Cross Sections for a LEND (Low Energy Nuclear Data)
0031 // LEND is Geant4 interface for GIDI (General Interaction Data Interface) 
0032 // which gives a discription of nuclear and atomic reactions, such as
0033 //    Binary collision cross sections
0034 //    Particle number multiplicity distributions of reaction products
0035 //    Energy and angular distributions of reaction products
0036 //    Derived calculational constants
0037 // GIDI is developped at Lawrence Livermore National Laboratory
0038 // Class Description - End
0039 
0040 // 071025 First implementation done by T. Koi (SLAC/SCCS)
0041 // 101118 Name modifications for release T. Koi (SLAC/PPA)
0042 
0043 #include "G4LENDHeader.hh"
0044 
0045 #include "G4LENDManager.hh"
0046 #include "G4LENDUsedTarget.hh"
0047 
0048 #include "G4VCrossSectionDataSet.hh"
0049 #include "G4DynamicParticle.hh"
0050 #include "G4Element.hh"
0051 #include "G4Material.hh"
0052 #include "G4ParticleDefinition.hh"
0053 #include "G4PhysicsTable.hh"
0054 
0055 #include <map>
0056 
0057 class G4LENDCrossSection : public G4VCrossSectionDataSet
0058 {
0059 
0060 //
0061 //
0062 //  G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A,
0063 //                         const G4Element* elm = 0,
0064 //                         const G4Material* mat = 0);
0065 // 
0066 //  G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A,
0067 //                              const G4Element* elm = 0,
0068 //                              const G4Material* mat = 0);
0069 //
0070 
0071    public:
0072    
0073       G4LENDCrossSection(const G4String name = "" );
0074    
0075       ~G4LENDCrossSection();
0076    
0077       //G4bool IsApplicable( const G4DynamicParticle* , const G4Element* );
0078 //TK110811
0079       G4bool IsIsoApplicable( const G4DynamicParticle* , G4int /*Z*/ , G4int /*A*/ , 
0080                               const G4Element* , const G4Material* );
0081       G4double GetIsoCrossSection( const G4DynamicParticle*, G4int /*Z*/, G4int /*A*/ ,
0082                                    const G4Isotope* , const G4Element* , const G4Material* );
0083 
0084       void BuildPhysicsTable( const G4ParticleDefinition& );
0085 
0086       void DumpPhysicsTable( const G4ParticleDefinition& );
0087       void DumpLENDTargetInfo( G4bool force = false );
0088 
0089 //TK110810
0090       //G4double GetCrossSection( const G4DynamicParticle* , const G4Element* , G4double aT );
0091       //G4double GetCrossSection(const G4DynamicParticle*, G4int , const G4Material* );
0092 
0093 
0094       void ChangeDefaultEvaluation( G4String name_tmp ){ default_evaluation = name_tmp; };
0095       void AllowNaturalAbundanceTarget(){ allow_nat = true; };
0096       void AllowAnyCandidateTarget(){ allow_any = true; };
0097 
0098       //Hadronic Framework still does not handle isotope in GPIL
0099       //G4VDiscreteProcess::PostStepGetPhysicalInteractionLenght()
0100       // G4HadronicProcess::GetMeanFreePath()
0101       //  G4double GetCrossSection(const G4DynamicParticle*, const G4Material*)
0102       //   G4double GetCrossSection(const G4DynamicParticle*, const G4Element*, G4double aTemperature); 
0103 
0104 //TK110810
0105       //G4bool IsIsoApplicable( const G4DynamicParticle* , G4int /*ZZ*/, G4int /*AA*/) { return true; }
0106       //G4bool IsZandAApplicable( const G4DynamicParticle* , G4int /*ZZ*/, G4int /*AA*/, const G4Element* , const G4Material* ) { return true; }
0107 //TK110810
0108       //G4double GetZandACrossSection(const G4DynamicParticle* , G4int /*Z*/, G4int /*A*/, G4double aTemperature);
0109       //G4double GetZandACrossSection(const G4DynamicParticle* , G4int /*Z*/, G4int /*A*/,  const G4Material* mat);
0110 //TK110810
0111       //G4double GetIsoCrossSection( const G4DynamicParticle* , const G4Isotope* , G4double );
0112 //      G4double GetIsoCrossSection( const G4DynamicParticle* , const G4Isotope* , const G4Material* mat );
0113 
0114    private:
0115 
0116       std::map< G4int , G4LENDUsedTarget* > usedTarget_map;
0117 
0118       G4String default_evaluation;
0119       G4bool allow_nat;
0120       G4bool allow_any;
0121 
0122       G4LENDManager* lend_manager;
0123       void recreate_used_target_map();
0124 
0125    protected :
0126       //G4String name;
0127       G4ParticleDefinition* proj;
0128       void create_used_target_map();
0129       G4GIDI_target* get_target_from_map( G4int nuclear_code );
0130       virtual G4double getLENDCrossSection( G4GIDI_target* , G4double , G4double ) { return 0.0; };
0131       //                                         elow       ehigh     xs_elow   xs_ehigh    ke (<elow)
0132       G4double GetUltraLowEnergyExtrapolatedXS( G4double , G4double , G4double , G4double , G4double );
0133 };
0134 #endif