cmpSDEvent.cc

Go to the documentation of this file.

#include "cmpSDEvent.h"

#include <iostream>
#include <iomanip>
#include <sstream>
#include <string>
#include <map>

#include <adst/DetectorGeometry.h>
#include <adst/FileInfo.h>
#include <adst/RecEvent.h>
#include <adst/SDEvent.h>
#include <adst/SdRecShower.h>

#include "testlib.h"
#include "adstReading.h"

using namespace std;

void
cmpSDEvent(const SDEvent& ref, const SDEvent& test)
{
  isDiag(test.GetEventId(), ref.GetEventId(), "SDEvent.GetEventId");
  isDiag(test.GetYYMMDD(), ref.GetYYMMDD(), "SDEvent.GetYYMMDD");
  isDiag(test.GetHHMMSS(), ref.GetHHMMSS(), "SDEvent.GetHHMMSS");
  isDiag(test.GetGPSSecond(), ref.GetGPSSecond(), "SDEvent.GetGPSSecond");
  isDiag(test.GetGPSNanoSecond(), ref.GetGPSNanoSecond(), "SDEvent.GetGPSNanoSecond");
  isDiag(test.HasAxis(), ref.HasAxis(), "SDEvent.HasAxis");
  isDiag(test.HasLDF(), ref.HasLDF(), "SDEvent.HasLDF");
  isDiag(test.HasCurvature(), ref.HasCurvature(), "SDEvent.HasCurvature");
  isDiag(test.IsT5(), ref.IsT5(), "SDEvent.IsT5");
  isDiag(test.Is6T5(), ref.Is6T5(), "SDEvent.Is6T5");
  isDiag(test.Is5T5(), ref.Is5T5(), "SDEvent.Is5T5");
  isDiag(test.HasStations(), ref.HasStations(), "SDEvent.HasStations");
  isDiag(test.IsSaturated(), ref.IsSaturated(), "SDEvent.IsSaturated");
  isDiag(test.GetRecLevel(), ref.GetRecLevel(), "SDEvent.GetRecLevel");
  isDiag(test.GetNoOfAccidentalStations(), ref.GetNoOfAccidentalStations(), "SDEvent.GetNoOfAccidentalStations");
  isDiag(test.GetNumberOfAccidentalStations(), ref.GetNumberOfAccidentalStations(), "SDEvent.GetNumberOfAccidentalStations");
  isDiag(test.GetNoOfAccidentalStations(), ref.GetNoOfAccidentalStations(), "SDEvent.GetNoOfAccidentalStations");
  isDiag(test.GetNoOfRemovedStations(), ref.GetNoOfRemovedStations(), "SDEvent.GetNoOfRemovedStations");
  isDiag(test.GetNumberOfCandidates(), ref.GetNumberOfCandidates(), "SDEvent.GetNumberOfCandidates");
  isDiag(test.GetT4Trigger(), ref.GetT4Trigger(), "SDEvent.GetT4Trigger");
  isDiag(test.GetT5Trigger(), ref.GetT5Trigger(), "SDEvent.GetT5Trigger");
  isDiag(test.GetTriggerSender(), ref.GetTriggerSender(), "SDEvent.GetTriggerSender");
  isDiag(test.GetTriggerAlgorithm(), ref.GetTriggerAlgorithm(), "SDEvent.GetTriggerAlgorithm");
  isDiag(test.IsLightning(), ref.IsLightning(), "SDEvent.IsLightning");
  isDiag(test.ReturnT4Name(), ref.ReturnT4Name(), "SDEvent.ReturnT4Name");
  isDiag(test.ReturnT5Name(), ref.ReturnT5Name(), "SDEvent.ReturnT5Name");
  isDiag(test.HasVEMTraces(), ref.HasVEMTraces(), "SDEvent.HasVEMTraces");
  isDiag(test.HasAllTraces(), ref.HasAllTraces(), "SDEvent.HasAllTraces");
  isDiag(test.HasGenStations(), ref.HasGenStations(), "SDEvent.HasGenStations");
  isDiag(test.HasParticles(), ref.HasParticles(), "SDEvent.HasParticles");
  isDiag(test.HasPEDistributions(), ref.HasPEDistributions(), "SDEvent.HasPEDistributions");
  isDiag(test.GetBadPeriodId(), ref.GetBadPeriodId(), "SDEvent.GetBadPeriodId");

// [a-z0-9A-Z]+
//  isDiag(test.GetSdRecShower().\&(), ref.GetSdRecShower().\&(), "sdRecShower.\&");

// GetRadiusPoint
// GetLogMPDHistogram
  isDiag(test.GetSdRecShower().GetAngleChi2(), ref.GetSdRecShower().GetAngleChi2(), "sdRecShower.GetAngleChi2");
  isDiag(test.GetSdRecShower().GetAngleNDoF(), ref.GetSdRecShower().GetAngleNDoF(), "sdRecShower.GetAngleNDoF");

  if (test.HasCurvature() && ref.HasCurvature()) {
     isRelEpsDiag(test.GetSdRecShower().GetCurvature(), ref.GetSdRecShower().GetCurvature(), 1.e-7, "sdRecShower.GetCurvature");
     isRelEpsDiag(test.GetSdRecShower().GetCurvatureError(), ref.GetSdRecShower().GetCurvatureError(), 1.e-7, "sdRecShower.GetCurvatureError");
     isRelEpsDiag(test.GetSdRecShower().GetRadiusOfCurvature(), ref.GetSdRecShower().GetRadiusOfCurvature(), 1.e-7, "sdRecShower.GetRadiusOfCurvature");
     isRelEpsDiag(test.GetSdRecShower().GetRadiusOfCurvatureError(), ref.GetSdRecShower().GetRadiusOfCurvatureError(), 1.e-7, "sdRecShower.GetRadiusOfCurvatureError");
  }

  if (test.HasLDF() && ref.HasLDF()) {
     isRelEpsDiag(test.GetSdRecShower().GetROpt(), ref.GetSdRecShower().GetROpt(), 1e-7, "sdRecShower.GetROpt");
     isRelEpsDiag(test.GetSdRecShower().GetROptError(), ref.GetSdRecShower().GetROptError(), 1e-7, "sdRecShower.GetROptError");
     isRelEpsDiag(test.GetSdRecShower().GetTimeResidualMean(), ref.GetSdRecShower().GetTimeResidualMean(), 1e-7, "sdRecShower.GetTimeResidualMean");
     isRelEpsDiag(test.GetSdRecShower().GetTimeResidualSpread(), ref.GetSdRecShower().GetTimeResidualSpread(), 1e-7, "sdRecShower.GetTimeResidualSpread");
     isRelEpsDiag(test.GetSdRecShower().GetS1000(), ref.GetSdRecShower().GetS1000(), 1e-7, "sdRecShower.GetS1000");
     isRelEpsDiag(test.GetSdRecShower().GetS1000Error(), ref.GetSdRecShower().GetS1000Error(), 1e-7, "sdRecShower.GetS1000Error");
     isRelEpsDiag(test.GetSdRecShower().GetS1000TotalError(), ref.GetSdRecShower().GetS1000TotalError(), 1e-7, "sdRecShower.GetS1000TotalError");
     isRelEpsDiag(test.GetSdRecShower().GetS1000BetaSys(), ref.GetSdRecShower().GetS1000BetaSys(), 1e-7, "sdRecShower.GetS1000BetaSys");
     isRelEpsDiag(test.GetSdRecShower().GetBeta(), ref.GetSdRecShower().GetBeta(), 1e-7, "sdRecShower.GetBeta");
     isRelEpsDiag(test.GetSdRecShower().GetBetaError(), ref.GetSdRecShower().GetBetaError(), 1e-7, "sdRecShower.GetBetaError");
     isRelEpsDiag(test.GetSdRecShower().GetGamma(), ref.GetSdRecShower().GetGamma(), 1e-7, "sdRecShower.GetGamma");
     isRelEpsDiag(test.GetSdRecShower().GetGammaError(), ref.GetSdRecShower().GetGammaError(), 1e-7, "sdRecShower.GetGammaError");
     isRelEpsDiag(test.GetSdRecShower().GetLDFChi2(), ref.GetSdRecShower().GetLDFChi2(), 1e-7, "sdRecShower.GetLDFChi2");
     isDiag(test.GetSdRecShower().GetLDFNdof(), ref.GetSdRecShower().GetLDFNdof(), "sdRecShower.GetLDFNdof");
     isRelEpsDiag(test.GetSdRecShower().GetLDFLikelihood(), ref.GetSdRecShower().GetLDFLikelihood(), 1e-7, "sdRecShower.GetLDFLikelihood");
  }

  isDiag(test.GetSdRecShower().GetLDFStatus(), ref.GetSdRecShower().GetLDFStatus(), "sdRecShower.GetLDFStatus");

  isRelEpsDiag(test.GetSdRecShower().SignalUncertainty(), ref.GetSdRecShower().SignalUncertainty(), 1.e-7, "sdRecShower.SignalUncertainty");
  isDiag(test.GetSdRecShower().HasLogMPDHistogram(), ref.GetSdRecShower().HasLogMPDHistogram(), "sdRecShower.HasLogMPDHistogram");

// [a-z0-9A-Z]+
//  isDiag(test.GetSdRecShower().GetLDF().\&(), ref.GetSdRecShower().GetLDF().\&(), "ldf.\&");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetS1000(), ref.GetSdRecShower().GetLDF().GetS1000(), 1e-7, "ldf.GetS1000");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetS1000Error(), ref.GetSdRecShower().GetLDF().GetS1000Error(), 1e-7, "ldf.GetS1000Error");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetS1000BetaSys(), ref.GetSdRecShower().GetLDF().GetS1000BetaSys(), 1e-7, "ldf.GetS1000BetaSys");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetS1000TotalError(), ref.GetSdRecShower().GetLDF().GetS1000TotalError(), 1e-7, "ldf.GetS1000TotalError");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetBeta(), ref.GetSdRecShower().GetLDF().GetBeta(), 1e-7, "ldf.GetBeta");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetBetaError(), ref.GetSdRecShower().GetLDF().GetBetaError(), 1e-7, "ldf.GetBetaError");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetBetaSys(), ref.GetSdRecShower().GetLDF().GetBetaSys(), 1e-7, "ldf.GetBetaSys");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetGamma(), ref.GetSdRecShower().GetLDF().GetGamma(), 1e-7, "ldf.GetGamma");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetGammaError(), ref.GetSdRecShower().GetLDF().GetGammaError(), 1e-7, "ldf.GetGammaError");
  {
    ToDo t("It seems this is not set in the reference (this is ne in trunk)");
    isRelEpsDiag(test.GetSdRecShower().GetLDF().GetNKGFermiMu(), ref.GetSdRecShower().GetLDF().GetNKGFermiMu(), 1e-7, "ldf.GetNKGFermiMu");
    isRelEpsDiag(test.GetSdRecShower().GetLDF().GetNKGFermiTau(), ref.GetSdRecShower().GetLDF().GetNKGFermiTau(), 1e-7, "ldf.GetNKGFermiTau");
  }
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetLDFChi2(), ref.GetSdRecShower().GetLDF().GetLDFChi2(), 1e-7, "ldf.GetLDFChi2");
  isDiag(test.GetSdRecShower().GetLDF().GetLDFNdof(), ref.GetSdRecShower().GetLDF().GetLDFNdof(), "ldf.GetLDFNdof");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetLDFLikelihood(), ref.GetSdRecShower().GetLDF().GetLDFLikelihood(), 1e-7, "ldf.GetLDFLikelihood");
  isRelEpsDiag(test.GetSdRecShower().GetLDF().GetLDFStatus(), ref.GetSdRecShower().GetLDF().GetLDFStatus(), 1e-7, "ldf.GetLDFStatus");

