Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/cuts/include/G4ProductionCutsTable.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 // G4ProductionCutsTable
 27 //
 28 // Class description:
 29 //
 30 // G4ProductionCutsTable is a singleton class for a table of
 31 // G4ProductionCuts objects. This class manages tables of production
 32 // cuts and energy cuts for each particle type.
 33 
 34 // Author: M.Asai, 5 October 2002 - First implementation
 35 // Modifications: H.Kurashige, 2004-2008
 36 // --------------------------------------------------------------------
 37 #ifndef G4ProductionCutsTable_hh 
 38 #define G4ProductionCutsTable_hh 1
 39 
 40 #include <cmath>
 41 #include <vector>
 42 
 43 #include "globals.hh"
 44 #include "G4ios.hh"
 45 #include "G4MaterialCutsCouple.hh"
 46 #include "G4MCCIndexConversionTable.hh"
 47 #include "G4Region.hh"
 48 
 49 class G4RegionStore;
 50 class G4VRangeToEnergyConverter;
 51 class G4LogicalVolume;
 52 class G4VPhysicalVolume;
 53 class G4ProductionCuts;
 54 class G4ProductionCutsTableMessenger;
 55 
 56 class G4ProductionCutsTable  
 57 {
 58   public:
 59 
 60     static G4ProductionCutsTable* GetProductionCutsTable();
 61       // This static method returns the singleton pointer of this class object.
 62       // At first invocation, the singleton object is instantiated
 63 
 64     G4ProductionCutsTable(const G4ProductionCutsTable&) = delete;
 65     G4ProductionCutsTable& operator=(const G4ProductionCutsTable&) = delete;
 66 
 67     virtual ~G4ProductionCutsTable();
 68 
 69     void CreateCoupleTables();
 70       // Creates material cuts couples table and allocate the other tables
 71 
 72     void UpdateCoupleTable(G4VPhysicalVolume* currentWorld);
 73       // Triggers an update of the table of G4ProductionCuts objects
 74 
 75     void SetEnergyRange(G4double lowedge, G4double highedge);
 76       // Sets the limits of energy cuts for all particles
 77 
 78     G4double GetLowEdgeEnergy() const;
 79     G4double GetHighEdgeEnergy() const;
 80       // Get the limits of energy cuts for all particles
 81 
 82     G4double GetMaxEnergyCut();
 83     void SetMaxEnergyCut(G4double value);
 84       // Get/set max cut energy of RangeToEnergy converter
 85       // for all particle types
 86 
 87     void DumpCouples() const;
 88       // Displays a list of registered couples
 89 
 90     const G4MCCIndexConversionTable* GetMCCIndexConversionTable() const;
 91       // Gives the pointer to the MCCIndexConversionTable
 92 
 93     const std::vector<G4double>* GetRangeCutsVector(std::size_t pcIdx) const;
 94     const std::vector<G4double>* GetEnergyCutsVector(std::size_t pcIdx) const;
 95   
 96     std::size_t GetTableSize() const;
 97       // Returns the size of the couple table
 98 
 99     const G4MaterialCutsCouple* GetMaterialCutsCouple(G4int i) const;
100       // Returns the pointer to the couple
101 
102     const G4MaterialCutsCouple* GetMaterialCutsCouple(const G4Material* aMat,
103                                             const G4ProductionCuts* aCut) const;
104      // Returns the pointer to the couple
105 
106     G4int GetCoupleIndex(const G4MaterialCutsCouple* aCouple) const;
107     G4int GetCoupleIndex(const G4Material* aMat,
108                          const G4ProductionCuts* aCut) const;
109       // Return the index of the couple.
110       // -1 is returned if index is not found
111 
112     G4bool IsModified() const;
113       // Returns TRUE if at least one production cut value is modified
114   
115     void PhysicsTableUpdated();
116       // Resets the status of IsModified(). This method must be exclusively
117       // used by the RunManager when physics tables are built
118 
119     G4ProductionCuts* GetDefaultProductionCuts() const;
120       // Returns the default production cuts
121  
122     G4double ConvertRangeToEnergy(const G4ParticleDefinition* particle,
123                                   const G4Material* material, 
124                                         G4double range);
125       // Gives energy corresponding to range value.
126       // -1 is returned if particle or material is not found
127 
128     void ResetConverters();
129       // Resets all range to energy converters  
130 
131     G4bool StoreCutsTable(const G4String& directory, 
132                           G4bool ascii = false);
133       // Stores cuts and material information in files under the
134       // the specified directory
135 
136     G4bool RetrieveCutsTable(const G4String& directory,
137                              G4bool ascii = false);
138       // Retrieve material cut couple information 
139       // in files under the specified directory
140 
141     G4bool CheckForRetrieveCutsTable(const G4String& directory, 
142                                      G4bool ascii = false);
143       // Checks stored material and cut values are consistent
144       // with the current detector setup
145 
146     G4double* GetRangeCutsDoubleVector(std::size_t pcIdx) const;
147     G4double* GetEnergyCutsDoubleVector(std::size_t pcIdx) const;
148       // Methods for backward compatibility
149 
150     void SetEnergyCutVector(const std::vector<G4double>& cutE, std::size_t idx);
151       // User defined cut vectors (idx < 4) range cut should be defined
152       // to avoid inconsistency in physics
153 
154     void  SetVerboseLevel(G4int value);
155     G4int GetVerboseLevel() const;
156       // Control flag for output message
157       //  0: Silent
158       //  1: Warning message
159       //  2: More
160 
161   protected:
162 
163     G4ProductionCutsTable();
164 
165     virtual G4bool StoreMaterialInfo(const G4String& directory, 
166                                      G4bool ascii = false);
167       // Stores material information in files under the specified directory
168 
169     virtual G4bool CheckMaterialInfo(const G4String& directory, 
170                                      G4bool ascii = false);
171       // Checks stored material is consistent with the current detector setup
172 
173     virtual G4bool StoreMaterialCutsCoupleInfo(const G4String& directory, 
174                                                G4bool ascii = false);
175       // Stores materialCutsCouple information in files under the
176       // specified directory
177 
178     virtual G4bool CheckMaterialCutsCoupleInfo(const G4String& directory, 
179                                                G4bool ascii = false);
180       // Checks stored materialCutsCouple is consistent with
181       // the current detector setup
182 
183     virtual G4bool StoreCutsInfo(const G4String& directory, 
184                                  G4bool ascii = false);
185       // Stores cut values information in files under the specified directory
186 
187     virtual G4bool  RetrieveCutsInfo(const G4String& directory,
188                                      G4bool ascii = false);
189       // Retrieves cut values information in files under the
190       // specified directory
191 
192   private:
193 
194     void ScanAndSetCouple(G4LogicalVolume* aLV,
195                           G4MaterialCutsCouple* aCouple,
196                           G4Region* aRegion);
197 
198     G4bool IsCoupleUsedInTheRegion(const G4MaterialCutsCouple* aCouple,
199                                    const G4Region* aRegion) const;
200 
201 
202   private:
203 
204     static G4ProductionCutsTable* fProductionCutsTable;
205 
206     std::vector<G4MaterialCutsCouple*> coupleTable;
207     std::vector<std::vector<G4double>*> rangeCutTable;
208     std::vector<std::vector<G4double>*> energyCutTable;
209 
210     std::vector<G4double>* userEnergyCuts[4] = {nullptr, nullptr, nullptr, nullptr};
211 
212     G4RegionStore* fG4RegionStore = nullptr;
213     G4VRangeToEnergyConverter* converters[NumberOfG4CutIndex]; 
214 
215     G4ProductionCuts* defaultProductionCuts = nullptr;
216 
217     G4MCCIndexConversionTable mccConversionTable;
218 
219     // These two vectors are for backward compatibility
220     G4double* rangeDoubleVector[NumberOfG4CutIndex];
221     G4double* energyDoubleVector[NumberOfG4CutIndex];
222 
223     enum { FixedStringLengthForStore = 32 }; 
224 
225     G4ProductionCutsTableMessenger* fMessenger = nullptr;
226     G4int verboseLevel = 1;
227     G4bool firstUse = true;
228 };
229 
230 // ------------------
231 // Inline methods
232 // ------------------
233 
234 inline 
235 const std::vector<G4double>*
236 G4ProductionCutsTable::GetRangeCutsVector(std::size_t pcIdx) const
237 { 
238   return rangeCutTable[pcIdx]; 
239 }
240 
241 inline 
242 const std::vector<G4double>*
243 G4ProductionCutsTable::GetEnergyCutsVector(std::size_t pcIdx) const
244 {
245   return energyCutTable[pcIdx]; 
246 }
247 
248 inline 
249 std::size_t G4ProductionCutsTable::GetTableSize() const
250 {
251   return coupleTable.size(); 
252 }
253 
254 inline 
255 const G4MaterialCutsCouple*
256 G4ProductionCutsTable::GetMaterialCutsCouple(G4int i) const
257 { 
258   return coupleTable[std::size_t(i)]; 
259 }
260 
261 inline 
262 G4bool G4ProductionCutsTable::IsModified() const
263 { 
264   if(firstUse) return true;
265   for(auto itr=coupleTable.cbegin(); itr!=coupleTable.cend(); ++itr)
266   { 
267     if((*itr)->IsRecalcNeeded())
268     {
269       return true; 
270     }
271   }
272   return false;
273 }
274   
275 inline 
276 void G4ProductionCutsTable::PhysicsTableUpdated()
277 { 
278   for(auto itr=coupleTable.cbegin(); itr!=coupleTable.cend(); ++itr)
279   { 
280     (*itr)->PhysicsTableUpdated(); 
281   }
282 }
283 
284 inline
285 G4double*
286 G4ProductionCutsTable::GetRangeCutsDoubleVector(std::size_t pcIdx) const
287 {
288   return rangeDoubleVector[pcIdx];
289 }
290 
291 inline
292 G4double*
293 G4ProductionCutsTable::GetEnergyCutsDoubleVector(std::size_t pcIdx) const
294 {
295   return energyDoubleVector[pcIdx];
296 }
297 
298 inline
299 G4ProductionCuts* G4ProductionCutsTable::GetDefaultProductionCuts() const
300 {
301   return defaultProductionCuts;
302 }
303 
304 inline
305 G4bool G4ProductionCutsTable::IsCoupleUsedInTheRegion(
306                                  const G4MaterialCutsCouple* aCouple,
307                                  const G4Region* aRegion) const
308 {
309   G4ProductionCuts* fProductionCut = aRegion->GetProductionCuts();
310   auto mItr = aRegion->GetMaterialIterator();
311   std::size_t nMaterial = aRegion->GetNumberOfMaterials();
312   for(std::size_t iMate=0;iMate<nMaterial; ++iMate, ++mItr)
313   {
314     if(aCouple->GetMaterial()==(*mItr) &&
315        aCouple->GetProductionCuts()==fProductionCut)
316     {
317       return true;
318     }
319   }
320   return false;
321 }
322 
323 inline
324 const G4MaterialCutsCouple* 
325 G4ProductionCutsTable::GetMaterialCutsCouple(const G4Material* aMat,
326                                              const G4ProductionCuts* aCut) const
327 {
328   for(auto cItr=coupleTable.cbegin(); cItr!=coupleTable.cend(); ++cItr)
329   { 
330     if((*cItr)->GetMaterial()!=aMat) continue;
331     if((*cItr)->GetProductionCuts()==aCut) return (*cItr);
332   }
333   return nullptr;
334 }
335 
336 inline
337 G4int
338 G4ProductionCutsTable::GetCoupleIndex(const G4MaterialCutsCouple* aCouple) const
339 { 
340   G4int idx = 0;
341   for(auto cItr=coupleTable.cbegin(); cItr!=coupleTable.cend(); ++cItr)
342   {
343     if((*cItr)==aCouple) return idx;
344     ++idx;
345   }
346   return -1;
347 }
348 
349 inline
350 G4int G4ProductionCutsTable::GetCoupleIndex(const G4Material* aMat,
351                                             const G4ProductionCuts* aCut) const
352 {
353   const G4MaterialCutsCouple* aCouple = GetMaterialCutsCouple(aMat,aCut);
354   return GetCoupleIndex(aCouple);
355 }
356 
357 inline 
358 G4int G4ProductionCutsTable::GetVerboseLevel() const
359 {
360   return verboseLevel;
361 }
362 
363 inline
364 const G4MCCIndexConversionTable* 
365 G4ProductionCutsTable::GetMCCIndexConversionTable() const
366 {
367   return &mccConversionTable;
368 }
369 
370 #endif
371