include/Geant4/G4PenelopeIonisationCrossSection.hh

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 //
0027 // Author: Luciano Pandola
0028 //
0029 // History:
0030 // -----------
0031 // 14 Mar 2012   L. Pandola   1st implementation.
0032 //
0033 // -------------------------------------------------------------------
0034 //
0035 //! Class description:
0036 //!  Access to Penelope ionisation cross sections for e- and e+. To be
0037 //!  registered as a specific model in G4UAtomicDeexcitation.
0038 //!
0039 //! NOTICE: working only for e- at the moment (no interface available for
0040 //!         e+)
0041 //!
0042 // -------------------------------------------------------------------
0043
0044 #ifndef G4PENELOPEIONISATIONCROSSSECTION_HH
0045 #define G4PENELOPEIONISATIONCROSSSECTION_HH 1
0046
0047 #include "globals.hh"
0048 #include "G4VhShellCrossSection.hh"
0049 #include <map>
0050
0051 class G4AtomicTransitionManager;
0052 class G4PenelopeOscillatorManager;
0053 class G4PenelopeOscillator;
0054 class G4PenelopeCrossSection;
0055 class G4PenelopeIonisationXSHandler;
0056
0057 class G4PenelopeIonisationCrossSection : public G4VhShellCrossSection
0058 {
0059 public:
0060   //! Constructor.
0061   explicit G4PenelopeIonisationCrossSection();
0062
0063   //! Destructor. Clean all tables.
0064   ~G4PenelopeIonisationCrossSection();
0065
0066   //! Purely virtual method from the base interface. Returns the cross
0067   //! section for all levels of element Z in material mat at the
0068   //! given energy
0069   std::vector<G4double> GetCrossSection(G4int Z,
0070                     G4double incidentEnergy,
0071                     G4double mass,
0072                     G4double deltaEnergy,
0073                     const G4Material* mat) override;
0074
0075   //! Purely virtual method from the base interface. Returns the
0076   //! cross section for the given shell in the element Z of material
0077   //! mat at the specified energy
0078   G4double CrossSection(G4int Z,
0079             G4AtomicShellEnumerator shell,
0080             G4double incidentEnergy,
0081             G4double mass,
0082             const G4Material* mat) override;
0083
0084   //! Purely virtual method from the base interface. Returns the
0085   //! shell ionisation probabilities for the given Z in the
0086   //! material mat at the specified energy.
0087   std::vector<G4double> Probabilities(G4int Z,
0088                       G4double incidentEnergy,
0089                       G4double mass,
0090                       G4double deltaEnergy,
0091                       const G4Material* mat) override;
0092   //! Getter/setter for the verbosity level
0093   void SetVerbosityLevel(G4int vl){fVerboseLevel = vl;};
0094   G4int GetVerbosityLevel(){return fVerboseLevel;};
0095
0096   G4PenelopeIonisationCrossSection & operator=(const G4PenelopeIonisationCrossSection &right)
0097   = delete;
0098   G4PenelopeIonisationCrossSection(const G4PenelopeIonisationCrossSection&) = delete;
0099
0100 private:
0101   //!The shells in Penelope are organized per *material*, rather than per
0102   //!element, so given a material one has to find the proper index for the
0103   //!given Z and shellID. An appropriate look-up table is used to avoid
0104   //!recalculation.
0105   G4int FindShellIDIndex(const G4Material* mat,G4int Z,G4AtomicShellEnumerator shell);
0106
0107   std::map< std::pair<const G4Material*,G4int>, G4DataVector*> *fShellIDTable;
0108
0109   //Oscillator manager
0110   G4PenelopeOscillatorManager* fOscManager;
0111   G4PenelopeIonisationXSHandler* fCrossSectionHandler;
0112   const G4AtomicTransitionManager* fTransitionManager;
0113
0114   G4double fLowEnergyLimit;
0115   G4double fHighEnergyLimit;
0116   G4int fVerboseLevel;
0117   G4int fNMaxLevels;
0118 };
0119
0120 #endif
0121