#ifndef __PNDDRCOPTMATLIF_H__
#define __PNDDRCOPTMATLIF_H__

/*! 
  \brief Class for optical material representation.

  In this case lithium fluoride. See 
  http://ircamera.as.arizona.edu/nircam/materials/SPIE_Aug05/5904-25_Leviton.pdf

*/

class PndDrcOptMatLiF : public PndDrcOptMatAbs
{
 protected:

  //!                                   Sellmeier coefficient for refractive index.
  double  m_B1;                         
  //!                                   Sellmeier coefficient for refractive index.  
  double  m_B2;                         
  //!                                   Sellmeier coefficient for refractive index.
  double  m_B3;                         
  //!                                   Sellmeier coefficient for refractive index.
  double  m_C1;                         
  //!                                   Sellmeier coefficient for refractive index.  
  double  m_C2;                         
  //!                                   Sellmeier coefficient for refractive index.
  double  m_C3;                         
  //!                                   Random generator.
  mutable TRandom m_ran;                

 private:
  /*! \brief Auxiliary function for assignment operator and copy constructor..
    \param s The object to copy.
  */
  void copy(const PndDrcOptMatLiF& s);         

  /*! \brief calculate Sellmeier coefficients which are here given as function 
    of temperature. Source is the paper 
    http://ircamera.as.arizona.edu/nircam/materials/SPIE_Aug05/5904-25_Leviton.pdf

    \param b1 Sellmeier coefficient.
    \param b2 Sellmeier coefficient.
    \param b3 Sellmeier coefficient.
    \param c1 Sellmeier coefficient [nm2].
    \param c2 Sellmeier coefficient [nm2].
    \param c3 Sellmeier coefficient [nm2].
    \param temp temperature in kelvin.
  */
  void calcSellmeier(double& b1, double& b2, double& b3,
		     double& c1, double& c2, double& c3,
		     double temp);
 public:

  //!                                             Empty constructor.
  PndDrcOptMatLiF();                

  /*! Copy constructor.
    /param Object to copy.
  */ 
  PndDrcOptMatLiF(const PndDrcOptMatLiF& mat);

  /*! Assignment operator.
    /param Object to assign.
  */ 
  PndDrcOptMatLiF& operator=(const PndDrcOptMatLiF& mat);

  // Implementation of pure virtual functions from base class.
  PndDrcOptMatLiF* clone() const;
  double refIndex(const double lambda) const;
  double refIndexDeriv(const double lambda) const;
  bool   absorptionFlag(double lambda, double length) const;  //!< not yet implemented.

 
};
#endif