Implementation of histogram source morphology. More...

#include <JMorphology.hh>

Inheritance diagram for JASTRONOMY::JMorphologyHistogram:

Public Member Functions
virtual	~JMorphologyHistogram ()
	Virtual destructor.

virtual std::istream &	read (std::istream &in) override
	Read source morphology from input stream.

virtual std::ostream &	write (std::ostream &out) const override
	Write source morphology to output stream.

virtual JSourceLocation	get () const override
	Get location on sky.

const JSourceLocation &	getSourceLocation () const
	Get source location.

	operator JAngle3D () const
	Type conversion operator.

double	getDeclination () const
	Get declination.

double	getRightAscension () const
	Get right ascension.

double	getDot (const JSourceLocation &location) const
	Dot product.

void	set (const angle_type_deg &angle)
	Convert angle.

bool	equals (const angle_type &angle, const double precision=std::numeric_limits< double >::min()) const
	Check equality.

Static Public Member Functions
static std::istream &	read (std::istream &in, JSourceLocation &source)
	Read source location in degrees from input stream.

static std::ostream &	write (std::ostream &out, const JSourceLocation &source)
	Write source location in degrees to output stream.

Public Attributes
std::string	filename

std::string	histname

Protected Member Functions
void	load ()
	Load histogram from file.

Protected Attributes
TFile *	in = NULL

TH2D *	h2 = NULL

double	_theta_

double	_phi_

Detailed Description

Implementation of histogram source morphology.

Note that all input angles are in degrees.

Definition at line 383 of file JMorphology.hh.

Constructor & Destructor Documentation

◆ ~JMorphologyHistogram()

virtual JASTRONOMY::JMorphologyHistogram::~JMorphologyHistogram ( )

inlinevirtual

Virtual destructor.

Definition at line 389 of file JMorphology.hh.

    {
      if (in != NULL) {
        in->Close();
      }
    }

Member Function Documentation

◆ read() [1/2]

virtual std::istream & JASTRONOMY::JMorphologyHistogram::read ( std::istream & in )

inlineoverridevirtual

Read source morphology from input stream.

Parameters

in	input stream

Returns: input stream

Implements JASTRONOMY::JMorphology.

Definition at line 403 of file JMorphology.hh.

    {
      JMorphology::read(in, *this)
        >> this->filename
        >> this->histname;
 
      load();
 
      return in;
    }

◆ write() [1/2]

virtual std::ostream & JASTRONOMY::JMorphologyHistogram::write ( std::ostream & out ) const

inlineoverridevirtual

Write source morphology to output stream.

Parameters

out	output stream

Returns: output stream

Implements JASTRONOMY::JMorphology.

Definition at line 421 of file JMorphology.hh.

    {
      return JMorphology::write(out, *this)
        << ' ' << this->filename
        << ' ' << this->histname;
    }

◆ get()

virtual JSourceLocation JASTRONOMY::JMorphologyHistogram::get ( ) const

inlineoverridevirtual

Get location on sky.

Returns: location on sky

Implements JASTRONOMY::JMorphology.

Definition at line 434 of file JMorphology.hh.

    {
      using namespace JPP;
 
      Double_t x;
      Double_t y;
 
      h2->GetRandom2(x, y, gRandom);
 
      x = sin(getRadians(x));
      y = sin(getRadians(y));
        
      JDirection3D u(x, y, sqrt(1.0 - x*x - y*y));
 
      const JRotation3D Rs(this->getSourceLocation());
 
      u.rotate_back(Rs);
 
      return JSourceLocation(u);
    }

◆ load()

void JASTRONOMY::JMorphologyHistogram::load ( )

inlineprotected

Load histogram from file.

Definition at line 462 of file JMorphology.hh.

    {
      in = TFile::Open(filename.c_str(), "exist");
 
      if (in == NULL || !in->IsOpen()) {
        THROW(JFileOpenException, "File: " << filename << " not opened.");
      }
 
      h2 = dynamic_cast<TH2D*>(in->Get(histname.c_str()));
 
      if (h2 == NULL) {
        THROW(JValueOutOfRange, "Histogram: " << histname << " not found.");
      }
    }

◆ read() [2/2]

static std::istream & JASTRONOMY::JMorphology::read	(	std::istream &	in,
		JSourceLocation &	source )

inlinestaticinherited

