ShowerSRecData_ROOT.cc

Go to the documentation of this file.

#include <evt/ShowerScintillatorRecData.h>

#include <fwk/LocalCoordinateSystem.h>
#include <utl/LameShadowPtr_imp.h>
#include <io/ShowerSRecData_ROOT.h>
#include <io/TabulatedFunctionErrors_ROOT.h>
#include <io/RiseTime1000_ROOT.h>
#include <io/SdFootprintData_ROOT.h>
#include <io/StreamerUtilities.h>

using namespace io;


ClassImp(ShowerSRecData_ROOT)


ShowerSRecData_ROOT&
ShowerSRecData_ROOT::operator=(const ShowerSRecData_ROOT& s)
{
  if (this != &s) {
    DeepCopy(fLDF, s.fLDF);
    fBarycenter = s.fBarycenter;
    fBarytime = s.fBarytime;
    fCore = s.fCore;
    fCoreError = s.fCoreError;
    fCorrelationXY = s.fCorrelationXY;
    fCoreTime = s.fCoreTime;
    fCoreTimeError = s.fCoreTimeError;
    fAxis = s.fAxis;
    fBeta = s.fBeta;
    fBetaError = s.fBetaError;
    fBetaSystematics = s.fBetaSystematics;
    fGamma = s.fGamma;
    fGammaError = s.fGammaError;
    fNKGFermiMu = s.fNKGFermiMu;
    fNKGFermiTau = s.fNKGFermiTau;
    fCurvature = s.fCurvature;
    fCurvatureError = s.fCurvatureError;
    fCurvatureTimeOffset = s.fCurvatureTimeOffset;
    fCurvatureTimeOffsetError = s.fCurvatureTimeOffsetError;
    fShowerSizeType = s.fShowerSizeType;
    fShowerSize = s.fShowerSize;
    fShowerSizeError = s.fShowerSizeError;
    fShowerSizeSystematics = s.fShowerSizeSystematics;
    fShowerSizeGeomagneticCorr = s.fShowerSizeGeomagneticCorr;
    fShowerSizeGeomagneticCorrError = s.fShowerSizeGeomagneticCorrError;
    fShowerSizeAtmosphericCorr = s.fShowerSizeAtmosphericCorr;
    fShowerSizeAtmosphericCorrError = s.fShowerSizeAtmosphericCorrError;
    fLDFOptimumDistance = s.fLDFOptimumDistance;
    fLDFChi2 = s.fLDFChi2;
    fLDFNdof = s.fLDFNdof;
    fThetaError = s.fThetaError;
    fPhiError = s.fPhiError;
    fCorrelationThetaPhi = s.fCorrelationThetaPhi;
    fAngleChi2 = s.fAngleChi2;
    fAngleNdof = s.fAngleNdof;
    fTimeResidualMean = s.fTimeResidualMean;
    fTimeResidualSpread = s.fTimeResidualSpread;
    fEnergy = s.fEnergy;
    fEnergyError = s.fEnergyError;
    fEnergyLDFSystematics = s.fEnergyLDFSystematics;
    fNumberOfActiveStations = s.fNumberOfActiveStations;
    fT4Trigger = s.fT4Trigger;
    fT5Trigger = s.fT5Trigger;
    fT5HottestStation = s.fT5HottestStation;
    fT5PriorActiveNeighbors = s.fT5PriorActiveNeighbors;
    fT5PosteriorActiveNeighbors = s.fT5PosteriorActiveNeighbors;
    fT5PosteriorCoreTriangle = s.fT5PosteriorCoreTriangle;
    fIsTOTdMoPSSilent = s.fIsTOTdMoPSSilent;
    std::copy(s.fReconstructionSeed, s.fReconstructionSeed+3, fReconstructionSeed);
    fLDFRecStage = s.fLDFRecStage;
    fSeedAxis = s.fSeedAxis;
    fSeedBarycenter = s.fSeedBarycenter;
    fBadPeriodId = s.fBadPeriodId;
    DeepCopy(fRiseTime1000, s.fRiseTime1000);
    DeepCopy(fFootprintData, s.fFootprintData);
  }
  return *this;
}


