GaisserHillas4Parameter.h

Go to the documentation of this file.

#ifndef _evt_GaisserHillas4Parameter_h_
#define _evt_GaisserHillas4Parameter_h_

#include <evt/VGaisserHillasParameter.h>
#include <evt/GaisserHillasTypes.h>
#include <utl/TabulatedFunction.h>
#include <iostream>
#include <vector>
#include <utility>


namespace io {
  class GaisserHillasParameters_ROOT;
}

namespace evt {

  class GaisserHillas4Parameter : public VGaisserHillasParameter {

  public:
    typedef std::pair<gh::EShapeParameter, double> ShapeParameter;

    GaisserHillas4Parameter(const gh::EFunctionType functionType = gh::eClassic);
    GaisserHillas4Parameter(const double nmax, const double xmax,
                            const ShapeParameter par1, ShapeParameter par2,
                            const gh::EFunctionType functionType = gh::eClassic);
    virtual ~GaisserHillas4Parameter() { }

    GaisserHillas4Parameter* Create() const
    { return new GaisserHillas4Parameter; }

    GaisserHillas4Parameter* Clone() const
    { return new GaisserHillas4Parameter(*this); }

    gh::EFunctionType GetFunctionType() const { return fFunctionType; }

    double GetShapeParameter(const gh::EShapeParameter par) const;
    double GetShapeParameterError(const gh::EShapeParameter par) const;

    double GetCorrelationXMaxShapeParameter(const gh::EShapeParameter par) const;
    double GetCorrelationNMaxShapeParameter(const gh::EShapeParameter par) const;
    double GetCorrelationShapeParameters() const { return fRhoShapePar1Par2; }

    bool IsInternal(const gh::EShapeParameter par) const;

    double Eval(const double depth) const;
    void Dump(std::ostream& os = std::cout) const;
    double GetIntegral() const;
    double GetIntegralError() const;

    double InverseLeft(const double h) const;
    double InverseRight(const double h) const;


    void SetShapeParameter(const gh::EShapeParameter par,
                           const double value, const double error);
    void SetCorrelationXMaxShapeParameter(const gh::EShapeParameter par,
                                          const double rho);
    void SetCorrelationNMaxShapeParameter(const gh::EShapeParameter par,
                                          const double rho);
    void SetCorrelationShapeParameters(const double rho) { fRhoShapePar1Par2 = rho; }

    // \\deprecated
    double GetXZero() const;
    double GetXZeroError() const;
    double GetLambda() const;
    double GetLambdaError() const;
    double GetNMaxLambdaCorrelation() const;
    double GetNMaxX0Correlation() const;
    double GetXMaxLambdaCorrelation() const;
    double GetXMaxX0Correlation() const;
    double GetLambdaX0Correlation() const;

  private:
    double Inverse(const double h, const int branch) const;
    enum EInternalShapeParameter {
      eFirst = 0,
      eLast
    };

    EInternalShapeParameter ExternalToInternal(const gh::EShapeParameter par) const;
    gh::EShapeParameter InternalToExternal(const EInternalShapeParameter) const;
    static gh::EShapeParameter InternalToExternal(const gh::EFunctionType,
                                                  const EInternalShapeParameter);

    double CalculateR(const double asym) const;

    gh::EFunctionType fFunctionType;

    double fShapePar[eLast+1];
    double fShapeParError[eLast+1];
    double fRhoNMaxShapePar[eLast+1];
    double fRhoXMaxShapePar[eLast+1];
    double fRhoShapePar1Par2;

    static utl::TabulatedFunction fRvsAsymTable;

    friend class io::GaisserHillasParameters_ROOT;

  };

}


#endif