collider.h

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

// TRENTO: Reduced Thickness Event-by-event Nuclear Topology
// Copyright 2015 Jonah E. Bernhard, J. Scott Moreland
// TRENTO3D: Three-dimensional extension of TRENTO by Weiyao Ke
// MIT License

#ifndef COLLIDER_H
#define COLLIDER_H

#include <memory>

#include "fwd_decl.h"
#include "event.h"
#include "nucleon.h"
#include "output.h"

namespace trento {


struct records {
        int i;
        double b;
        double npart;
        double mult;
};

class Collider {
 public:
  explicit Collider(const VarMap& var_map);

  ~Collider();

  void run_events();
  const std::vector<records> & all_records() const{
        return all_records_;
  }
  const Event & expose_event() const{
    return event_;
  }

 private:
  // Most of these are pretty self-explanatory...

  double sample_impact_param();
  std::unique_ptr<Nucleus> nucleusA_, nucleusB_;

  NucleonProfile nucleon_profile_;

  const int nevents_;

  int ntrys_;

  const double bmin_, bmax_;

  const int npartmin_, npartmax_;

  const double stotmin_, stotmax_;

  const double asymmetry_;

  Event event_;

  Output output_;

  bool with_ncoll_;

  // take down id, b, npart, ncoll, mult for each event.
  std::vector<records> all_records_;
};

}  // namespace trento

#endif  // COLLIDER_H