In This Package:

DaqRunInfoTestAlg.cc

Go to the documentation of this file.

#include "DaqRunInfoTestAlg.h"
#include "DaqRunInfoSvc/DaqRunInfo.h"
#include "DaqRunInfoSvc/DaqCalibRunInfo.h"
#include "DatabaseSvc/IDatabaseSvc.h"
#include "DaqRunInfoSvc/IDaqRunInfoSvc.h"

#include "Conventions/Site.h"
#include "Conventions/SimFlag.h"
#include "Context/TimeStamp.h"
#include "Context/ContextRange.h"
#include "DatabaseInterface/DbiResultPtr.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include "DatabaseInterface/DbiSqlContext.h"

#include <fstream>

DaqRunInfoTestAlg::DaqRunInfoTestAlg(const string& name, ISvcLocator* pSvcLocator)
:GaudiAlgorithm(name, pSvcLocator), m_daqRunInfoSvc(0)
{
  declareProperty("daqRunInfoSvcName", m_daqRunInfoSvcName = "DaqRunInfoSvc",
                  "Name of the service for access run information");
  declareProperty("fileName", m_fileName = "", "Name of the raw data file");
  declareProperty("outputFile", m_outputFile = "runConfig.txt", "Name of the output file holding run configuration");
  declareProperty("runNo", m_runNo = 400, "run number");
  declareProperty("verbose", m_verbose = "simple", "Verbose to print run information, simple or full");
}

DaqRunInfoTestAlg::~DaqRunInfoTestAlg()
{
}

StatusCode DaqRunInfoTestAlg::initialize()
{
  // Initialize the necessary services
  StatusCode sc = this->service(m_daqRunInfoSvcName, m_daqRunInfoSvc, true);
  if(sc.isFailure()){
      error() << "Failed to get service: " << m_daqRunInfoSvcName << endreq;
      return sc;
  }
  return StatusCode::SUCCESS;
}

StatusCode DaqRunInfoTestAlg::execute()
{
  return StatusCode::SUCCESS;
}

StatusCode DaqRunInfoTestAlg::finalize()
{
  if ( m_fileName == ""){
    warning() << "No raw data file name is provided!"<<endreq;
    warning() << "Use run number to access the run information."<<endreq;
  }
  else{ 
    m_runNo = m_daqRunInfoSvc->runNumber(m_fileName);
    info() << "The run number of "<<m_fileName<<" is "<<m_runNo<<endreq; 
  }
  
  if(m_runNo == 0)
    warning() << "No run number is provided!"<<endreq;
  else 
    printRunInfo(m_runNo);

  m_daqRunInfoSvc->release();
  return StatusCode::SUCCESS;
}

