Station_ROOT.cc

Go to the documentation of this file.

#include <io/Station_ROOT.h>
#include <io/StationSimData_ROOT.h>
#include <io/StationRecData_ROOT.h>
#include <io/StationCalibData_ROOT.h>
#include <io/StationTriggerData_ROOT.h>
#include <io/StationGPSData_ROOT.h>
#include <io/sevtScintillator_ROOT.h>
#include <io/SmallPMTData_ROOT.h>

using namespace io;


ClassImp(Station_ROOT)


Station_ROOT::~Station_ROOT()
{
  delete fSimData;
  delete fRecData;
  delete fCalibData;
  delete fTriggerData;
  delete fGPSData;
  delete fVemTrace;
  delete fScintillator;
  delete fSmallPMTData;
}


Station_ROOT&
Station_ROOT::operator=(const Station_ROOT& s)
{
  if (this != &s) {
    fId = s.fId;
    fPMTs = s.fPMTs;
    DeepCopy(fSimData, s.fSimData);
    DeepCopy(fRecData, s.fRecData);
    DeepCopy(fCalibData, s.fCalibData);
    DeepCopy(fTriggerData, s.fTriggerData);
    DeepCopy(fGPSData, s.fGPSData);
    fTraceStart = s.fTraceStart;
    DeepCopy(fVemTrace, s.fVemTrace);
    fReconstructionStatus = s.fReconstructionStatus;
    fSaturation[0] = s.fSaturation[0];
    fSaturation[1] = s.fSaturation[1];
    fSmallPMTSaturation = s.fSmallPMTSaturation;
    fRejectionStatus = s.fRejectionStatus;
    fT2Life = s.fT2Life;
    DeepCopy(fScintillator, s.fScintillator);
    DeepCopy(fSmallPMTData, s.fSmallPMTData);
  }
  return *this;
}


Station_ROOT::Station_ROOT(const sevt::Station& s) :
  fId(s.GetId()),
  fPMTs(s.PMTsBegin(sdet::PMTConstants::eAnyType), s.PMTsEnd(sdet::PMTConstants::eAnyType)),
  fTraceStart(s.GetTraceStartTime()),
  fSaturation{s.IsLowGainSaturation(), s.IsHighGainSaturation()},
  fSmallPMTSaturation(s.IsSmallPMTSaturation()),
  fRejectionStatus(s.GetRejectionStatus()),
  fT2Life(s.GetT2Life())
{
  if (s.HasSimData())
    ConvertCopy(fSimData, s.GetSimData());

  if (s.HasRecData())
    ConvertCopy(fRecData, s.GetRecData());

  if (s.HasCalibData())
    ConvertCopy(fCalibData, s.GetCalibData());

  if (s.HasTriggerData())
    ConvertCopy(fTriggerData, s.GetTriggerData());

  if (s.HasGPSData())
    ConvertCopy(fGPSData, s.GetGPSData());

  if (s.HasVEMTrace())
    fVemTrace = new MultiTraceD_ROOT(s.VEMTracesBegin(), s.VEMTracesEnd());

  if (s.IsSilent())
    fReconstructionStatus = sevt::StationConstants::eSilent;
  else if (s.IsRejected())
    fReconstructionStatus = sevt::StationConstants::eRejected;
  else if (s.IsCandidate())
    fReconstructionStatus = sevt::StationConstants::eCandidate;

  if (s.HasScintillator())
    ConvertCopy(fScintillator, s.GetScintillator());

  if (s.HasSmallPMTData())
    ConvertCopy(fSmallPMTData, s.GetSmallPMTData());
}


void
Station_ROOT::operator>>(sevt::Station& s)
  const
{
  if (fId != s.GetId())
    throw utl::IOFailureException("Station id in source and destination are not equal!");

  for (const auto& p : fPMTs) {
     const int pid = p.fId;
     p >> s.GetPMT(pid);
  }

  if (fSimData) {
    s.MakeSimData();
    *fSimData >> s.GetSimData();
  }

  if (fRecData) {
    s.MakeRecData();
    *fRecData >> s.GetRecData();
  }

  if (fCalibData) {
    s.MakeCalibData();
    *fCalibData >> s.GetCalibData();
  }

  if (fTriggerData) {
    s.MakeTriggerData();
    *fTriggerData >> s.GetTriggerData();
  }

  if (fGPSData) {
    s.MakeGPSData();
    *fGPSData >> s.GetGPSData();
  }

  utl::TimeStamp ts;
  fTraceStart >> ts;
  s.SetTraceStartTime(ts);

  if (fVemTrace) {
    for (const auto& st : fVemTrace->fTraces) {
      const auto source = static_cast<sevt::StationConstants::SignalComponent>(st.first);
      if (!s.HasVEMTrace(source))
        s.MakeVEMTrace(source);
      st.second >> s.GetVEMTrace(source);
    }
  }

  if (fReconstructionStatus == sevt::StationConstants::eSilent)
    s.SetSilent();
  else if (fReconstructionStatus == sevt::StationConstants::eRejected)
    s.SetRejected(fRejectionStatus);
  else if (fReconstructionStatus == sevt::StationConstants::eCandidate)
    s.SetCandidate();

  s.SetLowGainSaturation(fSaturation[0]);
  s.SetHighGainSaturation(fSaturation[1]);
  s.SetSmallPMTSaturation(fSmallPMTSaturation);

  s.SetT2Life(fT2Life);

  if (fScintillator)
    *fScintillator >> s.GetScintillator();

  if (fSmallPMTData)
    *fSmallPMTData >> s.GetSmallPMTData();
}