ConexInterface.h

Go to the documentation of this file.

#ifndef _ConexShowerGeneratorKG_ConexInterface_h_
#define _ConexShowerGeneratorKG_ConexInterface_h_

#include <ConexDynamicInterface.h>
#include <io/CONEXStructures.h>

#include <string>
#include <vector>


namespace io {
  class CONEXHeader;
  class CONEXShower;
  class CONEXLeadingParticles;
}


class LeadingInteractionsData {
public:
  LeadingInteractionsData() : pId(0), pEnergy(0), energyCM(0), kinel(0), mult(0),
                              matg(0), Depth(0), Height(0) {}
  LeadingInteractionsData(int _pid, double _kinel, double _pEnergy, double _energyCM,
                          int _mult, int _matg, double _depth, double _height) :
    pId(_pid), pEnergy(_pEnergy), energyCM(_energyCM), kinel(_kinel), mult(_mult),
    matg(_matg), Depth(_depth),
    Height(_height) {}
  int pId;           // parent id
  double pEnergy;    // parent energy
  double energyCM;   // center-of-mass energy
  double kinel;      // inelasticity
  int mult;          // multiplicity
  int matg;          // target nucleus mass
  double Depth;      // depth
  double Height;     // height
};


class LeadingInteractionsParticleData {
public:
  LeadingInteractionsParticleData(double e, double m, double px, double py, double pz,
                                  int id, int nInt)
  {fEnergy=e; fMass=m; fpx=px; fpy=py; fpz=pz; fID=id; fnInt=nInt;}
  double fEnergy;
  double fMass;
  double fpx;
  double fpy;
  double fpz;
  int fID;
  int fnInt;
private:
  LeadingInteractionsParticleData();
};


namespace ConexShowerGeneratorKG {

  class ConexInterface : public ConexDynamicInterface {

  public:
    ConexInterface(const EHEModel model,
                   const std::string& parameterPathName,
                   const int StoredGenerations,
                   const int seed);
    ~ConexInterface();

    void SimulateShower(int particleType,
                        double energyInGeV,
                        double zenith,
                        double azimuth,
                        double altitude,
                        int storedParticles);

    static bool ParticleSort(const LeadingInteractionsParticleData& lhs,
                             const LeadingInteractionsParticleData& rhs);

    const io::CONEXHeader& GetHeader() const { return *fHeader; }
    const io::CONEXShower& GetShower() const { return *fShower; }
    const std::vector<io::CONEXLeadingParticles>& GetLPvector() const { return fLeadingParticles; }

  private:
    ConexInterface();

    int fRandomSeeds[3] = { 0 };
    io::CONEXHeader* fHeader = nullptr;
    io::CONEXShower* fShower = nullptr;
    std::vector<io::CONEXLeadingParticles> fLeadingParticles;
  };

}


#endif