ShowerSRecData_ROOT::~ShowerSRecData_ROOT()
{
  delete fLDF;
  delete fRiseTime1000;
  delete fFootprintData;
}


ShowerSRecData_ROOT::ShowerSRecData_ROOT(const evt::ShowerSRecData& s) :
  fBarycenter(s.GetBarycenter()),
  fBarytime(s.GetBarytime()),
  fCore(s.GetCorePosition()),
  fCoreError(s.GetCoreError()),
  fCorrelationXY(s.GetCorrelationXY()),
  fCoreTime(s.GetCoreTime()),
  fCoreTimeError(s.GetCoreTimeError().GetInterval()),
  fAxis(s.GetAxis()),
  fBeta(s.GetBeta()),
  fBetaError(s.GetBetaError()),
  fBetaSystematics(s.GetBetaSystematics()),
  fGamma(s.GetGamma()),
  fGammaError(s.GetGammaError()),
  fNKGFermiMu(s.GetNKGFermiMu()),
  fNKGFermiTau(s.GetNKGFermiTau()),
  fCurvature(s.GetCurvature()),
  fCurvatureError(s.GetCurvatureError()),
  fCurvatureTimeOffset(s.GetCurvatureTimeOffset()),
  fCurvatureTimeOffsetError(s.GetCurvatureTimeOffsetError()),
  fShowerSizeType(s.GetShowerSizeType()),
  fShowerSize(s.GetShowerSize()),
  fShowerSizeError(s.GetShowerSizeError()),
  fShowerSizeSystematics(s.GetShowerSizeSystematics()),
  fShowerSizeGeomagneticCorr(s.GetSizeGeomagneticCorr()),
  fShowerSizeGeomagneticCorrError(s.GetSizeGeomagneticCorrError()),
  fShowerSizeAtmosphericCorr(s.GetSizeAtmosphericCorr()),
  fShowerSizeAtmosphericCorrError(s.GetSizeAtmosphericCorrError()),
  fLDFOptimumDistance(s.GetLDFOptimumDistance()),
  fLDFChi2(s.GetLDFChi2()),
  fLDFNdof(s.GetLDFNdof()),
  fLDFLikelihood(s.GetLDFLikelihood()),
  fThetaError(s.GetThetaError()),
  fPhiError(s.GetPhiError()),
  fCorrelationThetaPhi(s.GetCorrelationThetaPhi()),
  fAngleChi2(s.GetAngleChi2()),
  fAngleNdof(s.GetAngleNdof()),
  fTimeResidualMean(s.GetTimeResidualMean()),
  fTimeResidualSpread(s.GetTimeResidualSpread()),
  fEnergy(s.GetEnergy()),
  fEnergyError(s.GetEnergyError()),
  fEnergyLDFSystematics(s.GetEnergyLDFSystematics()),
  fNumberOfActiveStations(s.GetNumberOfActiveStations()),
  fT4Trigger(s.GetT4Trigger()),
  fT5Trigger(s.GetT5Trigger()),
  fT5HottestStation(s.GetT5HottestStation()),
  fT5PriorActiveNeighbors(s.GetT5PriorActiveNeighbors()),
  fT5PosteriorActiveNeighbors(s.GetT5PosteriorActiveNeighbors()),
  fT5PosteriorCoreTriangle(s.GetT5PosteriorCoreTriangle()),
  fIsTOTdMoPSSilent(s.IsTOTdMoPSSilent()),
  fLDFRecStage(s.GetLDFRecStage()),
  fSeedAxis(s.GetSeedAxis()),
  fSeedBarycenter(s.GetSeedBarycenter()),
  fBadPeriodId(s.GetBadPeriodId())
{
  if (s.HasLDF())
    ConvertCopy(fLDF, s.GetLDF());

  const std::vector<int>& seed = s.GetReconstructionSeed();
  if (seed.size() == 3)
    std::copy(seed.begin(), seed.end(), fReconstructionSeed);

  if (s.HasRiseTime1000())
    ConvertCopy(fRiseTime1000, s.GetRiseTime1000());
  if (s.HasFootprintData())
    ConvertCopy(fFootprintData, s.GetFootprintData());
}


