GenHeavyIon.h

// -*- C++ -*-
//
// This file is part of HepMC
// Copyright (C) 2014-2023 The HepMC collaboration (see AUTHORS for details)
//
#ifndef HEPMC3_HEAVYION_H
#define HEPMC3_HEAVYION_H
#include <iostream>
#include <map>
#include "HepMC3/Attribute.h"
 
namespace HepMC3 {
 
class GenHeavyIon : public Attribute {
 
public:

    GenHeavyIon()
        : Ncoll_hard(-1), Npart_proj(-1), Npart_targ(-1), Ncoll(-1),
#ifndef HEPMC3_NO_DEPRECATED
          spectator_neutrons(-1), spectator_protons(-1),
#endif
          N_Nwounded_collisions(-1), Nwounded_N_collisions(-1),
          Nwounded_Nwounded_collisions(-1), impact_parameter(-1.0),
          event_plane_angle(-1.0),
#ifndef HEPMC3_NO_DEPRECATED
          eccentricity(-1.0),
#endif
          sigma_inel_NN(-1.0), centrality(-1.0), user_cent_estimate(-1.0),
          Nspec_proj_n(-1), Nspec_targ_n(-1),
          Nspec_proj_p(-1), Nspec_targ_p(-1), forceoldformat(false) {}
 
//
// Fields
//
public:

    int    Ncoll_hard;

    int    Npart_proj;

    int    Npart_targ;

    int    Ncoll;
 
#ifndef HEPMC3_NO_DEPRECATED
    int    spectator_neutrons;

    int    spectator_protons;
#endif

    int    N_Nwounded_collisions;

    int    Nwounded_N_collisions;

    int    Nwounded_Nwounded_collisions;

    double impact_parameter;

    double event_plane_angle;
 
#ifndef HEPMC3_NO_DEPRECATED
    double eccentricity;
#endif

    double sigma_inel_NN;

    double centrality;

    double user_cent_estimate;
 

    int Nspec_proj_n;

    int Nspec_targ_n;

    int Nspec_proj_p;

    int Nspec_targ_p;

    std::map<int,double> participant_plane_angles;

    std::map<int,double> eccentricities;
 
//
// Functions
//
public:

    bool from_string(const std::string &att) override;

    bool to_string(std::string &att) const  override;
 
#ifndef HEPMC3_NO_DEPRECATED

    bool operator==( const GenHeavyIon& ) const;
    bool operator!=( const GenHeavyIon& ) const;


    void set( const int&nh, const int&np, const int&nt, const int&nc, const int&ns, const int&nsp,
              const int&nnw=0, const int&nwn=0, const int&nwnw=0,
              const double& im=0., const double& pl=0., const double& ec=0., const double& s=0., const double& cent=0., const double& usrcent=0. );
    bool is_valid() const;

    bool forceoldformat;
 
#endif
 
};
 
#ifndef HEPMC3_NO_DEPRECATED
typedef GenHeavyIon HeavyIon; 
#endif
 
 
} // namespace HepMC3
 
#endif