BDSFieldLoaderBDSIM.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) BDSIM Collaboration, 2001 - 2026.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSFIELDLOADERBDSIM_H
#define BDSFIELDLOADERBDSIM_H
 
#include "BDSDimensionType.hh"
#include "BDSFieldValue.hh"
 
#include "G4String.hh"
 
#include <map>
#include <string>
#include <vector>
 
#ifdef USE_GZSTREAM
#include "src-external/gzstream/gzstream.h"
#endif
 
class BDSArray4DCoords;
class BDSArray3DCoords;
class BDSArray2DCoords;
class BDSArray1DCoords;
 
template <class T>
class BDSFieldLoaderBDSIM
{
public:
  BDSFieldLoaderBDSIM();
  ~BDSFieldLoaderBDSIM();
  
  BDSArray4DCoords* Load4D(const G4String& fileName); 
  BDSArray3DCoords* Load3D(const G4String& fileName); 
  BDSArray2DCoords* Load2D(const G4String& fileName); 
  BDSArray1DCoords* Load1D(const G4String& fileName); 
 
private:
  void CleanUp();
 
  void Terminate(const G4String& message = "");
 
  void Load(const G4String& fileName,
            const unsigned int nDim);
 
  void ProcessData(const std::string& line,
                   const unsigned long xIndex,
                   const unsigned long yIndex,
                   const unsigned long zIndex,
                   unsigned long long currentLineNumber);
 
  T file;
  
  unsigned long                nColumns;  
  std::vector<G4float>         lineData;  
  BDSFieldValue                fv;        
  std::map<G4String, double>   header;    
  BDSArray4DCoords*            result;    
  G4String                     loopOrder; 
  const std::vector<G4String>  headerMustBePositiveKeys;
  unsigned long indexOfFirstFieldValue;
 
  struct DimensionKeys
  {
    G4String number;
    G4String min;
    G4String max;
  };
  std::map<BDSDimensionType, DimensionKeys> dimKeyMap;
  std::map<BDSDimensionType, G4double> dimUnitsMap;
};
 
#endif