StatusCode DaqRunInfoTestAlg::printRunInfo(int runNo)
{
  const DaqRunInfo* row;
  row = m_daqRunInfoSvc->GetRunInfo(runNo);
  if(!row){
    warning() << "runNo: "<<runNo<<" is not found in table DaqRunInfo."<<endreq;
    return StatusCode::FAILURE;
  }
 
  DbiResult* res = row->GetOwner();
  const DbiValidityRec vr = res->GetValidityRec();
  const ContextRange cr = vr.GetContextRange();
  info() <<"runNo: "<<row->GetRunNo()<<endreq;
  info() <<"timeStart: "<<cr.GetTimeStart()<<endreq;
  info() <<"timeEnd: "<<cr.GetTimeEnd()<<endreq;
  info() <<"triggerType: "<<row->GetTriggerType()<<endreq;
  info() <<"runType: "<<row->GetRunType()<<endreq;
  info() <<"detectorMask: "<<row->GetDetectorMask()<<endreq;
  info() <<"partitionName: "<<row->GetPartitionName()<<endreq;
  info() <<"schemaVersion: "<<row->GetSchemaVersion()<<endreq;
  info() <<"dataVersion: "<<row->GetDataVersion()<<endreq;
  info() <<"baseVersion: "<<row->GetBaseVersion()<<endreq;

  string runType = row->GetRunType();  
  if(runType == "ADCalib"){
    info() <<"runNo: "<<runNo<<" is a AD calibration run."<<endreq;
  }

  vector<DayaBay::DaqRunConfig> vec;
  if(m_verbose == "simple")
    m_daqRunInfoSvc->simpRunConfigList(runNo, vec);
  else if (m_verbose == "full")
    m_daqRunInfoSvc->fullRunConfigList(runNo, vec);
  else
    m_daqRunInfoSvc->simpRunConfigList(runNo, vec);
    
  info() <<"runNo: "<<runNo<<" has "<<vec.size()<<" records of DaqRunConfig "<<endreq;
  ofstream output;
  output.open(m_outputFile.c_str());
  output <<"------- Overview of run information ----------"<<endl;
  output <<"runNo: "<<row->GetRunNo()<<endl;
  output <<"timeStart: "<<cr.GetTimeStart()<<endl;
  output <<"timeEnd: "<<cr.GetTimeEnd()<<endl;
  output <<"triggerType: "<<row->GetTriggerType()<<endl;
  output <<"runType: "<<row->GetRunType()<<endl;
  output <<"detectorMask: "<<row->GetDetectorMask()<<endl;
  output <<"partitionName: "<<row->GetPartitionName()<<endl;
  output <<"schemaVersion: "<<row->GetSchemaVersion()<<endl;
  output <<"dataVersion: "<<row->GetDataVersion()<<endl;
  output <<"baseVersion: "<<row->GetBaseVersion()<<endl;

  if(runType == "ADCalib"){
    output <<"\n"<<"------- AD calibration run information ----------"<<endl;
    const DaqCalibRunInfo* calib = m_daqRunInfoSvc->GetCalibRunInfo(runNo);
    if(!calib){
      warning() << "runNo: "<<runNo<<" is not found in table DaqCalibRunInfo."<<endreq;
    }
    else{
      output<<"AdNo: "<<calib->GetAdNo()<<endl;
      output<<"DetectorId: "<<calib->GetDetectorId()<<endl;
      output<<"SourceIdA: "<<calib->GetSourceIdA()<<endl;
      output<<"ZPositionA: "<<calib->GetZPositionA()<<" cm"<<endl;
      output<<"SourceIdB: "<<calib->GetSourceIdB()<<endl;
      output<<"ZPositionB: "<<calib->GetZPositionB()<<" cm"<<endl;
      output<<"SourceIdC: "<<calib->GetSourceIdC()<<endl;
      output<<"ZPositionC: "<<calib->GetZPositionC()<<" cm"<<endl;
      output<<"Duration: "<<calib->GetDuration()<<" seconds"<<endl;
      output<<"LedNumber1: "<<calib->GetLedNumber1()<<endl;
      output<<"LedNumber2: "<<calib->GetLedNumber2()<<endl;
      output<<"LedVoltage1: "<<calib->GetLedVoltage1()<<" mV"<<endl;
      output<<"LedVoltage2: "<<calib->GetLedVoltage2()<<" mV"<<endl;
      output<<"LedFreq: "<<calib->GetLedFreq()<<" Hz"<<endl;
      output<<"LedPulseSep: "<<calib->GetLedPulseSep()<<endl;
      output<<"LtbMode: "<<calib->GetLtbMode()<<endl;
    }        
  }

  if(vec.size() != 0){ //print run configuration into file
    // print hardware lists
    vector<string> dets;
    vector<string> crates;
    vector<string> ltbs; 
    vector<string> fees;
    vector<string>::iterator itr;
    m_daqRunInfoSvc->findHWList(vec, dets, crates, ltbs, fees);
    output <<"\n"<<"------ Hardware setting --------------"<<endl;
    output<<"Detector:";
    for( itr = dets.begin(); itr != dets.end(); itr++)
      output << " "<<(*itr);
    output << endl;
    output <<"Crates:";
    for( itr = crates.begin(); itr != crates.end(); itr++)
      output << " "<<(*itr);
    output << endl;
    output <<"LTBs:";
    for( itr = ltbs.begin(); itr != ltbs.end(); itr++)
      output <<" "<<(*itr);
    output << endl;
    output <<"FEEs:";
    for( itr = fees.begin(); itr != fees.end(); itr++)
      output << " "<<(*itr);
    output << endl;
    
    // print FEE thresholds
    map<string, string> th_obj;
    map<string, bool> isuni;
    map<string, int> thrsuni;
    map<string, map<string, int> > thrs;
    m_daqRunInfoSvc->findFEEThreshold(runType, vec, fees, th_obj, isuni, thrsuni, thrs);
    map<string, string>::iterator obj;
    for( obj = th_obj.begin(); obj != th_obj.end(); obj++){
      if( isuni[obj->first]){
        output<<"Uniform thresholds for channels of "<<obj->first<<" ["<<obj->second<<"]: "<<thrsuni[obj->first]<<endl;
      }
      else
      {
        //output<<"Please check the threshold for each channel of "<<obj->first<<" ["<<obj->second<<"] at the table below."<<endl; 
        output<<"Thresholds for channels of "<<obj->first<<" ["<<obj->second<<"] is as following. "<<endl;
        map<string, int> amap;
        map<string, int>::iterator itt;
        amap = thrs[obj->first];
        for (itt = amap.begin(); itt != amap.end(); itt++){
          output<<"Threshold for one chaneel: "<<itt->first<<" = "<<itt->second<<endl;
        }
      }
    }
    
    // print LTB trigger setting
    map<string, string> mode;
    map<string, int> trig_src;
    map<string, string> trig_name;
    map<string, map<string, int> > trig_thres;
    m_daqRunInfoSvc->findLtbTrigger(runType, vec, ltbs, mode, trig_src, trig_name, trig_thres); 

    map<string, string>::iterator mapit;
    for( mapit = mode.begin(); mapit != mode.end(); mapit++){
      output<<"LTB object: "<<mapit->first<<endl;
      output<<"LTB work mode: "<<mapit->second<<endl;
      output<<"LTB_triggerSource: "<<trig_src[mapit->first]<<endl;
      output<<"LTB trigger setting: "<<trig_name[mapit->first]<<endl;
      map<string, int> thres;
      thres = trig_thres[mapit->first];
      map<string, int>::iterator th;
      for(th = thres.begin(); th != thres.end(); th++){
        if(th->first.find("MULTIPLICITY") != string::npos)
          output<<th->first<<" threshold >= "<<th->second + 1<<endl;
        else
          output<<th->first<<" threshold > "<<th->second<<endl;
      }
    }

    // print the list of run configuration
    output<<"\n"<<"------ Table of run configuration ------------"<<endl;
    output<<"dataVersion  "<<"creationVersion  "<<"className  "<<"objectId  "<<"name  "<<"parentPosition  "<<"intValue  "<<"floatValue  "<<"stringValue  "<<endl;
    vector<DayaBay::DaqRunConfig>::iterator it;
    for( it = vec.begin(); it != vec.end(); it++){
      output<<it->m_DataVersion<<"  "
            <<it->m_CreationVersion<<"  "
            <<it->m_ClassName<<"  "
            <<it->m_ObjectId<<"  "
            <<it->m_Name<<"  "
            <<it->m_ParentPosition<<"  "
            <<it->m_IntValue<<"  "
            <<it->m_FloatValue<<"  "
            <<it->m_StringValue<<endl;
    }
    output.close();
  }

  return StatusCode::SUCCESS;
}

---