// GetRiseTimeResults
// [a-z0-9A-Z]+
//  isDiag(test.GetSdRecShower().GetRiseTimeResults().\&(), ref.GetSdRecShower().GetRiseTimeResults().\&(), "risetime.\&");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetRiseTime1000(), ref.GetSdRecShower().GetRiseTimeResults().GetRiseTime1000(), 1e-7, "risetime.GetRiseTime1000");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetRiseTime1000Error(), ref.GetSdRecShower().GetRiseTimeResults().GetRiseTime1000Error(), 1e-7, "risetime.GetRiseTime1000Error");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetRiseTimeNDF(), ref.GetSdRecShower().GetRiseTimeResults().GetRiseTimeNDF(), 1e-7, "risetime.GetRiseTimeNDF");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetRiseTimeChi2(), ref.GetSdRecShower().GetRiseTimeResults().GetRiseTimeChi2(), 1e-7, "risetime.GetRiseTimeChi2");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetAlpha(), ref.GetSdRecShower().GetRiseTimeResults().GetAlpha(), 1e-7, "risetime.GetAlpha");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetBeta(), ref.GetSdRecShower().GetRiseTimeResults().GetBeta(), 1e-7, "risetime.GetBeta");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetXmax(), ref.GetSdRecShower().GetRiseTimeResults().GetXmax(), 1e-7, "risetime.GetXmax");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetXmaxErrorUp(), ref.GetSdRecShower().GetRiseTimeResults().GetXmaxErrorUp(), 1e-7, "risetime.GetXmaxErrorUp");
  isRelEpsDiag(test.GetSdRecShower().GetRiseTimeResults().GetXmaxErrorDown(), ref.GetSdRecShower().GetRiseTimeResults().GetXmaxErrorDown(), 1e-7, "risetime.GetXmaxErrorDown");




}