src/data_quality/MDQ_FADCBufferOverflow.cpp

Go to the documentation of this file.

/* Standard includes */
#include <stdio.h>
#include <stdlib.h>
#include <string>
#include <map>
#include <utility>

/* MIDAS includes */
#include "midas.h"

/* ROOT includes */
#include <TH1.h>
#include <TH2.h>
#include <TDirectory.h>

/* AlCap includes */
#include "TOctalFADCIsland.h"
#include "TOctalFADCBankReader.h"
#include "TGlobalData.h"
#include "TSetupData.h"

using std::string;
using std::map;
using std::vector;
using std::pair;

/*-- Module declaration --------------------------------------------*/
INT  MDQ_FADCBufferOverflow_init(void);
INT MDQ_FADCBufferOverflow_BOR(INT run_number);
INT MDQ_FADCBufferOverflow_EOR(INT run_number);
INT  MDQ_FADCBufferOverflow(EVENT_HEADER*, void*);

extern HNDLE hDB;
extern TGlobalData* gData;
extern TSetupData* gSetup;

static TH1* hDQ_FADCBufferOverflow_Total;
static TH1* hDQ_FADCBufferOverflow_Fraction;
static TH2* hDQ_FADCBufferOverflow_TotalByEvent;
static int n_total_midas_events;


ANA_MODULE MDQ_FADCBufferOverflow_module =
{
  "MDQ_FADCBufferOverflow",        /* module name           */
  "Andrew Edmonds",              /* author                */
  MDQ_FADCBufferOverflow,          /* event routine         */
  MDQ_FADCBufferOverflow_BOR,      /* BOR routine           */
  MDQ_FADCBufferOverflow_EOR,      /* EOR routine           */
  MDQ_FADCBufferOverflow_init,     /* init routine          */
  NULL,                          /* exit routine          */
  NULL,                          /* parameter structure   */
  0,                             /* structure size        */
  NULL,                          /* initial parameters    */
};

INT MDQ_FADCBufferOverflow_init()
{
  // See if the DataQuality_LowLevel/ directory already exists
  if (!gDirectory->Cd("DataQuality_LowLevel")) {
    
    std::string dir_name("DataQuality_LowLevel/");
    gDirectory->mkdir(dir_name.c_str());
    gDirectory->Cd(dir_name.c_str());
  }

  // Get run number
  char key_name[80];
  int size,run_number;
  sprintf(key_name, "Runinfo/Run number");
  size=sizeof(int);
  db_get_value(hDB,0,key_name, &run_number, &size, TID_INT,1);

  // Create the histograms
  hDQ_FADCBufferOverflow_Total = new TH1F(
    "hDQ_FADCBufferOverflow_Total",
    Form("Total number of MIDAS Events with buffer overflows per board, run %d",run_number),
    4,128, 132);
  hDQ_FADCBufferOverflow_Total->SetBit(TH1::kCanRebin);
  hDQ_FADCBufferOverflow_Total->GetXaxis()->SetTitle("FADC Board Number");
  hDQ_FADCBufferOverflow_Total->GetYaxis()->SetTitle("Number of events with buffer overflows");

  hDQ_FADCBufferOverflow_Fraction = new TH1F(
    "hDQ_FADCBufferOverflow_Fraction",
    Form("Fraction of MIDAS Events with FADC buffer overflows by board, run %d",run_number),
    4,128, 132);
  hDQ_FADCBufferOverflow_Fraction->SetBit(TH1::kCanRebin);
  hDQ_FADCBufferOverflow_Fraction->GetXaxis()->SetTitle("FADC Board Number");
  hDQ_FADCBufferOverflow_Fraction->GetYaxis()->SetTitle("Fraction of events with buffer overflows");

  hDQ_FADCBufferOverflow_TotalByEvent = new TH2F(
    "hDQ_FADCBufferOverflow_TotalByEvent",
    Form("FADC buffer overflows by board and event, run %d",run_number),
    4,128, 132,5000,0,5000);
  hDQ_FADCBufferOverflow_TotalByEvent->SetBit(TH1::kCanRebin);
  hDQ_FADCBufferOverflow_TotalByEvent->GetXaxis()->SetTitle("FADC Board Number");
  hDQ_FADCBufferOverflow_TotalByEvent->GetYaxis()->SetTitle("MIDAS event");

  gDirectory->Cd("/MidasHists/");

  return SUCCESS;
}

INT MDQ_FADCBufferOverflow_BOR(INT run_number)
{
  n_total_midas_events = 1;
  
  return SUCCESS;
}

INT MDQ_FADCBufferOverflow(EVENT_HEADER *pheader, void *pevent)
{
  // Get the event number
  int midas_event_number = pheader->serial_number;

  // Get the NBUF bank and get a pointer to the raw data
  unsigned int* raw; // points at the raw data
  int bankSize = bk_locate(pevent,"NBUF",&raw);

  int board_number = 0;
  int n_buffers_overflow = 0;

  bool Boards[256]; // We can have maximally 256 FADC boards due to the address space
  for(int i=0; i<256; ++i) Boards[i] = false;

  for (int i = 0; i < bankSize; i++) {
    board_number = *(raw+i);

    if(*(raw+i) >= 0 && *(raw+i) < 255) Boards[i] = true;

    n_buffers_overflow++;
    //printf("Event #%d: Board %d n buffers overflowed %d\n", midas_event_number, board_number,n_buffers_overflow);

    // Fill Diagnostic histogram
    hDQ_FADCBufferOverflow_Total->Fill(board_number, 1);   
    hDQ_FADCBufferOverflow_Fraction->Fill(board_number, 1);
    hDQ_FADCBufferOverflow_TotalByEvent->Fill(board_number,midas_event_number, 1);

  }

  /*
  for(int i=0; i<256; ++i){
    if(Boards[i] == true){
      hDQ_FADCBufferOverflow_Fraction->Scale(1.0 / n_total_midas_events);
    }
  }
  */  

  ++n_total_midas_events;

  return SUCCESS;
}

INT MDQ_FADCBufferOverflow_EOR(INT run_number)
{

  //for(int i=0; i<256; ++i){
    //if(Boards[i] == true){
      hDQ_FADCBufferOverflow_Fraction->Scale(1.0 / n_total_midas_events);
    //}
  //}

  return SUCCESS;
}

---