sersic_source.h

/*
 * sersic_source.h
 *
 *  Created on: Mar 6, 2010
 *      Author: R.B. Metcalf
 */
#ifndef SERSIC_SOURCE_H_
#define SERSIC_SOURCE_H_
 
#include "source.h"
 
 
class SourceSersic : public SourceColored
{
public:
  SourceSersic();
  SourceSersic(
               double my_mag            
               ,double my_Reff          
               ,double my_PA            
               ,double my_index         
               ,double my_q             
               ,double my_z             
               ,double my_zeropoint     
               //,const double *my_theta=0  /// optional angular position on the sky
  );
  SourceSersic(
                double my_mag            
                ,Band band               
                ,double my_zeropoint     
                ,double my_Reff          
                ,double my_PA            
                ,double my_index         
                ,double my_q             
                ,double my_z             
  );
 ~SourceSersic();
  
  SourceSersic(const SourceSersic &p);
  SourceSersic & operator=(const SourceSersic &p);
 
 
  void ReSet(PosType mag,Band band,PosType Reff,PosType PA,PosType my_index,PosType my_q,PosType my_z,const PosType *theta=0);
 
 inline PosType FractionRadius (PosType f) {return Reff*pow(-log (f)/bn,index);}
 
 inline PosType getSersicIndex() const { return index; }
 inline PosType getAxesRatio() const { return q; }
 inline PosType getReff() const { return Reff/arcsecTOradians; }
 //inline PosType getMag() const { return mag; }
 inline PosType getPA() const { return PA; }
 
 inline void setSersicIndex(PosType x)
 {
  index = x;
  bn = 1.9992*index - 0.3271; // approximation valid for 0.5 < n < 8
  I_n = 1./index*pow(bn,2*index)/tgamma(2*index);
  updateRadius();
 }
 
 inline void setAxesRatio(PosType x)
 {
  q = x;
  I_q = 1./q;
 }
 
 inline void setReff(PosType x // in arcseconds
  )
 {
  Reff = x*arcsecTOradians;
  I_r = 1./2./PI/Reff/Reff;
  updateRadius();
 }
 
 inline void setPA(PosType x)
 {
  PA = x;
    cosPA = cos(x);
    sinPA = sin(x);
 }
  
  void rotate(PosType theta){
    PA = PA + theta;
    cosPA = cos(PA);
    sinPA = sin(PA);
  }
 
 //inline PosType getTotalFlux() const { return flux_total; }
 
 PosType SurfaceBrightness(const PosType *x) const;
 void printSource();
 
private:
 inline void updateRadius()
 {
  // radius in Source
  // approximation of the radius that encloses 99% of the flux
  setRadius((3.73 - 0.926*index + 1.164*index*index)*Reff);
 }
 
 void assignParams(InputParams& params);
                      
 PosType Reff;  // in radians
 //PosType mag;
 PosType PA;
 PosType index;
 PosType bn;
 PosType q;
 //PosType flux;
 PosType I_r, I_n, I_q;
  PosType cosPA,sinPA;
};
 
#endif /* GALAXIES_SERSIC_H_ */