DataQualDB.cxx

Go to the documentation of this file.

#include "DataQualDB.h"

#include "RunQuality/RunQualityFinder.h"
#include "DcsUser/HvStatusFinder.h"
#include "DcsUser/CoilTools.h"
#include "SpillTiming/SpillServerMonFinder.h"

// #include "CandNtupleSR/NtpSRDataQuality.h"

Bool_t DataUtil::IsGoodData(const NtpStRecord *st) {

  // far detector data quality
  // =========================
  // number of crates enabled
  // (superseded by DbuFarRunQuality table)
  
  // Int_t cratemask = 0;
  // const NtpSRDataQuality *ntpDataQual;
  // ntpDataQual = &(st->dataquality);
  // cratemask = ntpDataQual->cratemask;
  // if( cratemask!=16 ) return 0;
 
  // get validity context
  // ====================
  const RecCandHeader* Header = &(st->GetHeader());
  VldContext cx = Header->GetVldContext();

  return IsGoodData(cx);
}

Bool_t DataUtil::IsGoodData(const VldContext& cx)
{   
  // If it's not real data, return kOK
  // =================================
  if( cx.GetSimFlag() != SimFlag::kData ){
    return 1;
  }

  // Far Detector Data Selection
  // ===========================
  if( cx.GetDetector() == Detector::kFar ) { 
  
    // requirements for good far detector data: 
    //  (i)   good run - IsGoodDataRUN(cx)
    //  (ii)  good hv - IsGoodDataHV(cx)
    //  (iii) good coil - IsGoodDataCOIL(cx)
    //  (iv)  good gps error - IsGoodDataGPS(cx)

    if( IsGoodDataRUN(cx) && IsGoodDataHV(cx)
     && IsGoodDataCOIL(cx) && IsGoodDataGPS(cx) ) return 1;
    else return 0;
  }
  
  // Near Detector Data Selection
  // ============================
  if( cx.GetDetector() == Detector::kNear ) { 
  
    // requirements for good near detector data: 
    //  (i)  good run - IsGoodDataRUN(cx)
    //  (ii) good coil - IsGoodDataCOIL(cx)

    if( IsGoodDataRUN(cx) 
     && IsGoodDataCOIL(cx) ) return 1;
    else return 0;
  }

  return 1;
}
  
Bool_t DataUtil::IsGoodDataRUN(const VldContext& cx)
{
  // If it's not real data, return kOK
  // =================================
  if( cx.GetSimFlag() != SimFlag::kData ){
    return 1;
  }

  // Far Detector Data - Check Run Status
  // ====================================
  if( cx.GetDetector() == Detector::kFar ) {
    return RunQualityFinder::Instance().IsOK(cx);
  }

  // Near Detector Data - Check Run Status
  // =====================================
  if( cx.GetDetector() == Detector::kNear ) {
    return RunQualityFinder::Instance().IsOK(cx);
  }

  return 1;
}
  
Bool_t DataUtil::IsGoodDataHV(const VldContext& cx)
{
  // If it's not real data, return kOK
  // =================================
  if( cx.GetSimFlag() != SimFlag::kData ){
    return 1;
  }

  // Far Detector Data - Check HV Status
  // ===================================
  // check HV status within 60 secs of snarl
  // (N.B: arguments given to HvStatusFinder: 
  //   validity context, time window, task number)
  if( cx.GetDetector() == Detector::kFar ) {
    HvStatus::HvStatus_t hv_ok = 
      HvStatusFinder::Instance().GetHvStatus(cx,60,1);
    return HvStatus::Good(hv_ok);
  }

  // Near Detector Data, return kOK
  // ==============================
  if( cx.GetDetector() == Detector::kNear ) {
    return 1;
  }

  return 1;
}
  
Bool_t DataUtil::IsGoodDataCOIL(const VldContext& cx)
{
  // If it's not real data, return kOK
  // =================================
  if( cx.GetSimFlag() != SimFlag::kData ){
    return 1;
  }

  // Far Detector Data - Check Coil Status
  // =====================================
  if( cx.GetDetector() == Detector::kFar ) {
    return CoilTools::IsOK(cx);     
  }
  
  // Near Detector Data - Check Coil Status
  // ======================================
  if( cx.GetDetector() == Detector::kNear ) {
    return CoilTools::IsOK(cx);     
  }

  return 1;
}
  
Bool_t DataUtil::IsGoodDataGPS(const VldContext& cx)
{
  // If it's not real data, return kOK
  // =================================
  if( cx.GetSimFlag() != SimFlag::kData ){
    return 1;
  }

  // Far Detector Data - Check GPS Status
  // ====================================
  // get spillserver monitoring block for this snarl,
  // get GPS worst case uncertainty (nsec)
  // spill must be real and within 5 minutes
  if( cx.GetDetector() == Detector::kFar ) {
  
    SpillServerMonFinder& smon = SpillServerMonFinder::Instance();
    const SpillServerMon& spill_near = smon.GetNearestSpill(cx);
    VldTimeStamp dt = spill_near.GetSpillTime()-cx.GetTimeStamp();

    Int_t dt_sec = abs(dt.GetSec());
    Int_t gps_error = spill_near.GetSpillTimeError();
    Int_t spill_type = spill_near.GetSpillType();

    if( spill_type==1 && dt_sec<360 && gps_error>1000 ) return 0;
    else return 1;
  }

  return 1;
}  

Bool_t DataUtil::IsGoodFDData(const NtpStRecord* st)
{ 
  return IsGoodData(st); 
}
  
Bool_t DataUtil::IsGoodFDData(const VldContext& cx)
{ 
  return IsGoodData(cx); 
}

Bool_t DataUtil::IsGoodFDDataRUN(const VldContext& cx)
{ 
  return IsGoodDataRUN(cx); 
}
  
Bool_t DataUtil::IsGoodFDDataHV(const VldContext& cx)
{ 
  return IsGoodDataHV(cx); 
}
  
Bool_t DataUtil::IsGoodFDDataCOIL(const VldContext& cx)
{ 
  return IsGoodDataCOIL(cx); 
}
  
Bool_t DataUtil::IsGoodFDDataGPS(const VldContext& cx)
{ 
  return IsGoodDataGPS(cx); 
}

Bool_t DataUtil::IsGoodNDData(const NtpStRecord* st)
{ 
  return IsGoodData(st); 
}
  
Bool_t DataUtil::IsGoodNDData(const VldContext& cx)
{ 
  return IsGoodData(cx); 
}

Bool_t DataUtil::IsGoodNDDataRUN(const VldContext& cx)
{ 
  return IsGoodDataRUN(cx); 
}
  
Bool_t DataUtil::IsGoodNDDataCOIL(const VldContext& cx)
{ 
  return IsGoodDataCOIL(cx); 
}

---