GaisserHillasParameters_ROOT.cc

Go to the documentation of this file.

#include <io/GaisserHillasParameters_ROOT.h>
#include <evt/VGaisserHillasParameter.h>
#include <evt/GaisserHillas2Parameter.h>
#include <evt/GaisserHillas4Parameter.h>
#include <evt/GaisserHillas6Parameter.h>
#include <evt/GaisserHillasTypes.h>

using namespace io;


ClassImp(GaisserHillasParameters_ROOT)


GaisserHillasParameters_ROOT::GaisserHillasParameters_ROOT(const evt::VGaisserHillasParameter& g) :
  fXMax(g.GetXMax()),
  fXMaxError(g.GetXMaxError()),
  fNMax(g.GetNMax()),
  fNMaxError(g.GetNMaxError()),
  fRhoNMaxXMax(g.GetNMaxXMaxCorrelation()),
  fChiSqr(g.GetChiSquare()),
  fNdof(g.GetNdof()),
  fGHtype(eUndefined),
  fGHFunctionType(evt::gh::eClassic),
  fIsdEdXProfile(g.IsdEdXProfile())
{
  if (dynamic_cast<const evt::GaisserHillas2Parameter*>(&g)) {
    fGHtype = e2Parameter;
  } else {
    const evt::GaisserHillas4Parameter* const gh4 =
      dynamic_cast<const evt::GaisserHillas4Parameter*>(&g);
    if (gh4) {
      fGHFunctionType = gh4->GetFunctionType();
      fGHtype = e4Parameter;

      const evt::gh::EShapeParameter lambda =
        evt::GaisserHillas4Parameter::InternalToExternal(gh4->GetFunctionType(),
                                                         evt::GaisserHillas4Parameter::eLast);
      const evt::gh::EShapeParameter x0 =
        evt::GaisserHillas4Parameter::InternalToExternal(gh4->GetFunctionType(),
                                                         evt::GaisserHillas4Parameter::eFirst);

      fXZero = gh4->GetShapeParameter(x0);
      fXZeroError = gh4->GetShapeParameterError(x0);
      fA = gh4->GetShapeParameter(lambda);
      fAError = gh4->GetShapeParameterError(lambda);

      fRhoNMaxX0 = gh4->GetCorrelationNMaxShapeParameter(x0);
      fRhoNMaxA = gh4->GetCorrelationNMaxShapeParameter(lambda);
      fRhoXMaxX0 = gh4->GetCorrelationXMaxShapeParameter(x0);
      fRhoXMaxA = gh4->GetCorrelationXMaxShapeParameter(lambda);
      fRhoAX0 = gh4->GetCorrelationShapeParameters();
    } else {
      const evt::GaisserHillas6Parameter* const gh6 =
        dynamic_cast<const evt::GaisserHillas6Parameter*>(&g);
      if (gh6) {
        fGHtype = e6Parameter;
        fXZero = gh6->GetXZero();
        fXZeroError = gh6->GetXZeroError();
        fA = gh6->GetA();
        fAError = gh6->GetAError();
        fB = gh6->GetB();
        fBError = gh6->GetBError();
        fC = gh6->GetC();
        fCError = gh6->GetCError();
      }
    }
  }
}


void
GaisserHillasParameters_ROOT::operator>>(evt::VGaisserHillasParameter*& g)
  const
{
  switch (fGHtype) {
  case e2Parameter:
    g = new evt::GaisserHillas2Parameter;
    break;
  case e4Parameter:
    {
      evt::GaisserHillas4Parameter* const gh4 =
        new evt::GaisserHillas4Parameter(evt::gh::EFunctionType(fGHFunctionType));

      const evt::gh::EShapeParameter lambda =
        evt::GaisserHillas4Parameter::InternalToExternal(gh4->GetFunctionType(),
                                                         evt::GaisserHillas4Parameter::eLast);
      const evt::gh::EShapeParameter x0 =
        evt::GaisserHillas4Parameter::InternalToExternal(gh4->GetFunctionType(),
                                                         evt::GaisserHillas4Parameter::eFirst);
      gh4->SetShapeParameter(x0, fXZero, fXZeroError);
      gh4->SetShapeParameter(lambda, fA, fAError);
      gh4->SetCorrelationNMaxShapeParameter(x0, fRhoNMaxX0);
      gh4->SetCorrelationNMaxShapeParameter(lambda, fRhoNMaxA);
      gh4->SetCorrelationNMaxShapeParameter(x0, fRhoXMaxX0);
      gh4->SetCorrelationNMaxShapeParameter(lambda, fRhoXMaxA);
      gh4->SetCorrelationShapeParameters(fRhoAX0);

      g = gh4;
      break;
    }
  case e6Parameter:
    {
      evt::GaisserHillas6Parameter* const gh6 = new evt::GaisserHillas6Parameter();
      gh6->SetXZero(fXZero, fXZeroError);
      gh6->SetA(fA, fAError);
      gh6->SetB(fB, fBError);
      gh6->SetC(fC, fCError);
      g = gh6;
      break;
    }
  default:
    g = new evt::GaisserHillas2Parameter();
  }

  g->SetXMax(fXMax, fXMaxError);
  g->SetNMax(fNMax, fNMaxError, fIsdEdXProfile);
  g->SetChiSquare(fChiSqr, fNdof);
  g->SetNMaxXMaxCorrelation(fRhoNMaxXMax);
}