Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/run/include/G4RNGHelper.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 // G4TemplateRNGHelper
 27 //
 28 // Class description:
 29 //
 30 // Helper class for RNG Engine seeds.
 31 // Used in MT builds to guarantee event reproducibility.
 32 // The function of this class is to return a RNG Engine seed given its index.
 33 // It is a simple templated container that allows to add seeds (AddOneSeed)
 34 // and retrieve a seed (GetSeed) by index.
 35 //
 36 // The normal use is with G4RNGHelper where each element of the container
 37 // represents a seed. To enforce strong-reproducibility the variant with
 38 // the RNG Engine status file names is available.
 39 
 40 // Author: A.Dotti (SLAC), 5 July 2013
 41 // --------------------------------------------------------------------
 42 #ifndef G4TemplateRNGHelper_hh
 43 #define G4TemplateRNGHelper_hh 1
 44 
 45 #include "globals.hh"
 46 
 47 #include <queue>
 48 #include <vector>
 49 
 50 template<class T>
 51 class G4TemplateRNGHelper
 52 {
 53   public:
 54     using SeedsQueue = std::vector<T>;
 55     using SeedsQueueSize_type = typename SeedsQueue::size_type;
 56 
 57   public:
 58     // The container is modeled as a (shared) singleton
 59     static G4TemplateRNGHelper<T>* GetInstance();
 60     static G4TemplateRNGHelper<T>* GetInstanceIfExist();
 61     virtual ~G4TemplateRNGHelper();
 62 
 63     // Returns seed given id
 64     virtual const T GetSeed(const G4int& sdId)
 65     {
 66       G4int seedId = sdId - 2 * offset;
 67       if (seedId < static_cast<G4int>(seeds.size())) {
 68         T& seed = seeds[seedId];
 69         return seed;
 70       }
 71       G4ExceptionDescription msg;
 72       msg << "No seed number " << seedId << "(" << seeds.size() << " available)\n"
 73           << " Original seed number " << sdId << " filled so far " << offset;
 74       G4Exception("G4RNGHelper::GetSeed", "Run0115", FatalException, msg);
 75       return T();
 76     }
 77 
 78     // Adds one seed to the collection
 79     void AddOneSeed(const T& seed) { seeds.push_back(seed); }
 80 
 81     // Fills N primary seed pairs
 82     void Fill(G4double* dbl, G4int nev, G4int nev_tot, G4int nrpe)
 83     {
 84       seeds.clear();
 85       for (G4int i = 0; i < nrpe * nev; ++i) {
 86         seeds.push_back((G4long)(100000000L * dbl[i]));
 87       }
 88       offset = 0;
 89       nev_filled = nev;
 90       nev_total = nev_tot;
 91       nRandParEvent = nrpe;
 92     }
 93 
 94     void Refill(G4double* dbl, G4int nev)
 95     {
 96       if (nev == 0) return;
 97       seeds.clear();
 98       for (G4int i = 0; i < nRandParEvent * nev; ++i) {
 99         seeds.push_back((G4long)(100000000L * dbl[i]));
100       }
101       offset += nev_filled;
102       nev_filled = nev;
103     }
104 
105     // Number of available seeds
106     const SeedsQueueSize_type GetNumberSeeds() const { return seeds.size(); }
107 
108     // Empty the seeds container
109     virtual void Clear() { seeds.clear(); }
110 
111   protected:
112     SeedsQueue seeds;
113     // Note: following numbers are number of events.
114     //       seeds are generated for nRandParEvent times n_event
115     G4int offset = 0;
116     G4int nev_filled = 0;
117     G4int nev_total = 0;
118     G4int nRandParEvent = 0;
119 
120   private:
121     G4TemplateRNGHelper() = default;
122 
123   private:
124     static G4TemplateRNGHelper<T>* instance;
125 };
126 
127 using G4RNGHelper = G4TemplateRNGHelper<G4long>;
128 using G4StringRNGHelper = G4TemplateRNGHelper<G4String>;
129 using G4SeedsQueue = std::queue<G4long>;
130 
131 #endif
132