Read source location in degrees from input stream.

Parameters

in	input stream
source	source

Returns: input stream

Definition at line 62 of file JMorphology.hh.

    {
      angle_type_deg angle;
 
      in >> angle;
 
      source.set(angle);
 
      return in;
    }

◆ write() [2/2]

static std::ostream & JASTRONOMY::JMorphology::write	(	std::ostream &	out,
		const JSourceLocation &	source )

inlinestaticinherited

Write source location in degrees to output stream.

Parameters

out	output stream
source	source

Returns: output stream

Definition at line 81 of file JMorphology.hh.

    {
      return out << angle_type_deg(source);
    }

◆ getSourceLocation()

const JSourceLocation & JASTRONOMY::JSourceLocation::getSourceLocation ( ) const

inlineinherited

Get source location.

Returns: source location

Definition at line 485 of file JAstronomy.hh.

    {
      return static_cast<const JSourceLocation&>(*this);
    }

◆ operator JAngle3D()

JASTRONOMY::JSourceLocation::operator JAngle3D ( ) const

inlineinherited

Type conversion operator.

Returns: polar angles [rad]

Definition at line 496 of file JAstronomy.hh.

    {
      return JAngle3D(PI/2 - _theta_, _phi_);
    }

◆ getDeclination()

double JASTRONOMY::JSourceLocation::getDeclination ( ) const

inlineinherited

Get declination.

Definition at line 502 of file JAstronomy.hh.

◆ getRightAscension()

double JASTRONOMY::JSourceLocation::getRightAscension ( ) const

inlineinherited

Get right ascension.

Definition at line 503 of file JAstronomy.hh.

◆ getDot()

double JASTRONOMY::JSourceLocation::getDot ( const JSourceLocation & location ) const

inlineinherited

Dot product.

Parameters

location source location

Returns: dot product

Definition at line 512 of file JAstronomy.hh.

    {
      return 
        cos(this->_theta_) * cos(location._theta_) * cos(this->_phi_ - location._phi_) +
        sin(this->_theta_) * sin(location._theta_);
    }

◆ set()

void JASTRONOMY::angle_type_rad::set ( const angle_type_deg & angle )

inlineinherited

Convert angle.

Parameters

angle angle [rad]

Definition at line 293 of file JAstronomy.hh.

    {
      static_cast<angle_type_rad&>(*this) = angle_type_rad(angle);
    }

◆ equals()

bool JASTRONOMY::angle_type::equals	(	const angle_type &	angle,
		const double	precision = std::numeric_limits<double>::min() ) const

inlineinherited

Check equality.

Parameters

angle	pair of angles
precision	precision

Returns: true if angles are equal; else false

Definition at line 194 of file JAstronomy.hh.

    {
      return (fabs(this->_theta_ - angle._theta_) <= precision &&
              fabs(this->_phi_   - angle._phi_)   <= precision);
    }

Member Data Documentation

◆ filename

std::string JASTRONOMY::JMorphologyHistogram::filename

Definition at line 455 of file JMorphology.hh.

◆ histname

std::string JASTRONOMY::JMorphologyHistogram::histname

Definition at line 456 of file JMorphology.hh.

◆ in

TFile* JASTRONOMY::JMorphologyHistogram::in = NULL

mutableprotected

Definition at line 477 of file JMorphology.hh.

◆ h2

TH2D* JASTRONOMY::JMorphologyHistogram::h2 = NULL

mutableprotected

Definition at line 478 of file JMorphology.hh.

◆ _theta_

double JASTRONOMY::angle_type::_theta_

protectedinherited

Definition at line 250 of file JAstronomy.hh.

◆ _phi_

double JASTRONOMY::angle_type::_phi_

protectedinherited

Definition at line 251 of file JAstronomy.hh.

The documentation for this struct was generated from the following file:

JMorphology.hh

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ~JMorphologyHistogram()

Member Function Documentation

◆ read() [1/2]

◆ write() [1/2]

◆ get()

◆ load()

◆ read() [2/2]

◆ write() [2/2]

◆ getSourceLocation()

◆ operator JAngle3D()

◆ getDeclination()

◆ getRightAscension()

◆ getDot()

◆ set()

◆ equals()

Member Data Documentation

◆ filename

◆ histname

◆ in

◆ h2

◆ _theta_

◆ _phi_