HepMCHeavyIon.hh

// -*- C++ -*-
#ifndef RIVET_HepMCHeavyIon_HH
#define RIVET_HepMCHeavyIon_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Tools/RivetHepMC.hh"
 
namespace Rivet {
 
 
  class HepMCHeavyIon : public Projection {
  public:


    HepMCHeavyIon();

    RIVET_DEFAULT_PROJ_CLONE(HepMCHeavyIon);


    using Projection::operator=;
 
 
  protected:

    void project(const Event& e);

    //int compare(const Projection& p) const;
    // Taken from Thrust.hh
    CmpState compare(const Projection&) const {
      return CmpState::EQ;
    }
 
  public:

    bool ok() const {
      return _hi != nullptr;
    }

    int Ncoll_hard() const;

    int Npart_proj() const;

    int Npart_targ() const;

    int Ncoll() const;
    int N_Nwounded_collisions() const;

    int Nwounded_N_collisions() const;

    int Nwounded_Nwounded_collisions() const;

    double impact_parameter() const;

    double event_plane_angle() const;
    double sigma_inel_NN() const;

    double centrality() const;

    double user_cent_estimate() const;

    int Nspec_proj_n() const;

    int Nspec_targ_n() const;

    int Nspec_proj_p() const;

    int Nspec_targ_p() const;

    map<int, double> participant_plane_angles() const;

    map<int, double> eccentricities() const;
 
 
  protected:

    ConstGenHeavyIonPtr _hi;
  };
}
 
 
#endif