include/Geant4/G4FermiNucleiProperties.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 // G4FermiBreakUpAN alternative FermiBreakUp model
0028 // by A. Novikov (January 2025)
0029 //
0030
0031 #ifndef G4FERMINUCLEIPROPERTIES_HH
0032 #define G4FERMINUCLEIPROPERTIES_HH
0033
0034 #include "G4FermiDataTypes.hh"
0035 #include "globals.hh"
0036 #include <vector>
0037
0038 // Caches values from larger G4NucleiProperties(x5-10 speed boost)
0039 class G4FermiNucleiProperties
0040 {
0041   public:
0042     void Initialize() { *this = G4FermiNucleiProperties(); }
0043
0044     template<typename DataSource>
0045     void Initialize(const DataSource& dataSource)
0046     {
0047       Initialize(dataSource.begin(), dataSource.end());
0048     }
0049
0050     template<typename Iter>
0051     void Initialize(Iter begin, Iter end)
0052     {
0053       nucleiMasses_.clear();
0054       static_assert(
0055         std::is_same_v<typename Iter::value_type, std::pair<const G4FermiNucleiData, G4double>>,
0056         "invalid iterator");
0057       for (auto it = begin; it != end; ++it) {
0058         InsertNuclei(it->first.atomicMass, it->first.chargeNumber, it->second);
0059       }
0060     }
0061
0062     static G4double GetNuclearMass(G4FermiAtomicMass atomicMass, G4FermiChargeNumber chargeNumber)
0063     {
0064       return Instance().GetNuclearMassImpl(atomicMass, chargeNumber);
0065     }
0066
0067     static G4bool IsStable(G4FermiAtomicMass atomicMass, G4FermiChargeNumber chargeNumber)
0068     {
0069       return Instance().IsStableImpl(atomicMass, chargeNumber);
0070     }
0071
0072     void InsertNuclei(G4FermiAtomicMass atomicMass, G4FermiChargeNumber chargeNumber, G4double mass,
0073                       G4bool isStable = true);
0074
0075     static G4FermiNucleiProperties& Instance()
0076     {
0077       static G4FermiNucleiProperties properties;
0078       return properties;
0079     }
0080
0081   private:
0082     G4FermiNucleiProperties();
0083
0084     G4double GetNuclearMassImpl(G4FermiAtomicMass atomicMass,
0085                                 G4FermiChargeNumber chargeNumber) const;
0086
0087     G4bool IsStableImpl(G4FermiAtomicMass atomicMass, G4FermiChargeNumber chargeNumber) const;
0088
0089     struct G4FermiMassData
0090     {
0091         G4double mass;
0092
0093         G4bool isStable = false;  // is nuclei stable
0094
0095         G4bool isCached = false;  // value has been inserted earlier
0096     };
0097
0098     mutable std::vector<G4FermiMassData> nucleiMasses_;
0099 };
0100
0101 #endif  // G4FERMINUCLEIPROPERTIES_HH