#ifndef HKalMixture_h
#define HKalMixture_h

// from ROOT
#include "TMixture.h"


class HKalMixture : public TMixture {

private:
    Float_t  fExciteEner;     // Mean excitation energy of mixture in MeV.

    Float_t *fDensityMixt;    // Stores densities in g/mm^3 of components.
    Float_t *fExciteEnerMixt; // Stores mean excitation energy in MeV of components.
    Float_t *fRadLengthMixt;  // Stores radiation length in mm of components.

public:

    HKalMixture();
    //HKalMixture(const char *name, const char *title, const Float_t mat[5], Float_t inter=0.);
    //HKalMixture(const char *name, const char *title, Float_t a, Float_t z, Float_t dense, Float_t exener, Float_t radl, Float_t inter=0.);
    HKalMixture(const char *name, const char *title, Int_t nmixt);
    virtual ~HKalMixture();

    static  void     calcMassFracFromVolFrac(Float_t w[], Int_t n, const Float_t rho[], const Float_t vol[]);
    virtual void     calcMixture      ();
    virtual void     defineElement    (Int_t n, const HKalMixture &mat, Float_t w);
    virtual void     defineElement    (Int_t n, const Float_t mat[5], Float_t w);
    virtual void     defineElement    (Int_t n, Float_t a, Float_t z, Float_t dense, Float_t exener, Float_t radl, Float_t w);

    virtual Float_t* GetDensityMixt   () const { return fDensityMixt;    }
    virtual Float_t  GetExciteEner    () const { return fExciteEner;     }
    virtual Float_t* GetExciteEnerMixt() const { return fExciteEnerMixt; }
    virtual Float_t* GetRadLengthMixt () const { return fRadLengthMixt;  }

    virtual void     print            (const char *opt="") const;

    ClassDef(HKalMixture, 0)
};

#endif // HKalMixture_h