SdStationChecker.cc

Go to the documentation of this file.

#include "SdStationChecker.h"
#include <evt/Event.h>
#include <sevt/SEvent.h>
#include <sevt/EventTrigger.h>
#include <utl/String.h>

using namespace fwk;
using namespace sevt;
using namespace utl;
using namespace std;


namespace SdStationCheckerOG {

  void
  MakeInfo(const vector<int>& ids, const string& message)
  {
    if (ids.empty())
      return;
    ostringstream info;
    info << "station" << String::Plural(ids) << ' ' << String::OfSortedIds(ids) << ' ' << message;
    INFO(info);
  }


  VModule::ResultFlag
  SdStationChecker::Run(evt::Event& event)
  {
    if (!event.HasSEvent())
      return eSuccess;
    auto& sEvent = event.GetSEvent();

    vector<int> randomTrigger;
    vector<int> badCompression;
    int noTrigger = 0;

    const bool isCommsCrisis = kCommsCrisis.IsInRange(det::Detector::GetInstance().GetTime());

    // loop on stations
    for (auto& station : sEvent.StationsRange()) {

      if (!station.HasTriggerData()) {
        station.SetRejected(StationConstants::eNoTrigger);
        ++noTrigger;
        continue;
      }

      // T2Life() has to be checked before SetSilent() is applied!
      if (!station.IsT2Life())
        station.SetRejected(StationConstants::eNotAliveT2);
      if (isCommsCrisis && !station.IsT2Life120())
        station.SetRejected(StationConstants::eNotAliveT120);

      const auto& trig = station.GetTriggerData();
      if (trig.GetErrorCode() & 0xff) { // for UUBs errorcodes are above 256.
        if (trig.IsSilent() && !station.IsRejected())
          station.SetSilent();
        else
          station.SetRejected(StationConstants::eErrorCode);
        continue;
      }

      if (trig.IsRandom()) {
        randomTrigger.push_back(station.GetId());
        station.SetRejected(StationConstants::eRandom);
        continue;
      }

      if (!station.HasCalibData()) {
        station.SetRejected(StationConstants::eNoCalibData);
        continue;
      }

      // exclude FADCs with Patrick's data
      if (station.GetCalibData().GetVersion() > 32000) {
        ostringstream warn;
        warn << "Station " << station.GetId() << " has LS calibration version "
             << station.GetCalibData().GetVersion() << '!';
        WARNING(warn);
        station.SetRejected(StationConstants::eNoCalibData);
        continue;
      }

      if (!station.HasGPSData()) {
        station.SetRejected(StationConstants::eNoGPSData);
        continue;
      }

      ApplyTimeCorrection(station);

      // check for "bad compressed data"
      if (sEvent.HasTrigger()) {
        const int trigSecond = sEvent.GetTrigger().GetTime().GetGPSSecond();
        const int timeDiff = int(station.GetGPSData().GetSecond()) - trigSecond;
        if (abs(timeDiff) > 1) {
          station.GetTriggerData().SetErrorCode(StationTriggerData::eBadCompress);
          station.SetRejected(StationConstants::eBadCompress);
          const int sId = station.GetId();
          badCompression.push_back(sId);
          ostringstream info;
          info << "Bad compress data: station " << sId << " has " << timeDiff
               << " sec of time difference to event trigger.";
          INFO(info);
          continue;
        }
      }

    } // end of loop over stations

    MakeInfo(randomTrigger, "with random trigger rejected.");
    MakeInfo(badCompression, "with bad compression data rejected.");

    if (noTrigger) {
      ostringstream info;
      info << noTrigger << " station" << String::Plural(noTrigger)
           << String::OfSortedIds(badCompression) << " without trigger data rejected.";
      INFO(info);
    }

    return eSuccess;
  }


  void
  SdStationChecker::ApplyTimeCorrection(Station& station)
  {
    auto& gpsData = station.GetGPSData();

    // NEW : TAP 26/04/2003 -> From CDAS v1r2 : taking into account Offsets...
    // Warning, part of the field is used for the tick offset:
    // GPS Offset = 0.01*(short)(gps->Offset & 0xffff)
    // Tick Offset = (short)(gps->Offset>>16)
    // New: taking into account 100ns jumps
    // From Moulin Rouge and Dia Noche we found that the TickFall-Tick
    // can be 0, 9, 10, 11 or a big number. The big number could be
    // understood if it is the trigger of another event. It was found
    // that if the dt is 0, there is a 100ns jump in the event, and not
    // in any other case, including big values. Hence this empiric
    // correction
    //
    // This is the code from IoSd v1r2 :
    // gps->NanoSecond =
    // (unsigned int)((gps->Tick*(1000000000.0 + gps->NextST - gps->CurrentST)
    // /gps->Next100) + gps->CurrentST + 0.01*(short)(gps->Offset & 0xffff))
    // -100*(gps->TickFall == gps->Tick);

    const unsigned int tick = gpsData.GetTick();
    const int currentST = gpsData.GetCurrentST();
    const int next100 = gpsData.GetNext100();
    const int nextST = gpsData.GetNextST();

#ifndef IOSD_V1R0
    const unsigned int tickFall = gpsData.GetTickFall();
    const int offset = gpsData.GetOffset();

    const unsigned int nanosecond =
      (unsigned int)((tick * (1e9 + nextST - currentST) / next100) + currentST +
                     0.01 * short(offset & 0xffff)) - 100 * (tickFall == tick);
#else
    const unsigned int nanosecond =
      (unsigned int)((tick * (1e9 + nextST - currentST) / next100) + currentST);
#endif

    gpsData.SetCorrectedNanosecond(nanosecond);
  }

}