DataQualityInterface.cxx

Go to the documentation of this file.

#include "DataQualityInterface.h"
  
#include "MessageService/MsgService.h"

#include "BeamDataUtil/BDSpillAccessor.h"
#include "BeamDataUtil/BeamMonSpill.h"
#include "BeamDataUtil/BMSpillAna.h"
#include "SpillTiming/SpillTimeFinder.h"
#include "SpillTiming/SpillServerMonFinder.h"
 
#include "DcsUser/HvStatus.h"
#include "DcsUser/HvStatusFinder.h"
#include "DcsUser/CoilTools.h"

#include <cmath>

//
// $Log $
//
 
ClassImp(DataQualityInterface)
 
CVSID("$Id: DataQualityInterface.cxx,v 1.9 2008/08/25 13:09:18 blake Exp $");
 
DataQualityInterface::DataQualityInterface() :
  fBeamTime(-1),
  fBeamStatus(-1),
  fNearDetStatus(-1),
  fBeamType(-1),
  fProtonsOnTarget(0.0),
  fDetectorTime(-1),
  fDetectorStatus(-1),
  fCoilStatus(-1),
  fHvStatus(-1),
  fGpsStatus(-1),
  fGpsError(-999),
  fDataQualityLookup(0)
{

}
  
DataQualityInterface::DataQualityInterface(const DataQualityInterface& rhs) :
  TObject(rhs)
{

}
  
DataQualityInterface::~DataQualityInterface()
{
  if(fDataQualityLookup) delete fDataQualityLookup;
}
 
Int_t DataQualityInterface::GetSpillStatus(VldContext vldc)
{
  this->ProcessBeamStatus(vldc,1);
  return fNearDetStatus;
}

Int_t DataQualityInterface::GetBeamStatus(VldContext vldc)
{
  this->ProcessBeamStatus(vldc,1);
  return fBeamStatus;
}
   
Int_t DataQualityInterface::GetBeamType(VldContext vldc)
{
  this->ProcessBeamStatus(vldc,1);
  return fBeamType;
}

Double_t DataQualityInterface::GetProtonsOnTarget(VldContext vldc)
{
  this->ProcessBeamStatus(vldc,1);
  return fProtonsOnTarget;
}

Int_t DataQualityInterface::GetDetectorStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,1);
  return fDetectorStatus;
}

Int_t DataQualityInterface::GetHvStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,1);
  return fHvStatus;
}

Int_t DataQualityInterface::GetCoilStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,1);
  return fCoilStatus;
}

Int_t DataQualityInterface::GetGpsStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,1);
  return fGpsStatus;
}

Int_t DataQualityInterface::LookupDetectorStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,0);
  return fDetectorStatus;
}

Int_t DataQualityInterface::LookupHvStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,0);
  return fHvStatus;
}

Int_t DataQualityInterface::LookupCoilStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,0);
  return fCoilStatus;
}

Int_t DataQualityInterface::LookupGpsStatus(VldContext vldc)
{
  this->ProcessDetectorStatus(vldc,0);
  return fGpsStatus;
}

void DataQualityInterface::ProcessBeamStatus(VldContext vldc, Bool_t useDB)
{
  if( vldc.GetTimeStamp().GetSec()!=fBeamTime ){
    MSG("DataQualityFilter",Msg::kVerbose) << " DataQualityInterface::ProcessBeamStatus(VldContext) " << endl;

    fBeamTime=vldc.GetTimeStamp().GetSec();
    fBeamStatus=-1;
    fNearDetStatus=-1;
    fBeamMonStatus=-1;
    fNearDetSpill=0;
    fBeamMonSpill=0;
    fBeamType=-1;
    fProtonsOnTarget=0.0;

    // USE BEAM DATABASE TABLES
    // ========================
    if( useDB ){

      // Check NearDet Spill Times
      SpillTimeFinder& NDaccess = SpillTimeFinder::Instance();
      fNearDetStatus = NDaccess.DataIsAvailable(vldc);
      const SpillTimeND& ndspill = NDaccess.GetNearestSpill(vldc);
      VldTimeStamp NDtime = ndspill.GetTimeStamp();

      Double_t nd_subtime = NDaccess.GetTimeToNearestSpill(vldc);

      if( nd_subtime>=0.0 && nd_subtime<1.0 ){
        fNearDetSpill = 1;
      }

      MSG("DataQualityFilter",Msg::kVerbose) << " NearDet [status=" << fNearDetStatus << ", spill=" << fNearDetSpill << ", time=" << nd_subtime << "]" << endl;
      
      // Check Beam Monitoring
      BDSpillAccessor& BMaccess = BDSpillAccessor::Get();
      const BeamMonSpill& BMspill = *(BMaccess.LoadSpill(NDtime));
      VldTimeStamp BMtime = BMspill.SpillTime();

      Double_t bm_subtime = ( BMtime.GetSec()-NDtime.GetSec() )
                   + 1.0e-9*( BMtime.GetNanoSec()-NDtime.GetNanoSec() );

      BMSpillAna BMana;
      BMana.UseDatabaseCuts();
      BMana.SetSpill(BMspill); 
      BMana.SetTimeDiff(bm_subtime);

      Double_t pots = 0.0;
      if(BMspill.fTortgt !=0.)      pots = BMspill.fTortgt;
      else if(BMspill.fTrtgtd !=0.) pots = BMspill.fTrtgtd;
      else if(BMspill.fTor101 !=0.) pots = BMspill.fTor101;

      if( fNearDetSpill
       && abs(bm_subtime)<1.0 ){
        fBeamMonSpill = 1;
        fProtonsOnTarget = pots;
        fBeamType = BMspill.BeamType();
        fBeamMonStatus = BMana.SelectSpill();
      }

      MSG("DataQualityFilter",Msg::kVerbose) << " BeamMon [spill=" << fBeamMonSpill << ", isgood=" << fBeamMonStatus << ", time=" << nd_subtime+bm_subtime << ", type=" << fBeamType << ", PoTs=" << fProtonsOnTarget << "]" << endl;

    }

    // USE BEAM LOOKUP TABLES
    // ======================
    else{

    }

    // APPLY BEAM QUALITY CUTS
    // =======================
    // check the beam monitoring
    //    1 = beam monitoring is good
    //    0 = beam monitoring is bad
    //   -1 = no beam monitoring information 

    if( fBeamMonStatus==1 && fNearDetStatus==1 ){
      fBeamStatus=1;
    }
    if( fBeamMonStatus==0 || fNearDetStatus==0 ){
      fBeamStatus=0;
    }

  }

  return;
}
  
