HepMCHeavyIon.hh

// -*- C++ -*-
#ifndef RIVET_HepMCHeavyIon_HH
#define RIVET_HepMCHeavyIon_HH
 
#include "Rivet/Projection.hh"
#include "Rivet/Tools/RivetHepMC.hh"
#include "Rivet/Event.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