ScintillatorSimData_ROOT.cc

Go to the documentation of this file.

#include <io/ScintillatorSimData_ROOT.h>
#include <mevt/ScintillatorSimData.h>
#include <io/StreamerUtilities.h>

using namespace io;


ClassImp(ScintillatorSimData_ROOT)


SPEPulse_ROOT::SPEPulse_ROOT(const mevt::ScintillatorSimData::SPEPulse& p) :
  fMu(p.fMu),
  fSigma(p.fSigma),
  fAmplitude(p.fAmplitude),
  fDestPixelId(p.fDestPixelId)
{ }


SPEPulse_ROOT&
SPEPulse_ROOT::operator>>(mevt::ScintillatorSimData::SPEPulse& p)
{
  p.fMu = fMu;
  p.fSigma = fSigma;
  p.fAmplitude = fAmplitude;
  p.fDestPixelId = fDestPixelId;
  return *this;
}


PEPulse_ROOT::PEPulse_ROOT(const mevt::ScintillatorSimData::PEPulse& p) :
  fT0(p.fT0),
  fAmplitude1(p.fAmplitude1),
  fAmplitude2(p.fAmplitude2),
  fAmplitude3(p.fAmplitude3),
  fRiseTime(p.fRiseTime),
  fFallTime1(p.fFallTime1),
  fFallTime2(p.fFallTime2),
  fFallTime3(p.fFallTime3),
  fDestPixelId(p.fDestPixelId)
{ }


PEPulse_ROOT&
PEPulse_ROOT::operator>>(mevt::ScintillatorSimData::PEPulse& p)
{
  p.fT0 = fT0;
  p.fAmplitude1 = fAmplitude1;
  p.fAmplitude2 = fAmplitude2;
  p.fAmplitude3 = fAmplitude3;
  p.fRiseTime = fRiseTime;
  p.fFallTime1 = fFallTime1;
  p.fFallTime2 = fFallTime2;
  p.fFallTime3 = fFallTime3;
  p.fDestPixelId = fDestPixelId;
  return *this;
}


PhotonTime_ROOT::PhotonTime_ROOT(const mevt::ScintillatorSimData::PhotonTime& t) :
  fTime(t.fTime)
{ }


PhotonTime_ROOT&
PhotonTime_ROOT::operator>>(mevt::ScintillatorSimData::PhotonTime& t)
{
  t.fTime = fTime;
  return *this;
}


ScintillatorSimData_ROOT::ScintillatorSimData_ROOT(const mevt::ScintillatorSimData& s) :
  fNParticles(s.GetNumberOfParticles()),
  fNumberOfInjectedMuons(s.GetNumberOfInjectedMuons()),
  fNumberOfCornerClippingMuons(s.GetNumberOfCornerClippingMuons()),
  fNumberOfElectrons(s.GetNumberOfElectrons()),
  fEnergyDeposit(s.GetEnergyDeposit()),
  fEnergyDepositMuons(s.GetEnergyDepositMuons()),
  fParticles(s.ParticlesBegin(), s.ParticlesEnd()),
  fAnalogicTrace(s.HasAnalogicTrace() ? new TraceD_ROOT(s.GetAnalogicTrace()) : nullptr),
  fSPEPulses(s.SPEPulsesBegin(), s.SPEPulsesEnd()),
  fPEPulses(s.PEPulsesBegin(), s.PEPulsesEnd()),
  fPhotonTimes(s.PhotonTimesBegin(), s.PhotonTimesEnd())
{ }


ScintillatorSimData_ROOT::~ScintillatorSimData_ROOT()
{
  delete fAnalogicTrace;
  
}


ScintillatorSimData_ROOT&
ScintillatorSimData_ROOT::operator=(const ScintillatorSimData_ROOT& s)
{
  if (this != &s) {
    fNParticles = s.fNParticles;
    fNumberOfInjectedMuons = s.fNumberOfInjectedMuons;
    fNumberOfCornerClippingMuons = s.fNumberOfCornerClippingMuons;
    fNumberOfElectrons = s.fNumberOfElectrons;
    fEnergyDeposit = s.fEnergyDeposit;
    fEnergyDepositMuons = s.fEnergyDepositMuons;
    fParticles = s.fParticles;
    DeepCopy(fAnalogicTrace, s.fAnalogicTrace);
    fSPEPulses = s.fSPEPulses;
    fPEPulses = s.fPEPulses;
    fPhotonTimes = s.fPhotonTimes;
  }
  return *this;
}


void
ScintillatorSimData_ROOT::operator>>(mevt::ScintillatorSimData& s)
  const
{
  // fNParticles will be set by AddParticle
  s.SetNumberOfInjectedMuons(fNumberOfInjectedMuons);
  s.SetNumberOfCornerClippingMuons(fNumberOfCornerClippingMuons);
  s.SetNumberOfElectrons(fNumberOfElectrons);
  s.SetEnergyDeposit(fEnergyDeposit);
  s.SetEnergyDepositMuons(fEnergyDepositMuons);

  for (const auto& p : fParticles)
    s.AddParticle(p.Convert());

  if (fAnalogicTrace) {
    if (!s.HasAnalogicTrace())
      s.MakeAnalogicTrace();
    *fAnalogicTrace >> s.GetAnalogicTrace();
  }


  for (const auto& p : fSPEPulses)
    s.AddSPEPulse(p.fMu, p.fSigma, p.fAmplitude, p.fDestPixelId);

  for (const auto& p : fPEPulses)
    s.AddPEPulse(p.fT0,
                 p.fAmplitude1, p.fAmplitude2, p.fAmplitude3,
                 p.fRiseTime,
                 p.fFallTime1, p.fFallTime2, p.fFallTime3,
                 p.fDestPixelId);

  for (const auto& p : fPhotonTimes)
     s.AddPhotonTime(p.fTime);
}