InitialFromFile.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file InitialFromFile.h

/*******************************************************************************
 * Copyright (c) The JETSCAPE Collaboration, 2018
 *
 * Modular, task-based framework for simulating all aspects of heavy-ion collisions
 * 
 * For the list of contributors see AUTHORS.
 *
 * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
 *
 * or via email to bugs.jetscape@gmail.com
 *
 * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
 * See COPYING for details.
 ******************************************************************************/

#ifndef INITIALFROMFILE_H
#define INITIALFROMFILE_H

#include <iostream>
#include <string>
#include <sstream>
#include <cmath>
#include "hdf5.h"
#include "Hydroinfo_h5.h"
#include "JetScapeModuleBase.h"
#include "InitialState.h"
#include "JetScapeLogger.h"

using namespace Jetscape;

class InitialFromFile : public Jetscape::InitialState {
  // this is wrapper class to read external files that
  // stores initial number of binary collisions and corresponding
  // configurations
public:
  InitialFromFile();
  ~InitialFromFile();

  //void Init();

  void Exec();

  void Clear();

  void InitTask();

  virtual void Write(weak_ptr<JetScapeWriter> w);

  double GetNpart() { return npart; };

  double GetNcoll() { return ncoll; };

  double GetTotalEntropy() { return totalentropy; };

private:
  // the hdf5 file pointer, e.g. *.hdf5
  hid_t H5file_ptr_;

  // the hdf5/group pointer, e.g. /event0
  hid_t H5group_ptr_;

  void ReadConfigs();

  void ReadNbcDist();

  void ReadEntropyDist();

  HydroinfoH5 *h5_helper_;

  int dim_x_, dim_y_;

  double npart = -1;
  double ncoll = -1;
  double totalentropy = -1;

  // Allows the registration of the module so that it is available to be used by the Jetscape framework.
  static RegisterJetScapeModule<InitialFromFile> reg;
};

#endif // INITIALFROMFILE_H