ShowerSimData_ROOT.cc

Go to the documentation of this file.

#include <io/ShowerSimData_ROOT.h>
#include <io/Point_ROOT.h>
#include <io/CoordinateSystem_ROOT.h>
#include <io/Vector_ROOT.h>
#include <io/TimeStamp_ROOT.h>
#include <io/GaisserHillasParameters_ROOT.h>
#include <io/LaserData_ROOT.h>
#include <io/AtmosphereParameters_ROOT.h>
#include <io/StreamerUtilities.h>
#include <io/GenParticle_ROOT.h>
#include <io/RadioSimulation_ROOT.h>

using namespace io;


ClassImp(ShowerSimData_ROOT)


ShowerSimData_ROOT::~ShowerSimData_ROOT()
{
  delete fParticleCS;
  delete fShowerCS;
  delete fLocalCS;
  delete fTimeStamp;
  delete fGHParameters;
  delete fLaserData;
  delete fAtmosphereParameters;
  delete fGenParticle;
  delete fRadioSimulation;
}


ShowerSimData_ROOT&
ShowerSimData_ROOT::operator=(const ShowerSimData_ROOT& s)
{
  if (this != &s) {
    fPrimaryParticle = s.fPrimaryParticle;
    fEnergy = s.fEnergy;
    fMuonNumber = s.fMuonNumber;
    fMuonWeightScale = s.fMuonWeightScale;
    fCalorimetricEnergy = s.fCalorimetricEnergy;
    fElectromagneticEnergy = s.fElectromagneticEnergy;
    fGroundParticleCoordinateSystemAzimuth = s.fGroundParticleCoordinateSystemAzimuth;
    fGroundParticleCoordinateSystemZenith = s.fGroundParticleCoordinateSystemZenith;
    fEnergyCutoff = s.fEnergyCutoff;
    fEnergyCutoffMuon = s.fEnergyCutoffMuon;
    fMinRadiusCut = s.fMinRadiusCut;
    fMaxRadiusCut = s.fMaxRadiusCut;
    fMinCherenkovWavelength = s.fMinCherenkovWavelength;
    fMaxCherenkovWavelength = s.fMaxCherenkovWavelength;
    fMagneticFieldInclination = s.fMagneticFieldInclination;
    fMagneticFieldDeclination = s.fMagneticFieldDeclination;
    fMagneticFieldStrength = s.fMagneticFieldStrength;
    fXFirst = s.fXFirst;
    fXInject = s.fXInject;
    fDistanceOfShowerMaximum = s.fDistanceOfShowerMaximum;
    fNmu = s.fNmu;
    fXmaxMu = s.fXmaxMu;
    fShowerNumber = s.fShowerNumber;
    fShowerRunId = s.fShowerRunId;
    DeepCopy(fParticleCS, s.fParticleCS);
    DeepCopy(fShowerCS, s.fShowerCS);
    DeepCopy(fLocalCS, s.fLocalCS);    
    DeepCopy(fTimeStamp, s.fTimeStamp);
    fLongitudinalProfile = s.fLongitudinalProfile;
    fFluorescencePhotons = s.fFluorescencePhotons;
    fCherenkovPhotons = s.fCherenkovPhotons;
    fCherenkovBeamPhotons = s.fCherenkovBeamPhotons;
    fCherenkovBeamProductionPhotons = s.fCherenkovBeamProductionPhotons;
    DeepCopy(fGHParameters, s.fGHParameters);
    DeepCopy(fLaserData, s.fLaserData);
    DeepCopy(fAtmosphereParameters, s.fAtmosphereParameters);
    DeepCopy(fGenParticle, s.fGenParticle);
    DeepCopy(fRadioSimulation, s.fRadioSimulation);
  }
  return *this;
}


