JOSCPROB::JOscillogram Struct Reference

Auxiliary class for creating oscillograms. More...

#include <JOscillogram.hh>

Public Member Functions
	JOscillogram (const JOscillogramAxis &abscissa, const JOscillogramAxis &ordinate, const JOscChannel &channel, const JOscProbInterpolatorInterface *pInterpolator)
	Constructor. More...
	JOscillogram (const std::string &abscissaName, const JGrid< double > abscissaBinning, const std::string &ordinateName, const JGrid< double > ordinateBinning, const JOscChannel &channel, const JOscProbInterpolatorInterface *pInterpolator)
	Constructor. More...
double	getEnergy (const int i, const int j) const
	Get energy corrresponding to the given bin indices. More...
double	getCosth (const int i, const int j) const
	Get cosine zenith angle corrresponding to the given bin indices. More...
double	getP (const int i, const int j) const
	Get oscillation probability for given bin indices. More...

Private Member Functions
std::pair< double, double >	getTransformation (const int i, const int j) const
	Get energy and cosine zenith angle corresponding to the given bin indices. More...

Private Attributes
JOscillogramAxis	abscissa
	Abscissa axis. More...
JOscillogramAxis	ordinate
	Ordinate axis. More...
JOscChannel	channel
	Oscillation channel. More...
const JOscProbInterpolatorInterface *	pInterpolator
	Pointer to oscillation probability interpolator. More...

Detailed Description

Auxiliary class for creating oscillograms.

Definition at line 78 of file JOscillogram.hh.

Constructor & Destructor Documentation

JOSCPROB::JOscillogram::JOscillogram ( const JOscillogramAxis &  abscissa, const JOscillogramAxis &  ordinate, const JOscChannel &  channel, const JOscProbInterpolatorInterface *  pInterpolator  )

inline

Constructor.

Parameters

abscissa	oscillogram abscissa axis definition
ordinate	oscillogram ordinate axis definition
channel	oscillation channel
pInterpolator	pointer to oscillation probability interpolator

Definition at line 88 of file JOscillogram.hh.

                                                                     :
      abscissa     (abscissa),
      ordinate     (ordinate),
      channel      (channel),
      pInterpolator(pInterpolator)
    {
      if (pInterpolator == NULL) {
        THROW(JNullPointerException, "JOscillogram::JOscillogram(...): Oscillation probability interpolator is not set");
      }
    }

JOSCPROB::JOscillogram::JOscillogram ( const std::string &  abscissaName, const JGrid< double >  abscissaBinning, const std::string &  ordinateName, const JGrid< double >  ordinateBinning, const JOscChannel &  channel, const JOscProbInterpolatorInterface *  pInterpolator  )

inline

Constructor.

Parameters

abscissaName	oscillogram abscissa variable name
abscissaBinning	oscillogram abscissa binning
ordinateName	oscillogram ordinate variable name
ordinateBinning	oscillogram ordinate binning
channel	oscillation channel
pInterpolator	pointer to oscillation probability interpolator

Definition at line 113 of file JOscillogram.hh.

                                                                     :
      abscissa     (JOscVars::getType(abscissaName), abscissaBinning),
      ordinate     (JOscVars::getType(ordinateName), ordinateBinning),
      channel      (channel),
      pInterpolator(pInterpolator)
    {
      if (pInterpolator == NULL) {
        THROW(JNullPointerException, "JOscillogram::JOscillogram(...): Oscillation probability interpolator is not set");
      }
    }

Member Function Documentation

double JOSCPROB::JOscillogram::getEnergy ( const int  i, const int  j  ) const

inline

Get energy corrresponding to the given bin indices.

Parameters

i	abscissa bin index
j	ordinate bin index

Returns

energy [GeV]

Definition at line 137 of file JOscillogram.hh.

    {
      const std::pair<double, double>& result = getTransformation(i, j);

      return result.first;
    }

double JOSCPROB::JOscillogram::getCosth ( const int  i, const int  j  ) const

inline

Get cosine zenith angle corrresponding to the given bin indices.

Parameters

i	abscissa bin index
j	ordinate bin index

Returns

energy [GeV]

Definition at line 152 of file JOscillogram.hh.

    {
      const std::pair<double, double>& result = getTransformation(i, j);

      return result.second;
    }

double JOSCPROB::JOscillogram::getP ( const int  i, const int  j  ) const

inline

Get oscillation probability for given bin indices.

Parameters

i	abscissa bin index
j	ordinate bin index

Returns

oscillation probability

Definition at line 167 of file JOscillogram.hh.

    {
      const double& energy = getEnergy(i, j);
      const double& costh  = getCosth (i, j);
      
      return (*pInterpolator)(channel, energy, costh);
    }

std::pair<double, double> JOSCPROB::JOscillogram::getTransformation ( const int  i, const int  j  ) const

inlineprivate

Get energy and cosine zenith angle corresponding to the given bin indices.

Parameters

i	abscissa bin index
j	ordinate bin index

Returns

transformed coordinate (E [GeV], costh)

Definition at line 186 of file JOscillogram.hh.

    {
      using namespace std;
      
      static const JBaselineCalculator& baselineCalculator = pInterpolator->getBaselineCalculator();
      
      static int ix = -1;
      static int iy = -1;

      static pair<double, double> result = {0.0, 0.0}; // Cache result

      if (i != ix || j != iy) {

        ix = i;
        iy = j;
        
        const double x = abscissa.getX(i);
        const double y = ordinate.getX(j);

        switch (ordinate.type) {
        case (int) JOscVars::COSTH:
          result.second = y;
          break;
        case (int) JOscVars::SINTH:
          result.second = sqrt((1 + y) * (1 - y));
          break;
        case (int) JOscVars::BASELINE:
          result.second = baselineCalculator.getCosth(y);
          break;
        default:
          THROW(JValueOutOfRange, "JOscillogram::getCosth(const double): Invalid ordinate type " << ordinate.type);
        }

        switch (abscissa.type) {
        case (int) JOscVars::ENERGY:
          result.first = x;
          break;
        case (int) JOscVars::LOG10E:
          result.first = pow(10.0, x);
          break;
        case (int) JOscVars::LOE:
          result.first = (x > 0.0 ? baselineCalculator.getBaseline(result.second) / x : 0.0);
          break;
        default:
          THROW(JValueOutOfRange, "JOscillogram::getEnergy(const double, const double): Invalid abscissa type " << abscissa.type);      
        }
      }

      return result;
    }

Member Data Documentation

JOscillogramAxis JOSCPROB::JOscillogram::abscissa

private

Abscissa axis.

Definition at line 238 of file JOscillogram.hh.

JOscillogramAxis JOSCPROB::JOscillogram::ordinate

private

Ordinate axis.

Definition at line 239 of file JOscillogram.hh.

JOscChannel JOSCPROB::JOscillogram::channel

private

Oscillation channel.

Definition at line 241 of file JOscillogram.hh.

const JOscProbInterpolatorInterface* JOSCPROB::JOscillogram::pInterpolator

private

Pointer to oscillation probability interpolator.

Definition at line 243 of file JOscillogram.hh.

---

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

• JOscillogram.hh