TrentoInitial.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file TrentoInitial.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 TRENTOINITIAL_H
#define TRENTOINITIAL_H

#include <tuple>
#include <memory>
#include <iostream>
#include <boost/filesystem.hpp>
#include <boost/filesystem/fstream.hpp>
#include <boost/program_options.hpp>

#include "fwd_decl.h"
#include "JetScapeModuleBase.h"
#include "collider.h"
#include "InitialState.h"
#include "JetScapeLogger.h"

using OptDesc = po::options_description;
using VarMap = po::variables_map;
using namespace trento;

namespace Jetscape {

typedef struct {
  double impact_parameter;
  double num_participant;
  double num_binary_collisions;
  double total_entropy;
  std::map<int, double> ecc; // order, eccentricity
  std::map<int, double> psi; // order, participant_plane
  double xmid, ymid;
} EventInfo;


class TrentoInitial : public InitialState {
public:
  // Initialize from XML configuration
  TrentoInitial();
  ~TrentoInitial();

  //void Init();
  void Exec();
  void Clear();
  void InitTask();

  struct RangeFailure : public std::runtime_error {
    using std::runtime_error::runtime_error;
  };
  EventInfo info_;

private:
  std::shared_ptr<trento::Collider> TrentoGen_;
  std::pair<double, double> GenCenTab(std::string proj, std::string targ,
                                      VarMap var_map, int cL, int cH);
  // Output output_;

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

} // end namespace Jetscape

#endif // TRENTOINITIAL_H