ShowerFRecData_ROOT.cc

Go to the documentation of this file.

#include <io/ShowerFRecData_ROOT.h>
#include <io/StreamerUtilities.h>

using namespace io;


ClassImp(ShowerFRecData_ROOT)


ShowerFRecData_ROOT::ShowerFRecData_ROOT(const evt::ShowerFRecData& sh) :
  fEmEnergy(sh.GetEmEnergy()),
  fEmEnergyError(sh.GetEmEnergyError()),
  fTotalEnergy(sh.GetTotalEnergy()),
  fTotalEnergyError(sh.GetTotalEnergyError()),
  fEnergyCutoff(sh.GetEnergyCutoff()),
  fEmEnergyErrors(sh.GetEmEnergyErrorsBegin(), sh.GetEmEnergyErrorsEnd()),
  fTotalEnergyErrors(sh.GetTotalEnergyErrorsBegin(), sh.GetTotalEnergyErrorsEnd()),
  fXmaxErrors(sh.GetXmaxErrorsBegin(), sh.GetXmaxErrorsEnd()),
  fAxis(sh.GetAxis()),
  fAxisError(sh.GetAxisError()),
  fCore(sh.GetCorePosition()),
  fCoreError(sh.GetCoreError()),
  fCoreTime(sh.GetCoreTime()),
  fCoreTimeError(sh.GetCoreTimeError().GetInterval()),
  fCorrelationXY(sh.GetCorrelationXY()),
  fCorrelationThetaPhi(sh.GetCorrelationThetaPhi()),
  fSDTimeResidual(sh.GetSDTimeResidual()),
  fStationIds(sh.GetStationIds())
{
  if (sh.HasEnergyDeposit())
    fEnergyDeposit = sh.GetEnergyDeposit();

  if (sh.HasLongitudinalProfile())
    fLongitudinalProfile = sh.GetLongitudinalProfile();

  if (sh.HasGHParameters())
    fGHParameters = sh.GetGHParameters();

  if (sh.HasMultipleGHParameters())
    fMultipleGHParameters = sh.GetMultipleGHParameters();

  if (sh.HasFluorescencePhotons())
    fFluorescencePhotons = sh.GetFluorescencePhotons();

  if (sh.HasCherenkovPhotons())
    fCherenkovPhotons = sh.GetCherenkovPhotons();
}


void
ShowerFRecData_ROOT::operator>>(evt::ShowerFRecData& sh)
  const
{
  sh.SetEmEnergy(fEmEnergy, fEmEnergyError);
  sh.SetTotalEnergy(fTotalEnergy, fTotalEnergyError);
  sh.SetEnergyCutoff(fEnergyCutoff);

  for (const auto& te : fEmEnergyErrors)
    sh.SetEmEnergyError(te.second, evt::ShowerFRecData::EUncertaintyType(te.first));
  for (const auto& te : fTotalEnergyErrors)
    sh.SetTotalEnergyError(te.second, evt::ShowerFRecData::EUncertaintyType(te.first));
  for (const auto& te : fXmaxErrors)
    sh.SetXmaxError(te.second, evt::ShowerFRecData::EUncertaintyType(te.first));

  sh.SetAxis(fAxis.Convert());
  sh.SetAxisError(fAxisError.Convert());
  sh.SetCorePosition(fCore.Convert());
  sh.SetCoreError(fCoreError.Convert());
  sh.SetCoreTime(fCoreTime.Convert(), utl::TimeInterval(fCoreTimeError));

  if (fEnergyDeposit) {
    sh.MakeEnergyDeposit();
    fEnergyDeposit >> sh.GetEnergyDeposit();
  }

  if (fLongitudinalProfile) {
    sh.MakeLongitudinalProfile();
    fLongitudinalProfile >> sh.GetLongitudinalProfile();
  }

  if (fFluorescencePhotons) {
    sh.MakeFluorescencePhotons();
    fFluorescencePhotons >> sh.GetFluorescencePhotons();
  }

  if (fCherenkovPhotons) {
    sh.MakeCherenkovPhotons();
    fCherenkovPhotons >> sh.GetCherenkovPhotons();
  }

  if (fGHParameters) {
    evt::VGaisserHillasParameter* gh = 0;
    fGHParameters >> gh;
    sh.MakeGHParameters(*gh);
  }

  if (fMultipleGHParameters)
    fMultipleGHParameters >> sh.GetMultipleGHParameters();

  sh.SetCorrelationXY(fCorrelationXY);
  sh.SetCorrelationThetaPhi(fCorrelationThetaPhi);
  sh.SetSDTimeResidual(fSDTimeResidual);

  for (const auto& id : fStationIds)
    sh.AddStationId(id);
}