CalibEvtQual.cc

Go to the documentation of this file.

#include <interface_main/SQEvent.h>
#include <fun4all/Fun4AllReturnCodes.h>
#include <fun4all/Fun4AllHistoManager.h>
#include <phool/PHNodeIterator.h>
#include <phool/PHIODataNode.h>
#include <phool/getClass.h>
#include "CalibEvtQual.h"
using namespace std;

CalibEvtQual::CalibEvtQual(const std::string& name) : SubsysReco(name)
{
  ;
}

int CalibEvtQual::Init(PHCompositeNode* topNode)
{
  return Fun4AllReturnCodes::EVENT_OK;
}

int CalibEvtQual::InitRun(PHCompositeNode* topNode)
{
  return Fun4AllReturnCodes::EVENT_OK;
}

int CalibEvtQual::process_event(PHCompositeNode* topNode)
{
  SQEvent* event = findNode::getClass<SQEvent>(topNode, "SQEvent");
  if (! event) return Fun4AllReturnCodes::ABORTEVENT;

  int run_id = event->get_run_id();
  int qual   = event->get_data_quality();

  int n_tdc; // expected number of tdc info
  if (run_id >= 28664) n_tdc = 101; // seen in run 28700
  else                 n_tdc = 108; // seen in run 28000
  // This run range was confirmed by checking n_tdc in the following runs;
  //  * 101 in runs 28680, 28670, 28665, 28664
  //  * 108 in runs 28500, 28600, 28650, 28660, 28661
  // Note that runs 28662 & 28663 contain no event.
  // Is the change of n_tdc related to elog #17385??
  // https://e906-gat6.fnal.gov:8080/SeaQuest/17385
  if      (event->get_n_board_taiwan() != n_tdc) qual |= ERR_N_TDC;
  if      (event->get_n_board_v1495      () < 2) qual |= ERR_N_V1495_0;
  else if (event->get_n_board_v1495      () > 2) qual |= ERR_N_V1495_2;
  if      (event->get_n_board_trig_bit   () < 1) qual |= ERR_N_TRIGB_0;
  else if (event->get_n_board_trig_bit   () > 1) qual |= ERR_N_TRIGB_2;
  if      (event->get_n_board_trig_count () < 1) qual |= ERR_N_TRIGC_0;
  else if (event->get_n_board_trig_count () > 1) qual |= ERR_N_TRIGC_2;
  if      (event->get_n_board_qie        () < 1) qual |= ERR_N_QIE_0;
  else if (event->get_n_board_qie        () > 1) qual |= ERR_N_QIE_2;

  event->set_data_quality(qual);
  //if (qual != 0) {
  //  cout << "  N! " << evt_id << " | " << ed->n_qie << " " << ed->n_v1495 << " " << ed->n_tdc << " " << ed->n_trig_b << " " << ed->n_trig_c << endl;
  //  }
  //}
  
  return Fun4AllReturnCodes::EVENT_OK;
}

int CalibEvtQual::End(PHCompositeNode* topNode)
{
  return Fun4AllReturnCodes::EVENT_OK;
}

void CalibEvtQual::PrintEvent(SQEvent* evt)
{
  cout << "SQEvent:  "
       << "  " << evt->get_run_id       ()
       << "  " << evt->get_spill_id     ()
       << "  " << evt->get_event_id     ()
       << "  " << evt->get_coda_event_id()
       << "  " << evt->get_data_quality ()
       << "  " << evt->get_vme_time     ()
       << "\n";
}