#ifndef HGEOMMEDIUM_H
#define HGEOMMEDIUM_H

using namespace std;
#include "TNamed.h"
#include <fstream>
#include <iomanip>

class HGeomMedium  : public TNamed {
  Int_t     medId;       // Unique number for the material and the medium
  Int_t     autoflag;    // Flag for default values of madfld .. minstep
  Int_t     nComponents; // Number of material components
  Int_t     weightFac;   // Factor for weights (1: relative w., -1: w. by number of atoms) 
  Double_t* ca;          // Atomic weights for the components
  Double_t* cz;          // Atomic numbers for the components
  Double_t* cw;          // Weights of the components in a mixture
  Double_t  density;     // Density in g cm(**-3)
  Double_t  radLen;      // Radiation lenght
  Int_t     sensFlag;    // Sensitivity flag
  Int_t     fldFlag;     // Field flag
  Double_t  fld;         // Maximum field value in kilogauss
  Double_t  epsil;       // Boundary crossing precision
  Double_t  madfld;      // Maximum angular deviation due to field
  Double_t  maxstep;     // Maximum step permitted
  Double_t  maxde;       // Maximum fractional energy loss
  Double_t  minstep;     // Minimum value for step
  Int_t     npckov;      // Number of optical parameters for Cerenkov
  Double_t* ppckov;      // Photon momentum
  Double_t* absco;       // Absoption length
  Double_t* effic;       // Detection efficiency
  Double_t* rindex;      // Refraction index
public:
  HGeomMedium(const Char_t* name="");
  ~HGeomMedium();
  void setMediumIndex(Int_t i) {medId=i;}
  Int_t getMediumIndex() {return medId;}
  void setAutoFlag(Int_t i) {autoflag=i;}
  Int_t getAutoFlag() {return autoflag;}
  void setNComponents(Int_t);
  Int_t getNComponents() {return nComponents;}
  Bool_t setComponent(Int_t,Double_t,Double_t,Double_t w=1.);
  void getComponent(Int_t,Double_t*);
  Int_t getWeightFac() {return weightFac;} 
  void setDensity(Double_t d) {density=d;}
  Double_t getDensity() {return density;}
  void setRadiationLength(Double_t l) {radLen=l;}
  Double_t getRadiationLength() {return radLen;} 
  void setNpckov(Int_t);
  Int_t getNpckov() {return npckov;}
  Bool_t setCerenkovPar(Int_t,Double_t,Double_t,Double_t,Double_t);
  void getCerenkovPar(Int_t,Double_t*);  
  void setMediumPar(Int_t,Int_t,Double_t,Double_t,Double_t maxDeviation=0.,
                    Double_t maxStep=0.,Double_t maxDE=0.,
                    Double_t minStepDouble_t=0.);
  void getMediumPar(Double_t*);
  inline Bool_t isSensitive();
  Int_t getSensitivityFlag() { return sensFlag;}
  Int_t getFieldFlag() {return fldFlag;}
  Double_t getField() {return fld;}
  Double_t getEpsil() {return epsil;}
  void read(fstream&, Int_t autoflag );
  void print();
  void write (fstream&);
  Bool_t calcRadiationLength();
private:
  ClassDef(HGeomMedium,0) // Class for tracking medium 
};

inline Bool_t HGeomMedium::isSensitive() {
  if (sensFlag) return kTRUE;
  else return kFALSE; 
}

#endif /* !HGEOMMEDIUM_H */