void
ShowerSRecData_ROOT::operator>>(evt::ShowerSRecData& s)
  const
{
  if (fLDF) {
    if (!s.HasLDF())
      s.MakeLDF();
    *fLDF >> s.GetLDF();
  }

  s.SetBarycenter(fBarycenter.Convert());
  s.SetBarytime(fBarytime.Convert());

  s.SetCorePosition(fCore.Convert());
  s.SetCoreError(fCoreError.Convert());
  s.SetCorrelationXY(fCorrelationXY);
  s.SetCoreTime(fCoreTime.Convert(), utl::TimeInterval(fCoreTimeError));
  s.SetAxis(fAxis.Convert());

  s.SetBeta(fBeta, fBetaError);
  s.SetBetaSystematics(fBetaSystematics);
  s.SetGamma(fGamma, fGammaError);
  s.SetNKGFermiParameters(fNKGFermiMu, fNKGFermiTau);
  s.SetCurvature(fCurvature, fCurvatureError);
  s.SetCurvatureTimeOffset(fCurvatureTimeOffset, fCurvatureTimeOffsetError);
  s.SetShowerSizeType(static_cast<evt::ShowerSRecData::ShowerSizeType>(fShowerSizeType));
  s.SetShowerSize(fShowerSize, fShowerSizeError);
  s.SetShowerSizeSystematics(fShowerSizeSystematics);
  s.SetSizeGeomagneticCorr(fShowerSizeGeomagneticCorr, fShowerSizeGeomagneticCorrError);
  s.SetSizeAtmosphericCorr(fShowerSizeAtmosphericCorr, fShowerSizeAtmosphericCorrError);

  s.SetLDFOptimumDistance(fLDFOptimumDistance);
  s.SetLDFChi2(fLDFChi2, fLDFNdof);
  s.SetLDFLikelihood(fLDFLikelihood);

  s.SetThetaError(fThetaError);
  s.SetPhiError(fPhiError);
  s.SetCorrelationThetaPhi(fCorrelationThetaPhi);
  s.SetAngleChi2(fAngleChi2, fAngleNdof);
  s.SetTimeResidualMean(fTimeResidualMean);
  s.SetTimeResidualSpread(fTimeResidualSpread);
  s.SetEnergy(fEnergy, fEnergyError);
  s.SetEnergyLDFSystematics(fEnergyLDFSystematics);

  s.SetNumberOfActiveStations(fNumberOfActiveStations);

  s.SetT4Trigger(fT4Trigger);
  s.SetT5Trigger(fT5Trigger);
  s.SetT5HottestStation(fT5HottestStation);
  s.SetT5PriorActiveNeighbors(fT5PriorActiveNeighbors);
  s.SetT5PosteriorActiveNeighbors(fT5PosteriorActiveNeighbors);
  s.SetIsTOTdMoPSSilent(fIsTOTdMoPSSilent);

  if (fReconstructionSeed[0])
    s.SetReconstructionSeed(std::vector<int>(fReconstructionSeed, fReconstructionSeed+3));

  s.SetLDFRecStage(fLDFRecStage);

  s.SetSeedAxis(fSeedAxis.Convert());
  s.SetSeedBarycenter(fSeedBarycenter.Convert());

  s.SetBadPeriodId(fBadPeriodId);

  if (fRiseTime1000) {
    if (!s.HasRiseTime1000())
      s.MakeRiseTime1000();
    *fRiseTime1000 >> s.GetRiseTime1000();
  }

  if (fFootprintData) {
    if (!s.HasFootprintData())
      s.MakeFootprintData();
    *fFootprintData >> s.GetFootprintData();
  }
}