ShowerSimData_ROOT::ShowerSimData_ROOT(const evt::ShowerSimData& s) :
  fPrimaryParticle(s.GetPrimaryParticle()),
  fEnergy(s.GetEnergy()),
  fMuonNumber(s.GetMuonNumber()),
  fMuonWeightScale(s.GetMuonWeightScale()),
  fCalorimetricEnergy(s.GetCalorimetricEnergy()),
  fElectromagneticEnergy(s.GetElectromagneticEnergy()),
  fGroundParticleCoordinateSystemAzimuth(s.fGroundParticleCoordinateSystemAzimuth),
  fGroundParticleCoordinateSystemZenith(s.fGroundParticleCoordinateSystemZenith),
  fEnergyCutoff(s.GetEnergyCutoff()),
  fEnergyCutoffMuon(s.GetEnergyCutoff(evt::ShowerSimData::eMuon)),
  fMinRadiusCut(s.GetMinRadiusCut()),
  fMaxRadiusCut(s.GetMaxRadiusCut()),
  fMinCherenkovWavelength(s.GetMinCherenkovWavelength()),
  fMaxCherenkovWavelength(s.GetMaxCherenkovWavelength()),
  fMagneticFieldInclination(s.GetMagneticFieldInclination()),
  fMagneticFieldDeclination(s.GetMagneticFieldDeclination()),
  fMagneticFieldStrength(s.GetMagneticFieldStrength()),
  fShowerNumber(s.GetShowerNumber()),
  fShowerRunId(s.GetShowerRunId()),
  fXFirst(s.GetXFirst()),
  fXInject(s.GetXInject()),
  fDistanceOfShowerMaximum(s.GetDistanceOfShowerMaximum()),
  fNmu(s.GetNmu()),
  fXmaxMu(s.GetXmaxMu()),
  fLongitudinalProfile(s.fLongitudinalProfile),
  fFluorescencePhotons(s.fFluorescencePhotons),
  fCherenkovPhotons(s.fCherenkovPhotons),
  fCherenkovBeamPhotons(s.fCherenkovBeamPhotons),
  fCherenkovBeamProductionPhotons(s.fCherenkovBeamProductionPhotons)
{
  if (s.HasGeometry()) {
    ConvertCopy(fParticleCS, s.GetGroundParticleCoordinateSystem());
    ConvertCopy(fShowerCS, s.GetShowerCoordinateSystem());
    ConvertCopy(fLocalCS, s.GetLocalCoordinateSystem());
  }
  if (s.HasTimeStamp())
    ConvertCopy(fTimeStamp, s.GetTimeStamp());
  if (s.HasGHParameters())
    ConvertCopy(fGHParameters, s.GetGHParameters());
  if (s.HasLaserData())
    ConvertCopy(fLaserData, s.GetLaserData());
  if (s.HasAtmosphereParameters())
    ConvertCopy(fAtmosphereParameters, s.GetAtmosphereParameters());
  if (s.HasParticleTree()) {
    ConvertCopy(fGenParticle, s.GetParticleTree());
  }
  if (s.HasRadioSimulation()) {
    ConvertCopy(fRadioSimulation, s.GetRadioSimulation());
  }
}


void
ShowerSimData_ROOT::operator>>(evt::ShowerSimData& s)
  const
{
  s.SetPrimaryParticle(fPrimaryParticle);
  s.SetEnergy(fEnergy);
  s.SetMuonNumber(fMuonNumber);
  s.SetMuonWeightScale(fMuonWeightScale);
  s.SetCalorimetricEnergy(fCalorimetricEnergy);
  s.SetElectromagneticEnergy(fElectromagneticEnergy);
  s.SetGroundParticleCoordinateSystemAzimuth(fGroundParticleCoordinateSystemAzimuth);
  s.SetGroundParticleCoordinateSystemZenith(fGroundParticleCoordinateSystemZenith);
  s.SetEnergyCutoff(fEnergyCutoff);
  s.SetEnergyCutoff(fEnergyCutoffMuon, evt::ShowerSimData::eMuon);
  s.SetMinRadiusCut(fMinRadiusCut);
  s.SetMaxRadiusCut(fMaxRadiusCut);
  s.SetMinCherenkovWavelength(fMinCherenkovWavelength);
  s.SetMaxCherenkovWavelength(fMaxCherenkovWavelength);
  s.SetMagneticFieldInclination(fMagneticFieldInclination);
  s.SetMagneticFieldDeclination(fMagneticFieldDeclination);
  s.SetMagneticFieldStrength(fMagneticFieldStrength);
  s.SetShowerNumber(fShowerNumber);
  s.SetShowerRunId(fShowerRunId);
  s.SetXFirst(fXFirst);
  s.SetXInject(fXInject);
  s.SetDistanceOfShowerMaximum(fDistanceOfShowerMaximum);

  s.SetNmu(fNmu);
  s.SetXmaxMu(fXmaxMu);

  if (fParticleCS)
    *fParticleCS >> s.fParticleCS;

  if (fShowerCS)
    *fShowerCS >> s.fShowerCS;

  if (fLocalCS)
    *fLocalCS >> s.fLocalCS;

  if (fTimeStamp)
    s.MakeTimeStamp(fTimeStamp->Convert());

  fLongitudinalProfile >> s.fLongitudinalProfile;
  fFluorescencePhotons >> s.fFluorescencePhotons;
  fCherenkovPhotons >> s.fCherenkovPhotons;
  fCherenkovBeamPhotons >> s.fCherenkovBeamPhotons;
  fCherenkovBeamProductionPhotons >> s.fCherenkovBeamProductionPhotons;

  if (fGHParameters) {
    evt::VGaisserHillasParameter* ghp = nullptr;
    *fGHParameters >> ghp;
    s.MakeGHParameters(*ghp);
  }

  if (fLaserData) {
    if (!s.HasLaserData())
      s.MakeLaserData();
    *fLaserData >> *s.fLaserData;
  }

  if (fAtmosphereParameters)
    s.MakeAtmosphereParameters(fAtmosphereParameters->Convert());

  if (fGenParticle)
    s.SetParticleTree(fGenParticle->Convert());

  if (fRadioSimulation) {
    if (!s.HasRadioSimulation())
      s.MakeRadioSimulation();
    *fRadioSimulation >> *s.fRadioSimulation;
  }
}