void DataQualityInterface::ProcessDetectorStatus(VldContext vldc, Bool_t useDB)
{
  if( vldc.GetTimeStamp().GetSec()!=fDetectorTime ){
    MSG("DataQualityFilter",Msg::kVerbose) << " DataQualityInterface::ProcessDetectorStatus(VldContext) " << endl;

    fDetectorTime=vldc.GetTimeStamp().GetSec();
    fDetectorStatus=-1;    
    fCoilStatus=-1;
    fHvStatus=-1;
    fGpsStatus=1;
    fGpsError=-999;

    // USE DETECTOR DATABASE TABLES
    // ============================
    if( useDB ){

      // Check Coil Current
      if( CoilTools::IsOK(vldc) ) fCoilStatus=1; else fCoilStatus=0; 
      MSG("DataQualityFilter",Msg::kVerbose) << " CoilCurrentStatus [status=" << fCoilStatus << "]" << endl;
      MSG("DataQualityFilter",Msg::kVerbose) << " CoilCurrent [lookup=0, status=" << fCoilStatus << "]" << endl;

      // Check High Voltage
      HvStatus::HvStatus_t status = HvStatusFinder::Instance().GetHvStatus(vldc);
      if( (status&HvStatus::kOK)==(HvStatus::kOK) ) fHvStatus=1;
      if( (status&HvStatus::kBad)==(HvStatus::kBad) ) fHvStatus=0;
      MSG("DataQualityFilter",Msg::kVerbose) << " HighVoltageStatus [statusword=" << status << ", status=" << fHvStatus << "]" << endl;
      MSG("DataQualityFilter",Msg::kVerbose) << " HighVoltage [lookup=0, status=" << fHvStatus << "]" << endl;      

      // Check GPS Error
      SpillServerMonFinder& GPSaccess = SpillServerMonFinder::Instance();
      if( GPSaccess.DataIsAvailable(vldc) ){
        const SpillServerMon& spill_nearest = GPSaccess.GetNearestSpill(vldc);

        Double_t gps_subtime = ( spill_nearest.GetTimeStamp().GetSec()-vldc.GetTimeStamp().GetSec() )
                      + 1.0e-9*( spill_nearest.GetTimeStamp().GetNanoSec()-vldc.GetTimeStamp().GetNanoSec() );

        if( gps_subtime>=0.0 && gps_subtime<1.0 ){
          fGpsError = 0;
          if( spill_nearest.GetSpillType()==1 ){
            fGpsError = spill_nearest.GetSpillTimeError();
            if( fGpsError>1000 ) fGpsStatus = 0;
          }
        }
      }
      MSG("DataQualityFilter",Msg::kVerbose) << " GpsErrorStatus [error=" << fGpsError << ", status=" << fGpsStatus << "]" << endl;
      MSG("DataQualityFilter",Msg::kVerbose) << " GpsError [lookup=0, status=" << fGpsStatus << "]" << endl;

    }

    // USE DETECTOR LOOKUP TABLES
    // ==========================
    else{
      if(!fDataQualityLookup) fDataQualityLookup = new DataQualityLookup();
  
      // Check Coil Current
      fCoilStatus=fDataQualityLookup->GetCoilStatus(vldc);
      MSG("DataQualityFilter",Msg::kVerbose) << " CoilCurrent [lookup=1, status=" << fCoilStatus << "]" << endl;

      // Check High Voltage
      fHvStatus=fDataQualityLookup->GetHvStatus(vldc);
      MSG("DataQualityFilter",Msg::kVerbose) << " HighVoltage [lookup=1, status=" << fHvStatus << "]" << endl;

      // Check GPS Error;
      fGpsStatus=1;
      MSG("DataQualityFilter",Msg::kVerbose) << " GpsError [lookup=1, status=" << fGpsStatus << "]" << endl;

    }

    // APPLY DETECTOR QUALITY CUTS
    // ===========================
    //  1 = coil and HV are good
    //  0 = coil or HV are bad or unknown
    // -1 = no detector information

    if( fCoilStatus==1 && fHvStatus==1 && fGpsStatus==1 ){
      fDetectorStatus=1;
    }
    if( fCoilStatus==0 || fHvStatus==0 || fGpsStatus==0 ){
      fDetectorStatus=0;
    }

  }

  return;
}

---